资源描述：
《图像分析与处理 图象恢复(二).ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、图象恢复的滤波方法逆滤波维纳滤波等功率谱滤波逆滤波假定退化图象遵从以下模型在不考虑噪声的情况下写成逆滤波该恢复方法取名为逆滤波。逆滤波逆滤波模型H(u,v)P(u,v)F(u,v)F(u,v)G(u,v)逆滤波实际应用时的缺点：（1）无噪声情况若在频谱平面对图象信号有决定影响的点或区域上，H(u,v)的值为零，故不能确定这些频率处的F(u,v)值，也就难以恢复原始图象f(x,y)。逆滤波（2）有噪声情况G(u,v)=F(u,v)H(u,v)+N(u,v)仍采用逆滤波器P(u,v)=1/H(u,v)作恢复滤波器(a)H(u,v)=0,没有定义。(b)H(u,v)=0
2、附近，H(u,v)较小，N(u,v)/H(u,v)会非常大，结果与大不相同，就不再象。维纳滤波器逆滤波复原方法对噪声极为敏感，要求信噪比较高，通常不满足该条件。因此希望找到一种方法，在有噪声条件下，从退化图像g(x,y)复原出f(x,y)的估计值，该估计值符合一定的准则。用向量f,g,n来表示f(x,y),g(x,y),n(x,y)，Q为对f的线性算子，在约束条件下求Qf的最小化而得到f的最佳估计。基本原理用拉格朗日法求微分，γ可以用来调节以满足约束条件。设Rf和Rn为f和n的相关矩阵：它们是对称矩阵。对于大多数图像而言，相邻象素之间相关性很强，在20个象素之外，
3、趋于零。在此条件下，Rf和Rn可以近似为分块循环矩阵：其中A和B为对角阵，W为酋阵。若QTQ用Rf－1Rn来代替当D为对角阵，分块循环矩阵因此：写成频域形式为：其中Sff(u,v),Snn(u,v)分别是f(x,y)和n(x,y)的谱密度。等功率谱滤波先验假设：图象和噪声均属均匀随机场，噪声的均值为零，且与图象不相关。令Sff(u,v)为信号的功率谱密度，它是信号的自相关函数的傅立叶变换。等功率谱滤波方法：（1）退化模型等功率谱滤波等功率谱滤波图象恢复的代数方法伪逆法投影迭代法伪逆法图象退化的离散模型：当没有噪声时，在平均误差最小的意义上，f的最佳估计为：H+是H
4、的伪逆。伪逆法当没有噪声时，若H是方阵，且存在逆矩阵，则：若H不是方阵，不能求逆。但根据伪逆矩阵理论：伪逆法H+称为H的伪逆矩阵。H+由下列四个条件唯一确定：（1）HH+H=H（2）H+HH+=H+（3）(HH+)T=HH+（4）(H+H)T=H+H这四个条件也可看成伪逆矩阵的定义。伪逆法１．对Ｍ＊Ｎ的矩阵，如Ｍ＞Ｎ，即线性方程组中方程组的数目多于未知量的个数，属“超定”（overdetermined）伪逆法２、对Ｎ＞Ｍ，称为“欠定”（underdetermined）３、对Ｍ＝Ｎ时，伪逆＝逆矩阵。投影迭代法不考虑噪声时，g=Hf投影迭代法投影迭代法表示成下列方程组
5、：投影迭代法令f(0)为f的初始值，下一个估值f(1)是f(0)在超平面a11f1+a12f2+……+a1NfN=g1上的投影。即投影迭代法f(2)是f(1)在超平面a21f1+a22f2+……+a2NfN=g2上的投影，依次求得f(3)……直至f(M)。图象恢复的非线性方法两种有代表性的非线性方法：（一）最大后验法：考虑了图象记录过程所具有的非线性并且把图象各点的集合平均看作是非平稳的随机过程。（二）最大熵法：考虑图象形成时的非线性，并且能保证图象函数有非负值。最大后验法若把图象记录的非线性考虑进去，图象退化模型为最大后验法若把图象记录的非线性考虑进去，图象退化
6、模型为最大后验法约束PSF解卷法h(x,y)-----成象系统的PSF。p(x,y)-----恢复滤波器的冲激响应或PSF。等效PSF：c(x,y)=h(x,y)*p(x,y)。(1)无噪声情况沉鱼落雁约束PSF解卷法(2)有噪声影响约束PSF解卷法约束PSF解卷法约束PSF解卷法框图h(x,y)+p(x,y)f(x,y)f(x,y)+n0(x,y)n(x,y)约束PSF解卷法方法TermsImagerestoration:图象恢复Degrade:退化Autocorrelation:自相关Convolution:卷积Rowstacking:堆叠Pad:补零Peri
7、od:周期Linearalgebra:线性代数TermsCirculantmatrix:循环矩阵Transpose:转置Blockmatrix:分块矩阵Blockcirculantmatrix:分块循环矩阵Impulse:冲激函数Diracdeltafunction:狄拉克函数Impulseresponse:冲激响应Pointspreadfunction(PSF):点扩展函数TermsPseudoinverse:伪逆
