如何利用共享数据手搓一份历年降雨分月数据

"前段时间做高标准农田项目，其中涉及水资源平衡 。评审的专家要求我们收集的降雨资料要更新到2023年，在网上地方气象局怎么也找不到数据，你说去当地气象局申请吧，这事情得部门打报告、搞保密协议、找领导签字，项目要这么搞GG的快。下面，你看小编如何利用网上数据，手搓出这份历年降雨分月数据的。 "

资源下载

数据是大神分享在国家青藏高原科学数据中心： 《中国1km分辨率逐月降水量数据集（1901-2023）》 。无需注册，网站为：

https://data.tpdc.ac.cn/zh-hans/data/faae7605-a0f2-4d18-b28f-5cee413766a2

数据下载，因为数据量挺大，20GB左右吧，网站也给出了下载方式:FTP下载。

在这里我采用的是Filezilla下载，只需下载一个FileZilla Client到电脑上，本机跟上面的主机建立联系即可下载数据到本地（我挂电脑一晚上就下载好了） 。

在文件——站点管理器——新站点，将前面图上的主机、端号、用户名、密码填上，即可完成下载。

一个晚上：1901-2023年的中国逐月降水量数据就下载好了，文件格式为.nc格式。

数据格式转化

接下来就是将nc格式文件转换成我们常用的tiff格式，这个操作在Arcmap里完成。步骤：ArcToolbox——多维工具——创建NetCDF栅格图层。输入netCDF文件为下载的文件，以2023为例子。

波段维度选择时间，其他默认，确认。

接下来 ArcToolbox——数据管理工具——图层和表视图——创建栅格图层 。每个波段代表一个月份。

导出成tif数据备用：

03

提取栅格数据导入到Excel中

将栅格进行投影，因为行政边界数据投影坐标系是：

WGS_1984_World_Mercator。

这里将上面生成的tif定义成同样的投影。

下面就只裁剪你需要的行政区的tif。 ArcToolbox——数据管理工具——栅格——栅格处理——裁剪 。

提取数据， ArcToolbox——Spatial Analyst工具——区域分析 ——以表格显示分区统计 。（涉及到扩展模块，还需要扩展模块的 Spatial Analyst打上勾，在自定义工具栏下面）。

区域字段：选择是数字的字段FID（字符串会报错001143），输出表路径采用英文（中文会报错999999），统计类型：使用mean。

由于上个步骤用的是FID字段，还得在生成的表格中添加字段，然后连接shp文件的区县字段给区县赋值。

导出表格，在excel中添加表格数据，即可获得区县当月的平均降雨量。