Come caricare l'output di memoria dall'elaborazione QGIS?

Molti algoritmi di elaborazione hanno un'opzione per salvare l'output come memoria. Se eseguito dalla casella degli strumenti, funziona correttamente, perché è possibile selezionare la casella "Apri file di output dopo l'esecuzione dell'algoritmo".

Guardando il codice sorgente per QGIS, il codice per caricare il livello di memoria sembra essere definito nella funzione Postprocessing.handleAlgorithmResults. La funzione utilizza l'accesso a un alg.outputselenco e procede come segue:

for out in alg.outputs:
    progress.setPercentage(100 * i / float(len(alg.outputs)))
    if out.hidden or not out.open:
        continue
    if isinstance(out, (OutputRaster, OutputVector, OutputTable)):
        try:
            if out.value.startswith('memory:'):
                layer = out.memoryLayer                                # NOTE!!
                QgsMapLayerRegistry.instance().addMapLayers([layer])   # NOTE!!
            else:
                # ... 

Quando esegui algoritmi di elaborazione dalla console, c'è un modo per caricare il layer senza accedere a questo oggetto? posso correre

processing.runalg("qgis:intersection", layer1, layer2, "memory:")

o anche

processing.runalg("qgis:intersection", layer1, layer2, "memory:myLayerName")

Non riesco comunque a trovare un modo per afferrare l'output risultante.

Aaaand l'ho trovato. Usa processing.runandload, che carica il livello di output nel sommario dopo aver eseguito l'algoritmo.

processing.runandload("qgis:intersection", layer1, layer2, "memory:myLayerName")
layer = QgsMapLayerRegistry.instance().mapLayersByName("memory:myLayerName")[0]
# Should do error checking as well, but this works!!

è il modo corretto, è spiegato nella documentazione http://docs.qgis.org/2.14/es/docs/user_manual/processing/console.html

il codice successivo funziona con in memoria tutto tranne l'ultimo che viene caricato

MDT=path/mdt.tif
drain=processing.runalg("grass:r.drain",MDT,"",(pun),False,False,False,"%f,%f,%f,%f"% (xmin, xmax, ymin, ymax),0,-1,0.00100,None)
vect=processing.runalg("grass:r.to.vect",drain['output'],0,False,"%f,%f,%f,%f"% (xmin, xmax, ymin, ymax),0,None)
bu=processing.runalg("qgis:fixeddistancebuffer",vect['output'],Metros_afecta,1,False,None)
buf=bu['OUTPUT']
bufe= QgsVectorLayer(buf,"area", "ogr")
#the last load the layer 
QgsMapLayerRegistry.instance().addMapLayers([bufe])

la processing.runalgrestituisce un dizionario in questo caso bu['OUTPUT']. OUTPUTÈ LA CHIAVE e il valore è un percorso temporaneo. Puoi vedere la chiave con processeing.alghelp ("elaborazione del nome") come elaborazione, alghelp ("grass: r.drain")

ritorno

processing.alghelp("grass:r.drain")
ALGORITHM: r.drain - Traces a flow through an elevation model on a raster map.
input <ParameterRaster>
coordinate <ParameterString>
vector_points <ParameterMultipleInput>
-c <ParameterBoolean>
-a <ParameterBoolean>
-n <ParameterBoolean>
GRASS_REGION_PARAMETER <ParameterExtent>
GRASS_REGION_CELLSIZE_PARAMETER <ParameterNumber>
GRASS_SNAP_TOLERANCE_PARAMETER <ParameterNumber>
GRASS_MIN_AREA_PARAMETER <ParameterNumber>
output <OutputRaster>

in questo caso la chiave viene emessa , fai attenzione con le lettere maiuscole che devi scrivere in maiuscolo o senza maiuscolo, in questo caso non maiuscolo.

Non sono sicuro se questo ti aiuterà nel tuo contesto (vuoi eseguire il tuo codice come autonomo o all'interno di QGIS dopo che un algoritmo di elaborazione è terminato?). In quest'ultimo caso, è possibile eseguire una query su oggetti caricati QGsVector-e facilmente QGsRasterlayeresaminando il registro di QGIS MapLayer.

# General function to retrieve a layer with a given name
def getLayerByName( layerName ):
  layerMap = QgsMapLayerRegistry.instance().mapLayers()
  for name, layer in layerMap.iteritems():
    if layer.name() == layerName:
        if layer.isValid():
          return layer
        else:
          return None

In generale, tutti i layer devono avere una sorgente da qualche parte sul disco rigido anche se il risultato viene appena aggiunto dopo che l'elaborazione ha terminato i suoi calcoli. Se guardi nei metadati dei livelli puoi trovare l'origine effettiva di un livello (se ha origine dall'elaborare normalmente da qualche parte in una cartella temporanea)