构造视线 (3D Analyst)
摘要
创建表示视线（从一个或多个视点到目标要素类的要素）的线要素。
用法
-
如果目标要素是线或面，则对周长进行统一采样。
-
输入要素类必须是点要素类。多点无效。
-
如果为观察点和目标要素指定了高度源，将生成三维输出。
使用连接字段为给定观察点指定一个或多个目标。如果未使用连接字段，则将所有点连接到所有目标。
观察点和目标要素的高度源默认为以下列表中遇到的第一个字段名称：
- Shape.Z（只对启用 Z 值的要素可用）
- Spot
- Z
- Z_Value
- Height
- Elev
- Elevation
- Contour
如果不存在适合的高度字段，则默认情况下使用 <无> 关键字来指示要素不含 Z 值。
如果所需高度字段在默认字段选择中没有更高优先级，将需要显式指定所需字段。同样，如果某个高度字段不是所需字段但要素类包含默认选择列表中的字段之一，则需要指定 <无> 关键字。
应按输出要素类的 x,y 单位给出采样距离单位。
以下字段将会添加到包含视线的输出要素类：
- OID_OBSERV - 观察点的 OID
- OID_TARGET - 目标要素的 OID
- DIST_ALONG - 当目标要素为线或面时沿此目标要素的距离
语法
| 参数 | 说明 | 数据类型 |
in_observer_points |
表示观察点的单点要素。不支持多点要素。 | Feature Layer |
in_target_features |
目标要素（点、多点、线和面）。 | Feature Layer |
out_line_feature_class |
包含视线的输出要素类。 | Feature Class |
observer_height_field (可选) |
从观察点属性表获得的观察点高度值的源。 从按优先级顺序列出的以下选项中选择默认的观察者高度字段字段。如果存在多个字段并且所需字段在默认字段选择上没有更高优先级，将需要指定所需字段。如果不存在适合的高度字段，将使用 <无> 关键字。同样，如果某个高度字段不是所需字段但要素类具有下列字段中的一个，则需要指定 <无> 关键字。
| String |
target_height_field (可选) |
目标的高度字段。 从按优先级顺序列出的以下选项中选择默认的目标高度字段字段。如果存在多个字段并且所需字段在默认字段选择上没有更高优先级，将需要指定所需字段。如果不存在适合的高度字段，将使用 <无> 关键字。同样，如果某个高度字段不是所需字段但要素类具有下列字段中的一个，则需要指定 <无> 关键字。如果不存在适合的高度字段，则默认情况下将使用 <无> 关键字。
| String |
join_field (可选) |
使用连接字段将观察点与特定目标相匹配。 | String |
sample_distance (可选) |
目标为线或面要素类时采样之间的距离。按输出要素类的 XY 单位解释采样距离单位。 | Double |
代码实例
下面的示例演示了如何在 Python 窗口中使用此工具：
import arcpy
from arcpy import env
arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'
arcpy.ConstructSightLines_3d('observer_pt.shp', 'target.shp',
'sightlines.shp', 'BASEHEIGHT',
'TOP_HEIGHT', 'NAME')
下面的示例演示了如何在独立 Python 脚本中使用此工具：
'''*********************************************************************
Name: Sight Line Visibility
Description: This script demonstrates how to calculate visibility
for sight lines against the obstructions presented by
terrain elevation & building models in a multipatch.
*********************************************************************'''
# Import system modules
import arcpy
import exceptions, sys, traceback
from arcpy import env
try:
arcpy.CheckOutExtension('3D')
# Set Local Variables
env.workspace = 'C:/data'
obs_pts = "Observers.shp"
target = "Observation_Targets.shp"
sight_lines = "in_memory/sightlines"
surface = "sample.gdb/elevation/terrain"
buildings = "city_buildings.shp"
outLOS = arcpy.CreateUniqueName("Line_of_Sight.shp")
obstruction_pts = arcpy.CreateUniqueName("Obstruction_Points.shp")
arcpy.AddMessage("Constructing sight lines...")
arcpy.ddd.ConstructSightLines(obs_pts, target, sight_lines)
arcpy.AddMessage("Calculating line of sight...")
arcpy.ddd.LineOfSight(surface, sight_lines, outLOS, obstruction_pts,
"CURVATURE", "REFRACTION", 0.35, 0, buildings)
arcpy.GetMessages(0)
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)