4. 热力图制作工具与环境搭建

热力图制作并不是单一软件能包打天下的任务。从原始监测数据的清洗、空间插值计算，到最终地图渲染与交互发布，不同环节适合不同工具。本节围绕三类主流工具——Python 数据科学生态、桌面 GIS 软件、WebGIS 前端库——说明它们各自的能力边界，并给出可落地的开发环境配置方案，帮助你在动手制作第一张热力图之前把"工具箱"准备好。

三类工具的分工与选型

理解工具分工，关键是看它在热力图生产链条中承担哪一段。热力图制作的典型链条是：数据准备 → 空间插值与核密度计算 → 栅格/矢量结果生成 → 地图渲染与交互 → 发布共享。

• Python 生态：擅长数据清洗、批量计算与自动化。常用组合为 pandas（表格处理）、geopandas（矢量空间数据）、numpy/scipy（数值计算）、pykrige（克里金插值）、matplotlib 或 folium（快速出图）。适合监测站点多、需要重复出图或做参数对比的场景。
• 桌面 GIS 软件：QGIS（开源）与 ArcGIS Pro（商业）提供成熟的插值工具、符号化系统和制图排版。适合需要精细制图、投影管理、成果出图交付的场景，交互式调参直观。
• WebGIS 库：Leaflet、OpenLayers、Mapbox GL JS、deck.gl 等，负责在浏览器中渲染可缩放、可交互的热力图。适合成果发布到网页、供多人查看的场景。

下图展示了三类工具在热力图生产链条中的位置与数据衔接关系，可作为环境搭建时的整体参照。

Python 环境配置实操

环境搭建最容易踩坑的是依赖冲突。推荐用 conda 创建独立环境，避免与系统 Python 混用。以下命令以 Windows/Linux 通用写法给出，示例版本号可按需调整：

# 1. 创建独立环境（Python 3.11）
conda create -n heatmap python=3.11 -y
conda activate heatmap

# 2. 安装核心依赖（conda 优先，地理库更稳）
conda install -c conda-forge geopandas rasterio pykrige -y
pip install pandas numpy matplotlib folium

# 3. 验证安装
python -c "import geopandas, pykrige, folium; print('环境就绪')"

安装完成后，用一段最小示例验证整条链路是否通畅：读取站点数据、做反距离权重（IDW）插值、输出一张预览图。示例数据为假设的 5 个监测点，仅用于验证环境，不代表真实监测结果。

import numpy as np, pandas as pd
from scipy.interpolate import griddata
import matplotlib.pyplot as plt

# 假设站点：x,y 为投影坐标(米)，pm25 为浓度(μg/m³)
df = pd.DataFrame({
    "x": [0, 1000, 2000, 500, 1500],
    "y": [0, 500, 0, 1500, 1200],
    "pm25": [35, 52, 41, 28, 60]
})

gx, gy = np.meshgrid(np.linspace(0, 2000, 120), np.linspace(0, 1500, 90))
gz = griddata((df.x, df.y), df.pm25, (gx, gy), method="cubic")

plt.contourf(gx, gy, gz, levels=15, cmap="YlOrRd")
plt.colorbar(label="PM2.5 (μg/m³)")
plt.scatter(df.x, df.y, c="k", s=18)
plt.title("IDW/插值环境验证示例")
plt.savefig("check.png", dpi=120)

若 check.png 正常生成且色带连续、站点位置合理，说明数值计算与绘图链路可用。注意 griddata 的 method 可选 nearest、linear、cubic，不同方法在边缘外推时表现差异明显，正式制作时应结合监测网络密度选择。

桌面 GIS 与 WebGIS 环境要点

QGIS 安装后需确认两点：一是"处理工具箱"中插值算法（如 IDW 插值、热力图 KDE）可用；二是工程坐标系设置正确——监测数据若为经纬度（EPSG:4326），插值前应投影到以米为单位的坐标系（如 CGCS2000 高斯投影或 UTM），否则距离权重会失真。ArcGIS Pro 对应使用"空间分析"扩展中的 IDW、克里金工具，需确保许可已激活。

WebGIS 环境相对轻量。以 Leaflet 为例，只需一个 HTML 文件引入 CDN 资源即可起步；若使用 leaflet.heat 插件渲染热力层，数据格式为 [纬度, 经度, 权重] 数组。生产环境建议用本地服务器（如 python -m http.server）打开页面，避免 file:// 协议下的跨域限制。Mapbox GL JS 与 deck.gl 需要申请访问令牌（token），配置时注意密钥不要硬编码进公开仓库。

环境自检清单