[GRASS-user] OT: importing .grd files (GEBCO)

[antonio rodriguez]

Hi,

I've recently faced this issue and with the help of some members of the 
list I got a solution to do this. My problem began when I tried to 
import some bathymetric data from the GEBCO initiative . This data is 
delivered in .grd format (GMT) Below I describe the methodology to 
accomplish the task (First in english, then in spanish)

Debian/Linux. GRASS 6.3. GMT package

1) In the console with the command grd2xyz we convert the .grd file to 
an .xyz one:

    grd2xyz file.grd > file.xyz

2) With the help of 'awk' we delete columns 1 & 2 (corresponding to 
lat/lon) and leave column '3' (bathimetry/altitude):

    cat file.xyz | awk '{print $3}' > baty.xyz

3) With the help of grdinfo we search for the dimensions and resolution 
of the .grd file

4) Knowing this we introduce a header in the baty.xyz file, i.e:
   
    north:   38
    south:   34
    east:    -2
    west:    -9
    rows:    241
    cols:    421

5) Now in GRASS, we create a new location with the above data. Datum and 
Projecion will depend on your area.

6) Then:

    r.in.ascii in=baty.xyz out=bathimetry

Hope this helps,

Antonio


------------------------------------------------------------------------------------------------------------------------------

1) Con la herramienta GMT grd2xyz convertimos el archivo .grd en .xyz

2) Mediante 'awk' eliminamos las columnas (1 y 2) correspondientes a los 
datos de lat y lon

    cat archivo.xyz | awk '{print $3}' > baty.xyz

3) Averiguamos las dimensiones (x e y) del área en cuestión (grdinfo 
archivo.grd)

4) Como sabemos los valores de lat y lon (grdinfo archivo.grd)los 
introducimos como cabecera en bati.xyz junto con las columnas (x) y las 
filas (y), tal que:
   north:   38
   south:   34
   east:    -2
   west:    -9
   rows:    241
   cols:    421

5) En GRASS nos creamos una nueva localidad (bat) con coordenadas 
lat(38n:34n)/lon(9w:2w), datum (eur50), proj = 5. Resolución (de grdinfo 
archivo.grd)

6) Con r.in.ascii importamos el archivo .xyz:

    r.in.ascii in=baty.xyz out=batimetria

Saludos,

Antonio