genrs.m

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

function [R,S]=genRS(N,r)
%generates matrices R and S for the calculation of image
%moments based on Fluss and Tuk's method described in 
%"On the Calculation of Moments"

%N-expected image size
%r-maximum index of moments needed

pm=r+1;
R1=ones(pm,1)*linspace(1,N,N);
R2=linspace(0,r,pm)'*ones(1,N);
R=R1.^R2;
S=R*triu(ones(N));