function?M=Caculatemask（zN）
%計(jì)算Zerniek模板
%M:模板
%z:徑向多項(xiàng)式
%N:模板大小
%例
%?syms?x?y;
%?M00=caculatemask（x7）;
%?M11=caculatemask（x-y*i7）;
%?M20=caculatemask（2*x^2+2*y^2-17）;
syms?x?y
M=zeros（NN）;
[XY]=meshgrid（1:N1:N）;
xx=（2.*X-N-1）/N;
yy=（N+1-2.*Y）/N;
inc=1/N;
%中間
for?row=2:N-1
????for?col=2:N-1
????????x1=xx（rowcol）-inc;
????????x2=xx（rowcol）+inc;
????????y1=yy（rowcol）-inc;
????????y2=yy（rowcol）+inc;
????????M（rowcol）=vpa（int（int（zyy1y2）xx1x2）4）;
????end
end
%上5
for?col=2:N-1
????x1=xx（1col）-inc;
????x2=xx（1col）+inc;
????y1=yy（1col）-inc;
????y2=sqrt（1-x^2）;
????M（1col）=vpa（int（int（zyy1y2）xx1x2）4）;
end
%下5
for?col=2:N-1
????x1=xx（Ncol）-inc;
????x2=xx（Ncol）+inc;
????y1=-sqrt（1-x^2）;
????y2=yy（Ncol）+inc;
????M（Ncol）=vpa（int（int（zyy1y2）xx1x2）4）;
end
%左5
for?row=2:N-1
????x1=-sqrt（1-y^2）;
????x2=xx（row1）+inc;
????y1=yy（row1）-inc;
????y2=yy（row1）+inc;
????M（row1）=vpa（int（int（zxx1x2）yy1y2）4）;
end
%右5
for?row=2:N-1
????x1=xx（rowN）-inc;
????x2=sqrt（1-y^2）;
????y1=yy（rowN）-inc;
????y2=yy（rowN）+inc;
????M（rowN）=vpa（int（int（zxx1x2）yy1y2）4）;
end

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

?????文件???????1187??2016-01-29?12:59??EdgeDetecion\Caculatemask.m

?????文件???????2327??2016-01-29?16:44??EdgeDetecion\Edgedectect.m

?????文件?????263224??2016-01-26?10:13??EdgeDetecion\Lena.bmp

?????文件????????272??2016-01-29?13:07??EdgeDetecion\Radialpoly.m

?????文件????????141??2016-01-29?21:37??EdgeDetecion\test.m

?????文件????????459??2016-01-29?13:07??EdgeDetecion\Zernikmoment.m

?????目錄??????????0??2016-01-29?21:40??EdgeDetecion

-----------?---------??----------?-----??----

???????????????267610????????????????????7