数字图像处理实验六图像复原

数字图像处理学习记录
实验内容
（1）基于均值滤波器的图像复原。

img=imread('cameraman.tif');%读取图片
subplot(2,3,1);
imshow(img); %显示图像
title(''),xlabel('原图');img_noise=double(imnoise(img,'gaussian',0.06)); %对图像附加高斯噪声
subplot(2,3,2);
imshow(img_noise,[]); %显示加噪图像
title(''),xlabel('高斯噪声污染的图像');img_mean=imfilter(img_noise,fspecial('average',3));%对附加有高斯噪声的图像实行均值滤波器
subplot(2,3,3);
imshow(img_mean,[]);%显示算术均值滤波后的图像
title(''),xlabel('用3x3算数均值滤波器滤波');img_mean=exp(imfilter(log(img_noise),fspecial('average',3)));%对附加有高斯噪声的图像实施几何均值滤波
subplot(2,3,4);
imshow(img_mean,[]);
title(''),xlabel('3x3的几何均值滤波器滤波');Q=-1.5;  %对高斯噪声图像进行Q取负值的逆谐波滤波
img_mean=imfilter(img_noise.^(Q+1),fspecial('average',3))./imfilter(img_noise.^Q,fspecial(('average'),3));
subplot(2,3,5);
imshow(img_mean,[]);
title(''),xlabel('Q=-1.5的逆谐波滤波器滤波');Q=1.5;  %对高斯噪声图像进行Q取正数的逆谐波滤波
img_mean=imfilter(img_noise.^(Q+1),fspecial('average',3))./imfilter(img_noise.^Q,fspecial(('average'),3));
subplot(2,3,6);
imshow(img_mean,[]);
title(''),xlabel('Q=1.5的逆谐波滤波器滤波');

（2）基于顺序统计滤波器的图像复原。

img=imread('cameraman.tif');%读取图片
subplot(2,3,1);
imshow(img); %显示图像
title(''),xlabel('原图');img_noise=double(imnoise(img,'salt & pepper',0.06));%加“椒盐”噪声
subplot(2,3,2);
imshow(img_noise,[]); %显示图像
title(''),xlabel('椒盐噪声污染的图像');img_mean=imfilter(img_noise,fspecial('average',5));%对附加“椒盐”噪声的图像实行算术均值滤波
subplot(2,3,3);
imshow(img_mean,[]); %显示图像
title(''),xlabel('均值滤波结果');img_median=medfilt2(img_noise);%对附加“椒盐”噪声的图像实行中值滤波
subplot(2,3,4);
imshow(img_median,[]); %显示图像
title(''),xlabel('中值滤波结果');img_median2=medfilt2(img_noise);
subplot(2,3,5);
imshow(img_median2,[]); %显示图像
title(''),xlabel('再次中值滤波结果');

（3）基于带阻滤波器的图像复原。

img=imread('lenazhouqizaosheng.jpg');%读取图片
[M,N]=size(img);%得到图像的高度和宽度
subplot(2,2,1);
imshow(img); %显示图像
title(''),xlabel('含有周期噪声的原图');IF=fftshift(fft2(img));  %对加噪图像进行傅里叶变换，并将原点移至中心点
IFV=log(1+abs(IF));  %加噪图像的频谱
subplot(2,2,2);
imshow(IFV,[]); %显示加噪图像的频谱
title(''),xlabel('加噪图像的频谱');freg=50;  %设置带阻滤波器的中心频率
width=5;  %设置带阻滤波器的频带宽度
ff=ones(M,N);
for i=1:Mfor j=1:Nff(i,j)=1-exp(-0.5*((((i-M/2)^2+(j-N/2)^2)-freg^2)/(sqrt((i-M/2)^2+(j-N/2)^2)*width))^2);%高斯带阻滤波器end
end
subplot(2,2,3);
imshow(ff,[]); %显示高斯带阻滤波器
title(''),xlabel('高斯带阻滤波器');out=IF.*ff;    %矩阵点乘实现频域滤波
out=ifftshift(out);  %原点移回左上角
out=ifft2(out);  %傅里叶变换
out=abs(out);    %取绝对值
out=out/max(out(:));  %归一化
subplot(2,2,4);
imshow(out,[]); %显示滤波结果
title(''),xlabel('滤波效果图');

数字图像处理实验六图像复原相关推荐

数字图像处理实验5图像复原
一.实验目的 (1)了解图像复原的目的及意义,加深对图像复原理论的认识. (2)掌握维纳滤波复原基本原理. (3)掌握约束最小二乘方复原方法. (4)掌握盲解卷积复原方法二.实验内容 (1)维纳滤 ...
数字图像处理实验六：方块编码（BTC）
基于Matlab的图像方块编码(BTC) DIP实验6:方块编码(BTC) 实验目的实验内容参考代码实验结果 DIP实验6:方块编码(BTC) 实验目的掌握方块编码的基本方法及压缩性能. 实验 ...
数字图像处理matlab实验对图像复原,数字图像处理实验07图像的复原处理
数字图像处理实验一.数字图像处理实验实验七图像的复原处理一.实验目的熟悉几种在实际应用中比较重要的图像复原技术,学会用MATLAB复原函数对退化图像进行复原处理. 二.实验内容 1．用点扩散 ...
数字图像处理——第六章彩色图像处理
数字图像处理--第六章彩色图像处理文章目录数字图像处理--第六章彩色图像处理 1 彩色模型 1.1 RGB彩色模型 1.2 CMY 和CMYK彩色模型 1.3 HSI彩色模型 2 伪彩色图像处 ...
数字图像处理第六章——彩色图像处理（上）
数字图像处理第六章数字图像处理---彩色图像处理 (一) 在 MATLAB 中彩色图像的表示 1.1 RGB图像 1.2 索引图像 1.3 处理 RGB图像和索引 (二)彩色空间之间的转换 2.1 ...
matlab数字图像实验报告,数字图像处理实验报告(matlab)
数字图像处理实验报告(matlab) 学院:自动化学院班级:电081班姓名:李林树学号:40850099 2011年10月实验一直方图均衡化一. 实验目的: 1. 熟悉图像数据在计算机中的 ...
图像处理边缘增强matlab,数字图像处理实验 matlab 图像增强边缘检测图像操作.doc...
数字图像处理实验 matlab 图像增强边缘检测图像操作实验1 点运算和直方图处理实验目的 1. 掌握利用Matlab图像工具箱显示直方图的方法 2. 掌握运用点操作进行图像处理的基本原理. ...
数字图像处理第六章 ——彩色图像处理（下）
数字图像处理第六章数字图像处理---彩色图像处理 (五) 彩色图像的空间滤波 5.1 彩色图像的平滑处理 5.2 彩色图像的锐化处理 (六) 直接在 RGB 矢量空间中处理 6.1 使用梯度的彩色边 ...
数字图像处理实验（总计23个）汇总
以下这些实验中的代码全部是我自己编写调试通过的,到此,最后进行一下汇总. 数字图像处理实验(1):PROJECT 02-01, Image Printing Program Based on Half ...
数字图像处理实验——Python语言实现
数字图像处理实验--Python语言实现实验一:数字图像处理入门实验二:直方图均衡实验三:线性平滑和锐化--掩模法实验四:非线性平滑--中值滤波实验五:非线性锐化--梯度法 GitHub地址 ...

数字图像处理实验六图像复原

数字图像处理实验六图像复原相关推荐

最新文章

热门文章

数字图像处理实验六 图像复原

数字图像处理实验六 图像复原相关推荐

最新文章

热门文章

数字图像处理实验六图像复原

数字图像处理实验六图像复原相关推荐