PostGIS-Interpolations- / Triangulationsoptionen

12

Ich bin auf der Suche nach einem Open-Source-Workflow zur Automatisierung der DEM-Erstellung. Wir haben eine Reihe von Sites, die in ein LIDAR-Dataset fallen. Wir möchten standortspezifische DEMs für jeden Standort erstellen und den Prozess automatisieren.

Bisher haben wir automatisiert:

  • Laden von Daten in PostGIS (LIDAR und Standortpunkte)
  • Erstellung von Standortgrenzen (Kombination aus st_buffer und st_envelope)

Jetzt suchen wir nach Optionen zum Interpolieren der Datenteilmengen an jedem Standort und zum Exportieren in Surfer7-Rasterdateien.

Derzeit werden die Daten in der PostGIS-Datenbank mit QuantumGIS abgefragt, in CSV-Dateien exportiert und anschließend manuell in Surfer7 importiert und in Grids umgewandelt. Hoffentlich können wir das auch automatisieren.

Für diesen Anwendungsfall möchten wir unsere LIDAR-Daten in ein DEM triangulieren. Andere mögliche Fälle, die wir uns vorstellen können, sind unterschiedliche Interpolationsmethoden. Wenn es also Optionen für inversen Abstand und Kriging gibt, sind wir interessiert!

Dies ist für uns eine Lernübung - wir arbeiten in kleinen Schritten!

tomtomnz
quelle
2
Ich würde sagen, dass Sie dem Entwicklungsteam voraus sind, das über die pointcloudEinführung des Datentyps in PostGIS nachgedacht hat, und hoffentlich bald daran arbeiten werden . In der Zwischenzeit gibt es ein points2grid- Dienstprogramm für PDAL , das interessant aussieht, das ich jedoch nicht verwendet habe.
Mike T

Antworten:

2

Sie können diese Automatisierung mit SAGA unter Windows oder Linux durchführen. SAGA kann große Punktdatensätze verarbeiten, die im SAGA-Format als pointcloud-Format (spc) importiert wurden. Dann haben Sie Funktionen zur Interpolation oder Triangulation und exportieren auch nach Surfer. Die Automatisierung erfolgt mit dem Befehlstool SAGA .

Niculita Mihai
quelle
In CVS importieren, können Sie dies mit dem psql-Befehl COPY (SELECT foo, bar FROM whatever) in bashscript setzen. TO '/tmp/ouput.csv' DELIMITER ',' CSV HEADER; Damit können Sie QGIS aus der
Toolkette entfernen
2

Wenn Sie die ursprünglichen LIDAR-Daten in .las / laz haben, können Sie lastools verwenden, http://www.cs.unc.edu/~isenburg/lastools/ , um mit großen LiDAR-Daten umzugehen. Es gibt ein las2tin und ein las2dem, die so aussehen, als würden sie direkt funktionieren, aber ich habe sie nicht benutzt. Ich habe das Folgende verwendet, um einen Satz von .laz-Dateien in einem Raster-DEM in GRASS zusammenzusetzen:

cat *.laz |~/Downloads/SRC/lastools/bin/las2txt -stdin  -parse xyztc -header pound -otxt |invproj -v -f %.9f +init=epsg:2284 |sed -e 's// /' | r.in.xyz input=- output=lidarJunk method=min type=FCELL x=1 y=2 z=3 fs=\    zscale=1.0 percent=100 --overwrite  
# the sed and fs=\ bits are for dealing with reformatting spaces and tabs between invproj and r.in.xyz

Lastools lässt sich gut automatisieren und sieht so aus, als könnten die Daten für ein DEM oder eine TIN an Grenzen Ihrer Wahl gekürzt werden.

Dave X.
quelle