前言
之前写过一篇维纳滤波在图像复原中的作用,讲述了图像退化模型以及维纳滤波的作用。维纳滤波使用的前提是知道信号和噪声的功率谱,但在实际应用中较难得到,只能根据先验知识进行估计。
本文介绍盲去卷积复原算法,并在MATLAB中进行实验,和维纳滤波的复原效果进行一个对比。盲去卷积的方法有多种,本文主要介绍由fish提出的基于露西-理查德森(Richardson-Lucy)的盲去卷积算法。
盲去卷积原理
露西-理查德森算法属于图像复原中的非线性算法,与维纳滤波这种较为直接的算法不同,该算法使用非线性迭代技术,在计算量、性能方面都有了一定提升。
露西-理查德森算法是由贝叶斯公式推导而来,因为使用了条件概率(即算法考虑了信号的固有波动,因此具有复原噪声图像的能力)。贝叶斯公式如下:
结合图像退化/复原模型,可以得到迭代函数:
其中 fi 就是第i轮迭代复原图像,对应贝叶斯公式中的p(x),g是退化函数,对应贝叶斯公式的p(y|x),c为退化图像(c(y)dy意为在退化图像上积分),如果满足等晕条件,即图像各区域的模糊函数相同,则迭代公式可化简如下:
<