comp.lang.idl-pvwave archive: archive » Re: rebinning data on new time samples without loops?

Home » Public Forums » archive » Re: rebinning data on new time samples without loops?

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

Re: rebinning data on new time samples without loops? [message #37285 is a reply to message #37270]

Fri, 12 December 2003 06:31

K. Bowman
Messages: 330
Registered: May 2000

Senior Member

In article <d916d48a.0312111543.8425684@posting.google.com>,
jbattat@cfa.harvard.edu (James) wrote:

> Hi there,
>
> I have a time series of data, y, sampled on an irregular grid, t1.
> I would like to rebin y onto an existing time series (also irregular),
> t2.
>
> for loop over t1
> dT = t2 - replicate(t1[i],nt2)
> ID = where(dT EQ min(dT))
> dump y1[i] data into y2[ID] if min(dT) < dTMax
> increment counter[ID]
> endfor loop
>
> then divide y2/counter to get the rebinned value.

Since your data is already sorted by time, you can use a binary search
(VALUE_LOCATE). It is *much* faster than WHERE on long vectors.

Loop over t2 (not t1), find the points in t1 closest to t2-dt and t2+dt
using VALUE_LOCATE. Then average as you are doing.

Ken bowman

Report message to a moderator

[Message index]

		Re: rebinning data on new time samples without loops? By: jbattat on Fri, 12 December 2003 13:11
		Re: rebinning data on new time samples without loops? By: K. Bowman on Fri, 12 December 2003 06:31
		Re: rebinning data on new time samples without loops? By: Chris Lee on Fri, 12 December 2003 01:18
		Re: rebinning data on new time samples without loops? By: Chris Lee on Sat, 13 December 2003 09:41

Previous Topic:	Permanently changing shell
Next Topic:	Re: Subscripting multidimensional arrays

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sun Nov 30 15:13:29 PST 2025

Total time taken to generate the page: 0.21669 seconds