数字图像处理实验(8):PROJECT 04-04,Highpass Filtering Using a Lowpass Image 标签: 图像处理MATLAB 2017-05-25 0
2024-09-04 23:34:44
实验要求:
高通滤波器可以通过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('高斯高通滤波后的图像');
实验结果:
最新文章
- 机器学习理论知识部分--偏差方差平衡(bias-variance tradeoff)
- UITableView heightForHeaderInSection遇到的坑
- pandas应用之分组因子暴露和分位数分析
- HTML5+JS 《五子飞》游戏实现(八)人机对战
- [四校联考P3] 区间颜色众数 (主席树)
- 【Tree 1】树形结构数据呈现的递归算法实现
- [原创]-IIS7.5优化,支持同时10万个请求
- [算法] 快速排序 Quick Sort
- shortcut switch in terminal start pos & end pos
- Git分支(3/5) -- 禁用 Fast Forward 合并
- ORACLE里锁有以下几种模式,v$locked_object,locked_mode
- WebGL three.js学习笔记 使用粒子系统模拟时空隧道(虫洞)
- Python网络爬虫实战:根据天猫胸罩销售数据分析中国女性胸部大小分布
- PAT 1008 数组元素循环右移问题
- 永久关闭selinux
- KMP与AC自动机模板
- linux 升级python2.7
- 回归JavaScript基础(十)
- 2016-2017-2 20155309 南皓芯java第六周学习总结
- anaconda-ks.cfg详解
热门文章
- UIView+PYJExtension
- js 预解析
- 【java规则引擎】一个基于drools规则引擎实现的数学计算例子
- LOJ10034 图书管理
- 利用git bash和git gui向git远程仓库提交文件
- spring--Autowired setter 方法
- 三个Telnet链接类----来自网络
- [Client] looks like we got no XML document in....
- (转)Android中的页面切换动画
- Win10的Hosts文件修改后无法保存的问题解决方法,实测可以