地表温度日数据转每月数据

根据https://www.jianshu.com/writer#/notebooks/8716507/notes/82102474处理得到的每日地表温度数据其单位并不是℃,这个影像的数值需要*0.02-273.15得到单位为℃的地表温度图层(详见https://www.docin.com/p-1638868980.html手册),但是仔细查看下载下来的地表温度数据,你会发现有很多地方是0值,如果是0值转换成℃图层的话,那这些地区对应的就是-273.15℃,这肯定是不合理的,所以对于每日的地表温度tif数据,我们还需要先剔除异常值,通过将0值区赋值为空值,然后再采用focalstatistics工具进行处理,即根据领域像元的数值计算空值区数值,具体实现方式在arcgis中实现,实现流程:

先在arcgis中点击spatial analysis tools->map algebra->raster calculator

输入SetNull("MOD11A2.A2018001.LST_Day_1km.tif" == 0,"MOD11A2.A2018001.LST_Day_1km.tif"),将温度图层中为0的数据设置为nodata;

第二步,根据空值像元周围数据填充空值区像元地表温度,实现途径:

spatial analysis tools->neighborhood->focal statistic,计算以5单位的圆形像元计算领域平均值赋值给空值区域,


地表温度日数据转每月数据_第1张图片
默认勾选ignore nodata in calculation,计算范围设置为比研究区大的矩形区域

此过程大致运行了9次完成了空值区域的填充,其他日期的地表温度数据按照同样的步骤剔除异常值,然后2018年1月份的地表温度平均值就是2018年第1天,第9天,第17天和第25天温度数据求平均。

地表温度日数据转每月数据_第2张图片
2018年1月份地表温度

计算结果就如图所示,1月份地表平均温度最低有零下36℃,数据显示黑龙江地区最低温度可达到零下53℃,因而计算结果无误。

你可能感兴趣的:(地表温度日数据转每月数据)