PostGIS插值/三角剖分选项


12

我正在寻找一个开源的工作流程来自动化DEM的构建。我们有一系列属于LIDAR数据集的站点。我们想要为每个站点创建站点特定的DEM,并且正在寻求使过程自动化。

到目前为止,我们已经实现了自动化:

  • 数据加载到PostGIS中(激光雷达和站点点)
  • 创建站点边界(st_buffer和st_envelope的组合)

现在,我们正在寻找在每个站点上插值数据子集并将其导出到Surfer7网格文件的选项。

当前,我们正在使用QuantumGIS查询PostGIS数据库中的数据,导出到csv文件,然后在Surfer7中手动导入和网格化数据。希望我们也可以使它自动化。

因此,对于此用例,我们希望将激光雷达数据三角剖分为DEM。我们可能想到的其他潜在情况涉及不同的插值方法-因此,如果有用于反距离和克里金法的选项-我们很感兴趣!

对我们来说,这是一个非常学习的练习-我们正在逐步进行中!


2
我要说的是,您在开发团队的领导之下,他们一直在讨论将pointcloud数据类型引入PostGIS,并希望能很快进行开发。同时,PDAL有一个points2grid实用程序看起来很有趣,但我没有使用过。
Mike T

Answers:


2

您可以在Windows或Linux中使用SAGA进行自动化。SAGA可以处理以SAGA格式作为点云格式(spc)导入的大点数据集。然后,您将具有用于插值或三角剖分的功能,并且还可以导出到Surfer。自动化是使用SAGA命令工具执行的。


导入到CVS中,您可以使用psql命令COPY(SELECT foo,bar FROM what)将其放入bashscript中。TO'/tmp/ouput.csv'DELIMITER',CSV HEADER; 这样您就可以从工具链中删除QGIS
simplexio 2013年

2

如果您在.las / laz中拥有原始LIDAR数据,则可以使用lastools(http://www.cs.unc.edu/~isenburg/lastools/)处理大型LiDAR数据。看起来像las2tin和las2dem可以直接工作,但是我没有使用它们。我已使用以下内容将一组.laz文件组装到GRASS中的栅格DEM中:

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可以很好地实现自动化,并且看起来可以将DEM或TIN的数据裁剪到您选择的范围。

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.