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

Home » Public Forums » archive » Re: efficient kernel or masking algorithm ?
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: efficient kernel or masking algorithm ? [message #22643 is a reply to message #22642] Wed, 29 November 2000 00:00 Go to previous messageGo to previous message
Richard Tyc is currently offline  Richard Tyc
Messages: 69
Registered: June 1999
Member
WOW, I need to look at these equations over about a dozen times to see what
is going on ?

I have been struggling with the variance of an nxn window of data, INCLUDING
central pixel

;mean of the neighboring pixels (including central)
mean=smooth(arr,n)
;square deviation from that mean
sqdev=(arr-mean)^2
;variance of an nxn window of data, INCLUDING central pixel
var=(smooth(sqdev,n)*n^2-sqdev)/(n^2-1)

This doesn't seem correct with test samples ? (Only difference is mean and
division by n^2-1 ??)

Thanks JD

Rich


J.D. Smith <jdsmith@astro.cornell.edu> wrote in message
news:3A25758E.A83B10CA@astro.cornell.edu...
>
> Oh my this is a common topic lately. See my recent posts in a thread
> with title "Array Manipulations". Here's the good stuff:
>
> ; the nxn window total
> total=smooth(arr,n)*n^2
> ; the nxn window total not including central pixel
> neighbors=smooth(arr,n)*n^2-arr
> ; the mean of the neighboring pixels (excluding central)
> neighmean=(smooth(arr,n)*n^2-arr)/(n^2-1)
> ; the square deviation from that mean
> sqdev=(arr-neighmean)^2
> ; the variance of an nxn window of data, excluding central pixel
> imvar=(smooth(sqdev,n)*n^2-sqdev)/(n^2-2)
>
> Take a look at the "EDGE*" keywords too, if you care about what happens
> near the borders.
>
> JD
>
> --
> J.D. Smith | WORK: (607) 255-6263
> Cornell Dept. of Astronomy | (607) 255-5842
> 304 Space Sciences Bldg. | FAX: (607) 255-5875
> Ithaca, NY 14853 |
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Previous Topic: Re: Contours driving me mad
Next Topic: How Computers Represent Floats

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

Current Time: Sat Oct 11 15:08:47 PDT 2025

Total time taken to generate the page: 1.60187 seconds