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

Home » Public Forums » archive » STANDARD DEVIATON
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: STANDARD DEVIATON [message #21052 is a reply to message #20952] Fri, 04 August 2000 00:00 Go to previous messageGo to previous message
Struan Gray is currently offline  Struan Gray
Messages: 178
Registered: December 1995
Senior Member
Ben Marriage, ben@met.ed.ac.uk writes:

> I believe that this gives different results to
> the inbuilt IDL routines.

It's possible I'm missing something, but the result of my function
should be the same as yours, but multiplied by sqrt(8/9). You and IDL
are dividing by (N-1) in the definition of variance, I am dividing by
N. Statisticians love to argue about which is appropriate and where,
but for an image of the s.d. it doesn't matter much. The box smooth
always divides by N, but you can always add a renormalisation step and
it'll still be faster than nested loops.

I use this as part of a routine to do so-called statistical
differencing of images, which brings out fine detail on a varying
background (it's like an unsharp mask weighted by the local image
statistics). In this case the the s.d. image is multiplied by an
arbitrary, user-selectable factor, so the difference between N and N-1
is irrelevant.


Struan
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Dear Santa: PVWAVE Style Date/Time Axes in IDL
Next Topic: Re: object newbie

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

Current Time: Wed Oct 08 19:03:36 PDT 2025

Total time taken to generate the page: 0.00398 seconds