PostGIS enterpolasyon / nirengi seçenekleri


12

DEM yapısını otomatikleştirmek için açık kaynaklı bir iş akışı arıyorum. LIDAR veri kümesine giren bir dizi sitemiz var. Her site için siteye özgü DEM'ler oluşturmak istiyoruz ve süreci otomatikleştirmek istiyoruz.

Şimdiye kadar otomatikleştirdik:

  • PostGIS'e veri yükleme (LIDAR ve site noktaları)
  • site sınırlarının oluşturulması (st_buffer ve st_envelope kombinasyonu)

Şimdi her sitedeki veri alt kümelerini enterpolasyona sokma ve bunları Surfer7 ızgara dosyalarına aktarma seçeneklerini arıyoruz.

Şu anda PostGIS veritabanındaki verileri QuantumGIS kullanarak sorgulıyoruz, csv dosyalarına dışa aktarıyoruz, sonra Surfer7'deki verileri manuel olarak içe aktarıyoruz ve ızgaralıyoruz. Umarım bunu da otomatikleştirebiliriz.

Bu nedenle, bu kullanım durumu için, lidar verilerimizi bir DEM'ye üçgenlemek istiyoruz. Düşünebileceğimiz diğer potansiyel durumlar farklı enterpolasyon yöntemlerini içerir - bu yüzden ters mesafe ve kriging için seçenekler varsa - ilgileniyoruz!

Bu bizim için bir öğrenme egzersizi - bebek adımlarıyla çalışıyoruz!


2
pointcloudVeri türünü PostGIS'e getirmeyi tartışan ve umarım yakında üzerinde çalışacak olan geliştirme ekibinin önünde olduğunuzu söyleyebilirim . Bu süre içinde, orada olduğu points2grid için yarar Pdal görünüyor ilginç olduğunu, ama bunu kullanmadım.
Mike T

Yanıtlar:


2

Bu otomasyonu SAGA ile Windows veya Linux'ta yapabilirsiniz. SAGA, SAGA biçiminde pointcloud biçimi (spc) olarak içe aktarılan büyük nokta veri kümelerini işleyebilir. Daha sonra enterpolasyon veya nirengi için fonksiyonlarınız var ve ayrıca Surfer'e ihracat yapıyorsunuz. Otomasyon, SAGA komut aracı kullanılarak gerçekleştirilir .


CVS'ye aktarın, bunu psql komutu COPY (SELECT foo, FROM olursa olsun) ile bashscript'e '/tmp/ouput.csv' DELIMITER ',' CSV HEADER; Bununla QGIS'i alet zincirinden düşürebilirsiniz
simplexio

2

Orijinal LIDAR verileriniz .las / laz'da varsa, büyük LiDAR verileriyle başa çıkmak için lastools, http://www.cs.unc.edu/~isenburg/lastools/ kullanabilirsiniz . Doğrudan çalışacak gibi görünen bir las2tin ve las2dem var, ama onları kullanmadım. GRASS içinde bir raster DEM içine bir .laz dosyaları kümesi oluşturmak için aşağıdakileri kullandım:

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 iyi bir otomasyona sahip olacak ve bir DEM veya TIN için verileri seçtiğiniz sınırlara kırpabilecek gibi görünüyor.

Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.