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

*To*: jdsmith(at)astro.cornell.edu*Subject*: Re: filtering problem*From*: Dave Brennan <9147261b(at)clinmed.gla.ac.uk>*Date*: Tue, 21 Nov 2000 09:46:03 +0000*Newsgroups*: comp.lang.idl-pvwave*Organization*: University of Glasgow*References*: <3A13B236.9A321113@clinmed.gla.ac.uk> <3A158C96.7205F78B@vni.com> <3A15B24E.F0158A4@astro.cornell.edu>*Xref*: news.doit.wisc.edu comp.lang.idl-pvwave:22280

Once again thanks for the help. However, life is not quite as simple as I imagined. Previously I said : >In detail: 'the algorithm should correct the pixel value by a >multiplicative factor found by dividing the global mean by the window >mean' >A further problem is I want the ability to set a threshold where data >below the threshold are not included in the statistics and not corrected >by the algorithm. I have later found that this is not entirely correct. 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. After scratching my head for a while I can't think of a simple solution without going back to my clunky method where I manually scan through the data and threshold the window each time. There must be a better way!! Any help would be greatly appreciated!! Cheers Dave Brennan

**Follow-Ups**:**Re: filtering problem***From:*Jaco van Gorkom

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

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

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

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