Arcgis基于格网的道路密度分析

前言

在ArcGIS中，计算研究区域内各个格网的道路密度主要有以下两种方法：

• 密度插值法：先通过线密度分析、核密度分析计算区域内的道路密度，再通过建立格网和值提取至点等操作将密度值关联到格网上。
• 计数法： 通过计数每个格网内道路的数量，然后与格网面积相除得到道路密度。

朴素和可靠的计数法：

为什么说计数法是朴素和可靠的呢？
朴素在于它的计算方法是统计各个格网内的道路数量，那么对于没有道路经过的格网，其密度就是0。
这种朴素是相对于密度插值法而言的，核密度分析法中的密度是基于核函数计算的，其函数是个光滑的曲面，因此得到的函数值是经过插值法得到，因此对于没有道路经过的格网而言，它们也会有一定的密度值。
可靠在于这种计算方式很符合我们对密度计算的认识，表征为单位面积上的道路数量。其不需要经过复杂的算法，一切都看似十分可控。

那么扯了这么久，下面正式来通过Arcgis来实现一下。

实现过程

原始数据

行政区划：粤港澳大湾区
道路数据：来自于osm的2018年道路数据

环境准备

这里的准备包括以下几项：

• 数据框坐标系
  • 其坐标系默认与第一个导入的数据的坐标系保持一致
• 数据坐标系
  • 如果需要计算面积、长度的话尽量采用投影坐标系来计算，因为地理坐标系是球面坐标，计算起来误差有点大

横轴墨卡托投影带计算方法：
$$X=$$