由激光雷达点描绘数据区

发送测量到的激光雷达或摄影测量数据时,通常不带详细的数据区边界。在通常情况下,测量区的 x 和 y 范围都由覆盖感兴趣区域的切片系统定义,而且这些切片都由数据来填充。下图展示了一个投影的激光雷达数据切片。这些切片的范围大致近似于实际研究区域的边界。

研究区域切片系统

通常会出现这样的情况,即实际的激光雷达数据未能完全覆盖投影区域周界上各切片的范围。如下图所示,数据仅能覆盖极小的范围,而且没有明确或绝对的边界,只能进行推断。此图形的中心区域便是一个这样的切片。

激光雷达点数据的覆盖范围

无论采取哪种方式,覆盖区域通常都不会是完全填充的矩形。

问题

如果在没有预先声明数据区(也就是说,在定义 terrain 数据集或 TIN 时包括剪切多边形)的情况下创建表面,则周界周围一些实际上是空白的区域也将被视为数据区。这些区域内的分析结果是不可靠的,因为估计的高度可能是基于很远处的样本得出的。

左下侧的图展示了密集激光雷达点集合(绿色区域)。内部的空隙是水体(在这些区域中激光雷达点通常被忽略)。很容易看出数据边界的形状不规则,但是除非以剪切多边形的形式提供明确的范围,否则 TIN 和 terrain 数据集的相关工具会填充其中的空白,这样将使实际的数据范围变得过分简单化。

激光雷达点没有剪切多边形的表面

如您所知,数据集范围之外的区域应排除在表面之外。问题就是要提供一个可精确表示此范围的多边形。

解决办法

解决办法是使用这些可以在表面中强制加入适当插值区的点来合成数据边界。下方,左侧的图片显示了激光雷达点。中间图片所显示的多边形是由点合成的边界。右侧的图片是使用激光雷达点和剪切多边形构成的表面。

激光雷达点剪切多边形经边界强化的表面

点间距是确定数据区时使用的主要变量。测量通常以明确的最小点间距进行,以便为插值器提供控制。不符合密度要求的区域是特例。这些特例通常可分为以下几个类别:水体、被遮盖的区域和休假(将遮盖区域的数据发送回数据提供者以进行修整)。绝大多数数据都应符合样本密度规范。点间距通常显示在元数据中。如果不知道激光雷达数据的点间距,则可使用 3D 分析的点文件信息地理处理工具。或者,使用 ArcMap 中的测量工具估算点间距。要了解点间距的详细信息,请参阅平均点间距

由激光雷达点描绘数据区

在了解激光雷达数据的点间距后,可按以下步骤描绘数据区:

步骤:
  1. 使用点转栅格地理处理工具栅格化激光雷达点。

    激光雷达点的栅格化有助于聚合激光雷达点所覆盖的区域。这便为后续步骤的处理奠定了的良好数据结构基础。操作过程只需告知地理处理工具要使用的单元分配类型和输出单元大小即可。使用 COUNT 作为单元分配类型值。将单元大小指定为比激光雷达数据的平均点间距大数倍的值。否则,因为点的间距不等将得到很多噪点。从提高处理效率和减少噪点的角度来看,使用的单元大小越大,效果便越好;但是使用的单元太大,结果拟合的紧密度又会受到影响。因此,该值的设定要适当。最好是从四倍平均点间距开始进行调整。

    点转栅格
  2. 使用条件函数地理处理工具为所有数据单元分配同一个值。

    在此工作流程中使用条件函数地理处理工具只是将所有栅格的数据单元转换为具有同一个值的单元。此值可定义将在第三步时扩展的栅格区域。所需要的仅是从点转栅格工具中取得输出并为正表达提供一个常数值。所有非零值单元都将被视为正值并将为其指定常数值。由于 COUNT 在栅格化过程中被用作单元分配方法,所以任何其中具有点的单元都必须具有大于零的值。

    “条件函数”工具
  3. 使用扩展地理处理工具填充小型 NoData 区域。

    除非使用的单元大小比平均点间距粗略很多,否则很有可能仍然有 NoData 单元。多数单元都可使用扩展地理处理工具加以消除。用户可能会希望将这些单元尽可能去除,以使后续步骤的矢量化过程中生成的多边形不会具有过多孔洞。但这样会导致不必要的高昂代价。

    “扩展”工具

    扩展工具可使感兴趣的区域向外扩展。在这种情况下,该区域的范围将包括所有编码值为 1 的数据单元。这将有效地减少内部的小缝隙。

    白色为 No Data 单元扩展的结果

    左图显示了许多单独的 NoData 单元和一些 NoData 单元的小型团簇(白色区域)。右侧图像则显示了使用扩展地理处理工具后的结果,其中的 NoData 单元(白色)大多都已移除。输出中仍然保留有一些孤立的 NoData 区域是可以接受的。剩下 NoData 单元将在最后的步骤加以处理。

  4. 使用收缩地理处理工具缩小数据单元的整体范围。

    在使用扩展消除孤立的 NoData 单元时,同时也会使数据区向外扩展。因此在实际上需要略微缩小。剪切多边形需要比实际的点范围小,这样在 terrain 数据集或 TIN 试图沿多边形边界估计 z 值时,在两侧均可找到点。这对获取估算精确的 z 值十分必要。要缩小栅格的数据边界,可使用下面所示的收缩地理处理工具。

    略微缩小该范围,以使步骤 5 中产生的多边形小于点的实际数据范围。这使软件在多边形边界上可以估算出更精确的 z 值。

    “收缩”工具

    此时,栅格相对更加整齐,其中的数据单元范围也略微小于激光雷达点的范围。

  5. 使用栅格转面地理处理工具使栅格矢量化。

    栅格转面地理处理工具可将栅格转换为面要素类。确保简化面选项已选中。如果并未选中,则输出将呈不平滑的锯齿状,而且将包含过多不必要的折点。

    栅格转面工具将输出面要素类。该结果代表过程开始时使用的点的数据范围。

    栅格转面工具

    此时,这一过程便基本完成。需要检查输出的正确性。某些情况下还需要一个额外的步骤,即移除剪切多边形中仍然存在的孔洞。

  6. 使用消除面部分地理处理工具移除所有仍存在的小孔洞。

    消除面部分地理处理工具可移除所有的内部环,仅留下外部的边界。

此时便得到了剪切多边形,可使用它定义 terrain 数据集或 TIN。剪切多边形应基本与激光雷达点的数据范围一致,但略微小于该范围。左图显示了得到的剪切多边形。右图是放大的视图,其中显示了面范围与源点范围的大小关系。注意剪切多边形范围略微小于源点边界的具体程度。

生成的剪切多边形剪切多边形和激光雷达点

7/10/2012