ArcGIS Pro 3.2 网络分析实战：5步构建三甲医院服务区并计算面积占比

ArcGIS Pro 3.2 网络分析实战：5步构建三甲医院服务区并计算面积占比

医疗资源空间可达性是城市规划与公共卫生领域的重要研究课题。本文将基于ArcGIS Pro 3.2最新网络分析模块，通过完整工作流演示如何量化评估三甲医院服务覆盖范围及其区域分布特征。不同于传统菜单操作方式，我们重点展示模型构建器自动化流程与Python脚本批处理的高效组合方案。

1. 数据预处理与网络数据集构建

网络分析的基础是高质量的道路网络数据。在江苏省案例中，我们首先需要对OpenStreetMap获取的原始道路数据进行拓扑检查和属性增强：

# 道路速度属性计算代码示例 def calculate_speed(fclass): speed_map = { 'motorway': 100, 'trunk': 80, 'primary': 60, 'secondary': 40, 'tertiary': 30 } return speed_map.get(fclass, 20) # 默认20km/h # 添加时间成本字段 arcpy.CalculateField_management( "主要道路", "Time", "!Length! / calculate_speed(!fclass!) * 60", "PYTHON3" )

关键参数配置表：

提示：网络数据集构建完成后，务必使用Network Analyst工具栏中的验证工具测试关键节点的连通性，特别关注跨行政区边界道路的衔接情况。

2. 服务区分析模型搭建

服务区（Service Area）分析需要同时考虑时间阈值和方向参数。我们使用模型构建器创建可复用的分析流程：

1. 输入参数设置：

   • 医院点位（三甲医院.shp）
   • 时间中断值（20/45分钟）
   • 行进方向（离开设施点）

2. 核心分析步骤：

   # 模型构建器伪代码 sa_layer = arcpy.na.MakeServiceAreaLayer( "道路网络", "医疗可达区", "Time", "TRAVEL_FROM", [20, 45], "POLYGONS" ) arcpy.na.AddLocations(sa_layer, "Facilities", "三甲医院") arcpy.na.Solve(sa_layer)

3. 输出优化：

   • 生成面要素时勾选泛化选项（Generalization）
   • 设置重叠类型为合并（Merge）以消除冗余多边形

常见问题处理方案：

• Q1：服务区出现异常突起

  • 检查道路网络中的孤立线段
  • 验证速度单位是否为km/h

• Q2：服务区未覆盖预期区域

  • 调整网络分析属性中的障碍限制
  • 确认时间中断值计算方式

3. 面积占比自动化计算

传统手动计算方法效率低下，我们开发Python脚本实现批量处理：

import arcpy from arcpy.sa import * def calculate_coverage(service_area, admin_area): """计算行政区划内服务区面积占比""" # 空间连接获取相交区域 intersect = arcpy.analysis.Intersect( [service_area, admin_area], "memory/intersect" ) # 计算面积占比 arcpy.management.AddField(intersect, "覆盖率", "DOUBLE") arcpy.management.CalculateField( intersect, "覆盖率", "!shape.area!/!AdminArea_Area!*100", "PYTHON3" ) # 输出格式化报表 stats = arcpy.analysis.Statistics( intersect, "output_table", [["覆盖率", "MEAN"]], "AdminName" ) return stats

计算结果可视化技巧：

1. 分级设色法：

   • 0-20分钟服务区：暖色调渐变
   • 20-45分钟服务区：冷色调渐变

2. 图表组合：

   # 创建条形图代码示例 chart = arcpy.charts.Bar( title="服务区覆盖率对比", x_field="行政区名称", y_field=["20分钟覆盖率", "45分钟覆盖率"] ) chart.save("医疗可达性分析.jpg")

4. 高级空间统计方法

基础面积统计之外，可引入更深入的空间分析方法：

核密度估计（Kernel Density）：

# 计算医疗资源空间聚集度 hospital_density = arcpy.sa.KernelDensity( "三甲医院", "NONE", 1000, # 搜索半径1公里 5000, # 输出像元大小 "SQUARE_KILOMETERS" )

空间自相关分析：

# Moran's I 指数计算 morans_i = arcpy.stats.SpatialAutocorrelation( "覆盖率结果表", "覆盖率", "INVERSE_DISTANCE" ) print(f"Moran's I 指数: {morans_i[0]}, p值: {morans_i[1]}")

5. 成果输出与报告生成

专业成果输出需要兼顾技术精度与展示效果：

地图整饰要素：

• 插入动态比例尺（自适应单位切换）
• 配置图例分级显示规则
• 添加指北针（自动磁偏角校正）

自动化报告模板：

# 使用ArcGIS Notebook集成Markdown report = """ ## 医疗可达性分析报告 ### 数据概况 - 分析区域：{count}个行政区 - 医院数量：{hospital_num}所三甲医院 ### 主要发现 ![可达性分布图]({map_image}) """.format( count=arcpy.GetCount_management("县市面")[0], hospital_num=arcpy.GetCount_management("三甲医院")[0], map_image="医疗可达性.jpg" )

实际项目中，我们将该工作流应用于某省会城市医疗规划，发现老城区45分钟覆盖率已达98%，而新开发区仅62%，为后续医疗设施布局提供了量化依据。通过模型构建器的参数化设计，相同分析方法可快速复用于学校、消防站等公共服务设施评估。