定焦评价函数程序.doc

定焦评价函数程序一边沿检测算子1.Robert算子clear all;close all;clc;N=11;for i=1:Nj=imread( int2str(i) .tif);c=222 272 300 270 221 194; r=21 21 75 121 121 75; bw=roipoly(j,c,r); h=fspecial(unsharp); f=roifilt2(h,j,bw);h1=0 -1;1 0;h2=-1 0;0 1;ff=filter2(h1,f)+ filter2(h2,f);f=sqrt(double(ff.2);f_rob(i)=sum(sum(f);fprintf(f-rob(%d)=%fn,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_rob maxindex=max(f_rob)2.Sobel算子clear all;close all;clc;N=11;for i=1:Nj=imread( int2str(i) .tif);c=222 272 300 270 221 194; r=21 21 75 121 121 75; bw=roipoly(j,c,r); h=fspecial(unsharp); f=roifilt2(h,j,bw);h1=-1 0 -1;-2 0 -2;-1 0 -1;h2=-1 -2 -1;0 0 0;1 2 1;ff=filter2(h1,f)+ filter2(h2,f);f=sqrt(double(ff.2);f_sob(i)=sum(sum(f);fprintf(f_sob(%d)=%fn,i,f_sob(i);endv_max=max(f_sob); f_sob=f_sob/v_max;figure;plot(f_sob ,-k+);xlabel(每个20um取一幅图像); ylabel(归一化定焦函数评价值); title(sobel算子做定价评价函数);maxf_sob maxindex=max(f_sob)3. Laplacion算子4领域clear all;close all;clc;N=11;for i=1:Nj=imread( int2str(i) .tif);c=222 272 300 270 221 194; r=21 21 75 121 121 75; bw=roipoly(j,c,r); h=fspecial(unsharp); f=roifilt2(h,j,bw);h=0 1 0;1 -4 1; 0 1 0;lf=filter2(h,f);f_lap(i)=sum(sum(lf.2);fprintf(f_lap(%f)=%fn,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_lap maxindex=max(f_lap)Laplacion算子8领域clear all;close all;clc;N=11;for i=1:Nj=imread( int2str(i) .tif);c=222 272 300 270 221 194; r=21 21 75 121 121 75; bw=roipoly(j,c,r); h=fspecial(unsharp); f=roifilt2(h,j,bw);h=1 1 1;1 -8 1; 1 1 1;lf=filter2(h,f);f_lap(i)=sum(sum(lf.2);fprintf(f_lap(%f)=%fn,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_lap maxindex=max(f_lap)4.log算子clear all;close all;clc;N=37;for i=1:Nj=imread( int2str(i) .bmp);c=222 272 300 270 221 194; r=21 21 75 121 121 75; bw=roipoly(j,c,r); h=fspecial(unsharp); f=roifilt2(h,j,bw);h=-2 -4 -4 -4 -2;-4 0 8 0 -4;-4 8 24 8 -4;-4 0 8 0 -4;-2 -4 -4 -4 -2;lf=filter2(h,f);f_log(i)=sum(sum(lf.2);fprintf(f_log(%d)=%fn,i,f_log(i);endv_max=max(f_log); f_log=f_log/v_max;figure;plot(f_log, -k+);xlabel(每个20um取一幅图像); ylabel(归一化定焦函数评价值); title(log算子做定价评价函数);maxf_log maxindex=max(f_log)二灰度梯度函数clear all;close all;clc;N=11;for i=1:Nf=imread( int2str(i) .tif);c=222 272 300 270 221 194; r=21 21 75 121 121 75; bw=roipoly(j,c,r); h=fspecial(unsharp); f=roifilt2(h,j,bw); M,N=size(f);F(i)=sum(sum(sqrt(double(f(1:M-1,2:N)-f(1:M-1,1:N-1).2+. (f(2:M,1:N-1)-f(1:M-1,1:N-1).2);fprintf(F(%d)=%fn,i,F(i);endv_max=max(F); F=F/v_max;figure;plot(F, -k+);xlabel(每隔0.1mm取一幅图像); ylabel(归一化定焦函数评价值); title(灰度梯度做定焦评价函数);F_max index=max(F)三 频域函数构造评价函数clear all;close all;clc;N=9;for i=1:Nj=imread( int2str(i) .tif);c=222 272 300 270 221 194; r=21 21 75 121 121 75; bw=roipoly(j,c,r); h=fspecial(unsharp); f=roifilt2(h,j,bw); h=fft2(f); h1=real(h)