图像处理第四章

本文档详细介绍了使用MATLAB进行图像噪声模拟的实验,包括高斯噪声、均匀噪声、盐和胡椒噪声、对数正态分布、瑞利分布、指数分布和厄朗分布等不同类型噪声的生成,并通过直方图展示其分布。同时,实验还涉及正弦噪声的周期频谱和模式,以及不同滤波技术,如反调和滤波器、最大滤波器、最小滤波器、中值滤波器和维纳滤波器在椒盐噪声处理中的应用。此外,研究了Radon变换和FanBeam变换在图像重建中的效果,探讨了不同滤波器参数对重建质量的影响。
摘要由CSDN通过智能技术生成

4.1题目:
在这里插入图片描述
代码:clc
clear

r = imnoise2('gaussian',100000,1,0,1);

bins = 100;
hist(r,bins)
title('gaussian')

r = imnoise2('uniform',100000,1,0,1);
bins = 100;
figure,hist(r,bins)
title('uniform')

r = imnoise2('salt & pepper',1000,1,0.1,0.27);
bins = 100;
figure,hist(r,bins)
title('salt & pepper')

r = imnoise2('lognormal',100000,1);
bins = 100;
figure,hist(r,bins)
title('lognormal')

r = imnoise2('rayleigh',100000,1,0,1);
bins = 100;
figure,hist(r,bins)
title('rayleigh')

r = imnoise2('exponential',100000,1);
bins = 100;
figure,hist(r,bins)
title('exponential')

r = imnoise2('erlang',100000,1);
bins = 100;
figure,hist(r,bins)
title('erlang')
运行结果:
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
4.2题目:
在这里插入图片描述
代码:
clc
clear
C = [0 64; 0 128; 32 32; 64 0; 128 0; -32 32];
[r,R,S] = imnoise3(512, 512, C);
imshow(S,[])
title('[6个]指定冲击的正弦噪声周期频谱[1]')

imshow(r,[])
title('[6个]相应的正弦噪声周期模式[1]')
% S1 = fftshift(S);
% imshowMy(S1,[])
% figure,mesh(S)

C1 = C/2;
[r1,R1,S1] = imnoise3(512, 512, C1);
imshow(S1,[]),title('[6个]指定冲击的正弦噪声周期频谱[2]')
imshow(r1,[]),title('[6个]相应的正弦噪声周期模式[2]')

C2 = [6 32; -2 2];
[r2,R2,S2] = imnoise3(512, 512, C2);
imshow(S2,[]),title('[2个]指定冲击的正弦噪声周期频谱[3]')
imshow(r2,[]),title('[2个]相应的正弦噪声周期模式[3]')

A = [1 5];
[r3,R3,S3] = imnoise3(512, 512, C2, A);
imshow(S3,[]) %有两个不清楚的点,因为其振幅较小
title('[2个][使用非默认的不同振幅]指定冲击的正弦噪声周期频谱[4]')
imshow(r3,[])
title('[2个][使用非默认的不同振幅]相应的正弦噪声周期模式[4]')

subplot(3, 2, 1), imshow(S,[]),title('[6个]指定冲击的正弦噪声周期频谱[1]')
subplot(3, 2, 2), imshow(r,[]),title('[6个]相应的正弦噪声周期模式[1]')
subplot(3, 2, 3), imshow(S1,[]),title('[6个]指定冲击的正弦噪声周期频谱[2]')
subplot(3, 2, 4), imshow(r1,[]),title('[6个]相应的正弦噪声周期模式[2]')
subplot(3, 2, 5), imshow(r2,[]),title('[2个][使用非默认的不同振幅]指定冲击的正弦噪声周期频谱[4]')
subplot(3, 2, 6), imshow(r3,[]),title('[2个][使用非默认的不同振幅]相应的正弦噪声周期模式[4]');
运行结果:
在这里插入图片描述
4.3题目:
在这里插入图片描述
代码:
clc
clear

r = imnoise2('gaussian',100000,1,0,1);
bins = 100;
hist(r,bins)
title('gaussian')

r = imnoise2('uniform',100000,1,0,1);
bins = 100;
figure,hist(r,bins)
title('uniform')

r = imnoise2('salt & pepper',1000,1,0.1,0.27);
bins = 100;
figure,hist(r,bins)
title('salt & pepper&

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值