zernike 矩不变量

矩不变量，适用于图像匹配，数字水印，图像重建

function	[A_nmzmlistcidxV_nm]	= zernike（imgnm）
if nargin>0
if nargin==1
	n	= 0;
end
d		= size（img）;			
img		= double（img）;
xstep		= 2/（d（1）-1）;
ystep		= 2/（d（2）-1）;
[xy]		= meshgrid（-1:xstep:1-1:ystep:1）;
circle1		= x.^2 + y.^2;
inside		= find（circle1<=1）;
mask		= zeros（d）;
mask（inside）	= ones（size（inside））;
[cimgcidx]	= clipimg（imgmask）;
z		= clipimg（x+i*ymask）;
p		= 0.9*abs（z）;   
theta		= angle（z）;
c	= 1;
for order=1:length（n）
	n1	= n（order）;
	if nargin<3
		m	= zpossible（n1）;		
	end
	for r=1:length（m）
		V_nmt		= zpoly（n1m（r）ptheta）;
		zprod		= cimg.*conj（V_nmt）;
		A_nm（c）		= （n1+1）*sum（sum（zprod））/pi;
		zmlist（c1:2）	= [n1 m（r）];
		if nargout==4
			V_nm（:c）	= V_nmt;
		end
		c		= c+1;
	end
end
else
end



function [cimgcindexdim]	= 

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

     文件     705654  2005-09-23 14:16  zernike\1.bmp

     文件     705654  2005-09-23 14:34  zernike\10.bmp

     文件       8254  2006-08-30 11:09  zernike\100.bmp

     文件      56614  2010-11-12 11:06  zernike\1000.bmp

     文件       8254  2006-08-28 09:56  zernike\101.bmp

     文件       8254  2006-08-28 09:13  zernike\102.bmp

     文件       8254  2006-08-28 09:14  zernike\103.bmp

     文件       8254  2006-08-28 09:56  zernike\104.bmp

     文件       8254  2006-08-28 09:59  zernike\105.bmp

     文件       8254  2006-08-28 09:58  zernike\106.bmp

     文件       8254  2006-08-28 09:53  zernike\107.bmp

     文件       8254  2006-09-02 22:37  zernike\108.bmp

     文件       8254  2006-08-28 10:01  zernike\109.bmp

     文件       8254  2006-08-28 09:51  zernike\110.bmp

     文件       8254  2006-08-28 09:50  zernike\111.bmp

     文件       8254  2006-08-28 09:50  zernike\112.bmp

     文件       8254  2006-08-28 10:04  zernike\113.bmp

     文件       8254  2006-08-28 10:07  zernike\114.bmp

     文件       8254  2006-08-28 10:08  zernike\115.bmp

     文件       8254  2006-08-28 10:09  zernike\116.bmp

     文件     705654  2005-09-23 14:44  zernike\15.bmp

     文件     705654  2005-09-23 14:48  zernike\16.bmp

     文件     705654  2005-09-23 14:58  zernike\17.bmp

     文件     705654  2005-09-23 15:00  zernike\18.bmp

     文件       8254  2006-08-30 10:11  zernike\201.bmp

     文件       8254  2006-08-30 21:02  zernike\202.bmp

     文件       8254  2006-08-30 21:03  zernike\203.bmp

     文件       8254  2006-08-30 21:04  zernike\204.bmp

     文件     705654  2005-09-23 15:08  zernike\22.bmp

     文件     705654  2005-09-23 15:27  zernike\23.bmp

............此处省略47个文件信息