QGIS + Python Teil 1

――Da ich gerade angefangen habe, GIS (QGIS) zu studieren, gibt es viele Zweifel an meinem Verständnis.

Umgebung zur Überprüfung der Ausführungsoperation

• Windows 7 64bit
• QGIS 2.8.11

Vorbereitung

--Wählen Sie "Python Console" aus "Plugins" in der QGIS-Menüleiste. --Wählen Sie die Editoranzeige über das Python-Konsolensymbol

• Verwenden Sie danach nicht die interaktive Konsole, sondern schreiben Sie den Code und führen Sie ihn im Editor aus.

Überprüfen Sie die Python-Version

test.py

import platform
print(platform.python_version())

• Das Ausführungsergebnis wird auf der Konsole angezeigt → 2.7.5 ――Es scheint Python der 2.x-Serie zu sein

Überprüfen Sie den Python-Ausführungspfad

test.py

import os
print(os.getcwd())

• Ausführungsergebnis → C: \ PROGRA ~ 1 \ QGIS2 ~ 1.18 \ bin
  • C:\Program Files\QGIS 2.18\bin

Vektorebene hinzufügen

test.py

filePath='C:/Users/xxxx/Documents/GIS DataBase/PyTest/lay01.shp'
iface.addVectorLayer(filePath,'test','ogr')

• Erstes Argument: Dateipfad
• Der Dateipfad kann auch mit "Schrägstrich" anstelle von "" beschrieben werden.
• Achten Sie bei Verwendung von "" auf Escape-Sequenzen.
• Zweites Argument: Der Name, den Sie der Ebene geben möchten
• Im Fall des obigen Codes lautet der Layername in QGIS "test lay01 Polygon".
• Zweites Argument + SP + Dateiname + SP + Typ
  • 
• Drittes Argument: Unbekannt, lassen Sie es als "ogr". providerKey
• Mach es nicht "org".

Holen Sie sich aktive Ebene

--Wählen Sie die Ebene im Voraus mit einer Maus oder ähnlichem.

test.py

layer=iface.activeLayer()
print('\n'.join(dir(layer)))

• Mit dir (xxx) können Sie eine Liste der Mitglieder von xxx-Objekten (Liste der verfügbaren Methoden usw.) abrufen
• geometryType()、getFeatures()
• Weitere Informationen finden Sie unter QgsVectorLayer-Klassenreferenz.

Untersuchen Sie den Geometrietyp einer Ebene

――Es scheint, dass Punkte (Punkte), Linien (Polylinien, Linien) und Polygone (Flächen) nicht in der Vektorebene von QGIS gemischt werden können.

• Daher kann der Typ der Vektorebene (wahrscheinlich) Punkt, Linie oder Polygon sein.
  • 

--Punkte, Linien und Polygone werden gemeinsam als Geometrie bezeichnet.

test.py

layer=iface.activeLayer()
type = layer.geometryType()
if type == QGis.Point :
    print('Point')
elif type == QGis.Line :
    print('Line')
elif type == QGis.Polygon:
    print('Polygon')

Ausgabe der Attribute des ausgewählten Features (ausgewähltes Polygon)

• Angenommen, das Feature verfügt über eine Attributspalte mit dem Namen "id".

  • 

• Angenommen, eine Funktion wird einzeln oder mehrmals ausgewählt.

test.py

layer=iface.activeLayer()
features = layer.selectedFeatures()

for feature in features:
    print(feature['id'])

--Verwenden Sie "layer.getFeatures ()", wenn Sie alle Features auf dem aktiven Layer mit oder ohne Auswahl als Ziel festlegen möchten.

test.py

layer=iface.activeLayer()
features = layer.getFeatures()

for feature in features:
    print(feature['id'])

Geben Sie die Scheitelpunktkoordinaten des Features (Polygon) aus.

test.py

# -*- coding: utf-8 -*- 

layer=iface.activeLayer()

if layer is None :
    print(u'Es gibt keine aktive Ebene.')

if layer.geometryType() != QGis.Polygon:
    print(u'Der Geometrietyp der aktiven Ebene ist kein Polygon.')

print(u'Aktive Schicht{0}Führt den Prozess für aus.'.format(layer.name()))
print('')

for feature in layer.getFeatures():

    print(u'ID={0}'.format(feature['id']))
    if feature.geometry().isMultipart():
        polygons = feature.geometry().asMultiPolygon()
    else:
        polygons = [ feature.geometry().asPolygon() ]
                
    for polygon in polygons:
        for vertices in polygon:
            for vertex in vertices:
                print('X={0:.9f}  Y={1:.9f}'.format(vertex.x(),vertex.y()))
 
    print('')

• Der obige "Vertex" ist eine Instanz von "class'qgis._core.QgsPoint"
• Einzelheiten finden Sie unter QgsPoint-Klassenreferenz.

Holen Sie sich den Begrenzungsrahmen einer Ebene

• Holen Sie sich den Begrenzungsrahmen für alle Features in der Ebene
• Der Begrenzungsrahmen ist die Klasse "QgsRectangle"
• Einzelheiten finden Sie unter QgsRectangle-Klassenreferenz.

test.py

# -*- coding: utf-8 -*- 
layer=iface.activeLayer()

if layer is None :
    print(u'Es gibt keine aktive Ebene.')

box = layer.extent()
print('xMinimum={0:.9f}'.format(box.xMinimum()))
print('xMaximum={0:.9f}'.format(box.xMaximum()))
print('yMinimum={0:.9f}'.format(box.yMinimum()))
print('yMaximum={0:.9f}'.format(box.yMaximum()))
print('width={0:.9f}'.format(box.width()))
print('height={0:.9f}'.format(box.height()))

Referenzmaterial

• Einführung in die QGIS-Programmierung 2016 Osaka Edition
• GIS Stack Exchange
• QGIS API Documentation 2.18
• Getting Started With Python Programming