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:
Publicar un comentario