Muestrear Datos Ráster usando Puntos o Polígonos (QGIS3)

Varios conjuntos de datos científicos y ambientales vienen como rásters grid. Los datos de elevación (DEM) son también distribuídos como archivos rásters. En estos archivos rásters, el parámetro que está siendo representado está codificado como los valores píxeles del ráster. A menudo, uno necesita extraer los valores píxeles en ciertas ubicaciones o agregarlos sobre un área. Esta funcionalidad está disponible en QGIS mediante algoritmos de procesamiento. Muestrear valores de ráster para capas punto y Estadísticas Zonales para capas polígono.

Vista general de la tarea

Dada una cuadrícula ráster de temperatura máxima en los EE.UU. continental, necesitamos extraer la temperatura en una capa punto de todas las áreas urbanas y calcular la temperatura promedio para una capa polígono de cada condado en los EE.UU.

Otras habilidades que aprenderá

  • Seleccionar y quitar múltiples capas de la Tabla de Contenido QGIS.

Obtener los datos

El Centro de Predicción del Clima de NOAA provee datos SIG relacionados con temperatura y precipitación en los EE.UU. Descargue el último archivo grid para temperaturas máximas. El archivo tendrá como nombre us.tmax_nohads_ll_{YYYYMMDD}_float.tif

Usaremos un archivo CSV del 2018 US Gazetteer que representa las áreas urbanas en los EE.UU. Descargue el Archivo Urban Areas Gazetteer.

US Census Bureau proporciona archivos shape TIGER/Line. Puede visitar el sitio FTP y descargar Archivo shape Census Tracts.

Para su comodidad, puede descargar directamente una copia de los conjuntos de datos de los enlaces abajo:

us.tmax_nohads_ll_20190501_float.tif

2018_Gaz_ua_national.zip

tl_2018_us_county.zip

Fuentes de Datos: [NOAACPC], [USGAZETTEER] [TIGER]

Procedimiento

  1. Descomprima y extraiga tanto 2018_Gaz_ua_national.zip como tl_2018_us_county.zip a una carpeta en su computadora. Abra QGIS y localice el archivo us.tmax_nohads_ll_20190501_float.tif. En el Explorador QGIS, arrástrelo a la pantalla.

../../_images/1154.png
  1. Verá una nueva capa ráster us.tmax_nohads_ll_20190501_float cargada en el panel Capas. Esta capa ráster contiene la temperatura máxima registrada en cada píxel en grados Celsius. A continuación cargaremos el archivo punto de áreas urbanas. Este archivo viene como un archivo texto en el formato Valores Separados por Tabulación (TSV). Clic el botón Abrir Administrador de Fuente de Datos en la Barra de Herramientas Fuente de Datos.

../../_images/2129.png
  1. Cambie a la pestaña Texto Delimitado. Clic el botón junto a Nombre de archivo y especifique la ruta al archivo texto que descargó. En la sección Formato de archivo, seleccione Delimitadores personalizados y marque Tabulación. Seleccione INTPTLONG como el campo X y INTPTLAT como el campo Y. Clic en Añadir y luego Cerrar.

../../_images/367.png
  1. Una nueva capa punto 2018_Gaz_ua_national se cargará en el panel Capas. Ahora estamos listos para extraer los valores de la capa ráster en esos puntos. Vaya a Procesos ‣ Caja de herramientas.

../../_images/438.png
  1. Busque y localice el algoritmo menuselection:Análisis ráster –> Muestrear valores ráster. Haga doble-clic para iniciarlo.

../../_images/534.png
  1. Seleccione 2018_Gaz_ua_national como la Capa de Punto de Entrada. Seleccione us.tmax_nohads_ll_20190501_float  como la Capa Ráster a Muestrear. Expanda los Parámetros avanzados e ingrese tmax como el Prefijo de columna de salida. Clic en Ejecutar. Una vez que termine el procesamiento, clic en Cerrar.

../../_images/633.png
  1. Una nueva capa Sampled Points será cargada en el panel Capas. Seleccione la herramienta Identificar en la Barra de herramientas Atributos y haga clic en cualquier punto. Verá los atributos mostrados en el panel Resultados de identificación. Verá un nuevo atributo llamado tmax_1 agregado a cada objeto. Éste es el valor píxel de la capa ráster extraído en la ubicación del punto. El 1 representa el número de banda del ráster. Si la capa ráster tenía múltiples bandas, vería múltiples columnas nuevas en la capa de salida.

../../_images/733.png
  1. La primera parte de nuestro análisis está completo. Quitemos las capas innecesarias. Mantenga presionada la tecla Mayús y seleccione las capas Sampled Points` y ``2018_Gaz_ua_national. Clic-derecho y seleccion Borrar para eliminarlas de QGIS. Cuando se le pregunta Eliminar 2 entradas de leyenda?, seleccione Aceptar.

../../_images/832.png
  1. Ahora usaremos la capa de condados para muestrear el ráster y calcular la temperatura promedio para cada condado. Localice el archivo tl_2018_us_county.shp en el Explorador QGIS y arrástrelo a la pantalla.

../../_images/933.png
  1. Se cargará una nueva capa tl_2018_us_county al panel Capas. Vaya a Procesos ‣ Caja de herramientas.

../../_images/1039.png
  1. Busque y localice el algoritmo Análisis ráster ‣ Estadísticas zonales y haga doble-clic para iniciarlo.

../../_images/1155.png
  1. Seleccione us.tmax_nohads_ll_20190501_float como la Capa ráster y tl_2018_us_county como la Capa vectorial que contiene zonas. Ingrese tmax_ como el Prefijo de columna salida. Clic en junto a Estadísticas a calcular.

../../_images/1242.png
  1. Seleccione sólo el valor Mean y clic en Aceptar.

../../_images/1340.png
  1. Ahora, haga clic en junto a :guilabel: Estadísticas zonales y seleccione la opción Guardar en archivo para guardar el resultado.

    ../../_images/1436.png
  2. Vaya a la carpeta donde desea guardar la salida. Guarde la capa de salida como us_county_tmax.shp. Haga clic en :guilabel: Ejecutar para iniciar el proceso. El algoritmo puede tardar unos minutos en completarse y verá la capa de salida añadida tan pronto como finalice el procesamiento. Haga clic en Cerrar.

../../_images/1532.png
  1. Haga clic con el botón derecho en la capa us_county_tmax y seleccione :guilabel: Abrir tabla de atributos.

../../_images/1629.png
  1. Usted verá una nueva columna llamada tmax_mean añadido en la tabla de atributos. Contiene el valor medio de temperatura extraído sobre el polígono para cada característica. Hay algunos valores nulos porque esos condados (pertenecientes a Alaska, Hawai y Puerto Rico) están fuera de la extensión de la capa ráster.

../../_images/1731.png

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