Estou tentando descobrir como funciona a deconvolução. Entendo a idéia por trás dela, mas quero entender alguns dos algoritmos reais que a implementam - algoritmos que tomam como entrada uma imagem borrada com sua função de amostra de ponto (blur kernel) e produzem como saída a imagem latent

Até agora eu encontrei Richardson – Lucy algoritmo em que a matemática não parece ser tão difícil, mas não consigo entender como o algoritmo funciona. Na Wikipedia, diz:

Isso leva a uma equação para a qual pode ser resolvida iterativamente de acordo ...

como ele não mostra o loop real. Alguém pode me apontar para um recurso onde o algoritmo real é explicado. No Google, só consigo encontrar métodos que usam Richardson – Lucy como uma de suas etapas, mas não o algoritmo atual de Richardson – Lucy.

algoritmo em qualquer idioma ou pseudo-código seria bom, no entanto, se um estiver disponível em Python, isso seria incríve

gradecemos antecipadament

Edita

Essencialmente, o que eu quero descobrir é a imagem borrada (nxm):

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

e o kernel (ixj) que foi usado para obter a imagem borrada:

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

Quais são as etapas exatas do algoritmo Richardson – Lucy para descobrir a imagem origina