En este tutorial vamos a explicar paso a paso:
- Cómo trabajar con archivos .asc en QGIS.
- Cómo convertir en otros formatos.
- Cómo crear curvas de nivel y sombreados
- Cómo crear visualizaciones en 3D.
- Por último te enseñaremos cómo crear un perfil topográfico en QGIS.
QGIS tiene fantásticas herramientas y complementos para el procesamiento de datos del terreno o elevaciones y para su visualización. Ante la pregunta:
Tengo un archivo .asc y quiero importarlo a QGIS. ¿Cómo lo hago?
La respuesta es que no necesitas realizar ningún tipo de conversión previa. QGIS lee perfectamente los archivos .asc (Arc/Info ASCII Grid) gracias al uso de la librería GDAL para leer y escribir formatos de datos raster, incluyendo Arc/Info Binary Grid, GeoTIFF, Erdas Imagine y muchos más. La librería GDAL soporta más de 100 formatos raster.
Índice
Obtener datos del terreno
Desde el Centro de Descargas del IGN podemos descargar un MDT200 (Modelo digital del terreno con paso de malla de 200 m.) de la provincia de León. Utilizaremos este ráster para realizar los ejercicios planteados en este tutorial.
Para poder realizar la descarga es necesario registrarse en el sitio. Una vez realizada la descarga tendremos un archivo llamado PNOA-MDT200-ETRS89-LEON-H30.asc
Añadiendo datos raster en QGIS
Podemos añadir el ráster arrastrando el archivo directamente a QGIS o a través del menú superior Capa → Añadir capa raster…
Una vez añadida la capa veremos los datos del terreno en QGIS. Cada píxel del raster representa la elevación media en metros en ese lugar. Los píxeles oscuros representan las zonas de baja altitud y píxeles más claros representan áreas con alta altitud.
Convertir a otro formato raster
Para convertir un archivo raster a otro formato raster en QGIS, debemos ir al menú superior Ráster > Conversión > Traducir (Convertir formato)
Haciendo uso de la librería GDAL, QGIS es capaz de convertir prácticamente cualquier formato en otro, la lista es muy amplia: .SID, . TIFF, . IMG, HDR.ADF, .DT0, .DT1, .DT2, . JPG, . PNG, .BMP, .TER, .ASC, .SQLITE, .E00, .MBTILES, etc. Consulta la lista de formatos raster de GDAL.
Extraer las curvas de nivel
Para extraer las curvas de nivel del modelo digital del terreno seleccionamos la herramienta Curvas de nivel del menú superior Ráster > Extracción > Curvas de nivel.
- En el cuadro de diálogo Curvas de nivel, seleccionamos MDT200-LEON-H30 como archivo de entrada.
- Como nombre del archivo de salida de las curvas de nivel mdt200_leon_curvas100.shp.
- Vamos a generar las curvas de nivel en intervalos de 100 metros, por lo que ponemos 100 como el intervalo entre las curvas de nivel.
Clic en Aceptar.
Una vez finalizado el proceso, verá las curvas de nivel cargadas en QGIS. Cada línea de esta capa representa una elevación. Todos los puntos a lo largo de estas líneas estarían a la misma altura. Y cuanto más próximas estén las líneas más pronunciada será la pendiente.
Crear un hillshade (mapa de sombras)
Ahora vamos a crear un sombreado del relieve del ráster en QGIS. Seleccionamos Ráster > Análisis > Mapa de sobras (Hillshade).
Como puedes ver en la imagen anterior, la herramienta Raster > Análisis también nos permite crear mapas de:
- Pendientes y orientación (tutorial aquí).
- Irregularidad (escabrosidad).
La creación de un hillshade es muy sencillo, en la ventana Mapa de sombras (Hillshade), seleccionamos MDT200-LEON-H30 como archivo de entrada y como nombre del archivo de salida del sombreado MDT200-LEON_hillshade:
Clic en Aceptar.
Jugando con la transparencia del hillshade y cambiando la rampa de colores del MDT original este puede ser el resultado:
Mostrando datos del terreno en 3D
En QGIS 3 ya podemos crear fácilmente visualizaciones en 3D de nuestros datos geográficos.
Simplemente desde el menú superior Ver > Nueva vista de mapa 3D.
La configuración 3D podemos modificar la escala vertical para exagerar la elevación y hacer la visualización más atractiva:
Crear perfiles topográficos con QGIS
Para crear perfiles topográficos en QGIS necesitamos tenemos dos opciones: utilizar la herramienta integrada en QGIS o utilizar algún complemento como como Profile Tool. En esta entrada te mostramos cómo hacerlo de las dos formas:
La herramienta integrada
El complemento Profile Tool
Cómo utilizar la calculadora ráster en QGIS
La calculadora ráster es una herramienta que permite realizar operaciones matemáticas sobre los valores de los píxeles existentes en un ráster, lo que es útil para la conversión y manipulación de dichos datos.
Podemos dividir la calculadora ráster en tres partes:
- Bandas ráster. Contiene todas las capas ráster cargadas en la interfaz de QGIS y que pueden ser utilizadas para la realización de cálculos sobre las mismas.
- Capa resultado. Lugar donde definimos la capa de salida. Además, podemos indicar la extensión de la zona de cálculo para establecer la resolución de la capa de salida.
- Contiene todos los operadores disponibles para añadir a la caja de expresiones de la herramienta. Cuenta con operadores matemáticos y funciones trigonométricas.
Los resultados de dichas operaciones se muestran en una nueva capa ráster en un formato admitido por la librería GDAL.
En este tutorial te mostramos cómo modificar los valores de un ráster a partir de una máscara con la calculadora ráster de QGIS.
Let’s connect!
Date de alta en nuestra newsletter y te enviaremos GRATIS el ebook que te ayudará a impulsar tu perfil GIS:
Vitaminas MappingGIS
Tan solo una vez al mes recibirás las últimas novedades del sector GIS y de nuestros cursos
Licenciado en Geografía. Máster en Sistemas de Información Geográfica. Consultor GIS desde el año 2004. En MappingGIS desde el año 2012 para ayudarte a impulsar tu perfil GIS y diferenciarte de la competencia. Echa un vistazo a todos nuestros cursos de SIG online.
Buenas tardes Aurelio:
Estoy intentando hacer un estudio y no lo consigo, te comento….
Tengo una serie de puntoS en una capa .shp en concreto estaciones de medición, por otro lado una capa raster con la altura z.
Me gustaría saber como vincular la altura del raster con la capa shp y obtener la altura a la que se encuentra capa punto extraída del raster.
Gracias y enhorabuena por los cursos y la divulgación en GIS
Hola RAfael,
Gracias por tu comentario.
Para resolver tu problema hay varias opciones: desde utilizar herramientas de SAGA como «Add grid values to points» o utilizar el plugin de QGIS «Point Sampling Tool», hasta convertir el raster a vectorial y realizar una intersección posterior entre las dos capas. Saludos!
En QGIS ¿se pueden hacer curvas de nivel maestras y secundarias?
Saludos
Hola Alberto,
En QGIS puedes generar todas las curvas de nivel, después puedes crear una columna nueva en la tabla de atributos para indicar qué curvas quieres que sean maestras y cuales secundarias. Con este paso previo ya puedes realizar una simbología categorizada, para representar las curvas maestras de un color y las secundarias de otro. Un saludo!
Buen día Amigo. Una duda para obtener las Pendientes: Cuando genero mi archivo (en formato .tif), luego de que quiero dar formato a éste clasificando las pendientes que se supone serán generadas, éste me indica que el único dato obtenido es de 89.9 (tanto min. como max.), por lo que no puedo hacer ninguna clasificación y mi mapa aparece de un solo color. Estaré haciendo algún proceso mal?
Hola Alonso, Me temo que el TIFF que tienes no se ha generado bien ya que los píxeles no contienen información. Prueba a generar de nuevo un TIFF con valores en sus celdillas. Saludos
Como descargo un Dem de la zona
Hola Duban,
Puedes descargar el DEM desde el Centro de Descargas del IGN o visitar esta página: 10 Fuentes de datos GIS gratis: Los mejores datos globales raster y vectoriales. Saludos!
hola a todos,
alguien sabe como se puede imprimir a escalas las secciones topográficas realizadas con las herramientas de qgis (qprof, VoGis profile tool). Necesito para realizar secciones geológicas.
gracias
Hola Aurelio, estoy empezando a usar Q GIS y necesito corregir un DEM que debido a la presencia de vegetación me impide realizar buenos análisis hidrológicos, vos que me recomendarias?. Por ultimo me gustaría poder utilizar la herramienta para escabado de cauces de SAGA pero me dice el programa que las capas no tienen la misma extensión de cuadricula, a que hace referencia? Desde ya muchas gracias!
Hola Ramiro, Como paso previo a los análisis hidrológicos será necesario corregir el DEM. Hay que realizar un quemado (burn) con una capa de ríos para indicar por donde irá el curso de agua. Puedes utilizar SAGA. Lo que debes hacer es tener amabas capas, la raster y la vectorial en el mismo sistema de referencia espacial.
Un saludo!
Hola,
Lo primero, enhorabuena por la web y gracias por compartir tanta información. Estoy utilizando Qgis y he realizado un mapa de relieve, a partir de un MDT y sus curvas de nivel, y he modificado la rampa de colores por tramos de altitud (de 0 a 100, 100 a 200 m., etc) para que quede lo más visual posible.
Mi duda es si hay alguna forma de poder mostrar dichos tramos por separado, y así poder mostrar sólo los tramos de altitud que interesen, en vez de todo en conjunto asociado al ráster de origen, como si fueran capas independientes. Lo he intentado recortando el ráster, pero necesito un shape poligonal para poder hacer el recorte, la capa de curvas de nivel es lineal y al intentar transformarla a poligonal, me da error.
Quizá sea alguna opción desde el mismo archivo ráster que permita visualizar por separado, pero de momento no doy con ello.
Gracias de antemano y un saludo.
Hola Antonio,
Muchas gracias por el comentario!
Creo que lo que necesitas es reclasificar el raster, tienes un algoritmo de GDAL que se denomina r.reclass que te permite realizarlo. Te creará una capa nueva con los valores categorizados, así te permitirá hacer filtros… Un saludo!
Gracias Aurelio, lo que pasa es que soy nuevo en qgis y no he trabajado con GDAL ni algoritmos, sabes o alguien conoce de algún tutorial para que pueda ver los pasos a seguir?
gracias de nuevo
Hola Antonio, No quiero hacer promoción, pero quizá te venga bien realizar el curso que tenemos de QGIS inicial-intermedio: http://mappinggis.com/cursos/qgis/ En el que se trabaja con algoritmos y tendrás a un tutor a tu disposición para resolver cualquier duda. Saludos!
No está mal, le echaré un vistazo. Al final lo he resuelto con la calculadora de ráster, y combinando los resultantes posteriores.
saludos!
Podrían decir como suavizar curvas de nivel, ya que al realizarlas salen con muchas esquinas en qgis, gracias un saludo
Hola Carlos,
Para suavizar curvas de nivel lo mejor es utilizar el algoritmo de GRASS > v.generalizar. Y elegir el método chaiken. Pruébalo y nos cuentas. Un saludo
Hola,
Yo soy estudiante senegales, uso Qgis en mis investigaciones y Me gustaria saber si hay una posibilidad de crear un MDT a partir de datos X, Y Z que tengo en una BDD Excel o .txt
Gracias por ayudarme
Hola,
No soy capaz de encontrar una herramienta que me de el siguiente resultado:
Tomando como base un MDT y dos puntos de coordenadas, necesitaría una herramienta que me diera la línea con mínima pendiente entre ambas coordenadas. Todo ello enfocado al diseño de senderos y caminos en orografía acusada.
Alguien podría orientarme?
Gracias y un saludo.
PD:Soy usuario tanto de ArcGIS como QGIS y Global Mapper
Hola, Aurelio. Una duda respecto a QGIS y su Profile Tools:
Como bien has indicado en tu entrada,
una vez creado el perfil longitudinal, podemos exportarlo a varios formatos ( PDF, PNG, SVG o PS ), pero ¿y si lo que nos interesa es obtener un listado, CSV por ejemplo, con la distancia y altitud Z del perfil?
El objetivo de todo esto es obtener un listado que me aporte la distancia (abscisa) y la altitud (ordenada) de una serie de senderos, partiendo de lo que tenemos: el track de los senderos y un MDE. Cómo obtener la representación de ese perfil está claro, como hemos visto y tú mismo comentas en tu entrada, pero lo que no logro resolver es cómo obtener desde QGIS ese listado en formato texto para poder manejarlo en otras aplicaciones.
Lo tengo resuelto vía aplicaciones externas, por ejemplo con TCX Converter, que me permite a partir del track calcular la Z y guardar el resultado en un fichero tipo csv. Pero lo que me gustaria es hacerlo desde el propio GIS, QGIS, gvSIG o ArcGIS.
Por tanto, la pregunta es: ¿Conoces alguna forma de obtener desde QGIS, gvSIG o ArcGIS el listado Dis/Z (Distancia/Altitud) de un eje, partiendo del TRACK (XY) y de un MDE (que es lo que nos aporta la Z)?
Gracias y un saludo
Resuelto: Con el algoritmos de SAGA «Profiles from lines», accesible desde la Caja de Herramientas de Procesado de QGIS, se puede obtener la salida tabulada que buscaba teniendo como datos de partida el TRACK y el MDE
Hola Francisco, muchas gracias por compartir con todos cómo lo has resuelto, la verdad es que yo no sabía como resolver este «problema».
Eres un fenómeno!