估算森林冠层密度和高度
许多环境应用都将森林冠层密度及高度用作变量。这些应用包括对生物量、森林范围和条件以及生物多样性的估算。冠层密度(又称为冠层盖度)是从空中观测到的植被地面比。冠层高度表示冠层顶端相对于地面的高度。激光雷达可用于确定这两个变量。
根据激光雷达点计算冠层密度和高度的步骤如下。首先,需要被归类为地面回波(裸露地表)与非地面回波的激光雷达信号。这种类型的点分类任务通常由数据提供程序执行。其次,需要考虑激光雷达数据的采集时间以及研究区域中的植被类型。如果存在许多落叶树并且是在秋季(落叶期)执行数据采集,则密度计算将不准确。
创建 LAS 数据集和图层
LAS 数据集可用来计算冠层密度和高度。使用创建 LAS 数据集地理处理工具来引用 LAS 格式的激光雷达数据。如果正在处理 ASCII 数据或要处理地理数据库中基于 LAS 的点,则请参阅本主题结尾的替代的工作流部分。创建 LAS 数据集后,使用图层来定义过滤器属性,这种属性将控制在 LAS 文件中使用哪些点。要创建图层,则添加数据集到地图或运行创建 LAS 数据集图层地理处理工具。定义图层时,选择所有回波,然后便可以处理类代码了。基本上,我们将处理两组点:地面点和地上点。下表包含了按 LAS 1.1 标准定义的 LAS 分类代码:
分类代码 | 分类类型 |
---|---|
0 | 从不分类 |
1 | 未分配 |
2 | 地面 |
3 | 低植被 |
4 | 中等植被 |
5 | 高植被 |
6 | 建筑物 |
7 | 噪点 |
8 | 模型关键 |
9 | 水域 |
10 | 为 ASPRS 定义预留 |
11 | 为 ASPRS 定义预留 |
12 | 重叠 |
13–31 | 为 ASPRS 定义预留 |
![注 注](rsrc/note.png)
如果点已分类,则在最近几年中创建的任何 LAS 文件均应使用这些代码。不幸的是,仍然存在一些不明确的点。例如,第 2 类是地面,但第 8 类也是地面。第 8 类(模型关键)点是一组特殊的地面点,用于等值线或其他需要一组已细化地面点的应用。是否具有此类点取决于数据的处理方式。如果不确定,可以同时指定这两类。即使发现不存在任何模型关键点,也不会有负面影响。植被也存在类似的问题。有时,供应商会将高于地面的所有事物都归为第 1 类,因为他们没有对这些事物进行更详细的分类。因此,如果不确定数据分类的详细信息,请使用第 1、3、4、5 类选择非地面点。这是一种获取全部植被点的合理方法。请注意,如果第 1 类中存在建筑物或其他人造非地面要素,则同样会获取到这些要素,并且它们会对结果的准确性产生细微影响。
计算密度
确定冠层密度最有效的方法是通过栅格化将研究区域划分为许多大小相等的小像元。在每个栅格像元中,将地上点的数量与点的总数相比。
此处要牢记的一项重要技巧是确定适合于此分析的像元大小。必须至少为平均点间距的四倍。可以使用更大的像元,但不能小于该值。
- 通过 RETURN_ COUNT 选项对地上点使用用作栅格数据的 LAS 点统计地理处理工具。
- 将生成的任意 NoData 像元转换为 0,这样后续操作便会将没有点的像元视为 0。可使用为空地理处理工具配合条件函数地理处理工具完成此操作。
- 对激光雷达地面多点重复执行步骤 1 和 2。
- 使用加地理处理工具将地上栅格和裸露地表栅格加在一起以获得每个像元的总数。
- 目前所创建的所有栅格都是长整型数据类型。您还需要一个浮点型栅格以便获取来自除地理处理工具(将在步骤 6 中使用)的浮点型输出。要生成浮点型栅格,请将加地理处理工具的输出栅格用作转为浮点型地理处理工具的输入。
- 现在使用除地理处理工具求出地上栅格数和浮点型栅格总数之比。您将得到一个介于 0.0 和 1.0 之间的比值,其中 0.0 表示无冠层,1.0 表示冠层浓密。
![“用作栅格数据的 LAS 点”地理处理工具 “用作栅格数据的 LAS 点”地理处理工具](015w/GUID-BBEB7855-0A09-41FE-B231-60D8EE7CE456-web.png)
![为空地理处理工具 为空地理处理工具](015w/GUID-0A42B608-7482-4E78-8058-91481DBC0CBA-web.png)
![条件函数地理处理工具 条件函数地理处理工具](015w/GUID-4DFFFFC6-8825-4D70-AA42-39E588BB5BAB-web.png)
![加地理处理工具 加地理处理工具](015w/GUID-F7643171-58EF-4543-A6FB-EC792C3E9667-web.png)
![“转为浮点型”地理处理工具 “转为浮点型”地理处理工具](015w/GUID-57354CFD-0B03-452D-A571-F0F6B9E51DEF-web.png)
![必填 必填](015w/GUID-77DBA5CC-93C7-4323-84BE-D60549B7D1C2-web.png)
下图表示了冠层密度。浅色区域表示基本没有植被。在这些区域中,大部分激光雷达波能够“看到”地面。深绿色区域表示植被冠层很茂密,激光雷达无法穿透到达地面。
![密度分析结果 密度分析结果](015w/GUID-994A185E-CF90-4ED7-9845-38FBE37E53B1-web.png)
计算高度
要确定冠层高度,需要用第一个回波表面 (DSM) 减去裸露地表 (DEM)。请按照从大型激光雷达点集合创建栅格 DEM 和 DSM 主题中的步骤操作生成这两个表面。
- 生成第一个回波栅格和裸露地表栅格后,使用减地理处理工具确定这两个栅格数据集的差值。对森林而言,差值结果即为冠层高度。
![“减”工具 “减”工具](015w/GUID-0C8D6728-A2E1-42BD-8449-6FE367C275AD-web.png)
下图表示了距地面的高度。它的范围为从蓝色(基本没有高度)到橙色(最高)。
![高度结果 高度结果](015w/GUID-68D7F290-C5C9-4E29-B34A-4B78E7EFE987-web.png)
激光雷达可用于计算植被的密度和高度。这具有多种用途,其中包括生物量估算、碳估算以及森林管理。
替代的工作流
如果激光雷达点是 ASCII 格式,则需要拥有两组文件:仅包含地面点的文件,以及仅包含非地面点的文件。如果文件以这种方式进行组织,则应使用 3D ASCII 文件转要素类地理处理工具将点加载到地理数据库的多点要素类中。将此工具运行两次,每组文件各运行一次。然后运行点转栅格地理处理工具来替换上述基于 LAS 数据集的工作流中使用的“LAS 点统计”工具。
如果您具有 LAS 格式的数据,但更希望仅将感兴趣点加载到地理数据库中,并处理这些点,则需要使用 LAS 转多点地理处理工具两次。第一次是为了加载地面点,第二次是为了加载非地面点。然后运行点转栅格地理处理工具来替换上述基于 LAS 数据集的工作流中使用的“LAS 点统计”工具。