clear?all
clc
I1=imread（‘test1.bmp‘）;
%A=[];
tn=36;
for?r=1:tn
????k=int2str（r）;
c=strcat（‘test‘k‘.BMP‘）;
I1=imread（c）;
if?r==1
????im=figure;
?[A1rect]=imcrop（I1）;
end
A1=imcrop（I1rect）;
I3=im2double（A1）;?
I=I3;
[mn]=size（I）;
b（r）=0;
for?i=2:（m-1）
???for?j=2:（n-1）
??????a1=[I（i+1j-1）+I（i+1j）+I（i+1j+1）]-[I（i-1j-1）+I（i-1j）+I（i-1j+1）];
?????a2=[I（i-1j+1）+I（ij+1）+I（i+1j+1）]-[I（i-1j-1）+I（ij-1）+I（i+1j-1）];
????????b（r）=b（r）+a1*a1+a2*a2;
???end
end
end
A=[b（1）b（2）b（3）b（4）b（5）b（6）b（7）b（8）b（9）b（10）b（11）b（12）b（13）b（14）b（15）b（16）b（17）b（18）b（19）b（20）b（21）b（22）b（23）b（24）b（25）b（26）b（27）b（28）b（29）b（30）b（31）b（32）b（33）b（34）b（35）b（36）];
x=1:tn;
plot（xA）;
[ps]=polyfit（xA20）;
y1=polyval（px）;
plot（xA‘or‘xy1）;
grid?on;
xlabel（‘鏡頭位置‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????820??2013-05-28?13:17??brenner.m