comp.lang.idl-pvwave archive: archive » Re: what is the best way to do a surface (or 2D) interpolation?

Home » Public Forums » archive » Re: what is the best way to do a surface (or 2D) interpolation?

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

Re: what is the best way to do a surface (or 2D) interpolation? [message #62704 is a reply to message #62702]

Tue, 23 September 2008 12:46

Vince Hradil
Messages: 574
Registered: December 1999

Senior Member

On Sep 23, 12:25 pm, Brian Larsen <balar...@gmail.com> wrote:
>> Here's a way to get verts:
>
>> sz = size(array)
>> nx = sz[0]
>> ny = sz[1]
>> nz = sz[2]
>> ns = sz[sz[0]+2]
>> verts = findgen(ns)
>> verts = transpose([ [verts mod nx], [verts/nx mod ny], [verts/nx/
>> ny] ])
>
>> BTW, I'd like to find a faster way, if there is one.
>
> This looks like the right thing but doesn't seem to give the right
> answer (or am I using it wrong?)
>
> ;; this is my data
> array=findgen(15,3)
> ;; and get the verts
> sz = size(array)
> nx = sz[0]
> ny = sz[1]
> nz = sz[2]
> ns = sz[sz[0]+2]
> verts = findgen(ns)
> verts = transpose([ [verts mod nx], [verts/nx mod ny], [verts/nx/
> ny] ])
>
> IDL> print, verts
> 0.00000 0.00000 0.00000
> 1.00000 0.500000 0.0333333
> 0.00000 1.00000 0.0666667
> 1.00000 1.50000 0.100000
> 0.00000 2.00000 0.133333
> 1.00000 2.50000 0.166667
>
> The z values that are here are not in my original array...
>
> Sorry to hijack your post Paula,
>
> Brian
>
> ------------------------------------------------------------ --------------
> Brian Larsen
> Boston University
> Center for Space Physicshttp://people.bu.edu/balarsen/Home/IDL

My example is for 3d arrays only! For 2d arrays, just use:
nx = sz[1]
ny = sz[2]
verts = lindgen(nx*ny)
verts = transpose( [ [verts mod nx], [verts/nx] ] )

There was a typo in the original, too. nx=sz[1] not sz[0], and so on.

Yes, it should be lindgen...

Report message to a moderator

[Message index]

		Re: what is the best way to do a surface (or 2D) interpolation? By: paulartcoelho on Wed, 24 September 2008 04:01
		Re: what is the best way to do a surface (or 2D) interpolation? By: Brian Larsen on Tue, 23 September 2008 13:26
		Re: what is the best way to do a surface (or 2D) interpolation? By: Vince Hradil on Tue, 23 September 2008 12:46
		Re: what is the best way to do a surface (or 2D) interpolation? By: Jean H. on Tue, 23 September 2008 10:53
		Re: what is the best way to do a surface (or 2D) interpolation? By: Brian Larsen on Tue, 23 September 2008 10:25
		Re: what is the best way to do a surface (or 2D) interpolation? By: Brian Larsen on Tue, 23 September 2008 10:15
		Re: what is the best way to do a surface (or 2D) interpolation? By: Vince Hradil on Tue, 23 September 2008 10:14
		Re: what is the best way to do a surface (or 2D) interpolation? By: Brian Larsen on Tue, 23 September 2008 09:54

Previous Topic:	Re: Getting array indices
Next Topic:	giving input to fortran program using spawn in IDL

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Thu Oct 09 17:23:00 PDT 2025

Total time taken to generate the page: 0.63643 seconds