Opciones de interpolación / triangulación de PostGIS

Question

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 MDEs específicos para cada sitio, y estamos buscando automatizar el proceso.

Hasta ahora, hemos automatizado:

• carga de datos en PostGIS (LIDAR y puntos del 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 de Surfer7.

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

Así que para este caso de uso, nos gustaría triangular nuestros datos lidar a un DEM. Otros casos potenciales en los que podemos pensar implican diferentes métodos de interpolación - así que si hay opciones para la distancia inversa y el kriging - ¡estamos interesados!

Esto es un ejercicio de aprendizaje para nosotros, ¡estamos trabajando con pasos de bebé!

Answer 1

Puedes hacer esta automatización con SAGA en Windows o en Linux. SAGA puede manejar grandes conjuntos de datos de puntos, importados en formato SAGA como formato de nube de puntos (spc). Luego tiene funciones para la interpolación o la triangulación, y también la exportación a Surfer. La automatización se realiza mediante Herramienta de mando SAGA .

Answer 2

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