ArcGIS缓冲区分析入门：学校 500 米服务范围，距离单位依赖坐标系

一、为什么“500 米”会算错：坐标系是第一性原理

1.1 一个被低估的入门陷阱

在 ArcGIS Pro 或 ArcMap 里，Buffer 工具的参数面板极其简洁：输入要素、输出要素类、距离、融合类型。很多人第一次做学校服务区时，直接拖入一个 WGS84 经纬度的学校点图层，距离填 500，单位选 Meters，点运行，得到一个看起来“圆圆的”多边形，任务就算完成了。

问题恰恰出在这里。当输入数据的坐标系是地理坐标系（GCS）时，ArcGIS 的缓冲距离单位默认按输入坐标系的角度单位解释，也就是“度”。你填的 500，在软件内部可能被当作 500 度来处理，或者被自动换算成一个与真实地面距离严重不符的值。结果就是缓冲区要么小得几乎看不见，要么大得离谱。这不是软件 bug，而是坐标系语义没有被正确理解。

Esri 官方文档在 Buffer 工具说明中明确指出：距离参数的单位取决于输入要素的坐标系；对于地理坐标系数据，建议先投影到合适的投影坐标系再进行距离分析（Esri, ArcGIS Pro Tool Reference, Buffer）。本文评述：这条建议看似平淡，实则是整个缓冲区分析的第一性原理——距离分析必须在“距离有意义”的空间中进行。

1.2 本文的分析主线

为了让全文逻辑不散，笔者确立一条贯穿始终的主线：

距离单位依赖坐标系 → 坐标系选择依赖分析目的 → 分析目的决定缓冲模型（欧氏/网络/三维）。

这条主线把“入门操作”和“进阶建模”串成一条链：先解决“500 米到底怎么算”，再解决“500 米算得对不对”，最后解决“500 米是不是我们真正想表达的服务范围”。

1.3 服务范围分析的现实意义

学校服务范围（school catchment / service area）是教育资源配置、学区划分、通学安全评估中的基础空间单元。国内《城市普通中小学校校舍建设标准》等文件对学校服务半径有原则性要求，各地在“15 分钟生活圈”“就近入学”等政策实践中也大量使用服务半径概念。国际上，OECD 与 UNESCO 的相关研究长期关注学校可达性与教育公平的关系（OECD, 2023; UNESCO, 2022）。

但政策文件里的“500 米”通常是沿道路的步行距离，而不是直线距离。这就埋下了第二个问题：欧氏缓冲区给出的是直线半径圆，而真实服务范围受路网、河流、围墙、铁路等阻隔影响。笔者认为，入门阶段用欧氏缓冲建立直觉是合理的，但必须清楚它的边界——它是“直线可达性”的近似，不是“步行可达性”的等价物。

二、地理坐标系 vs 投影坐标系：距离单位从哪里来

2.1 地理坐标系：用角度描述位置

地理坐标系（Geographic Coordinate System, GCS）以椭球体为基准，用经纬度（度、分、秒或十进制度）描述地球表面位置。常见的有 WGS84（EPSG:4326）、CGCS2000（EPSG:4490）、北京54、西安80等。它的坐标单位是角度，不是长度。

关键点在于：1 度经度对应的地面距离随纬度变化。在赤道附近约 111.32 公里，在纬度 30° 处约 96.5 公里，在纬度 60° 处约 55.8 公里。1 度纬度相对稳定，约 110.57 公里（因椭球扁率略有变化）。这意味着，同样“1 度”的缓冲距离，在哈尔滨和广州代表的地面长度完全不同。

国家大地坐标系 CGCS2000 自 2008 年起在我国全面启用，其椭球参数与 WGS84 极为接近，但在厘米级应用中仍需注意两者差异（中国国家标准化管理委员会, GB/T 22021-2008）。本文评述：对缓冲区分析而言，GCS 的正确用途是“存储与交换”，而不是“测量与缓冲”。

2.2 投影坐标系：把球面摊平成平面

投影坐标系（Projected Coordinate System, PCS）= 地理坐标系 + 投影方法 + 线性单位。它把椭球面投影到平面上，坐标单位变为米、英尺等长度单位。只有在这种坐标系下，“500 米”才是一个可以直接用于缓冲的、有明确物理意义的量。

我国常用的投影坐标系包括：

• 高斯-克吕格投影（Gauss-Krüger）：我国大比例尺地形图的标准投影，按 3° 或 6° 分带，如 CGCS2000 3-degree Gauss-Kruger zone 39（EPSG:4547）。
• UTM（通用横轴墨卡托）：国际通用，6° 分带，如 WGS84 UTM Zone 50N（EPSG:32650）。
• Web Mercator（EPSG:3857）：Web 地图事实标准，但面积和距离变形随纬度急剧增大，不适合做距离分析。
• Albers 等积投影：适合大范围面积统计，但距离保持性一般。

投影的核心矛盾是：没有任何一种投影能同时保持形状、面积、距离、方向全部正确。等距投影（equidistant）只在特定方向上保持距离正确，横轴墨卡托在中央经线上距离无变形，离中央经线越远变形越大。这就是为什么选投影必须结合研究区位置和分析目的。

2.3 距离单位依赖坐标系：三个层次

把“距离单位依赖坐标系”拆开看，至少有三个层次：

本文评述：很多教程只讲到第一层“要投影”，却忽略了第二、三层。真正专业的做法是：先明确距离的物理含义，再选择能最小化该含义误差的坐标系与工具。

2.4 投影变形量级：用数据说话

以高斯-克吕格 3° 带为例，中央经线处长度比恒为 1，离中央经线 1.5°（约 150 公里）处长度变形约为 1/10000 量级，即 500 米距离的变形约 5 厘米。这个量级对学校服务区分析完全可以忽略。

但如果误用 Web Mercator，情况就完全不同。Web Mercator 在纬度 φ 处的长度变形因子约为 1/cos(φ)。在纬度 30° 处约 1.155，在纬度 40° 处约 1.305，在纬度 60° 处约 2.0。也就是说，在北京（约北纬 40°）用 Web Mercator 做 500 米缓冲，实际地面距离可能只有约 383 米，误差超过 23%。

下表为模拟计算，用于说明量级（基于球面近似公式，非实测）：

本文评述：这张表是本文最想强调的“工程直觉”之一——用错坐标系，500 米可能变成 350 米，服务范围直接缩水三成。对教育资源配置这类涉及公平性的分析，这种误差是不可接受的。

三、学校 500 米服务区：完整操作路径

3.1 数据准备与检查

假设我们有一份某市小学点位数据（点要素）和一份行政区划面数据。数据来源可以是教育部门公开的学校名录、OpenStreetMap 的 amenity=school 要素，或高德/百度 POI（需注意坐标系偏移问题）。

第一步不是急着缓冲，而是检查坐标系。在 ArcGIS Pro 中，右键图层 → Properties → Source，查看 Spatial Reference。在 ArcMap 中，右键 → Properties → Source。也可以用 Python 快速批量检查：

import arcpy

fc = r"C:\data\schools.shp"
desc = arcpy.Describe(fc)
sr = desc.spatialReference
print("名称:", sr.name)
print("类型:", sr.type)          # Geographic 或 Projected
print("线性单位:", sr.linearUnitName)
print("角度单位:", sr.angularUnitName)
print("WKID:", sr.factoryCode)

如果输出 type 为 Geographic，说明数据是经纬度，不能直接做米制缓冲。如果 linearUnitName 为空而 angularUnitName 为 Degree，更要警惕。

国内数据还有一个特殊坑：高德地图使用 GCJ-02 坐标系，百度使用 BD-09，与 WGS84/CGCS2000 存在几十到上百米的偏移。如果学校点位来自高德 POI，而底图是 CGCS2000，直接缓冲会产生系统性位置偏差。需要先做坐标纠偏（GCJ-02 → WGS84），再进入分析流程。相关纠偏算法可参考开源项目 coordTransform_py。

3.2 选择投影坐标系

选择投影的原则是：研究区尽量靠近中央经线，投影变形尽量小，单位是米。对单个城市尺度的学校服务区分析，推荐：

• 国内数据：CGCS2000 3-degree Gauss-Kruger CM 117E（EPSG:4547）或对应带号，视城市经度而定。
• 国际数据：WGS84 UTM Zone 50N（EPSG:32650，覆盖东经 114°–120°）等。
• 如果研究区跨带，可考虑 Albers 等积投影自定义中央经线，或使用本地工程坐标系。

ArcGIS Pro 中投影操作：Analysis → Tools → Data Management → Projections and Transformations → Project。输入要素、输出坐标系选上述 EPSG，地理变换（Geographic Transformation）按区域选择，如 CGCS2000 与 WGS84 之间可选 WGS_1984_(ITRF00)_To_CGCS2000 等。

本文评述：投影不是“一次性设置”，而是分析流程中的关键决策点。建议把投影步骤显式写进工作流文档，而不是依赖软件默认。

3.3 生成 500 米缓冲区

投影完成后，使用 Analysis Tools → Proximity → Buffer（ArcGIS Pro 中为 Pairwise Buffer 或 Buffer）。参数设置：

Python 调用示例：

import arcpy

arcpy.analysis.Buffer(
    in_features=r"C:\data\schools_projected.shp",
    out_feature_class=r"C:\data\schools_500m.shp",
    buffer_distance_or_field="500 Meters",
    line_side="FULL",
    line_end_type="ROUND",
    dissolve_option="ALL"
)

注意 "500 Meters" 这个写法：显式带单位可以避免软件按坐标系默认单位解释，是工程上更稳妥的习惯。

3.4 结果校核

缓冲生成后，不能直接出图，至少做三项校核：

1. 量测校核：用 Measure 工具量取缓冲区边缘到学校点的距离，确认约等于 500 米。
2. 叠置校核：与行政区划、路网叠置，检查是否有明显不合理的跨界或穿越障碍。
3. 面积校核：单个 500 米圆的理论面积约 π×500² ≈ 785,398 平方米。若缓冲区面积远大于或远小于此值，需排查坐标系或单位问题。

本文评述：校核这一步最容易被跳过，但它是把“操作”变成“可信结果”的分水岭。建议把面积校核写成自动化脚本，作为每次缓冲分析的固定检查项。

3.5 可视化与出图

出图时，缓冲区图层建议用半透明填充（如紫色 #7c3aed，透明度 40%），学校点用醒目符号，底图用浅色。若需展示重叠情况，可用 Dissolve 后的合并面与未合并面叠加对比。

Esri 官方提供了多个缓冲区分析教程，例如 ArcGIS Learn: Analyze Crime Hot Spots 和 Buffer 工具官方文档。视频教程可参考 Esri 官方 YouTube 频道的 “How to create buffers in ArcGIS Pro”。

四、距离偏差的六大来源与诊断方法

4.1 来源一：坐标系类型误用

最典型的是在 GCS 数据上直接做米制缓冲。诊断方法：检查图层 Spatial Reference 的 type 字段。若为 Geographic，立即停止缓冲，先投影。

4.2 来源二：投影变形过大

研究区远离中央经线，或使用了 Web Mercator 等不适合距离分析的投影。诊断方法：计算研究区到中央经线的经度差，估算长度变形。对高斯-克吕格 3° 带，经度差超过 1.5° 就应考虑换带或自定义中央经线。

4.3 来源三：单位设置错误

距离填了 500 但单位选了 Kilometers，或软件默认单位与预期不符。诊断方法：显式指定单位，如 "500 Meters"。

4.4 来源四：数据偏移（GCJ-02/BD-09）

国内 POI 数据未纠偏，导致学校点位整体偏移。诊断方法：与已知控制点或高精度底图比对，偏移量通常在 50–500 米。处理方法：使用坐标纠偏算法转换到 WGS84/CGCS2000。

4.5 来源五：几何错误

学校点存在重复、空几何、无效几何，导致缓冲结果异常。诊断方法：使用 Check Geometry 和 Repair Geometry 工具。

4.6 来源六：模型选择错误

用欧氏缓冲代替网络服务区，忽略了路网阻隔。诊断方法：对比欧氏缓冲与网络服务区结果，若差异显著（如隔河、隔铁路），说明需要网络分析。

本文评述：这六类来源中，前两类是坐标系问题，后四类分别涉及数据、参数、几何与模型。把这张表当作缓冲区分析的“故障排查清单”，能覆盖绝大多数入门到中级场景的错误。

五、进阶：网络服务区、多环缓冲与三维缓冲

5.1 网络服务区（Service Area）

欧氏缓冲的局限在于忽略路网。ArcGIS Network Analyst 的 Service Area 工具可以基于道路网络计算 500 米步行或车行可达范围。核心步骤：构建网络数据集 → 新建 Service Area 图层 → 设置设施点（学校）→ 设置中断值（500 米）→ 求解。

网络服务区的结果通常是不规则多边形，受路网密度、单行道、过街天桥、河流阻隔影响。国内外大量研究用网络服务区替代欧氏缓冲评估学校可达性（如 Talen & Anselin, 1998 的经典可达性框架；近年国内学者在城市规划领域也有大量应用）。本文评述：网络服务区更接近真实步行体验，但对路网数据质量要求高，且计算成本显著上升。入门阶段先用欧氏缓冲建立直觉，进阶再上网络分析，是合理的学习路径。

5.2 多环缓冲（Multiple Ring Buffer）

实际分析中常需要 200 米、500 米、1000 米多级服务范围。ArcGIS 提供 Multiple Ring Buffer 工具，可一次生成多个距离环。参数中需注意 Dissolve Option 的选择：ALL 会合并所有环，NONE 保留独立环，LIST 按字段合并。

5.3 三维缓冲与地形影响

在山区，平面 500 米缓冲与真实地面 500 米步行距离差异巨大。ArcGIS 3D Analyst 提供 Buffer 3D 工具，可基于 TIN 或地形表面生成三维缓冲。相关研究指出，坡度对步行速度的影响可用 Tobler 步行函数等模型量化（Tobler, 1993）。本文评述：三维缓冲是欧氏缓冲在地形维度上的自然延伸，但计算复杂度和数据要求更高，适合山地城市或精细化通学分析。

5.4 椭球距离与大地线

对于跨区域、大范围分析，平面投影变形不可忽略，此时可用椭球距离（Geodesic Distance）工具，直接在椭球面上计算距离。ArcGIS Pro 的 Buffer 工具在地理坐标系下也支持 GEODESIC 方法。本文评述：椭球距离适合大尺度、高精度场景，但对城市级 500 米服务区，投影坐标系仍是更简单可靠的选择。

六、工程实践中的坐标系治理清单

6.1 数据入库阶段

• 统一元数据：记录每个图层的坐标系、来源、精度。
• 统一基准：国内项目建议统一到 CGCS2000。
• 统一投影：按研究区选择标准分带投影，避免混用。

6.2 分析流程阶段

• 显式投影：不依赖软件默认，每一步都明确坐标系。
• 显式单位：距离参数带单位。
• 显式校核：面积、量测、叠置三重检查。

6.3 成果输出阶段

• 成果图层附带坐标系说明。
• 出图时标注投影与距离含义（直线/网络）。
• 保留中间数据，便于复现与审计。

本文评述：坐标系治理不是一次性任务，而是贯穿数据全生命周期的纪律。把“距离单位依赖坐标系”这句话贴在工位上，比记住任何工具按钮都重要。

七、前沿与预判：从欧氏缓冲到可达性建模

7.1 从几何缓冲到行为建模

近年学校可达性研究正从“几何缓冲”走向“行为建模”。例如，基于手机信令、公交刷卡、共享单车轨迹的实证研究，可以反推真实通学距离分布。这类研究对“500 米”这一政策阈值提出了实证检验：不同城市、不同年龄段学生的实际步行耐受距离差异显著。

国际期刊如 Journal of Transport Geography、Computers, Environment and Urban Systems 近三年持续发表相关研究。本文评述：几何缓冲是规范性的（normative），行为数据是实证性的（positive），两者结合才能支撑科学决策。

7.2 2SFCA 与供需视角

两步移动搜索法（2SFCA）及其改进形式（如 Gaussian 2SFCA、3SFCA）把服务范围从“有没有覆盖”推进到“供需匹配程度”。学校学位供给与学龄人口需求的比值，比单纯的 500 米覆盖更能反映教育资源配置的公平性。相关方法在 Health & Place、Annals of the AAG 等期刊有大量讨论。

7.3 实时与动态服务范围

随着实时交通数据、天气数据、事件数据的可获取性提升，静态 500 米缓冲正在向动态服务范围演进。例如，考虑早晚高峰步行速度变化、雨天绕行、施工封路等因素，服务范围会随时间波动。ArcGIS Velocity、GeoEvent Server 等实时 GIS 平台为这类分析提供了基础设施。

本文评述：未来的服务范围分析，坐标系仍然是底座，但“距离”将从静态几何量变成动态时空量。掌握坐标系原理，是理解这一切变化的前提。

八、参考文献与数据来源