1.本发明涉及一种预警算法，尤其是一种基于真实孔径雷达的配准和预警算法。
背景技术：
2.矿山开采会产生滑坡，滑坡深部为煤矿采空区，煤矿开采导致岩土体应力改变，产生变形，这种危害危险性极大，为防止这种灾害现今均是利用边坡孔径雷达对矿区数据采集并数据，以期提前预警，现场雷达采集数据后，通过接口调用将监测数据上传至系统服务端，服务端接收到数据后进行数据存储以及数据分析，用户通过客户端对历史数据或实时数据进行查询，此时就需要一种预警算法能够将孔径雷达的配准结果在客户端进行展示，同时通过点选择或区域选择对选择的数据进行曲线分析，计算预计滑坡时间。
技术实现要素：
3.本发明提供了一种基于真实孔径雷达的配准和预警算法，解决了现有技术中存在的缺乏此种预警算法的问题。
4.本发明为解决上述技术问题所采用的技术方案是：以雷达扫描的数据为核心，通过对宽、高、型变量、位置偏移量的计算来确定该雷达扫描区域的预警等级，具体步骤如下：
5.1)在地形上选择一个区域，设选择的图形上横坐标是x，纵坐标是y；
6.2)设定相对于雷达采集数据的横轴偏移量是offx，设定相对于雷达采集数据的纵轴偏移量是offy；
7.3)设定分析的相对宽度是w，设定分析的相对高度是h，雷达扫描总数据量为length，设定地形的相对偏移量为i，设t为雷达扫描选中地形的平均形变量；
8.4)设a为实际的数据位置横向坐标，通过a＝(x-offx)/w计算；
9.5)设d为实际的数据位置纵向坐标，通过d＝(y-offy)/h计算；
10.6)设岩层密度占用比例为avg，通过avg＝(a/2+d/2)/(i+length)计算；
11.7)设预警等级g，通过g＝f(a,d,i)计算；
12.8)设定预警值：
13.黄色预警：s1＝(30+avg
×
length)
×214.橙色预警：s2＝(60+avg
×
length)
×215.红色预警：s3＝(90+avg
×
length)
×216.无预警：d
×
100+a＞length
17.9)用t和s1、s2、s3比较判断预警等级如下：
18.当r＝0时，未屏蔽，判断关系：
19.t＜s1，输出为0，即g＝f(a,d,i)＝0，则预警等级为无；
20.t≥s1，t＜s2，输出为1，即g＝f(a,d,i)＝1，则预警等级为黄色预警；
21.t≥s2，t＜s3，输出为2，即g＝f(a,d,i)＝2，则预警等级为橙色预警；
22.t≥s3，输出为3，即g＝f(a,d,i)＝3，则预警等级为红色预警；
23.d
×
100+a＞length,则预警等级为无，无预警；
24.当雷达扫描被屏蔽，判断关系：
25.输出为0，则预警等级为无。
26.其中函数值为x∈n，y∈n，offx∈[0,200],offy∈[0,200],w∈n，h∈n，length∈[1,+∞]，i∈n，a∈n，d∈n，t∈n，g∈{0，1，2，3}，其中n为自然数。
[0027]
优选的，根据输出的g数值判断对应的预警等级，判断危险等级，预警等级g数值越小，预警等级越低，预警等级g数值越大，预警等级越高。
[0028]
优选的，所选雷达扫描区域的区域面积一定时，雷达扫描总数据量和地形的相对偏移量和预警值成正比关系。
[0029]
本发明采用上述结构，具有以下的优点：可以通过运算把握矿山预警等级，客观的反映矿山的地形变化，以便于采取有针对性的措施；雷达监测站点定时向服务器发送数据，通过运算把握该地区的预警等级，客观的反应当前地区的地形变化，准确度较高，对及时有效采取防范措施有重大意义。
具体实施方式
[0030]
为能清楚说明本方案的技术特点，下面通过具体实施方式，对本发明进行详细阐述。
[0031]
以雷达扫描的数据为核心，通过对宽、高、型变量、位置偏移量的计算来确定该雷达扫描区域的预警等级，具体步骤如下：
[0032]
1)在地形上选择一个区域，设选择的图形上横坐标是x，纵坐标是y；
[0033]
2)设定相对于雷达采集数据的横轴偏移量是offx，设定相对于雷达采集数据的纵轴偏移量是offy；
[0034]
3)设定分析的相对宽度是w，设定分析的相对高度是h，雷达扫描总数据量为length，设定地形的相对偏移量为i，设t为雷达扫描选中地形的平均形变量；
[0035]
4)设a为实际的数据位置横向坐标，通过a＝(x-offx)/w计算；
[0036]
5)设d为实际的数据位置纵向坐标，通过d＝(y-offy)/h计算；
[0037]
6)设岩层密度占用比例为avg，通过avg＝(a/2+d/2)/(i+length)计算；
[0038]
7)设预警等级g，通过g＝f(a,d,i)计算；
[0039]
8)设定预警值：
[0040]
黄色预警：s1＝(30+avg
×
length)
×2[0041]
橙色预警：s2＝(60+avg
×
length)
×2[0042]
红色预警：s3＝(90+avg
×
length)
×2[0043]
无预警：d
×
100+a＞length
[0044]
9)用t和s1、s2、s3比较判断预警等级如下：
[0045]
当r＝0时，未屏蔽，判断关系：
[0046]
t＜s1，输出为0，即g＝f(a,d,i)＝0，则预警等级为无；
[0047]
t≥s1，t＜s2，输出为1，即g＝f(a,d,i)＝1，则预警等级为黄色预警；
[0048]
t≥s2，t＜s3，输出为2，即g＝f(a,d,i)＝2，则预警等级为橙色预警；
[0049]
t≥s3，输出为3，即g＝f(a,d,i)＝3，则预警等级为红色预警；
[0050]
d
×
100+a＞length,则预警等级为无，无预警；
[0051]
当雷达扫描被屏蔽，判断关系：
[0052]
输出为0，则预警等级为无。
[0053]
其中函数值为x∈n，y∈n，offx∈[0,200],offy∈[0,200],w∈n，h∈n，length∈[1,+∞]，i∈n，a∈n，d∈n，t∈n，g∈{0，1，2，3}，其中n为自然数。
[0054]
根据输出的g数值判断对应的预警等级，判断危险等级，预警等级g数值越小，预警等级越低，预警等级g数值越大，预警等级越高。
[0055]
所选雷达扫描区域的区域面积一定时，雷达扫描总数据量和地形的相对偏移量和预警值成正比关系。
[0056]
实施例：
[0057]
将雷达对矿区采集的数据带入预警等级函数g＝f(a,d,i)中,
[0058]
第一组数据输入：x＝100，y＝50，offx＝10，offy＝20
[0059]
w＝10，h＝10，length＝1000，i＝100，t＝163
[0060]
计算过程：
[0061]
中间带入计算：
[0062]
a＝(x-offx)/w＝(100-10)/10＝9
[0063]
d＝(y-offy)/h＝(50-20)/10＝3
[0064]
avg＝(a/2+d/2)/(i+length)＝(9/2+3/2)/(100+1000)＝0.005454
[0065]
计算预警值：
[0066]
无预警：(d
×
100+a＝3
×
100+9＝309)＜(length＝1000)
[0067]
黄色预警：s1＝(30+avg
×
length)
×
2＝(30+6/1100
×
1000)
×
2＝70.908
[0068]
橙色预警：s2＝(60+avg
×
length)
×
2＝(60+avg
×
length)
×
2＝130.908
[0069]
红色预警：s3＝(90+avg
×
length)
×
2＝(90+avg
×
length)
×
2＝190.908
[0070]
t≥s2，t＜s3输出为2，即g＝f(a,d,i)＝f(9,3,100)＝2
[0071]
输出的结果：g＝2
[0072]
预警等级为二级，预警状态为橙色预警。
[0073]
第二组数据输入：
[0074]
x＝90，y＝60，offx＝30，offy＝10，w＝20，h＝10
[0075]
length＝200，i＝10，t＝60
[0076]
计算过程：
[0077]
中间带入计算：
[0078]
a＝(x-offx)/w＝(90-30)/20＝3
[0079]
d＝(y-offy)/h＝(60-10)/10＝3
[0080]
avg＝(a/2+d/2)/(i+length)＝(3/2+3/2)/(100+1000)＝0.0142
[0081]
计算预警值：
[0082]
无预警：d
×
100+a＝3
×
100+3＝303)＞(length＝200)
[0083]
输出的结果：g＝0
[0084]
预警等级为无，无预警。
[0085]
第三组数据输入：
[0086]
x＝110，y＝65，offx＝10，offy＝5，w＝5，h＝30
[0087]
length＝400，i＝40，t＝260
[0088]
计算过程：
[0089]
中间带入计算：
[0090]
a＝(x-offx)/w＝(110-10)/5＝20
[0091]
d＝(y-offy)/h＝(65-5)/30＝2
[0092]
avg＝(a/2+d/2)/(i+length)＝(20/2+2/2)/(40+400)＝0.025
[0093]
计算预警值：
[0094]
无预警：(d
×
100+a＝2
×
100+20＝220)＜(length＝400)
[0095]
黄色预警：s1＝(30+avg
×
length)
×
2＝(30+0.025
×
400)
×
2＝80
[0096]
橙色预警：s2＝(60+avg
×
length)
×
2＝(60+0.025
×
400)
×
2＝140
[0097]
红色预警：s3＝(90+avg
×
length)
×
2＝(90+0.025
×
400)
×
2＝200
[0098]
t≥s3,输出为3，即g＝f(a,d,i)＝f(20,2,40)＝3
[0099]
输出的结果：g＝3
[0100]
预警等级为三级，预警状态为红色预警。
[0101]
第四组数据输入：
[0102]
x＝44，y＝66，offx＝2，offy＝4，w＝2，h＝2
[0103]
length＝400，i＝5，t＝270
[0104]
计算过程：
[0105]
中间带入计算：
[0106]
a＝(x-offx)/w＝(44-2)/2＝21
[0107]
d＝(y-offy)/h＝(66-4)/2＝31
[0108]
avg＝(a/2+d/2)/(i+length)＝(21/2+31/2)/(5+400)＝0.0641
[0109]
计算预警值：
[0110]
无预警：(d
×
100+a＝31
×
100+21＝331)＜(length＝400)
[0111]
黄色预警：s1＝(30+avg
×
length)
×
2＝(30+0.0641
×
400)
×
2＝111.28
[0112]
橙色预警：s2＝(60+avg
×
length)
×
2＝(60+0.0641
×
400)
×
2＝145.64
[0113]
红色预警：s3＝(90+avg
×
length)
×
2＝(90+0.0641
×
400)
×
2＝231.28
[0114]
t≥s3,输出为3，即g＝f(a,d,i)＝f(21,31,5)＝3
[0115]
输出的结果：g＝3
[0116]
预警等级为三级，预警状态为红色预警。
[0117]
第五组数据输入：
[0118]
x＝33，y＝55，offx＝6，offy＝2，w＝10，h＝10
[0119]
length＝1000，i＝20，t＝146
[0120]
计算过程：
[0121]
中间带入计算：
[0122]
a＝(x-offx)/w＝(33-6)/10＝2.7
[0123]
d＝(y-offy)/h＝(55-2)/10＝5.3
[0124]
计算比例：
[0125]
avg＝(a/2+d/2)/(i+length)＝(2.7/2+5.3/2)/(20+1000)＝0.0039
[0126]
计算预警值：
[0127]
无预警：(d
×
100+a＝5.3
×
100+2.7＝532.7)＜(length＝1000)
[0128]
黄色预警：s1＝(30+avg
×
length)
×
2＝(30+0.0039
×
1000)
×
2＝37.8
[0129]
橙色预警：s2＝(60+avg
×
length)
×
2＝(60+0.0039
×
1000)
×
2＝127.8
[0130]
红色预警：s3＝(90+avg
×
length)
×
2＝(90+0.0039
×
1000)
×
2＝187.8
[0131]
t≥s2,并且t＜s3,输出为2，即g＝f(a,d,i)＝f(2.7,5.3,20)＝2
[0132]
输出的结果：g＝2
[0133]
预警等级为二级，预警状态为橙色预警。
[0134]
此时，可以根据输出的g数值判断对应的预警等级，然后根据预警等级计算预计滑坡时间，做出相应的预防措施。
[0135]
其中，雷达和雷达扫描技术均是现有技术。
[0136]
上述具体实施方式不能作为对本发明保护范围的限制，对于本技术领域的技术人员来说，对本发明实施方式所做出的任何替代改进或变换均落在本发明的保护范围内。
[0137]
本发明未详述之处，均为本技术领域技术人员的公知技术。