基于matlab圖像頻域分析與仿真-資料下載頁
2025-11-03 15:27本頁面
【導(dǎo)讀】實驗證明，低通濾波器和高通濾波器都能較好的改善圖像質(zhì)量。X射線穿過身體并由位于對面環(huán)中的相應(yīng)檢測器收集起來。當物體沿垂直于檢測器的方向運動時就產(chǎn)生一系列的切片，這些切片組成了物體內(nèi)部的再現(xiàn)圖像。許多能獲得三維圖像的設(shè)備和分析處理三維圖像的系統(tǒng)已經(jīng)研制成功了，圖像處理技術(shù)得到了廣泛的應(yīng)用。地理學(xué)用相同或相似的技術(shù)從航空和衛(wèi)星圖像中研究污染模式。在文章中Kim提出了保持亮度特性的直方圖均衡算法。Kim的改進算法提出后，引起了許多學(xué)者的關(guān)注。
【正文】 fftshift(F1)。 subplot(2,2,1),imshow(F)。 subplot(2,2,2),imshow(F2)。 subplot(2,2,3),imshow(log(1+abs(F3)),[])。 2）圖像的快速傅立葉變換，僅對相位部分做傅立葉反變換 F=imread( 39。E:\應(yīng)用軟件\matlab\bin\139。,39。jpeg39。)。 F1=fft2(F)。 F2=log(1+abs(F1))。 F3=fftshift(F1)。 subplot(2,2,1),imshow(F)。 subplot(2,2,2),imshow(F2)。 subplot(2,2,3),imshow(log(1+abs(F3)),[])。 3）低通濾波器的應(yīng)用實例： h=imread(39。E:\應(yīng)用軟件\matlab\bin\139。,39。jpeg39。)。 i=rgb2gray(h)。 subplot(2,2,1)。imshow(i)。 f=double(i)。 g=fft2(f)。 g=fftshift(g)。 [M,N]=size(g)。 d0=15。 m=fix(M/2)。n=fix(N/2)。 for i=1:M。 for j=1:N。 d=sqrt((im)^2+(jn)^2)。 if(d=d0)。 h=1。 else h=0。 end result(i,j)=h*g(i,j)。 end end result=ifftshift(result)。 j1=ifft2(result)。 j2=uint8(real(j1))。 subplot(2,2,2)。imshow(j2)。 d0=30。 m=fix(M/2)。n=fix(N/2)。 for i=1:M。 for j=1:N。 d=sqrt((im)^2+(jn)^2)。 if(d=d0)。 h=1。 else h=0。 end result(i,j)=h*g(i,j)。 end end result=ifftshift(result)。 j1=ifft2(result)。 j2=uint8(real(j1))。 subplot(2,2,3)。imshow(j2)。 d0=100。 m=fix(M/2)。n=fix(N/2)。 for i=1:M。 for j=1:N。 d=sqrt((im)^2+(jn)^2)。 if(d=d0)。 h=1。 else h=0。 end result(i,j)=h*g(i,j)。 end end result=ifftshift(result)。 j1=ifft2(result)。 j2=uint8(real(j1))。 subplot(2,2,4)。imshow(j2)。4） 巴特沃斯高通濾波器應(yīng)用實例： i=imread(39。E:\應(yīng)用軟件\matlab\bin\139。,39。jpeg39。)。 i=rgb2gray(i)。 subplot(2,2,1)。imshow(i)。 f=double(i)。 g=fft2(f)。 g=fftshift(g)。 [M,N]=size(g)。 nn=2。 d0=15。 m=fix(M/2)。n=fix(N/2)。 for i=1:M for j=1:N d=sqrt((im)^2+(jn)^2)。 if(d==0) h=0。 else h=1/(1+*(d0/d)^(2*nn))。 end result(i,j)=h*g(i,j)。 end end result=ifftshift(result)。 j1=ifft2(result)。 j2=uint8(real(j1))。 subplot(2,2,2)。imshow(j2)。 d0=25。 m=fix(M/2)。n=fix(N/2)。 for i=1:M for j=1:N d=sqrt((im)^2+(jn)^2)。 if(d==0) h=0。 else h=1/(1+*(d0/d)^(2*nn))。 end result(i,j)=h*g(i,j)。 end end result=ifftshift(result)。 j1=ifft2(result)。 j2=uint8(real(j1))。 subplot(2,2,3)。imshow(j2)。 d0=80。 m=fix(M/2)。n=fix(N/2)。 for i=1:M for j=1:N d=sqrt((im)^2+(jn)^2)。 if(d==0) h=0。 else h=1/(1+*(d0/d)^(2*nn))。 end result(i,j)=h*g(i,j)。 end end result=ifftshift(result)。 j1=ifft2(result)。 j2=uint8(real(j1))。 subplot(2,2,4)。imshow(j2)。21
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1