7.3. Mapa Gravimétrico. Anomalía de Bouguer.

Se explica como crear un mapa gravimétrico de anomalías de Bouguer a partir de una tabla de datos ASCII del ICGEM.

El script permite crear un mapa de anomalías de Bouguer a partir de la tabla de datos del modelo Eigen-6c4 descargados de la página del ICGEM (International Centre for Global Earth Model).

0. Descargar datos: Para este ejemplo utilizo la tabla de datos "eigen-64c-76909.gdf". El archivo tiene un encabezado de 37 líneas (hasta "end of head") y luego una tabla de datos en 3 columnas (longitud, latitud y anomalía de Bouguer). Ver Anexo 3 para más información.

1. Crear grilla: Para crear una grilla a partir de una tabla de datos ASCII uniformemente distribuidos (eigen-64c-76909.gdf) se utiliza xyz2grd. Con -G se indica que se utilice el nombre definido por la variable %CUT% para guardar la grilla. -I indica el espaciado entre datos ("Grid Step") definido al descargar la grilla (ver Anexo 3) y escrito en la línea 24 del archivo. -R define la región geográfica de la grilla (debe ser igual o menor a la abarcada por los datos).

Nota: Para datos distribuidos aleatoriamente es necesario primero utilizar algún módulo que los distribuya uniformemente (por ej.: surface,  triangulate, nearneighbor).

gmt xyz2grd -R -G%CUT% "eigen-6c4*.gdf" -h37 -I0.1

2. Hacer Mapa: En los siguientes pasos se utilizan las herramientas makecpt, grd2cpt, grdgradient, grdimage, psscale, etc explicados previamente (pasos 5 a 8 en entradas 5.1 y 7.1). El uso de makecpt y grd2cpt permiten ajustar l distinta

gmt makecpt -T-550/450/50 -Z >%color%
gmt grd2cpt %CUT% -Z >%color%
gmt grd2cpt %CUT% -Z -L-150/150 -D >%color%

makecpt -T-550/450/50

grd2cpt %CUT% -Z

grd2cpt %CUT% -Z -L-150/150 -D

3. Resolución: En caso de elegir una resolución de grilla (Grid Step) menor a la que corresponde (calculada en la página del ICGEM), se observarán líneas grises (que corresponden a zonas sin datos) en el mapa.

gmt xyz2grd -R -G%CUT% "eigen-6c4*.gdf" -h37 -I0.095

xyx2grd -I0.095