function CtrM = CentralM(p,q,I)
%
% CentralM 计算二维p+q阶中心矩
%
% I(x, y)是图像在点(x, y)的强度(灰度值)。
[m,n] = size(I);
CtrM = 0;
m10 = Moment2D(1,0,I);
m01 = Moment2D(0,1,I);
m00 = Moment2D(0,0,I);
xaverage = m10/m00;
yaverage = m01/m00;
for x = 1:m
for y= 1:n
CtrM = CtrM+(x-xaverage)^p * (y-yaverage)^q * I(x,y);
end;
end;
%CtrM ,p,q ,xaverage,yaverage%