因自己不间断需要用到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格式文件后续处理
如有错误敬请指正,谢谢。