Opciones de interpolación / triangulación PostGIS

12

Estoy buscando un flujo de trabajo de código abierto para automatizar la construcción de DEM. Tenemos una serie de sitios que se encuentran dentro de un conjunto de datos LIDAR. Queremos crear DEM específicos del sitio para cada sitio, y estamos buscando automatizar el proceso.

Hasta ahora, hemos automatizado:

  • carga de datos en PostGIS (LIDAR y puntos de sitio)
  • creación de límites del sitio (combinación de st_buffer y st_envelope)

Ahora estamos buscando opciones para interpolar los subconjuntos de datos en cada sitio y exportarlos a archivos de cuadrícula Surfer7.

Actualmente, estamos consultando los datos en la base de datos PostGIS usando QuantumGIS, exportando a archivos csv, luego importando y cuadriculando manualmente los datos en Surfer7. Esperemos que podamos automatizar esto también.

Entonces, para este caso de uso, nos gustaría triangular nuestros datos LIDAR a un DEM. Otros posibles casos en los que podemos pensar involucran diferentes métodos de interpolación, así que si hay opciones para distancia inversa y kriging, ¡estamos interesados!

Este es un ejercicio de aprendizaje para nosotros, ¡estamos trabajando en pasos pequeños!

postgis dem interpolation lidar tomtomnz
fuente
2
Diría que está por delante del equipo de desarrollo, que ha estado discutiendo la pointcloudposibilidad de llevar el tipo de datos a PostGIS, y con suerte trabajará en él pronto . Mientras tanto, hay una utilidad points2grid para PDAL que parece interesante, pero no la he usado.
Mike T

Respuestas:

2

Podrías hacer esta automatización con SAGA , en Windows o en Linux. SAGA puede manejar conjuntos de datos de puntos grandes, importados en formato SAGA como formato pointcloud (spc). Luego tiene funciones de interpolación o triangulación, y también exporta a Surfer. La automatización se realiza utilizando la herramienta de comando SAGA .

Niculita Mihai
fuente
Importar a CVS, puede poner esto en bashscript con el comando psql COPIAR (SELECCIONAR foo, barra DESDE lo que sea) A '/tmp/ouput.csv' DELIMITER ',' CSV HEADER; Con eso, puede eliminar QGIS de la cadena de herramientas
simplexio
2

Si tiene los datos LIDAR originales en .las / laz, puede usar lastools, http://www.cs.unc.edu/~isenburg/lastools/ para manejar grandes datos LiDAR. Hay un las2tin y las2dem que parecen funcionar directamente, pero no los he usado. He usado lo siguiente para ensamblar un conjunto de archivos .laz en un DEM ráster dentro de 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 se automatizaría bien y parecería que podrían recortar los datos para un DEM o TIN a los límites de su elección.

Dave X
fuente