comp.lang.idl-pvwave archive: archive

Home » Public Forums » archive » Re: GRIDDATA woes

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

Re: GRIDDATA woes [message #59053 is a reply to message #59044]

Mon, 03 March 2008 06:24

Kenneth P. Bowman
Messages: 585
Registered: May 2000

Senior Member

In article
<57808cc6-8454-45f1-a104-50e465ef294c@v3g2000hsc.googlegroups.com>,
"ben.bighair" <ben.bighair@gmail.com> wrote:

> I have seen a number of messages on the newsgroup about interpolation
> from an irregular grid to a regular one. None appear to address the
> issues around gridding on a sphere. I don't think I can use anything
> as simple as INTERPOLATE since the input array is sampled at irregular
> intervals.
>
> So how is this kind of interpolation supposed to be done?

If your grid is rectangular and separable (in the sense that all the
longitudes in each "column" of data are the same and all of the
latitudes in each "row" of data are same), even if the coordinates
are not regularly spaced, then it is actually quite
easy to interpolate to any set of points (regular or irregular) using
INTERPOLATE. This should be much faster than triangulating.

This problem looks just like the one David Fanning was working
on recently, and here is an outline of the solution

> Assuming that your data is 2-D (x = longitude and y = latitude), create
> the grids that you want to interpolate to
> nx = 360
> ny = 181
> x �= FINDGEN(nx)
> y �= -90.0 + FINDGEN(ny)
> Compute the "interpolation coordinates" from the original grid
> j �= VALUE_LOCATE(y_original, y)
> yj = j + (y - y_original[j])/(y_original[j+1] - y_original[j])
> Since the input and output grids are the same in the x-direction, you
> don't need to do anything with x. �Expand x and yi into 2-D arrays
> xx = REBIN(x, nx, ny, /SAMPLE)
> yy = REBIN(REFORM(yi, 1, ny), nx, ny, /SAMPLE)
> Then interpolate
> new = INTERPOLATE(original, xx, yy)

By happy chance, the interpolation chapter from my book is the sample
that is posted online here

http://idl.tamu.edu/Book.html

Ken Bowman

Report message to a moderator

[Message index]

		Re: GRIDDATA woes By: David Fanning on Tue, 04 March 2008 13:14
		Re: GRIDDATA woes By: David Fanning on Tue, 04 March 2008 11:49
		Re: GRIDDATA woes By: David Fanning on Tue, 04 March 2008 09:54
		Re: GRIDDATA woes By: wgallery on Tue, 04 March 2008 09:43
		Re: GRIDDATA woes By: David Fanning on Tue, 04 March 2008 09:14
		Re: GRIDDATA woes By: David Fanning on Tue, 04 March 2008 09:05
		Re: GRIDDATA woes By: Brian Larsen on Tue, 04 March 2008 08:23
		Re: GRIDDATA woes By: ben.bighair on Tue, 04 March 2008 08:05
		Re: GRIDDATA woes By: David Fanning on Tue, 04 March 2008 08:03
		Re: GRIDDATA woes By: David Fanning on Tue, 04 March 2008 07:35
		Re: GRIDDATA woes By: Kenneth P. Bowman on Mon, 03 March 2008 06:24
		Re: GRIDDATA woes By: wgallery on Mon, 03 March 2008 11:00
		Re: GRIDDATA woes By: Kenneth P. Bowman on Tue, 04 March 2008 19:04
		Re: GRIDDATA woes By: Kenneth P. Bowman on Tue, 04 March 2008 19:02
		Re: GRIDDATA woes By: wgallery on Tue, 04 March 2008 14:17
		Re: GRIDDATA woes By: Kenneth Bowman on Wed, 05 March 2008 06:14
		Re: GRIDDATA woes By: Kenneth P. Bowman on Tue, 04 March 2008 18:55
		Re: GRIDDATA woes By: ben.bighair on Tue, 04 March 2008 14:46
		Re: GRIDDATA woes By: pgrigis on Tue, 04 March 2008 14:39

Previous Topic:	Re: Speeding up a call to the where function
Next Topic:	GRIDDATA woes

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Oct 11 06:12:28 PDT 2025

Total time taken to generate the page: 1.20095 seconds