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