viernes, 10 de noviembre de 2017

Reproyectando un ráster con processing gdalogr:warpreproject

Una forma expedita de reproyectar un ráster con PyQGIS es mediante la herramienta gdalogr:warpreproject de Processing.

La información sobre los parámetros a ser empleados se encuentran en la imagen siguiente:


El código de ejemplo se dispone a continuación:

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
import processing

uri = "/home/zeito/pyqgis_data/utah_demUTM12.tif"
layer = processing.getObjectFromUri(uri)
SOURCE_SRS = layer.crs().authid()
DEST_SRS = 4326

processing.runandload("gdalogr:warpreproject",
                     {"INPUT":layer,
                      "SOURCE_SRS":SOURCE_SRS,
                      "DEST_SRS":"EPSG:"+ str(DEST_SRS),
                      "NO_DATA":"0",
                      "METHOD":0,
                      "RTYPE":5,
                      "COMPRESS":0,
                      "BIGTIFF":2,
                      "OUTPUT":None})

que se empleó en la Python Console para reproyectar un ráster desde EPSG:32612 (proyección leída automáticamente en el código) a EPSG:4326 y cuyo resultado se visualiza en la imagen siguiente:


No hay comentarios: