Options d'interpolation / triangulation PostGIS

12

Je recherche un flux de travail open source pour automatiser la construction de DEM. Nous avons une série de sites qui font partie d'un ensemble de données LIDAR. Nous voulons créer des DEM spécifiques au site pour chaque site, et nous cherchons à automatiser le processus.

Jusqu'à présent, nous avons automatisé:

  • chargement des données dans PostGIS (LIDAR et points de site)
  • création de limites de site (combinaison de st_buffer et st_envelope)

Nous recherchons maintenant des options pour interpoler les sous-ensembles de données sur chaque site et les exporter vers des fichiers de grille Surfer7.

Actuellement, nous interrogeons les données dans la base de données PostGIS à l'aide de QuantumGIS, exportons vers des fichiers csv, puis importons et quadrillons manuellement les données dans Surfer7. J'espère que nous pourrons également automatiser cela.

Donc, pour ce cas d'utilisation, nous aimerions trianguler nos données lidar en DEM. D'autres cas potentiels auxquels nous pouvons penser impliquent différentes méthodes d'interpolation - donc s'il existe des options pour la distance inverse et le krigeage - nous sommes intéressés!

C'est vraiment un exercice d'apprentissage pour nous - nous travaillons par petits pas!

tomtomnz
la source
2
Je dirais que vous êtes en avance sur l'équipe de développement, qui a discuté de la pointcloudpossibilité d' apporter le type de données à PostGIS, et nous espérons y travailler bientôt . En attendant, il existe un utilitaire points2grid pour PDAL qui semble intéressant, mais je ne l'ai pas utilisé.
Mike T

Réponses:

2

Vous pouvez faire cette automatisation avec SAGA , sous Windows ou sous Linux. SAGA peut gérer des ensembles de données de gros points, importés au format SAGA au format pointcloud (spc). Ensuite, vous avez des fonctions d'interpolation ou de triangulation, et exportez également vers Surfer. L'automatisation est effectuée à l'aide de l' outil de commande SAGA .

Niculita Mihai
la source
Importez dans CVS, vous pouvez le mettre en bashscript avec la commande psql COPY (SELECT foo, bar FROM any) TO '/tmp/ouput.csv' DELIMITER ',' CSV HEADER; Avec cela, vous pouvez supprimer QGIS de la chaîne d'outils
simplexio
2

Si vous avez les données LIDAR originales en .las / laz, vous pouvez utiliser lastools, http://www.cs.unc.edu/~isenburg/lastools/ pour traiter les gros volumes de données LiDAR. Il y a un las2tin et un las2dem qui semblent fonctionner directement, mais je ne les ai pas utilisés. J'ai utilisé ce qui suit pour assembler un ensemble de fichiers .laz dans un DEM raster dans 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 s'automatiserait bien et donnerait l'impression de pouvoir couper les données d'un DEM ou d'un TIN aux limites de votre choix.

Dave X
la source