function?zernike7（I）

I=imread（‘2.jpg‘）;

%?7*7Zernike模板
M00=...
[
???????0	0.0287	0.0686	0.0807	0.0686	0.0287	0
????0.0287	0.0815	0.0816	0.0816	0.0816	0.0815	0.0287
????0.0686	0.0816	0.0816	0.0816	0.0816	0.0816	0.0686
????0.0807	0.0816	0.0816	0.0816	0.0816	0.0816	0.0807
????0.0686	0.0816	0.0816	0.0816	0.0816	0.0816	0.0686
????0.0287	0.0815	0.0816	0.0816	0.0816	0.0815	0.0287
???????0	0.0287	0.0686	0.0807	0.0686	0.0287	0
];
M11R=...
[
???????0	-0.015	-0.019	0	0.019	0.015	??0
????-0.0224	-0.0466	-0.0233	0	0.0233	0.0466	0.0224
????-0.0573	-0.0466	-0.0233	0	0.0233	0.0466	0.0573
????-0.069	-0.0466	-0.0233	0	0.0233	0.0466	0.069
????-0.0573	-0.0466	-0.0233	0	0.0233	0.0466	0.0573
????-0.0224	-0.0466	-0.0233	0	0.0233	0.0466	0.0224
???????0	-0.015	-0.019	0	0.019	0.015	??0
];
M11I=...
[
???????0	-0.0224	-0.0573	-0.069	-0.0573	-0.0224	??0
????-0.015	-0.0466	-0.0466	-0.0466	-0.0466	-0.0466	-0.015
????-0.019	-0.0233	-0.0233	-0.0233	-0.0233	-0.0233	-0.019
???????0	???0	???0	???0	???0	???0	??0
????0.019	0.0233	0.0233	0.0233	0.0233	0.0233	0.019
????0.015	0.0466	0.0466	0.0466	0.0466	0.0466	0.015
???????0	0.0224	0.0573	0.069	0.0573	0.0224	??0
];
M20=...
[
???????0	0.0225	0.0394	0.0396	0.0394	0.0225	0
????0.0225	0.0271	-0.0128	-0.0261	-0.0128	0.0271	0.0225
????0.0394	-0.0128	-0.0528	-0.0661	-0.0528	-0.0128	0.0394
????0.0396	-0.0261	-0.0661	-0.0794	-0.0661	-0.0261	0.0396
????0.0394	-0.0128	-0.0528	-0.0661	-0.0528	-0.0128	0.0394
????0.0225	0.0271	-0.0128	-0.0261	-0.0128	0.0271	0.0225
???????0	0.0225	0.0394	0.0396	0.0394	0.0225	0
];


if?length（size（I））==3?I=rgb2gray（I）;?end
I=im2bw（I0.6）;
K=double（I）;
[m?n]=size（K）;
xs=double（zeros（mn））;
ys=double（zeros（mn））;

%?卷積運算
A11I=conv2（M11IK）;
A11R=conv2（M11RK）;
A20=conv2（M20K）;

%?截掉多余部分
A11I=A11I（4:end-34:end-3）;
A11R=A11R（4:end-34:end-3）;
A20=A20（4:end-34:end-3）;

J=zeros（size（K））;
boundary=J;
theta=atan2（A11IA11R）;%計算theta

%計算邊緣的三個參數
A11C=A11R.*cos（theta）+A11I.*sin（theta）;
l=A20./A11C;
k=1.5*A11C./（（1-l.^2）.^1.5）;
e=abs（l）>1/3.5;
k（e）=0;

%邊緣判斷條件
a=abs（l）<1/sqrt（2）*2/7;
b=abs（k）>max（I（:））/10;
%?ab分別為距離和邊緣強度判斷結果
J（a&b）=1;
%將圖像的最邊緣去除
%?boundary（2:end-12:end-1）=1;
%?J（~boundary）=0;
format?short
%?[xy]=find（J==1）;%邊緣的像素級坐標
%?O=[x?y];
%?Z=[x+l（find（J==1））.*cos（theta（find（J==1）））?y+l（find（J==1））.*sin（theta（find（J==1）））];%亞像素坐標
%?%?fprintf（‘%.4f?%.4f\n‘Z‘）;

[Lnum]=bwlabel（J8）;%對二值圖像進行標記
%自動化搜索連通域
s=zeros（1num）;
for?i=1:num
???s（i）=size（find（L==i）1）;
end
[bwLlabel]=sort（s‘descend‘）;
if?label（1）????index1=label（1）;
????index2=label（2）;
else
????index1=label（2）;
????index2=label（1）;
end

%計算左邊探針的最前端坐標
[r1c1]=find（L==index1）;
A1=[r1?c1];
y1=max（A1（:2））;%該連通域中y最大值為針尖處
x1=max（A1（find（A1（:2）==y1）1））;
x1sub=x1+3.5*l（x1y1）*cos（theta（x1y1））;
y1sub=y1+3.5*l（x1y1）*sin（theta（x1y1））;

%計算最右邊探針的最前端坐標
[r2c2]=find（L==index2）;
A2=[r2?c2];
y2=min（A2（:2））;%該連通域中y最小為連通域
x2=max（A2（find（A2（:2）==

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????83415??2015-01-13?19:34??2.jpg
?????文件????????3585??2015-01-19?09:21??zernike7.m