噪声图像的复原

数字图像处理名称：西征学编号：实验4噪声图像恢复一、实验目的1、熟悉和掌握MATLAB图像处理工具箱的使用；2、了解和掌握常用映像的恢复和增强技术。二、相关知识第一维纳滤波用法：j=deconvwner (I，PSF，nsr)使用维纳滤波算法恢复图I的图像，然后返回到图像j。I是n维数组。PSF是点扩展函数的卷积。NSP是附加噪声的噪声与信号的功率比。例如：I=im 2 double(im read(camera man . TIF)；IMS how(I)；Title(原始图像)；%模拟运动模糊Matlab中文论坛LEN=21theta=11；PSF=fspecial(motion，LEN，THETA)；Blured=imfilter (I，PSF，conv，circular)；Figure，im show (blured)%恢复映像www.iLoveMWnr2=deconvwner (blured _ noisy，PSF)；Figure，imshow(wnr2)Title(中断的存储)双频域滤波器频域滤波器分为低通滤波器和高通滤波器两类。频域低通滤波器的基本思路：G(u，v)=F(u，v)H(u，v)F(u，v)是需要钝化图像的傅立叶变换形式，H(u，v)是选定的低通滤波器变换函数，G(u，v)是通过H(u，v)减少F(u，v)的高频部分而得到的结果，并且使用傅立叶变换功能实现了钝化理想的遍历过滤器(ILPF)具有传递函数。其中是指定的非负值，(u，v)是到过滤器中心的距离。点的轨迹是圆。n阶巴特沃兹低通滤波器(BLPF)(距离原点处基于频率的发生)的传递函数如下与理想通过滤波器不同，巴特沃兹率通过滤波器的传递函数在这里突然不连续。高斯低通滤波器(GLPF)的传递函数包括其中是标准偏差。相应的高通滤波器还包括理想的高通滤波器、n阶巴特沃兹高通滤波器、高斯高通滤波器。给出了低通滤波器的传递函数，利用以下简单关系，得到了相应高通滤波器的传递函数。三、实验阶段1、生成运动模糊图像，使用维纳过滤器恢复图像，显示结果。I=imread(f:1.jpg)I=rgb2gray(s)I=im 2 double(I)；%模拟运动模糊LEN=21theta=11；PSF=fspecial(motion，LEN，THETA)；Blured=imfilter (I，PSF，conv，circular)；还原% imageWnr2=deconvwnr(blurred，PSF)；Subplot(1，2，1)；IMS how(blurred)；Title(运动模糊图像)Subplot(1，2，2)；IMS how(wnr 2)；Title(还原映像)2、利用MATLAB实现频域滤波程序f=im read(room . TIF)；f=FFT 2(f)；%图像的傅里叶变换%转换后，对图像进行对数变化，然后变换其坐标以使其居中s=FFT shift(log(1 ABS(F)；S=g scale(S)；% 0-256范围内的光谱图像缩放Imshow(S)%光谱图像显示h=fs pecial(Sobel)；%生成空间“Sobel”模板Freqz2(h)%查看相应频域过滤器的图像PQ=padded size(size(f)；生成%过滤器时所需大小的矩阵H=freqz2(h，PQ(1)，PQ(2)；%生成频率区域中的“Sobel”过滤器h1=IFFT shift(H)；重新排列数据序列，使%原点位于频率矩阵的左上角Imshow(abs(H)，)%以图形方式显示过滤器Figure，imshow(abs(H1)，)Gs=imfilter(双(f)，h)；使用%模板h过滤空间Gf=dftfilt(f，H1)；使用%过滤器过滤图像的频域Figure，imshow(gs，)Figure，imshow(gf，)Figure、imshow(abs(gs)、)Figure、imshow(abs(gf)、)f=im read(number . TIF)；阅读%图片PQ=padded size(size(f)；生成%过滤器时所需大小的矩阵d0=0.05 * PQ(1)；% gauss设置高通过滤器的阈值H=hpfilter(高斯、pq (1)、pq (2)、d0)；% gauss生成高通滤波器G=dftfilt(f，H)；过滤%图像Figure，imshow(f)%显示原始图像Figure，imshow(g，)%显示过滤后的图像四、实验设备1.电脑2.MATLAB程序；移动内存(u盘等)；4.记录用笔，纸。五、实验报告内容1.叙述实验过程。提交实验的原始图像和结果图像。六、实验报告要求要说明实验的基本阶段，使用数据和图片给出每个阶段获得的实验结果，并进行必要的讨论，必须包括原始图像和相应的计算/处理后图像。七、考试问题1.结合实验评价维纳滤波和频域滤波有什么优点？2.在映像还原过程中需要注意哪些事项？1、生成运动模糊图像，使用维纳过滤器恢复图像，显示结果。I=I=im read(foot ball . jpg)；I=RGB 2灰色(I)；I=im 2 double(I)；LEN=21theta=11；PSF=fspecial(motion，LEN，THETA)；Blured=imfilter (I，PSF，conv，circular)；Wnr2=deconvwnr(blurred，PSF)；Subplot(1，2，1)；IMS how(blurred)；Title(运动模糊图像)Subplot(1，2，2)；IMS how(wnr 2)；Title(还原映像)2、利用MATLAB实现频域滤波ClcClear all，close allf=im read(camera man . TIF)；f=FFT 2(f)；s=FFT shift(log(1 ABS(F)；S=g scale(S)；Subplot(3，3，1)、imshow(S)和title(原始光谱)。h=fs pecial(Sobel)；Freqz2(h)PQ=padded size(size(f)；H=freqz2(h，PQ(1)，PQ(2)；h1=IFFT shift(H)；Subplot (3，3，2)、im show (ABS (h)、 title(频域中的“Sobel”过滤器)；Subplot (3，3，3)、im show (ABS (h1)、和title(重新排列数据序列后的“Sobel”筛选器)；Gs=imfilter(双(f)，h)；Gf=dftfilt(f，H1)；Subplot (3，3，4)、im show (GS，)、title(使用模板h过滤图像空间)；Subplot (3，3，5)、im show (gf，)、title(使用滤镜过滤图像的频域)；Subplot (3，3，6)、im show (ABS (GS)、title(使用模板h空域过滤图像的模式值)；Subplot (3，3，7)、im show (ABS (gf)、title(图像过滤器的频域过滤后的模式值)；f=im read(camera man . TIF)；PQ=padded size(size(f)；d0=0.05 * PQ(1)；H=hpfilter(高斯、pq (1)、pq (2)、d0)；G=dftfilt(f，H)；Subplot(3，