Wichtig

Übersetzen ist eine Gemeinschaftsleistung Sie können mitmachen. Diese Seite ist aktuell zu 100.00% übersetzt.

11.3. Georeferenzierung

Der georefRun Georeferenzierer ist ein Werkzeug zum Ausrichten von nicht referenzierten Raster- oder Vektorlayern auf bekannte Koordinatensysteme unter Verwendung von sogenannten Passpunkte (oder auch „GCP“ von eng.: „Ground Control Points“). Er unterstützt den Export transformierter Raster (GeoTIFF) und Vektoren (Shapefiles, Geopakete usw.) oder das Schreiben begleitender Weltdateien (nur für Raster). Der grundlegende Ansatz für die Georeferenzierung eines Layers besteht darin, Punkte festzulegen, für die Sie die Koordinaten möglichst genau bestimmen können.

Funktionen

Icon

Zweck

Icon

Zweck

addRasterLayer

Raster öffnen

addOgrLayer

Vektor öffnen

start

Georeferenzierung starten

gdalScript

GDAL-Skript erzeugen

loadGCPpoints

Passpunkte laden

saveGCPPointsAs

Passpunkte speichern

transformSettings

Transformationseinstellungen

addGCPPoint

GCP-Punkt hinzufügen

deleteGCPPoint

GCP-Punkt löschen

moveGCPPoint

Passpunkt verschieben

pan

Verschieben

zoomIn

Hinein zoomen

zoomOut

Heraus zoomen

zoomToLayer

Auf den Layer zoomen

zoomLast

Zoom zurück

zoomNext

Zoom vor

linkGeorefToQGis

Georeferenzierung mit QGIS verbinden

linkQGisToGeoref

QGIS mit Georeferenzierung verbinden

fullHistogramStretch

Volle Histogrammstreckung

localHistogramStretch

Lokale Histogrammstreckung

Table Georeferenzierung 1: Georeferenzierfunktionen

11.3.1. Georeferenzierung von Rasterlayern

Als X- und Y-Koordinaten (Grad, Minuten, Sekunden (dd mm ss.ss), Dezimalgrad (dd.dd) oder projizierte Koordinaten (mmmm.mm)), die dem ausgewählten Punkt auf dem Bild entsprechen, können zwei alternative Verfahren verwendet werden:

  • Das Raster selbst liefert manchmal Kreuze mit Koordinaten, die auf das Bild „geschrieben“ sind. In diesem Fall können Sie die Koordinaten manuell eingeben.

  • Verwendung von bereits georeferenzierten Layern. Dies können entweder Vektor- oder Rasterdaten sein, die dieselben Objekte/Features enthalten, die Sie auf dem Bild haben, das Sie georeferenzieren möchten, und mit der Projektion, die Sie für Ihr Bild wünschen. In diesem Fall können Sie die Koordinaten eingeben, indem Sie auf den in der QGIS Kartenansicht geladenen Referenzdatensatz klicken.

Das übliche Verfahren zur Georeferenzierung eines Bildes umfasst die Auswahl mehrerer Punkte auf dem Raster, die Angabe ihrer Koordinaten und die Auswahl eines entsprechenden Transformationstyps. Auf der Grundlage der eingegebenen Parameter und Daten berechnet der Georeferencer die Parameter der Weltdatei. Je mehr Koordinaten Sie angeben, desto besser wird das Ergebnis sein.

Wählen Sie aus der Menüleiste Layer ► georefRun Georeferenzierung. Der Georeferenzierung-Dialog erscheint wie in Abb. 11.34 gezeigt.

Für dieses Beispiel verwenden wir ein Rasterkarte von South Dakota von SDGS. Es kann später zusammen mit den Daten aus der GRASS spearfish60 Location visualisiert werden. Sie können das Topo-Blatt hier herunterladen: https://grass.osgeo.org/sampledata/spearfish_toposheet.tar.gz.

../../../_images/georef.png

Abb. 11.34 Georeferenzierung Dialog

11.3.1.1. Eingabe von Passpunkten (GCPs)

  1. Um mit der Georeferenzierung eines nicht referenzierten Rasters zu beginnen, müssen wir es über die Schaltfläche addRasterLayer Raster öffnen… laden. Das Raster wird im Hauptarbeitsbereich des Dialogs angezeigt. Sobald das Raster geladen ist, können wir mit der Eingabe von Referenzpunkten beginnen.

  2. Fügen Sie mit der Schaltfläche addGCPPoint GCP-Punkt hinzufügen Punkte zum Hauptarbeitsbereich hinzu und geben Sie deren Koordinaten ein (siehe Abbildung Abb. 11.35). Für diesen Vorgang stehen Ihnen folgende Optionen zur Verfügung:

    • Klicken Sie auf einen Punkt im Rasterbild und geben Sie die X- und Y-Koordinaten zusammen mit dem KBS des Punktes manuell ein.

    • Klicken Sie auf einen Punkt im Rasterbild und wählen Sie die Schaltfläche pencil Aus Kartenansicht, um die X- und Y-Koordinaten mithilfe einer bereits in die QGIS Kartenansicht geladenen georeferenzierten Karte hinzuzufügen. Das KBS wird automatisch gesetzt.

    • Bei der Eingabe von GCPs aus der Kartenansicht haben Sie die Möglichkeit, das Georeferenzierungsfenster auszublenden, während Sie Punkte aus der Kartenansicht auswählen. Wenn das Kontrollkästchen checkbox Georeferenzierungsfensters automatisch verbergen aktiviert ist, wird nach dem Klicken auf pencil Aus Kartenansicht das Hauptfenster des Georeferenzierers ausgeblendet, bis ein Punkt auf der Kartenansicht hinzugefügt wird. Der Dialog Kartenkoordinaten eingeben bleibt geöffnet. Ist das Kontrollkästchen nicht markiert, bleiben beide Fenster geöffnet, während ein Punkt auf der Kartenansicht ausgewählt wird. Diese Option ist nur wirksam, wenn das Georeferenzierungsfenster nicht an der Hauptschnittstelle angedockt ist.

  3. Fahren Sie mit der Eingabe der Punkte fort. Sie sollten mindestens vier Punkte haben, und je mehr Koordinaten Sie angeben können, desto besser wird das Ergebnis sein. Es gibt zusätzliche Werkzeuge zum Zoomen und Schwenken des Arbeitsbereichs, um einen relevanten Satz von GCP-Punkten zu finden.

    Tipp

    Um nicht ständig zwischen den Schaltflächen pan Verschieben, addGCPPoint GCP-Punkt hinzufügen und moveGCPPoint GCP-Punkt verschieben wechseln zu müssen, können Sie die Pfeiltasten der Tastatur zum Verschieben und das Mausrad zum Skalieren der georeferenzierten Karte verwenden.

  4. Nachdem Sie einige Punkte angegeben haben, können Sie die Schaltflächen linkQGisToGeoref QGIS mit Georeferenzierung verbinden und/oder linkGeorefToQGis Georeferenzierung mit QGIS verbinden verwenden, um den Kartenausschnitt des QGIS-Hauptfensters an die aktuelle Ansicht in der Georeferenzierung anzupassen und/oder umgekehrt.

  5. Mit dem Werkzeug moveGCPPoint Passpunkt verschieben können Sie die GCPs sowohl in der Kartenansicht als auch im Georeferenzierungsfenster verschieben, wenn Sie sie korrigieren müssen.

../../../_images/choose_points.png

Abb. 11.35 Punkte zum Rasterbild hinzufügen

Die zur Karte hinzugefügten Punkte werden in einer separaten Textdatei ([Dateiname].points) gespeichert, normalerweise zusammen mit dem Rasterbild. Dies ermöglicht es uns, den Georeferencer zu einem späteren Zeitpunkt erneut zu öffnen und neue Punkte hinzuzufügen oder bestehende zu löschen, um das Ergebnis zu optimieren. Die Punktedatei enthält Werte in der Form: mapX, mapY, pixelX, pixelY. Sie können die Schaltflächen loadGCPpoints GCP-Punkte laden und saveGCPPointsAs GCP-Punkte speichern als verwenden, um die Dateien zu verwalten.

11.3.1.2. Definieren der Transformationseinstellungen

Nachdem Sie Ihre Passpunkte zum Rasterbild hinzugefügt haben, müssen Sie die Transformationseinstellungen für den Georeferenzierungsprozess festlegen.

../../../_images/transformation_settings.png

Abb. 11.36 Festlegen der Einstellungen für die Transformation der Georeferenzierung

Verfügbare Transformations-Algorithmen

Es stehen mehrere Transformationsalgorithmen zur Verfügung, die von der Art und Qualität der Eingabedaten, der Art und dem Ausmaß der geometrischen Verzerrung, die Sie für das Endergebnis in Kauf nehmen, und der Anzahl der Passpunkte (GCPs) abhängen.

Derzeit sind die folgenden Transformationstypen verfügbar:

  • Der Lineare Algorithmus wird zur Erstellung einer Weltdatei verwendet und unterscheidet sich von den anderen Algorithmen, da er die Rasterpixel nicht wirklich transformiert. Er ermöglicht die Positionierung (Verschiebung) des Bildes und eine einheitliche Skalierung, aber keine Drehung oder andere Transformationen. Dieser Algorithmus ist am besten geeignet, wenn Ihr Bild eine qualitativ gute Rasterkarte in einem bekannten KBS ist, aber nur die Georeferenzierungsinformationen fehlen. Es werden mindestens 2 GCPs benötigt.

  • Die Helmert-Transformation ermöglicht auch eine Drehung. Sie ist besonders nützlich, wenn es sich bei Ihrem Raster um eine qualitativ gute lokale Karte oder ein orthorektifiziertes Luftbild handelt, das jedoch nicht mit dem Raster in Ihrem KBS ausgerichtet ist. Es werden mindestens 2 GCPs benötigt.

  • Der Polynomial 1-Algorithmus ermöglicht eine allgemeinere affine Transformation, insbesondere auch eine gleichmäßige Scherung. Gerade Linien bleiben gerade (d. h. kollineare Punkte bleiben kollinear) und parallele Linien bleiben parallel. Dies ist besonders nützlich für die Georeferenzierung von Datenkartogrammen, die möglicherweise mit unterschiedlichen Bodenpixelgrößen in verschiedenen Richtungen aufgezeichnet (oder erfasst) wurden. Es sind mindestens 3 GCPs erforderlich.

  • Die Polynomial-Algorithmen 2-3 verwenden allgemeinere Polynome 2. oder 3. Grades statt nur affiner Transformationen. Dadurch können sie Krümmungen oder andere systematische Verzerrungen des Bildes berücksichtigen, z. B. fotografierte Karten mit gekrümmten Kanten. Es werden mindestens 6 (bzw. 10) GCPs benötigt. Winkel und lokaler Maßstab werden nicht beibehalten oder gleichmäßig über das Bild behandelt. Insbesondere können gerade Linien gekrümmt werden, und es kann zu erheblichen Verzerrungen an den Rändern oder weit entfernt von den GCPs kommen, wenn die an die Daten angepassten Polynome zu weit extrapoliert werden.

  • Der Projektiv Algorithmus verallgemeinert das Polynom 1 auf eine andere Art und Weise und erlaubt Transformationen, die eine zentrale Projektion zwischen zwei nicht parallelen Ebenen, dem Bild und der Kartenansicht, darstellen. Gerade Linien bleiben gerade, aber die Parallelität wird nicht beibehalten und der Maßstab des Bildes ändert sich entsprechend der Änderung der Perspektive. Diese Transformationsart eignet sich am besten für die Georeferenzierung von schrägen Fotos (und nicht von flachen Scans) von Karten guter Qualität oder von schrägen Luftbildern. Ein Minimum von 4 GCPs ist erforderlich.

  • Der Algorithmus Thin Plate Spline (TPS) schließlich „ verformt“ das Raster unter Verwendung mehrerer lokaler Polynome, um es an die angegebenen GCPs anzupassen, wobei die Oberflächenkrümmung insgesamt minimiert wird. Bereiche, die von den GCPs entfernt sind, werden in der Ausgabe verschoben, um die GCP-Anpassung zu ermöglichen, werden aber ansonsten nur minimal lokal deformiert. TPS ist am nützlichsten für die Georeferenzierung beschädigter, deformierter oder anderweitig leicht ungenauer Karten oder schlecht orthorektifizierter Luftbilder. Es ist auch nützlich für die ungefähre Georeferenzierung und implizite Reprojektion von Karten mit unbekanntem Projektionstyp oder unbekannten Parametern, bei denen jedoch ein regelmäßiges Raster oder ein dichter Satz von Ad-hoc-GCPs mit einer Referenzkarten-Layer abgeglichen werden kann. Technisch gesehen sind mindestens 10 GCPs erforderlich, in der Regel jedoch mehr, um erfolgreich zu sein.

Wenn bei allen Algorithmen mit Ausnahme von TPS mehr als die minimalen GCPs angegeben sind, werden die Parameter so angepasst, dass der Gesamtrestfehler minimiert wird. Dies ist hilfreich, um die Auswirkungen von Registrierungsfehlern, d. h. leichten Ungenauigkeiten bei Zeigerklicks oder getippten Koordinaten, oder anderen kleinen lokalen Bildverformungen zu minimieren. Ohne weitere GCPs zum Ausgleich könnten solche Fehler oder Verformungen zu erheblichen Verzerrungen führen, insbesondere in der Nähe der Ränder des georeferenzierten Bildes.␣ Wenn jedoch mehr als die minimalen GCPs angegeben werden, stimmen sie in der Ausgabe nur ungefähr überein. Im Gegensatz dazu passt TPS alle angegebenen GCPs genau an, kann aber zu erheblichen Verformungen zwischen nahegelegenen GCPs mit Registrierungsfehlern führen.

Definieren der Resampling-Methode

Welche Art der Neuabtastung Sie wählen, hängt wahrscheinlich von Ihren Eingabedaten und dem letztendlichen Ziel der Arbeit ab. Wenn Sie die statistischen Daten des Rasters nicht verändern wollen (abgesehen von der ungleichmäßigen geometrischen Skalierung bei Verwendung anderer Transformationen als Linear, Helmert oder Polynom 1), sollten Sie die Option „Nächster Nachbar“ wählen. Im Gegensatz dazu führt z. B. das „kubische Resampling“ in der Regel zu einem visuell glatteren Ergebnis.

Es ist möglich, zwischen fünf verschiedenen Resampling-Methoden zu wählen:

  1. Nächster Nachbar

  2. Bilinear (2x2 Kern)

  3. Kubisch (4x4 Kern)

  4. Kubisch B-Spline (4x4 Kern)

  5. Lanczos (6x6 Kern)

Definieren der Transformationseinstellungen

Für das georeferenzierte Ausgaberaster müssen mehrere Optionen festgelegt werden.

  • Das Kontrollkästchen checkbox Nur World-Datei erzeugen ist nur verfügbar, wenn Sie sich für die lineare Transformation entscheiden, da dies bedeutet, dass das Rasterbild tatsächlich nicht transformiert wird. In diesem Fall ist das Feld Ausgabedatei nicht aktiviert, da nur eine neue Weltdatei erstellt wird.

  • Für alle anderen Transformationsarten müssen Sie eine Ausgabedatei definieren. Standardmäßig wird eine neue Datei ([Dateiname]_modifiziert) im selben Ordner zusammen mit dem ursprünglichen Rasterbild erstellt.

  • In einem nächsten Schritt müssen Sie das Ziel-KBS (Koordinatenreferenzsystem) für das georeferenzierte Raster definieren (siehe Arbeiten mit Projektionen).

  • Bei Bedarf können Sie PDF-Karte erzeugen und auch PDF-Bericht erzeugen aktivieren. Der Bericht enthält Informationen über die verwendeten Transformationsparameter, ein Bild der Residuen und eine Liste mit allen GCPs und ihren RMS-Fehlern.

  • Außerdem können Sie das Kontrollkästchen checkbox Zielauflösung setzen aktivieren und die Pixelauflösung des Ausgaberasters festlegen. Die Standardeinstellung für die horizontale und vertikale Auflösung ist 1.

  • Die checkbox Falls nötig 0 für Transparenz verwenden-Option kann aktiviert werden, wenn Pixel mit dem Wert 0 transparent dargestellt werden sollen. In unserem Beispiel-Toposheet würden alle weißen Flächen transparent sein.

  • Die checkbox Passpunkte speichern speichert die GCP Punkte in einer Datei neben dem Ausgaberaster.

  • Schließlich lädt checkbox Nach Fertigstellung ins Projekt laden das Ausgaberaster automatisch in die QGIS Kartenansicht, wenn die Transformation abgeschlossen ist.

11.3.1.3. Durchführung der Transformation

Nachdem alle GCPs gesammelt und alle Transformationseinstellungen definiert wurden, drücken Sie einfach die Schaltfläche start :sup:` Georeferenzierung starten`, um das neue georeferenzierte Raster zu erstellen.

11.3.2. Georeferenzierung des Vektor-Layers

Die Georeferenzierung von Vektor-Layern funktioniert ähnlich wie die Raster-Georeferenzierung, jedoch werden statt Bildpixeln Vektor-Geometrien (Punkte, Linien oder Polygone) mit bekannten räumlichen Referenzen abgeglichen.

Das Standardverfahren beginnt genauso wie bei der Raster-Georeferenzierung: Öffnen Sie QGIS und fügen Sie der Kartenfläche einen Layer hinzu, der als Referenz dienen soll. Dies kann ein georeferenzierter Raster- oder Vektor-Layer oder ein WMS-Layer sein.

Öffnen Sie das Dialogfeld „“Georeferenzierung“ über Layer ► georefRun Georeferenzierung.

Beginnen Sie mit der Georeferenzierung, indem Sie die folgenden Schritte ausführen (in diesem Beispiel verwenden wir eine nicht referenzierte Datei „alaska.shp”):

../../../_images/vector_georeferencer_dialog.png

Abb. 11.37 Dialogfeld Vektor- Georeferenzierung

  1. Laden Sie den nicht referenzierten Vektor-Layer mit der Schaltfläche addOgrLayer. Der Vektor-Layer wird im Hauptarbeitsbereich des Dialogfelds angezeigt.

  2. Verwenden Sie die Schaltfläche addGCPPoint GCP-Punkt hinzufügen, um einen Punkt zum Arbeitsbereich hinzuzufügen. Geben Sie dessen Koordinaten manuell ein und legen Sie das KBS fest, oder klicken Sie auf die Schaltfläche pencil Aus Kartenansicht, um die Koordinaten aus einem georeferenzierten Layers in der Hauptkartenansicht von QGIS auszuwählen. In diesem Fall wird das KBS automatisch festgelegt

  3. Definieren Sie die Transformationseinstellungen:

    • Wählen Sie den Transformationstyp aus. Die Transformationsalgorithmen sind dieselben wie bei der Raster-Georeferenzierung. Weitere Informationen finden Sie unter Definieren der Transformationseinstellungen.

    • Definieren Sie das Ziel-KBS ( Koordinatenbezugssystem) für den georeferenzierten Vektor (siehe Arbeiten mit Projektionen).

    • Legen Sie das Ausgabeformat und den Pfad der Datei fest (z. B. GeoPackage, Shapefile). Standardmäßig wird eine neue Datei mit der Endung „_modified“ im selben Ordner wie die ursprüngliche Vektordatei erstellt.

    • Optional können Sie PDF-Karte erstellen und PDF-Bericht erstellen aktivieren. Der Bericht enthält Transformationsparameter, GCP-Residuen und eine Zusammenfassung der RMS-Fehler.

    • Aktivieren Sie das Kontrollkästchen „GCP-Punkte speichern“, um die GCPs zusammen mit dem Ausgabevektor-Layer in einer Datei zu speichern.

    • Aktivieren Sie das Kontrollkästchen „Nach Fertigstellung ins Projekt laden“, um das Ergebnis direkt zur Kartenfläche hinzuzufügen.

  4. Klicken Sie auf start Georeferenzierung starten, um die Transformation auszuführen und den georeferenzierten Vektor-Layer zu generieren.

11.3.3. Layer-Eigenschaften anzeigen und anpassen

Durch Klicken auf die Option Quelleneigenschaften im Menü Einstellungen wird je nach Art des Layers, den Sie georeferenzieren, entweder Rasterlayer-Eigenschaften oder Vektorlayer-Eigenschaften geöffnet.

11.3.4. Konfiguration des Georeferenzierers

Sie können das Verhalten der Georeferenzierung unter Einstellungen ► Georeferenzierung konfigurieren anpassen (oder das Tastaturkürzel Strg+P verwenden).

  • Unter Punktspitze können Sie die Anzeige von GCP-IDs und X/Y-Koordinaten sowohl im Georeferenzierungsfenster als auch in der Kartenansicht umschalten.

  • Resteinheiten steuert, ob Resteinheiten in Pixeln oder Karteneinheiten angegeben werden

  • PDF-Bericht ermöglicht es Ihnen, die Randgröße in mm für den Berichtsexport festzulegen

  • PDF-Karte ermöglicht die Auswahl eines Papierformats für den Export der Karte

  • Schließlich können Sie das Kontrollkästchen Georeferenzierungsfenster docken aktivieren. Dadurch wird das Georeferenzierungsfenster an das Hauptfenster von QGIS angedockt, anstatt es als separates Fenster anzuzeigen, das minimiert werden kann.