Remuestreo de imágenes GeoTIFF a la misma resolución utilizando QGIS

Question

Tengo dos imágenes geotiff que me gustaría remuestrear a la misma resolución.

Si utilizo ArcGIS Desktop esto es fácil de hacer, pero ¿cómo lo hago utilizando QGIS?

Answer 1

Esto también es fácil en QGIS, aunque un poco menos obvio. Hay un par de maneras de hacerlo:

• Calculadora de Raster - simplemente utilice la calculadora de raster y puede establecer la resolución y la extensión allí y puede hacer que coincidan con otro raster seleccionando la banda de raster que quiere hacer coincidir en la lista de Bandas de Raster y luego haciendo clic en el botón "Extensión de la capa actual". Los campos de columnas y filas le permitirán establecer la resolución. Sin embargo, este método no le permite controlar el método de remuestreo.

• Utilizando GDAL_Warp - esta herramienta le permite establecer la resolución de salida, ya sea especificando la anchura y la altura de la trama de salida o especificando el interruptor -tr (véase la documentación). Puede acceder a la herramienta GDAL_warp yendo a Raster->Projections->Warp (¡he dicho que no era obvio desde el punto de vista del remuestreo!).

  • (v2.x) Si quieres usar el interruptor -tr, rellena todas las casillas para la trama de entrada y la salida, etc. (tus valores SRS de origen y destino serán presumiblemente los mismos en este caso - aunque no tienen que serlo si también estás reproyectando). A continuación, haga clic en el pequeño icono del lápiz en la parte inferior y edite la línea de comandos gdal-warp generada automáticamente para incluir su interruptor -tr. Gdal_wrap le permite especificar el algoritmo que desea utilizar para el remuestreo y por lo tanto es un poco menos de un instrumento contundente que el uso de la calculadora de trama.

  • (v3.x) El parámetro -tr se habilita utilizando la opción Resolución del archivo de salida en unidades georreferenciadas de destino caja. Por ejemplo, para reducir la muestra de un MDE de 1 m a un MDE de 2 m, puede introducir 2 en ese campo. Sin embargo, no existe la opción de pasar dos argumentos diferentes para los píxeles no cuadrados. Digamos que su tamaño de píxel objetivo es 0.3125,0.25 , lo que significa que el xres es 0.3125 y el yres es 0.25 . Si ahora pasas el valor 0.3125 en esa casilla, se pondrá -tr 0.3125 0.3125 en el comando. Para contrarrestar esta limitación, simplemente copie el código, péguelo en la línea de comandos, edite la bandera -tr y ejecútelo. Por ejemplo:

    gdalwarp -t_srs EPSG:4326 -tr 0.3125 0.25 -r near -te 71.40625 24.875 84.21875 34.375 -te_srs EPSG:4326 -of GTiff foo.tiff bar.tiff

    (dependiendo de tu instalación y de las variables de entorno, puede que también tengas que indicar explícitamente la ruta de acceso a gdalwarp).

Answer 2

Normalmente utilizo la herramienta de Alineación de Tramas.

Está localizado en la pestaña de rasterización y se puede utilizar para establecer rasterizaciones en la misma medida también. También puede realizar algunas operaciones al remuestrear a una nueva resolución (como la media, el valor mínimo y máximo, etc.).

Aquí hay un enlace con algo de información al respecto: https://docs.qgis.org/2.18/en/docs/user_manual/working_with_raster/raster_analysis.html?highlight=raster%20alignment#id3

Yo uso la versión 2.18.13. No sé si la herramienta está disponible en versiones anteriores.

Answer 3

En QGIS esto se puede hacer utilizando r.rsamp.interp de la caja de herramientas de GRASS. Permite aumentar o disminuir la resolución del raster.

Esta herramienta permite seleccionar muchos Métodos de interpolación por muestreo (más cercano, bilineal, bicúbico). También permite recortar o enmascarar una imagen durante el proceso de remuestreo con el botón Ampliación de la región GRASS . Para seleccionar el área de recorte es posible utilizar una vista vectorial, de trama, de lienzo o simplemente seleccionando un área con el ratón.