Vektorová data v QGIS Python API

Programování v GIS 4

Jan Caha

2025-09-30

Objekty pro práci s vektorovými daty

  • QgsVectorLayer - reprezentace vrstvy v QGIS
  • QgsVectorDataProvider- abstrakce čtení a zápisu (pokud je povolen) dat do zdroje dat
  • QgsFeature - reprezentace prvku ve vrstvě, přístup jeho geometrii a atributům
  • QgsGeometry - reprezentace geometrie prvku
    • QgsAbstractGeometry - abstraktní třída implementace geometrie - součást QgsGeometry
  • QgsFeatureRequest - objekt pro získání prvku z datového zdroje
  • QgsCoordinateReferenceSystem - reprezentace souřadnicového systému
  • QgsFields - seznam atributů vrstvy
  • QgsField- atribut a jeho definice

QgsVectorLayer

  • vytváření jako QgsVectorLayer(uri, name, driver)
    • uri- cesta k souboru, nebo např. připojení do DB
    • name - název vrstvy v QGIS
    • driver - obvykle ogr nebo postgres, memory ale i další
  • přístup k prvkům pomocí funkce getFeatures() vrací QgsFeatureIterator
    • lze vyspecifikovat buď pomocí QgsFeatureRequest nebo QgsRectangle

Tvorba QgsField

from qgis.PyQt.QtCore import QMetaType
from qgis.core import QgsField

field = QgsField("field", QMetaType.Type.Double)

QgsFeature

  • přístup k atributům a geometrii
  • do vrstvy se vkládá skrze QgsVectorDataProvider metodami addFeature() nebo addFeatures()
  • editace vrstvy pomocí funkcí changeAttributeValue(), changeGeometry(), changeAttributeValues()
  • vrstvu je nutné nejdřív otevřít pro editaci pomocí startEditing() následně uzavřít pomocí commitChanges(), nebo změny zrušit pomocí rollBack()

QgsGeometry

  • hlavní třída pro práci s geometrií
  • obsahuje řadu funkcí (i statické)
  • existují ale dilčí třídy odvozené od QgsAbstractGeometry pro konkrétní typy geometrie
    • QgsPoint, QgsLineString, QgsPolygon
    • QgsMultiPoint, QgsMultiLineString, QgsMultiPolygon
    • QgsGeometryCollection
    • QgsSurface, QgsCurve, QgsCompoundCurve, QgsTriangle

Tvorba QgGeometry

from qgis.core import QgsGeometry, QgsPointXY, QgsTriangle

t = QgsTriangle(
    QgsPointXY(0, 0),
    QgsPointXY(1, 0),
    QgsPointXY(0.5, 1),
)

g = QgsGeometry(t)

print(g.asWkt())

Tvorba vrstvy s memory driverem

  • data uložená v RAM paměti, po dobu běhu QGIS, nebo skriptu
from qgis.core import Qgis, QgsFields, QgsMemoryProviderUtils

layer = QgsMemoryProviderUtils.createMemoryLayer(
    "Sampling Table",
    fields=QgsFields(),
    geometryType=Qgis.WkbType.NoGeometry,
)

Tvorba vrstvy s memory driverem

from qgis.core import QgsVectorLayer

table = QgsVectorLayer(
    "NoGeometry?"
    "field=a:double&"
    "field=b:int",
    "source",
    "memory",
)

Řešení problému s potenciálním návratem None

  • velká část funkcí qgis může vracet buď objekt případně None, pokud operace neproběhne korektně
  • např. QgsVectorLayer() pokud se nepodaří načíst data ze zdroje
  • problém je, že mypy pak označuje velkou část kódu jako potenciálně problematickou
  • řešení je několik, nejsnazší a nejméně narušující je použití cast(...) z modulu typing

Ukázka

from qgis.core import QgsVectorLayer
from typing import cast

layer = QgsVectorLayer("cesta/k/souboru.gpkg", "vrstva", "ogr")

print(layer.featureCount())  # mypy hlásí chybu, layer může být None

layer = cast(QgsVectorLayer, layer)  

print(layer.featureCount())  # zde mypy už nehlásí chybu

Dotazy?