Análisis Avanzado de Ráster (QGIS3)

En el tutorial anterior Estilos y Análisis Básicos Ráster (QGIS3), aprendió a realizar álgebra de ráster con la Calculadora de raster. Este tutorial se basa en estas técnicas y le muestra cómo utilizar otras herramientas de análisis ráster de la caja de herramientas de procesamiento. Aprenderá a procesar con los rásteres de Uso y Cobertura de la Tierra (UCT) en QGIS para extraer ciertos tipos de clases de cobertura del suelo y mapear los cambios.

Vista general de la tarea

Utilizaremos el conjunto de datos South African National Land Cover para identificar y extraer los asentamientos informales en la ciudad de Johannesburgo (Sudáfrica). También utilizaremos un conjunto de datos de evaluación de cambios para identificar los patrones de crecimiento urbano en la ciudad desde 2014 hasta 2018.

Otras habilidades que aprenderá

  • Cómo reproyectar datos ráster a otra proyección.

  • Cómo cargar un archivo Excel en QGIS.

  • Cómo ajustar la transparencia de una capa ráster en QGIS.

Obtener los datos

Para este tutorial descargaremos los siguientes conjuntos de datos

  1. El conjunto de datos South African National Land Cover 2018: El nuevo conjunto de datos South African National Land-Cover 2018 se ha generado a partir de imágenes satelitales Sentinel 2 multitemporales de 20 metros y contiene 73 clases de cobertura de la tierra.

  2. Evaluaciones de cambio de la cobertura de la tierra nacional de Sudáfrica 2018: Este conjunto de datos compara el cambio en 20 clases de cobertura de la tierra desde 2014 hasta 2018.

  3. Límite de la ciudad de Johannesburgo: Un archivo shape de los límites de la ciudad de Johannesburgo, Sudáfrica.

El portal Environmental Geographical Information Systems (E-GIS) proporciona acceso a datos geoespaciales medioambientales de Sudáfrica. De este portal descargaremos los rásteres de South African National Land-Cover (SANLC).

  1. Visite la página de Descarga de datos GIS. Haga clic en Estoy de acuerdo para aceptar las Condiciones de uso y continuar.

../../_images/data1.png
  1. Deberá crear una cuenta gratuita para descargar el conjunto de datos. Haga clic en Deseo crear una cuenta y siga las instrucciones para crear una cuenta.

../../_images/data2.png
  1. Una vez iniciada la sesión, busque South African National Land Cover (SANLC) 2018 Computer Automated Land Cover (CALC). Este conjunto de datos se ofrece en dos proyecciones diferentes. Para este tutorial, descargaremos el conjunto de datos ALBERS. Haga clic en él para descargar el archivo SA_NLC_2018_Albers_CALC_data.zip.

../../_images/data3.png
  1. A continuación, busque New South African National Land Cover (SANLC) 2014 and 2018 Change Assessment Datasets Computer Automated Land Cover (CALC) y haga clic en SA_NLC_2014_2018_CLASS_CHANGE_CALC (DATASET AND REPORT) para descargar el archivo SA_NLC_2014_2018_CLASS_CHANGE_ALBERS_CALC.tif.vat.zip.

../../_images/data4.png

La Ciudad de Johannesburgo publica conjuntos de datos espaciales como parte del Marco de Desarrollo Espacial 2040 (SDF) para Johannesburgo data. De este sitio descargaremos el archivo shape de los límites.

  1. Haga clic en el enlace de descarga http://bit.ly/joburg-sdf-16.

../../_images/data5.png
  1. Haga clic en el directorio SDF Shapefiles.

../../_images/data6.png
  1. Descargue el archivo SDF Shapefiles.zip y descomprímalo en una carpeta.

../../_images/data7.png

Para su comodidad, en el siguiente enlace encontrará una versión recortada de las capas necesarias:

SA_NLC_2018_ALBERS_CALC.tif

SA_NLC_2014_2018_CLASS_CHANGE_ALBERS_CALC.tif

COJ_Boundary.zip

Fuente de Datos: [SANLC] [COJ]

Procedimiento

  1. Busque la carpeta descargada y descomprimida en el navegador. Expándala y arrastra y suelta el SA_NLC_2018_ALBERS_CALC.tif en el lienzo.

../../_images/011.png
  1. Una vez cargada la capa, podrá observar que el SRC estará configurado como SRC Desconocido en la parte inferior derecha. Haga doble clic en él para abrir el cuadro de diálogo Propiedades del proyecto - SRC.

../../_images/021.png
  1. En la parte inferior, verá una vista previa de la extensión de la proyección. Esta SRC desconocida es una proyección de área equitativa de Lambert personalizada definida para el país de Sudáfrica. Más adelante reproyectaremos esta capa a otra proyección. Haga clic en Aceptar.

../../_images/031.png
  1. Cargue las otras dos capas SA_NLC_2014_2018_CLASS_CHANGE_ALBERS_CALC.tif y COJ_Boundary. Verá que la capa rasterizada cubre todo el país. Para nuestro análisis, sólo nos interesa el área cubierta por la capa COJ_Boundary. Ahora recortaremos la capa rasterizada a esta región. Vaya a Procesamiento ‣ Caja de herramientas ‣ GDAL ‣ Extracción de rásteres ‣ Recortar ráster por capa de máscara. Haga doble clic para abrirla.

../../_images/041.png
  1. En la Cortar ráster por capa de máscara, seleccionamos SA_NLC_2018_ALBERS_CALC como Capa de entrada, y luego COJ_Boundary como Capa de máscara. También tenemos una opción para reproyectar los datos a otra proyección. Es una buena práctica mantener todas las capas de datos en la misma proyección. Reproyectaremos los rásteres para que el SRC coincida con el de la capa COJ_Boundary. Seleccione EPSG:4326 - WGS 84 como SRC objetivo`.

../../_images/051.png
  1. The default output data format is GeoTiff. GeoTiff files can get very large if they are not compressed. A good practice is to always apply a loss-less compression when creating new raster layers. Expand Advanced Parameters and choose High Compression as the Profile. Next, click the ... button next to Clipped (mask) and select Save to file… to enter layer name as SA_NLC_2018_Clipped. Click Run.

../../_images/061.png
  1. Una vez que el algoritmo termine, no cierre la ventana. Vamos a aplicar la misma operación a la otra capa rasterizada. Cambie a la pestaña Parámetros y cambie la Capa de entrada a SA_NLC_2014_2018_CLASS_CHANGE_ALBERS_CALC. Mantenga todas las demás opciones pero cambie el nombre de la capa de salida a SA_NLC_2014_2018_CLASS_CHANGE_Clipped. Haga clic en Ejecutar.

../../_images/071.png
  1. Las dos capas recortadas se cargarán ahora en el lienzo. Seleccione la capa original y haga clic en Quitar capa para eliminarlas.

../../_images/081.png
  1. Las tres capas restantes están ahora en el mismo SRC. Ahora podemos cambiar el SRC del proyecto al SRC de las capas. Haga clic con el botón derecho en cualquier capa recortada y elija Capa SRC ‣ Establecer SRC del proyecto desde la capa.

../../_images/091.png
  1. Now the project CRS will be set to EPSG:4326. Bring the SA_NLC_2018_Clipped layer to top.

../../_images/101.png
  1. Haga clic en SA_NLC_2018_Clipped y utilice la herramienta identificar de la barra de herramientas Atributos para hacer clic en la imagen e inspeccionar los valores de los píxeles. Verá que los valores de los píxeles van de 1 a 73. Estos valores representan una clase distinta de uso/cobertura del suelo.

../../_images/111.png
  1. Las clases del conjunto de datos se describen en la Presentación de SANLC 2018, que puede descargarse del ``Portal de EGI <https://egis.environment.gov.za/sa_national_land_cover_datasets>`_. Para este ejercicio, nos interesan los asentamientos informales representados por los números de clase 51 a 54.

../../_images/121.png
  1. Vamos a extraer los píxeles que pertenecen a estas clases. Vaya a la herramienta Procesamiento ‣ Caja de herramientas ‣ Análisis rasterizado ‣ Calculadora rasterizada. Haga doble clic para abrirla.

../../_images/13.png
  1. Click the ... button next to Input layers. Select the SA_NLC_2018_Clipped layer and click OK. Next, click on the ε button to open the expression window.

../../_images/14.png
  1. Enter the following expression to select pixels from class 51-54. The source image has only 1 band. The @1 suffix indicates the band number. Click OK to close the expression window.

"SA_NLC_2018_Clipped@1" >= 51 AND "SA_NLC_2018_Clipped@1" <= 54
../../_images/15.png
  1. Next, click the ... button next to Calculated and select Save to File….

../../_images/16.png
  1. Nombre el archivo de salida residencial_informal.tif y haga clic en Ejecutar.

../../_images/17.png
  1. Una vez finalizado el proceso, se añadirá a QGIS una nueva capa, residencial_informal. Esta capa rasterizada sólo tiene valores de dos píxeles: 1 donde nuestra expresión evaluó true y 0 donde fue false. Los píxeles que aparecen en blanco son los que pertenecen a las clases de asentamientos informales. Estilizaremos mejor esta capa para poder ver claramente los asentamientos informales. Haga clic en el botón Abrir el panel de estilo de la capa.

../../_images/18.png
  1. Select the residential_informal layer and change the renderer to be Paletted/Unique values. Click the Add values manually (+) button.

../../_images/19.png
  1. Cambie el Valor a 1 y escriba Residencial Informal como Etiqueta. Seleccione un color de su elección.

../../_images/20.png
  1. Ahora podemos ver todos los asentamientos informales de la ciudad de Johannesburgo. Sería útil verlos en contexto con un mapa base. Tenemos acceso a una variedad de mapas base desde el complemento QuickMapServices. Una vez instalado el complemento, ve a Web ‣ QuickMapServices ‣ OSM ‣ OSM Standard para añadir la capa de OpenStreetMap.

../../_images/21.png
  1. Ahora puede identificar y verificar fácilmente si nuestro análisis ha identificado correctamente los asentamientos informales. Puede seleccionar la capa residencial_informal y cambiar a la pestaña Transparencia en el panel Estilos de capa`. Puede reducir la Opacidad global para ver juntos los píxeles extraídos y el mapa base.

../../_images/22.png
  1. Ya ha completado la primera parte del tutorial. Ahora utilizaremos la capa ráster SA_NLC_2014_2018_CHANGE_Clipped para identificar las regiones que se urbanizaron entre 2014 y 2018. Desactive todas las capas excepto SA_NLC_2014_2018_CHANGE_Clipped, y haga clic en el botón Abrir el panel de estilo de capas`. Cambie a la pestaña Transparencia e introduzca 0 en Valor adicional sin datos`. Esto hará que los píxeles con valor 0 sean transparentes.

../../_images/23.png
  1. Utilice la herramienta Identificar de la barra de herramientas Atributos para hacer clic en la imagen e inspeccionar los valores de los píxeles. Verá que los valores de los píxeles van de 21 a 420. Cada valor indica una transición de una de las 73 clases de origen a otra clase.

../../_images/24.png
  1. Su descarga de datos viene con una hoja de cálculo llamada lcccodes.xlsx. Esta hoja tiene una hoja 03 urban_change_codes que da más detalles sobre cada valor de píxel. Nos interesan todos los valores de los píxeles en los que cualquier clase de 2014 cambió a una clase edificada de 2018. En la imagen de abajo, estos están resaltados en azul.

../../_images/25.png
  1. Nuestro objetivo es mapear los cambios en la clase construida. Aplicaremos una transformación en la capa SA_NLC_2014_2018_CHANGE_Clipped para que todos los valores de los píxeles sean mapeados desde sus valores originales a cualquiera de los siguientes valores.

1

Todos los píxeles que fueron una clase construida tanto en 2014 como en 2018

2

Todos los píxeles que cambiaron de clase no construida en 2014 a clase construida en 2018.

0

Todos los píxeles que quedan

  1. Para ello, necesitamos crear una tabla que especifique estas reglas. Como QGIS es capaz de leer hojas de cálculo directamente, es el método más conveniente para crear esta tabla. Nuestra hoja de cálculo debe tener 3 columnas, MIN, MAX, y OUTPUT. Cada fila debe ser el rango de valores ráster de entrada al que se le debe asignar un valor de salida. Cree una hoja de cálculo como la que se muestra a continuación y guárdela en su ordenador como reclass.xlsx. También puede descargar una copia lista para usar desde este enlace - reclass.xlsx

../../_images/27.png
  1. Localice el archivo reclass.xlsx en el navegador. Arrástrelo y suéltelo en la ventana principal.

../../_images/28.png
  1. A new layer Sheet1 will be added to the Layers panel. Right-click on it and select Open Attribute Table. Verify that the sheet was imported correctly and you have 3 columns named MIN, MAX and OUTPUT. Open the Processing Toolbox ‣ Reclassify by layer tool.

../../_images/29.png
  1. In the Reclassify by layer dialog, select SA_NLC_2014_2018_CHANGE_Clipped as the Raster layer. Select Sheet1 as the Layer containing class breaks. Select MIN, MAX and OUTPUT fields for their respective fields.

../../_images/30.png
  1. Amplíe la sección Parámetros avanzados. Cambie la etiqueta Límites de rango a min <= value <= max. Haga clic en el botón de Recalificación de la trama e introduzca el nombre del archivo de salida como builtup_change.tif. Haga clic en Ejecutar.

../../_images/31.png
  1. Una vez finalizado el procesamiento, se añadirá al lienzo una nueva capa builtup_change con valores de píxel 0-2. En el panel de Estilos de capa, seleccione Valores únicos y haga clic en Clasificar.

../../_images/32.png
  1. Elija el color de su elección para cada categoría y etiquete los valores de los píxeles «0», «1» y «2» como «no edificado», «edificado existente» y «edificado nuevo».

../../_images/33.png
  1. Ahora, en la pestaña Transparencia, reduzca la Opacidad globa, y active la capa OSM Standard para ver tanto los píxeles de buildup_change como el mapa base juntos.

../../_images/34.png

If you want to give feedback or share your experience with this tutorial, please comment below. (requires GitHub account)