插图
使用方法
使用自然邻域插值法时,请考虑指定采样距离,该距离应等于或大于表面中数据点平均点间距的一半。
在使用仅插值折点选项时,不会输出折点落在表面数据区域外的要素,除非输入表面为栅格并且正在使用最邻近插值方法。
语法
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 |
in_feature_class | 待处理的输入要素。 | Feature Layer |
out_feature_class | 将由此工具生成的要素类。 | Feature Class |
sample_distance (可选) | 用于内插 z 值的间距。默认情况下,该参数是栅格数据集的像元大小或三角化网格面的自然增密。 | Double |
z_factor (可选) | Z 值将乘上的系数。此值通常用于转换 Z 线性单位来匹配 XY 线性单位。默认值为 1,此时高程值保持不变。如果输入表面的空间参考具有已指定线性单位的 Z 基准,则将禁用此参数。 | Double |
method (可选) | 用于确定输出要素的高程值的插值方法。可用选项取决于正在使用的表面类型。
| String |
vertices_only (可选) | 指定是否仅沿输入要素的折点进行插值,从而忽略采样距离选项。
| Boolean |
pyramid_level_resolution (可选) | 此工具将使用 terrain 金字塔等级的 z 容差或窗口大小分辨率。默认值为 0(z 容差),或全分辨率(窗口大小)。 | Double |
preserve_features (可选) | 指定是否将在输出中保留一个或多个折点落在栅格数据区域范围之外的要素。此参数仅当输入表面为栅格并且使用最邻近插值法时可用。
| Boolean |
代码示例
下面的示例演示了如何在 Python 窗口中使用此工具。
arcpy.env.workspace = "C:/data"
arcpy.InterpolateShape_3d("my_tin", "roads.shp", "roads_interp.shp")
下面的示例演示了如何在独立 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