从Excel到地图:手把手教你用ArcGIS 10.2处理气象数据,搞定气温分布图

从Excel到地图:手把手教你用ArcGIS 10.2处理气象数据,搞定气温分布图

从Excel到地图:手把手教你用ArcGIS 10.2处理气象数据,搞定气温分布图

在环境监测、农业规划和气候研究中,气温空间分布图是最基础也最重要的可视化工具之一。许多初学者面对杂乱的气象站原始数据时,往往卡在数据预处理阶段——那些躺在Excel里的数字如何变成专业地图上的等温线?本文将带你完整走通这条数据处理流水线:从Excel表格清洗到ArcGIS空间分析,最终输出可直接用于科研报告的气温分布图。

1. 气象数据的前期处理:Excel中的关键四步

气象数据通常以TXT或CSV格式从观测站导出,直接导入ArcGIS会导致字段混乱。我们需要先在Excel中完成数据标准化。以下是实测有效的操作流程:

1.1 文本分列:结构化原始数据

打开从中国气象数据网下载的TXT文件时,所有信息可能挤在一列中。使用Excel的文本分列向导可将其拆分为多列:

  1. 选中数据列 → 点击"数据"选项卡 → 选择"分列"
  2. 在向导中选择"分隔符号" → 勾选"空格"或"制表符"(根据实际数据格式)
  3. 预览分列效果 → 完成

注意:若数据包含非标准分隔符,建议先用记事本查看原始格式

1.2 数据关联:VLOOKUP实战

气象数据通常分散在多个文件,比如站点气温和站点坐标分开存储。使用VLOOKUP函数可以快速关联:

=VLOOKUP(查找值, 表格区域, 返回列序号, FALSE)

具体操作示例:

  • 在气温表新增"经度"、"纬度"列
  • 输入公式:=VLOOKUP(A2, 坐标表!A:C, 2, FALSE)获取经度
  • 同理获取纬度值

提示:FALSE参数确保精确匹配,避免数据错位

1.3 数据质量检查

在导出到ArcGIS前,必须检查:

  • 所有站点是否有完整坐标
  • 气温值是否存在异常(如<-50℃或>50℃)
  • 日期格式是否统一

推荐使用Excel的条件格式快速定位异常值:

  1. 选中气温列 → 开始 → 条件格式 → 色阶
  2. 观察极端值颜色标记

1.4 导出为CSV

最终保存时选择"CSV UTF-8"格式,确保中文字符正常显示。这是ArcGIS最兼容的表格格式。

2. ArcGIS中的数据导入与坐标转换

2.1 创建空间数据

在ArcMap 10.2中导入CSV文件后:

  1. 右键表格 → 显示XY数据
  2. X字段选经度,Y字段选纬度
  3. 坐标系选择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 插值结果优化

生成预览后,需调整三类参数:

  1. 统计检验:查看交叉验证的RMSE(应<1℃)
  2. 可视化分级:右键图层属性 → 符号系统 → 分类
    • 方法:自然间断点(Jenks)
    • 类别数:7-9级为宜
  3. 掩膜提取:使用省界裁剪多余区域
    arcpy.gp.ExtractByMask_sa("temp_surface", "hubei_boundary", "final_map")

4. 专业地图制作与输出

4.1 地图元素排版

切换到布局视图,添加:

  • 比例尺(公里/英里双单位)
  • 指北针(简洁箭头样式)
  • 图例(调整色带文字说明)
  • 数据来源声明(必填)

进阶技巧

  • 插入动态文本(如"制图日期:<%d-%m-%Y%>")
  • 使用格网标注经纬度

4.2 多格式输出

根据使用场景选择输出格式:

格式用途DPI建议
PDF印刷出版300+
PNG网页展示150
SVG矢量编辑-
AI进一步美化-

最后保存地图文档(.mxd)时,使用相对路径存储,方便项目迁移。

常见问题解决方案

Q1 插值结果出现异常高温圈?

  • 检查原始数据是否有录入错误
  • 尝试调整半变异函数范围(Range参数)

Q2 部分站点未显示在地图?

  • 确认坐标是否为十进制度数(非度分秒)
  • 检查数据框坐标系是否与图层一致

Q3 图例色带不连续?

  • 右键图层 → 属性 → 符号系统
  • 取消勾选"仅显示分类范围内的值"

实际操作中,建议保存关键步骤的中间文件(如清洗后的Excel、投影前的Shapefile等),方便回溯修改。记得定期使用ArcGIS的"文档压缩"功能清理冗余数据,保持项目文件轻量化。