[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

*Subject*: Re: filtering problem*From*: "Jaco van Gorkom" <gorkom(at)rijnh.nl>*Date*: Tue, 21 Nov 2000 13:11:29 -0000*Newsgroups*: comp.lang.idl-pvwave*Organization*: UUNET WorldCom server (post doesn't reflect views of UUNET WorldCom*References*: <3A13B236.9A321113@clinmed.gla.ac.uk> <3A158C96.7205F78B@vni.com> <3A15B24E.F0158A4@astro.cornell.edu> <3A1A445B.14A912A3@clinmed.gla.ac.uk>*Xref*: news.doit.wisc.edu comp.lang.idl-pvwave:22283

"Dave Brennan" <9147261b@clinmed.gla.ac.uk> wrote in message 3A1A445B.14A912A3@clinmed.gla.ac.uk">news:3A1A445B.14A912A3@clinmed.gla.ac.uk... ... > Although any value below the thershold should not be included in the > statistics, pixels below the threshold should be corrected by the algorithm. > > Therefore is it possible to change the code: > > pro thresh, a, n, t > m = a ge t > wh = where(m,cnt) > if cnt eq 0 then return > a[wh] = a[wh] * mean(a[wh]) * (smooth(float(m),n,/EDGE) / $ > (smooth(a*m,n,/EDGE)+1.e-30))[wh] > end > > to allow this. If I understand you correctly, you want to correct the below-threshold pixels by the factor calculated from the valid pixels surrounding it. Couln't you just leave out some of the specific subscripting [wh] in the final assignment? a = a * mean(a[wh]) * (smooth(float(m),n,/EDGE) / $ (smooth(a*m,n,/EDGE)+1.e-30)) Correct me if I'm wrong here... Jaco ---------------- Jaco van Gorkom e-mail: gorkom@rijnh.nl FOM-Instituut voor Plasmafysica "Rijnhuizen", The Netherlands

**Follow-Ups**:**Re: filtering problem***From:*Dave Brennan

**References**:**filtering problem***From:*Dave Brennan

**Re: filtering problem***From:*Margaret Journey

**Re: filtering problem***From:*J.D. Smith

**Re: filtering problem***From:*Dave Brennan

- Prev by Date:
**Re: Spline Fitting and FWHM** - Next by Date:
**Re: filtering problem** - Prev by thread:
**Re: filtering problem** - Next by thread:
**Re: filtering problem** - Index(es):