我经常遇到带有高程数据的LiDAR(.LAS)文件。当我这样做时,我总是很难找到一种将它们转换为DEM进行水文分析的方法。我已经开发了几个工作流,即使它们通常(但并非总是如此!)输出我需要的东西,但它们似乎并没有非常优化。另外,我希望有一个可以依靠的工作流,而不必每次都尝试三种不同的工作。
这些是我当前的工作流: PointCloudCombiner-> RasterDEMGenerator-> Geotiff编写器。 (参考时间5分钟)。
ArcMap:


LAS到多点->多点到地形->地形到栅格(10分钟)。
LAS到多点->创建TIN(因为Topo到栅格无法读取Point.Z值)-> TIN到栅格(15分钟)。
创建新的Mosaic数据集并将LAS添加为栅格。导出到ESRI GRID。

LAStools:
LAStoTXT->制作XY事件层->到Shapefile->地形到光栅(3-4小时)知道要栅格化的LAS数据集和要TIN的LAS数据集,但是他们通常无法识别我的LAS文件。据我了解,Arcmap中的#1是ESRI首选方法(?)。您将使用什么?
我可以访问ArcGIS Standard 10.1(即将于10.2启用),3D分析师,空间分析师。欢迎使用开放源代码解决方案,并在需要时使用脚本编写(Python)。

评论

另请参阅:gis.stackexchange.com/questions/13618/…

您是否尝试过LAStools LASGRID?使用LAStoTXT并制作事件层听起来并不正确。由于LAStoTXT会生成ASCII txt文件-您应使用ASCII转RASTER工具箱工具来创建栅格。或者,您可以尝试将LASGRID从LAS直接转到栅格。这将为您提供LiDAR DEM,您可以从中生成轮廓并将其用作TopoToRaster工具的输入。至于水文分析,您无需详细说明,但我不能确信仅LiDAR是合适的。您需要测量流并将其刻录到栅格中,进行归一化,减少噪声等。

我应该检查一下该工具。我确实有溪流,湖泊,水槽等,但这是供以后处理。在这个阶段,我只是在寻找一种方法,可以将LAS文件转换为DEM,并且不会花费太多时间并且可以生成有用的高程模型。

对于LAStools,您应该使用las2dem或blast2dem。

#1 楼

将LiDAR数据集转换为DEM时,您要获取一组离散数据点并将其转换为单个连续数据集。假设您的.las文件包含X(纬度),Y(经度)和Z(海拔)值,平均分辨率约为1米。此处的分辨率非常重要-我们仅讨论平均值,因此几乎找不到整个数据集的〜1米分辨率。取而代之的是,我们将找到该分辨率的“估计”值之内的值。因此,您可以将这些点转换为栅格DEM或TIN。 X和Y值应该看到最小的失真,但是您会注意到Z值可能不是您所期望的。这是因为计算机不知道未落在LiDAR点之一上的单元格中正确的Z值是多少。在LiDAR点之间,已应用插值算法来估计合理的Z值。选择相对于分析目标正确的插值方法是从LiDAR到DEM的关键部分。在输出DEM上设置适当的分辨率很重要-始终将分辨率设置为低于LiDAR数据集的分辨率。因此,对于约1米的分辨率,我将为DEM设置3米的分辨率,以尽量减少失真。

我有使用来自LiDAR的DEM研究滑坡和泥石流的经验。滑坡和泥石流是非常线性的特征,发生在地形中的其他线性特征附近。因此,当我从LiDAR转换为DEM时,我想要一种最能强调线性特征的插值方法。这恰好是TIN(不规则三角网)。您说的目的是进行水文分析。也许您应该尝试样条插值方法来构建您的DEM。样条插值的工作原理是在所有数据点上绘制连续的重叠线,以创建非常平滑的栅格表面。确定您的水槽,填满它们,绘制计数,然后重复。错误的问题。您应该询问使用哪种插值方法,而不是要求使用软件工作流程来构建水文正确的DEM。如果您是我,我会尝试样条曲线插值方法。

在软件方面,处理LiDAR数据需要占用大量CPU / RAM。如果您有> 6 GB的RAM,我建议您使用GRASS GIS。他们拥有我曾经使用过的最好的LiDAR处理软件(它是FOSS),但是您必须做一些内存分配。否则,我建议您坚持使用ArcGIS。他们在他们的网站上有一些很棒的文档,说明如何做您想做的事情。

评论


我应该尝试另一种插值方法。我的计算机上有40 GB的RAM,但是我在内存分配方面没有经验。但是,ArcMap在处理大型las数据集时有点慢(我通常有50至1亿个未过滤的点),因此仍然值得一试。

–马丁
2013年12月6日下午6:38

如果您的PC上有40 GB RAM,那么您甚至不必担心使用GRASS时会修改内存分配。这是处理GRASS中LiDAR点的重要指南:grasswiki.osgeo.org/wiki/LIDAR。我喜欢这个软件的地方在于,它可以让您完全控制数据处理的每个步骤,这与其他专有软件不同,它可以在幕后为您做出许多通用决策。

–asonnenschein
2013年12月6日15:27

极好的答案。提出传感器分辨率(激光雷达点间距)和表面分辨率之间的区别是很大的。但是,出于水文目的,ToopoToRaster插值远远优于样条线方法。如果您没有或希望使用Esri工具,请访问源代码并获取ANUDEM(gis.stackexchange.com/a/131870/108)。

–马特·威尔基
17年11月2日在20:47

如果您使用TopoToRaster以外的方法(该方法是为轮廓而不是激光雷达数据设计的,是否不确定是否合适?),则应在DEM上进行一些后处理坑填充/裂化,以确保水文正确性。您可以使用Whitebox GAT(出色的工具,但无法处理庞大的数据集),Taudem,ArcHydro对ArcGIS,SAGA的免费扩展或Grass r.hydrodem进行扩展。所有这些都是开源的。对于大型数据集,我建议使用Grass或Taudem。

–fionag11
19年7月17日在17:24



#2 楼

如果只需要执行一次,则可能要考虑下载Merrick&Company的MARS 30天评估版。完整的软件套件非常昂贵(11995美元),但是我认为使用评估软件,您可以使用已经存在的水体数据集来强制围绕多边形的恒定高程。

评论


那确实是一笔不菲的价格!我不时地这样做,所以恐怕试验不会成功。但是,无论如何我都应该下载它,看看它能做什么!

–马丁
2013年12月6日下午6:33