gdal_grid – Online in der Cloud

PROGRAMM:

NAME/FUNKTION

gdal_grid – gdal_grid erstellt ein regelmäßiges Raster aus den verstreuten Daten

ZUSAMMENFASSUNG

gdal_grid [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/
CInt16/CInt32/CFloat32/CFloat64}]
[-of format] [-co "NAME=VALUE"]
[-zfield Feldname] [-z_increase Erhöhungswert] [-z_multiply Multiplikationswert]
[-a_srs srs_def] [-spat xmin ymin xmax ymax]
[-clipsrc |WKT|datasource|spat_extent]
[-clipsrcsql sql_statement] [-clipsrclayer Ebene]
[-clipsrcwhere Ausdruck]
[-l Layername]* [-where Ausdruck] [-sql select_statement]
[-txe xmin xmax] [-tye ymin ymax] [-outsize xsize ysize]
[-a Algorithmus[:parameter1=wert1]*] [-q]

BESCHREIBUNG

Dieses Programm erstellt ein regelmäßiges Raster (Raster) aus den verstreuten Daten, die aus dem OGR gelesen werden
Datenquelle. Eingabedaten werden interpoliert, um Rasterknoten mit Werten zu füllen, die Sie auswählen können
aus verschiedenen Interpolationsmethoden.

Ab GDAL 1.10 ist es möglich, Folgendes festzulegen GDAL_NUM_THREADS Konfigurationsmöglichkeit
um die Verarbeitung zu parallelisieren. Der anzugebende Wert ist die Anzahl der Arbeitsthreads oder
ALL_CPUS um alle Kerne/CPUs des Computers zu nutzen.

-ot tippe:
Damit die Ausgabebänder den angegebenen Datentyp haben.

-oder Format:
Wählen Sie das Ausgabeformat aus. Der Standardwert ist GeoTIFF (GTiff). Verwenden Sie den kurzen Formatnamen.

-txe xMin xmax:
Legen Sie georeferenzierte X-Ausdehnungen der zu erstellenden Ausgabedatei fest.

-tye ymin ymax:
Legen Sie georeferenzierte Y-Ausdehnungen der zu erstellenden Ausgabedatei fest.

-Übergröße xgröße ygröße:
Legen Sie die Größe der Ausgabedatei in Pixel und Zeilen fest.

-a_srs srs_def:
Überschreiben Sie die Projektion für die Ausgabedatei. Der srs_def kann irgendein Übliches sein
GDAL/OGR-Formulare, vollständiges WKT, PROJ.4, EPSG:n oder eine Datei, die das WKT enthält.

-zfield Feldname:
Identifiziert ein Attributfeld für die Features, das zum Abrufen eines Z-Werts verwendet werden soll. Das
Der Wert überschreibt den aus dem Feature-Geometriedatensatz gelesenen Z-Wert (natürlich, wenn Sie einen Z-Wert haben).
Geben Sie einen Wert in der Geometrie an, andernfalls haben Sie keine Wahl und sollten einen Feldnamen angeben
enthält den Z-Wert).

-z_increase erhöhen_Wert:
Ergänzung zum Attributfeld für die Features, aus denen ein Z-Wert ermittelt werden soll. Der
Die Addition sollte dieselbe Einheit wie der Z-Wert haben. Der Ergebniswert ist Z-Wert + Z
Wert steigern. Der Standardwert ist 0.

-z_multiply multiply_value:
Dies ist das Multiplikationsverhältnis für das Z-Feld. Dies kann z. B. für die Umstellung vom Fuß auf den Fuß verwendet werden
Meter oder von der Höhe bis zur Tiefe. Der Ergebniswert ist (Z-Wert + Z-Erhöhung).
Wert) * Z-Multiplikationswert. Der Standardwert ist 1.

-a [Algorithmus[:Parameter1=Wert1][:Parameter2=Wert2]...]:
Legen Sie den Namen des Interpolationsalgorithmus oder der Datenmetrik und (optional) seine Parameter fest.
Weitere Informationen finden Sie auch in den INTERPOLATION ALGORITHMEN und DATEN METRIK Abschnitte zur weiteren Diskussion
Verfügbare Optionen.

-spuckte xMin ymin xmax ymax:
Fügt einen räumlichen Filter hinzu, um nur Features auszuwählen, die im Begrenzungsrahmen enthalten sind
beschrieben durch (xmin, ymin) - (xmax, ymax).

-clipsrc [xmin ymin xmax ymax]|WKT|datasource|spat_extent:
Fügt einen räumlichen Filter hinzu, um nur Features auszuwählen, die innerhalb der angegebenen Begrenzung enthalten sind
Box (ausgedrückt im Quell-SRS), WKT-Geometrie (POLYGON oder MULTIPOLYGON), aus a
Datenquelle oder auf die räumliche Ausdehnung der -spuckte Option, wenn Sie die verwenden spat_extent
Stichwort. Wenn Sie eine Datenquelle angeben, möchten Sie diese im Allgemeinen verwenden
Kombination der -clipsrclayer, -clipsrcwhere or -clipsrcsql nach.

-clipsrcsql sql_statement:
Wählen Sie die gewünschten Geometrien stattdessen mithilfe einer SQL-Abfrage aus.

-clipsrclayer Layername:
Wählen Sie die benannte Ebene aus der Quellclip-Datenquelle aus.

-clipsrcwhere Ausdruck:
Beschränken Sie gewünschte Geometrien basierend auf der Attributabfrage.

-l Layername:
Gibt die Layer aus der Datenquelle an, die für Eingabe-Features verwendet werden. Mai
mehrfach angegeben werden, jedoch mindestens ein Layername oder ein -sql Option muss sein
spezifiziert.

-wo Ausdruck:
Ein optionaler SQL-WHERE-Abfrageausdruck, der auf die Auswahl von Features angewendet werden soll
Prozess aus der/den Eingabeschicht(en).

-sql select_statement:
Eine SQL-Anweisung, die anhand der Datenquelle ausgewertet werden soll, um eine virtuelle Schicht davon zu erstellen
zu verarbeitende Merkmale.

-co 'NAME=WERT':
Übergibt eine Erstellungsoption an den Ausgabeformattreiber. Mehrere -co Optionen können sein
aufgeführt. Informationen zu rechtlichen Erstellungsoptionen für jedes Format finden Sie in der formatspezifischen Dokumentation.

-Q:
Unterdrücken Sie den Fortschrittsmonitor und andere fehlerfreie Ausgaben.

src_datasource:
Jede von OGR unterstützte lesbare Datenquelle.

dst_filename:
Die von GDAL unterstützte Ausgabedatei.

INTERPOLATION ALGORITHMEN

Es stehen zahlreiche Interpolationsalgorithmen zur Auswahl.

Invdist
Inverse Distanz zu einer Potenz. Dies ist der Standardalgorithmus. Es hat folgende Parameter:

Werkzeuge:
Gewichtungsleistung (Standard 2.0).

Glättung:
Glättungsparameter (Standard 0.0).

radius1:
Der erste Radius (X-Achse, wenn der Drehwinkel 0 ist) der Suchellipse. Stellen Sie diesen Parameter ein
auf Null, um ein ganzes Punktarray zu verwenden. Der Standardwert ist 0.0.

radius2:
Der zweite Radius (Y-Achse, wenn der Drehwinkel 0 ist) der Suchellipse. Stellen Sie dies ein
Parameter auf Null setzen, um ein ganzes Punktarray zu verwenden. Der Standardwert ist 0.0.

Winkel:
Winkel der Suchellipsendrehung in Grad (gegen den Uhrzeigersinn, Standard 0.0).

max_points:
Maximale Anzahl der zu verwendenden Datenpunkte. Suchen Sie nicht nach mehr Punkten als dieser Anzahl.
Dies wird nur verwendet, wenn die Suchellipse eingestellt ist (beide Radien sind ungleich Null). Null bedeutet das
Alle gefundenen Punkte sollten verwendet werden. Der Standardwert ist 0.

min_points:
Mindestanzahl der zu verwendenden Datenpunkte. Wenn weniger Punkte vorhanden sind, wurde der Gitterknoten gefunden
gilt als leer und wird mit der NODATA-Markierung gefüllt. Dies wird nur bei der Suche verwendet
Ellipse ist eingestellt (beide Radien sind ungleich Null). Der Standardwert ist 0.

keine Daten:
NODATA-Marker zum Füllen leerer Punkte (Standard 0.0).

durchschnittlich
Algorithmus für gleitenden Durchschnitt. Es hat folgende Parameter:

radius1:
Der erste Radius (X-Achse, wenn der Drehwinkel 0 ist) der Suchellipse. Stellen Sie diesen Parameter ein
auf Null, um ein ganzes Punktarray zu verwenden. Der Standardwert ist 0.0.

radius2:
Der zweite Radius (Y-Achse, wenn der Drehwinkel 0 ist) der Suchellipse. Stellen Sie dies ein
Parameter auf Null setzen, um ein ganzes Punktarray zu verwenden. Der Standardwert ist 0.0.

Winkel:
Winkel der Suchellipsendrehung in Grad (gegen den Uhrzeigersinn, Standard 0.0).

min_points:
Mindestanzahl der zu verwendenden Datenpunkte. Wenn weniger Punkte vorhanden sind, wurde der Gitterknoten gefunden
gilt als leer und wird mit der NODATA-Markierung gefüllt. Der Standardwert ist 0.

keine Daten:
NODATA-Marker zum Füllen leerer Punkte (Standard 0.0).

Beachten Sie, dass es wichtig ist, die Suchellipse für die Methode des gleitenden Durchschnitts festzulegen. Es ist ein Fenster
Dies wird bei der Berechnung der Gitterknotenwerte gemittelt.

nächste
Algorithmus für den nächsten Nachbarn. Es hat folgende Parameter:

radius1:
Der erste Radius (X-Achse, wenn der Drehwinkel 0 ist) der Suchellipse. Stellen Sie diesen Parameter ein
auf Null, um ein ganzes Punktarray zu verwenden. Der Standardwert ist 0.0.

radius2:
Der zweite Radius (Y-Achse, wenn der Drehwinkel 0 ist) der Suchellipse. Stellen Sie dies ein
Parameter auf Null setzen, um ein ganzes Punktarray zu verwenden. Der Standardwert ist 0.0.

Winkel:
Winkel der Suchellipsendrehung in Grad (gegen den Uhrzeigersinn, Standard 0.0).

keine Daten:
NODATA-Marker zum Füllen leerer Punkte (Standard 0.0).

DATEN METRIK

Neben der Interpolationsfunktionalität gdal_grid kann zur Berechnung einiger Datenmetriken verwendet werden
unter Verwendung der angegebenen Fenster- und Ausgabegittergeometrie. Diese Kennzahlen sind:

Minimum:
Minimaler Wert, der in der Suchellipse des Gitterknotens gefunden wurde.

maximal:
Maximaler Wert, der in der Suchellipse des Gitterknotens gefunden wurde.

Angebot:
Ein Unterschied zwischen den minimalen und maximalen Werten, die in der Gitterknoten-Suchellipse gefunden wurden.

zählen:
Eine Reihe von Datenpunkten, die in der Suchellipse des Gitterknotens gefunden wurden.

durchschnittliche Entfernung:
Ein durchschnittlicher Abstand zwischen dem Gitterknoten (Mittelpunkt der Suchellipse) und allen
Die in der Gitterknoten-Suchellipse gefundenen Datenpunkte.

durchschnittliche_Entfernung_Punkte:
Ein durchschnittlicher Abstand zwischen den in der Gitterknoten-Suchellipse gefundenen Datenpunkten. Der
Der Abstand zwischen jedem Punktpaar innerhalb der Ellipse wird berechnet und der Durchschnitt aller Punkte gebildet
Die Entfernungen werden als Gitterknotenwert festgelegt.

Alle Metriken verfügen über die gleichen Optionen:

radius1:
Der erste Radius (X-Achse, wenn der Drehwinkel 0 ist) der Suchellipse. Stellen Sie diesen Parameter ein
auf Null, um ein ganzes Punktarray zu verwenden. Der Standardwert ist 0.0.

radius2:
Der zweite Radius (Y-Achse, wenn der Drehwinkel 0 ist) der Suchellipse. Stellen Sie dies ein
Parameter auf Null setzen, um ein ganzes Punktarray zu verwenden. Der Standardwert ist 0.0.

Winkel:
Winkel der Suchellipsendrehung in Grad (gegen den Uhrzeigersinn, Standard 0.0).

min_points:
Mindestanzahl der zu verwendenden Datenpunkte. Wenn weniger Punkte vorhanden sind, wurde der Gitterknoten gefunden
gilt als leer und wird mit der NODATA-Markierung gefüllt. Dies wird nur bei der Suche verwendet
Ellipse ist eingestellt (beide Radien sind ungleich Null). Der Standardwert ist 0.

keine Daten:
NODATA-Marker zum Füllen leerer Punkte (Standard 0.0).

LESEN ABSATZ GETRENNT WERTE

Oft haben Sie eine Textdatei mit einer Liste von durch Kommas getrennten XYZ-Werten, mit denen Sie arbeiten können (so
sogenannte CSV-Datei). Sie können diese Art von Datenquelle problemlos verwenden gdal_grid. Alles was du brauchst
Erstellen Sie einen virtuellen Datensatz-Header (VRT) für Ihre CSV-Datei und verwenden Sie ihn als Eingabedatenquelle
für gdal_grid. Einzelheiten zum VRT-Format finden Sie auf der Beschreibungsseite zum virtuellen Format.

Hier ist ein kleines Beispiel. Lassen Sie uns eine CSV-Datei namens haben dem.csv mit

Ostwert, Nordwert, Höhe
86943.4,891957,139.13
87124.3,892075,135.01
86962.4,892321,182.04
87077.6,891995,135.01

Für die oben genannten Daten erstellen wir dem.vrt Header mit folgendem Inhalt:



dem.csv
wkbPoint




Diese Beschreibung spezifiziert die sogenannte 2.5D-Geometrie mit den drei Koordinaten X, Y und Z. Z
Der Wert wird für die Interpolation verwendet. Jetzt können Sie es verwenden dem.vrt mit allen OGR-Programmen (Start
mit ogrinfo um zu testen, ob alles gut funktioniert). Die Datenquelle enthält eine einzelne Ebene
namens 'dem' gefüllt mit Punktmerkmalen, die aus Werten in der CSV-Datei erstellt wurden. Benutze dies
Mit dieser Technik können Sie CSV-Dateien mit mehr als drei Spalten verarbeiten, Spalten wechseln usw.

Wenn Ihre CSV-Datei keine Spaltenüberschriften enthält, kann dies wie folgt behandelt werden
Weg:



Die Beschreibungsseite für durch Kommas getrennte Werte enthält Details zum CSV-Format, das von unterstützt wird
GDAL/OGR.

BEISPIEL

Im Folgenden wird eine Raster-TIFF-Datei aus der in beschriebenen VRT-Datenquelle erstellt LESEN ABSATZ
GETRENNT WERTE Abschnitt unter Verwendung der inversen Distanz-zu-einer-Potenz-Methode. Werte zu
Die Interpolation wird aus dem Z-Wert des Geometriedatensatzes gelesen.

gdal_grid -a invdist:power=2.0:smoothing=1.0 -txe 85000 89000 -tye 894000 890000 -outsize 400 400 -of GTiff -ot Float64 -l dem dem.vrt dem.tiff

Der nächste Befehl macht dasselbe wie der vorherige, liest jedoch zu interpolierende Werte
aus dem mit angegebenen Attributfeld -zfield Option statt Geometriedatensatz. Also rein
In diesem Fall werden die X- und Y-Koordinaten aus der Geometrie und die Z-Koordinaten aus der Geometrie übernommen
'Elevation' Feld. GDAL_NUM_THREADS wird ebenfalls festgelegt, um die Berechnung zu parallelisieren.

gdal_grid -zfield "Elevation" -a invdist:power=2.0:smoothing=1.0 -txe 85000 89000 -tye 894000 890000 -outsize 400 400 -of GTiff -ot Float64 -l dem dem.vrt dem.tiff --config GDAL_NUM_THREADS ALL_CPUS

AUTOREN

Andrey Kisseljow dron@ak4719.spb.edu

Verwenden Sie gdal_grid online über die Dienste von onworks.net