基于Matlab的红外图像处理算法研究 附录：源程序.docVIP

源程序： 1、线性变换增加对比度： X1=imread(pout.tif); figure,imshow(X1) f0=0;g0=0; f1=70;g1=30; f2=180;g2=230; f3=255;g3=255; figure,plot([f0,f1,f2,f3],[g0,g1,g2,g3]) axis tight,xlabel(f),ylabel(g) title(intensity transformation) r1=(g1-g0)/(f1-f0); b1=g0-r1*f0; r2=(g2-g1)/(f2-f1); b2=g1-r2*f1; r3=(g3-g2)/(f3-f2); b3=g2-r3*f2; [m,n]=size(X1); X2=double(X1); for i=1:m for j=1:n f=X2(i,j); g(i,j)=0; if (f=0)(f=f1) g(i,j)=r1*f+b1; elseif (f=f1)(f=f2) g(i,j)=r2*f+b2; elseif (f=f2)(f=f3) g(i,j)=r3*f+b3; end end end figure,imshow(mat2gray(g)) f0=0,g0=0; f1=10,g1=5; f2=160,g2=240 f3=255,g3=255 figure,plot([f0,f1,f2,f3],[g0,g1,g2,g3]) axis tight,xlabel(f),ylabel(g) title(intensity transformation) 运用Matlab函数 X1=imread(pout.tif); figure,imshow(X1) J=imadjust(X1,[0.25,0.6],[],1.2322) figure,imshow(J) 2、求反 X1=imread(pout.tif); figure,imshow(X1) f1=200; g1=256; figure,plot([0,f1],[g1,0]) axis tight,xlabel(f),ylabel(g) title(intensity transformation) k=g1/f1; [m,n]=size(X1); X2=double(X1); for i=1:m for j=1:n f=X2(i,j); g(i,j)=0; if(f=0)(f=f1) g(i,j)=g1-k*f; else g(i,j)=0; end end end figure,imshow(mat2gray(g)) 3、对数变换处理 X1=imread(liftingbody.png); figure,imshow(X1); c=255/log(256); x=0:255; y=c*log(1+x); figure,plot(x,y) axis tight,xlabel(f),ylabel(g) title(intensity transformation) [m,n]=size(X1); X2=double(X1); for i=1:m for j=1:n g(i,j)=c*log(X2(i,j)+1); end end figure,imshow(mat2gray(g)) 运用Matlab函数 X1=imread(circuit.tif); figure,imshow(X1) J=imadjust(X1,[],[0.1 0.8],0.35287) figure,imshow(J) 4、直方图均衡化 I=imread(tire.tif); J=histeq(I); imshow(I) figure,imshow(J) figure,imhist(J,64) figure,imhist(I,64) 5、直方图规定化 I=imread(tire.tif); hgram=0:255 J=histeq(I,hgram); imshow(I) figure,imshow(J) figure,imhist(I,64) figure,imhist(J,64) 6、平滑滤波 均值平滑滤波 I=imread(moon.tif); J=imnoise(I,salt,0.02); imshow(I) figure,imshow(J) K1=f