I'm trying to use GRIDDATA, but I find it incredibly slow and running out of memory very easily.
This is the command:
grid = GRIDDATA(new_easting, new_northing, new_array, /GRID, xout=xout, yout=yout)
- new_easting is a 211*179 long array with x coordinates of my data set
- new_northing is a 211*179 long array with y coordinates of my data set
- new_array is a 211*179 integer array with the values of my data set
- xout is a 2736 long array with the x coordinates I want my data gridded and interpolated to
- yout is a 2736 long array with the y coordinates I want my data gridded and interpolated to
I'm running IDL 7.1.1, 64-bit version, using Windows 7 on a 3GHz Intel dual processor. Running the grid command causes the software to stop - IDL never moves on to the next line. If I manually stop the program, then continue it, it moves on. But not without manual intervention. Strange.
In addition, if I try to make my xout and yout ~27 million long arrays (as I originally wanted), I run out of memory. As far as I can tell, this shouldn't be the case when I'm running a 64-bit version of IDL.
Can anybody help?
Anders
|