基于点云进行插值 (数据管理)

摘要

从点云对数字地形模型 (DTM) 或数字表面模型 (DSM) 进行插值。

使用情况

  • 点云的格式可以是 LAS 文件或解决方案点表。

参数

标注说明数据类型
输入 LAS 文件夹或点表

输入文件、文件夹或要素图层的路径和名称。 输入可以是 LAS 文件的文件夹或正射映射工具的解决方案点表。 对于云存储,请提供云存储路径,例如 C:\Temp\Cloud.acs\lasfolder

LAS 文件可以是生成点云工具的输出,其中 LAS 点分为地面点和地上点两种类型。 解决方案点表为计算区域网平差工具或计算相机模型工具的输出。

Folder; File; Feature Class; Feature Layer
输出栅格

输出栅格数据集的位置、名称和文件扩展名。 您也可以提供云存储路径(例如 C:\Temp\Cloud.acs\lasfolder)以保存输出栅格数据集。

可以以大多数可写栅格格式(例如 TIFF、CRF 或 IMG)创建输出。

Raster Dataset
像元大小

输出栅格数据集的像元大小。

Double
插值方法
(可选)

指定将用于从点云对输出栅格数据集进行插值的方法。

  • TIN 线性插值法将使用三角测量方法。 亦称作不规则三角网 (TIN) 线性插值法,适用于不规则分布的稀疏点,例如从区域网平差计算中得到的解决方案点。
  • TIN 自然邻域插值法将使用自然邻域法。 其与三角测量类似,但是生成的表面更为平滑且运算量更大。
  • 反距离权重平均插值法将使用反距离权重 (IDW) 平均方法。 此方法适用于规则分布的密集点,例如通过生成点云工具得到的点云 LAS 文件。 IDW 搜索半径将根据平均点密度自动进行计算。
String
平滑方法
(可选)

指定将用于平滑输出栅格数据集的滤波器。

  • 高斯 3 x 3将使用具有 3 x 3 窗口的高斯滤波器。
  • 高斯 5 x 5将使用具有 5 x 5 窗口的高斯滤波器。
  • 高斯 7 x 7将使用具有 7 x 7 窗口的高斯滤波器。
  • 高斯 9 x 9将使用具有 9 x 9 窗口的高斯滤波器。
  • 无平滑将不使用任何平滑滤波器。
String
表面类型
(可选)

指定将创建数字地形模型还是数字表面模型。

  • 数字地形模型将通过仅插值地面点来创建数字地形模型。
  • 数字表面模型将通过插值所有点来创建数字表面模型。
String
输入填充 DEM
(可选)

用于填充 NoData 区域的 DEM 栅格输入。 NoData 区域可以存在于像素未从输入获得足够生成值所需信息的位置。

Raster Dataset; Raster Layer; Mosaic Dataset; Mosaic Layer
分类地面选项
(可选)

从输入 LAS 数据中提取地面点。

表面类型参数设置为数字地形模型时,该选项将启用。

注:

地理处理窗格中设置选项时,输入用于检测地面点的方法名称,并在列表框中选择相应的值。

  • Classify - 根据不同地形类型,采用多种方法分类地面点。 选项不区分大小写。
    • standard - 该方法对坡度变化具有容差,因此可以捕获到使用 conservative 方法无法捕获到的地面地形上的平缓波动,但是该方法无法捕获到通过 aggressive 方法捕获到的坡度起伏较大的地貌波动。 这是默认设置。
    • conservative - 与其他方法相比,该方法对地面坡度的变化使用了更为严格的限制,因此该方法可将地面与草地和灌木丛等低地植被区分开来。 该方法最适合最小曲率的地形。
    • aggressive - 对于 standard 方法容易忽略的山脊和山顶等地貌起伏较大的地面区域,此方法可检测到这些内容。 该方法最好在 ReuseGround 选项设置为 1 的情况下用于该工具的第二个迭代。 该方法不宜用于城市区域或地势平坦的乡村区域,因为可能会将地势较高的对象(例如,发电塔、植被和建筑物局部)错误分类为地面。
  • LowNoise - 用于将点分类为低噪声点的地面以下距离。 单位为米。 默认值为 0.25 米。
  • HighNoise - 用于将点分类为高噪声点的地面以上距离。 单位为米。 默认值是 100 米。
  • ReuseGround - 指定现有地面点是否会重新分类或重新使用。 0 表示重新分类,1 表示重用。 默认值为 0。
  • ReuseLowNoise - 指定是否重用或重新分类现有低噪声点。 0 表示重新分类,1 表示重用。 默认值为 0。
  • ReuseHighNoise - 指定是否重用或重新分类现有高噪声点。 0 表示重新分类,1 表示重用。 默认值为 0。
Value Table

arcpy.management.InterpolateFromPointCloud(in_container, out_raster, cell_size, {interpolation_method}, {smooth_method}, {surface_type}, {fill_dem}, {options})
名称说明数据类型
in_container

输入文件、文件夹或要素图层的路径和名称。 输入可以是 LAS 文件的文件夹或正射映射工具的解决方案点表。 对于云存储,请提供云存储路径,例如 C:\Temp\Cloud.acs\lasfolder

LAS 文件可以是生成点云工具的输出,其中 LAS 点分为地面点和地上点两种类型。 解决方案点表为计算区域网平差工具或计算相机模型工具的输出。

Folder; File; Feature Class; Feature Layer
out_raster

输出栅格数据集的位置、名称和文件扩展名。 您也可以提供云存储路径(例如 C:\Temp\Cloud.acs\lasfolder)以保存输出栅格数据集。

可以以大多数可写栅格格式(例如 TIFF、CRF 或 IMG)创建输出。

Raster Dataset
cell_size

输出栅格数据集的像元大小。

Double
interpolation_method
(可选)

指定将用于从点云对输出栅格数据集进行插值的方法。

  • TRIANGULATION将使用三角测量方法。 亦称作不规则三角网 (TIN) 线性插值法,适用于不规则分布的稀疏点,例如从区域网平差计算中得到的解决方案点。
  • NATURAL_NEIGHBOR将使用自然邻域法。 其与三角测量类似,但是生成的表面更为平滑且运算量更大。
  • IDW将使用反距离权重 (IDW) 平均方法。 此方法适用于规则分布的密集点,例如通过生成点云工具得到的点云 LAS 文件。 IDW 搜索半径将根据平均点密度自动进行计算。
String
smooth_method
(可选)

指定将用于平滑输出栅格数据集的滤波器。

  • GAUSS3x3将使用具有 3 x 3 窗口的高斯滤波器。
  • GAUSS5x5将使用具有 5 x 5 窗口的高斯滤波器。
  • GAUSS7x7将使用具有 7 x 7 窗口的高斯滤波器。
  • GAUSS9x9将使用具有 9 x 9 窗口的高斯滤波器。
  • NONE将不使用任何平滑滤波器。
String
surface_type
(可选)

指定将创建数字地形模型还是数字表面模型。

  • DTM将通过仅插值地面点来创建数字地形模型。
  • DSM将通过插值所有点来创建数字表面模型。
String
fill_dem
(可选)

用于填充 NoData 区域的 DEM 栅格输入。 NoData 区域可以存在于像素未从输入获得足够生成值所需信息的位置。

Raster Dataset; Raster Layer; Mosaic Dataset; Mosaic Layer
options
[[name, value],...]
(可选)

从输入 LAS 数据中提取地面点。

surface_type 参数设置为 DTM 时,此参数处于活动状态。

  • Classify - 根据不同地形类型,采用多种方法分类地面点。 选项不区分大小写。
    • standard - 该方法对坡度变化具有容差,因此可以捕获到使用 conservative 方法无法捕获到的地面地形上的平缓波动,但是该方法无法捕获到通过 aggressive 方法捕获到的坡度起伏较大的地貌波动。 这是默认设置。
    • conservative - 与其他方法相比,该方法对地面坡度的变化使用了更为严格的限制,因此该方法可将地面与草地和灌木丛等低地植被区分开来。 该方法最适合最小曲率的地形。
    • aggressive - 对于 standard 方法容易忽略的山脊和山顶等地貌起伏较大的地面区域,此方法可检测到这些内容。 该方法最好在 ReuseGround 选项设置为 1 的情况下用于该工具的第二个迭代。 该方法不宜用于城市区域或地势平坦的乡村区域,因为可能会将地势较高的对象(例如,发电塔、植被和建筑物局部)错误分类为地面。
  • LowNoise - 用于将点分类为低噪声点的地面以下距离。 单位为米。 默认值为 0.25 米。
  • HighNoise - 用于将点分类为高噪声点的地面以上距离。 单位为米。 默认值是 100 米。
  • ReuseGround - 指定现有地面点是否会重新分类或重新使用。 0 表示重新分类,1 表示重用。 默认值为 0。
  • ReuseLowNoise - 指定是否重用或重新分类现有低噪声点。 0 表示重新分类,1 表示重用。 默认值为 0。
  • ReuseHighNoise - 指定是否重用或重新分类现有高噪声点。 0 表示重新分类,1 表示重用。 默认值为 0。
Value Table

代码示例

InterpolateFromPointCloud 示例 1(Python 窗口)

这是一个 InterpolateFromPointCloud 函数的 Python 示例,其中通过重新分类地面点生成了 DTM。

# Import system modules 
import arcpy
 
# Execute 
arcpy.management.InterpolateFromPointCloud(in_container=r"C:\data\LASFoler", out_raster=r"C:\data\dtm.crf", cell_size=0.2, interpolation_method="IDW", smooth_method="GAUSS5x5", surface_type="DTM", fill_dem=None, options="Classify standard;LowNoise 0.25;HighNoise 110;ReuseGround 0;ReuseLowNoise 1;ReuseHighNoise 1")
InterpolateFromPointCloud 示例 2(独立脚本)

这是一个 InterpolateFromPointCloud 函数 Python 示例,其中生成了 DSM。

# Define input parameters 
import arcpy
in_container="C:/data/LASFoler"
out_raster="C:/data/dsm.crf"
# Execute 
arcpy.management.InterpolateFromPointCloud(in_container,out_raster, 0.2, "TRIANGULATION", "GAUSS5x5", "DSM")
InterpolateFromPointCloud 示例 3(独立脚本)

这是一个 InterpolateFromPointCloud 函数 Python 示例,其中通过云存储生成了 DSM。

# Define input parameters 
import arcpy
in_container="C:/data/LASFoler"
out_raster="C:/data/Azure.acs/ProductFolder/dsm.crf"
# Execute 
arcpy.management.InterpolateFromPointCloud(in_container,out_raster, 0.2, "TRIANGULATION", "GAUSS5x5", "DSM")

许可信息

  • Basic: 否
  • Standard: 需要 ArcGIS Reality for ArcGIS Pro
  • Advanced: 是

相关主题