理查森-露西算法是如何工作的？代码示例？

Question

我想弄清楚反褶积是如何工作的。我理解它背后的想法，但我想了解实现它的一些实际算法-算法，将带有点样本函数（模糊核）的模糊图像作为输入，并生成作为输出的潜在图像。

到目前为止，我发现了Richardson–Lucy算法，其中的数学似乎没有那么困难，但我无法理解实际的算法是如何工作的。维基百科上说：

This leads to an equation for which can be solved iteratively according...

但是它没有显示实际的循环。有谁能给我指一个解释实际算法的资源吗。在Google上，我只能找到使用Richardson-Lucy作为步骤之一的方法，而不是实际的Richardson-Lucy算法。

任何语言或伪代码中的算法都是不错的，但是如果Python中有算法，那就太棒了。

先走一步。

编辑

基本上我想知道的是模糊图像（nxm）：

x00 x01 x02 x03 .. x0n
x10 x11 x12 x13 .. x1n
...
xm0 xm1 xm2 xm3 .. xmn

以及用于获取模糊图像的内核（ixj）：

p00 p01 p02 .. p0i
p10 p11 p12 .. p1i
...
pj0 pj1 pj2 .. pji

为了找出原始图像，理查森-露西算法的具体步骤是什么。

Answer 1

下面是一个开源的Python实现： http://code.google.com/p/iocbio/wiki/IOCBioMicroscope