需要 Spatial Analyst 许可。
描述
根据落入每个单元一定半径范围内的折线要素计算每单位面积的量级。
插图
使用
计算密度时,将仅考虑邻域内部的线段。如果某特定像元的邻域范围内未落入任何线,则为该像元指定 NoData。
半径参数值越大,生成的密度栅格的概化程度便越高。值越小,生成的栅格所显示的信息越详细。
输出栅格的值将始终为浮点型。
输出像元大小可以通过数值进行定义,也可以从现有栅格数据集获取。如果没有将像元大小明确指定为参数值,则将从像元大小环境获取相应值(前提是已指定环境)。如果未指定参数像元大小和环境像元大小,但已设置捕捉栅格环境,则将使用捕捉栅格的像元大小。如果未指定任何内容,则像元大小会通过使用范围的宽度或高度中的较小值除以 250 来计算,其中范围位于在环境中指定的输出坐标系内。
如果使用数值指定像元大小,则工具会直接将其用于输出栅格。
如果使用栅格数据集指定像元大小,则该参数将显示栅格数据集的路径而不是像元大小的值。如果数据集的空间参考与输出空间参考相同,则栅格数据集的像元大小将直接用于分析。如果数据集的空间参考与输出空间参考不同,则将基于所选的像元大小投影方法进行投影。
如果面积单位比例因子的单位相对于要素(线段的长度)很小,则输出值可能会很小。要获取较大的值,请使用单位较大的面积单位比例因子(例如,平方千米与平方米)。
有关适用于此工具的地理处理环境的详细信息,请参阅分析环境和 Spatial Analyst。
语法
LineDensity(in_polyline_features, population_field, {cell_size}, {search_radius}, {area_unit_scale_factor})
参数 | 说明 | 数据类型 |
in_polyline_features | 要计算密度的输入线要素。 | Feature Layer |
population_field | 表示各折线总体值(线应被统计的次数)的数值字段。 Population 字段的值可以是整型或浮点型。 以下列出的是该字段的选项和默认特性。
| Field |
cell_size (可选) | 将创建的输出栅格的像元大小。 此参数可以通过数值进行定义,也可以从现有栅格数据集获取。如果未将像元大小明确指定为参数值,则将使用环境像元大小值(如果已指定);否则,将使用其他规则通过其他输出计算像元大小。有关详细信息,请参阅用法。 | Analysis Cell Size |
search_radius (可选) | 在其范围内计算密度的搜索半径。单位基于输出空间参考投影的线性单位。 例如,如果单位为米,若要包含一英里邻域内的所有要素,可将搜索半径设置为 1609.344(1 英里 = 1609.344 米)。 默认值为输出空间参考中输出范围的宽度或高度的最小值除以 30。 | Double |
area_unit_scale_factor (可选) | 输出密度值的所需面积单位。 基于输出空间参考的线性单位选择默认单位。若要转换密度输出,可将此单位更改为合适的单位。线密度值同时转换长度和面积单位。 未指定输出空间参考的情况下,输出空间参考与输入要素类相同。默认输出密度单位通过输出空间参考的线性单位确定,如下所示。如果输出线性单位是米,输出面积密度单位将设置为 SQUARE_KILOMETERS,输出平方千米(点要素)或千米每平方千米(折线要素)。如果输出线性单位是英尺,输出面积密度单位将设置为 SQUARE_MILES。 如果输出单位不是英尺和米,输出面积密度单位将设置为 SQUARE_MAP_UNITS。即,输出密度单位为输出空间参考的线性单位的平方。例如,如果输出线性单位是厘米,输出面积密度单位将是 SQUARE_MAP_UNITS,即平方厘米。如果输出线性单位是千米,输出面积密度单位将是 SQUARE_MAP_UNITS,即平方千米。 可用选项及相应的输出密度单位如下:
| String |
返回值
名称 | 说明 | 数据类型 |
out_raster | 输出线密度栅格。 总为浮点栅格。 | Raster |
代码示例
线密度分析 (LineDensity) 示例 1(Python 窗口)
本例将根据折线 shapefile 的长度字段计算密度栅格。
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outLDens = LineDensity("roads.shp", "LENGTH", 45, 1000, "SQUARE_MILES")
outLDens.save("C:/sapyexamples/output/ldensout")
线密度分析 (LineDensity) 脚本示例 2(独立脚本)
本例将根据折线 shapefile 的长度字段计算密度栅格。
# Name: LineDensity_Ex_02.py
# Description: Calculates a magnitude per unit area from polyline features
# that fall within a radius around each cell.
# 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
inPolylineFeatures = "roads.shp"
populationField = "length"
cellSize = 120
searchRadius = 1500
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute LineDensity
outLineDensity = LineDensity(inPolylineFeatures, populationField, cellSize,
searchRadius, "SQUARE_MILES")
# Save the output
outLineDensity.save("C:/sapyexamples/output/linedensity")
环境
许可信息
- Basic: 需要 Spatial Analyst
- Standard: 需要 Spatial Analyst
- Advanced: 需要 Spatial Analyst