直接采用循环实现：.doc

直接采用循环实现： ? %?计算图像平面上区域的矩%?参考：《图像工程（中册）图像分析》P240%?f???????????输入二维图像(函数)%?(p,q)????需要计算的阶数%?返回：%????????函数?f(x,y)?的?p+q?阶矩function?[moment]=CalculateMoment(f,p,q)[xSize,ySize]=size(f);????????%?计算横纵范围f=double(f);????????????????????????%?类型抓换moment=0.0;??tictesttimes=10000;????????????%?测试次数for?mm=1:testtimes????for?x=1:xSize????????for?y=1:ySize????????????moment=moment+(x^p)*(y^q)*f(x,y);??%?计算代码????????end????endendmoment=moment/testtimestoc 输出为： moment = 1.2885e+003 Elapsed time is 12.563000 seconds.