图像退化与复原Word文档格式.docx

图像退化与复原Word文档格式.docx

《图像退化与复原Word文档格式.docx》由会员分享，可在线阅读，更多相关《图像退化与复原Word文档格式.docx（12页珍藏版）》请在冰豆网上搜索。

来实现，要实现线性移不变退化图像的复原，必须知道退化系统的退化函数，即

点扩展函数h（x,y）。

在点扩展函数已知的情况下，常见图像复原方法有逆滤波

和维纳滤波两种。

在考虑噪声的情况下，逆滤波的原理可以表示如下:

F（"

）二^^=F（u,v）+^）

H（u,v）'

/H（u,v）

通常情况下，N（u,v）是未知的，因此即使知道退化模型也不能复原图像

此外，当H（u,v）的任何元素为零或者值很小时，N（u,v）/H（u,v）的比值决定

着复原的结果，从而导致图像复原结果出现畸变。

对丁这种情况，通常采用限制

滤波频率使其难以接近原点值，从而减少遇到零值的可能性。

维纳滤波则克服了逆滤波的缺点，其数学模型表示如下：

F（u,v）顼

|H*（u,v）「

2

（u,v）H（u,v）+Sh（u,v）/Sf（u,v）

G（u,v）

然而，为退化图像的功率谱很少是已知的，因此常常用下面表达式近似:

F（u,v）=

&

1|H*（u,v）|2

:

扩声|H（u,v）|2+k

因此，本实验的内容就是利用上述经典图像复原的原理，对降质退化图像进行复原。

四.实验步骤

本次实验主要包括光电图像的退化模型和复原方法实现两大部分内容。

（一）图像的退化图像

1、大气湍流的建模

1）湍流引起图像退化的光学传递函数（OTF生成。

已知湍流退化模型的OTF表达式如下：

H（u,v）=exp[-k（u2+v2）5/6]（7）

其中，k为一个常数，反映了大气湍流的严重程度。

（u,v）分别代表了（x,y）

方向的频率坐标。

为了生成中心化的OTF可以考虑将式（7）改写为：

22\5/6

H（u,v）=exp[-k（（u-M/2）+（v-N/2）2）]（8）

其中，M,N为图像的长和宽。

2）读入一幅灰度图像，设定式（8）中k=0.0025,进行退化试验。

分别显示原始图像、退化模型和退化图像。

3）设定k=0.001、0.00025重复上一步的试验。

读入原始图像

对原图形进行灰度处理

显示原图像

将上述结果进行fft处理

得到FP

设计退化湍流模型为H

显示传递函数

由FP与H进行相关处理,

得到退化图像

结束

图2大气湍流的退化过程

2、运动模糊的图像退化试验

1）匀速直线运动引起图像退化的光学传递函数（OTF）生成。

已知相机匀

速直线运动的OTF表达式如下:

H（u,v）=Tsin[p（ua+vb）]e-jp（ua+vb）

p（ua+vb）

其中，T为相机曝光时间，a,b分别表示（x,y）方向的速度；

（u,v）分别对

应（x,y）方向的频率坐标。

2）读入一幅灰度图像，设定式（9）中T=1.0,a=b=0.1,编写MATLAB代

码进行模糊退化试验。

要求分别显示原始图像、退化模型和退化图像。

3）设定不同的值，a,b值，重复上一步的试验。

4）利用数字显微镜或其它图像采集设备，移动物体过程中，采集图像。

（二）图像复原试验

1、逆滤波

设计运动模糊模型为H

图3运动模糊的图像退化

显示退化图像

1）根据试验

（一）设计一幅退化图像（包括噪声污染+模糊退化两部分），

其中模糊退化可选高斯模糊、大气湍流模糊或运动模糊（方向可任意指定，如

10度、20度、45度等），噪声模型可自行设定。

2）利用MATLAB编程实现利用全逆滤波方法对退化图像的复原。

要求在同

一个窗口下显示原始退化图像、复原结果及复原结果与理想图像的差值图共3

个图，并对复原结果进行必要的分析。

逆滤波复原公式如下：

其中，G（u,v）为退化图像的傅立叶变换，H（u,v）为退化系统的光学传递函数

（OTF。

3）伪逆滤波：

为了防止逆滤波中H（u,v）过小，使得复原后的图像数据过大

和放大噪声，可采用频谱半径（阈值）限制下的逆滤波方法，即

R

?

V

+

^

其中，R为中心化频谱H（u,v）中某点到原点（零频）的距离或半径。

另一种替代方法是直

其中，（T为一个阈值，用丁限制频谱的幅度值。

这种方法被称为伪逆滤波。

实验要求利用式（11）方式的伪逆滤波重复实验步骤内容2）所涉及的图像。

2、Wiener滤波

1）针对以上逆滤波设计的退化图，编程实现利用Wiener滤波对其进行复

原。

滤波原理如下：

（13）

F（u,v）=|5虫冬G（u,v）

U（u,v）H（u,v）+k

（OTF,k为一个与信噪比有关的调节因子。

要求在同一个窗口下显示理想图

像（退化前）、退化图像、复原结果等共3个图，并对复原结果进行必要的分析。

2）改变k值，重复试验内容1）。

以上应根据原理自行编写代码，不允许直接调用MATLAfi带的deconvwnr（）

函数

显示原始

退化图像

分别显示三种

方式下的复原

结果F1

图4全逆，伪逆，wiener滤波复原过程

五.实验结果及分析

分析：

由上述结果可知，大气湍流会使图像变得模糊，而k值越大，其模糊效果越明显。

由上述结果可知，随着a,b的值变大，图像模糊变得明显，人眼看起来好像是由丁运动速度过快造成的模糊。

3、图像复原试验

图10,图11,图12分别为全逆，伪逆，wiener对运动模糊（a=b=0.1）滤波的结果，从中可以看出wiener是三者中对运动模糊复原效果最好的滤波方式，且wiener中k值越小复原效果越好。

六.实验心得体会和建议

心得体会：

通过这次实验使我了解了图像退化的原因，以及相关的退化模

型，并学会以matlab为平■台利用退化模型对图像进行退化处理以及退化图像的

复原处理。

建议：

可以要求利用C或C+辿行图像的退化与复原处理。

七.程序源代码

%title:

atmosphere

%explain:

本程序利用大气湍流模型对理想图像进行退化f=imread（'

3.jpg'

）;

figure

（1）

subplot（131）,imshow（f）,title（'

f=rgb2gray（f）;

Fp=fft2（f）;

原始图彳'

）

[m,n]=size（f）;

%绘制网格点[v,u]=meshgrid（1:

n,1:

m）;

u=u-floor（m/2）;

v=v-floor（n/2）;

k=0.00025;

Duv=u.A2+v.A2;

H=exp（-k.*Duv.A（5/6））;

G=H.*fftshift（Fp）;

f1=abs（ifft2（G））;

nchar=num2str（k）;

ltext=strcat（'

k='

nchar）;

%

标题注释

subplot（132）,imshow（H）,title（['

subplot（133）,imshow（f1,[]）,title（'

传递函数'

ltext]）;

退化图像'

）；

move

本程序利用运动模糊模型对理想图像进行退化f=imread（'

[v,u]=meshgrid（1:

T=1.0;

a=0.3;

b=0.3;

%a=0.1,b=0.1;

%a=0.01,b=0.01;

z=pi*（u*a+v*b）+eps;

H=T./z.*sin（z）.*exp（-1j*z）;

nchar=num2str（a）;

a=b='

传递函数'

退化图像'

recovery

本程序利用运动模糊对加入椒盐噪声的理想图像进行退化，之后分

别用全逆滤波

%伪逆滤波，wiener滤波对设计的退化图进行处理，观察三种滤波的复原效果。

clc,clearall,closeall;

f=imread（'

4.jpg'

figure,imshow（f）,title（'

原始图像'

f1=imnoise（f,'

salt&

pepper'

0.02）;

%加入椒盐噪声后的图像

%画网格点

%====================================================================

%运动模糊

%T=1.0,a=0.1,b=0.1;

%T=1.0,a=0.01,b=0.01;

a=0.1;

b=0.1;

G=H.*fftshift（fft2（f1））;

%退化图像的频域

f2=real（ifft2（ifftshift（G）））;

%加了椒盐噪声的退化图像

G1=fftshift（fft2（f2））;

%求原始退化图像的频域表示

以逆滤波

F=G1./H;

F1=real（ifft2（ifftshift（F）））;

以反中心化，反傅立叶变换取实部得到复原结果图（包含噪声）

F2=F1-double（f）;

%复原结果与理想图像的差值

figure,

subplot（131）,imshow（f2,[]）,title（'

原始退化图像'

subplot（132）,imshow（F1,[]）;

title（'

复原结果'

subplot（133）,imshow（F2,[]）;

差值图'

%=======================

%»

逆滤波

%=======================Duv=u.A2+v.A2;

ifDuv.A（1/2）<

=1000

P=1./H;

else

P=0;

end

F=G1.*P;

原始退化图像'

复原结果'

差值图'

%====================================================================

%wiener滤波

%====================================================================k1=0.000025;

F=（（1./H）.*（abs（H）A2）./（（abs（H）A2）+k1））.*G1;

%反中心化，反傅立叶变换取实部得到复原结果

图（包含噪声）

八.思考题

1.简要叙述图像退化的原因。

答：

原因：

在图像的获取和传输过程中，由丁成像系统各种因素的影响，使得最终得到的图像相比丁理想图像质量降低。

2.图像复原与图像增强有什么异同？

表格1图像复原与图像增强的异同

图像复原

图像增强

目的

改善图像质量。

过程

及

评价

不需要考虑退化模型；

客观过程，恢复原始图像的最优估值。

主观过程，提供便于人眼观察或机器识别。

实现

手段

空域（反卷积）或频域滤波

空域（卷积）或频域滤波

3.简述同态滤波的设计思想和处理流程，说明该方法适合丁何种场合下的图像恢复。

设计思想：

通过把频率过滤和灰度变换相结合进行图像处理，依靠图像的照度/反射率模型作为频域处理的基础，利用压缩亮度范围和增强对比度来改善图像的质量。

处理流程：

1）对原图取对数得到如下两个加性分量：

lnf（x,y）=lnfi（x,y）+lnf「（x,y）

（14）

2）频域表小为:

FFTlnfx,y=FFTlnfix,y+FFTlnfrx,y（15）

3）设计一个滤波器H（u,v）进行滤波处理

4）傅立叶反变换，即IFFT。

5）求指数结果，得到复原结果。

适用场合：

作为非线性复原方法的同态滤波，它可以用来恢复在非均匀光照射下的图像。

报告评分:

指导教师签字: