数字图像处理实验(8):PROJECT 04-04,Highpass Filtering Using a Lowpass Image 标签: 图像处理MATLAB 2017-05-25 0
2017-05-25 09:48
856 查看
实验要求:
高通滤波器可以通过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:P for j = 1:Q H(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:P for y = 1:Q img_f(x, y) = img_fp(x, y) .* (-1)^(x+y); end end img_F = fft2(img_f); img_G = img_F .* H; img_g = real(ifft2(img_G)); % img_g = ones .* img_g; for x = 1:P for y = 1:Q img_g(x, y) = img_g(x, y) .* (-1)^(x+y); end end img_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('高斯高通滤波后的图像');
实验结果:
相关文章推荐
- 数字图像处理实验(9):PROJECT 04-05,Correlation in the Frequency Domain 标签: 图像处理MATLAB 2017-05-25 10:14
- 数字图像处理实验(4):PROJECT 02-04 [Multiple Uses],Zooming and Shrinking Images by Bilinear Interpolation 标签: 图像处理MATLAB
- 数字图像处理实验(7):PROJECT 04-03 , Lowpass Filtering 标签: 图像处理MATLAB 2017-05-25 09:30 109人
- 数字图像处理实验(17):PROJECT 06-04,Color Image Segmentation 标签: 图像处理MATLAB 2017-05-27 21:13
- 数字图像处理实验(10):PROJECT 05-01 [Multiple Uses],Noise Generators 标签: 图像处理MATLAB 2017-05-26 23:36
- 数字图像处理实验(6):PROJECT 04-02,Fourier Spectrum and Average Value 标签: 图像处理MATLABfft 2017-05-07 23:1
- 数字图像处理实验(13):PROJECT 05-04,Parametric Wiener Filter 标签: 图像处理MATLAB 2017-05-27 10:59
- 数字图像处理实验(1):PROJECT 02-01, Image Printing Program Based on Halftoning 标签: 图像处理MATLAB 2017-04-2
- 数字图像处理实验(17):PROJECT 06-04,Color Image Segmentation
- 数字图像处理实验(4):PROJECT 02-04 [Multiple Uses],Zooming and Shrinking Images by Bilinear Interpolation
- 数字图像处理实验(15):PROJECT 06-02,Pseudo-Color Image Processing 标签: 图像处理MATLAB 2017-05-27 20:53
- 数字图像处理实验(14):PROJECT 06-01,Web-Safe Colors 标签: 图像处理MATLAB 2017-05-27 20:45 116人阅读
- 数字图像处理实验(11):PROJECT 05-02,Noise Reduction Using a Median Filter 标签: 图像处理MATLAB 2017-05-26 23:
- 数字图像处理:基于MATLAB的车牌识别项目 标签: 图像处理matlab算法 2017-06-24 09:17 98人阅读 评论(0)
- 数字图像处理实验(总计23个)汇总 标签: 图像处理MATLAB 2017-05-31 10:30 175人阅读 评论(0)
- 数字图像处理实验(16):PROJECT 06-03,Color Image Enhancement by Histogram Processing 标签: 图像处理MATLAB 2017
- 数字图像处理实验(12):PROJECT 05-03,Periodic Noise Reduction Using a Notch Filter 标签: 图像处理MATLAB 2017-0
- 数字图像处理实验(2):PROJECT 02-02, Reducing the Number of Gray Levels in an Image 标签: 图像处理MATLAB 2017-
- 数字图像处理实验(12):PROJECT 05-03,Periodic Noise Reduction Using a Notch Filter
- 数字图像处理实验(5):Proj03-01 ~ Proj03-06 标签: 图像处理matlab 2017-04-30 10:39 184人阅读