comp.lang.idl-pvwave archive: archive » Re: 3-dimensional integration? [IDL 5.4]

Home » Public Forums » archive » Re: 3-dimensional integration? [IDL 5.4]

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: 3-dimensional integration? [IDL 5.4] [message #35545 is a reply to message #35533]

Mon, 23 June 2003 07:04

Timm Weitkamp
Messages: 66
Registered: August 2002

Member

Oops! When I just wrote,

> Circular regions take more than one line, but are easily doable too.
> Assume that r is the radius of the circle over which you want to
> integrate (in pixels). Then you should do something like this:
>
> dim = SIZE(data, /DIMENSIONS)
> xarr = findgen(dim[0]) # (1+fltarr(dim[1]))
> yarr = findgen(dim[1]) ## (1+fltarr(dim[0]))
> rarr = SQRT((xarr-x)^2 + (yarr-y)^2)
> idxList = WHERE(rarr LE r, npix)
> IF npix NE 0 THEN tot = TOTAL(data[idxList])

... I hit the "send" key too early. The code above works, but its last two
lines can (and should) be merged into the more elegant

tot = TOTAL( (rarr LE r) * data )

I think this way of directly using a logical expression for indexing,
which in many cases I find very handy, is sometimes ignored by the
worshipers of WHERE.

Timm

--
Timm Weitkamp <http://people.web.psi.ch/weitkamp>

Report message to a moderator

[Message index]

		Re: 3-dimensional integration? [IDL 5.4] By: Timm Weitkamp on Wed, 25 June 2003 01:08
		Re: 3-dimensional integration? [IDL 5.4] By: Ingo Salzmann on Tue, 24 June 2003 07:18
		Re: 3-dimensional integration? [IDL 5.4] By: Timm Weitkamp on Mon, 23 June 2003 07:04
		Re: 3-dimensional integration? [IDL 5.4] By: Timm Weitkamp on Mon, 23 June 2003 06:36

Previous Topic:	Re: sec : U Re: Sinc interpolation
Next Topic:	Re: How to do nested objects??

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Oct 08 18:59:43 PDT 2025

Total time taken to generate the page: 0.00393 seconds