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

Home » Public Forums » archive » Image warping in IDL
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: Image warping in IDL [message #51405 is a reply to message #51319] Mon, 20 November 2006 12:17 Go to previous messageGo to previous message
JD Smith is currently offline  JD Smith
Messages: 850
Registered: December 1999
Senior Member
On Tue, 14 Nov 2006 11:20:38 +0100, Wox wrote:

> On Fri, 10 Nov 2006 15:51:34 -0700, JD Smith <jdsmith@as.arizona.edu>
> wrote:
>
> <snip>
>
>> If you code this up, let us know whether it was faster.
>
> Wow, thanks! A nice opportunity to brush up my histogram skills :-).
>
> I'm almost there and I will post some benchmarks when finished.
>
> One question though: "Don't you need a second loop (a loop over the
> repeat counts in H4)?"

Yes, as I see you figured out. Nice implementation. As you found,
explicitly loop from 1 to omax in your histogram of repeat counts is
fine, and solves the problem without any monkeying of indices. In
fact, the snippet j=1,omax4-omin4 works only when omin4 is zero (which
it seems to be always for you). j=1,omax should work. If you want to
handle the j=1 case seperately for efficiency (as you've done), just do so
and start the loop at 2. Also, I couldn't quite understand the
rebin([3,2,3,2],4,npix) for selecting which 4 of the 9 output pixels
actually receive any data. It seems like those are fixed offsets, which
wouldn't work when the offset direction rotates around. Maybe something
about your mapping lets you get away with that.

JD
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Making movies
Next Topic: Re: Pair Counts in an Annulus, for large data sets

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

Current Time: Wed Oct 08 14:32:39 PDT 2025

Total time taken to generate the page: 0.00413 seconds