moment.m

A set of functions to calculate Hu s seven invariant moments. Good accuracy. Written for my Project

function [m] = moment(img, p, q)
% Calculate the regular moment of a binary image with order p, q
% img - binary image
% p, q - order
img = double(img);
[y x] = size(img);
m = zeros(size(p));
	for i = 1:y
		for j = 1:x
			m = m + img(i,j).*(i.^p).*(j.^q);
		end
	end
return