clc,clear,close all % 清理命令区、清理工作区、关闭显示图形
warning off % 消除警告
feature jit off % 加速代码运行
D0 = 20; % 阻止的频率点与频域中心的距离
n = 2; % 阶次
im = imread(‘coloredChips.png’); % 原图像
R = imnoise(im(:,:,1),‘gaussian’,0,0.01); % R + 白噪声
G = imnoise(im(:,:,2),‘gaussian’,0,0.01); % G + 白噪声
B = imnoise(im(:,:,3),‘gaussian’,0,0.01); % B + 白噪声
im = cat(3,R,G,B); % 原图像 + 白噪声
R1 = freqfilter_btw_lp(R,D0,n); % 巴特沃斯低通滤波器
G1 = freqfilter_btw_lp(G,D0,n); % 巴特沃斯低通滤波器
B1 = freqfilter_btw_lp(B,D0,n); % 巴特沃斯低通滤波器
im1 = cat(3,R1,G1,B1);
figure(‘color’,[1,1,1])
subplot(121),imshow(im,[]); title(‘原始图像’)
subplot(122),imshow(im1,[]); title(‘巴特沃斯低通滤波图像’);
function im5 = freqfilter_btw_lp(im,D0,n)
if ~isa(im,‘double’)
im1 = double(im)/255;
end
im2 = fft2(im1); % 傅里叶变换
im3 = fftshift(im2); % 中心化

[N1, N2] = size(im3);
n1 = fix(N1 / 2);
n2 = fix(N2 / 2);
for i = 1:N1
for j = 2:N2
d = sqrt((i-n1)2+(j-n2)2);
h = 1/(1 + 0.414 * (d / D0)^(2*n)); % 巴特沃斯低通滤波器
result(i,j) = h * im3(i,j);
end
end
result = ifftshift(result); % 反中心化
im4 = ifft2(result); % 反变换
im5 = im2uint8(real(im4)); % 滤波图像

end

MATLAB巴特沃斯低通滤波图像相关推荐

  1. 用matlab编程实现数字图像理想低通滤波、高斯低通滤波和巴特沃斯低通滤波去噪算法

    1 理想低通滤波 %理想低通 I = imread('fig.png'); I=rgb2gray(I); figure(1); subplot(221),imshow(I); title('原图像') ...

  2. c++ opencv数字图像处理:频率域滤波--低通滤波--巴特沃斯低通滤波

    文章目录 前言 一.巴特沃斯低通滤波器(BLPF) 二.代码 三.说明 前言 数字图像处理c++ opencv(VS2019 opencv4.53)持续更新 一.巴特沃斯低通滤波器(BLPF) D2( ...

  3. MATLAB巴特沃斯陷波滤波图像

    clc,clear,close all % 清理命令区.清理工作区.关闭显示图形 warning off % 消除警告 feature jit off % 加速代码运行 D0 = 4; % 阻止的频率 ...

  4. python理想低通滤波、巴特沃斯低通滤波、高斯低通滤波实现

    代码 代码如下(示例): import numpy as np import cv2 as cv image = cv.imread('2.PNG') # print(image.shape) ima ...

  5. PX4中的二阶巴特沃斯低通滤波

    源码 const float fr = sample_freq/_cuttof_freq; const float ohm = tanf(M_PI_F, fr); const float c = 1. ...

  6. Matlab实现 理想低通、巴特沃斯低通、高斯低通、理想高通、巴特沃斯高通、高斯高通(d=10,50,150)

    数字图像处理第二次编程课后作业 理想低通(d=10,50,150): close all; clear all;%% ---------Ideal Lowpass Filters (Fre. Doma ...

  7. matlab常见的图像增强技术(包括基于幂次变换,对图像进行均衡化处理,巴特沃斯低通,理想低通,梯形低通滤波, 均值滤波,中值滤波,最大,最小值滤波,修正后的阿尔法滤波器)

    1.基于幂次变换中的r值,比较不同r 值下图像增强的效果 代码 : I = imread('D:\图片\TH.JFIF');subplot (1,4,1);imshow(I);title('原始图像' ...

  8. matlab 对图像进行低通滤波,matlab图像处理实现低通滤波

    matlab图像处理实现低通滤波 matlab图像处理实现低通滤波 clc; clear all; img=imread('lena.jpg'); img_noise = imnoise(img, ' ...

  9. MATLAB高斯低通滤波图像

    clc,clear,close all % 清理命令区.清理工作区.关闭显示图形 warning off % 消除警告 feature jit off % 加速代码运行 D0 = 20; % 阻止的频 ...

最新文章

  1. PTA 基础程序编程集 7-2 然后是几点 C语言
  2. linux七个运行级如何设置,RHEL7 运行级别简介及切换操作
  3. 设计模式在Netty中的应用-装饰者模式源码举例
  4. 关于SAP CRM One Order状态(Status)和Status Profile的处理逻辑
  5. latex 分页_latex 图片跨页显示问题???
  6. C# DllImport用法和路径问题
  7. 征稿延期 | 2019亚洲语言处理国际大会(IALP2019)征稿延期
  8. python 反编译exe
  9. 新疆公需课继续教育答案 自动获取
  10. STM32F407VET6+cubemx+FSMC+ST7789
  11. 盗心贼歌曲用计算机多少数字,抖音上常见背景音乐歌词盗心的贼是那首歌?
  12. Redis技能—底层IO多路复用
  13. 采集百度搜寻结果,图片不显示的解决方法
  14. 手机时钟软件推荐,创意时钟APP介绍
  15. JAVA知识两百问(181~200)
  16. google搜索特殊字符
  17. EXCEL VBA 备忘录
  18. 【JDBC】idea添加mysql-jar包(很轻松)
  19. Springboot jar包外置教程
  20. linux学习课程从入门到精通:文件的归档和压缩

热门文章

  1. [logback]去掉kafka日志与additivity
  2. 原来这才是使用普洱熟茶减肥的正确方式!
  3. 一个普通程序员的副业历程
  4. 有哪些半入耳蓝牙耳机推荐,盘点学生党平价蓝牙耳机
  5. 心火肝火胃火肺火的症状区别及治疗方法
  6. 小程序api 分享scene_网课查题题库接口API-在线免费授权,微信小程序网课答案api接口...
  7. Java Student类
  8. 操作系统学习福利:600 条最强 Linux 命令总结
  9. Windows使用任务计划程序执行定时任务,开机自启任务,登录自启任务
  10. Delphi常用属性