图像复原_先建立退化的数学模型,然后根据该模型对退化图像进行

1、第六章 图像复原第六章 图像复原（Image Restoration）图像增强：用适当方式改善图像质量，增强图像的视觉效果，以适应人眼的视觉和心理，不用考虑增强处理后的图像是否符合原有图像，是否失真。图像复原：试图利用退化过程的先验知识，去恢复已被退化图像的本来面目。成像系统受各种因素的影响，导致了图像质量的降低，或者说是退化。由于获得图像的方法不同（光学、光电子或电子等），有多种退化形式，都使成像的分辨率和对比度退化，例如：传感器噪声摄像机聚焦不佳物体与摄像机之间的相对移动随机大气湍流光学系统的象差成像光源和射线的散射图像退化的主要表现形式：图像模糊图像受到干扰由于成像系统造成图像退化的典型

2、现象是模糊，所以图像复原的一个基本任务就是去模糊。第六章 图像复原由于图像复原是建立在比较严格的数学推导之上的，存在较复杂的数学运算是本章特点。图像复原的基本思路：先建立退化的数学模型，然后根据该模型对退化图像进行拟合。图像复原模型可以用连续数学和离散数学处理，处理项的实现可在空间域卷积，或在频域相乘。第六章 图像复原6.1 图像退化模型14设一成像系统的物像映射关系（或说退化图像）为g(x,y)=Tf(x,y)式中g(x,y)：输出的退化函数或称退化图像；f(x,y)：输入的图像函数；Ti：成像系统作用的运算符，可以认为是线性算子。按照线性系统的性质：复杂激励f(x,y)无数个函数激励之和f

3、(x,y)用二维(x,y)函数的二利用函数的选择性，可以把图像维卷积表示出来+f(x,y)=退化图像： f(,)(x,y)dd。+g(x,y)=Tf(,)(x,y)dd+=-f(,)T(x,y)dd+ f(,)h(x,y,)dd因为f(,)是加在(x,y)上的权重因子，按线性性质，符号Ti可移到积分号内而直接作用于函数，得系统脉冲响应（点扩散函数，PSF: Point Spread Function）：第六章 图像复原h(x,y,)=T(x,y)结论：首先求出系统对基元函数的响应（输出）表达式（即系统脉冲响应），再乘以相应权重因子后求和，就可以得到退化图像。脉冲响应只依赖(x)和(y)，是具有

4、移不变性质：h(x,y,)=h(x,y)据此，我们得到+g(x,y)=f(,)h(x,y)ddf(x,y)h(x,y) （卷积形式）正是点扩散函数h使图像退化。这就是近代光学中用成像概念来描述的退化图像（degraded image）的数学表达式。If h(x,y,)=h1(x,)h2(y,)， i.e. Ti是一个h可分离系统，计算时的二维问题就变成行和列的两次一维计算。6.1.1 空间域图像退化模型（Image Degraded Model）实际问题中，我们还要考虑图像噪声，if image noise is为加性的，得图像退化模型+n(x,y),且g(x,y)=f(,)h(x,y)dd+

5、n(x,y)如上所述：这个系统模型的特点是线性、移不变的（或图像领域常称空间不变）。第六章 图像复原冲响应），加上了信号噪声，输出退化图像。若已知h和n(x,y)，经过反演运算，可以得到一个近似于原图像也就是说，输入图像函数f(x,y)，通过系统T（这里h为系统T的脉(x,y)f(x,y)的最佳估计, 故求的复原图像（Restored image）f(x,y)的过程就是Image Restoration f(x,y)的最佳估计f(x,y)=T1g(x,y)n(x,y)f例 运动模糊的degraded Model运动模糊：成像系统与物体间相对运动造成的像模糊。设只有像函数f(x,y)相对系统的移

6、动。设x(t)：x方向的移动分量；y(t)：y方向的移动分量。胶片上的总曝光量为在快门开闭时间T的积分，即运动模糊成像表示为：g(x,y)=fxx(t),yy(t)dt0T且在曝光时间T内的再简化一次，令物体仅在x方向上作匀速直线运动，总位移量为a，物体沿x方向的变换分量为ax(t)=t, T第六章 图像复原有ag(x,y)=f(xt)dt=g(x)T0alet t1=t, thenTTg(x,y)=g(x)=f(xt1)dt1=f(x)h(x)a0where: h(x)=aT T, 0xa 为沿x方向造成运动模糊的点扩散a函数。即，运动物体沿x方向移动时的图像退化模型为：T 0xah(x)=

7、a。 0 elsewhere或因成像系统的离焦造成点像弥散成一个圆盘: 6第六章 图像复原hr0h(x,y)=00 elsewhereh0为点光强；r0为圆盘半径6.1.2 频域退化模型首先，将空间域退化模型写成卷积形式g(x,y)=f(x,y)h(x,y)+n(x,y)After Fourier Transform,G(u,v)=F(u,v)iH(u,v)+N(u,v)其中H(u,v)又称系统的传递函数。从频域看，由于系统传递函数的退化导致输出产生退化。1) 运动模糊的传递函数if PSF 为一矩形函数，且是沿x方向移动的运动模糊的点扩散函数，则系统点扩散函数的FT代表了系统因运动模糊的传递

8、函数 TjuasinuaeH(u)= ua第六章 图像复原2） 模糊的传递函数5当光学成像系统离焦时，点光源的像成圆盘，其传递函数J1(r0)H(u,v)=2 r0=, J1为一阶贝塞尔函数（振荡衰减函数），令wherex=r0，则1xk1J1(x)=()k!(nk)!+2k=02k+1xx3x5=+ 216384该传递函数在以原点为中心，r0的倍数为半径处存在零点，形成一些同心暗环，由离焦图像的频谱上估计出这些同心圆的半径，即可确立离焦模糊的传递函数。第六章 图像复原3）大气湍流模糊的传递函数5在长时间曝光时，大气湍流将使图像产生模糊，这种传递函数可近似为高斯函数22H(u,v)=expc(

9、u+v)式中c为与湍流性质相关的常数。6.1.3 能量域的退化模型相关定理卷积定理的一个特例6这个定理构成了函数在空域与频域的一种联系。两个函数相关互相关函数，这里+f(x) g(x)=f()g(x+)d，这里类似卷积计算，但不是g(x)， 不需要折叠操作，只需要直接使g(x)移动、与f(x)相乘并积分即可。对应Fourier 变换+f(x) g(x)e+j2uxdx=G(u)F(u)可证G(u)F(u)证明： g(s)ej2usds。 f()ej2ud第六章 图像复原+j2uxj2ux f( )g(x+)dedxf(x) g(x)edx=+j2ux =f()g(x+)edxd+=这里因为+f

10、()ej2uG(u)d=G(u)F(u)g(x+)ej2uxdx 中令s=x+ ，可得+g(s)ej2usej2uds=ej2uG(u)。把相关定理写成傅立叶变换对形式：f(x) g(x)G(u)F(u)反之f(x)g(x)F(u) G(u)若f(x)=g(x)，是同一函数，f(x) g(x)表示自相关函数 f(x) f(x)F(u)F(u)=F(u)2。对二维的情况，如果函数，相关定义为 f(x,y)andg(x,y)是连续变量x,y的第六章 图像复原+f(x,y) g(x,y)=相关定理可表示成： f(,)g(x+,y+)ddf(x,y) g(x,y)G(u,v)F*(u,v)f(x,y)

11、g(x,y)G(u,v) F(u,v)and(x,y)看成是大小分别为 g 把f(x,y) AB and CD的离散数列，在x,y方向上延拓这些数列为每个周期M,N。为了避免交叠效应的产生，令周期为：M =A+ C1N=B+D1此时，f(x,y)andg(x,y)有如下形式：f(x,y) 0xA1 and 0yB1fe(x,y)=0 AxM1 and ByN1g(x,y) 0xC1 and 0yD1 ge(x,y)=0 CxM1 and DyN1则二维离散相关可表示为：fe(x,y) ge(x,y)=fe(m,n)ge(x+m,y+n) m=0n=0M1N1Where x=0,1,2,.,M1

12、; y=0,1,2,.,N1.6.1.4 能量表示1,6如果把图像作为二维随机过程的一个样本，它的功率谱密度定义为 11第六章 图像复原其自相关的傅里叶变换，设Pg：已退化图像的功率谱Pf：理想图像（输入图像）的功率谱 ： Pf=F(u,v)2H(u,v)：系统传输函数按照线性系统的基本原理，对于具有系统传输函数H(u,v)的线性系统，其输出功率谱为Pg=PfH(u,v)考虑噪声影响，退化模型的能量表达式为 2Pg=PfH(u,v)+Pn Pn为噪声的功率谱如何得到Pf和Pn？ 2谱的统计特性由实验测定： 加性噪声的Pn谱特点：位于高频段、且幅值为常数的均匀分布（；第六章 图像复原 如果成像系

13、统是理想的（理想化到衍射现象都没有），点扩散函数被看成是函数，则传递函数H(u,v)1，所以实际得到的图像功率谱为 PgPf。由谱的统计特性知，它集中在低频段，频率从零到某一频率数值，幅值由高到低平滑下降。上图的虚线段2AB显示的是Pf； 实际上，H(u,v)1，PfH(u,v)的结果使Pf的趋势下降，而且下降不平滑，因而实际图像的功率曲线由ABC所示。现在，由AB段来估算AB段，作为Pf的近似，是可行的。6.1.5 MATLAB中的模糊与噪声7若不存在噪声，图像的模糊结果就完全来源于点扩散函数PSF的作用，通过精确描述失真的PSF对模糊图像进行卷积操作，就可实现去模糊。分别创建一个仿真运动模

14、糊PSF和一个均值滤波PSF并模糊原图像 % 使用函数fspecial创建点扩散函数PSFI=imread(girl.bmp);%I=I(10+1:256,100+1:256,:); % 裁剪图像subplot(1,2,1);imshow(I);title(original image)LEN=40; % 指定运动位移为40个像素THETA=30; 运动角度为30度PSF=fspecial(motion,LEN,THETA);Blurred=imfilter(I,PSF,circular,conv);第六章 图像复原subplot(1,2,2);imshow(Blurred); title(Blurred image)H=fspecial(disk,10); % 均值滤波PSF blurredfilter=imfilter(I,H); figure,imshow(blurredfilter); title(Filter Blurred image)第六章 图像复原本章参考文献1 黄贤武等，数字图像处理与压缩编码技术，电子科大版，2000，