Zernike不变矩 MATLAB实现

Zernike矩一种新的图像识别特征，可以用在很多方面

function [iimg]=guiyihua（fbeta）

bimg	= double（f）;

[rowcol]	= size（bimg）;
[xy]   	= meshgrid（1:col1:row）;


m00	= gmoment（bimg00）;
m10	= gmoment（bimg10）;
m01	= gmoment（bimg01）;
if  nargin==1
    inva=1;
else
    inva= 1/sqrt（beta/m00）;
end
%计算重心
thecen1= m10/m00;
thecen2= m01/m00;
xi	= x*inva+thecen1;
yi	= y*inva+thecen2;
%xi	= x*inva;
%yi	= y*inva;
%xi	= （x+thecen1）*inva;
%yi	= （y+thecen2）*inva;
iimg	= interp2（x+inva*col/2+.5y+inva*row/2+.5bimgxiyi‘*linear‘）;
wnan		= find（isnan（iimg）==1）;
iimg（wnan）	= zeros（size（wnan））;

function	[m_pq]=gmoment（imgpq）
if nargin~=3
	p=0;q=0;
end
[rowcol]	= size（img）;
[xy]   	= meshgrid（1:col1:row）;
x		= x.^p;
y		= y.^q;

gmom		= x.*y.*double（img）;
m_pq		= sum（sum（gmom））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        777  2010-06-23 09:49  Zernike矩\guiyihua.m

     文件       1039  2010-05-04 01:37  Zernike矩\zernike.m

     目录          0  2010-07-20 20:22  Zernike矩

----------- ---------  ---------- -----  ----

                 1816                    3