| what is the best way to do a surface (or 2D) interpolation? [message #62714] |
Tue, 23 September 2008 03:04 |
paulartcoelho
Messages: 30
Registered: March 2007
|
Member |
|
|
hi,
i have a set of 2D data, say:
IDL> help,vxjc
VXJC FLOAT = Array[6, 16]
IDL> print,vxjc
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
21.2471 0.00000 0.00000 9.55912 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
38.7544
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 0.00000
0.00000
0.00000 0.00000 0.00000 0.00000 30.4393
0.00000
it's really like this, full of zeros and some few bins with non-zero
valuers. i know what are the x- and y- coordinates for this data, say:
for columns:
IDL> print,vz
-1.62839 -1.23045 -0.628389 -0.327359 0.0483046
0.246672
for rows:
IDL> print,vage
1.00000e+08 1.40000e+08 2.00000e+08 2.80000e+08 4.00000e+08
5.60000e+08 7.90000e+08 1.12000e+09 1.58000e+09 2.24000e+09
3.16000e+09 4.47000e+09 6.31000e+09 8.91000e+09 1.25900e+10
1.58500e+10
but i need to represent this surface at some other coordinates, say
that i need x-axis to have values
newx = [-2.0, -1.5, -1.0, -0.5, 0.0, 0.5]
and y-axis a vector with 10 elements, spaced regularly in log.
what is the best way to rebin and/or interpolate my original 2D data
to the desired new coordinates?
i've been trying without success with the gridding/interpolation
routines, that i'm completely new at.
many thanks for any help,
paula
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
