需要 Spatial Analyst 许可。
摘要
通过托布勒函数的倒数定义垂直成本系数与垂直相对移动角度之间的关系,从而得出以小时为单位的徒步旅行时间。
插图
说明
VfHikingTime 对象用于距离累积和距离分配 Spatial Analyst 工具,以及路径距离、路径距离分配和路径距离回溯链接旧的距离工具。
垂直系数 (VF) 对象用于定义垂直成本系数和垂直相对移动角度 (VRMA) 之间的关系。
VF 定义从一个像元移至下一像元时所遇到的垂直阻力。
VRMA 标识“起始”像元或处理像元与“终止”像元之间的坡度角度。
当 VRMA 大于切削角上限且小于切削角下限时,在两个像元之间移动的 VF 设置为源自函数的值。 如果 VRMA 大于 highCutAngle 值或小于 lowCutAngle 值,则 VF 设置为无穷大。
托布勒徒步旅行函数的公式如下:
VfHikingTime 的公式(即托布勒徒步旅行函数的倒数)如下所示:
其中 d 是用于计算坡度的角度。
注:
有关托布勒徒步旅行函数的其他信息,请参阅地理处理工具帮助中的使用垂直系数调整遇到的距离和以下资源:
- Tobler, Waldo. 1993. 地理分析与模型三篇论文:非各向同性地理建模;地理学的几何特征的猜测;全球空间分析 (93-1)。 检索自 https://escholarship.org/uc/item/05r820mz
语法
VfHikingTime ({lowCutAngle}, {highCutAngle})
参数 | 说明 | 数据类型 |
lowCutAngle | The VRMA degree defining the lower threshold, below which (less than) the VFs are set to infinity. (默认值为 -70.0) | Double |
highCutAngle | The VRMA degree defining the upper threshold, beyond which (larger than) the VFs are set to infinity. (默认值为 70.0) | Double |
属性
属性 | 说明 | 数据类型 |
lowCutAngle (可读写) | 用于定义阈值下限的 VRMA 度数,如果低于(小于)该值,则将 VF 设置为无穷大。 | Double |
highCutAngle (可读写) | 用于定义阈值上限的 VRMA 度数,如果高于(大于)该值,则将 VF 设置为无穷大。 | Double |
代码示例
此代码示例演示了如何创建 VfHikingTime 类并在 Python 窗口中的距离累积工具中使用它。
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
myVerticalFactor = VfHikingTime(-30, 30)
outDistAccum = DistanceAccumulation("Source.shp", "", "elev.tif",
"cost.tif", "elev.tif",
myVerticalFactor)
outDistAccum.save("C:/sapyexamples/output/distAccumVfH.tif")
使用距离累积工具和 VfHikingTime 类执行距离累积分析。
# Name: VfHikingTime_Ex_02.py
# Description: Uses the VfHikingTime object to run the
# DistanceAccumulation tool
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inSourceData = "source.shp"
inCostRaster = "costraster.tif"
inElevation = "elev.tif"
# Create the VfHikingTime Object
lowCutAngle = -30
highCutAngle = 30
myVerticalFactor = VfHikingTime(lowCutAngle, highCutAngle)
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Run DistanceAccumulation
outDistAccum = DistanceAccumulation(inSourceData, "", inElevation,
inCostRaster, inElevation,
myVerticalFactor)
# Save the output
outDistAccum.save("C:/sapyexamples/output/distaccumvfH2.tif")