Tutorial de Rasterización de Shapefiles con Python y QGIS 3.0

Rasterize.png

QGIS 3.0 empezó a distribuirse este año y ha salido con variedad de novedades, incluye el cambio completo de su API a Python 3, el cual tiene diferencias ligeras de sintaxis con su predecesor, así mismo el código de la API de QGIS 3.0 ha variado drásticamente, incluyendo nuevos algoritmos y variando otros complemente, sin mencionar que su interfaz es completamente diferente.

En este tutorial se observará la rasterización de un shapefile con el complemento “GDAL” la cual es una librería muy potente que incluye variedad de algoritmos en los cuales se basa la mayoría de código de los SIG.

 

Video

 

Código de Python

Este es el código utilizado en este tutorial:

from qgis.utils import iface

import processing

from qgis.core import QgsProject


shape = r"C:\Users\Jonathan\FlopyShapes\flopy\QuilcaChili.shp"

layer = QgsVectorLayer(shape, 'cuenca')

#QgsProject.instance().addMapLayer(layer)

extents = layer.extent()

xmin = extents.xMinimum()

xmax = extents.xMaximum()

ymin = extents.yMinimum()

ymax = extents.yMaximum()

alg = 'gdal:rasterize'

params={"INPUT" : layer,"FIELD":None,"BURN":1,"UNITS":1,"WIDTH":100,"HEIGHT":100,"EXTENT" : "%f,%f,%f,%f"% (xmin, xmax, ymin, ymax),"NODATA":0,"DATA_TYPE":1,"INVERT":0,"OUTPUT" : r'C:\Users\Jonathan\Downloads\Raster.tif'}

test = processing.run(alg, params)

iface.addRasterLayer(r'C:\Users\Jonathan\Downloads\Raster.tif',"Raster")

Smiley face

Suscríbase a nuestro boletín gratuito para recibir noticias, datos interesantes y fechas de nuestros cursos en recursos hídricos.

Posted on August 30, 2018 and filed under TutorialPython, TutorialQGIS.