matlab 数字图像滤波,数字图像处理-Matlab实现图像的滤波

最新推荐文章于 2024-01-01 12:31:58 发布
最新推荐文章于 2024-01-01 12:31:58 发布
阅读量445 收藏 1
点赞数
文章标签： matlab 数字图像滤波

目录

1.空域高斯滤波器

2.中值滤波

3.频率域高斯滤波

4.滤波器测试脚本

5.测试结果

1.空域高斯滤波器

function [OutImage,ConvAarray] = FilterGauss(InputImage,Size)

%================================================================

% 功能：单通道图像高斯滤波器

% 参数：InputImage为输入单通道图像，Size为模板尺寸

% 返回值：OutImage为InputImage同维数组

% 主要思路：利用卷积核，模板对称，利用相关等于卷积

% 备注：如果图像为多通道则需要重复调用

% 调用方法：OutImage = FilterGauss(InputImage)

% 日期：2019.12.02

% 作者：Leetion

[iLimit,jLimit] = size(InputImage);

OutImage = InputImage;

switch (Size)

case 3

ConvAarray =double((1/16)*[1 2 1;2 4 2;1 2 1]);

for yIndex = 2:iLimit-1

for xIndex = 2:jLimit-1

NeighbourArray = [InputImage(yIndex-1,xIndex-1) InputImage(yIndex-1,xIndex) InputImage(yIndex-1,xIndex+1);...

InputImage(yIndex,xIndex-1) InputImage(yIndex,xIndex) InputImage(yIndex,xIndex+1);...

InputImage(yIndex+1,xIndex-1) InputImage(yIndex+1,xIndex) InputImage(yIndex+1,xIndex+1)];

OutImage(yIndex,xIndex) = uint8(sum(sum(ConvAarray.*double(NeighbourArray))));

end

end

case 5

ConvAarray =double((1/90)*[1 2 4 2 1;2 4 8 4 2;4 8 16 8 4;2 4 8 4 2; 1 2 4 2 1]);

for yIndex = 3:iLimit-2

for xIndex = 3:jLimit-2

NeighbourArray = [InputImage(yIndex-2,xIndex-2) InputImage(yIndex-2,xIndex-1) InputImage(yIndex-2,xIndex) InputImage(yIndex-2,xIndex+1) InputImage(yIndex-2,xIndex+2);...

InputImage(yIndex-1,xIndex-2) InputImage(yIndex-1,xIndex-1) InputImage(yIndex-1,xIndex) InputImage(yIndex-1,xIndex+1) InputImage(yIndex-1,xIndex+2);...

InputImage(yIndex,xIndex-2) InputImage(yIndex,xIndex-1) InputImage(yIndex,xIndex) InputImage(yIndex,xIndex+1) InputImage(yIndex,xIndex+2);...

InputImage(yIndex+1,xIndex-2) InputImage(yIndex+1,xIndex-1) InputImage(yIndex+1,xIndex) InputImage(yIndex+1,xIndex+1) InputImage(yIndex+1,xIndex+2);...

InputImage(yIndex+2,xIndex-2) InputImage(yIndex+2,xIndex-1) InputImage(yIndex+2,xIndex) InputImage(yIndex+2,xIndex+1) InputImage(yIndex+2,xIndex+2)];

OutImage(yIndex,xIndex) = uint8(sum(sum(ConvAarray.*double(NeighbourArray))));

end

end

otherwise

warning("只支持三阶和五阶");

end

end

2.中值滤波

function OutImage = FilterMiddle(InputImage,Size)

%================================================================

% 功能：单通道图像中值滤波器

% 参数：InputImage为输入单通道图像，Size为邻域尺寸

% 返回值：OutImage为InputImage同维数组

% 主要思路：求解邻域中值

% 备注：如果图像为多通道则需要重复调用

% 调用方法：OutImage = FilterMiddle(InputImage)

% 日期：2019.12.02

% 作者：Leetion

[iLimit,jLimit] = size(InputImage);

OutImage = InputImage;

if Size(1)==1 && Size(2)==3

flag = 1;

elseif Size(1)==1 && Size(2)==5

flag = 2;

elseif Size(1)==3 && Size(2)==1

flag = 3;

elseif Size(1)==5 && Size(2)==1

flag = 4;

elseif Size(1)==3 && Size(2)==3

flag = 5;

elseif Size(1)==5 && Size(2)==5

flag = 6;

else

warning("不支持的尺寸")

end

switch(flag)

case 1

for yIndex = 1:iLimit

for xIndex = 2:jLimit-1

NeighbourArray = [InputImage(yIndex,xIndex-1) InputImage(yIndex,xIndex) InputImage(yIndex,xIndex+1)];

OutImage(yIndex,xIndex) = median(NeighbourArray);

end

end

case 2

for yIndex = 1:iLimit

for xIndex = 3:jLimit-2

NeighbourArray = [InputImage(yIndex,xIndex-2) InputImage(yIndex,xIndex-1) InputImage(yIndex,xIndex) InputImage(yIndex,xIndex+1) InputImage(yIndex,xIndex+2)];

OutImage(yIndex,xIndex) = median(NeighbourArray);

end

end

case 3

for yIndex = 2:iLimit-1

for xIndex = 1:jLimit

NeighbourArray = [InputImage(yIndex-1,xIndex),InputImage(yIndex,xIndex),InputImage(yIndex+1,xIndex)];

OutImage(yIndex,xIndex) = median(NeighbourArray);

end

end

case 4

for yIndex = 3:iLimit-2

for xIndex = 1:jLimit

NeighbourArray = [InputImage(yIndex-2,xIndex) InputImage(yIndex-1,xIndex),InputImage(yIndex,xIndex),InputImage(yIndex+1,xIndex) InputImage(yIndex+2,xIndex)];

OutImage(yIndex,xIndex) = median(NeighbourArray);

end

end

case 5

for yIndex = 2:iLimit-1

for xIndex = 2:jLimit-1

NeighbourArray = [InputImage(yIndex-1,xIndex-1) InputImage(yIndex-1,xIndex) InputImage(yIndex-1,xIndex+1)...

InputImage(yIndex,xIndex-1) InputImage(yIndex,xIndex) InputImage(yIndex,xIndex+1)...

InputImage(yIndex+1,xIndex-1) InputImage(yIndex+1,xIndex) InputImage(yIndex+1,xIndex+1)];

OutImage(yIndex,xIndex) = median(NeighbourArray);

end

end

case 6

for yIndex = 3:iLimit-2

for xIndex = 3:jLimit-2

NeighbourArray = [InputImage(yIndex-2,xIndex-2) InputImage(yIndex-2,xIndex-1) InputImage(yIndex-2,xIndex) InputImage(yIndex-2,xIndex+1) InputImage(yIndex-2,xIndex+2)...

InputImage(yIndex-1,xIndex-2) InputImage(yIndex-1,xIndex-1) InputImage(yIndex-1,xIndex) InputImage(yIndex-1,xIndex+1) InputImage(yIndex-1,xIndex+2)...

InputImage(yIndex,xIndex-2) InputImage(yIndex,xIndex-1) InputImage(yIndex,xIndex) InputImage(yIndex,xIndex+1) InputImage(yIndex,xIndex+2)...

InputImage(yIndex+1,xIndex-2) InputImage(yIndex+1,xIndex-1) InputImage(yIndex+1,xIndex) InputImage(yIndex+1,xIndex+1) InputImage(yIndex+1,xIndex+2)...

InputImage(yIndex+2,xIndex-2) InputImage(yIndex+2,xIndex-1) InputImage(yIndex+2,xIndex) InputImage(yIndex+2,xIndex+1) InputImage(yIndex+2,xIndex+2)];

OutImage(yIndex,xIndex) = median(NeighbourArray);

end

end

end

end

3.频率域高斯滤波

function OutImage = FreFilterGauss(InputImage,D0)

%================================================================

% 功能：频域高斯滤波

% 参数：InputImage为输入单通道图像,D0为截止频率

% 返回值：OutImage为输出图像

% 主要思路：FFT转换到频率域，频率域乘积高斯传函，再IFFT2到时域

% 备注：使用库函数计算FFT和IFFT

% 调用方法：OutImage = FreFilterGauss(InputImage,D0)

% 日期：2019.12.02

% 作者：Leetion

[iLimit,jLimit] = size(InputImage);

FreInputImage = fft2(InputImage);

for vIndex = 1:iLimit

for uIndex = 1:jLimit

Distance = sqrt((vIndex-(1+iLimit)/2)^2+(uIndex-(1+jLimit)/2)^2);

H = exp(-(Distance^2/(2*D0^2)));

FreOutImage(vIndex,uIndex) = H*FreInputImage(vIndex,uIndex);

end

end

OutImage = uint8(ifft2(FreOutImage));

end

4.滤波器测试脚本

%================================================================

% 功能：滤波器测试

% 主要思路：

% 备注：

% 日期：2019.11.29

% 作者：李习银

clc,clear

close all

InputImage = imread('eight.tif');

figure(1)

subplot(2,3,1);

imshow(InputImage);

title("原图");

InputImage = imnoise(InputImage,'gaussian',0.001,0.001);

% InputImage = imnoise(InputImage,'salt & pepper');

subplot(2,3,2);

imshow(InputImage);

title("Gaussian Noise");

OutImage = FilterGauss(InputImage,3);

subplot(2,3,3);

imshow(OutImage);

title("Gaussian Filter");

subplot(2,3,4);

OutImage = FilterMiddle(InputImage,[5,5]);

imshow(OutImage);

title("Middle Filter");

subplot(2,3,5);

OutImage = FreFilterGauss(InputImage,1500);

imshow(OutImage);

title("FreGauss Filter");

figure(2)

subplot(2,3,1);

imshow(InputImage);

title("原图");

InputImage = imnoise(InputImage,'salt & pepper');

subplot(2,3,2);

imshow(InputImage);

title("Gaussian Noise");

OutImage = FilterGauss(InputImage,3);

subplot(2,3,3);

imshow(OutImage);

title("Gaussian Filter");

subplot(2,3,4);

OutImage = FilterMiddle(InputImage,[5,5]);

imshow(OutImage);

title("Middle Filter");

subplot(2,3,5);

OutImage = FreFilterGauss(InputImage,1500);

imshow(OutImage);

title("FreGauss Filter");

5.测试结果

decf1375789b53a584ef16443b4bae3b.png 高斯白噪声

2bf74845a6a48ba18d360fca3eb66d32.png 椒盐噪声

祝昊
关注
MATLAB实现图像滤波及噪声消除【图像处理实战】
04-23
项目详情请参见：https://handsome-man.blog.csdn.net/article/details/130327927 本资源使用MATLAB实现对添加了噪声的图像，使用滤波器对图像进行平滑处理，实现图像滤波及噪声消除。 项目代码可顺利编译运行~
MATLAB--数字图像处理--图像的几何变换
小席的学习日志
04-11 5649
图像的几何变换 图像的几何变换是将一幅图像中的坐标映射到另外一幅图像中的新坐标位置，它不改变图像的像素值，只是改变像素所在的几何位置，使原始图像按照需要产生位置、形状和大小的变换。 1、图像的平移 2、图像的镜像 3、图像的缩放 4、图像的转置 5、图像的旋转 6、图像的剪切 7、图像的空间变换 ...
MATLAB--数字图像处理 图像噪声与滤波处理
热门推荐
海到无边天作岸 山登绝顶我为峰
12-12 1万+
一、实验名称 图像的噪声与滤波处理 二、实验目的 1.熟悉MATLAB软件的使用。 2.掌握图像的噪声与滤波处理。 三、实验内容 1.对一张图片添加不同强度的高斯噪声、椒盐噪声 2.对图像进行不同模板的均值滤波、高斯加权滤波、中值滤波，对比结果并分析不同滤波方式的优劣 四、实验仪器与设备 Win10 64位电脑 MATLAB2017a 五、实验原理 高斯噪声    &...
数字图像处理-Matlab实现图像点运算
LLeetion的博客
12-04 4027
目录 图像反色 图像线性点运算 图像对数点运算 图像幂数点运算 图像灰度均衡化 以上函数测试脚本 1.线性，反色，对数，幂数 2.灰度均衡化 结果 图像反色 图像反色利用的公式为。Matlab本擅长矩阵运算可以直接运算，由于实验要求不能使用矩阵运算，使用循环完成，以下其他点运算也是一样的处理，不再解释。 function OutImage = InverseImageCo...
数字图像处理-Matlab实现图像的几何变换
LLeetion的博客
12-04 3492
目录 图像的平移 图像的旋转 图像的切变 图像的镜像 上述函数的测试脚本 结果 图像的平移 function OutImage = ImageShiftFun(InputImage,deltaX,deltaY,method) %================================================================ % 功能：图像的x...
Matlab实现数字图像处理——滤波
weixin_44616879的博客
05-29 1万+
Matlab实现数字图像处理——滤波同态滤波局部直方图均衡傅立叶变换自适应直方图均衡巴特沃斯滤波器均值滤波高斯滤波器半邻域滤波中值滤波多窗口中值滤波百分比滤波锐化拉普拉斯边界检测与增强噪声对拉普拉斯边界检测的影响 本文介绍了使用matlab代码实现数字图像处理中各种滤波方法。 同态滤波 I = imread('FeiWang.jpg'); G = rgb2gray(I); figure(); imshow(G), title('OriginalFigure'); homo_flt(G); function
数字图像均值滤波和高斯滤波分析_高斯滤波_matlab_均值滤波_数字图像处理_
09-30
总之，均值滤波和高斯滤波是数字图像处理中重要的噪声抑制工具，MATLAB提供便捷的接口实现这两种滤波。通过对这两种滤波方法的理解和实践，我们可以更好地掌握图像处理技术，并在处理实际问题时做出明智的选择。
matlab数字图像处理课程-频率域滤波matlab-源码
10-01
在本课程资源中，主题聚焦于使用MATLAB进行数字图像处理，特别是频率域滤波的实践。MATLAB是一款强大的数学计算软件，广泛应用于科学计算、数据分析以及图像处理等领域。在这个课程中，你将学习如何利用MATLAB的强大...
图像空间域滤波降噪_数字图像处理 matlab
07-15
数字图像处理课程的基础实验内容，图像空间域滤波降噪
图像处理的常见滤波方法，matlab代码
09-17
包含图像处理中常用的几种滤波方法，有matlab代码和ppt讲解，包含均值滤波，中值滤波，高斯滤波，双边滤波，有示例演示。
MATLAB 图像滤波器
tutouxiaotu的博客
09-26 8511
是一种典型的非线性滤波，是基于排序统计理论的一种能够有效抑制噪声的非线性信号处理技术，基本思想是用像素点邻域灰度值的中值来代替该像素点的灰度值，让周围的像素值接近真实的值从而消除孤立的噪声点。2．接着学习了三种不同的算子对图像的锐化处理，切身感受了不同的算子对图像锐化的侧重点不同，比如Sobel算子能够锐化的边缘信息较强，在今后应用中如果想提取图像的边缘信息最好选用Sobel算子。椒盐噪声又称脉冲噪声，它随机改变一些像素值，是由图像传感器，传输信道，解码处理等产生的黑白相间的亮暗点噪声。
MATLAB--数字图像处理 均值滤波
海到无边天作岸 山登绝顶我为峰
09-16 1万+
均值滤波 代码： t=imread('a1.jpg'); [m,n,z]=size(t); t1=imnoise(t,'gaussian',0,0.02); imshow(t1),title('加入高斯噪声后') t2=t; for i=1:m-2 for j=1:n-2 t2(i+1,j+1,:)=round(t2(i,j,:)/9+t2(i,j+1,:)/9+t2(i,j+2,:)/9+t2...
数字图像处理 中值滤波 MATLAB实验
Jasmine-Lily的博客
06-10 5725
一、原理_均值滤波 中值滤波的基本思想是将图像中每个像素的灰度值用其邻域内像素灰度的中值代替，它是一种非线性平滑滤波算法。 设加噪图像为 f(x,y) ，经中值滤波处理后的图像为g(x,y) ，则： 式中，S是(x,y)像素点的邻域。本实验分别选用3×3、5×5、7×7的中值滤波窗口对图像进行处理。 需要注意的是，当模板滑动到图像边缘时，模板的部分行或列就会处于图像之外，本实验可采用下面的任一种方法处理这种问题： ①先在图像边缘以外补上一行和一列灰度为零的像素点，再对原图像中的每个像素点进行模板处理。 ②
matlab的图像滤波,Matlab实现图像滤波
weixin_29062823的博客
03-24 2892
Matlab实现图像滤波一、实验目的1．实现图片添加高斯白噪声,椒盐噪声,泊松噪声,乘性噪声等并进行比较；2．实现图片进行均值滤波，中值滤波，高斯滤波，拉普拉斯滤波等处理；二、实验步骤1.读取图片image=imread(‘test2.bmp’); %读取图片subplot(5,2,1); %指定显示图片位置为5乘2矩阵第1个位置imshow(image),title(‘原图’); %显示所读取的...
数字图像处理---低高通滤波实验（MATLAB实现）
OpenCV学习笔记
04-29 1万+
数字图像处理MATLAB实现高低通滤波器代码以及实验原理和实验效果图
基于Matlab的各种图像滤波Filter算法（代码开源）
最新发布
black_sneak的博客
01-01 6128
基于Matlab的各种图像滤波算法，图像滤波算法包含：均值滤波、中值滤波、高斯滤波、双边滤波、引导滤波。滤波算法可调用API函数，方便快捷！
MATLAB图像的频域低通滤波（灰度图像滤波+彩色图像滤波）
m0_60004633的博客
03-29 3560
MATLAB图像的频域低通滤波（灰度图像滤波+彩色图像滤波）
基于Matlab的图像加噪滤波处理和图像边缘检测
qq_45156707的博客
10-27 3104
对一幅原始图像，加上高斯噪声，实现均值滤波，高斯滤波及中值滤波，比较滤波结果；对图像使用Sobel及Canny算子检测，比较边缘检测结果。
数字图像处理实验：空域滤波与MATLAB实现
"本次实验是关于数字图像处理的空域滤波，主要涉及图像的均值滤波和中值滤波。实验使用MATLAB作为工具，通过读取图像、添加高斯噪声和椒盐噪声，然后应用不同大小的模板进行滤波处理，对比不同滤波方法的效果。" 在...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值