В моем предыдущем ответе я написал реализацию для подмножества функции regionprops
.Цель состояла в том, чтобы найти ориентацию изображения, которая была получена из моментов изображения .Вот часть, относящаяся к вам:
function outmom = raw_moments(im,i,j)
outmom = sum(sum( ((1:size(im,1))'.^j * (1:size(im,2)).^i) .* im ));
end
function cmom = central_moments(im,i,j)
rawm00 = raw_moments(im,0,0);
centroids = [raw_moments(im,1,0)/rawm00 , raw_moments(im,0,1)/rawm00];
cmom = sum(sum( (([1:size(im,1)]-centroids(2))'.^j * ...
([1:size(im,2)]-centroids(1)).^i) .* im ));
end
Код соответствует уравнениям из статьи Википедии, поэтому никаких дополнительных объяснений не требуется ..