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

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")