Se explica como crear un mapa gravimétrico de anomalías de Bouguer a partir de una tabla de datos ASCII del ICGEM.
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%
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
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.
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 |
gmt xyz2grd -R -G%CUT% "eigen-6c4*.gdf" -h37 -I0.095
xyx2grd -I0.095 |
No hay comentarios.:
Publicar un comentario