本数据集分析自NASA MERRA2卫星数据集,原始数据为全球42个气压层次下的气温,以6小时/次的频率采集,数据集代号M2I6NPANA。本数据集的分析时间范围是2000.1.1至2021.6.30,分析的粒度包括了每天4个时间点,每个时间点一个数据表。
数据格式见预览图片。其中第一列id,是地级市的编号;第二列name,代表城市中文名,根据2020年最新的行政区划地图进行逆温值归属;第三列date,代表分析的日期;第四列和第五列,分别代表地面层至第2层(g2),以及地面层至第3层(g3)之间的逆温强度。这里的g2是常规定义的逆温强度,g3作为检验稳健性的变量,可以按需使用。
接下来介绍具体的分析流程:
第一步,将栅格数据转化为逆温矩阵。由于中国不同地区的海拔差异较大,离地面最近的气压层,其气压不尽相同,表现为高海拔地区在高气压层出现缺失值。因此,我们将每个地区(栅格)内温度非缺失值,且气压最高的一层,定义为该栅格的地面层(g),往上的第二层、第三层气压层分别定义为L2和L3。
第二步,对栅格进行双线性插值。原始的栅格分辨率是0.5 ° x 0.625 °,我们将分辨率通过插值提升到0.1° x0.1°,随后在插值矩阵上进行后续分析。
第三步,计算逆温强度矩阵。在定义层之后,我们计算层温度差,L2-g和L3-g。如果值为正数,说明存在逆温,保持原值;如果值为负数,说明不存在逆温,对逆温强度无贡献,被重置为0。
第四步,根据行政区划分配逆温强度。得到逆温强度矩阵后,我们对每个地区的逆温强度进行平均运算,从得到一个地区在特定时间粒度下的当日逆温强度。地图采用了2020年的行政区划,分析至地市级。地图文件附在数据集链接中,可以根据需要下载。
【私信可优惠】
本数据集可以分析至县级水平。如果需要任何定制分析,欢迎联系。
邮箱:comtrader@foxmail.com