从Excel到地图：手把手教你用ArcGIS 10.2处理气象数据，搞定气温分布图
从Excel到地图：手把手教你用ArcGIS 10.2处理气象数据，搞定气温分布图
在环境监测、农业规划和气候研究中，气温空间分布图是最基础也最重要的可视化工具之一。许多初学者面对杂乱的气象站原始数据时，往往卡在数据预处理阶段——那些躺在Excel里的数字如何变成专业地图上的等温线？本文将带你完整走通这条数据处理流水线：从Excel表格清洗到ArcGIS空间分析，最终输出可直接用于科研报告的气温分布图。
1. 气象数据的前期处理：Excel中的关键四步
气象数据通常以TXT或CSV格式从观测站导出，直接导入ArcGIS会导致字段混乱。我们需要先在Excel中完成数据标准化。以下是实测有效的操作流程：
1.1 文本分列：结构化原始数据
打开从中国气象数据网下载的TXT文件时，所有信息可能挤在一列中。使用Excel的文本分列向导可将其拆分为多列：
- 选中数据列 → 点击"数据"选项卡 → 选择"分列"
- 在向导中选择"分隔符号" → 勾选"空格"或"制表符"（根据实际数据格式）
- 预览分列效果 → 完成
注意：若数据包含非标准分隔符，建议先用记事本查看原始格式
1.2 数据关联：VLOOKUP实战
气象数据通常分散在多个文件，比如站点气温和站点坐标分开存储。使用VLOOKUP函数可以快速关联：
=VLOOKUP(查找值, 表格区域, 返回列序号, FALSE)具体操作示例：
- 在气温表新增"经度"、"纬度"列
- 输入公式：
=VLOOKUP(A2, 坐标表!A:C, 2, FALSE)获取经度 - 同理获取纬度值
提示：FALSE参数确保精确匹配，避免数据错位
1.3 数据质量检查
在导出到ArcGIS前，必须检查：
- 所有站点是否有完整坐标
- 气温值是否存在异常（如<-50℃或>50℃）
- 日期格式是否统一
推荐使用Excel的条件格式快速定位异常值：
- 选中气温列 → 开始 → 条件格式 → 色阶
- 观察极端值颜色标记
1.4 导出为CSV
最终保存时选择"CSV UTF-8"格式，确保中文字符正常显示。这是ArcGIS最兼容的表格格式。
2. ArcGIS中的数据导入与坐标转换
2.1 创建空间数据
在ArcMap 10.2中导入CSV文件后：
- 右键表格 → 显示XY数据
- X字段选经度，Y字段选纬度
- 坐标系选择WGS1984（EPSG:4326）
此时生成的是临时点图层，需右键导出为正式Shapefile：
- 输出位置选择专用工作文件夹
- 命名规范建议：
Hubei_Temp_Stations.shp
2.2 投影转换实战
省级分析通常需要投影坐标系。以湖北省为例，转换到兰伯特投影：
| 步骤 | 工具位置 | 关键参数 |
|---|---|---|
| 1 | 工具箱 → 数据管理 → 投影和变换 → 要素 → 投影 | 输出坐标系：Asia_Lambert_Conformal_Conic |
| 2 | 地理处理 → 环境 → 处理范围 | 设置为"湖北省边界图层" |
# ArcPy脚本实现（供高级用户参考） arcpy.Project_management("Hubei_Temp_Stations", "Hubei_Temp_Stations_Lambert", "PROJCS['Asia_Lambert_Conformal_Conic']")重要：始终检查转换后的点位是否仍落在省界内，防止投影参数错误
3. 气温空间插值：克里金法深度配置
3.1 插值方法选型
在Geostatistical Analyst工具栏中，主要方法比较：
| 方法 | 适用场景 | 本例选择 |
|---|---|---|
| 反距离权重 | 数据密集且均匀 | × |
| 克里金 | 考虑空间自相关性 | √ |
| 协同克里金 | 有辅助变量（如海拔） | 可选 |
克里金法参数设置技巧：
- 趋势移除：选择"一次"多项式（气温通常随纬度线性变化）
- 半变异函数：指数模型更适合气温连续变化
- 步长数：建议设为经度范围/2（湖北约24步）
- 扇区类型：8扇区提高各向异性分析
3.2 插值结果优化
生成预览后，需调整三类参数：
- 统计检验：查看交叉验证的RMSE（应<1℃）
- 可视化分级：右键图层属性 → 符号系统 → 分类
- 方法：自然间断点（Jenks）
- 类别数：7-9级为宜
- 掩膜提取：使用省界裁剪多余区域
arcpy.gp.ExtractByMask_sa("temp_surface", "hubei_boundary", "final_map")
4. 专业地图制作与输出
4.1 地图元素排版
切换到布局视图，添加：
- 比例尺（公里/英里双单位）
- 指北针（简洁箭头样式）
- 图例（调整色带文字说明）
- 数据来源声明（必填）
进阶技巧：
- 插入动态文本（如"制图日期：<%d-%m-%Y%>"）
- 使用格网标注经纬度
4.2 多格式输出
根据使用场景选择输出格式：
| 格式 | 用途 | DPI建议 |
|---|---|---|
| 印刷出版 | 300+ | |
| PNG | 网页展示 | 150 |
| SVG | 矢量编辑 | - |
| AI | 进一步美化 | - |
最后保存地图文档（.mxd）时，使用相对路径存储，方便项目迁移。
常见问题解决方案
Q1 插值结果出现异常高温圈？
- 检查原始数据是否有录入错误
- 尝试调整半变异函数范围（Range参数）
Q2 部分站点未显示在地图？
- 确认坐标是否为十进制度数（非度分秒）
- 检查数据框坐标系是否与图层一致
Q3 图例色带不连续？
- 右键图层 → 属性 → 符号系统
- 取消勾选"仅显示分类范围内的值"
实际操作中，建议保存关键步骤的中间文件（如清洗后的Excel、投影前的Shapefile等），方便回溯修改。记得定期使用ArcGIS的"文档压缩"功能清理冗余数据，保持项目文件轻量化。