Posteado por: rengifo | enero 8, 2016

Como crear un mapa de relieve coloreado usando Gdal y GRASS GIS


Estimados lectores,

Ante todo espero que hayan tenido un felíz año nuevo. En esta oportunidad les escribo brevemente para demostrar como crear un mapa de relieve coloreado de acuerdo a la altura. Básicamente se trata crear un RGB con relieve.

Lo primero que necesitan es tener instalado el GRASS GIS 6.4 que es  la version que uso en este caso , y la libreria  GDAL.También va a necesitar un MED o tif que contenga los valores de altura.

El  primero paso es que crear una escala que comprenda los nivel de altura  observados en el MED/tif.

En este caso usé un MED de 25 metros del Archipielago de Svalbard, archipielago localizado entre Noruega y el polo Norte.

De acuerdo al MED  de Svalvard el mismo tiene alturas máximas de  1700 msnm (sobre el nivel del mar). Asi que usando  notepad u otro editor de texto simple  cree una tabla de colores con los niveles de altura y los codigo de los colores deseados en código RGB. La misma la guardé   en formato txt. (tabla_colores.txt) en el mismo lugar donde tenía el MED/tif.

0 110 220 110
700 240 250 160
1200 230 220 170
1500 220 220 220
1700 250 250 250

Una vez hecho esto  abrí gdal en la consola de comando del dos (cmd) y navegé  al lugar donde tenía el  MED og Tif.

Dentro de la línea de comando escribí el siguiente  comando:

gdaldem color-relief svalvard.tif tabla_colores.txt svalvard_color_relieve.tif

Vamos descomponer el comando para que sea más fácil entenderlo:

“gdaldem color-relief”: es una utilidad  de  gdal y dice que tiene tomar el svalvard.tif(input), ver la altura y asignarle los colores definido en la tabla_colores.txt y  guardarlo en un nuevo tif llamado svalvard_color_relieve.tif.

Una vez hecho  esto pueden ver la imágen y ver si están satisfechos con la distribución de los colores con respecto a la altura. De otra manera pueden modificar  la tabla y correr el comando de nuevo el tif anterior será resescrito. (Ver imágen 1)

Posteriormente se importa  svalvard.tif(input) y svalvard_color_relieve.tif  en GRASS GIS. Al importar svalvard_color_relieve.tif  GRASS GIS descompone el .tif en  RGB (rojo, verde y azul). Usando  el comando r.composite en GRASS GIS  se combina los 3 canales  y se restituye como la imágen producida por gdal anteriormente. (ver imágen 1 )

 

 

sv

Imagen 1: Muestra el tif original a la izuierda y  el tif(RGB) creado a con  GDAL a la derecha.

Con el  .tif original se genera un relieve de sombras (shade) usando el comando r.relief y se  fusiona con el RGB usando el comando r.shade
El resultado de la fusion  se exporta usando el comando r.out.tiff (ver imagenes 2,3,4,5,6)

2016-01-08 10_44_10-

Imágen 2 : Ejecutando el comando r.relies desd GRASS GIS

 

2016-01-08 10_44_58-GRASS GIS 7.0.0 Map Display_ 1  - Location_ svalvard@PERMANENT

Imágen 3: Fusionando el RGB y el hillshade con el comando r.shade

 

 

2016-01-08 09_02_21-Untitled - ArcMap

Imágen  4: Exportando el   el relieve coloreado usando el comando r.out.tiff en GRASS GIS 6.4

 

2016-01-08 09_02_33-Untitled - ArcMap

Imágen 5: definiendo las opciones de exportacion, es importante * crear un world file de otra manera el tif pierde la proyección .

 

2016-01-08 12_44_24-GRASS GIS 6.4.4 Map Display_ 1  - Location_ svalvard

Imágen 6: Mapa de relieve coloreado

Es  importante  remarcar que en el caso de Svalvard  habria que ajustar los colores para reflejar la existencia de los glaciales, como se muestra en el mapa de GEBCO (ver imagen 7).

 

2016-01-08 11_13_25-QGIS 2.8.1-Wien

Imagen 7: Extraida del  servicio de mapas en linea (WMS) de GEBCO.

Espero  haberles despertado la curiosidad respecto al uso combinado de GDAL y GRASS GIS. Para personas que dominen el inglés se encuentra algo parecido en en este blog.

Me despido cordialmente hasta la proxima vez.

Ps: disculpa por la falta de acentos.. por alguna razón no funciona en wordpress hoy.

 

 

 

 

 

 

 

 

Anuncios

Responder

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s

Categorías

A %d blogueros les gusta esto: