资源描述：
《定焦评价函数程序》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、定焦评价函数程序一．边沿检测算子1.Robert算子clearall;closeall;clc;N=11;fori=1:Nj=imread([int2str(i)'.tif']);c=[222272300270221194];r=[21217512112175];bw=roipoly(j,c,r);h=fspecial('unsharp');f=roifilt2(h,j,bw);h1=[0-1;10];h2=[-10;01];ff=filter2(h1,f)+filter2(h2,f);f=sqrt(double(ff.^2));f_rob(i)=sum(sum(f))
2、;fprintf('f-rob(%d)=%f',i,f_rob(i));endv_max=max(f_rob);f_rob=f_rob/v_max;figure;plot(f_rob,'-k+');xlabel('每个20um取一幅图像');ylabel('归一化定焦函数评价值');title('roberts算子做定价评价函数');[maxf_robmaxindex]=max(f_rob)2.Sobel算子clearall;closeall;clc;N=11;fori=1:Nj=imread([int2str(i)'.tif']);c=[2222723002702
3、21194];r=[21217512112175];bw=roipoly(j,c,r);h=fspecial('unsharp');f=roifilt2(h,j,bw);h1=[-10-1;-20-2;-10-1];h2=[-1-2-1;000;121];ff=filter2(h1,f)+filter2(h2,f);f=sqrt(double(ff.^2));f_sob(i)=sum(sum(f));fprintf('f_sob(%d)=%f',i,f_sob(i));endv_max=max(f_sob);f_sob=f_sob/v_max;figure;plot
4、(f_sob,'-k+');xlabel('每个20um取一幅图像');ylabel('归一化定焦函数评价值');title('sobel算子做定价评价函数');[maxf_sobmaxindex]=max(f_sob)3.Laplacion算子4领域clearall;closeall;clc;N=11;fori=1:Nj=imread([int2str(i)'.tif']);c=[222272300270221194];r=[21217512112175];bw=roipoly(j,c,r);h=fspecial('unsharp');f=roifilt2(h,j,b
5、w);h=[010;1-41;010];lf=filter2(h,f);f_lap(i)=sum(sum(lf.^2));fprintf('f_lap(%f)=%f',i,f_lap(i));endv_max=max(f_lap);f_lap=f_lap/v_max;figure;plot(f_lap,'-k+');xlabel('每个20um取一幅图像');ylabel('归一化定焦函数评价值');title('laplance算子做定价评价函数');[maxf_lapmaxindex]=max(f_lap)Laplacion算子8领域clearall;close
6、all;clc;N=11;fori=1:Nj=imread([int2str(i)'.tif']);c=[222272300270221194];r=[21217512112175];bw=roipoly(j,c,r);h=fspecial('unsharp');f=roifilt2(h,j,bw);h=[111;1-81;111];lf=filter2(h,f);f_lap(i)=sum(sum(lf.^2));fprintf('f_lap(%f)=%f',i,f_lap(i));endv_max=max(f_lap);f_lap=f_lap/v_max;figu
7、re;plot(f_lap,'-k+');xlabel('每个20um取一幅图像');ylabel('归一化定焦函数评价值');title('laplance算子做定价评价函数');[maxf_lapmaxindex]=max(f_lap)4.log算子clearall;closeall;clc;N=37;fori=1:Nj=imread([int2str(i)'.bmp']);c=[222272300270221194];r=[21217512112175];bw=roipoly(j,c,r);h=fspecial('unsharp')
