CN112348960B - 一种适用于全球空间范围的空域冲突检测方法 - Google Patents

一种适用于全球空间范围的空域冲突检测方法 Download PDF

Info

Publication number
CN112348960B
CN112348960B CN202011354731.8A CN202011354731A CN112348960B CN 112348960 B CN112348960 B CN 112348960B CN 202011354731 A CN202011354731 A CN 202011354731A CN 112348960 B CN112348960 B CN 112348960B
Authority
CN
China
Prior art keywords
airspace
grid
level
conflict
code
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011354731.8A
Other languages
English (en)
Other versions
CN112348960A (zh
Inventor
万路军
徐鑫宇
黄阿倩
戴江斌
沈堤
霍丹
陈致远
蔡明�
高志周
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Air Force Engineering University of PLA
Original Assignee
Air Force Engineering University of PLA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Air Force Engineering University of PLA filed Critical Air Force Engineering University of PLA
Priority to CN202011354731.8A priority Critical patent/CN112348960B/zh
Publication of CN112348960A publication Critical patent/CN112348960A/zh
Application granted granted Critical
Publication of CN112348960B publication Critical patent/CN112348960B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • G06T17/205Re-meshing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T9/00Image coding
    • G06T9/001Model-based coding, e.g. wire frame

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Multimedia (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明提供一种适用于全球空间范围的空域冲突检测方法,包括:构建空域时间维编码,检测时间维冲突;构建空域高度范围编码,检测高度范围冲突;建立地球表面剖分网格体系;网格化空域;判断网格编码集合间有无交集;判断空域间水平距离;进行频率维冲突检测;本发明方法实现了空域在时间维、空间维和频率维的一体化冲突检测;提高空域表达的精度及冲突检测准确度;计算量小,提高了实时性;本发明方法基于网格编码,编码之间采用比较、异或、求并运算,而非传统的基于经纬度坐标的浮点计算,更符合计算机的计算逻辑,提高运算效率。

Description

一种适用于全球空间范围的空域冲突检测方法
技术领域
本发明属于空中交通管制技术领域,特别涉及一种空域冲突检测方法。
背景技术
随着国民经济的发展,各领域、各行业都有航空活动需求,除军用航空需要划设各类训练和演习空域、民用航空不断开辟新的航路航线之外,各类航空旅游、航空摄影等需求也日益激增,加上不断增多的需要单独划设空域的无人机飞行空域,使得不断增长的用空需求与空域资源的有限性的矛盾日益凸显,空域正变得越来越拥挤。各类空域相互交织,必然导致空域之间的冲突成倍增加,如何高效、快速、准确地进行空域冲突检测,对于提高飞行安全、防止危及飞行安全的情况发生至关重要。冲突检测可分为飞行冲突检测和空域冲突检测两种。前者的研究已有半个多世纪,已有大量的研究文献,适用范围是合作型航空器之间,通过遍历解算航迹的多元方程,在航迹预测的基础上进行冲突检测,此类方法计算量大、耗时长,故无法应用于飞行流量大的区域和军民航之间冲突检测。而空域冲突检测研究的是航空器活动所在的空域有无冲突,而非航空器本身,通过空域间隔将各类航空器的活动空域隔离开,适用范围更广。
空域的属性包括时间属性、空间属性和频率属性,空间属性又包括空域的高度范围和水平范围。因此进行空域冲突检测必须在时间维、高度范围、水平范围和频率维都要检测。
现有的空域冲突检测方法包括图形相交型和栅格型。图形求交型的原理是将空域看成形状各异的图形,判断图形之间线段有无相交进行冲突判定,如专利CN 104050374 B以及专利CN 108922250 A。该方法的不足是没有考虑空域的安全间隔,当空域所代表的图形没有相交时,也可能出现空域间隔小于安全间隔而存在冲突的情况,造成“漏检”。栅格型是构建空域的栅格模型,并赋予栅格以内在属性,利用栅格的内在属性判定空域有无冲突,如专利CN 111477034 A。该方法的不足是没有考虑空域大小与网格大小的关系,而统一用一个能包含空域大小的矩形网格表示空域,对空域的描述过于粗放,影响检测结果的准确性。此外,以上方法均没有考虑空域的使用频率之间的冲突。
发明内容
为了克服现有技术的不足,本发明提供一种适用于全球空间范围的空域冲突检测方法,包括:
步骤1:构建空域时间维编码,检测时间维冲突;
步骤2:构建空域高度范围编码,检测高度范围冲突;
步骤3:建立地球表面剖分网格体系;
步骤4:网格化空域;
步骤5:判断网格编码集合间有无交集;
步骤6:判断空域间水平距离;
步骤7:进行频率维冲突检测。
进一步的,步骤1根据空域的使用时间范围,构建格式统一的时间维编码,便于计算机识别和读取;若空域的时间维无冲突,直接判定为空域之间无冲突,不需要进行下面步骤的判定;若空域的时间维有冲突,则转入步骤2;步骤2对时间维存在冲突的空域,根据空域的高度范围,构建格式统一的高度范围编码;若空域的高度范围无冲突,则转入步骤7;若空域的高度范围有冲突,则转入步骤3;在步骤3中进行第1至7层级剖分与编码;而后转换经纬度坐标与网格编码;步骤4中根据每个空域的类型和大小,用不同层级的多个网格去填充空域;步骤5筛选出必然存在冲突的空域和可能存在冲突的空域;步骤6中将计算出的空域间距离与规定的安全间隔对比,若大于安全间隔,则判定为空域在水平范围上不存在冲突;否则,空域之间在水平范围上存在冲突;步骤7中对时间维有冲突而空间维无冲突的所有空域,进行频率维冲突检测。
进一步的,
步骤1包括:
步骤1-1:时间维编码;
步骤1-2:检测时间维冲突;
步骤2包括:
步骤2-1:高度范围编码;
步骤2-2:检测高度范围冲突;
步骤3包括:
步骤3-1:进行第1层级剖分与编码;
步骤3-2:进行第2~7层级剖分与编码;
步骤3-3:转换经纬度坐标与网格编码;
步骤4包括:
步骤4-1:计算空域的外包矩形;
步骤4-2:计算空域的最小外包网格;
步骤4-3:计算包含要进行水平范围冲突检测的多个空域的最小外包网格;
步骤4-4:计算空域边界所在的网格;
步骤4-5:计算空域内部所在的网格;
步骤4-6:网格聚合;
步骤6包括:
步骤6-1:判定空域间相对方位;
步骤6-2:确定两空域在相对方位上最边缘网格的行列坐标;
步骤6-3:通过坐标计算出空域间水平距离;
步骤7包括:
步骤7-1:统一频率维的编码格式;
步骤7-2:进行编码运算;
步骤7-3:根据步骤7-2的运算结果,进行冲突判定。
进一步的,步骤1-1中时间编码为一串16位的数字,前8位数字代表为空域的启用时间,后8位数字代表空域的结束使用时间,中间用“-”隔开;步骤1-2中两两空域之间进行时间维冲突检测,两个空域只有时间范围不重叠,才无冲突;步骤2-1中高度编码用大写英文字母H加上三位数字表示,数字的单位为百米;步骤2-2中两两空域之间进行高度范围冲突检测,两个空域只有高度范围不重叠,才无冲突;步骤3-1中选择本初子午线与赤道交点为剖分原点,用4°×4°的第一层级网格将纬度范围为[-90°,90°]、经度范围为[-180°,180°]的地球表面剖分为46×90份;步骤3-2中上一层级的网格为父网格,下一层级的网格为子网格;每1个上一层级父网格平均分成16份,得到下一层级16个子网格,对16个子网格按照皮亚诺空间填充曲线进行16进制编码,直到剖分至第7层级;步骤4-1中分别过纬度最小的点和纬度最大的点作两条正东正西的直线;分别过经度最小的点和经度最大的点作两条正南正北的直线,四条线的相交范围即该空域的外包矩形;步骤4-2中比较出网格面积大于步骤4-1中空域外包矩形的面积且与该面积最接近的网格层级;步骤4-3中选择一个最小外包网格的编码,取其编码均相等前多位编码,编码对应的网格即为多个空域的最小外包网格;步骤4-4中将步骤4-3得到的最小外包网格剖分至指定层级;步骤4-5中将空域边界网格行列坐标集合中行坐标相等的网格分在同一组,该组中列坐标最小为Pmin,最大为Pmax,按步长为1填充Pmin~Pmax中的列坐标,行坐标不变,按照纬度由低到高、经度由西向东的顺序逐行进行扫描,直到空域内部填充完毕,获得空域内部网格的行列坐标集合;步骤4-6中将同属于1个父网格的16个子网格聚合成1个父网格,用1个父网格编码代替16个子网格编码;步骤5中当空域对应的网格编码集合间有交集,即两空域包含同一个网格时,空域在水平范围内必然存在冲突,不再进行下一步冲突判定;只有空域对应的网格编码集合间无交集,即两空域不包含同一个网格时,空域在水平范围内才有可能无冲突,转入步骤6;步骤6中将计算出的距离与《间隔规定》和《中华人民共和国飞行基本规则》规定的安全间隔对比,若大于安全间隔,则判定为空域在水平范围上不存在冲突;否则,空域之间在水平范围上存在冲突。
进一步的,步骤1-1中在前8位和后8位数字中,第1-2位数字代表某月,第3-4位数字代表某月,第5-8位数字代表某时某分,采用24小时制,空域的时间维编码默认省略当年年份;若空域使用时间跨年,则在空域启用时间和结束使用时间编码前面各加上4位数字,代表年份编码;步骤1-2中在两个空域中任选一个空域,记为空域A,取空域A的前8位编码,先与另外一个空域B的后8位编码进行逐位比较,最终比较结果为Result-1,再取空域A的后8位编码与空域B的前8位编码进行逐位比较,最终比较结果为Result-2,对Result-1和Result-2进行异或运算,判断时间维是否冲突;步骤2-1中高度范围的编码共8位,其中前4位为空域的下限高度,后4位为空域的上限高度;步骤2-2中在时间维上存在冲突的两个空域中任选一个空域,记为空域M,取空域M高度范围编码的前4位,先与另外一个空域N的后4位编码进行逐位比较,最终比较结果为Result-3;再取空域A的后4位编码与空域N的前4位编码进行逐位比较,最终比较结果为Result-4,对Result-3和Result-4进行异或运算,判断高度范围是否冲突;步骤3-1中第一层级网格序号在纬度方向上用A~Y和a~y共48个英文字母从赤道按照纬度由低到高编码,编码跳过I、O、i和o,其中北纬大写、南纬小写;第一层级网格序号在经度方向上从本初子午线开始由西向东用00~89共90个数字编码;步骤3-3中某点经纬度坐标为(L,B),其中L为纬度,B为经度,L和B的表示形式均为度分秒,记作L=LD°LM′LS″,B=BD°BM′BS″;以步骤4-4中以最小外包网格内第N层级网格中左下网格为原点网格,其坐标为(0,0),建立直角坐标系,向东、向北从0开始,按步长为1编号,直到编号到最小外包网格内最东和最北端的第N层级网格,确定空域边界每条线段端点所在网格的行列坐标,端点之间网格的行列坐标集由插值计算得到,去掉空域角点重复的坐标,得到空域边界网格的行列坐标集合,其中列坐标在前,行坐标在后;步骤7-1中频率单位统一为KHz,单位前面的数值统一设定为9位,不足9位的在数值前面补0;步骤7-2中对两空域所使用的频率值先逐位进行异或运算,再对异或运算得到的9个数值逐位进行求并运算;步骤7-3中若最终求并运算的结果为0,则两空域所使用的频率相等,频率之间会存在相互干扰,空域间存在频率冲突;若结果为1,则两空域所使用的频率不相等,频率之间不存在相互干扰,空域间不存在频率冲突。
进一步的,步骤3-1中第1层级网格的边长为512千米;步骤3-2中第2~7级网格的大小分别为1°×1°、16′×16′、4′×4′、1′×1′、16″×16″、4″×4″,网格边长分别为128千米、32千米、8千米、2千米、512米、128米。
进一步的,步骤3-3中转换经纬度坐标与网格编码具体包括:
步骤3-3-1:计算第1层级网格编码;计算公式为:
LD/4°的商为M1、余数为N1,M1+1对应的英文字母为该点所在的第1层级网格的纬度方向编号,北纬大写,南纬小写;BD/4°商为m1、余数为n1,m1即为该点所在的第1层级网格的经度方向编号;两个编号组合成该点所在的第1层级网格编码;
步骤3-3-2:计算第2层级网格编码;计算公式为:
步骤1的余数N1除以1°得到数值M2,M2即为该点所在的第2层级网格的纬度方向编号;n1除以1°得到数值m2,m2即为该点所在的第2层级网格的经度方向编号,根据两个编号对应得出第2层级网格编码;
步骤3-3-3:计算第3层级网格编码,计算公式为:
LM/16′的商为M3、余数为N3,M2即为该点所在的第3层级网格的纬度方向编号;BM/16′商为m3、余数为n3,m3即为该点所在的第3层级网格的经度方向编号,根据两个编号对应出第3层级网格编码;
步骤3-3-4:计算出该点位于第4层级网格的编码,计算公式为:
N3/4′的商为M4、余数为N4,M4即为该点所在的第4层级网格的纬度方向编号;n3/4′商为m4、余数为n4,m4即为该点所在的第4层级网格的经度方向编号,根据两个编号对应出第4层级网格编码;
步骤3-3-5:计算出该点位于第5层级网格的编码,计算公式为:
N4/1′的商为M5,M5即为该点所在的第5层级网格的纬度方向编号;n4/1′商为m5、余数为n5,m5即为该点所在的第5层级网格的经度方向编号,根据两个编号对应出第5层级网格编码;
步骤3-3-6:计算出该点位于第6层级网格的编码,计算公式为:
LS/16″的商为M6、余数为N6,M6即为该点所在的第6层级网格的纬度方向编号;BS/16″商为m6、余数为n6,m6即为该点所在的第6层级网格的经度方向编号,根据两个编号对应出第6层级网格编码;
步骤3-3-7:计算出该点位于第7层级网格的编码,计算公式为:
N6/4″的商为M7,M7即为该点所在的第7层级网格的纬度方向编号;n6/4″商为m7,m7即为该点所在的第7层级网格的经度方向编号,根据两个编号对应出第7层级网格编码。
进一步的,步骤6-1中根据两空域的中心点经纬度坐标,判定两空域的方位关系;步骤6-2中根据表达精度对应的最小剖分层级N,确定两空域在相对方位上最边缘的第N级网格编码,取最后两位编码转化为一维二进制编码,再拆分为二维二进制编码,最后转化为二维十进制编码,以最小外包网格的左下角为坐标原点(0,0),得到两空域在相对方位上最边缘的第N级网格的坐标;步骤6-3中通过坐标计算两空域在相对方位上最边缘第N级网格之间的横向间距与纵向间距,使用勾股定理求得两空域间水平距离,若水平距离小于水平安全间隔,则两空域在水平范围存在冲突,判定两空域空间上有冲突;否则,两空域在水平范围不冲突,判定两空域在空间维上不冲突。
进一步的,步骤6-3中水平安全间隔包括横向安全间隔、纵向安全间隔,水平安全间隔为10公里。
本发明方法实现了空域在时间维、空间维和频率维的一体化冲突检测,提高了冲突检测的准确度和实时性。
附图说明
图1为本发明方法流程框图;
图2为本发明中地球表面剖分金字塔示意图;
图3为本发明中地球表面剖分网格体系和编码方案平面示意图;
图4为本发明中计算空域最小外包网格示意图;
图5(a)为本发明中填充空域内部网格时的扫描路径表达示意图;
图5(b)为本发明中16个子网格聚合成1个父网格示意图;
图6为本发明中判定空域间相对方位示意图;
图7为本发明中计算空域间水平距离示意图;
图8为本发明方法程序实现流程框图。
具体实施方式
为使本发明的目的、技术方案以及优势更加明晰,下面结合附图和实施例对本发明的内容进行进一步详细说明。
本发明的基本原理是通过在空域的时间维、高度范围和频率维进行格式统一的编码,便于计算机识别和读取,利用编码间比较、异或和求并运算,实现空域在时间维、高度范围和频率维上的冲突检测。通过构建地球表面剖分网格体系,并对网格进行编码,每个网格都对应唯一的编码,对空域进行网格化表达。将空域表达为网格集合,由网格与编码的唯一对应关系,得到空域的编码集合,通过编码集合求交集,对无交集的空域对,通过相对方位上最边缘网格的编码,计算出空域间水平间隔,进行水平范围上空域冲突检测。进而实现空域在时间维、高度范围、水平范围和频率维冲突检测的目的。
如图1所示,本发明方法包括以下步骤:
步骤1:构建空域时间维编码,检测时间维冲突。即根据空域的使用时间范围,构建格式统一的时间维编码,便于计算机识别和读取。若空域的时间维无冲突,直接判定为空域之间无冲突,不需要进行下面步骤的判定;若空域的时间维有冲突,则转入步骤2。
具体方法为:
步骤1-1:对时间维编码。时间编码为一串16位的数字,前8位数字代表为空域的启用时间,后8位数字代表空域的结束使用时间,中间用“-”隔开。在前8位和后8位数字中,第1-2位数字代表某月,第3-4位数字代表某月,第5-8位数字代表某时某分,采用24小时制。空域的时间维编码默认省略当年年份。例如今年是2020年,06100800-06111400代表该空域的使用期限为2020年6月10日早8点到2020年6月11日下午2点。若空域使用时间跨年,则在空域启用时间和结束使用时间编码前面各加上4位数字,代表年份编码,例如202012300800-202101052000,代表该空域的使用期限为2020年12月30日早8点到2021年1月5日晚8点。
步骤1-2:检测时间维冲突。两两空域之间进行时间维冲突检测,两个空域只有时间范围不重叠,才无冲突。因此在两个空域中任选一个空域,记为空域A,取空域A的前8位编码,先与另外一个空域B的后8位编码进行逐位比较,最终比较结果为Result-1。再取空域A的后8位编码与空域B的前8位编码进行逐位比较,最终比较结果为Result-2。若逐位比较结果为等于,相应位的计算结果为空,用符号表示;若逐位比较结果是大于,相应位的计算结果为0;若逐位比较结果是小于,相应位的计算结果为1。一旦数字1或0出现,则不再进行后面位的比较,最终比较结果Result-1和Result-2的取值即为1或0。对Result-1和Result-2进行异或运算,公式为/>其中符号/>代表异或运算。若运算结果为0,则两空域的时间范围不重叠,两空域在时间维上无冲突,直接判定为空域间无冲突,不再进行下一步的冲突检测;若运算结果为1,则两空域的时间范围重叠,两空域在时间维上有冲突,进行下一步高度范围冲突检测。
步骤2:构建空域高度范围编码,检测高度范围冲突。即对时间维存在冲突的空域,根据空域的高度范围,构建格式统一的高度范围编码。若空域的高度范围无冲突,则转入步骤7;若空域的高度范围有冲突,则转入步骤3。
具体方法为:
步骤2-1:高度范围编码。高度编码为了与时间编码隔开,用大写英文字母H加上三位数字表示,数字的单位为百米。为了保持高度编码长度一致,便于计算机识读,如果空域高度固定,则重复编码,高度范围的编码共8位,其中前4位为空域的下限高度,后4位为空域的上限高度。默认高度为标准大气压高度,如需表示真高则用H加下标表示。例如H008H008代表标压高800米,H080H110代表标压高8000米到11000米,H010H020代表真高1000米到2000米。
步骤2-2:检测高度范围冲突。方法与时间维冲突检测方法类似。在时间维上存在冲突的两个空域中任选一个空域,记为空域M,取空域M高度范围编码的前4位,先与另外一个空域N的后4位编码进行逐位比较,最终比较结果为Result-3。再取空域A的后4位编码与空域N的前4位编码进行逐位比较,最终比较结果为Result-4。若逐位比较结果为等于,相应位的计算结果为空;若逐位比较结果大于,相应位的计算结果为0;若逐位比较结果小于,相应位的计算结果为1。一旦数字1或0出现,则不再进行后面位的比较,最终比较结果Result-3和Result-4的取值即为1或0。对Result-3和Result-4进行异或运算,公式为 其中符号/>代表异或运算。若运算结果为0,则两空域的高度范围不重叠,两空域在高度范围上无冲突,转入步骤7;若运算结果为1,则两空域的高度范围重叠,两空域在高度维上有冲突,转入步骤3。
步骤3:建立地球表面剖分网格体系。对于在高度范围也存在冲突的空域,需要进行水平范围冲突检测,首先需要建立地球表面剖分网格体系。即将地球表面剖分为多层级、多尺度的网格,各层级网格有对应的面积,然后按照编码规则对各层级的网格进行编码。地球表面的任一网格都有唯一的编码,网格编码与经纬度坐标之间能够进行快速转换。
在步骤3中,建立地球表面剖分网格体系的具体方法为:
步骤3-1:进行第1层级剖分与编码。如图2所示,选择本初子午线与赤道交点为剖分原点,用4°×4°的第一层级网格将纬度范围为[-90°,90°]、经度范围为[-180°,180°]的地球表面剖分为46×90份。第一层级网格序号在纬度方向上用A~Y和a~y共48个英文字母从赤道按照纬度由低到高编码,为避免字母与数字1和0混淆,编码跳过I、O、i和o,其中北纬大写、南纬小写;第一层级网格序号在经度方向上从本初子午线开始由西向东用00~89共90个数字编码。第1层级网格的边长为512千米。
步骤3-2:进行第2~7层级剖分与编码。上一层级的网格为父网格,下一层级的网格为子网格。每1个上一层级父网格平均分成16份,得到下一层级16个子网格,对16个子网格按照皮亚诺空间填充曲线进行16进制编码,直到剖分至第7层级。所以子网格编码继承上一层级父网格的编码,并比上一层级父网格编码多1位。第2~7级网格的大小分别为1°×1°、16′×16′、4′×4′、1′×1′、16″×16″、4″×4″,网格边长分别为128千米、32千米、8千米、2千米、512米、128米。剖分层级选择为7级,是因为空域间的水平安全间隔为10公里,若空域间水平间隔小于10公里,则空域间在水平范围上有冲突。第7层级网格的边长为128米,如果再往下剖分到第8和第9层级,则第8和第9层级网格边长分别为32米和8米。不管是128米、32米还是8米,相比于10公里来说都已经很小,对于提高空域间水平范围冲突检测的精度来说意义并不大,反而过多的剖分层级会使得网格编码变得冗长,占用更多的计算机内部存储空间,降低编码间计算效率。因此,地球表面剖分网格体系将剖分层级设定为7级。例如,第3层级某16′×16′网格的编码为H07C2,将其平均分成16个4′×4′的第4层级网格,对这16个子网格按照图2中的编码方式进行编码,编码分别为{H07C20,H07C21,H07C22,H07C23,H07C24,H07C25,H07C26,H07C27,H07C28,H07C29,H07C2A,H07C2B,H07C2C,H07C2D,H07C2E,H07C2F}。经过步骤2-1和步骤2-2,将地球剖分成无缝隙、无重叠的网格,每个网格都有唯一的编码以及相对应的面积。各层级网格大小、尺度和数量见表1。
表1
步骤3-3:转换经纬度坐标与网格编码。设某点的经纬度坐标为(L,B),其中L为纬度,B为经度,L和B的表示形式均为度分秒,记作L=LD°LM′LS″,B=BD°BM′BS″。如图3所示,第1层级网格中纵坐标代表纬度方向编码,大写字母为北纬,小写字母为南纬,横坐标代表经度方向编码00~89,图3中只显示了部分第1层级网格及编码。第2~7层级网格中纵坐标代表纬度方向,编号为0、1、2、3,横坐标代表经度方向,编号为0、1、2、3,编号坐标顺序为先纬向,再经向,例如编号坐标(2,0)对应的编号为8。经纬度坐标转化位网格编码的详细步骤为:
步骤3-3-1:计算第1层级网格编码。计算公式为:
由于第1层级网格为4°×4°网格,故计算出LD/4°的商为M1、余数为N1,M1+1对应的英文字母即为该点所在的第1层级网格的纬度方向编号,北纬大写,南纬小写;BD/4°商为m1、余数为n1,m1即为该点所在的第1层级网格的经度方向编号。两个编号组合成该点所在的第1层级网格编码。例如西安市钟楼的经纬度坐标为(34°15′39″N,108°56′32″E),34°/4°=8余2°,108°/4°=27,位于第1层级中北纬向第9个、经向第28个4°×4°网格,所以西安钟楼所在的第1层级网格的编码为J27。
步骤3-3-2:计算第2层级网格编码。计算公式为:
由于第2层级网格为1°×1°网格,故步骤1的余数N1除以1°得到数值M2,M2即为该点所在的第2层级网格的纬度方向编号;n1除以1°得到数值m2,m2即为该点所在的第2层级网格的经度方向编号,根据两个编号对应得出第2层级网格编码。上例中,2°/1°=2,0/1°=0,坐标(2,0)对应的1°×1°网格编号为8,所以西安钟楼所在的第2层级网格的编码为J278。
步骤3-3-3:计算第3层级网格编码。计算公式为:
第3层级网格为16′×16′网格,故LM/16′的商为M3、余数为N3,M2即为该点所在的第3层级网格的纬度方向编号;BM/16′商为m3、余数为n3,m3即为该点所在的第3层级网格的经度方向编号。根据两个编号对应出第3层级网格编码。例如15′/16′=0余15′,56′/16′=3余8′,坐标(0,3)对应的16′×16′网格编码为5,所以西安钟楼所在的第3层级网格的编码为J2785。
步骤3-3-4:计算出该点位于第4层级网格的编码。计算公式为:
第4层级网格为4′×4′网格,故N3/4′的商为M4、余数为N4,M4即为该点所在的第4层级网格的纬度方向编号;n3/4′商为m4、余数为n4,m4即为该点所在的第4层级网格的经度方向编号。根据两个编号对应出第4层级网格编码。例如15′/4′=3余3′,8′/4′=2余0′,坐标(3,2)对应的4′×4′网格编码为E,所以西安钟楼所在的第4层级网格的编码为J2785E。
步骤3-3-5:计算出该点位于第5层级网格的编码。计算公式为:
第5层级网格为1′×1′网格,故N4/1′的商为M5,M5即为该点所在的第5层级网格的纬度方向编号;n4/1′商为m5、余数为n5,m5即为该点所在的第5层级网格的经度方向编号。根据两个编号对应出第5层级网格编码。例如3′/1′=3余0′,0′/1′=0余0′,坐标(3,0)对应的1′×1′网格编码为A,所以西安钟楼所在的第5层级网格的编码为J2785EA。
步骤3-3-6:计算出该点位于第6层级网格的编码。计算公式为:
第6层级网格为16″×16″网格,故LS/16″的商为M6、余数为N6,M6即为该点所在的第6层级网格的纬度方向编号;BS/16″商为m6、余数为n6,m6即为该点所在的第6层级网格的经度方向编号。根据两个编号对应出第6层级网格编码。例如39″/16″=2余7″,32″/16″=2余0″,坐标(2,2)对应的16″×16″网格编码为C,所以西安钟楼所在的第6层级网格的编码为J2785EAC。
步骤3-3-7:计算出该点位于第7层级网格的编码。计算公式为:
第7层级网格为4″×4″网格,故N6/4″的商为M7,M7即为该点所在的第7层级网格的纬度方向编号;n6/4″商为m7,m7即为该点所在的第7层级网格的经度方向编号。根据两个编号对应出第7层级网格编码。例如7″/4″=1余3″,0″/4″=0余0″,坐标(1,0)对应的4″×4″网格编码为2,所以西安钟楼所在的第7层级网格的编码为J2785EAC2。
步骤4:网格化空域。假设要进行水平范围冲突检测的空域共有K个,根据每个空域的类型和大小,用不同层级的多个网格去填充空域,提高空域表达的精度,进而提高冲突检测的准确度。
具体包括:
步骤4-1:计算空域的外包矩形。每个空域有m个角点P1~Pm(m的取值根据该空域实际角点个数确定),经纬度坐标分别为P1(L1,B1),P2(L2,B2)……,Pm(Lm,Bm),其中L为纬度,B为经度。由于各层级网格为正方形网格,方向正南正北、正东正西,故分别过纬度最小(Lmin)的点和纬度最大(Lmax)的点作两条正东正西的直线,距离记为d1;分别过经度最小(Bmin)的点和经度最大(Bmax)的点作两条正南正北的直线,距离记为d2。四条线的相交范围即空域的外包矩形。如图4所示,例如某空域有9个角点P1~P9,分别过纬度最小的点P5和纬度最大的点P1作两条正东正西的直线,距离记为d1;分别过经度最小的点P1和经度最大的点P4作两条正南正北的直线,距离记为d2。四条线的相交范围即该空域的外包矩形。采用上述方法,计算得到K个空域的外包矩形。
步骤4-2:计算空域的最小外包网格。如图4所示,计算出某空域外包矩形的面积S=d1×d2。第1~7层级网格都有固定的面积,比较出网格面积大于S且与S最接近的网格层级n。在空域的所有角点中,任意选择一点,按照步骤3-3的方法计算出该点所在的第n层级网格编码,该编码对应的网格即为最小外包网格。例如,若图4中d1=5千米,d2=7千米,则外包矩形面积为35平方千米。因为第4层级网格为64平方千米,大于35平方千米且与35平方千米最接近,所以该空域的最小外包网格对应的层级为4。假设任意选择的一个角点为P2,将P2的经纬度坐标按步骤3-3的方法转化为第4层级的网格编码K06E75,该编码对应的网格即为该空域的最小外包网格。采用上述方法,计算得到K个空域的最小外包网格。
步骤4-3:计算包含要进行水平范围冲突检测的K个空域的最小外包网格。这K个空域的最小外包网格的编码中前W位编码均相等,从第W+1位开始不相等,则任意选择一个最小外包网格的编码,取其前W位编码,即为所有K个空域的最小外包网格编码,该编码对应的网格即为K个空域的最小外包网格。
步骤4-4:计算空域边界所在的网格。如图5所示,将步骤4-3得到的最小外包网格剖分至指定层级N(n≤N≤7),以最小外包网格内第N层级网格中左下网格为原点网格,其坐标为(0,0),建立直角坐标系。向东、向北从0开始,按步长为1编号,直到编号到最小外包网格内最东和最北端的第N层级网格。确定空域边界每条线段端点所在网格的行列坐标,端点之间网格的行列坐标集由插值计算得到,并进行唯一性判断,去掉空域角点重复的坐标,得到最终的空域边界网格的行列坐标集合,其中列坐标在前,行坐标在后。例如,设图5中空域最小外包网格的编码为M24A,属于第2层级网格,将最小外包网格剖分到第4层级,以最小外包网格的左下点为原点,建立直角坐标系,向东、向北为正向,依次用01234…等数字编码。得到空域边界网格的行列坐标集合为{(2,0)、(2,1)、(2,2)、(2,3)、(2,4)、(2,5)、(2,6)、(2,7)、(2,8)、(2,9)、(2,10)、(3,0)、(4,0)、(5,0)、(6,0)、(7,1)、(7,2)、(7,3)、(7,4)、(7,5)、(7,6)、(7,7)、(7,8)、(7,9)、(7,10)、(3,11)、(4,11)、(5,11)、(6,11)}。
步骤4-5:计算空域内部所在的网格。将空域边界网格行列坐标集合中行坐标相等的网格分在同一组,该组中列坐标最小为Pmin,最大为Pmax,按步长为1填充Pmin~Pmax中的列坐标,行坐标不变。按照纬度由低到高、经度由西向东的顺序逐行进行扫描,直到空域内部填充完毕,获得空域内部网格的行列坐标集合,如图5(a)所示,Z字型曲线即为填充空域内部网格时的扫描路径。加上步骤4-4得到的边界坐标集,得到完整的空域所在网格的行列坐标集合。利用行列坐标与网格编码之间的对应关系,得到空域所在网格的编码集合。例如,在图5(a)中,行坐标为1的坐标为(2,1)和(7,1),按步长为1填充1-7中的列坐标,行坐标不变,得到该行中间的网格的坐标为(3,1)、(4,1)、(5,1)、(6,1)。依次类推,按照纬度由低到高、经度由西向东的顺序逐行进行扫描,直到空域内部填充完毕,获得空域内部网格的行列坐标集合,加上步骤4-4得到的边界坐标集,得到完整的空域所在网格的行列坐标集合为{(3,0)、(4,0)、(5,0)、(6,0)、(2,1)、(3,1)(4,1)、(5,1)、(6,1)、(7,1)、(2,2)、(3,2)(4,2)、(5,2)、(6,2)、(7,2)、(2,3)、(3,3)(4,3)、(5,3)、(6,3)、(7,3)、(2,4)、(3,4)(4,4)、(5,4)、(6,4)、(7,4)、(2,5)、(3,5)(4,5)、(5,5)、(6,5)、(7,5)、(2,6)、(3,6)(4,6)、(5,6)、(6,6)、(7,6)、(2,7)、(3,7)(4,7)、(5,7)、(6,7)、(7,7)、(2,8)、(3,8)(4,8)、(5,8)、(6,8)、(7,8)、(2,9)、(3,9)(4,9)、(5,9)、(6,9)、(7,9)、(2,10)、(3,10)(4,10)、(5,10)、(6,10)、(7,10)、(3,11)、(4,11)、(5,11)、(6,11)}。利用行列坐标与网格编码之间的对应关系,得到空域所在网格的编码集合为{M24A05、M24A10、M24A11、M24A14、M24A06、M24A07、M24A12、M24A13、M24A16、M24A17、M24A0C、M24A0D、M24A18、M24A19、M24A1C、M24A1D、M24A0E、M24A0F、M24A1A、M24A1B、M24A1E、M24A1F、M24A24、M24A25、M24A30、M24A31、M24A34、M24A35、M24A26、M24A27、M24A32、M24A33、M24A36、M24A37、M24A2C、M24A2D、M24A38、M24A39、M24A3C、M24A3D、M24A2E、M24A2F、M24A3A、M24A3B、M24A3E、M24A3F、M24A84、M24A85、M24A90、M24A91、M24A94、M24A95、M24A86、M24A87、M24A92、M24A93、M24A96、M24A97、M24A8C、M24A8D、M24A98、M24A99、M24A9C、M24A9D、M24A8F、M24A9A、M24A9B、M24A9C}。
步骤4-6:网格聚合。如图5(b)所示,将同属于1个父网格的16个子网格聚合成1个父网格,用1个父网格编码代替16个子网格编码。得到空域的最少网格编码集合。例如,在如图5(b)中,将同属于1个父网格的16个子网格聚合成1个父网格,用1个父网格编码代替16个子网格编码。得到空域的最少网格编码集合为{M24A05、M24A06、M24A07、M24A0C、M24A0D、M24A0E、M24A0F、M24A24、M24A25、M24A26、M24A27、M24A2C、M24A2D、M24A2E、M24A2F、M24A84、M24A85、M24A86、M24A87、M24A8C、M24A8D、M24A8F、M24A1、M24A3、M24A9}。
步骤5:判断网格编码集合间有无交集。此步的目的是筛选出必然存在冲突的空域和可能存在冲突的空域,达到减小计算量、提高实时性的目的。当空域对应的网格编码集合间有交集,即两空域包含同一个网格时,空域在水平范围内必然存在冲突,不再进行下一步冲突判定;只有空域对应的网格编码集合间无交集,即两空域不包含同一个网格时,空域在水平范围内才有可能无冲突,转入步骤6。如图7所示,对于空域A和空域D来说,编码集合存在交集,不再计算空域间水平距离。空域A与空域B、空域A与空域C、空域B与空域C、空域B与空域D、空域C与空域D的编码集合无交集,两空域不包含同一个网格,空域在水平范围内才有可能无冲突,转入步骤6,计算空域间水平距离。
步骤6:判断空域间水平距离。而对于步骤5中网格编码集合间无交集的空域对,才需要利用网格编码计算空域间距离。将计算出的距离与《间隔规定》和《中华人民共和国飞行基本规则》规定的安全间隔对比,若大于安全间隔,则判定为空域在水平范围上不存在冲突;否则,空域之间在水平范围上存在冲突。
具体包括以下步骤:
步骤6-1:判定空域间相对方位。如图6所示,即根据两空域的中心点经纬度坐标,判定两空域的方位关系。判定相对方位的方法如下:设空域1中心点的经纬度坐标为(L1,B1),空域2中心点的经纬度坐标为(L2,B2),其中L表示纬度,其中北纬为正,南纬为负;B表示经度,其中东经为正,西经为负。若满足L1<L2且B1<B2,则空域2位于空域1的东北方向;若满足L1>L2且B1>B2,则空域2位于空域1的西南方向;若满足L1<L2且B1>B2,则空域2位于空域1的西北方向;若满足L1>L2且B1<B2,则空域2位于空域1的东南方向;若满足L1<L2且B1=B2,则空域2位于空域1的正北方向;若满足L1>L2且B1=B2,则空域2位于空域1的正南方向;若满足L1=L2且B1>B2,则空域2位于空域1的正西方向;若满足L1=L2且B1<B2,则空域2位于空域1的正东方向。在图7中,
步骤6-2:确定两空域在相对方位上最边缘网格的行列坐标。如图7所示,根据表达精度对应的最小剖分层级N,确定两空域在相对方位上最边缘的第N级网格编码,取最后两位编码转化为一维二进制编码,再拆分为二维二进制编码,最后转化为二维十进制编码。以最小外包网格的左下角为坐标原点(0,0),得到两空域在相对方位上最边缘的第N级网格的坐标。例如,在图7中,采用步骤6-1的判定方法,可判定空域C位于空域A的东北反向。4个空域的最小外包网格为第3层级网格,所对应的编码为J2764,根据表达精度需要确定最小剖分层级为5级,故确定空域A东北角的第5级网格编码为J27643F,空域C西南角的第5级网格编码为J2764D8。将3F转化为一维二进制编码为00111111,进而拆分为二维二进制编码为0111和0111,并转化为二维十进制编码7和7,即空域A东北角第5级网格的行列坐标为(7,7);将D8转化为一维二进制编码为11011000,进而拆分为二维二进制编码为1010和1100,并转化为二维十进制编码10和12,即空域C西南角第5级网格的行列坐标为(12,10)。
步骤6-3:通过坐标计算出空域间水平距离。如图7所示,通过坐标计算两空域在相对方位上最边缘第N级网格之间的横向间距与纵向间距,使用勾股定理求得两空域间水平距离,若水平距离小于水平安全间隔(包括横向安全间隔、纵向安全间隔),则两空域在水平范围存在冲突,判定两空域空间上有冲突;否则,两空域在水平范围不冲突,判定两空域在空间维上不冲突。水平安全间隔为10公里。例如,图7中,通过坐标计算空域A和空域C在相对方位第5级网格之间的横向间距与纵向间距。12-7-1=4,10-7-1=2,即空域A和空域C横向之间由4个、纵向有2个第5级网格,故横向间距为8km,纵向间距为4km。勾股定理得到空域A和空域C的水平距离为8.9公里,而水平安全间隔为10公里,小于水平安全间隔,故空域A和空域C有冲突。
通过构建地球表面剖分网格体系,并对网格进行编码,每个网格都对应唯一的编码,对空域进行网格化表达。将空域表达为网格集合,由网格与编码的唯一对应关系,得到空域的编码集合,通过编码集合求交集,对无交集的空域对,通过相对方位上最边缘网格的编码,计算出空域间水平间隔,进行水平范围上空域冲突检测。
步骤7:对时间维有冲突而空间维无冲突的所有空域,进行频率维冲突检测。
具体包括:
步骤7-1:统一频率维的编码格式。频率单位有Hz、KHz、MHz、GHz、THz,对空通信频率都在KHz以上,故将频率单位统一为KHz,单位前面的数值统一设定为9位,不足9位的在数值前面补0。例如某空域1的使用频率为5MHz,则编码为000005000,空域2的使用频率为225KHz,则编码为000000225。
步骤7-2:进行编码运算。两空域所使用的频率值先逐位进行异或运算,异或运算的结果为0或1,方法为当数值a与b相等时,当数值a与b不相等时,/>其中符号/>代表异或运算。再对异或运算得到的9个数值逐位进行求并运算,公式为1∪0=1,1∪1=1,0∪1=1,0∪0=0,其中符号∪代表求并运算,目的是检查逐位异或运算得到的9个数值是否均为0。例如步骤7-1中两空域编码分别为000005000和000000225,逐位进行异或运算的结果为000001111,再逐位进行求并运算,第1位0与第二位0求并运算结果为0,该预算结果0与第三位0求并运算结果为0,依次类推,得到最终的计算结果为1。
步骤7-3:根据步骤7-2的运算结果,进行冲突判定。若最终求并运算的结果为0,则两空域所使用的频率相等,频率之间会存在相互干扰,空域间存在频率冲突;若结果为1,则两空域所使用的频率不相等,频率之间不存在相互干扰,空域间不存在频率冲突。
图8为本发明方法程序实现的详细流程图。结合实施例,具体说明本发明的具体实施步骤。如图7所示,某区域内要进行冲突检测的空域共4个,分别是空域A、空域B、空域C、空域D。其中,空域A为一个16km×16km的正方形空域,使用时间范围为2020年10月17日早上6点到2020年11月10日早上6点,高度范围为标压高5000米到10000米,使用频率为5MHz。空域B为半径为8km的圆形空域,使用时间范围为2020年10月20日早上8点到2020年10月30日早上8点,高度范围为标压高6000米到8000米,使用频率为225KHz。空域C为长边为16公里,短边为4km的跑道性空域,使用时间范围为2020年10月25日早上6点到2020年11月5日下午8点,高度范围为标压高6000米到7500米,使用频率为500KHz。空域D为长16km×8km的矩形空域,使用时间范围为2020年10月1日早上6点到2020年10月6日早上6点,高度范围为标压高5000米到10000米,使用频率为3MHz。
以上4个空域按照本发明方法进行冲突检测,首先进行步骤1,构建空域的时间维编码,进行时间维冲突检测。具体方法为:
步骤1-1:时间维编码。将空域A、B、C、D的使用时间范围分别编码为10170600-11100600、10200800-10300800、10250600-11052000、10010600-10060600。
步骤1-2:进行时间维冲突检测。两两空域之间进行时间维冲突检测,以空域A和空域D为例,空域A的前8位编码为10170600,先与空域B的后8位编码10060600进行逐位比较,最终比较结果Result-1为0;再进行11100600与10010600之间的运算,最终比较结果Result-2为0。Result-1和Result-2进行异或运算,运算结果为0,则空域A和D的时间范围不重叠,在时间维上无冲突,直接判定空域A和D无冲突,不再进行下一步的冲突检测。按照同样的方法,判定出空域A、B、C之间在时间维上存在冲突。
对步骤1中存在时间维冲突的空域A、B、C,进行步骤2,构建空域的高度范围编码,进行高度范围冲突检测。具体方法为:
步骤2-1:高度范围编码。空域A、B、C的高度范围编码分别为H050H100、H060H080、H060H075。
步骤2-2:进行高度范围冲突检测。两两空域之间进行高度范围冲突检测,以空域A和空域B为例,空域A的高度范围编码的前4位为H050,先与空域B的后4位编码H080进行逐位比较,最终比较结果为Result-3为1,再进行H100与H060之间的逐位比较,最终比较结果为Result-4为0,Result-3和Result-4进行异或运算,运算结果为1,则空域A和B的高度范围重叠,两空域在高度维上有冲突。按照同样的方法,判定出空域A、B、C之间在高度范围上存在冲突。转入步骤3。
由于空域A、B、C之间在高度范围上也存在冲突,故需要对空域A、B、C进行水平范围冲突检测。首先,进行步骤3,建立地球表面剖分网格体系。具体方法为:
步骤3-1:进行第1层级剖分与编码。如图2所示,选择本初子午线与赤道交点为剖分原点,用4°×4°的第一层级网格将纬度范围为[-90°,90°]、经度范围为[-180°,180°]的地球表面剖分为46×90份。第一层级网格序号在纬度方向上用A~Y和a~y共46个英文字母从赤道按照纬度由低到高编码,为避免字母与数字1和0混淆,编码跳过I、O、i和o,其中北纬大写、南纬小写;第一层级网格序号在经度方向上从本初子午线开始由西向东用00~89共90个数字编码。第1层级网格的边长为512千米。
步骤3-2:进行第2~7层级剖分与编码。上一层级的网格称为父网格,下一层级的网格称为子网格。每1个上一层级父网格平均分成16份,得到下一层级16个子网格,对16个子网格按照皮亚诺空间填充曲线进行16进制编码,直到剖分至第7层级。所以子网格编码继承上一层级父网格的编码,并比上一层级父网格编码多1位。第2~7级网格的大小分别为1°×1°、16′×16′、4′×4′、1′×1′、16″×16″、4″×4″,网格边长分别为128千米、32千米、8千米、2千米、512米、128米。例如,图7中,第3层级某16′×16′网格的编码为J2764,将其平均分成16个4′×4′的第4层级网格,对这16个子网格按照图2中的编码方式进行编码,编码分别为{J27640,J27641,J27642,J27643,J27644,J27645,J27646,J27647,J27648,J27649,J2764A,J2764B,J2764C,J2764D,J2764E,J2764F}。经过步骤2-1和步骤2-2,将地球剖分成无缝隙、无重叠的网格,每个网格都有唯一的编码以及相对应的面积。
步骤3-3:转换经纬度坐标与网格编码。西安市钟楼的经纬度坐标为(34°15′39″N,108°56′32″E),按照步骤3-3-1到步骤3-3-7的方法,将该经纬度坐标转化为网格编码作为步骤3-3的一个实施例。34°/4°=8余2°,108°/4°=27,位于第1层级中北纬向第9个、经向第28个4°×4°网格,所以西安钟楼所在的第1层级网格的编码为J27。2°/1°=2,0/1°=0,坐标(2,0)对应的1°×1°网格编号为8,所以西安钟楼所在的第2层级网格的编码为J278。15′/16′=0余15′,56′/16′=3余8′,坐标(0,3)对应的16′×16′网格编码为5,所以西安钟楼所在的第3层级网格的编码为J2785。15′/4′=3余3′,8′/4′=2余0′,坐标(3,2)对应的4′×4′网格编码为E,所以西安钟楼所在的第4层级网格的编码为J2785E。3′/1′=3余0′,0′/1′=0余0′,坐标(3,0)对应的1′×1′网格编码为A,所以西安钟楼所在的第5层级网格的编码为J2785EA。39″/16″=2余7″,32″/16″=2余0″,坐标(2,2)对应的16″×16″网格编码为C,所以西安钟楼所在的第6层级网格的编码为J2785EAC。7″/4″=1余3″,0″/4″=0余0″,坐标(1,0)对应的4″×4″网格编码为2,所以西安钟楼所在的第7层级网格的编码为J2785EAC2。
建立地球表面剖分网格体系之后,要对空域A、B、C按照步骤4进行空域网格化表达,包括:
步骤4-1:计算空域的外包矩形。在图7中,空域A的外包矩形为空域A本身。空域B的外包矩形为8km×8km的正方形,空域C的外包矩形为8km×2km的矩形。
步骤4-2:计算空域的最小外包网格。空域A的最小外包网格编码为J2764,空域B的最小外包网格编码为J27647,空域C的最小外包网格编码为J2764。
步骤4-3:计算包含空域A、B、C的最小外包网格。空域A、B、C的最小外包网格编码前5位均相等,因此包含空域A、B、C的最小外包网格的编码为J2764。
步骤4-4:计算空域边界所在的网格。如图7所示,根据表达精度的需要,将步骤4-3得到的最小外包网格剖分至第5层级,以左下角的第5层级网格坐标为(0,0),建立直角坐标系,并向东、向北依次编码。空域A的边界所在网格的行列坐标为{(0,0)、(1,0)、(2,0)、(3,0)、(4,0)、(5,0)、(6,0)、(7,0)、(0,1)、(0,2)、(0,3)、(0,4)、(0,5)、(0,6)、(0,7)、(1,7)、(2,7)、(3,7)、(4,7)、(5,7)、(6,7)、(7,7)、(7,1)、(7,2)、(7,3)、(7,4)、(7,5)、(7,6)},空域B的边界所在网格的行列坐标为{(12,4)、(13,4)、(14,4)、(15,4)、(12,5)、(12,6)、(12,7)、(13,7)、(14,7)、(15,7)、(15,5)、(15,6)},空域C的边界所在网格的行列坐标为{(12,10)、(12,14)}。
步骤4-5:计算空域内部所在的网格。经过逐行填充后,获得空域A、B、C内部网格的行列坐标集合,加上步骤4-4得到的边界坐标集,得到完整的空域所在网格的行列坐标集合。利用行列坐标与网格编码之间的对应关系,得到空域所在网格的编码集合。其中空域A所在的第5层级网格编码集合为{J276400、J276401、J276402、J276403、J276404、J276405、J276406、J276407、J276408、J276409、J27640A、J27640B、J27640C、J27640D、J27640E、J27640F、J276410、J276411、J276412、J276413、J276414、J276415、J276416、J276417、J276418、J276419、J27641A、J27641B、J27641C、J27641D、J27641E、J27641F、J276420、J276421、J276422、J276423、J276424、J276425、J276426、J276427、J276428、J276429、J27642A、J27642B、J27642C、J27642D、J27642E、J27642F、J276430、J276431、J276432、J276433、J276434、J276435、J276436、J276437、J276438、J276439、J27643A、J27643B、J27643C、J27643D、J27643E、J27643F},空域B所在的第5层级网格编码集合为{J276470、J276471、J276472、J276473、J276474、J276475、J276476、J276477、J276478、J276479、J27647A、J27647B、J27647C、J27647D、J27647E、J27647F},空域C所在的第5层级网格编码集合为{J2764D8、J2764DA、J2764F0、J2764F2、J2764F8}。
步骤4-6:网格聚合。将同属于1个父网格的16个子网格聚合成1个父网格,用1个父网格编码代替16个子网格编码。得到空域的最少网格编码集合。聚合后,空域A的最少网格编码集合为{J27640、J27641、J27642、J27643},空域B的最少网格编码集合为{J27647},空域C的最少网格编码集合为{J2764D8、J2764DA、J2764F0、J2764F2、J2764F8},
空域A、B、C经过步骤4的网格化表达后,进入步骤5,计算网格编码集合间有无交集。空域A与空域B、空域A与空域C、空域B与空域C的编码集合无交集,说明上述两空域间不包含同一个网格,空域在水平范围内才有可能无冲突,转入步骤6。
在步骤6中,计算空域间水平范围的距离,进行水平范围冲突判定。方法包括以下步骤:
步骤6-1:判定方位。在图6中,可判定空域B位于空域A的东北方向,空域C位于空域A的东北反向,空域B位于空域C的东南方向。
步骤6-2:确定两空域在相对方位上最边缘网格的行列坐标。最小剖分层级为5级,由于空域C位于空域A的东北反向,故确定空域A东北角的第5级网格编码为J27643F,空域C西南角的第5级网格编码为J2764D8。将3F转化为一维二进制编码为00111111,进而拆分为二维二进制编码为0111和0111,即空域1东北角第5级网格的行列坐标为(7,7);将D8转化为一维二进制编码为11011000,进而拆分为二维二进制编码为1010和1100,即空域3西南角第5级网格的行列坐标为(12,10)。
步骤6-3:通过坐标计算出空域间水平距离。通过坐标计算两空域相对方位第5级网格之间的横向间距与纵向间距。12-7-1=4,10-7-1=2,即空域A和空域C横向之间由4个、纵向有2个第5级网格,故横向间距为8km,纵向间距为4km。勾股定理得到空域A和空域C的水平距离为8.9公里,而水平安全间隔为10公里,小于水平安全间隔,故空域A和空域C在水平范围内有冲突。同理可得空域A和空域B、空域B和空域C之间的水平距离小于水平安全间隔,空域A和空域B、空域B和空域C之间在水平范围内也有冲突。加上步骤1判定的空域A、B、C在时间维上存在冲突,步骤2判定的空域A、B、C在高度范围上存在冲突,故最终判定为空域A、B、C之间存在冲突,冲突检测结束。
图7所代表的实施例中,4个空域在进行频率冲突检测前,就已经最终判定空域间存在冲突,不需要再进行步骤7。在本发明的另一个实施例中,空域E的使用频率维500KHz,空域F的使用频率维850KHz,经过步骤1-6,空域E和F之间在时间维有冲突,而在空间维无冲突,故空域E和F之间需要进行步骤7——频率维冲突检测。结合该实施例说明步骤7所用的方法。
步骤7-1:统一频率维的编码格式。空域E的使用频率编码为000000500,空域F的使用频率编码为000000850。
步骤7-2:进行编码运算。逐位进行异或运算的结果为000000110,再逐位进行求并运算,第1位0与第二位0求并运算结果为0,该预算结果0与第三位0求并运算结果为0,依次类推,得到最终的计算结果为1。
步骤7-3:根据最终的运算结果,进行冲突判定。最终的计算结果为1,则空域E和F所使用的频率不相等,频率之间不存在相互干扰,空域E和F间不存在频率冲突。最终判定为空域E和F之间无冲突,冲突检测结束。
本发明方法实现了空域在时间维、空间维和频率维的一体化冲突检测,特别是对于水平范围冲突检测,通过构建7层网格体系和对应的编码方案,将各类空域用网格表达,用不同层级的多个网格而不是用特定层级的一个网格去表达空域,提高空域表达的精度,进而提高冲突检测准确度;通过编码集合求交集,筛选出必然存在冲突的空域和可能存在冲突的空域,减小计算量,提高了实时性;对于可能存在冲突的空域,利用网格编码进行距离计算,与水平安全间隔对比,进而判定是否存在水平范围冲突;本发明方法的计算过程都是基于网格编码,编码之间采用比较、异或、求并运算,而非传统的基于经纬度坐标的浮点计算,更符合计算机的计算逻辑,提高计算机的运算效率。

Claims (5)

1.一种适用于全球空间范围的空域冲突检测方法,包括:
步骤1:构建空域时间维编码,检测时间维冲突;
步骤2:构建空域高度范围编码,检测高度范围冲突;
步骤3:建立地球表面剖分网格体系;
步骤4:网格化空域;
步骤5:判断网格编码集合间有无交集;
步骤6:判断空域间水平距离;
步骤7:进行频率维冲突检测;
步骤1根据空域的使用时间范围,构建格式统一的时间维编码,便于计算机识别和读取;若空域的时间维无冲突,直接判定为空域之间无冲突,不需要进行下面步骤的判定;若空域的时间维有冲突,则转入步骤2;步骤2对时间维存在冲突的空域,根据空域的高度范围,构建格式统一的高度范围编码;若空域的高度范围无冲突,则转入步骤7;若空域的高度范围有冲突,则转入步骤3;在步骤3中进行第1至7层级剖分与编码;而后转换经纬度坐标与网格编码;步骤4中根据每个空域的类型和大小,用不同层级的多个网格去填充空域;步骤5筛选出必然存在冲突的空域和可能存在冲突的空域;步骤6中将计算出的空域间距离与规定的安全间隔对比,若大于安全间隔,则判定为空域在水平范围上不存在冲突;否则,空域之间在水平范围上存在冲突;步骤7中对时间维有冲突而空间维无冲突的所有空域,进行频率维冲突检测;
步骤1包括:
步骤1-1:时间维编码;
步骤1-2:检测时间维冲突;
步骤2包括:
步骤2-1:高度范围编码;
步骤2-2:检测高度范围冲突;
步骤3包括:
步骤3-1:进行第1层级剖分与编码;
步骤3-2:进行第2~7层级剖分与编码;
步骤3-3:转换经纬度坐标与网格编码;
步骤4包括:
步骤4-1:计算空域的外包矩形;
步骤4-2:计算空域的最小外包网格;
步骤4-3:计算包含要进行水平范围冲突检测的多个空域的最小外包网格;
步骤4-4:计算空域边界所在的网格;
步骤4-5:计算空域内部所在的网格;
步骤4-6:网格聚合;
步骤6包括:
步骤6-1:判定空域间相对方位;
步骤6-2:确定两空域在相对方位上最边缘网格的行列坐标;
步骤6-3:通过坐标计算出空域间水平距离;
步骤7包括:
步骤7-1:统一频率维的编码格式;
步骤7-2:进行编码运算;
步骤7-3:根据步骤7-2的运算结果,进行冲突判定;
步骤1-1中时间编码为一串16位的数字,前8位数字代表为空域的启用时间,后8位数字代表空域的结束使用时间,中间用“-”隔开;步骤1-2中两两空域之间进行时间维冲突检测,两个空域只有时间范围不重叠,才无冲突;步骤2-1中高度编码用大写英文字母H加上三位数字表示,数字的单位为百米;步骤2-2中两两空域之间进行高度范围冲突检测,两个空域只有高度范围不重叠,才无冲突;步骤3-1中选择本初子午线与赤道交点为剖分原点,用4°×4°的第一层级网格将纬度范围为[-90°,90°]、经度范围为[-180°,180°]的地球表面剖分为46×90份;步骤3-2中上一层级的网格为父网格,下一层级的网格为子网格;每1个上一层级父网格平均分成16份,得到下一层级16个子网格,对16个子网格按照皮亚诺空间填充曲线进行16进制编码,直到剖分至第7层级;步骤4-1中分别过纬度最小的点和纬度最大的点作两条正东正西的直线;分别过经度最小的点和经度最大的点作两条正南正北的直线,四条线的相交范围即该空域的外包矩形;步骤4-2中比较出网格面积大于步骤4-1中空域外包矩形的面积且与该面积最接近的网格层级;
步骤4-3中选择一个最小外包网格的编码,取其编码均相等前多位编码,编码对应的网格即为多个空域的最小外包网格;步骤4-4中将步骤4-3得到的最小外包网格剖分至指定层级;步骤4-5中将空域边界网格行列坐标集合中行坐标相等的网格分在同一组,该组中列坐标最小为Pmin,最大为Pmax,按步长为1填充Pmin~Pmax中的列坐标,行坐标不变,按照纬度由低到高、经度由西向东的顺序逐行进行扫描,直到空域内部填充完毕,获得空域内部网格的行列坐标集合;步骤4-6中将同属于1个父网格的16个子网格聚合成1个父网格,用1个父网格编码代替16个子网格编码;步骤5中当空域对应的网格编码集合间有交集,即两空域包含同一个网格时,空域在水平范围内必然存在冲突,不再进行下一步冲突判定;只有空域对应的网格编码集合间无交集,即两空域不包含同一个网格时,空域在水平范围内才有可能无冲突,转入步骤6;步骤6中将计算出的距离与规定的安全间隔对比,若大于安全间隔,则判定为空域在水平范围上不存在冲突;否则,空域之间在水平范围上存在冲突;
步骤1-1中在前8位和后8位数字中,第1-2位数字代表某月,第3-4位数字代表某日,第5-8位数字代表某时某分,采用24小时制,空域的时间维编码默认省略当年年份;若空域使用时间跨年,则在空域启用时间和结束使用时间编码前面各加上4位数字,代表年份编码;步骤1-2中在两个空域中任选一个空域,记为空域A,取空域A的前8位编码,先与另外一个空域B的后8位编码进行逐位比较,最终比较结果为Result-1,再取空域A的后8位编码与空域B的前8位编码进行逐位比较,最终比较结果为Result-2,对Result-1和Result-2进行异或运算,判断时间维是否冲突;步骤2-1中高度范围的编码共8位,其中前4位为空域的下限高度,后4位为空域的上限高度;步骤2-2中在时间维上存在冲突的两个空域中任选一个空域,记为空域M,取空域M高度范围编码的前4位,先与另外一个空域N的后4位编码进行逐位比较,最终比较结果为Result-3;再取空域A的后4位编码与空域N的前4位编码进行逐位比较,最终比较结果为Result-4,对Result-3和Result-4进行异或运算,判断高度范围是否冲突;步骤3-1中第一层级网格序号在纬度方向上用A~Y和a~y共48个英文字母从赤道按照纬度由低到高编码,编码跳过I、O、i和o,其中北纬大写、南纬小写;第一层级网格序号在经度方向上从本初子午线开始由西向东用00~89共90个数字编码;步骤3-3中某点经纬度坐标为(L,B),其中L为纬度,B为经度,L和B的表示形式均为度分秒,记作L=LD°LM′LS″,B=BD°BM′BS″;以步骤4-4中以最小外包网格内第N层级网格中左下网格为原点网格,其坐标为(0,0),建立直角坐标系,向东、向北从0开始,按步长为1编号,直到编号到最小外包网格内最东和最北端的第N层级网格,确定空域边界每条线段端点所在网格的行列坐标,端点之间网格的行列坐标集由插值计算得到,去掉空域角点重复的坐标,得到空域边界网格的行列坐标集合,其中列坐标在前,行坐标在后;步骤7-1中频率单位统一为KHz,单位前面的数值统一设定为9位,不足9位的在数值前面补0;步骤7-2中对两空域所使用的频率值先逐位进行异或运算,再对异或运算得到的9个数值逐位进行求并运算;步骤7-3中若最终求并运算的结果为0,则两空域所使用的频率相等,频率之间会存在相互干扰,空域间存在频率冲突;若结果为1,则两空域所使用的频率不相等,频率之间不存在相互干扰,空域间不存在频率冲突。
2.如权利要求1所述的一种适用于全球空间范围的空域冲突检测方法,其特征在于:步骤3-1中第1层级网格的边长为512千米;步骤3-2中第2~7级网格的大小分别为1°×1°、16′×16′、4′×4′、1′×1′、16″×16″、4″×4″,网格边长分别为128千米、32千米、8千米、2千米、512米、128米。
3.如权利要求2所述的一种适用于全球空间范围的空域冲突检测方法,其特征在于:步骤3-3中转换经纬度坐标与网格编码具体包括:
步骤3-3-1:计算第1层级网格编码;计算公式为:
LD/4°的商为M1、余数为N1,M1+1对应的英文字母为该点所在的第1层级网格的纬度方向编号,北纬大写,南纬小写;BD/4°商为m1、余数为n1,m1即为该点所在的第1层级网格的经度方向编号;两个编号组合成该点所在的第1层级网格编码;
步骤3-3-2:计算第2层级网格编码;计算公式为:
步骤1的余数N1除以1°得到数值M2,M2即为该点所在的第2层级网格的纬度方向编号;n1除以1°得到数值m2,m2即为该点所在的第2层级网格的经度方向编号,根据两个编号对应得出第2层级网格编码;
步骤3-3-3:计算第3层级网格编码,计算公式为:
LM/16′的商为M3、余数为N3,M3即为该点所在的第3层级网格的纬度方向编号;BM/16′商为m3、余数为n3,m3即为该点所在的第3层级网格的经度方向编号,根据两个编号对应出第3层级网格编码;
步骤3-3-4:计算出该点位于第4层级网格的编码,计算公式为:
N3/4′的商为M4、余数为N4,M4即为该点所在的第4层级网格的纬度方向编号;n3/4′商为m4、余数为n4,m4即为该点所在的第4层级网格的经度方向编号,根据两个编号对应出第4层级网格编码;
步骤3-3-5:计算出该点位于第5层级网格的编码,计算公式为:
N4/1′的商为M5,M5即为该点所在的第5层级网格的纬度方向编号;n4/1′商为m5、余数为n5,m5即为该点所在的第5层级网格的经度方向编号,根据两个编号对应出第5层级网格编码;
步骤3-3-6:计算出该点位于第6层级网格的编码,计算公式为:
LS/16″的商为M6、余数为N6,M6即为该点所在的第6层级网格的纬度方向编号;BS/16″商为m6、余数为n6,m6即为该点所在的第6层级网格的经度方向编号,根据两个编号对应出第6层级网格编码;
步骤3-3-7:计算出该点位于第7层级网格的编码,计算公式为:
N6/4″的商为M7,M7即为该点所在的第7层级网格的纬度方向编号;n6/4″商为m7,m7即为该点所在的第7层级网格的经度方向编号,根据两个编号对应出第7层级网格编码。
4.如权利要求2所述的一种适用于全球空间范围的空域冲突检测方法,其特征在于:
步骤6-1中根据两空域的中心点经纬度坐标,判定两空域的方位关系;步骤6-2中根据表达精度对应的最小剖分层级N,确定两空域在相对方位上最边缘的第N级网格编码,取最后两位编码转化为一维二进制编码,再拆分为二维二进制编码,最后转化为二维十进制编码,以最小外包网格的左下角为坐标原点(0,0),得到两空域在相对方位上最边缘的第N级网格的坐标;步骤6-3中通过坐标计算两空域在相对方位上最边缘第N级网格之间的横向间距与纵向间距,使用勾股定理求得两空域间水平距离,若水平距离小于水平安全间隔,则两空域在水平范围存在冲突,判定两空域空间上有冲突;否则,两空域在水平范围不冲突,判定两空域在空间维上不冲突。
5.如权利要求4所述的一种适用于全球空间范围的空域冲突检测方法,其特征在于:步骤6-3中水平安全间隔包括横向安全间隔、纵向安全间隔;水平安全间隔为10公里。
CN202011354731.8A 2020-11-27 2020-11-27 一种适用于全球空间范围的空域冲突检测方法 Active CN112348960B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011354731.8A CN112348960B (zh) 2020-11-27 2020-11-27 一种适用于全球空间范围的空域冲突检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011354731.8A CN112348960B (zh) 2020-11-27 2020-11-27 一种适用于全球空间范围的空域冲突检测方法

Publications (2)

Publication Number Publication Date
CN112348960A CN112348960A (zh) 2021-02-09
CN112348960B true CN112348960B (zh) 2024-05-07

Family

ID=74365010

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011354731.8A Active CN112348960B (zh) 2020-11-27 2020-11-27 一种适用于全球空间范围的空域冲突检测方法

Country Status (1)

Country Link
CN (1) CN112348960B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113051633B (zh) * 2021-04-02 2024-03-22 中国人民解放军空军工程大学 一种航空器运行可视化方法
CN113505189B (zh) * 2021-09-07 2022-03-04 中国电子科技集团公司第二十八研究所 一种基于多层级全球格网的面状目标栅格化方法
CN113781846B (zh) * 2021-09-28 2022-12-16 中国人民解放军空军工程大学 一种水平维度空域冲突检测方法
CN114419281B (zh) * 2022-01-18 2022-12-13 中国人民解放军93209部队 一种计算空域网格的空间几何关系的方法
WO2024007256A1 (zh) * 2022-07-07 2024-01-11 南京航空航天大学 空域数字化栅格的无人机冲突探测方法、装置和存储介质

Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4970682A (en) * 1981-01-13 1990-11-13 Harris Corporation Digital map generator and display system
US5636123A (en) * 1994-07-15 1997-06-03 Rich; Richard S. Traffic alert and collision avoidance coding system
CN102609525A (zh) * 2012-02-10 2012-07-25 北京大学 一种统一现有经纬度剖分网格的方法
CN103108189A (zh) * 2013-03-01 2013-05-15 北京大学 一种用于北斗遥感影像数据传输方法
CN103136371A (zh) * 2013-03-21 2013-06-05 北京大学 多源空间数据的剖分标识生成方法以及数据检索方法
CN103150338A (zh) * 2013-02-04 2013-06-12 北京大学 一种用于影像数据提取的剖分预处理方法及数据提取方法
CN104050374A (zh) * 2014-06-19 2014-09-17 中国电子科技集团公司第二十八研究所 一种空域冲突的快速检测方法
US9307360B1 (en) * 2015-01-09 2016-04-05 NinthDecimal, Inc. Systems and methods to identify a predefined geographical region in which a mobile device is located
CN105894862A (zh) * 2016-05-05 2016-08-24 中国民用航空华东地区空中交通管理局 一种空中交通管制智能化指挥系统
CN107808552A (zh) * 2017-10-27 2018-03-16 中航信移动科技有限公司 航班行为检测方法和装置
CN109559567A (zh) * 2018-12-12 2019-04-02 四川九洲空管科技有限责任公司 一种双层网格化空域资源动态状态管理方法及系统
CN110532342A (zh) * 2019-08-16 2019-12-03 南京王师大数据有限公司 全球空间目标区域位置编码方法及装置
CN111210669A (zh) * 2020-01-11 2020-05-29 中国人民解放军93209部队 一种基于推演评估的空域冲突解脱规则创建方法及系统
CN111429757A (zh) * 2020-03-09 2020-07-17 中国电子科技集团公司第十五研究所 一种空域使用冲突自动检测方法及系统
CN111427978A (zh) * 2019-11-22 2020-07-17 中国电子科技集团公司电子科学研究院 基于城市形态学和r树的三维模型索引方法
CN111475597A (zh) * 2020-03-31 2020-07-31 中国人民解放军战略支援部队信息工程大学 非刚性网格编码、空间对象唯一标识、查询方法及装置
CN111477034A (zh) * 2020-03-16 2020-07-31 中国电子科技集团公司第二十八研究所 基于栅格模型的大规模空域使用计划冲突检测与解脱方法
KR102152879B1 (ko) * 2020-01-17 2020-09-07 주식회사 씨엠월드 3차원 격자를 이용한 공간 정보 저장 방법 및 그 장치
CN111898819A (zh) * 2020-07-24 2020-11-06 盈嘉互联(北京)科技有限公司 空间网格划分方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8606491B2 (en) * 2011-02-22 2013-12-10 General Electric Company Methods and systems for managing air traffic

Patent Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4970682A (en) * 1981-01-13 1990-11-13 Harris Corporation Digital map generator and display system
US5636123A (en) * 1994-07-15 1997-06-03 Rich; Richard S. Traffic alert and collision avoidance coding system
CN102609525A (zh) * 2012-02-10 2012-07-25 北京大学 一种统一现有经纬度剖分网格的方法
CN103150338A (zh) * 2013-02-04 2013-06-12 北京大学 一种用于影像数据提取的剖分预处理方法及数据提取方法
CN103108189A (zh) * 2013-03-01 2013-05-15 北京大学 一种用于北斗遥感影像数据传输方法
CN103136371A (zh) * 2013-03-21 2013-06-05 北京大学 多源空间数据的剖分标识生成方法以及数据检索方法
CN104050374A (zh) * 2014-06-19 2014-09-17 中国电子科技集团公司第二十八研究所 一种空域冲突的快速检测方法
US9307360B1 (en) * 2015-01-09 2016-04-05 NinthDecimal, Inc. Systems and methods to identify a predefined geographical region in which a mobile device is located
CN105894862A (zh) * 2016-05-05 2016-08-24 中国民用航空华东地区空中交通管理局 一种空中交通管制智能化指挥系统
CN107808552A (zh) * 2017-10-27 2018-03-16 中航信移动科技有限公司 航班行为检测方法和装置
CN109559567A (zh) * 2018-12-12 2019-04-02 四川九洲空管科技有限责任公司 一种双层网格化空域资源动态状态管理方法及系统
CN110532342A (zh) * 2019-08-16 2019-12-03 南京王师大数据有限公司 全球空间目标区域位置编码方法及装置
CN111427978A (zh) * 2019-11-22 2020-07-17 中国电子科技集团公司电子科学研究院 基于城市形态学和r树的三维模型索引方法
CN111210669A (zh) * 2020-01-11 2020-05-29 中国人民解放军93209部队 一种基于推演评估的空域冲突解脱规则创建方法及系统
KR102152879B1 (ko) * 2020-01-17 2020-09-07 주식회사 씨엠월드 3차원 격자를 이용한 공간 정보 저장 방법 및 그 장치
CN111429757A (zh) * 2020-03-09 2020-07-17 中国电子科技集团公司第十五研究所 一种空域使用冲突自动检测方法及系统
CN111477034A (zh) * 2020-03-16 2020-07-31 中国电子科技集团公司第二十八研究所 基于栅格模型的大规模空域使用计划冲突检测与解脱方法
CN111475597A (zh) * 2020-03-31 2020-07-31 中国人民解放军战略支援部队信息工程大学 非刚性网格编码、空间对象唯一标识、查询方法及装置
CN111898819A (zh) * 2020-07-24 2020-11-06 盈嘉互联(北京)科技有限公司 空间网格划分方法及装置

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
ATC系统中短期冲突探测方法的研究;程先峰;龚维强;;网络安全技术与应用;20100915(第09期);全文 *
一种改进的基于ADS-B信息的空域冲突检测算法;王啸;赵修斌;;科学技术与工程(第32期);全文 *
低空通航飞行计划预先调配方法研究;田宇;张洪海;钱晓鹏;许卫卫;王立超;;武汉理工大学学报(交通科学与工程版)(第03期);全文 *
基于4D航迹的飞机连续下降运行飞行引导技术;李广文;贾秋玲;齐林;张翔伦;章卫国;;中国科学:技术科学;20180319(第03期);全文 *
基于全球剖分网格的空间数据编码方法;金安;程承旗;;测绘科学技术学报(第03期);全文 *
基于多级信息网格的海量遥感数据存储管理研究;李爽;程承旗;童晓冲;陈波;翟卫欣;;测绘学报(第S1期);全文 *
多级地理空间网格框架及其关键技术初探;关丽;吕雪锋;;地理与地理信息科学(第03期);全文 *
混合空域中无人机飞行防相撞技术;沈华;陈金良;周志靖;张建峰;;指挥信息系统与技术(第06期);全文 *
虚拟地理环境时空建模及其作战计划推演应用;李锋;万刚;蒋秉川;安籽鹏;;测绘学报(第08期);全文 *

Also Published As

Publication number Publication date
CN112348960A (zh) 2021-02-09

Similar Documents

Publication Publication Date Title
CN112348960B (zh) 一种适用于全球空间范围的空域冲突检测方法
Cheng et al. Using big data to track marine oil transportation along the 21st-century Maritime Silk Road
CN111291776A (zh) 基于众源轨迹数据的航道信息提取方法
CN110633262B (zh) 基于Spark的地图相交区域面积计算方法及系统
EP4206611A1 (en) Data processing method and device
CN116069882B (zh) 空域网格图生成方法
CN113139760B (zh) 基于风雨大数据的台风风险综合评估方法及系统
CN104422451A (zh) 一种道路识别方法及装置
CN112862176A (zh) 一种公共服务设施选址方法及装置
CN110866015A (zh) 一种基于局部网格的动目标活动范围记录方法
CN112069371A (zh) 一种基于格网化邻域的海上船舶有效救援搜索方法
CN112486987A (zh) 基于经纬度的城市查询方法、装置、设备及存储介质
CN114238384B (zh) 区域定位方法、装置、设备和存储介质
CN110489511B (zh) 等高线接边高程错误修正方法、系统及电子设备和介质
Wang et al. Spatial-temporal characteristics and causes of changes to the county-level administrative toponyms cultural landscape in the eastern plains of China
Li et al. Road extraction from high spatial resolution remote sensing image based on multi-task key point constraints
Runtuh et al. Representing landslides as polygon (areal) or points? How different data types influence the accuracy of landslide susceptibility maps
CN114219911A (zh) 一种基于立体剖分框架的空域网格化建模方法
CN114116948A (zh) 地理矢量数据空间缓冲区分析方法、装置、设备及介质
He et al. The fractal or scaling perspective on progressively generated intra-urban clusters from street junctions
CN117312471B (zh) 海量位置点的海陆属性判断方法及装置
CN114637305B (zh) 一种无人机最短路径规划方法及装置
CN114139099B (zh) 基于轨迹密度匀化和层次分割的道路交叉口信息提取方法
Rahımı et al. Monitoring the Long Term Urban Expansion of Zanjan City Using Remote Sensing and Geographic Information Systems
CN116303846A (zh) 一种基于全球网格的区域快速生成方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant