利用Arcgis提取nc数据并计算区域年均温

因自己不间断需要用到nc数据，故写下这篇笔记，本篇包括过程的步骤和注意事项，注释详细，希望同时与更多人交流学习。

---

1.从国家青藏高原科学数据中心（https://data.tpdc.ac.cn/）获取得到中国1km分辨率逐月平均气温数据集（1901-2022）。

2.解压数据至不含中文的路径中。打开Arcmap，选择工具箱中的多维分析工具-创建NetCDF栅格图层，导入数据。（本例子导入数据为2021年nc数据：tmp_2021.nc）

注：导入数据后若系统未自动读取变量等参数，表明数据所在路径存在中文，需修改数据保存路径，后重新导入。

3.打开工具箱中的数据管理工具-图层和表视图-创建栅格图层，导入数据，对输出数据进行命名，得到月均温栅格数据。（本例子利用数据内容temperature和对应时间，将二者组合得到命名tmp202101，并选择对应波段1）

⌈nc数据可能为一年数据，也可能为多年组合数据，当nc数据为一年数据时，波段可选数字为1~12，表示当年的1月至12月；若nc数据为多年组合数据，如2006~2008年（3年），则波段可选数字为1~36，1表示该组合数据的第一个月，即2006年1月，2表示该组合数据的第一个月，即2006年2月，依此类推，13表示该组合数据的第13个月，即2007年1月。⌋

按顺序命名并选择对应波段，确认即可创建。

TIP：当波段较多时，也可以右键创建栅格图层，选择批处理。统一选择输入栅格，命名输出栅格，并选择对应波段。

4.得到每月对应的月均温数据后，选择工具箱-Spatial Analyst工具（空间分析工具）-地图代数-栅格计算器，计算全国一年年均温。（本例子得到2021年全国年均温）

计算式放在这里供大家参考（所有符号均为英文符号）：("tmp202101"+"tmp202102"+"tmp202103"+"tmp202104"+"tmp202105"++"tmp202106"+"tmp202107"++"tmp202108"+"tmp202109"++"tmp202110"+"tmp202111"+"tmp202112")*0.1/12

5.得到一年年均温后可将其导出为.tif格式数据，方便后期使用。（本例子命名与输出栅格保持一致）

至此，由nc数据提取并计算得到全国年均温数据的步骤已全部完成。

6.选择工具箱-Spatial Analyst工具（空间分析工具）-区域分析-以表格显示分区统计，得到分区域的气温数据。（本例子选择的矢量边界为省级边界，表示统计对象为矢量数据中的各省。统计类型选MEAN或查看ALL中的MEAN即可计算矢量数据中各省的年均温。）

注：区域字段应选择数字类字段，字符串型字段（中文）会出现乱码。

统计类型可点击显示帮助查看，如下图。

分区域年均温结果如下图所示：

---

主要参考资料：

①利用月均温数据计算年均温/NC文件如何用ArcGIS打开？

②ArcGIS软件对降水量数据与气温数据的nc格式文件后续处理 

③Arcmap转nc文件为TIFF格式

如有错误敬请指正，谢谢。