插值 Shape (3D Analyst)

需要 3D Analyst 许可。

需要 Spatial Analyst 许可。

摘要

通过从表面插入 Z 值创建 3D 要素。

了解有关“插值 Shape”工作原理的详细信息

插图

插值 Shape 工具图示

使用情况

  • 此工具使用从输入表面的重叠部分派生的高度值创建 3D 要素。 3D 面将仅在其周边存储 z 值,因为 3D 面的内部在渲染时将被随机定义。 因此,3D 面通常不适合表示非平面高度信息。 要生成表面的真实表示,请使用面插值为多面体工具。

  • 任何曲线或面线段都将根据采样距离参数值进行增密。 如果未定义采样距离,则该值将从输入表面派生。 对于栅格,默认采样大小将是栅格的像元大小。 对于 TIN、terrain 或 LAS 数据集,默认采样将基于三角化表面生成的边。 如果曲线短于采样大小,则曲线将使用起点和终点简化为两点线。

  • 使用自然邻域插值法时,请考虑指定采样距离,该距离应等于或大于表面中数据点平均点间距的一半。

  • 在使用仅插值折点参数时,不会输出折点落在表面数据区域外的要素,除非输入表面为栅格并且使用最邻近插值方法。

参数

标注说明数据类型
输入表面

将用于插入 Z 值的表面。

LAS Dataset Layer; Mosaic Layer; Raster Layer; Terrain Layer; TIN Layer; Image Service
输入要素

待处理的输入要素。

Feature Layer
输出要素类

将生成的要素类。

Feature Class
采样距离
(可选)

用于内插 z 值的间距。 默认情况下,该参数是栅格数据集的像元大小或三角化网格面的自然增密。

Double
Z 因子
(可选)

Z 值将乘上的系数。 此值通常用于转换 z 线性单位来匹配 x,y 线性单位。 默认值为 1,此时高程值保持不变。 如果输入表面的空间参考具有已指定线性单位的 z 基准,则此参数不可用。

Double
方法
(可选)

指定将用于确定输出要素的高程值的插值方法。 可用选项取决于表面类型。

  • 双线性将使用双线性插值法来确定查询点的值。 如果输入为栅格表面,则其为默认值。
  • 最邻近法将使用最邻近插值法来确定查询点的值。 如果使用此方法,则将仅针对输入要素的折点对表面值进行插值。 此选项仅适用于栅格表面。
  • 线性将根据由包含查询点 x,y 位置的三角形定义的平面获取高程值。 此为 TIN、terrain 和 LAS 数据集的默认插值方法。
  • 自然邻域法通过将基于区域的权重应用于查询点的自然邻域获取高程值。
  • 合并最小 Z 值将根据在查询点自然邻域中找到的最小 z 值获取高程值。
  • 合并最大 Z 值将根据在查询点自然邻域中找到的最大 z 值获取高程值。
  • 合并最近的 Z 值将根据在查询点自然邻域中找到的最邻近值获取高程值。
  • 合并最接近平均值的 z 值将根据距查询点所有自然邻域的平均值最近的 z 值获取高程值。
String
仅插值折点
(可选)

指定是否仅沿输入要素的折点进行插值,并忽略采样距离选项。 如果输入表面为栅格并且选择最邻近插值法,则只能在要素折点处插值 z 值。

  • 选中 - 插值将仅沿折点发生。
  • 未选中 - 将使用采样距离进行插值。 这是默认设置。
Boolean
金字塔等级分辨率
(可选)

将使用 terrain 金字塔等级的 z 容差或窗口大小分辨率。 默认值为 0,或全分辨率。

Double
保留部分位于表面外的要素
(可选)

指定是否将在输出中保留一个或多个折点落在栅格数据区域范围之外的要素。 此参数仅当输入表面为栅格并且使用最邻近插值法时可用。

  • 选中 - 落在栅格表面范围之外的每个折点都将具有各自的 z 值,此值派生自针对栅格表面内的折点计算的 z 值趋势,并将保留在输出中。
  • 未选中 - 将在输出中跳过至少一个折点落在栅格表面范围之外的要素。 这是默认设置。
Boolean

arcpy.ddd.InterpolateShape(in_surface, in_feature_class, out_feature_class, {sample_distance}, {z_factor}, {method}, {vertices_only}, {pyramid_level_resolution}, {preserve_features})
名称说明数据类型
in_surface

将用于插入 Z 值的表面。

LAS Dataset Layer; Mosaic Layer; Raster Layer; Terrain Layer; TIN Layer; Image Service
in_feature_class

待处理的输入要素。

Feature Layer
out_feature_class

将生成的要素类。

Feature Class
sample_distance
(可选)

用于内插 z 值的间距。 默认情况下,该参数是栅格数据集的像元大小或三角化网格面的自然增密。

Double
z_factor
(可选)

Z 值将乘上的系数。 此值通常用于转换 z 线性单位来匹配 x,y 线性单位。 默认值为 1,此时高程值保持不变。 如果输入表面的空间参考具有已指定线性单位的 z 基准,则此参数不可用。

Double
method
(可选)

指定将用于确定输出要素的高程值的插值方法。 可用选项取决于表面类型。

  • BILINEAR将使用双线性插值法来确定查询点的值。 如果输入为栅格表面,则其为默认值。
  • NEAREST将使用最邻近插值法来确定查询点的值。 如果使用此方法,则将仅针对输入要素的折点对表面值进行插值。 此选项仅适用于栅格表面。
  • LINEAR将根据由包含查询点 x,y 位置的三角形定义的平面获取高程值。 此为 TIN、terrain 和 LAS 数据集的默认插值方法。
  • NATURAL_NEIGHBORS通过将基于区域的权重应用于查询点的自然邻域获取高程值。
  • CONFLATE_ZMIN将根据在查询点自然邻域中找到的最小 z 值获取高程值。
  • CONFLATE_ZMAX将根据在查询点自然邻域中找到的最大 z 值获取高程值。
  • CONFLATE_NEAREST将根据在查询点自然邻域中找到的最邻近值获取高程值。
  • CONFLATE_CLOSEST_TO_MEAN将根据距查询点所有自然邻域的平均值最近的 z 值获取高程值。
String
vertices_only
(可选)

指定是否仅沿输入要素的折点进行插值,并忽略采样距离选项。

  • DENSIFY将使用采样距离进行插值。 这是默认设置。
  • VERTICES_ONLY插值将仅沿折点发生。
Boolean
pyramid_level_resolution
(可选)

将使用 terrain 金字塔等级的 z 容差或窗口大小分辨率。 默认值为 0,或全分辨率。

Double
preserve_features
(可选)

指定是否将在输出中保留一个或多个折点落在栅格数据区域范围之外的要素。 此参数仅当输入表面为栅格并且使用最邻近插值法时可用。

  • PRESERVE落在栅格表面范围之外的每个折点都将具有各自的 z 值,此值派生自针对栅格表面内的折点计算的 z 值趋势,并将保留在输出中。
  • EXCLUDE将在输出中跳过至少一个折点落在栅格表面范围之外的要素。 这是默认设置。
Boolean

代码示例

InterpolateShape 示例 1(Python 窗口)

下面的示例演示了如何在 Python 窗口中使用此工具。

arcpy.env.workspace = "C:/data"
arcpy.ddd.InterpolateShape("my_tin", "roads.shp", "roads_interp.shp")
InterpolateShape 示例 2(独立脚本)

下面的示例演示了如何在独立 Python 脚本中使用此工具。

'''*********************************************************************
Name: InterpolateShape Example
Description: This script demonstrates how to use InterpolateShape
             on all 2D features in a target workspace.
*********************************************************************'''
# Import system modules
import arcpy

# Set local variables
inWorkspace = arcpy.GetParameterAsText(0)
surface = arcpy.GetParameterAsText(1)

try:
    # Set default workspace
    arcpy.env.workspace = inWorkspace
    # Create list of feature classes in target workspace
    fcList = arcpy.ListFeatureClasses()
    if fcList:
        for fc in fcList:
            desc = arcpy.Describe(fc)
            # Find 2D features
            if not desc.hasZ:
                # Set Local Variables
                outFC = "{0}_3D.shp".format(desc.basename)
                method = "BILINEAR"
                # Execute InterpolateShape
                arcpy.ddd.InterpolateShape(surface, fc, outFC, 
                                           10, 1, method, True)
            else:
                print("{0} is not a 2D feature.".format(fc))
    else:
        print("No feature classes were found in {0}.".format(env.workspace))
    
except arcpy.ExecuteError:
    print(arcpy.GetMessages())
    
except Exception as err:
    print(err)

许可信息

  • Basic: 需要 3D Analyst
  • Standard: 需要 3D Analyst
  • Advanced: 需要 3D Analyst

相关主题