Opcje interpolacji / triangulacji PostGIS

12

Szukam przepływu pracy typu open source do automatyzacji budowy DEM. Mamy szereg witryn należących do zestawu danych LIDAR. Chcemy stworzyć DEM dla każdej witryny i chcemy zautomatyzować ten proces.

Do tej pory zautomatyzowaliśmy:

  • ładowanie danych do PostGIS (LIDAR i punkty strony)
  • tworzenie granic serwisu (połączenie st_buffer i st_envelope)

Teraz szukamy opcji interpolacji podzbiorów danych w każdej witrynie i eksportowania ich do plików siatki Surfer7.

Obecnie przeszukujemy dane w bazie danych PostGIS za pomocą QuantumGIS, eksportujemy do plików csv, a następnie ręcznie importujemy i gridujemy dane w Surfer7. Mamy nadzieję, że możemy to również zautomatyzować.

W tym przypadku chcielibyśmy triangulować nasze dane lidar do DEM. Inne potencjalne przypadki, o których możemy pomyśleć, obejmują różne metody interpolacji - więc jeśli istnieją opcje odwrotnej odległości i krigingu - jesteśmy zainteresowani!

To jest dla nas ćwiczenie edukacyjne - pracujemy krok po kroku!

tomtomnz
źródło
2
Powiedziałbym, że wyprzedzasz zespół programistów, który dyskutuje nad przeniesieniem pointcloudtypu danych do PostGIS i mam nadzieję, że wkrótce nad nim pracujesz . W międzyczasie istnieje narzędzie point2grid dla PDAL, które wygląda interesująco, ale go nie użyłem.
Mike T

Odpowiedzi:

2

Możesz wykonać tę automatyzację za pomocą SAGA , w systemie Windows lub Linux. SAGA może obsługiwać duże zbiory danych punktowych, importowane w formacie SAGA jako format pointcloud (spc). Następnie masz funkcje interpolacji lub triangulacji, a także eksport do Surfer. Automatyzacja odbywa się za pomocą narzędzia poleceń SAGA .

Niculita Mihai
źródło
Zaimportuj do CVS, możesz umieścić to w bashscript za pomocą polecenia psql COPY (WYBIERZ foo, słupek od czegokolwiek) DO '/tmp/ouput.csv' DELIMITER ',' CSV HEADER; Dzięki temu możesz usunąć QGIS z łańcucha narzędzi
simplexio
2

Jeśli masz oryginalne dane LIDAR w .las / laz, możesz użyć lastools, http://www.cs.unc.edu/~isenburg/lastools/ do obsługi dużych danych LiDAR. Jest las2tin i las2dem, które wyglądają, jakby działały bezpośrednio, ale ich nie użyłem. Użyłem poniżej, aby złożyć zestaw plików .laz w rastrowym DEM w 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 dobrze się automatyzuje i wygląda na to, że mogłyby przyciąć dane dla DEM lub TIN do wybranych granic.

Dave X
źródło