[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Q: about Richardson_Lucy Alg. (deconvolution)
- Subject: Re: Q: about Richardson_Lucy Alg. (deconvolution)
- From: edward.s.meinel(at)aero.org
- Date: Fri, 16 Jul 1999 13:06:39 GMT
- Newsgroups: comp.lang.idl-pvwave
- Organization: Deja.com - Share what you know. Learn what you don't.
- References: <7mlo5n$pps$1@sylvester.vcn.bc.ca>
- Xref: news.doit.wisc.edu comp.lang.idl-pvwave:15696
In article <7mlo5n$pps$1@sylvester.vcn.bc.ca>,
"Herbert H. Tsang" <tsang@vcn.bc.ca> wrote:
>
> Anyone is familiar w/ the Richardson_Lucy Deconvolution alg.? I have
found
> a IDL program (lucy_guess.pro) that claims to perform this
deconvolution.
>
> However, when examining the code I am not quite sure if it is correct.
in
> lucy_guess.pro line 89 it has this loop
> for k=0,n_iter do begin
> weight = convol(old,psfr)
> new = old * convol( (data/weight) ,psf)
> old = new
> endfor
>
> but shouldn't it be
> for k=0,n_iter do begin
> weight = convol(old,psfr,/edge_truncate)
> new = old * correlate( (data/weight) ,psf)
> old = new
> endfor
>
> Any idea??
You are correct about the correlation step: the inner operation is a
convolution used to simulate the imaging process; the outer operation is
a correlation. Whether to use /EDGE_TRUNCATE or not is up to the user.
Ed Meinel
Sent via Deja.com http://www.deja.com/
Share what you know. Learn what you don't.