%???感謝親親使用此代碼，此代碼解決您的問題了嗎~（@^_^@）~
%???沒解決的話告訴親親一個好消息，我這里可以1分錢幫助代碼改錯，還提供1分錢成品代碼（′▽‘〃）哦~
%???登錄淘寶店鋪“大成軟件工作室”便可領取
%???是的，親親真的沒有看錯，撓破頭皮的問題真的1分錢就可以解決了\（^o^）/YES!
%???小的這就把傳送門給您，記得要收藏好哦（づ￣3￣）づ╭～
%???傳送門：https://item.taobao.com/item.htm?spm=a1z10.1-c.w4004-15151018122.5.uwGoq5&id=538759553146
%???如果傳送門失效，親親可以來店鋪討要，客服MM等親親來騷擾哦~（*/ω╲*）
function?[oimgXY]=baseZernike1（img）
load?M2;?????%獲取模板7*7的
%卷積運算
Z00=conv2（M00img）;
Z11I=conv2（M11Iimg）;
Z11R=conv2（M11Rimg）;
Z20=conv2（M20img）;
Z31I=conv2（M31Iimg）;
Z31R=conv2（M31Rimg）;
Z40=conv2（M40img）;
%截掉多余部分
Z00=Z00（4:end-34:end-3）;
Z11I=Z11I（4:end-34:end-3）;
Z11R=Z11R（4:end-34:end-3）;
Z20=Z20（4:end-34:end-3）;
Z31I=Z31I（4:end-34:end-3）;
Z31R=Z31R（4:end-34:end-3）;
Z40=Z40（4:end-34:end-3）;
%設置k，oimg，theta，Zz11，l等的初始矩陣，以防后來矩陣運算大小不匹配。
k=zeros（size（img1）size（img2））;
oimg=k;
theta=k;
theta31=k;
Zz11=k;
Zz31=k;
z11=k;
z31=k;
Zz40=k;
Zz20=k;
h=k;
bounder=k;
l=100*ones（size（img1）size（img2））;
l1=100*ones（size（img1）size（img2））;
l2=100*ones（size（img1）size（img2））;
%?????????防止分子、分母均為0而產生的NAN
a=abs（Z11R）<0.0001;
b=abs（Z11I）<0.0001;
c=a&b;Z11R（c）=1;Z11I（c）=1;
theta=atan（Z11I./Z11R）;
g=abs（Z31R）<0.0001;
j=abs（Z31I）<0.0001;
f=g&j;Z31R（f）=1;Z31I（f）=1;
theta31=atan（Z31I./Z31R）;

Zz11=Z11R.*cos（theta）+Z11I.*sin（theta）;
Zz31=Z31R.*cos（theta31）+Z31I.*sin（theta31）;?
Zz20=Z20;
Zz40=Z40;
l1=sqrt（（5*Zz40+3*Zz20）./（8*Zz20））;
l2=sqrt（（5*Zz31+Zz11）./（6*Zz11））;
l=（l2+l1）/2;
?e=abs（l）>1/（3.5*sqrt（2））;
?k=1.5.*Zz11./（（1-l2.^2）.^1.5）;
?k（e）=0;
?h=（Z00-（k*pi）/2+k.*asin（l2）+k.*l2.*sqrt（1-l2.^2））./pi;
a=abs（l2-l1）<1;
?b=abs（k）>=1.7;
%b=abs（k）>max（img（:））/10;
%?amp=sqrt（Z11R.^2+Z11I.^2）;
%??t=amp>15;
%?ab分別為距離和邊緣強度判斷結果c為之前需要去掉的NAN部分
?%d=a&b&t&~c&~f;
d=a&b&~c&~f;
oimg（d）=1;
bounder（3:end-23:end-2）=1;
oimg（~bounder）=0;

figureimshow（oimg）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-09-23?05:36??baseZernike1\
?????文件????????2116??2016-09-23?10:11??baseZernike1\baseZernike1.m
?????文件????????1255??2009-03-13?09:40??baseZernike1\M2.mat
?????文件??????786872??2009-03-11?16:01??baseZernike1\基于Zernike正交矩的圖像亞像素邊緣檢測算法改進.pdf