图像处理与分析实验报告

1、 hghu学院实验预习报告班 级 电子51 姓 名 gughu 学 号 2595959 指导教师 杨 课 程 名 称 Matlab应用图像处理 成 绩 实 验 室 电子楼107 实 验 项 目 图像增强与复原 同组实验者 实 验 日 期 2015年12月11日 一、 实验目的1.了解灰度变换增强和空域滤波增强的Matlab实现方法2.掌握直方图灰度变换方法3.掌握噪声模拟和图像滤波函数的使用方法4.了解图像复原的Matlab实现方法 二、主要仪器设备及材料1.计算机2.MATLAB软件3.实验指导书三、 实验原理与步骤1. 灰度变换增强A) 线段上像素灰度分布读入灰度图像debye1.tif，

2、采用交互式操作，用improfile绘制一条线段的灰度值。读入RGB图像flowers.tif，显示所选线段上红、绿、蓝颜色分量的分布B) 直方图变换<i>直方图显示读入图像rice.tif，在一个窗口中显示灰度级n=64，128和256的图像直方<ii>直方图灰度调节利用函数imadjust调解图像灰度范围，观察变换后的图像及其直方图的变化。<iii>直方图均衡化分别对图像pout.tif和tire.tif进行直方图均衡化处理，比较处理前后图像及直方图分布的变化。2. 空域滤波增强A) 噪声模拟利用函数imnoise给图像eight.tif分别添加高斯(g

3、aussian)噪声和椒盐(salt & pepper)噪声。B) 空域滤波<i> 对上述噪声图像进行均值滤波和中值滤波，比较滤波效果。<ii> 总结均值滤波和中值滤波的特点及使用场合。均值滤波器是一种最常用的线性低通平滑滤波器。可抑制图像中的噪声,但同时也使图像变得模糊中值滤波器是一种最常用的非线性平滑滤波器。可消除图像中孤立的噪声点,又可产生较少的模糊<iii> *对图像'saturn.tif'采用'laplacian'高通滤波器进行锐化滤波。3. 图像复原A) 模糊与噪声<i> 对图像flowers.

4、tif分别采用运动PSF和均值滤波PSF进行模糊。<ii> 在上述模糊图像上再添加噪声B) 维纳滤波复原<i> 使用维纳滤波复原函数deconvwnr复原无噪声模糊图像。%非真实PSF %非真实PSF <ii> *使用维纳滤波复原函数deconvwnr复原模糊噪声图像。<iii> *设置信噪比和相关函数的维纳滤波复原。四、原始实验数据1. 灰度变换增强A) 线段上像素灰度分布读入灰度图像debye1.tif，采用交互式操作，用improfile绘制一条线段的灰度值。imshow('debye1.jpg')improfile读入RG

5、B图像flowers.tif，显示所选线段上红、绿、蓝颜色分量的分布imshow('flowers.tif')improfileB) 直方图变换<i>直方图显示读入图像rice.tif，在一个窗口中显示灰度级n=64，128和256的图像直方图。I=imread('flowers.tif');imshow(I)figure,imhist(I,64)figure,imhist(I,128)<ii>直方图灰度调节利用函数imadjust调解图像灰度范围，观察变换后的图像及其直方图的变化。I=imread('flowers.tif

6、9;);imshow(I)figure,imhist(I)J=imadjust(I,0.15 0.9,0 1);figure,imhist(J)figure,imshow(J)I=imread('cameraman.tif');imshow(I)figure,imhist(I)J=imadjust(I,0 0.2,0.5 1);figure,imhist(J)figure,imshow(J)<iii>直方图均衡化分别对图像pout.tif和tire.tif进行直方图均衡化处理，比较处理前后图像及直方图分布的变化。I=imread('pout.tif'

7、);imshow(I)figure,imhist(I)J=histeq(I);figure,imhist(J)figure,imshow(J)I=imread('tire.tif');imshow(I)figure,imhist(I)J=histeq(I);figure,imshow(J)figure,imhist(J)2. 空域滤波增强A) 噪声模拟利用函数imnoise给图像eight.tif分别添加高斯(gaussian)噪声和椒盐(salt & pepper)噪声。I=imread('eight.tif');imshow(I)I1=imnoise

8、(I,'gaussian',0,0.01);figure,imshow(I1)I2=imnoise(I,'salt & pepper');figure,imshow(I2)B) 空域滤波<i> 对上述噪声图像进行均值滤波和中值滤波，比较滤波效果。I=imread('eight.tif');J=imnoise(I,'salt & pepper',0.02);imshow(J)K1=medfilt2(J,3 3); % 3*3中值滤波figure,imshow(K1)K2=filter2(fspecial(

9、'average',5),J)/255;figure,imshow(K2)<ii>总结均值滤波和中值滤波的特点及使用场合。均值滤波器是一种最常用的线性低通平滑滤波器。可抑制图像中的噪声,但同时也使图像变得模糊中值滤波器是一种最常用的非线性平滑滤波器。可消除图像中孤立的噪声点,又可产生较少的模糊<iii> *对图像'saturn.tif'采用'laplacian'高通滤波器进行锐化滤波。I=imread('flowers.tif');imshow(I)h=fspecial('laplacian'

10、;); I2=filter2(h,I);figure,imshow(I2)3. 图像复原A) 模糊与噪声<i> 对图像flowers.tif分别采用运动PSF和均值滤波PSF进行模糊。I=imread('flowers.tif');I=I(10+1:256,222+1:256,:); %剪切图像imshow(I)len=30; %运动位移 theta=10; %运动角度PSF=fspecial('motion',len,theta);blurred=imfilter(I,PSF,'circular','conv');f

11、igure,imshow(blurred)I=imread('flowers.tif');imshow(I) H=fspecial('motion',50,45); %运动PSPmotionblur=imfilter(I,H);figure,imshow(motionblur)H=fspecial('disk',10);%均值滤波PSPaverageblur=imfilter(I,H);figure,imshow(averageblur)<ii> 在上述模糊图像上再添加噪声J=imnoise(motionblur,'salt

12、& pepper');figure,imshow(J)B) 维纳滤波复原<i> 使用维纳滤波复原函数deconvwnr复原无噪声模糊图像。len=30; theta=10; PSF=fspecial('motion',len,theta);blurred=imfilter(I,PSF,'circular','conv');figure,imshow(blurred)wnr1=deconvwnr(blurred,PSF); %真实PSFfigure,imshow(wnr1)%非真实PSF len=30; theta=10

13、; PSF=fspecial('motion',len,theta);blurred=imfilter(I,PSF,'circular','conv');figure,imshow(blurred)wnr2=deconvwnr(blurred,fspecial('motion',2*len,theta); figure,imshow(wnr2)%非真实PSF len=30; theta=10; PSF=fspecial('motion',len,theta);blurred=imfilter(I,PSF,'

14、circular','conv');figure,imshow(blurred)wnr3=deconvwnr(blurred,fspecial('motion',len,2*theta); figure,imshow(wnr3)<ii> *使用维纳滤波复原函数deconvwnr复原模糊噪声图像。PSF=fspecial('motion',len,theta);blurred=imfilter(I,PSF,'circular','conv');noise=0.1*randn(size(I);Blu

15、rredNoisy=imadd(blurred,im2uint8(noise);figure,imshow(BlurredNoisy)wnr4=deconvwnr(BlurredNoisy,PSF);figure,imshow(wnr4) <iii> *设置信噪比和相关函数的维纳滤波复原(p179图9.12)。NSR=sum(noise(:).2)/sum(im2double(I(:).2); %计算信噪比wnr5=deconvwnr(BlurredNoisy,PSF,NSR);figure,imshow(wnr5) %设置信噪比的复原结果NP=abs(fftn(noise).2;

16、 %噪声能量谱密度NCORR=fftshift(real(ifftn(NP); %噪声自相关函数IP=abs(fftn(im2double(I).2; %图像能量谱密度ICORR=fftshift(real(ifftn(IP); %图像自相关函数wnr6=deconvwnr(BlurredNoisy,PSF,NCORR,ICORR);figure,imshow(wnr6) %设置自相关函数的复原结果 指导教师签名 日期 年 月 日Huh 学院实验报告班 级 电子51 姓 名 gughu 学 号 2595959 指导教师 杨 课 程 名 称 Matlab应用图像处理 成 绩 实 验 室 电子楼1

17、07 实 验 项 目 图像增强与复原 同组实验者 实 验 日 期 2015年12月11日 一、实验目的1.了解灰度变换增强和空域滤波增强的Matlab实现方法2.掌握直方图灰度变换方法3.掌握噪声模拟和图像滤波函数的使用方法4.了解图像复原的Matlab实现方法二、主要仪器设备及材料1.计算机2.MATLAB软件3.实验指导书三、实验原理与步骤1. 灰度变换增强A) 线段上像素灰度分布读入灰度图像debye1.tif，采用交互式操作，用improfile绘制一条线段的灰度值。读入RGB图像flowers.tif，显示所选线段上红、绿、蓝颜色分量的分布B) 直方图变换<i>直方图显示

18、读入图像rice.tif，在一个窗口中显示灰度级n=64，128和256的图像直方<ii>直方图灰度调节利用函数imadjust调解图像灰度范围，观察变换后的图像及其直方图的变化。<iii>直方图均衡化分别对图像pout.tif和tire.tif进行直方图均衡化处理，比较处理前后图像及直方图分布的变化。2. 空域滤波增强A) 噪声模拟利用函数imnoise给图像eight.tif分别添加高斯(gaussian)噪声和椒盐(salt & pepper)噪声。B) 空域滤波<i> 对上述噪声图像进行均值滤波和中值滤波，比较滤波效果。<ii>

19、总结均值滤波和中值滤波的特点及使用场合。均值滤波器是一种最常用的线性低通平滑滤波器。可抑制图像中的噪声,但同时也使图像变得模糊中值滤波器是一种最常用的非线性平滑滤波器。可消除图像中孤立的噪声点,又可产生较少的模糊<iii> *对图像'saturn.tif'采用'laplacian'高通滤波器进行锐化滤波。3. 图像复原A) 模糊与噪声<i> 对图像flowers.tif分别采用运动PSF和均值滤波PSF进行模糊。<ii> 在上述模糊图像上再添加噪声B) 维纳滤波复原<i> 使用维纳滤波复原函数deconvwnr复原无噪声模糊图像。%非真