PyQGISを使用して多数のマップを生成しますか？

種の分布図を多数（数百）作成する必要があります。各種の分布を含むシェープファイルがあり、各種について、関係する種の名前、凡例（領域を区別するため）を含む画像（jpg、pngなど）としてマップを取得したい年間分布、繁殖、非繁殖など...）。

これを行うためにQGISを使用したいと思います。

同様の要件があり、すべての種の点の局所性を含むシェープファイルに基づいてマップを生成するQGISプラグインをまとめました（共通の識別子として属性テーブルの一意の分類名を想定しています）。私の要件はそれほど複雑ではありませんでした-私は季節的な情報、タイトル、または伝説を必要としませんでしたが、それはあなたにとって有用な出発点かもしれません。より複雑な側面については、マップコンポーザーを使用する必要があります。詳細については、PyQGISクックブックを参照してください。

プラグイン

プラグインは、マップの作成を自動化し、範囲、解像度、その他の側面を構成できます。グリッドオーバーレイと同じスタイルを出力に適用します。現在、QGISの開発バージョン（1.9以降）でのみ動作します。

• 短い紹介ページ
• QGISプラグインページ
• GitHubリポジトリ

Sextanteスクリプト

プラグインを作成する前に、SEXTANTEを使用してロジックを作成しました。このユーザースクリプトは1.8でも動作するはずです（テストしていません）。配布スタイルファイル（.qml）は、出力配布のスタイルです（配布オーバーレイのスタイルは無視されます）。現在、オペレーティングシステムのデフォルト（Linuxでは/ tmp、Windowsではさまざまな場所-TEMP環境変数で定義されています）に基づいて、出力マップを一時ディレクトリに配置します。ただし、コードで簡単に定義できます。また、コードの範囲と出力解像度（および海に別の色が必要な場合は背景色）を編集する必要があります。

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

今日は少し時間をかけてこれに取り組みました。だから私はあなたのスクリプトにいくつかの変更を加えました。.qmlファイルを使用し、Seasonalフィールドが同じシェープファイルにあるため、selectbyattributeおよびsaveselectedfeaturesステップを追加する必要はありません。以下に、私が行ったことを示します。

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

あなたがそれを改善するための発言やアドバイスがあれば、躊躇しないでください。

これを改善するには、範囲を選択する（たとえば、ヨーロッパの場合）ことが最善です。この範囲を使用して、この範囲内に含まれる種のみを選択します。これは、すべての種のマップを取得しているためです。たとえば、ヨーロッパ以外の種も含まれています（そのため、空のマップがたくさんあります）。それは可能だと思いますか？

乾杯、

Onesime