Варіанти інтерполяції / триангуляції після ГІС

12

Я шукаю робочий процес із відкритим кодом для автоматизації побудови DEM. У нас є низка сайтів, які підпадають під набір даних LIDAR. Ми хочемо створити конкретні DEM для кожного сайту, і ми прагнемо автоматизувати процес.

Поки ми автоматизували:

завантаження даних у PostGIS (LIDAR та точки сайту)
створення меж сайту (поєднання st_buffer і st_envelope)

Зараз ми шукаємо варіанти інтерполяції підмножини даних на кожному сайті та експортування їх у сіткові файли Surfer7.

В даний час ми запитуємо дані в базі даних PostGIS за допомогою QuantumGIS, експортуємо в файли csv, а потім вручну імпортуємо та передаємо їх у сітку в Surfer7. Сподіваємось, ми також можемо це автоматизувати.

Отже, для цього випадку ми хочемо тріангулювати наші дані про лідар до DEM. Інші потенційні випадки, про які ми можемо подумати, включають різні методи інтерполяції - тому, якщо є варіанти зворотної відстані та кригінгу - нас цікавить!

Це дуже багато навчальних вправ для нас - ми працюємо в дитячих кроках!

— tomtomnz
джерело

2

Я б сказав, що ви випередили команду з розробки, яка обговорювала питання про те, pointcloudщоб привести тип даних до PostGIS, і, сподіваємось, незабаром буде працювати над цим . Тим часом є PD2- утиліта для PDAL, яка виглядає цікаво, але я не використовував її.

— Майк Т

2

Ви можете зробити цю автоматизацію за допомогою SAGA , в Windows або Linux. SAGA може обробляти набори даних великих точок, імпортованих у форматі SAGA у форматі pointcloud (spc). Тоді ви маєте функції для інтерполяції чи тріангуляції, а також експортуєте до Surfer. Автоматизація виконується за допомогою командного інструменту SAGA .

— Нікуліта Міхай
джерело

Імпортуючи в CVS, ви можете поставити це до bashscript командою psql COPY (SELECT foo, bar OD Any) TO '/tmp/ouput.csv' DELIMITER ',' CSV HEADER; Завдяки цьому ви можете скинути QGIS з ланцюжка інструментів

— simplexio

2

Якщо у вас є оригінальні дані LIDAR в .las / laz, ви можете використовувати lastools, http://www.cs.unc.edu/~isenburg/lastools/, щоб мати справу з великими даними LiDAR. Є las2tin та las2dem, схожі на те, що вони працюватимуть безпосередньо, але я не використовував їх. Нижче я використовував для збирання набору файлів .laz у растровій DEM в межах 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 добре автоматизуватиметься і виглядатиме так, ніби вони могли вирізати дані для DEM або TIN до меж на ваш вибір.

— Дейв X
джерело