Opsi interpolasi / triangulasi PostGIS

12

Saya mencari alur kerja open source untuk mengotomatisasi konstruksi DEM. Kami memiliki serangkaian situs yang termasuk dalam dataset LIDAR. Kami ingin membuat DEM khusus situs untuk setiap situs, dan kami ingin mengotomatiskan prosesnya.

Sejauh ini, kami telah otomatis:

pemuatan data ke dalam PostGIS (LIDAR dan poin-poin situs)
pembuatan batas situs (kombinasi st_buffer dan st_envelope)

Sekarang kami sedang mencari opsi untuk menginterpolasi subset data di setiap situs dan mengekspornya ke file grid Surfer7.

Saat ini, kami sedang menanyakan data dalam database PostGIS menggunakan QuantumGIS, mengekspor ke file csv, kemudian secara manual mengimpor dan memasukkan data ke dalam Surfer7. Semoga kita dapat mengotomatisasi ini juga.

Jadi untuk kasus penggunaan ini, kami ingin melakukan triangulasi data LIDAR kami ke DEM. Kasus potensial lain yang dapat kita pikirkan melibatkan metode interpolasi yang berbeda - jadi jika ada opsi untuk jarak terbalik dan kriging - kami tertarik!

Ini merupakan latihan pembelajaran bagi kami - kami sedang mengerjakan langkah kecil!

— tomtomnz
sumber

2

Saya katakan Anda berada di depan tim pengembangan, yang telah mendiskusikan membawa pointclouddatatype ke PostGIS, dan semoga akan segera mengatasinya . Sementara itu, ada utilitas points2grid untuk PDAL yang terlihat menarik, tetapi saya belum menggunakannya.

— Mike T

2

Anda dapat melakukan otomatisasi ini dengan SAGA , di Windows atau di Linux. SAGA dapat menangani set data titik besar, yang diimpor dalam format SAGA sebagai format pointcloud (spc). Kemudian Anda memiliki fungsi untuk interpolasi atau triangulasi, dan juga ekspor ke Surfer. Otomasi dilakukan dengan menggunakan alat perintah SAGA .

— Niculita Mihai
sumber

Impor ke CVS, Anda dapat menempatkan ini ke bashscript dengan perintah psql COPY (PILIH foo, bilah DARI apa pun) KE '/tmp/ouput.csv' DELIMITER ',' CSV HEADER; Dengan itu Anda dapat menjatuhkan QGIS dari rantai alat

— simplexio

2

Jika Anda memiliki data LIDAR asli di .las / laz, Anda dapat menggunakan lastools, http://www.cs.unc.edu/~isenburg/lastools/ untuk menangani data LiDAR yang besar. Ada las2tin dan las2dem yang terlihat seperti mereka akan bekerja secara langsung, tetapi saya belum menggunakannya. Saya telah menggunakan di bawah ini untuk merakit satu set file .laz menjadi DEM raster dalam 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 akan terotomasi dengan baik, dan sepertinya mereka dapat memotong data untuk DEM atau TIN hingga batas pilihan Anda.

— Dave X
sumber