¿Cómo cargar múltiples capas vectoriales en el lienzo de mapas de QGIS utilizando scripts de python?

Question

Estoy tratando de escribir un script independiente en python para crear un lienzo de mapa QGIS. Las capas vectoriales que estoy tratando de cargar son Shapefiles.

Actualmente, puedo cargar cada capa vectorial, añadirlas al registro del mapa y mostrarlas en mi aplicación de lienzo de mapa. El problema que tengo es que necesito que todas estas capas vectoriales se muestren en el lienzo del mapa al mismo tiempo. Actualmente el lienzo del mapa de qgis sólo muestra la última capa vectorial que se estableció en "canvas.setExtent(vector_layer.extent())". He intentado crear una lista de todas mis capas y luego ponerla en canvas.setExtent(list_of_layers.extent()) pero las listas no tienen el método "extent.()".

He leído el cookbook de pyQGIS y entiendo que el método .extent() es para un determinado vector que se quiere ver... ¿Cómo puedo mostrar todas las capas vectoriales a la vez (sin usar la consola de Python)? He indagado en otros Q&As que tenían preguntas similares a esta; sin embargo, no estoy encontrando ninguna solución que realmente aborde mi problema.

Aquí está mi siguiente código:

import modules
from qgis.core import *
from qgis.gui import *
from PyQt4.QtCore import *
from PyQt4 import QtGui
import os, sys
import qgis

def main():

    # supply path to qgis isntall location
    QgsApplication.setPrefixPath("/usr", True)

    # create a reference to the QgsApplication, setting the 
    # second arguement to False disables the GUI
    app = QgsApplication([], True)

    # load providers
    app.initQgis()

    # create Qt widget
    canvas = QgsMapCanvas()
    canvas.setCanvasColor(Qt.white)

    # enable this for smooth rendering
    canvas.enableAntiAliasing(True)

    # not updated US6SP10M files from ENC_ROOT
    source_dir = "/home/cassandra/desktop/file_formats/Shapefiles"

    # load vector layers
    for files in os.listdir(source_dir):

        # load only the shapefiles
        if files.endswith(".shp"):
            vlayer = QgsVectorLayer(source_dir + "/" + files, files, "ogr")

            # add layer to the registry
            QgsMapLayerRegistry.instance().addMapLayer(vlayer)

            # set extent to the extent of our layer
            canvas.setExtent(vlayer.extent())

            # set the map canvas layer set
            canvas.setLayerSet([QgsMapCanvasLayer(vlayer)])

    # refresh canvas and show it        
    canvas.refresh()
    canvas.show()
    app.exec_()
    app.exitQgis()

main()

Answer 1

¡Lo he descubierto! Era algo muy pequeño. Básicamente, las capas que se quieren ver en el lienzo se proporcionan a través de canvas.setLayerSet(list_of_layers).

# total list of layers actually displayed on map canvas
canvas_layers = []

# load vector layers
for files in os.listdir(source_dir):

    # load only the shapefiles
    if files.endswith(".shp"):

        # create vector layer object
        vlayer = QgsVectorLayer(source_dir + "/" + files, files, "ogr")
        print(files)

        # add the layer to the registry
        QgsMapLayerRegistry.instance().addMapLayer(vlayer)

        # combine extent of the current vector layer with the extent of the created "extent" rectangle object
        extent.combineExtentWith(vlayer.extent())
        canvas_layers.append(QgsMapCanvasLayer(vlayer))

# set extent to the extent of a larger rectangle so we can see all geometries       
canvas.setExtent(extent)

# provide set of layers for display on the map canvas
canvas.setLayerSet(canvas_layers)