Nota

Este proyecto solo se incluye en el disco de máquina virtual OSGeoLive (VMDK)

Inicio Rápido GMT¶

GMT es una colección de herramientas que permiten a los usuarios manipular conjuntos de datos (x, y) y (x, y, z) (incluyendo el filtrado, ajuste de tendencia, retícula, proyección, etc.) y producir ilustraciones Archivo de PostScript Encapsulado (EPS) que van desde simples diagramas x-y a través de mapas de contorno hasta superficies iluminadas artificialmente y vistas en perspectiva 3D en blanco y negro, tonos grises, patrones ashurados y colores de 24 bits.

En este rápido tutorial GMT utilizamos un modelo digital del terreno (DTM) para generar un mapa de relieve sombreado.

Contenidos

Crear un mapa
Crear una paleta de colores
Definir el degradado
Mostrar la barra de color
Cosas para probar
¿Qué sigue?

Crear un mapa ¶

Como nuestro conjunto de datos de entrada utilizaremos un DTM del conjunto de datos de ejemplo GRASS GIS, al que GDAL puede acceder mediante el controlador GDAL-GRASS.

Desde un vistazo rápido a los metadatos DTM, abra una terminal y utilice el comando GDAL gdalinfo:

gdalinfo /home/user/grassdata/nc_basic_spm_grass7/PERMANENT/cellhd/elevation

Podemos ver que el DTM está en un sistema de coordenadas proyectado y el rango de datos es entre 50 y 160 metros.

Podemos convertir el DTM en un formato de archivo GMT netCDF compatible con gdal_translate como se muestra a continuación:

gdal_translate -of GMT /home/user/grassdata/nc_basic_spm_grass7/PERMANENT/cellhd/elevation elevation.gmt

Para nuestro ejemplo GMT primero transformamos el conjunto de datos en coordenadas geográficas longitud-latitud (EPSG:4326) utilizando el comando gdalwarp :

gdalwarp elevation.gmt geo_elevation.gmt -t_srs EPSG:4326

Crear una paleta de colores ¶

Para asignar un mapa de color a nuestra red podemos hacer una paleta de colores GMT tablas con el comando makecpt. Estableceremos el rango y el número de colores con la bandera -T y asignar un mapa de colores con la bandera -C.

gmt makecpt -Chaxby -T50/160/10 -Z > elevation.cpt

Puedes experimentar con otras tablas de color, vea man gmt makecpt para una lista completa.

Definir el degradado ¶

Para ejemplificar el relieve sombreado tenemos que especificar cómo se distribuye la intensidad del color a lo largo de la cuadrícula. Para ello, calculamos la direccional derivada, o gradiente, de nuestra red utilizando el comando grdgradient. En grdgradient usaremos el -Ne bandera para normalizar la salida usando una distribución de Laplace acumulativa y la -A bandera para la dirección Aimutal para un derivado direccional.

gmt grdgradient geo_elevation.gmt -Ne0.8 -A100 -fg -Ggradient.nc

El resultado es un archivo netCDF que luego puede ser utilizado como la intensidad de color mediante el comando grdimage como sigue:

gmt grdimage geo_elevation.gmt -Igradient.nc -JM6i -P -Ba -Celevation.cpt -K > GMT_tut_OSGeoLive.ps

Mostrar la barra de color ¶

Por último podemos añadir una colorbar con psscale:

gmt psscale -DjTC+w5i/0.25i+h+o0/-1i -Rgeo_elevation.gmt -J -Celevation.cpt -I0.4 -By+lm -O >> GMT_tut_OSGeoLive.ps

Aquí la bandera -D se utiliza para seleccionar la posición de la barra de color TC: Top Center, el tamaño 5x0.25 se expresan en pulgadas i la orientación por h: horizontal y la bandera -By para establecer una etiqueta m en el y axis.

El documento de salida ps debe parecerse: