Generujesz dużą liczbę map za pomocą PyQGIS?

Muszę wykonać dużą liczbę (setki) map rozmieszczenia gatunków. Mam plik kształtu, który zawiera rozkłady dla każdego gatunku, a dla każdego z nich chciałbym uzyskać mapę jako obraz (jpg, png lub inny), który zawiera nazwę danego gatunku, legendę (w celu rozróżnienia obszarów roczne dystrybucje, hodowla, hodowla itp.).

Chciałbym do tego użyć QGIS.

Miałem podobne wymaganie i stworzyłem wtyczkę QGIS do generowania map w oparciu o plik kształtu z lokacjami punktowymi dla wszystkich gatunków (przyjmuje on unikalną nazwę taksonu w tabeli atrybutów jako wspólny identyfikator). Moje wymagania nie były tak skomplikowane - nie potrzebowałem sezonowych informacji, tytułów ani legend, ale może to być przydatny punkt wyjścia. W przypadku bardziej złożonych aspektów konieczne będzie skorzystanie z narzędzia do tworzenia map. Więcej informacji na ten temat zawiera książka kucharska PyQGIS .

Podłącz

Wtyczka automatyzuje tworzenie map i pozwala konfigurować zakresy, rozdzielczość i inne aspekty. Stosuje ten sam styl do wydruku, co nakładka siatki. Obecnie działa tylko w wersji rozwojowej QGIS (1.9 lub nowszej).

• Krótka strona wstępna
• Strona wtyczki QGIS
• Repozytorium GitHub

Skrypt Sextante

Przed stworzeniem wtyczki opracowałem logikę za pomocą SEXTANTE. Ten skrypt użytkownika powinien również działać w wersji 1.8 (nie testowałem go). Plik stylu dystrybucji (.qml) to styl dystrybucji danych wyjściowych (ignoruje styl nakładki dystrybucji). Obecnie umieszcza mapy wyników w katalogu temp w oparciu o wartości domyślne systemu operacyjnego (/ tmp w systemie Linux i różne miejsca w systemie Windows - zdefiniowane przez zmienną środowiskową TEMP). Możesz jednak dość łatwo to zdefiniować w kodzie. Będziesz także musiał edytować zakres i rozdzielczość wyjściową w kodzie (i kolor tła, jeśli chcesz inny kolor dla morza).

#Definition of inputs and outputs
#==================================
##[Scratch]=group
##all_localities=vector
##taxon_field=field all_localities
##africa_map=vector
##sa_map=vector
##grid_layer=vector
##distribution_style_file=file

#Algorithm body
#==================================
from qgis.core import *
from PyQt4.QtCore import *
from PyQt4.QtGui import *
from sextante.core.QGisLayers import QGisLayers
from sextante.core.SextanteVectorWriter import SextanteVectorWriter
import tempfile
import os

def print_map(taxon,taxon_shp):
    #load taxon layer (necessary?)
    #QGisLayers.load(taxon_shp,name = "taxon",style = distribution_style_file)
    taxon_layer = QgsVectorLayer(taxon_shp,"taxon","ogr")
    QgsMapLayerRegistry.instance().addMapLayer(taxon_layer)
    taxon_layer.loadNamedStyle(distribution_style_file)

    # create image (dimensions 325x299)
    img = QImage(QSize(325,299), QImage.Format_ARGB32_Premultiplied)

    # set image's background color
    color = QColor(192,192,255)   # blue sea
    img.fill(color.rgb())

    # create painter
    p = QPainter()
    p.begin(img)
    p.setRenderHint(QPainter.Antialiasing)

    render = QgsMapRenderer()

    # create layer set
    africa_layer = QGisLayers.getObjectFromUri(africa_map)
    sa_layer = QGisLayers.getObjectFromUri(sa_map)
    #taxon_layer = QGisLayers.getObjectFromUri(taxon_shp)

    lst = []
    lst.append(taxon_layer.id())    
    lst.append(sa_layer.id())
    lst.append(africa_layer.id())

    render.setLayerSet(lst)

    # set extent (xmin,ymin,xmax,ymax)
    rect = QgsRectangle(14.75,-36.00,34.00,-21.00)
    render.setExtent(rect)

    # set output size
    render.setOutputSize(img.size(), img.logicalDpiX())

    # do the rendering
    render.render(p)
    p.end()

    # save image
    #outdir = os.path.dirname(os.path.abspath(output))
    tempdir = tempfile.gettempdir()
    img.save(os.path.join(tempdir,taxon+".png"),"png")

    # remove taxon layer from project
    QgsMapLayerRegistry.instance().removeMapLayers([taxon_layer.id()])

tempdir = tempfile.gettempdir()   
taxa = sextante.runalg('qgis:listuniquevalues', all_localities, taxon_field, None)['UNIQUE_VALUES'].split(";")
for taxon in taxa:
    sextante.runalg('qgis:selectbyattribute', all_localities, taxon_field, 0, taxon)
    sextante.runalg('qgis:selectbylocation', grid_layer, all_localities, 0)
    filename = os.path.join(tempdir,"taxon.shp")    #memory file better?
    sextante.runalg('qgis:saveselectedfeatures', grid_layer, filename)
    print_map(taxon,filename)

Pracowałem dziś nad tym trochę czasu. Wprowadziłem więc kilka zmian w twoim skrypcie. Nie muszę dodawać dodatkowego kroku selectbyattribute i saveselectedfeatures, ponieważ używam plików .qml, a pole Seasonal znajduje się w tym samym pliku kształtu. Poniżej możesz zobaczyć, co zrobiłem:

#Definition of inputs and outputs
#==================================
##[Scratch]=group
##all_localities=vector
##taxon_field=field all_localities
##seasonal_field=field all_localities
##countries_map=vector
##distribution_style_file=file
##output_folder=folder

#Algorithm body
#==================================
from qgis.core import *
from PyQt4.QtCore import *
from PyQt4.QtGui import *
from sextante.core.QGisLayers import QGisLayers
from sextante.core.SextanteVectorWriter import SextanteVectorWriter
import tempfile
import os

def print_map(taxon,taxon_shp):
#load taxon layer (necessary?)
#QGisLayers.load(taxon_shp,name = "taxon",style = distribution_style_file)
taxon_layer = QgsVectorLayer(taxon_shp,"taxon","ogr")
QgsMapLayerRegistry.instance().addMapLayer(taxon_layer)
taxon_layer.loadNamedStyle(distribution_style_file)

# create image (dimensions 325x299)
img = QImage(QSize(325,299), QImage.Format_ARGB32_Premultiplied)

# set image's background color
color = QColor(221,249,254)   # blue sea
img.fill(color.rgb())

# create painter
p = QPainter()
p.begin(img)
p.setRenderHint(QPainter.Antialiasing)

render = QgsMapRenderer()

# create layer set
countries_layer = QGisLayers.getObjectFromUri(countries_map)
taxon_layer = QGisLayers.getObjectFromUri(taxon_shp)

lst = []
lst.append(taxon_layer.id())    
lst.append(countries_layer.id())
render.setLayerSet(lst)

# set extent (xmin,ymin,xmax,ymax)
rect = QgsRectangle(-11,32,39,71)
render.setExtent(rect)
# set output size
render.setOutputSize(img.size(), img.logicalDpiX())

# do the rendering
render.render(p)
p.end()

#save image
#outdir = os.path.dirname(os.path.abspath(output))
tempdir = output_folder
img.save(os.path.join(tempdir,taxon+".png"),"png")

# remove taxon layer from project
QgsMapLayerRegistry.instance().removeMapLayers([taxon_layer.id()])

tempdir = tempfile.gettempdir()  

taxa = sextante.runalg('qgis:listuniquevalues', all_localities, taxon_field, None)        ['UNIQUE_VALUES'].split(";")

for taxon in taxa:
sextante.runalg('qgis:selectbyattribute', all_localities, taxon_field, 0, taxon)
filename = os.path.join(tempdir,"taxon.shp")    #memory file better?
sextante.runalg('qgis:saveselectedfeatures', all_localities, filename)
print_map(taxon,filename)

Jeśli masz jakieś uwagi lub porady, aby je poprawić, nie wahaj się.

Aby to poprawić, najlepiej byłoby, gdybyśmy wybrali zakres (na przykład w Europie), wykorzystują ten zakres, aby wybrać tylko gatunki objęte tym zakresem. To dlatego, że otrzymuję mapy dla wszystkich gatunków, nawet tych, które są na przykład poza Europą (więc mam wiele pustych map). Czy uważasz, że to możliwe?

Twoje zdrowie,

Onesime