核密度分析 (空间分析)

许可等级:BasicStandardAdvanced

摘要

使用核函数根据点或折线 (polyline) 要素计算每单位面积的量值以将各个点或折线 (polyline) 拟合为光滑锥状表面。

了解有关核密度分析工作原理的详细信息

插图

Kernel Density illustration
OutRas = KernelDensity(InPts, None, 30)

用法

语法

KernelDensity (in_features, population_field, {cell_size}, {search_radius}, {area_unit_scale_factor})
参数说明数据类型
in_features

要计算密度的输入要素(点或线)。

Feature Layer
population_field

表示各要素的 population 值的字段。Population 字段表示遍布于用来创建连续表面的景观内的计数或数量。

population 字段的值可以是整型或浮点型。

以下列出的是该字段的选项和默认特性。

  • 如果不使用任何项目或特殊值,则选择 None,这样每一要素就只计数一次。

  • 如果输入要素包含 Z,则可以使用 形状

  • 否则,默认字段为 POPULATION。以下条件同样适用。

    • 如没有 POPULATION 字段,但是存在 POPULATIONxxxx 字段,默认使用此字段。“xxxx”可以为任何有效字符,例如 POPULATION6POPULATION1974POPULATIONROADTYPE
    • 如没有 POPULATION 字段或 POPULATIONxxxx 字段,但是存在 POP 字段,默认使用此字段。
    • 如没有 POPULATIONPOPULATIONxxxxPOP 字段,但是存在 POPxxxx 字段,默认使用此字段。
    • 如没有 POPULATIONPOPULATIONxxxxPOPPOPxxxx 字段,则默认使用 NONE
Field
cell_size
(可选)

输出栅格数据集的像元大小。

如果专门进行设置,则是环境中的值。如果未设置环境,单元大小为输出空间参考中输出范围的宽度或高度较小值除以 250。

Analysis Cell Size
search_radius
(可选)

在其范围内计算密度的搜索半径。单位基于输出空间参考投影的线性单位。

例如,如果单位为米,若要包含一英里邻域内的所有要素,可将搜索半径设置为 1609.344(1 英里 = 1609.344 米)。

使用“Silverman 经验规则”的空间变量专为输入数据集计算默认搜索半径(带宽),该变量可有效避免空间异常值(即距离其余点太远的点)。有关该算法的描述,请参阅上述“使用提示”。

Double
area_unit_scale_factor
(可选)

输出密度值的所需面积单位。

基于输出空间参考投影的线性单位选择默认单位。若要转换密度输出,可将此单位更改为合适的单位。线密度值同时转换长度和面积单位。

例如,如果输入单位是米,则默认输出面积密度单位为平方米(点要素)或千米每平方千米(折线要素)。

基于输入要素单位的默认密度单位有:

  • SQUARE_MAP_UNITS 如果单位为未知,点或十进制度。
  • SQUARE_MILES 用于英尺、码、英里或海里。
  • SQUARE_KILOMETERS 用于米或千米。
  • SQUARE_INCHES 用于英寸。
  • SQUARE_CENTIMETERS 用于厘米。
  • SQUARE_MILLIMETERS 用于毫米。
String

返回值

名称说明数据类型
out_raster

输出核密度栅格。

其总为浮点栅格。

Raster

代码实例

核密度 (KernelDensity) 示例 1(Python 窗口)

本示例根据点 shapefile 来计算平滑的密度栅格。

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outKDens = KernelDensity("rec_sites.shp", "NONE", 45, 1200, "SQUARE_KILOMETERS")
outKDens.save("C:/sapyexamples/output/kdensout")
核密度 (KernelDensity) 示例 2(独立脚本)

本示例根据点 shapefile 来计算平滑的密度栅格。

# Name: KernelDensity_Ex_02.py
# Description: Calculates a magnitude per unit area from point or polyline 
#    features using a kernel function to fit a smoothly tapered 
#    surface to each point or polyline.
# 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
inFeatures = "rec_sites.shp"
populationField = "NONE"
cellSize = 60
searchRadius = 2500


# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")

# Execute KernelDensity
outKernelDensity = KernelDensity(inFeatures, populationField, cellSize,
                                 searchRadius, "SQUARE_KILOMETERS")

# Save the output 
outKernelDensity.save("C:/sapyexamples/output/kerneldout")

环境

相关主题

许可信息

ArcGIS for Desktop Basic: 需要 Spatial Analyst
ArcGIS for Desktop Standard: 需要 Spatial Analyst
ArcGIS for Desktop Advanced: 需要 Spatial Analyst
5/10/2014