comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » filtering problem
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: filtering problem [message #22621 is a reply to message #22539] Tue, 21 November 2000 00:00 Go to previous messageGo to previous message
Jaco van Gorkom is currently offline  Jaco van Gorkom
Messages: 97
Registered: November 2000
Member
"Dave Brennan" <9147261b@clinmed.gla.ac.uk> wrote in message
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
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Spline Fitting and FWHM
Next Topic: Getting BIG structures into IDL via CALL_EXTERNAL

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 10 09:41:21 PDT 2025

Total time taken to generate the page: 1.27984 seconds