[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Q: about Richardson_Lucy Alg. (deconvolution)
- Subject: Q: about Richardson_Lucy Alg. (deconvolution)
- From: "Herbert H. Tsang" <tsang(at)vcn.bc.ca>
- Date: 15 Jul 1999 22:43:03 GMT
- Newsgroups: comp.lang.idl-pvwave
- Organization: Vancouver CommunityNet
- User-Agent: tin/pre-1.4-19990517 ("Psychonaut") (UNIX) (SunOS/5.6 (sun4u))
- Xref: news.doit.wisc.edu comp.lang.idl-pvwave:15688
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??
-- Herbert
--
-------------------------------+-----------------------------
| Tel: 604.464.8025
H E R B E R T T S A N G | Email: hhpt@yahoo.com
| herbert@canada.com