前不久给大家分享了使用 R 语言爬取 GSOD 气象站点数据的方法：
×
并且给大家分享了 GSOD 气象站点数据：
×
最近我仿照之前处理气象数据的方法处理了 平均观测站气压与平均海平面气压 数据，也就是首先把气象站点的数据插值成 0.1˚x0.1˚ 分辨率的网格数据，然后把网格数据转换成栅格数据再分区域平均得到各省市区县的年度、月度和日度数据，具体处理细节如下：
- 因为之前的数据是使用 2020 年行政区划进行插值处理得到，为了方便对比，这次依然使用了 2020 年的行政区划
- 使用的网格是 0.1˚x0.1˚ 分辨率的网格，也就是每个格点长宽都是 0.1˚（纬度/经度），更加精细；
- 先把插值得到的网格数据构建成栅格数据数据，然后再分区域汇总；
- 从栅格数据分区域汇总平均即可得到省市区县的面板数据；
- 为了避免过度平滑的问题，进行插值的时候仅仅使用每个格点附近的 5 个站点数据；
- 考虑到过早年份数据中的气象站点过少，仅仅处理了 1973 年至今的数据。
以 2023 年 4 月 18 日平均观测站气压的处理为例，首先是该天所有气象站点的数据：
使用 IDW 法插值成 0.1˚x0.1˚ 分辨率的网格数据：
把插值得到的网格数据转换成栅格数据（这里展示的是 2022 年年度平均的）：
栅格数据分区域汇总就可以得到省市区县的数据了，例如区县的：
按照上面的处理思路循环处理 1973年1月1日～2023年4月18日的数据即可各省市区县分年、分月、逐日的平均观测站气压和海平面气压了，为了方便大家的使用，我把数据保存成了 xslx 格式的：
每份数据都包含各级行政区划代码，方便大家和其他数据集进行匹配，作为参考，下图展示了 2021 年 8 月 1 日各城市平均观测站气压数据：
2022 年 6 月各省份平均观测站气压：
由于这次处理过程中生成了栅格数据，所以附件中也有每天的栅格数据。
最后的这个图展示了 1973～2022 年中国各省市年平均观测站气压的变化（2023 年的数据不足全年）：
下图展示了这次处理的结果与之前版本的结果对比：
看起来不太一致，可能是数据的口径不同，无法比较。使用的时候大家也可以看看哪个好用用哪个。
之前版本的数据在这里（如果只需要 2020 年之前的数据更建议用这个）：
×
注意事项
- 由于压强可能受到地形、植被、水文、人类活动等因素的影响，因此使用上文所说的方法计算得到的各区域均值可能与实际值有所偏差；
- 虽然上文中的图只展示了平均观测站压强，实际上附件中的数据里还有平均海平面压强；
- 这里压强的单位是 hPa。
处理方法
感兴趣的小伙伴可以参考这个课程学习：
×
点击这里跳转到 RStata 短书平台获取附件：旧版本｜1973～2023年各省市区县平均观测站气压与海平面气压日度面板数据
评论