ArcGIS实战：手把手教你绘制土壤重金属污染分布图（以贵阳Cd镉为例）
ArcGIS实战：从采样数据到污染热图的完整工作流
清晨六点，实验室的离心机还在嗡嗡作响。这是环境监测站的王工连续第三天加班处理土壤样本数据了。他面前的Excel表格里记录着贵阳周边200多个采样点的镉含量数据，但如何将这些数字转化为直观的空间分布图？这个问题困扰着许多像王工这样的基层环保工作者。本文将用厨房烹饪般的步骤分解，带你掌握从原始数据到专业污染分布图的完整ArcGIS工作流。
1. 数据准备与预处理
任何GIS分析都始于干净规整的数据。打开你的Excel文件时，请确认至少包含三列关键数据：采样点编号、经纬度坐标（或平面坐标）、重金属浓度值。常见的初学者错误是直接使用"东经113°26'"这样的文本格式坐标，这会导致GIS无法识别。
数据清洗 checklist：
- 删除所有合并单元格和空行
- 浓度单位统一为mg/kg（避免部分数据用ppm）
- 经度范围验证（贵阳地区约在106-107°E之间）
- 异常值标注（如浓度超过1000mg/kg需核实）
# 示例：用Python快速检查数据范围 import pandas as pd df = pd.read_excel('soil_samples.xlsx') print(f"经度范围: {df['Longitude'].min()}~{df['Longitude'].max()}") print(f"镉浓度范围: {df['Cd'].min()}~{df['Cd'].max()}mg/kg")提示：建议保留原始数据副本，所有操作在复本上进行。GIS中的很多操作不可逆。
2. 空间数据导入与坐标系配置
在ArcGIS Pro中创建新项目后，通过"地图"选项卡→"添加数据"导入Excel文件。此时数据还只是属性表，需要右键选择"显示XY数据"，将经度、纬度字段分别指定为X、Y字段。
常见坐标系选择参考表：
| 适用范围 | 推荐坐标系 | EPSG代码 |
|---|---|---|
| 贵阳市区 | CGCS2000 / 3-degree Gauss-Kruger zone 35 | 4547 |
| 贵阳周边区县 | WGS 1984 UTM Zone 48N | 32648 |
| 全省范围 | China Geodetic Coordinate System 2000 | 4490 |
坐标系配置不当会导致后续插值结果偏移数公里。我曾见过某环保局报告因错选坐标系，整个污染区域向南平移了800米。完成坐标定义后，务必通过"导出数据"将点要素永久保存为shapefile或地理数据库要素类。
3. 空间插值方法深度解析
ArcGIS提供12种空间插值方法，针对重金属污染分析，这些方法各有优劣：
克里金法（Kriging）：
- 优势：考虑空间自相关，可生成预测标准差表面
- 参数要点：
- 半变异函数选择：球形(Spherical)适合局部变异
- 搜索半径：贵阳城区建议设置500-1000米
- 邻域点数：12-15个采样点为佳
# 半变异函数模型选择参考 variogram_models = { 'Spherical': '适合有明显变程的情况', 'Exponential': '适合渐变型空间变异', 'Gaussian': '适合非常平滑的连续表面' }反距离权重法（IDW）虽然简单，但容易在采样稀疏区产生"牛眼"效应。某次分析中，使用IDW导致一个采样点周围出现夸张的圆形污染晕，而实际调查发现该点是实验室标样误放入实测样本所致。
4. 专业制图与可视化技巧
得到插值结果后，右击图层选择"符号系统"，这才是展现专业性的关键阶段。不建议直接使用默认的"拉伸"渲染，重金属污染图更适用"分类"方法。
镉污染等级划分建议（根据GB 15618-2018）：
| 污染等级 | 镉含量(mg/kg) | 推荐色标 |
|---|---|---|
| 清洁区 | <0.3 | 深绿色 |
| 警戒线 | 0.3-0.6 | 浅绿色 |
| 轻度污染 | 0.6-1.0 | 黄色 |
| 中度污染 | 1.0-2.0 | 橙色 |
| 重度污染 | ≥2.0 | 红色 |
在布局视图中添加必备地图元素时，注意：
- 比例尺用线段式而非数字式
- 指北针避免使用复杂样式
- 图例需注明"基于克里金插值结果"
- 添加数据来源和采样时间说明
最后导出图片时，300dpi的TIFF格式适合报告印刷，而PNG格式更适合电子文档。记得在图层属性→元数据中填写完整的作者和版权信息。