插值 Shape (3D Analyst)
插图
![]() |
用法
-
在使用自然邻域法插值选项时,请确保指定了合理的采样距离。该距离通常是用于构建 TIN 或 terrain 数据集的数据的平均点间距的 0.5 倍至 1.0 倍。
-
在使用仅插值折点选项时,将忽略折点落在表面数据区域外的输入要素,并且不会输出它们。请先裁剪要素,再运行插值 Shape,以确保要素完全位于表面上。
语法
InterpolateShape_3d (in_surface, in_feature_class, out_feature_class, {sample_distance}, {z_factor}, {method}, {vertices_only}, {pyramid_level_resolution})
参数 | 说明 | 数据类型 |
in_surface |
用于内插 z 值的 LAS 数据集、栅格、TIN 或 terrain 表面。 | LAS Dataset Layer, Raster Layer; Terrain Layer; TIN Layer |
in_feature_class |
待输入的要素类。 | Feature Layer |
out_feature_class |
输出要素类。 | Feature Class |
sample_distance (可选) |
用于内插 z 值的间距。默认情况下,该参数是栅格的像元大小或 TIN 的自然增密。 | Double |
z_factor (可选) |
高程值将乘上的系数。此系数通常用于转换 Z 线性单位,以匹配 XY 线性单位的值。默认值为 1,此时高程值保持不变。 | Double |
method (可选) |
用于确定输出要素的高程值的插值方法。可用选项取决于正在使用的表面类型。BILINEAR 插值可用于栅格表面,查询点根据四个最邻近的像元中找到的值获取其高程。Terrain 和 TIN 数据集提供了以下选项:
| String |
vertices_only (可选) |
指定是否仅沿输入要素的折点进行插值,从而忽略采样距离选项。
| Boolean |
pyramid_level_resolution (可选) |
此工具将使用 terrain 金字塔等级的 z 容差或窗口大小分辨率。默认值为 0(z 容差),或全分辨率(窗口大小)。 | Double |
代码实例
插值 Shape (InterpolateShape) 示例 1(Python 窗口)
下面的示例演示了如何在 Python 窗口中使用此工具:
import arcpy
from arcpy import env
arcpy.CheckOutExtension("3D")
env.workspace = "C:/data"
arcpy.InterpolateShape_3d("my_tin", "roads.shp", "roads_interp.shp")
插值 Shape (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
from arcpy import env
import exceptions, sys, traceback
# Set local variables
inWorkspace = arcpy.GetParameterAsText(0)
surface = arcpy.GetParameterAsText(1)
try:
arcpy.CheckOutExtension("3D")
# Set default workspace
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)
arcpy.CheckInExtension('3D')
except arcpy.ExecuteError:
print arcpy.GetMessages()
except:
# Get the traceback object
tb = sys.exc_info()[2]
tbinfo = traceback.format_tb(tb)[0]
# Concatenate error information into message string
pymsg = 'PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}'\
.format(tbinfo, str(sys.exc_info()[1]))
msgs = 'ArcPy ERRORS:\n {0}\n'.format(arcpy.GetMessages(2))
# Return python error messages for script tool or Python Window
arcpy.AddError(pymsg)
arcpy.AddError(msgs)
相关主题
许可信息
ArcGIS for Desktop Basic:需要 3D Analyst
ArcGIS for Desktop Standard:需要 3D Analyst
ArcGIS for Desktop Advanced:需要 3D Analyst
9/15/2013