图像复原课件

第5章图像复原第5章图像复原1图像复原－又称为图像恢复与图像增强相似－都要得到在某种意义上改进的图像，或者说，希望要改进输入图像的视觉质量不同之处－图像增强技术一般要借助人的视觉系统的特性，以取得看起来好的视觉结果，而图像复原则认为图像是在某种情况下退化或恶化了(图像品质下降了)，现在需要根据相应的退化模型和知识重建或恢复原始的图像图像恢复技术是要将图像退化的过程模型化，并据此采取相反的过程以得到原始的图像图像复原－又称为图像恢复与图像增强相似－都要得到在某种意义上2图像恢复的内容退化模型和循环矩阵对角化复原的代数方法逆滤波最小二乘方滤波交互式恢复空间复原技术图像恢复的内容退化模型和循环矩阵对角化3退化模型和循环矩阵对角化退化模型产生原因光学系统中的衍射传感器非线性畸变光学系统的像差摄影胶片的非线性大气流的扰动效应图像运动造成的模糊几何畸变退化模型和循环矩阵对角化退化模型产生原因光学系统中的衍射传感4定义：f[x,y]:原始图像g[x,y]:退化图像n[x,y]:加性噪声g[x,y]=H{f[x,y]}+n[x,y]H{}:系统或操作图像恢复就是在给定g(x,y)和代表退化的H的基础上，得到对f(x,y)的某个近似的过程定义：f[x,y]:原始图像图像恢复就是在给定g(x,y)5H{.}+f(x,y)n(x,y)g(x,y)简单的通用退化模型H{.}+f(x,y)n(x,y)g(x,y)简单的通用退化6线性：

H{k1f1+k2f2}=k1H{f1}+k2H{f2}相加性：令k1=k2=1，则H{f1+f2}=H{f1}+H{f2}一致性：令f2=0，则H{k1f1}=k1H{f1}位置（空间）不变性：

H{f[x-a,y-b]}=g[x-a,y-b]H多具有的性质2幅图像常数图像任意位置的响应只与在该位置的输入值有关，而与位置本身无关线性：H多具有的性质2幅图像常数图像任意位置的响应只与7常见具体退化模型示例空间不变线性摄影胶片的冲洗过程非线性光学成像系统，由于孔径衍射产生的退化目标运动造成的模糊退化模糊退化随机噪声迭加，随机性的退化常见具体退化模型示例空间不变线性摄影胶片的冲洗过程非线性光学8退化模型的计算假设对2个函数f(x)和h(x)进行均匀采样，其结果放到尺寸为A和B地2个数组。对f(x)，x的取值范围是0,1,2…A-1；对h(x)，x的取值范围是0,1,2,….B-1。利用卷积计算g(x)。为了避免卷积的各个周期重叠，取M≥A+B-1，并将函数用0扩展补齐fe(x)和he(x)表示扩展函数，卷积为ge(x)=∑fe(m)

he(x-m)

x=0,1,…M-1矩阵表示g=Hfg和f是M维列矢量：fT=[f[0],f[1],…,f[M-1]]gT=[g[0],g[1],…,g[M-1]]退化模型的计算假设对2个函数f(x)和h(x)进行均匀采样，9H称为M×M循环矩阵H=g=Hf+n(1)考虑噪声H称为M×M循环矩阵H=g=Hf+n(1)10如果直接对式(1)进行计算求解f，计算量达，如M=N=512,则H的尺寸为262144×262144，可以通过对角化H来简化当k=0,1…M-1时，循环矩阵H(设为M×M)的特征矢量和特征值分别为循环矩阵对角化如果直接对式(1)进行计算求解f，计算量达，如M=N=51211将H的M个特征矢量组成1个M×M的矩阵WW=[w(0)w(1)w(2)…w(M-2)w(M-1)]H=WDW-1

D=W-1HW其中：D(k,k)=λ(k)，WW-1=W-1W=I将H的M个特征矢量组成1个M×M的矩阵WW=[w(0)12复原的代数方法图像复原的主要目的是当给定退化的图像g以及H和n的某种假设，估计出原始图像f代数复原方法的中心是寻找一个估计的f^，它使事先确定的某种优度准则为最小复原的代数方法图像复原的主要目的是当给定退化的图像g以及H和13无约束复原方法由退化模型可知，其噪声项为：n=g-Hf在并不知道n的情况下，希望找到一个f^，使得Hf^在最小二乘方意义上来说近似于g，也就是说，寻找一个f^，使得

||n||2=||g–Hf^||2

J(f^)无约束复原方法由退化模型可知，其噪声项为：||n||2=14nTn=(g–Hf^)T(g–Hf^)或实际上是求J(f^)的极小值问题，除了要求J(f^)为最小外，不受任何其它条件约束，因此称为无约束复原dJ(f^)/df^=0=-2HT(g–Hf^)即f^=(HTH)-1HTgM=N时，则有f^=H-1(HT)-1HTg=H-1g(2)nTn=(g–Hf^)T(g–Hf^)或实际上是15约束复原方法在最小二乘方复原处理中，为了在数学上更容易处理，常常附加某种约束条件。如令Q为f的线性算子，最小二乘复原问题可看成是使形式为||Qf^||2函数，服从约束条件||g–Hf^||2=||n||2的最小问题，这种带有附件条件的极值问题可用拉各朗日乘数法处理约束复原方法在最小二乘方复原处理中，为了在数学上更容易处理，16处理过程寻找一个f^，使下述准则函数为最小J(f^)=||Qf^||2+α{||g–Hf^||2-||n||2}拉各朗日系数dJ(f^)/df^=0f^=(HTH+γQTQ)–1HTgα=1/λ处理过程寻找一个f^，使下述准则函数为最小J(f^)=|17逆滤波基本原理在不考虑噪声的情况下，假设M=N，则根据前面的公式，有f^=H–1g=WD–1W–1g或W–1f^=D–1W–1gF^(u,v)=G(u,v)/H(u,v)傅立叶变换经过傅立叶反变换，可求得原始图像f(x,y)反向滤波的作用逆滤波基本原理在不考虑噪声的情况下，假设M=N，则根据前面的18在有噪声的情况下F^(u,v)=F(u,v)+N(u,v)/H(u,v)从上面两式可以看出，在进行复原处理时可能会发生下列情况：H(u,v)=0或H(u,v)非常小，在这种情况下，即使无噪声，也无法精确恢复f(x,y)在有噪声存在时，在H(u,v)的邻域内，H(u,v)的值可能比N(u,v)的值小的多，由上式得到的噪声项可能会非常大，不能使f(x,y)正确恢复在有噪声的情况下F^(u,v)=F(u,v)+N(u19一般说，逆滤波不能正确估计H(u,v)的零点实际中，不用1/H(u,v)，而用另外一个关于u,v的函数M(u,v)处理框图为：f(x,y)H(u,v)+M(u,v)G(u,v)N(u,v)F^(u,v)M(u,v)=1/H(u,v)u2+v2≤w201u2+v2＞w20一般说，逆滤波不能正确估计H(u,v)的零点处理框图为：f(20例子原始图像散焦模糊利用原始图像的一个邻域光谱面恢复利用大的邻域进行恢复例子原始图像散焦模糊利用原始图像的一个邻域光谱面恢复利用大的21消除匀速直线运动引起的模糊在获取图像时，由于景物和摄像机之间的相对运动，造成图像的模糊。假设对平面匀速运动的景物采集1幅图像f(x,y)，并设x0(t)和y0(t)分别是景物在x和y方向的运动分量，T是采集时间长度。实际采集到的由运动而造成的模糊图像g(x,y)为：g(x,y)=∫0Tf[x-x0(t),y-y0(t)]dtG(u,v)=F(u,v)

0Texp{-j2p[ux0(t)+vy0(t)]}dt=F(u,v)H(u,v)H(u,v)=

0Texp{-j2p[ux0(t)+vy0(t)]}dt如果知道运动分量x0(t)和y0(t)，从上式直接得到H(u,v)消除匀速直线运动引起的模糊在获取图像时，由于景物和摄像机之间22只考虑x方向的运动情况，即x0(t)=ct/T，y0(t)＝0，则H(u,v)=

0Texp{-j2πux0(t)}dt=(T/πuc)sin(πuc)e-jπuc当n为整数时，H在u=n/c处为0根据模糊图像g(x,y)，不考虑y的变化g(t)=

0Tf[x-x0(t)]dt=

0Tf[x-ct/T]dt=

x-cxf[t]dt,0

x

L其中：t=x–ct/T.因此，g’(x)=f(x)–f(x-c) or f(x)=g’(x)+f(x-c)只考虑x方向的运动情况，即x0(t)=ct/T，y0(t)＝23假设L=Kc,其中K使整数.变量x为： x=z+mc其中z在[0,c]范围，m是(x/a)的整数部分(m的取值为0,1,…,K-1).则, f(z+mc)=g’(z+mc)+f[z+(m-1)c]这个等式能通过递归f(z)解决，f(z)=f(z-a)定义为运动范围的一部分0

z<c. f(x)=Sk=0mg’(x-kc)+f(x-mc)由于g(x)为已知，问题就转化为估计f(x)假设L=Kc,其中K使整数.变量x为：24为了估计

f(x): f(x-mc)=f(x)–f^(x).对每个kc

x<(k+1)c计算，并将k=0,1,…,K-1的结果加起来，得到

Kf(x)=Sk=0K-1f(x+kc)–Sk=0K-1f^(x+kc),0

x<c f(x)

A–(1/K)Sk=0K-1f^(x+kc),0

x<c其中(1/K)Sk=0K-1f(x+kc)是未知的，但是当K很大时，其接近平均值，将其设为常数A。因此,

f(x-mc)

A–(1/K)Sk=0K-1f^(x+kc-mc),0xL A–(1/K)Sk=0K-1Sj=0kg’[x-mc+(k-j)c]

或对0

x,y

Lf(x,y)A–(1/K)Sk=0K-1Sj=0kg’[x-mc+(k-j)c,y]+Sj=0mg’[x-jc,y]为了估计f(x): f(x-mc)=f(x)–f25最小二乘方滤波最小二乘滤波也就是维纳滤波，它是使原始图像f(x,y)及其恢复图像f^(x,y)之间的均方误差最小的复原方法具体的数学公式推导过程忽略，直接给出公式Sf(u,v):为

f[x,y]的功率普，Sh(u,v)为n[x,y]的功率普最小二乘方滤波最小二乘滤波也就是维纳滤波，它是使原始图像f(26讨论一下上式的几种情况（1）如果s=1，方括号中的项就是维纳滤波器（2）如果s是变量，就称为参数维纳滤波器（3）当没有噪声时，Sn(u,v)=0，维纳滤波器就退化为理想的逆滤波器（4）当Sn(u,v)和Sf(u,v)未知时，用常数K可代替因此必须调节s以满足f^=(HTH+sQTQ)–1HTg讨论一下上式的几种情况（1）如果s=1，方括号中的项就是维纳27逆滤波和维纳滤波恢复比较110100SNR退化图像

傅立叶功率普逆滤波恢复维纳滤波恢复光谱图

逆滤波和维纳滤波恢复比较110100SNR退化图像28原始图像逆滤波恢复模糊和增加噪声约束的最小二乘滤波原始图像逆滤波恢复模糊和增加噪声约束的最小二乘滤波29交互式恢复前面讨论都是自动解析的恢复方法，在具体恢复工作中，常常需要人机结合，由人来控制恢复过程，以达到一些特殊的效果实际中，有时图像会被1种2-D的正弦干扰模式(也叫相关噪声)覆盖。令η(x,y)代表幅度为A，频率分量为(u0,v0)的正弦干扰模式，即：η(x,y)=Asin(u0x+v0y)傅立叶变换为：N(u,v)=-jA[δ(u-u0/(2π),v-v0/v(2π))-δ(u+u0/(2π),v+v0/(2π))]交互式恢复前面讨论都是自动解析的恢复方法，在具体恢复工作中，30这里退化仅由噪声造成，所以有：G(u,v)=F(u,v)+N(u,v)利用前面讲的带阻滤波器消除，以去掉正弦干扰模式影响正弦干扰傅立叶光谱恢复图像n(x,y)=Asin(u0x+v0y)这里退化仅由噪声造成，所以有：G(u,v)=F(u,v)+N31例子周期干扰退化图像傅立叶谱例子周期干扰退化图像傅立叶谱32干涉模式处理后图像干涉模式处理后图像33空间复原技术空间变换灰度插值空间复原技术空间变换34空间变换在图像的获取或显示过程中，产生几何失真，如成像系统有一定的几何非线性，因此会造成如图所示的枕形失真或桶形失真，另外，由于地球表面呈球形，摄取的平面图像也将会有较大的几何失真。对这些图像必须加以校正，以免影响分析精度原始图像枕形失真桶形失真空间变换在图像的获取或显示过程中，产生几何失真，如成像系统有35校正