Re: filtering problem

[Dave Brennan <9147261b(at)clinmed.gla.ac.uk>]

Thanks jaco. You just hit the nail on the head. I forgot to remove the [wh]
after the smooth operation (changed a[wh] to a and couldn't work out why it
wasn't working!!!), therefore I was getting array size errors. There were a
few further lines I had to add to do exactly what I wanted but that seems
to have solved my problem. Thanks folks. Now on to the next part of the
code.......

Cheers

Dave Brennan

>
> 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