太阳阴影体 (3D Analyst)

摘要

创建封闭体积,对每个要素在给定日期和时间使用阳光投射的阴影进行建模。

使用情况

  • 所有输入要素都应位于同一区域,因为相对太阳位置的计算基于第一个要素在第一个要素类中的位置。

  • 起始日期和时间结束日期和时间参数中分别提供一个日期,即可确定阴影模拟的日出和日落情况。 如果太阳在给定的日期和时间不可见,或者如果太阳的相对位置与输入要素成 90 度垂直角,则不会生成阴影体。

  • 通过向太阳光的方向拉伸输入要素,阴影将被模拟成一个闭合多面体。 光线均视为平行光线并且按太阳的相对位置计算出的方向进行传播。每个阴影体都在一个垂直于太阳光线水平投影的垂直平面开始和结束。

  • 以下字段将归因于阴影体要素:

    • SOURCE - 投射阴影体的要素类的名称。
    • SOURCE_ID - 投射阴影体的要素的唯一 ID。
    • DATE_TIME - 用于计算太阳位置的当地日期和时间。
    • AZIMUTH - 正北方向与太阳和地平线的相对位置的垂直投影之间的角度(以度为单位)。 其值的范围为 0 到 360 之间。
    • VERT_ANGLE - 地平线与太阳相对位置之间的角度(以度为单位),其中 0 度表示太阳处于地平线的位置,90 度表示太阳位于地平线的正上方。
    注:

    通常,每个阴影体都应紧靠或紧贴其投影源要素。 如果无法以这种方式生成阴影,将会从要素的外部范围边界处进行创建。 如果以这种方式至少创建了一个阴影,则会包括一个名为 HUGS_FEATR 的字段,该字段会指示哪些阴影紧靠其对应要素。

参数

标注说明数据类型
输入要素

用于模拟阴影的多面体要素。

Feature Layer
起始日期和时间

为阴影建模计算日光轨迹的日期和时间。

Date
输出要素类

将存储生成的阴影体的多面体要素类。

Feature Class
调整为夏令时
(可选)

指定是否针对夏令时 (DST) 调整时间值。

  • 未选中 - 不遵守 DST。 这是默认设置。
  • 选中 - 遵守 DST。
Boolean
时区
(可选)

参与输入所在的时区。 默认设置是操作系统设置的时区。

String
结束日期和时间
(可选)

用于计算太阳位置的最终日期和时间。 如果仅提供日期,则最终时间假定为日落。

Date
迭代间隔
(可选)

该值用于定义从开始日期起的迭代时间。

Double
迭代单位
(可选)

定义应用到起始日期和时间的迭代值的单位。

  • 迭代值将表示天数。 这是默认设置。
  • 小时迭代值将代表一个或多个小时。
  • 迭代值将代表一个或多个分钟。
String

arcpy.ddd.SunShadowVolume(in_features, start_date_and_time, out_feature_class, {adjusted_for_dst}, {time_zone}, {end_date_and_time}, {iteration_interval}, {iteration_unit})
名称说明数据类型
in_features
[in_features,...]

用于模拟阴影的多面体要素。

Feature Layer
start_date_and_time

为阴影建模计算日光轨迹的日期和时间。

Date
out_feature_class

将存储生成的阴影体的多面体要素类。

Feature Class
adjusted_for_dst
(可选)

指定是否针对夏令时 (DST) 调整时间值。

  • ADJUSTED_FOR_DST遵守 DST。
  • NOT_ADJUSTED_FOR_DST不遵守 DST。 这是默认设置。
Boolean
time_zone
(可选)

参与输入所在的时区。 默认设置是操作系统设置的时区。

String
end_date_and_time
(可选)

用于计算太阳位置的最终日期和时间。 如果仅提供日期,则最终时间假定为日落。

Date
iteration_interval
(可选)

该值用于定义从开始日期起的迭代时间。

Double
iteration_unit
(可选)

定义应用到起始日期和时间的迭代值的单位。

  • DAYS迭代值将表示天数。 这是默认设置。
  • HOURS迭代值将代表一个或多个小时。
  • MINUTES迭代值将代表一个或多个分钟。
String

代码示例

SunShadowVolume 示例 1(Python 窗口)

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

arcpy.env.workspace = 'C:/data'
arcpy.ddd.SunShadowVolume('sample.fgdb/buildings', 
                         start_date_and_time='12/25/2011 10:00 AM', 
                         out_feature_class='shadows_dec25.shp', 
                         adjusted_for_dst='ADJUSTED_FOR_DST', 
                         time_zone='Eastern_Standard_Time',
                         end_date_and_time='12/25/2011 3:00 PM', 
                         iteration_interval='HOURS', iteration_unit=1)
SunShadowVolume 示例 2(独立脚本)

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

'''*********************************************************************
Name: Model Shadows For GeoVRML Models
Description: Creates a model of the shadows cast by GeoVRML models
             imported to a multipatch feature class for a range of dates
             and times. A range of times from the start time and end
             time can also be specified by setting the EnforceTimes
             Boolean to True. This sample is designed to be used in a
             script tool.
*********************************************************************'''
# Import system modules
import arcpy
from datetime import datetime, time, timedelta

#*************************  Script Variables  **************************
inFiles = arcpy.GetParameterAsText(0) # list of input features
spatialRef = arcpy.GetParameterAsText(1) # list of GeoVRML files
outFC = arcpy.GetParameterAsText(2) # multipatch from 3D files
inTimeZone = arcpy.GetParameterAsText(3) # time zone
startDate = arcpy.GetParameter(4) # starting date as datetime
endDate = arcpy.GetParameter(5) # ending date as datetime
dayInterval = arcpy.GetParameter(6) # day interval as long (0-365)
minInterval = arcpy.GetParameter(7) # minute interval as long (0-60)
enforceTime = arcpy.GetParameter(8) # minute interval as Boolean
outShadows = arcpy.GetParameterAsText(9) # output shadow models
outIntersection = arcpy.GetParameterAsText(10) # shadow & bldg intersection

# Function to find all possible date/time intervals for shadow modelling
def time_list():
    dt_result = [startDate]
    if dayInterval:
        if endDate: #Defines behavior when end date is supplied
            while startDate < endDate:
                startDate += timedelta(days=dayInterval)
                dt_result.append(startDate)
            dt_result.append(endDate)
        else: # Behavior when end date is not given
            daymonthyear = datetime.date(startDate)
            while startDate <= datetime(daymonthyear.year, 12, 31, 23, 59):
                startDate += timedelta(days=dayInterval)
                dt_result.append(startDate)
    return dt_result


importFC = arcpy.CreateUniqueName('geovrml_import', 'in_memory')

# Import GeoVRML files to in-memory feature
arcpy.ddd.Import3DFiles(inFiles, importFC, 'ONE_FILE_ONE_FEATURE',
                        spatialRef, 'Z_IS_UP', 'wrl')

# Ensure that building models are closed
arcpy.ddd.EncloseMultiPatch(importFC, outFC, 0.05)

# Discard in-memory feature
arcpy.management.Delete(importFC)
dt_result = time_list()
for dt in dt_result:
    if dt == dt_result[0]:
        shadows = outShadows
    else:
        shadows = arcpy.CreateUniqueName('shadow', 'in_memory')
    arcpy.ddd.SunShadowVolume(outFC, dt, shadows, 'ADJUST_FOR_DST',
                              inTimeZone, '', minInterval, 'MINUTES')
    if dt is not dt_result[0]:
        arcpy.management.Append(shadows, outShadows)
        arcpy.management.Delete(shadows)
arcpy.ddd.Intersect3D(outFC, outIntersection, outShadows, 'SOLID')

许可信息

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

相关主题