核密度分析的工作原理

需要 Spatial Analyst 许可。

核密度分析工具用于计算要素在其周围邻域中的密度。此工具既可计算点要素的密度，也可计算线要素的密度。

可能的用途包括针对社区规划分析房屋密度或犯罪情况，或探索道路或公共设施管线如何影响野生动物栖息地。可使用 population 字段赋予某些要素比其他要素更大的权重，该字段还允许使用一个点表示多个观察对象。例如，一个地址可以表示一栋六单元的公寓，或者在确定总体犯罪率时可赋予某些罪行比其他罪行更大的权重。对于线要素，分车道高速公路可能比狭窄的土路具有更大的影响。

如何计算核密度

对于点要素

核密度分析用于计算每个输出栅格像元周围的点要素的密度。

概念上，每个点上方均覆盖着一个平滑曲面。在点所在位置处表面值最高，随着与点的距离的增大表面值逐渐减小，在与点的距离等于搜索半径的位置处表面值为零。仅允许使用圆形邻域。曲面与下方的平面所围成的空间的体积等于此点的 Population 字段值，如果将此字段值指定为 NONE 则体积为 1。每个输出栅格像元的密度均为叠加在栅格像元中心的所有核表面的值之和。核函数以 Silverman 的著作（1986 年版，第 76 页，equation 4.5）中描述的四次核函数为基础。

如果 population 字段设置使用的是除 NONE 之外的值，则每项的值用于确定点被计数的次数。例如，值 3 会导致点被算作三个点。值可以为整型也可以为浮点型。

默认情况下，单位是根据输入点要素数据的投影定义的线性单位进行选择的，或是在输出坐标系环境设置中以其他方式指定的。如果选择的是面积单位，则计算所得的像元密度将乘以相应因子，然后写入到输出栅格。

例如，如果输入单位为米，则输出面积单位将默认为平方千米。将以米和千米为单位的单位比例因子进行比较，将得到相差 1,000,000（1,000 米 x 1,000 米）倍的值。

对于线要素

核密度分析还可用于计算每个输出栅格像元的邻域内的线状要素的密度。

概念上，每条线上方均覆盖着一个平滑曲面。其值在线所在位置处最大，随着与线的距离的增大此值逐渐减小，在与线的距离等于指定的搜索半径的位置处此值为零。由于定义了曲面，因此曲面与下方的平面所围成的空间的体积等于线长度与 Population 字段值的乘积。每个输出栅格像元的密度均为叠加在栅格像元中心的所有核表面的值之和。用于线的核函数是根据 Silverman 著作中所述的用于计算点密度的四次核函数改编的。

一条线段与覆盖在其上方的核表面。

上图显示的是一条线段与覆盖在其上方的核表面。线段对密度的影响等同于栅格像元中心处核表面的值对密度的影响。

默认情况下，单位是根据输入折线要素数据的投影定义的线性单位进行选择的，或是在输出坐标系环境设置中以其他方式指定的。

指定输出面积单位因子后，它会转换长度单位和面积单位。例如，如果线性单位是米，则输出面积单位将默认为 SQUARE_KILOMETERS ，而所得到的线密度单位将转换为千米/平方千米。将以米和千米为单位的面积比例因子相比较，最终结果将是相差 1,000 倍的密度值。

您可以通过手动选择相应的因子来控制密度单位。要将密度单位设置为米/平方米（而不是默认的千米/平方千米），请将面积单位设置为 SQUARE_METERS。同样，若要将输出的密度单位设置为英里/平方英里，请将面积单位设置为 SQUARE_MILES。

有关特定距离单位的详细信息，请参阅以下主题。

• 了解有关 ArcGIS 中单位和转换常数的详细信息

如果 population 字段使用的是除 NONE 之外的值，则线的长度将由线的实际长度乘以此线的 population 字段的值而得出。

计算核密度的公式

以下公式定义如何计算点的核密度以及如何在核密度公式中确定默认搜索半径。

预测点的密度

新 (x, y) 位置的预测密度由以下公式确定：

其中：

• i = 1,…，n 是输入点。如果它们位于 (x, y) 位置的半径距离内，则仅包括总和中的点。
• pop_i 是 I 点的 population 字段值，它是一个可选参数。
• dist_i 是点 i 和 (x, y) 位置之间的距离。

然后将计算的密度乘以点数，或者乘以 population 字段的总和（如果有）。这种校正使空间配额等于点数（或总和或 population 字段），而非总等于 1。该实现使用四次核（Silverman，1986）。需要为要估算密度的每个位置分别计算公式。由于正在创建栅格，因此计算将应用于输出栅格中每个像元的中心。

默认搜索半径（带宽）

在 ArcGIS 10.2.1 中，计算默认搜索半径的方法已得到了改善。

用于确定默认搜索半径（同时也被称为带宽）的算法如下所示：

1. 计算输入点的平均中心。如果提供了 Population 字段，则此字段和所有以下计算将按该字段中的值加权。
2. 计算所有点的（加权）平均中心的距离。
3. 计算这些距离的（加权）中值 D_m。
4. 计算（加权）标准距离 SD。

   有关此内容的详细信息，请参阅标准距离空间统计工具。

5. 使用以下公式计算带宽。

其中：

• D_m 是（加权）平均中心的（加权）中值距离。
• n 是没有使用 population 字段的点数，如果提供了 population 字段，则 n 是 population 字段值的总和。
• SD 是标准距离。

请注意，方程的 min 部分表示将使用两选项（SD 或 ）中计算得出值较小者。

有两种计算标准距离的方法，即未加权和加权。

未加权距离

其中：

• x _i 、y _i 和 z _i 是要素 i 的坐标
• {x̄, ȳ, z̄} 表示要素的平均中心
• N 等于要素总数。

加权距离

其中：

• w_i 是要素 i 的权重
• {x _w, y _w, z _w} 表示加权平均中心

方法

这种用于选择搜索半径的方法基于“Silverman 经验规则”带宽估算公式，但该公式经调整已适用于两个维度。这种计算默认半径的方法通常可避免稀疏数据集中经常出现的“点周围的圆环”现象，并且防止空间异常值 - 即远离其余点的几个点。

重新创建来自 ArcGIS 先前版本的结果

在 ArcGIS 10.2.1 之前的版本中，默认搜索半径是通过将输入范围的宽度或高度中的较小者除以 30 来计算的。如果您希望重现这些结果，请将默认值输出像元大小参数乘以 25/3 或 8.333333。

参考文献

Silverman, B. W. Density Estimation for Statistics and Data Analysis. 纽约：Chapman and Hall，1986 年