数字图像处理实验(8):PROJECT 04-04,Highpass Filtering Using a Lowpass Image
生活随笔
收集整理的這篇文章主要介紹了
数字图像处理实验(8):PROJECT 04-04,Highpass Filtering Using a Lowpass Image
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
實驗要求:
高通濾波器可以通過1減去低通濾波器的傳遞函數得到。
使用公式 計算可以的得到 。
實驗代碼:
% PROJECT 04-04 Highpass Filtering Using a Lowpass Image close all; clc; clear all;% img = imread('Fig4.11(a).jpg'); img = mat2gray(img); figure; subplot(2,2,1); imshow(img); title('原圖像');% 產生濾波函數 [M, N] = size(img); P = 2 * M; Q = 2 * N;alf = 50; H = zeros(P, Q); for i = 1:Pfor j = 1:QH(i, j) = exp(-((i-P/2)^2 + (j-Q/2)^2) / (2 * alf^2));end end% H = ones(P, Q); subplot(2,2,2); imshow(H); title('濾波函數');% % 圖像填充 [M, N] = size(img); P = 2 * M; Q = 2 * N;img_fp = zeros(P, Q); img_fp(1:M, 1:N) = img(1:M, 1:N);% [X, Y] = meshgrid(1:P, 1:Q); % ones = (-1)^(X+Y);% img_f = ones .* img_fp; img_f = zeros(P, Q); for x = 1:Pfor y = 1:Qimg_f(x, y) = img_fp(x, y) .* (-1)^(x+y);end endimg_F = fft2(img_f);img_G = img_F .* H; img_g = real(ifft2(img_G));% img_g = ones .* img_g;for x = 1:Pfor y = 1:Qimg_g(x, y) = img_g(x, y) .* (-1)^(x+y);end endimg_o = img_g(1:M, 1:N);subplot(2,2,3); imshow(img_o, []); title('高斯低通濾波后的圖像');% img_gmask = img_o - img; subplot(2,2,4); imshow(img_gmask, []); title('高斯高通濾波后的圖像');前面的代碼與此前的實驗相同,下面是最后的高通濾波處理。img_o 是前面得到的高斯低通濾波處理后的圖像,使用原圖像減去它,得到的就是高斯高通濾波后的圖像了。
img_gmask = img_o - img; subplot(2,2,4); imshow(img_gmask, []); title('高斯高通濾波后的圖像');實驗結果:
總結
以上是生活随笔為你收集整理的数字图像处理实验(8):PROJECT 04-04,Highpass Filtering Using a Lowpass Image的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 数字图像处理实验(7):PROJECT
- 下一篇: 数字图像处理实验(9):PROJECT