圖像處理程序

注意：在用這些程序的時(shí)候最好更改一下變量跟圖片，更改圖片的方法用第一個(gè)例子來(lái)說(shuō)明：I1=imread(‘tire.tif’);這個(gè)句子中的tire.tif改成你所用圖片的地址就行了，比如我改成：D:\我的文檔\圖片\MyPictures\haha.tif既用上面圖片的地址替換掉tire.tif即可，不用輸入地址是因?yàn)閠ire.tif這個(gè)圖片在matlab默認(rèn)的圖片收藏文件夾下。如果你把你自己的圖片放到里面也可以直接用你的圖片名稱(chēng)來(lái)代替。1.對(duì)圖像pout.tif加入椒鹽噪聲后，實(shí)現(xiàn)butterworth低通濾波。clear;I1=imread(‘tire.tif’);subplot(2,2,1);imshow(I1);title('原始圖像');I2=imnoise(I1,'salt&pepper',0.02); %加入椒鹽噪聲subplot(2,2,2);imshow(I2);title('噪聲圖像');f=double(I2);g=fft2(f) %傅里葉變換g=fftshift(g) %轉(zhuǎn)換數(shù)據(jù)矩陣[N1,N2]=size(g);n=2;d0=50;n1=fix(N1/2);n2=fix(N2/2);fori=1:N1 %對(duì)頻域中循環(huán)濾波 forj=1:N2 d=sqrt((i-n1)^2+(j-n2)^2); %計(jì)算butterworth低通濾波轉(zhuǎn)換函數(shù) h=1/(1+0.414*(d/d0)^(2*n)); result1(i,j)=h*g(i,j); if(g(i,j)>50) result2(i,j)=0; else result2(i,j)=g(i,j); end endendresult1=ifftshift(result1); %進(jìn)行反變換result2=ifftshift(result2); %進(jìn)行反變換X2=ifft2(result1);X3=uint8(real(X2));subplot(2,2,3);imshow(X3);title('butterworth濾波圖像');2.高通濾波clear;I1=imread(‘tire.tif’);subplot(2,2,1);imshow(I1);title('原始圖像');I2=imnoise(I1,'salt&pepper',0.02); %加入椒鹽噪聲subplot(2,2,2);imshow(I2);title('噪聲圖像');f=double(I2);g=fft2(f) %傅里葉變換g=fftshift(g) %轉(zhuǎn)換數(shù)據(jù)矩陣[N1,N2]=size(g);n=2;d0=50;n1=fix(N1/2);n2=fix(N2/2);fori=1:N1 %對(duì)頻域中循環(huán)濾波 forj=1:N2 d=sqrt((i-n1)^2+(j-n2)^2); %計(jì)算butterworth低通濾波轉(zhuǎn)換函數(shù) ifd==0h=0;elseh=1/(1+(d0/d)^(2*n));endresult1(i,j)=h*g(i,j); if(g(i,j)<50) result2(i,j)=0; else result2(i,j)=g(i,j); end endendresult1=ifftshift(result1); %進(jìn)行反變換result2=ifftshift(result2); %進(jìn)行反變換X2=ifft2(result1);X3=uint8(real(X2));subplot(2,2,3);imshow(X3);title('butterworth濾波圖像')X4=ifft2(result2);X5=uint8(real(X4));subplot(2,2,4)imshow(X5);title('理想高通濾波圖像');3．線(xiàn)性銳化濾波clear;I=imread(‘tire.tif’);h=fspecial('laplacian');I2=filter2(h,I);subplot(1,2,1);imshow(I);title('原始圖像');subplot(1,2,2);imshow(I2);title('濾波后的圖像');4．直方圖規(guī)定劃clear;I=imread(‘tire.tif’);hgram=0:255;J=histeq(I,hgram);subplot(2,2,1);imshow(I);title('原始圖像');subplot(2,2,2);imshow(J);title('直方圖規(guī)定化后圖像');subplot(2,2,3);imhist(I);title('原始圖像直方圖');subplot(2,2,4);imhist(J);title('規(guī)定化后圖像直方圖');5．直方圖均衡化clear;I=imread(‘tire.tif’);J=histeq(I);subplot(2,2,1);imshow(I);title('原始圖像');subplot(2,2,2);imshow(J);title('直方圖均衡化圖像');subplot(2,2,3);imhist(I);title('原始圖像直方圖');subplot(2,2,4);imhist(J);title('均衡化圖像直方圖');I=imread(‘tire.tif’);[J,T]=histeq(I);6．均值濾波clear;I=imread(‘tire.tif’);J=imnoise(I,'salt&pepper',0.02);subplot(2,2,1);imshow(J);title('噪聲圖像');K1=filter2(fspecial('average',3),J)/255;K2=filter2(fspecial('average',5),J)/255K3=filter2(fspecial('average',7),J)/255subplot(2,2,2);imshow(K1);title('3*3模板均值濾波');subplot(2,2,3);imshow(K2);title('5*5模板均值濾波');subplot(2,2,4);imshow(K3);title('7*7模板均值濾波');7．中值濾波clear;I=imread(‘tire.tif’);J=imnoise(I,'salt&pepper',0.02);subplot(2,2,1);imshow(J);title('噪聲圖像');K1=medfilt2(J,[3,3]);K2=medfilt2(J,[5,5]);K3=me