Importing interpolated data

GRASS GRASS user
1

From: Henk Coetzee <henkc@geoscience.org.za>
Subject: Importing interpolated data

I am looking for ways of importing already interpolated data from
Generic Mapping Tools (gmt) which uses the NetCDF data format. Grass's
interpolation tools are good, but I find them extremely slow, compared
to the tools from gmt. The nature of the gmt algorithms is such that I
could produce perfectly matching regions. Is there a way to do this? An
alternative would be to write the data back to x,y,z values from the
grid, bring it in as a site file and write that to a raster. Any ideas??

How about using the GMT command "grdreformat" to create a Sun standard
8-bit rasterfile and then trying to import that instead of the grd2xyz
approach you mention above. (GMT-3.3.1 is current, but I think grdreformat
goes back to GMT-3.1.)

Ben Horner-Johnson
ben@earth.nwu.edu

2

Hi

I am working on streamlining the importation of data from GMT. There are
a number of options, all of which seem to work well. The problem that I
am finding is that I have raster files with a large number of NaN
values. These seem to convert to zeros, which is not what I want. Any
suggestions?

Thanks

Henk

------------------------------------------------
Henk Coetzee
Council for Geoscience
Private Bag X112, Pretoria, 0001, South Africa
Tel: +27-12-841-1192 Fax: +27-12-841-1424
------------------------------------------------