用matlab制作数字图像高频噪声

高频噪声的制作思想可参考https://blog.csdn.net/qq_39148922/article/details/86609629

本文介绍数字图像高频噪声的制作,原始噪声用正态分布随机数生成

源代码如下

%数字图像高频噪声的制作

%运行参数
mLength=800;%噪声图像的宽度
nLength=600;%噪声图像的高度
ampi=0.46;%噪声的幅度
freqCut=0.34*sqrt(mLength^2+nLength^2);%高频噪声的截止频率
filterStep=3;%滤波器的阶数

%生成图像的中心距离矩阵
[u,v]=dftuv(mLength,nLength);
%滤波器的频域函数
Huv=1./(1+(freqCut^2./(u.^2+v.^2)).^filterStep);
%白噪声
whiteNoise=sqrt(ampi)*randn(mLength,nLength);
%噪声通过滤波器,得到高频噪声
highFreqNoise=dftfilt(whiteNoise,Huv);

%结果展示
figure
imshow(whiteNoise+0.5);
title('原始噪声')

figure
imshow(highFreqNoise+0.5)
title('高频噪声')

 

[u,v]=dftuv(mLength,nLength);

dftuv这个函数如果你的matlab工作环境中没有,可以从这里获得它的源代码:

https://blog.csdn.net/qq_39148922/article/details/86625377

在当前的工作目录里面新建一个函数,粘贴进去即可使用

 

这个函数的返回值[u,v]是为了滤波器频域函数的表达方便

u的图形是这样:

v的图形是这样:

 

从图形中也可以看出,它们分别表示横纵方向上矩阵到中心的距离

这样就可以利用两个距离矩阵表达出所需的滤波器频域函数

Huv=1./(1+(freqCut^2./(u.^2+v.^2)).^filterStep);

接下来的函数

highFreqNoise=dftfilt(whiteNoise,Huv);

函数dftfilt的行为是将参数1进行二维傅里叶变换得到频谱,再将得到的频谱和参数2相乘,相乘后进行逆傅里叶变换得到返回值

函数的源代码可从这里获得:https://blog.csdn.net/qq_39148922/article/details/86625526

至此即可获得高频噪声

 

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值