—种基于栅格数据快速识别极端高/低温事件的方法与流程

—
种基于栅格数据快速识别极端高/低温事件的方法
技术领域
1.本发明涉及极端气候识别技术领域，具体是指—种基于栅格数据快速识别极端高/低温事件的方法。


背景技术：

2.现有的识别极端低温事件的算法，大都基于站点数据，且处理步骤繁琐，运行时间久。鉴于此，本方案提出一种基于栅格数据的极端低温事件快速识别方法，运行速度快，识别精准。
3.本方案弹性大，可根据自己的研究需要进行改进，一方面，对于计算极端高低温阈值部分，可根据需要使用不同方法界定阈值，然后将数据在后续步骤中替换即可；另一方面，对于设定的最小影响面积阈值a0，可根据研究需要进行更改，方便快捷；对于识别极端低温事件，id_tmin、id_tmax函数已逐步注释，易于理解和更改，可根据需要改成计算持续多天的极端高低温事件。
4.本方法仅对极端高低温事件的发生位置(纬度
‑
lat、经度
‑
lon)、日最低或最高气温、开始时间、结束时间和影响的格点面积进行统计；因为进行极端高低温事件时空变化研究的方法很多，本方案并未提出具体计算强度、频次等方法，因此，可根据自己研究需要在统计出来的各极端高低温事件参数的基础上进行进一步计算分析。


技术实现要素：

5.本发明要解决的技术问题是以上所述的技术问题，提供一种运行速度快，识别精准的基于栅格数据的极端低温事件快速识别方法。
6.为解决上述技术问题，本发明提供的技术方案为：—种基于栅格数据快速识别极端高/低温事件的方法，所述的识别极端高低温事件的方法为：
7.步骤1：导入逐日最低或最高气温栅格数据，如果是点数据(行为时间，列为格点)，需运行matlab软件函数raw2matrix将其转为面数据(lat
×
lon
×
time)；若导入数据为面数据，则直接进行步骤2；
8.步骤2：1)、计算极端低温阈值，剔除高温地区(格点年均日最低气温>研究区平均日最低气温)；选取逐年冬季(12
‑
2月)日最低气温序列的第5百分位数的多年平均值作为极端低温的阈值，若界定极端低温阈值的方法为绝对阈值法或研究区无高温地区，则可无需剔除高温地区步骤直接进行步骤3；2)、计算极端高温阈值，剔除高寒地区(格点年均日最高气温<研究区平均最高气温)；选取逐年夏季(6
‑
8月)日最高气温序列的第95百分位数的多年平均值作为极端高温的阈值，若界定极端高温阈值的方法为绝对阈值法或研究区无高寒地区，则可无需剔除高寒地区步骤直接进行步骤3；
9.步骤3：设定极端低温事件或极端高温事件最小影响面积阈值a0，剔除影响面积小于a0的事件。设定最小面积阈值a0，若一次低温或高温事件的影响面积小于a0则进行剔除。因此，此步骤首先对整个时间尺度进行扫描，剔除小于a0的极端低温事件。最小面积a0阈值
可变，可根据研究区的大小来进行设定，最小面积a0阈值设为25000km2，若进行中国尺度研究，可设为150000km2，若进行全球尺度研究，可设为500000km2；
10.步骤4：通过matlab软件的核心函数id_tmin或id_tmax识别极端低温或高温事件并输出相关信息。
11.作为改进，所述的matlab软件函数raw2matrix表达式为：
12.[lon,lat,obs]＝raw2matrix(lon,lat,data)。
[0013]
作为改进，所述的步骤4中的核心函数id_tmin、id_tmax为：
[0014]
[tmin_event,day,days]＝id_tmin(tmin,lat,lon,p5,time,d,r)；
[0015]
[tmax_event,day,days]＝id_tmax(tmax,lat,lon,p95,time,d,r)；
[0016]
其中，输入变量：tmin为逐日最低气温栅格数据(面数据，lat
×
lon
×
time)，tmax为逐日最高气温栅格数据(面数据，lat
×
lon
×
time)，lat、lon分别为步骤1中采用函数raw2matrix生成的经纬度面数据，p5为步骤2中计算得到的第5百分位阈值，p95为步骤2中计算得到的第95百分位阈值，time为日最低或最高气温数据对应的时间，d为识别极端高温或极端低温事件的持续天数，r为数据的空间分辨率；
[0017]
输出变量：tmin_event为该时间段内所有发生的极端低温事件，tmax_event为该时间段内所有发生的极端高温事件，共有5行；每一列对应的为一次极端低温或极端高温事件，每一行分别对应每次极端低温或极端高温事件的：纬度(lat)、经度(lon)、时间(page/time)、日最低气温(tmin)、持续天数。day和days分别对应的为持续d天及持续d天以上的极端低温事件，共有8行；每一列对应的为一次极端低温事件，每一行分别对应每次极端低温事件的：纬度(lat)、经度(lon)、日最低气温(tmin)或日最高气温(tmax)、开始时间、结束时间、格点影响面积、总影响面积。
[0018]
本发明与现有技术相比的优点在于：通过点数据转面数据函数raw2matrix，对数据进行批量处理，大大优化并减少了程序的运行时间，识别极端高低温事件函数id_tmin、id_tmax，此函数主要采用了matlab内置函数bwconncomp，以连通分量的形式来识别极端低温或高温事件，大大减少了程序的处理时间，能够快速识别极端低温或极端高温事件并统计不同持续事件内极端低温或极端高温事件的发生位置(纬度
‑
lat、经度
‑
lon)、日最低或最高气温、开始时间、结束时间和影响的格点面积，效率高、精度高。
附图说明
[0019]
图1是本发明—种基于栅格数据快速识别极端高/低温事件的方法极端高低温事件过程的示意图。
[0020]
图2是本发明—种基于栅格数据快速识别极端高/低温事件的方法识别极端低温事件工作的流程图。
[0021]
图3是本发明—种基于栅格数据快速识别极端高/低温事件的方法识别极端高温事件工作的流程图。
具体实施方式
[0022]
下面结合附图对本发明做进一步的详细说明。
[0023]
结合附图1～2，—种基于栅格数据快速识别极端高/低温事件的方法，所述的识别
极端高低温事件的方法为：
[0024]
步骤1：导入逐日最低或最高气温栅格数据，如果是点数据(行为时间，列为格点)，需运行matlab软件函数raw2matrix将其转为面数据(lat
×
lon
×
time)；若导入数据为面数据，则直接进行步骤2；
[0025]
步骤2：1)、计算极端低温阈值，剔除高温地区(格点年均日最低气温>研究区平均日最低气温)；选取逐年冬季(12
‑
2月)日最低气温序列的第5百分位数的多年平均值作为极端低温的阈值，若界定极端低温阈值的方法为绝对阈值法或研究区无高温地区，则可无需剔除高温地区步骤直接进行步骤3，极端低温阈值计算方法可变，可根据研究需要界定极端低温阈值并把结果导入本程序进行识别极端低温事件；2)、计算极端高温阈值，剔除高寒地区(格点年均日最高气温<研究区平均最高气温)；选取逐年夏季(6
‑
8月)日最高气温序列的第95百分位数的多年平均值作为极端高温的阈值，若界定极端高温阈值的方法为绝对阈值法或研究区无高寒地区，则可无需剔除高寒地区步骤直接进行步骤3，极端高温阈值计算方法可变，可根据研究需要界定极端高温阈值并把结果导入本程序进行识别极端高温事件；
[0026]
步骤3：设定极端低温事件或极端高温事件最小影响面积阈值a0，剔除影响面积小于a0的事件。设定最小面积阈值a0，若一次低温或高温事件的影响面积小于a0则进行剔除。因此，此步骤首先对整个时间尺度进行扫描，剔除小于a0的极端低温事件。最小面积a0阈值可变，可根据研究区的大小来进行设定，最小面积a0阈值设为25000km2，若进行中国尺度研究，可设为150000km2，若进行全球尺度研究，可设为500000km2；
[0027]
步骤4：通过matlab软件的核心函数id_tmin或id_tmax识别极端低温或高温事件并输出相关信息，一次极端低温或极端高温事件应包括发生、发展与消亡的过程。其中，发展过程包括极端低温或极端高温的移动、分裂与重新组合。由于同一极端低温或极端高温事件在短时间内移动距离较短，在发展这一过程中空间范围存在重叠，因此，对于相邻时间极端低温或极端高温对象，若其空间范围拓扑相交，则认为极端低温对象或极端高温对象属于同一次极端低温或极端高温事件，如图1所示。
[0028]
所述的matlab软件函数raw2matrix表达式为：
[0029]
[lon,lat,obs]＝raw2matrix(lon,lat,data)，将点数据转换成面数据。
[0030]
所述的步骤4中的核心函数id_tmin、id_tmax为：
[0031]
[tmin_event,day,days]＝id_tmin(tmin,lat,lon,p5,time,d,r)；
[0032]
[tmax_event,day,days]＝id_tmax(tmax,lat,lon,p95,time,d,r)；
[0033]
其中，输入变量：tmin为逐日最低气温栅格数据(面数据，lat
×
lon
×
time)，tmax为逐日最高气温栅格数据(面数据，lat
×
lon
×
time)，lat、lon分别为步骤1中采用函数raw2matrix生成的经纬度面数据，p5为步骤2中计算得到的第5百分位阈值，p95为步骤2中计算得到的第95百分位阈值，time为日最低或最高气温数据对应的时间，d为识别极端高温或极端低温事件的持续天数，r为数据的空间分辨率；
[0034]
输出变量：tmin_event为该时间段内所有发生的极端低温事件，tmax_event为该时间段内所有发生的极端高温事件，共有5行；每一列对应的为一次极端低温或极端高温事件，每一行分别对应每次极端低温或极端高温事件的：纬度(lat)、经度(lon)、时间(page/time)、日最低气温(tmin)、持续天数。day和days分别对应的为持续d天及持续d天以上的极端低温事件，共有8行；每一列对应的为一次极端低温事件，每一行分别对应每次极端低温
事件的：纬度(lat)、经度(lon)、日最低气温(tmin)或日最高气温(tmax)、开始时间、结束时间、格点影响面积、总影响面积。
[0035]
以上对本发明及其实施方式进行了描述，这种描述没有限制性，附图中所示的也只是本发明的实施方式之一，实际的结构并不局限于此。总而言之如果本领域的普通技术人员受其启示，在不脱离本发明创造宗旨的情况下，不经创造性的设计出与该技术方案相似的结构方式及实施例，均应属于本发明的保护范围。