ベクターデータの読み込みと図形へのアクセス

ベクターデータを読み込みます。すでにQGISで読み込んでいる場合と、ファイルを直接指定して読み込む場合があります。また、読み込んだレイヤの図形への順次アクセスを学びます。

Menu

• 読み込み済みレイヤへのアクセス
• ベクターデータの新規読み込み
• アクティブなレイヤへのアクセス
• レイヤの表示/非表示の切り替え
• 図形への順次アクセス

読み込み済みレイヤへのアクセス

QGISで読み込み済みのレイヤへアクセスします。

読み込み済みのレイヤ一覧にはQgsProjectクラスのインスタンスからアクセス出来ます。シングルトンクラスになっているので、QgsProject.instance()を使用します。 まずはデータの読み込みがされていない状態実行してみましょう。何も読み込まれていないので空の辞書が返されます。

>>>QgsProject.instance().mapLayers()
{}

QGISのユーザーインターフェイス上で、データソースマネージャから、実習用にダウンロードしたne_50m_admin_0_countries.shpを追加します（データの取得方法は READMEを参照してください。）。先ほどと同様にデータにアクセスすると、読み込んだ画像のファイル名が辞書で返されます。複数のファイルを読み込んでいる場合には全てのファイル名が辞書に格納されます。

>>>QgsProject.instance().mapLayers()
{'ne_50m_admin_0_countries_8d55f8f5_ddff_4585_888a_caa06235d30a': <qgis._core.QgsVectorLayer object at 0x0000023676C76288>}

ベクターデータの新規読み込み

addVecterLayer()を使ってコンソールからレイヤを追加してみましょう。C:\work\ne_50m_populated_places.shpを読み込みます。最初の引き数は読み込むファイル名、2番目はレイヤ名称として使用されます。3番目の引き数は使用するQgisDataProviderへのキーになります。データの保存先ディレクトリは適宜、ご自分の環境に読み替えてください。 QgisDataProviderはデータソースからの読み込み、書き込み、ハンドリング全般を行うクラスです。扱うデータの形式に応じてQgisDataProviderを指定します。ここではベクターデータを扱っているのでogrを使用します。

>>>iface.addVectorLayer(r'C:\work\ne_50m_populated_places.shp', 'cities', 'ogr')
<qgis._core.QgsVectorLayer object at 0x0000023676C76318>

都市のポイントデータがマップキャンバスに追加されました。

アクティブなレイヤへのアクセス

現在選択されている(操作対象となっている)レイヤはifaceを通して取得します。

>>>layer = iface.activeLayer()

レイヤの表示/非表示の切り替え

QgsProject.instance()から、レイヤツリーにアクセスし、レイヤツリーを通して各レイヤの表示/非表示を切り替えます。

>>>QgsProject.instance().layerTreeRoot().findLayer(layer).setItemVisibilityChecked(False)
>>>QgsProject.instance().layerTreeRoot().findLayer(layer).setItemVisibilityChecked(True)

図形への順次アクセス

QGISでne_50m_admin_0_countries.shpのレイヤを選択し、各要素にアクセスしてみます。

name()で読み込んだレイヤの名称を表示します。

>>>layer = iface.activeLayer()
>>>layer.name()
'ne_50m_admin_0_countries'

QgsFeatureIteratorクラスを通して各図形に順次アクセスします。

>>>features = layer.getFeatures()
>>>for feature in features:
... feature.geometry()
<QgsGeometry: MultiPolygon (((31.28789062500001705 -22.40205078125001137, 31.197265625 -22.34492187500001137,
(以下略)

field()を使ってアクティブなベクターレイヤの全フィールド名を取得します。

>>>fields = layer.fields()
>>>for field in fields:
... print(field.name())
featurecla
scalerank
LABELRANK
(以下略)

for文の中で属性'NAME'を指定し、全地物の'NAME'を取得します。

>>>features = layer.getFeatures()
>>>for feature in features:
... feature['NAME']
'Zimbabwe'
'Zambia'
'Yemen'
(以下略)

属性'NAME'が'Japan'のベクトル要素を選択して、その要素の範囲にズームしてみましょう。

>>>features = layer.getFeatures()
>>>for feature in features:
... if feature['NAME'] == 'Japan':
...     break
>>>canvas = iface.mapCanvas()
>>>canvas.setExtent(feature.geometry().boundingBox())
>>>canvas.refresh()

教材の利用に関するアンケート