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
%計(jì)算重心
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））;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????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