2022.01.07，今天是服務外包競賽：隨便拿個獎隊的專案進行的第七天，今天根據專案要求繼續學習matlab數字影象處理

實驗七 影象增強—頻域濾波

一、 實驗目的

1．掌握怎樣利用傅立葉變換進行頻域濾波

2．掌握頻域濾波的概念及方法

3．熟練掌握頻域空間的各類濾波器

4．利用MATLAB程式進行頻域濾波

二、實驗原理

頻域濾波分為低通濾波和高通濾波兩類，對應的濾波器分別為低通濾波器和高通濾波器。頻域低通過濾的基本思想：

G(u,v)=F(u,v)H(u,v)

F(u,v)是需要鈍化影象的傅立葉變換形式，H(u,v)是選取的一個低通過濾器變換函式，G(u,v)是通過H(u,v)減少F(u,v)的高頻部分來得到的結果，運用傅立葉逆變換得到鈍化後的影象。

n階巴特沃茲低通濾波器(BLPF)(在距離原點處出現截至頻率)的傳遞函式為

與理想地通濾波器不同的是，巴特沃茲率通濾波器的傳遞函式並不是在處突然不連續。

高斯低通濾波器(GLPF)的傳遞函式為

其中，為標準差。

二、 實驗步驟：

1．調入並顯示所需的圖片；

2．利用MATLAB提供的低通濾波器實現影象訊號的濾波運算，並與空間濾波進行比較。

3．利用MATLAB提供的高通濾波器對影象進行處理。

4．記錄和整理實驗報告。

f1=imread('C:\Users\Gateway\Desktop\ͼƬ\camera.jpg');

f=rgb2gray(f1);

F=fft2(f);

S=fftshift(log(1+abs(F)));

figure;imshow(S,[]); %²úÉúƵÆ×ͼÏñ

h=fspecial('sobel');

figure;freqz2(h);

PQ=paddedsize(size(f));

H=freqz2(h,PQ(1),PQ(2));

H1=ifftshift(H);

figure;imshow(abs(H),[]);

Sobel邊緣提取

f=imread('C:\Users\Gateway\Desktop\ͼƬ\dabao.jpg');

f=rgb2gray(f);

f=im2double(f);

[VSFAT]=edge(f,'sobel','vertical');

figure;imshow(f);

figure;imshow(VSFAT);