盲去卷积原理及在图像复原的应用

前言

       之前写过一篇维纳滤波在图像复原中的作用，讲述了图像退化模型以及维纳滤波的作用。维纳滤波使用的前提是知道信号和噪声的功率谱，但在实际应用中较难得到，只能根据先验知识进行估计。

       本文介绍盲去卷积复原算法，并在MATLAB中进行实验，和维纳滤波的复原效果进行一个对比。盲去卷积的方法有多种，本文主要介绍由fish提出的基于露西-理查德森(Richardson-Lucy)的盲去卷积算法。

盲去卷积原理

     露西-理查德森算法属于图像复原中的非线性算法，与维纳滤波这种较为直接的算法不同，该算法使用非线性迭代技术，在计算量、性能方面都有了一定提升。

     露西-理查德森算法是由贝叶斯公式推导而来，因为使用了条件概率（即算法考虑了信号的固有波动，因此具有复原噪声图像的能力）。贝叶斯公式如下：

      结合图像退化/复原模型，可以得到迭代函数：

  

    其中 fi 就是第i轮迭代复原图像，对应贝叶斯公式中的p(x)，g是退化函数，对应贝叶斯公式的p(y|x)，c为退化图像（c(y)dy意为在退化图像上积分），如果满足等晕条件，即图像各区域的模糊函数相同，则迭代公式可化简如下：

<