[QGIS-it-user] sviluppo di uno script processing, importa layer in numpy

[Giacomo Fontanelli]

Salve

sto scrivendo un mio script processing

nella prima parte del processing vero e proprio (scusate il gioco di parole"
faccio una rasterizzazione di un layer vettoriale poligonale

        processPar = {
            "INPUT":parameters[self.INPUT],
            "FIELD":"id",
            "BURN":1,
            "UNITS":0,
            "WIDTH":nCol,
            "HEIGHT":nRow,
            "EXTENT": parameters[self.INPUT],
            "NODATA":"NA",
            "OPTIONS":"",
            "DATA_TYPE":0,
            "INIT":0,
            "INVERT":False,
            "EXTRA":"",
            "OUTPUT":parameters[self.OUTPUT]}

        rasterizedResult = processing.run(
            "gdal:rasterize",
            processPar,
            is_child_algorithm=True,
            context=context,
            feedback=feedback)

nella seconda parte vorrei poter importare il file raster prodotto prima
(in teoria)
da gdal:rasterize in numpy.

Per fare questo credo di dover passare prima da gdal, quindi faccio

        rasterGDAL = gdal.Open(rasterizedResult[OUTPUT])

e poi

  band = rasterGDAL.GetRasterBand(1)
        tempArray = band.ReadAsArray()

Purtroppo però alla linea

        rasterGDAL = gdal.Open(rasterizedResult[OUTPUT])

ottengo un errore

File "<string>", line 218, in processAlgorithm
NameError: name 'OUTPUT' is not defined


Un piccolo aiuto?

Grazie
-------------- parte successiva --------------
Un allegato HTML è stato rimosso...
URL: <http://lists.osgeo.org/pipermail/qgis-it-user/attachments/20201117/2c2d1f3a/attachment.html>