Opções de interpolação / triangulação PostGIS

12

Estou procurando um fluxo de trabalho de código aberto para automatizar a construção do DEM. Temos uma série de sites que se enquadram em um conjunto de dados LIDAR. Queremos criar DEMs específicos para cada site e pretendemos automatizar o processo.

Até agora, automatizamos:

  • carregamento de dados no PostGIS (LIDAR e pontos do site)
  • criação de limites do site (combinação de st_buffer e st_envelope)

Agora, estamos procurando opções para interpolar os subconjuntos de dados em cada site e exportá-los para os arquivos de grade do Surfer7.

Atualmente, estamos consultando os dados no banco de dados PostGIS usando o QuantumGIS, exportando para arquivos csv e depois importando e agrupando manualmente os dados no Surfer7. Espero que possamos automatizar isso também.

Portanto, neste caso de uso, gostaríamos de triangular nossos dados do lidar para um DEM. Outros casos potenciais em que podemos pensar envolvem diferentes métodos de interpolação - por isso, se existem opções para distância inversa e krigagem - estamos interessados!

Este é um exercício de aprendizado para nós - estamos trabalhando em etapas de bebê!

tomtomnz
fonte
2
Eu diria que você está à frente da equipe de desenvolvimento, que vem discutindo a possibilidade de levar o pointcloudtipo de dados para o PostGIS, e esperamos trabalhar em breve . Nesse meio tempo, há uma points2grid utilidade para PDAL que parece interessante, mas eu não tê-lo usado.
Mike T

Respostas:

2

Você poderia fazer essa automação com SAGA , no Windows ou no Linux. O SAGA pode lidar com grandes conjuntos de dados pontuais, importados no formato SAGA como formato pointcloud (spc). Então você tem funções para interpolação ou triangulação e também exporta para o Surfista. A automação é realizada usando a ferramenta de comando SAGA .

Niculita Mihai
fonte
Importar para o CVS, você pode colocar isso em bashscript com o comando psql COPY (SELECT foo, bar FROM qualquer) TO '/tmp/ouput.csv' DELIMITER ',' CSV HEADER; Com isso você pode deixar cair QGIS da cadeia de ferramentas
simplexio
2

Se você possui os dados originais do LIDAR em .las / laz, pode usar lastools, http://www.cs.unc.edu/~isenburg/lastools/ para lidar com grandes dados do LiDAR. Existe um las2tin e um las2dem que parecem funcionar diretamente, mas eu não os usei. Eu usei o abaixo para montar um conjunto de arquivos .laz em um DEM raster no 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

As ferramentas de barra automatizariam bem e pareceriam que poderiam cortar os dados de um DEM ou TIN nos limites de sua escolha.

Dave X
fonte