Opzioni di interpolazione / triangolazione PostGIS

12

Sto cercando un flusso di lavoro open source per automatizzare la costruzione di DEM. Abbiamo una serie di siti che rientrano in un set di dati LIDAR. Vogliamo creare DEM specifici per ogni sito e stiamo cercando di automatizzare il processo.

Finora abbiamo automatizzato:

caricamento dei dati in PostGIS (LIDAR e punti del sito)
creazione dei confini del sito (combinazione di st_buffer e st_envelope)

Ora stiamo cercando opzioni per interpolare i sottoinsiemi di dati in ciascun sito ed esportarli in file di griglia Surfer7.

Attualmente stiamo eseguendo una query sui dati nel database PostGIS utilizzando QuantumGIS, esportando in file CSV, quindi importando e eseguendo manualmente la griglia dei dati in Surfer7. Speriamo di poter automatizzare anche questo.

Quindi, per questo caso d'uso, vorremmo triangolare i nostri dati lidar su un DEM. Altri potenziali casi a cui possiamo pensare riguardano diversi metodi di interpolazione - quindi se ci sono opzioni per la distanza inversa e il kriging - siamo interessati!

Questo è per noi un esercizio di apprendimento: stiamo lavorando a piccoli passi!

— tomtomnz
fonte

2

Direi che sei in vantaggio rispetto al team di sviluppo, che ha discusso di portare il pointcloudtipo di dati su PostGIS e speriamo che ci lavorerà presto . Nel frattempo, c'è un points2grid utility per PDAL che sembra interessante, ma io non l'ho usato.

— Mike T,

2

È possibile eseguire questa automazione con SAGA , in Windows o in Linux. SAGA è in grado di gestire set di dati di grandi dimensioni, importati in formato SAGA come formato pointcloud (spc). Quindi hai funzioni per l'interpolazione o la triangolazione, ed esporta anche su Surfer. L'automazione viene eseguita utilizzando lo strumento di comando SAGA .

— Niculita Mihai
fonte

Importa in CVS, puoi metterlo su bashscript con il comando psql COPIA (SELEZIONA foo, barra DA qualunque cosa) A '/tmp/ouput.csv' DELIMITER ',' CSV HEADER; Con ciò puoi eliminare QGIS dalla catena degli strumenti

— simplexio il

2

Se disponi dei dati LIDAR originali in .las / laz, potresti utilizzare lastools, http://www.cs.unc.edu/~isenburg/lastools/ per gestire i grandi dati LiDAR. C'è un las2tin e las2dem che sembrano funzionare direttamente, ma non li ho usati. Ho usato il seguito per assemblare un set di file .laz in un DEM raster in GRASS:

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 si automatizzerebbe bene e sembrerebbe che potrebbe tagliare i dati per un DEM o TIN ai limiti della tua scelta.

— Dave X
fonte