momentbyrs.m

fast calculation of image moments based on Fluss and Tuk s method described in "On the Calculat

function m=momentbyRS(B,p,q,R,S)
%fast calculation of image
%moments based on Fluss and Tuk's method described in 
%"On the Calculation of Moments"

%B - binary image matrix
%R,S - previously generated matrices

d=filter2([1,-1],B);
[i,j]=find(d==-1);
[k,l]=find(d==1);
m=S(p+1,l)*R(q+1,k)'-S(p+1,j)*R(q+1,i)';