[GRASSLIST:6589] Re: Who wants contour lines labeled "300.000000"?

[Jose Gomez-Dans]

On 4/21/05, Dylan Beaudette <dylan at iici.no-ip.org> wrote:
> On Wednesday 20 April 2005 09:48 pm, Gordon Keith wrote:
> > On Thu, 21 Apr 2005 08:19, Dylan Beaudette wrote:
> > > I have found GMT to be an excellent solution for creating cartographic
> > > products directly based on GRASS rasters.
> >
> > Would you have any pointers on going the other way?
> >
> > I have some GMT grid files I would like to analyse in GRASS but haven't
> > been able to find a way to read them in yet. I've tried both r.in.bin for
> > type 1 grids and r.in.gdal and both create rasters with all values NaN.
> >
> > Regards
> > Gordon
> 
> 
> Hmm.. I haven't tried going the other way: from GMT-> GRASS ... r.in.bin
> would 
> be my initial guess, as the NetCDF support in r.in.gdal has never worked for
> 
> me. 

You can do that easily enough. First, note that you want grid
registered GMT grids. You can use grdsample with the -T option if
that's not the case. Then, convert the grids to a binary file using
grd2xyz -ZTLf <grid_file> > grid.bin, and read them into GRASS with
r.in.bin. This works fine for me here, but I don't know whether GDAL
might be an option. My method is a bit long-winded, but "it works for
me" :)

Greetings,
Jose