CN116541799B - 航道结构数据精度提升方法、***、电子设备及存储介质 - Google Patents
航道结构数据精度提升方法、***、电子设备及存储介质 Download PDFInfo
- Publication number
- CN116541799B CN116541799B CN202310506127.XA CN202310506127A CN116541799B CN 116541799 B CN116541799 B CN 116541799B CN 202310506127 A CN202310506127 A CN 202310506127A CN 116541799 B CN116541799 B CN 116541799B
- Authority
- CN
- China
- Prior art keywords
- data
- modulation degree
- degree difference
- sampling point
- ddm
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 70
- 238000000547 structure data Methods 0.000 title claims abstract description 45
- 238000005070 sampling Methods 0.000 claims abstract description 143
- 238000012545 processing Methods 0.000 claims abstract description 18
- 238000012937 correction Methods 0.000 claims abstract description 13
- 238000007499 fusion processing Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000010354 integration Effects 0.000 claims description 12
- 238000013480 data collection Methods 0.000 claims description 11
- 230000001133 acceleration Effects 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 6
- 230000004927 fusion Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 4
- 238000011835 investigation Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010006 flight Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/251—Fusion techniques of input or preprocessed data
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Navigation (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种航道结构数据精度提升方法、***、电子设备及存储介质,其方法包括:S1、采集QAR数据并提取航空器飞行相关的航迹数据及飞行性能数据;S2、将航班着陆跑道入口、航向台位置的经纬度坐标数据及航迹经纬度坐标数据经过横轴墨卡托投影转换成平面直角坐标系;S3、对航迹经纬度坐标数据的相邻两个采样点坐标进行插值处理;S4、在平面直角坐标系上计算采样点至跑道延长线的垂直距离di,以及采样点与跑道延长线垂直交点至航向台的距离Di;S5、将调制度差DDMi与采集模块提取到的调制度差DDM进行数据对应比较并修正融合处理。本发明经过对航迹数据和航道结构数据进行处理与修正后,极大地提高了航道结构数据的精度和准确性。
Description
技术领域
本发明涉及航道结构数据处理领域,尤其涉及一种航道结构数据精度提升方法、***、电子设备及存储介质。
背景技术
仪表着陆***是机场最后进近着陆阶段最重要的导航设备,其运行状态直接影响飞机着陆和机场运行标准,对机场的运行和安全有重要作用;其由航向台、下滑台和测距仪台组成。近年来,大型繁忙机场航班运行中飞行员反映航向台信号不稳定的事件时有发生,且由于航向台信号不稳定导致中止进近或复飞等不安全事件,为机场运行人员带来巨大压力。一旦机场发生中止进近或复飞等不安全事件后,机场运行人员需进行不安全事件的追溯和原因分析。
目前不安全追溯事件和原因排查主要采用的方式有设备故障的检查、周边干扰源的排查或通过实际的飞行校验。事件原因排查时间较长,且会花费大量的人力和物力,同时也降低了原因排查的效率,为机场安全运行带来了巨大的风险和威胁。由于航班QAR数据(QAR,英文全称Quick Access Recorder,是一种快速存取记录器,通过机载设备采集记录飞机运行中的各类数据,记录容量为128MB,可连续记录时间达600h,可以同时采集数百个数据,其包含了各种参量、飞机位置信息和机组操作等与飞机运行密切相关的数据信息,经由地面译码软件解析后,可对飞机飞行的状态和情况进行分析)具有数据内容丰富、数据具有高维时序特性且采样间隔规律等优点,可通过对航班QAR数据解析获取航空器在飞行过程中机载设备接收到的航向台信号相关数据信息。由于飞行干扰与扰动、机载设备精度与数据采集模式,OAR记录设备本身性能和精度等多方面的原因,航班OAR数据精度就相对不高,以至于无法得到高精度的航道结构数据以影响航道结构数据结果的分析与判断。基于此,如何实现航道结构数据(核心为调制度差数据)的精度提升,便于得到精度较高且修正后的调制度差数据,可有效提高了航空器飞行过程中不安全事件的追溯和原因查找的效率和准确性。
发明内容
本发明的目的在于克服背景技术所指出的技术问题,提供一种航道结构数据精度提升方法、***、电子设备及存储介质,通过提取航班在飞行进近阶段计算接收端记录的航向台相关的数据,经过对航迹数据和航道结构数据进行处理与修正后,极大地提高了航道结构数据的精度和准确性,为机场相关人员在事件排查过程中提供了强有力的精准数据支撑。
本发明的目的通过下述技术方案实现:
一种航道结构数据精度提升方法,其方法包括:
S1、采集QAR数据进行解析、译码、滤波处理后,提取航空器飞行相关的航迹数据及飞行性能数据,航迹数据包括航迹经纬度坐标数据,飞行性能数据包括飞机地速、空速、水平加速度、真航向:
S2、采集航班着陆跑道入口的经纬度坐标数据和航向台位置的经纬度坐标数据,将航班着陆跑道入口的经纬度坐标数据、航向台位置的经纬度坐标数据、航迹经纬度坐标数据经过横轴墨卡托投影转换成平面直角坐标系,在平面直角坐标系上绘制航道扇区的边缘线,航班着陆跑道入口的经纬度坐标数据记为(x0,y0),航向台位置的经纬度坐标数据记为(xL,yL),航迹经纬度坐标数据记为(xj,yj),i=1,2,3,……,N;
S3、对航迹经纬度坐标数据的相邻两个采样点坐标进行插值处理:
根据飞行性能数据通过二重积分的方法插值得到新的采样点方法如下:
其中(x′j,y′j)表示以第j个采样点的坐标经过二重积分方法插值得到新的采样点坐标,(xa,ya)表示已知A点的坐标,已知A点为一个已知采样点或航班着陆跑道入口处或航向台位置;va表示A点的地速,aj表示第j个采样点的水平加速度,θj表示第j个采样点的真航向,tj表示第j个采样点的时刻,tj+1表示第j+1个采样点的时刻;
在所有相邻两个采样点坐标之间依次***N0个采样点并组合成新的航迹经纬度坐标数据,同时将新的采样点经过横轴墨卡托投影转换成平面直角坐标系;
S4、在平面直角坐标系上计算各个采样点至跑道延长线的垂直距离di,以及各个采样点与跑道延长线垂直交点至航向台的距离Di,计算公式如下:
其中k0表示跑道延长线的斜率,c0表示在跑道延长线上航班着陆跑道入口到航向台位置直线段的截距;(Xi,Yi)表示新的航迹经纬度坐标数据中的采样点i的坐标;若采样点i位于跑道延长线左侧,则垂直距离di为正值,并记为d′i,d′i=di;若采样点i位于跑道延长线右侧,则垂直距离di为负值,并记为d′i,d′i=-di;
其中ki表示采样点i与跑道延长线垂直交点所在直线的斜率,ci表示采样点i与跑道延长线垂直交点之间直线段的截距;
将采样点i的垂直距离di与距离Di对应关联存储,按照如下公式计算得到采样点i的调制度差DDMi:
其中Wi/2表示采样点i与跑道延长线垂直交点至航道扇区边缘线的距离;
由此计算得到所有采样点随距离Di的调制度差DDMi,由此得到航道结构数据的调制度差DDMi数据集合;
S5、从OAR数据中提取调制度差DDM,设定阈值F,按照如下方法进行数据比较并修正融合处理;
S51、若调制度差DDMi数据集合中的调制度差DDMi与提取到的调制度差DDM存在数据正负差异,则保留调制度差DDMi,并删除提取到的调制度差DDM;
S52、若调制度差DDMi数据集合中的调制度差DDMi与提取到的调制度差DDM差值小于或等于阈值F,则采用提取到的调制度差DDM替换所对应的调制度差DDMi;
S53、若调制度差DDMi数据集合中的调制度差DDMi与提取到的调制度差DDM差值大于阈值F,则采保留调制度差DDMi,并删除提取到的调制度差DDM;
由此得到修正后的调制度差DDMi数据集合。
为了更好地实现本发明航道结构数据精度提升方法,QAR数据中航迹数据的采样时间间隔为T0;在步骤S3中,若相邻两个采样点坐标的间隔时间小于或等于T0,则在相邻两个采样点坐标之间依据飞行性能数据通过二重积分方法依次***n0个采样点;若相邻两个采样点坐标的间隔时间大于T0,则在相邻两个采样点坐标之间依据飞行性能数据通过二重积分方法依次***至少2n0个采样点。
本发明航道结构数据精度提升方法优选的技术方案是:在步骤S4中,跑道延长线的斜率k0通过如下公式计算得到:
在跑道延长线上航班着陆跑道入口到航向台位置直线段的截距c0通过如下公式计算得到:
本发明航道结构数据精度提升方法优选的技术方案是:采样点i与跑道延长线垂直交点所在直线的斜率ki通过如下公式计算得到:
采样点i与跑道延长线垂直交点之间直线段的截距ci通过如下公式计算得到:
一种航道结构数据精度提升***,包括:
采集模块,用于采集OAR数据进行解析、译码、滤波处理后,提取航空器飞行相关的航迹数据及飞行性能数据,航迹数据包括航迹经纬度坐标数据,飞行性能数据包括飞机地速、空速、水平加速度、真航向:采集模块还用于从OAR数据中提取调制度差DDM以及采集航班着陆跑道入口的经纬度坐标数据和航向台位置的经纬度坐标数据;
计算处理***,包括插值计算模块、调制度差计算模块和调制度差修正融合模块,计算处理***内部构建有成平面直角坐标系并将航班着陆跑道入口的经纬度坐标数据、航向台位置的经纬度坐标数据、航迹经纬度坐标数据经过横轴墨卡托投影转换成平面直角坐标系,在平面直角坐标系上绘制航道扇区的边缘线;插值计算模块对航迹经纬度坐标数据的相邻两个采样点坐标进行插值处理,在所有相邻两个采样点坐标之间依次***N0个采样点并组合成新的航迹经纬度坐标数据,同时将新的采样点经过横轴墨卡托投影转换成平面直角坐标系;调制度差计算模块计算得到所有采样点随距离Di的调制度差DDMi并得到航道结构数据的调制度差DDMi数据集合;调制度差修正融合模块设定阈值F,将调制度差DDMi与采集模块提取到的调制度差DDM进行数据对应比较并修正融合处理,由此得到修正后的调制度差DDMi数据集合;
输出模块,用于向外输出修正后的调制度差DDMi数据。
一种电子设备,包括:至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行本发明航道结构数据精度提升方法的步骤。
一种存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现本发明航道结构数据精度提升方法的步骤。
本发明较现有技术相比,具有以下优点及有益效果:
(1)本发明通过提取航班在飞行进近阶段计算接收端记录的航向台相关的数据,经过对航迹数据和航道结构数据进行处理与修正后,极大地提高了航道结构数据的精度和准确性,能够得到高精度的调制度差DDMi数据,为机场相关人员在事件排查过程中提供了强有力的精准数据支撑,便于实现仪表着陆***的信号评估、事件有效排查。
(2)本发明将航迹经纬度数据转换成平面直角坐标系,创造性地通过二重积分对航迹数据的采样点进行插值处理,求得采样点(或称航迹点)距跑道中线的垂直距离以及航向台至两条直线交点的距离等距离数据,方法简单,结果准确,可快速实现对航道结构数据进行修正处理,提高了航道结构数据的精度。
(3)本发明可直接通过航班OAR数据来分析航向台的航道结构数据,节约了人工信号测量的时间成本,提高了不安全事件排查的效率;经过对QAR数据中航向台航道结构数据相关的数据进行修正,提高了航道结构数据的精度,准确性高;为航向台信号不稳定性原因排查与分析提供了精准数据支撑,也为机场航空器的安全运行提供了技术保障,具有数据生产效率高、质量高、数据可靠等优点。
附图说明
图1为本发明航道结构数据精度提升方法的方法流程示意图;
图2为本发明步骤S5中数据比较并修正融合处理的方法流程示意图;
图3为实施例中平面直角坐标系以采样点i为例的关系示意图;
图4为实施例中航道结构数据修正前后的对比示意图;
图5为本发明航道结构数据精度提升***的功能模块框图。
具体实施方式
下面结合实施例对本发明作进一步地详细说明:
实施例
如图1~图4所示,一种航道结构数据精度提升方法,其方法包括:
S1、采集OAR数据(采集的OAR数据包括航班号、飞行日期、经纬度、飞行高度、水平加速度、地速、空速、真航向、磁航向、航向台DDM值、航向台工作频率和航班着陆跑道号等相关数据;由于航班OAR数据是原码数据,需经对OAR数据进行解析并译码后转为可读取数据;OAR数据中航迹数据的采样时间间隔为T0;本实施例所使用的QAR数据采样时间间隔为1秒)进行解析、译码、滤波处理后,提取航空器飞行相关的航迹数据及飞行性能数据,航迹数据包括航迹经纬度坐标数据,飞行性能数据包括飞机地速、空速、水平加速度、真航向:
S2、采集航班着陆跑道入口的经纬度坐标数据和航向台位置的经纬度坐标数据,将航班着陆跑道入口的经纬度坐标数据、航向台位置的经纬度坐标数据、航迹经纬度坐标数据经过横轴墨卡托投影转换成平面直角坐标系,在平面直角坐标系上绘制航道扇区的边缘线(从航道线向两边扩展到DDM为0.155(150μA)的各点轨迹所限制的区域;在跑道入口两边105m处的DDM为0.155(150μA),例如as为航道扇区角,根据航向台至跑道末端的距离、跑道长度和跑道入口两边的宽度可计算出航道扇区角的值),航班着陆跑道入口的经纬度坐标数据记为(x0,y0),参见图3的示例;航向台位置的经纬度坐标数据记为(xL,yL),参见图3的示例;航迹经纬度坐标数据记为(xj,yj),i=1,2,3,……,N。
S3、对航迹经纬度坐标数据的相邻两个采样点坐标进行插值处理。
根据飞行性能数据通过二重积分的方法插值得到新的采样点方法如下:
其中(x′j,y′j)表示以第j个采样点的坐标经过二重积分方法插值得到新的采样点坐标,(xa,ya)表示已知A点的坐标,已知A点为一个已知采样点或航班着陆跑道入口处或航向台位置;va表示A点的地速,aj表示第j个采样点的水平加速度,θj表示第j个采样点的真航向,tj表示第j个采样点的时刻,tj+1表示第j+1个采样点的时刻。
在所有相邻两个采样点坐标之间依次***N0个采样点并组合成新的航迹经纬度坐标数据,同时将新的采样点经过横轴墨卡托投影转换成平面直角坐标系。正常情况下,QAR数据中采集的航迹数据会按照采样时间间隔T0进行航迹采样点的采集,即航迹数据中任两个采样点坐标的间隔时间都等于T0。由于原始数据在译码解析过程会存在部分数字段的错位、信息缺失、参数取值超出理论取值范围和出现不合逻辑的跳变等异常情况,需结合异常数据所处时间点附近一段时间内的飞机状态相关的参数数据,对异常数据进行识别与删除。
若相邻两个采样点坐标的间隔时间小于或等于T0,则在相邻两个采样点坐标之间依据飞行性能数据通过二重积分方法依次***n0个采样点(即N0=n0,相邻两个采样点可以为航迹数据中两个相邻的采样点,也可以为航迹数据中一个采样点与相邻***新的采样点,也可以为两个相邻***新的采样点)。若相邻两个采样点坐标的间隔时间大于T0,则在相邻两个采样点坐标之间依据飞行性能数据通过二重积分方法依次***至少2n0个采样点(即N0=2n0)。上述所提到的采样点或***采样点也即航迹数据中的航迹点。
S4、如图3所示,在平面直角坐标系上计算各个采样点至跑道延长线的垂直距离di(参见图3以采样点i为例的距离示例),以及各个采样点与跑道延长线垂直交点至航向台的距离Di(参见图3以采样点i为例的距离示例),计算公式如下:
其中k0表示跑道延长线的斜率,c0表示在跑道延长线上航班着陆跑道入口到航向台位置直线段的截距;(Xi,Yi)表示新的航迹经纬度坐标数据中的采样点i的坐标;若采样点i位于跑道延长线左侧,则垂直距离di为正值,并记为d′i,d′i=di;若采样点i位于跑道延长线右侧,则垂直距离di为负值,并记为d′i,d′i=-di。
跑道延长线的斜率k0通过如下公式计算得到:
在跑道延长线上航班着陆跑道入口到航向台位置直线段的截距c0通过如下公式计算得到:
其中ki表示采样点i与跑道延长线垂直交点所在直线的斜率,ci表示采样点i与跑道延长线垂直交点之间直线段的截距。跑道延长线在水平面内最靠近跑道中心的DDM为零的点集合,也即跑道中心线的延长线。
采样点i与跑道延长线垂直交点所在直线的斜率ki通过如下公式计算得到:
采样点i与跑道延长线垂直交点之间直线段的截距ci通过如下公式计算得到:
将采样点i的垂直距离di与距离Di对应关联存储,按照如下公式计算得到采样点i的调制度差DDMi:
其中Wi/2表示采样点i与跑道延长线垂直交点至航道扇区边缘线的距离;
由此计算得到所有采样点随距离Di的调制度差DDMi,由此得到航道结构数据的调制度差DDMi数据集合;
S5、从OAR数据中提取调制度差DDM(即OAR数据所提取出的调制度差DDM,具体为150Hz和90Hz两个音频信号对射频的调制度百分数的差值),设定阈值F,按照如下方法(参见图2)进行数据比较并修正融合处理;
S51、若调制度差DDMi数据集合中的调制度差DDMi与提取到的调制度差DDM存在数据正负差异,则保留调制度差DDMi,并删除提取到的调制度差DDM;
S52、若调制度差DDMi数据集合中的调制度差DDMi与提取到的调制度差DDM差值小于或等于阈值F,则采用提取到的调制度差DDM替换所对应的调制度差DDMi;
S53、若调制度差DDMi数据集合中的调制度差DDMi与提取到的调制度差DDM差值大于阈值F,则采保留调制度差DDMi,并删除提取到的调制度差DDM;
由此得到修正后的调制度差DDMi数据集合。以某机场I类运行的航向台航道结构数据和该机场的着陆航班OAR数据为例,修正后的调制度差DDMi数据和OAR数据中直接读取的DDM值进行对比,对比曲线如图4所示,从图中可与看出,修正后的调制度差DDMi数据(也即本申请航道结构数据)精度明显得到了较高提升。
如图5所示,一种航道结构数据精度提升***,包括:
采集模块,用于采集OAR数据进行解析、译码、滤波处理后,提取航空器飞行相关的航迹数据及飞行性能数据,航迹数据包括航迹经纬度坐标数据,飞行性能数据包括飞机地速、空速、水平加速度、真航向:采集模块还用于从OAR数据中提取调制度差DDM以及采集航班着陆跑道入口的经纬度坐标数据和航向台位置的经纬度坐标数据;
计算处理***,包括插值计算模块、调制度差计算模块和调制度差修正融合模块,计算处理***内部构建有成平面直角坐标系并将航班着陆跑道入口的经纬度坐标数据、航向台位置的经纬度坐标数据、航迹经纬度坐标数据经过横轴墨卡托投影转换成平面直角坐标系,在平面直角坐标系上绘制航道扇区的边缘线;插值计算模块对航迹经纬度坐标数据的相邻两个采样点坐标进行插值处理,在所有相邻两个采样点坐标之间依次***N0个采样点并组合成新的航迹经纬度坐标数据,同时将新的采样点经过横轴墨卡托投影转换成平面直角坐标系;调制度差计算模块计算得到所有采样点随距离Di的调制度差DDMi并得到航道结构数据的调制度差DDMi数据集合;调制度差修正融合模块设定阈值F,将调制度差DDMi与采集模块提取到的调制度差DDM进行数据对应比较并修正融合处理,由此得到修正后的调制度差DDMi数据集合;
输出模块,用于向外输出修正后的调制度差DDMi数据。
一种电子设备,包括:至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行本发明航道结构数据精度提升方法的步骤。
一种存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现本发明航道结构数据精度提升方法的步骤。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种航道结构数据精度提升方法,其特征在于:其方法包括:
S1、采集QAR数据进行解析、译码、滤波处理后,提取航空器飞行相关的航迹数据及飞行性能数据,航迹数据包括航迹经纬度坐标数据,飞行性能数据包括飞机地速、空速、水平加速度、真航向:
S2、采集航班着陆跑道入口的经纬度坐标数据和航向台位置的经纬度坐标数据,将航班着陆跑道入口的经纬度坐标数据、航向台位置的经纬度坐标数据、航迹经纬度坐标数据经过横轴墨卡托投影转换成平面直角坐标系,在平面直角坐标系上绘制航道扇区的边缘线,航班着陆跑道入口的经纬度坐标数据记为(x0,y0),航向台位置的经纬度坐标数据记为(xL,yL),航迹经纬度坐标数据记为(xj,yj),i=1,2,3,……,N;
S3、对航迹经纬度坐标数据的相邻两个采样点坐标进行插值处理:
根据飞行性能数据通过二重积分的方法插值得到新的采样点方法如下:
其中(x′j,y′j)表示以第j个采样点的坐标经过二重积分方法插值得到新的采样点坐标,(xa,ya)表示已知A点的坐标,已知A点为一个已知采样点或航班着陆跑道入口处或航向台位置;va表示A点的地速,aj表示第j个采样点的水平加速度,θj表示第j个采样点的真航向,tj表示第j个采样点的时刻,tj+1表示第j+1个采样点的时刻;
在所有相邻两个采样点坐标之间依次***N0个采样点并组合成新的航迹经纬度坐标数据,同时将新的采样点经过横轴墨卡托投影转换成平面直角坐标系;
S4、在平面直角坐标系上计算各个采样点至跑道延长线的垂直距离di,以及各个采样点与跑道延长线垂直交点至航向台的距离Di,计算公式如下:
其中k0表示跑道延长线的斜率,c0表示在跑道延长线上航班着陆跑道入口到航向台位置直线段的截距;(Xi,Yi)表示新的航迹经纬度坐标数据中的采样点i的坐标;若采样点i位于跑道延长线左侧,则垂直距离di为正值,并记为d′i,d′i=di;若采样点i位于跑道延长线右侧,则垂直距离di为负值,并记为d′i,d′i=-di;
其中ki表示采样点i与跑道延长线垂直交点所在直线的斜率,ci表示采样点i与跑道延长线垂直交点之间直线段的截距;
将采样点i的垂直距离di与距离Di对应关联存储,按照如下公式计算得到采样点i的调制度差DDMi:
其中Wi/2表示采样点i与跑道延长线垂直交点至航道扇区边缘线的距离;
由此计算得到所有采样点随距离Di的调制度差DDMi,由此得到航道结构数据的调制度差DDMi数据集合;
S5、从QAR数据中提取调制度差DDM,设定阈值F,按照如下方法进行数据比较并修正融合处理;
S51、若调制度差DDM数据集合中的调制度差DDMi与提取到的调制度差DDM存在数据正负差异,则保留调制度差DDMi,并删除提取到的调制度差DDM;
S52、若调制度差DDMi数据集合中的调制度差DDMi与提取到的调制度差DDM差值小于或等于阈值F,则采用提取到的调制度差DDM替换所对应的调制度差DDMi;
S53、若调制度差DDMi数据集合中的调制度差DDMi与提取到的调制度差DDM差值大于阈值F,则采保留调制度差DDMi,并删除提取到的调制度差DDM;
由此得到修正后的调制度差DDMi数据集合。
2.按照权利要求1所述的航道结构数据精度提升方法,其特征在于:QAR数据中航迹数据的采样时间间隔为T0;在步骤S3中,若相邻两个采样点坐标的间隔时间小于或等于T0,则在相邻两个采样点坐标之间依据飞行性能数据通过二重积分方法依次***n0个采样点;若相邻两个采样点坐标的间隔时间大于T0,则在相邻两个采样点坐标之间依据飞行性能数据通过二重积分方法依次***至少2n0个采样点。
3.按照权利要求1所述的航道结构数据精度提升方法,其特征在于:在步骤S4中,跑道延长线的斜率k0通过如下公式计算得到:
在跑道延长线上航班着陆跑道入口到航向台位置直线段的截距c0通过如下公式计算得到:
4.按照权利要求1所述的航道结构数据精度提升方法,其特征在于:采样点i与跑道延长线垂直交点所在直线的斜率ki通过如下公式计算得到:
采样点i与跑道延长线垂直交点之间直线段的截距ci通过如下公式计算得到:
5.一种航道结构数据精度提升***,其特征在于:包括:
采集模块,用于采集QAR数据进行解析、译码、滤波处理后,提取航空器飞行相关的航迹数据及飞行性能数据,航迹数据包括航迹经纬度坐标数据,飞行性能数据包括飞机地速、空速、水平加速度、真航向:采集模块还用于从QAR数据中提取调制度差DDM以及采集航班着陆跑道入口的经纬度坐标数据和航向台位置的经纬度坐标数据;
计算处理***,包括插值计算模块、调制度差计算模块和调制度差修正融合模块,计算处理***内部构建有成平面直角坐标系并将航班着陆跑道入口的经纬度坐标数据、航向台位置的经纬度坐标数据、航迹经纬度坐标数据经过横轴墨卡托投影转换成平面直角坐标系,在平面直角坐标系上绘制航道扇区的边缘线;插值计算模块对航迹经纬度坐标数据的相邻两个采样点坐标进行插值处理,在所有相邻两个采样点坐标之间依次***N0个采样点并组合成新的航迹经纬度坐标数据,同时将新的采样点经过横轴墨卡托投影转换成平面直角坐标系;调制度差计算模块计算得到所有采样点随距离Di的调制度差DDMi并得到航道结构数据的调制度差DDMi数据集合,采样点随距离Di的调制度差DDMi获得方法如下:
在平面直角坐标系上计算各个采样点至跑道延长线的垂直距离di,以及各个采样点与跑道延长线垂直交点至航向台的距离Di,计算公式如下:
其中k0表示跑道延长线的斜率,c0表示在跑道延长线上航班着陆跑道入口到航向台位置直线段的截距;(Xi,Yi)表示新的航迹经纬度坐标数据中的采样点i的坐标;若采样点i位于跑道延长线左侧,则垂直距离di为正值,并记为d′i,d′i=di;若采样点i位于跑道延长线右侧,则垂直距离di为负值,并记为d′i,d′i=-di;
其中ki表示采样点i与跑道延长线垂直交点所在直线的斜率,ci表示采样点i与跑道延长线垂直交点之间直线段的截距;
将采样点i的垂直距离di与距离Di对应关联存储,按照如下公式计算得到采样点i的调制度差DDMi:
其中Wi/2表示采样点i与跑道延长线垂直交点至航道扇区边缘线的距离;调制度差修正融合模块设定阈值F,将调制度差DDMi与采集模块提取到的调制度差DDM进行数据对应比较并修正融合处理,由此得到修正后的调制度差DDMi数据集合;
输出模块,用于向外输出修正后的调制度差DDMi数据。
6.一种电子设备,其特征在于:包括:至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器执行如权利要求1-4任一所述的方法的步骤。
7.一种存储介质,其上存储有计算机程序,其特征在于:所述计算机程序被处理器执行时实现如权利要求1-4中任一项所述的方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310506127.XA CN116541799B (zh) | 2023-05-06 | 2023-05-06 | 航道结构数据精度提升方法、***、电子设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310506127.XA CN116541799B (zh) | 2023-05-06 | 2023-05-06 | 航道结构数据精度提升方法、***、电子设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116541799A CN116541799A (zh) | 2023-08-04 |
CN116541799B true CN116541799B (zh) | 2023-11-07 |
Family
ID=87453756
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310506127.XA Active CN116541799B (zh) | 2023-05-06 | 2023-05-06 | 航道结构数据精度提升方法、***、电子设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116541799B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117407477B (zh) * | 2023-10-26 | 2024-05-14 | 航科院中宇(北京)新技术发展有限公司 | 地理信息数据演变识别处理方法、***及存储介质 |
CN118049995B (zh) * | 2024-04-16 | 2024-07-05 | 暨南大学 | 一种基于卫惯组合技术的船载智能导航方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110456655A (zh) * | 2019-06-29 | 2019-11-15 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 微波着陆仿真模拟*** |
CN114754775A (zh) * | 2022-04-07 | 2022-07-15 | 广州嘉航通信科技有限公司 | 航道变化分析预警方法、***、计算机设备、及存储介质 |
WO2022257620A1 (zh) * | 2021-06-08 | 2022-12-15 | 东航技术应用研发中心有限公司 | 着陆空中距离的测量方法和*** |
CN115798265A (zh) * | 2022-11-01 | 2023-03-14 | 航科院中宇(北京)新技术发展有限公司 | 一种基于数字孪生技术的数字塔台构建方法及其实现*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE112020007038T5 (de) * | 2020-06-11 | 2023-01-26 | Mitsubishi Electric Corporation | Radareinrichtung und radarbild-generierungsverfahren |
-
2023
- 2023-05-06 CN CN202310506127.XA patent/CN116541799B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110456655A (zh) * | 2019-06-29 | 2019-11-15 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 微波着陆仿真模拟*** |
WO2022257620A1 (zh) * | 2021-06-08 | 2022-12-15 | 东航技术应用研发中心有限公司 | 着陆空中距离的测量方法和*** |
CN114754775A (zh) * | 2022-04-07 | 2022-07-15 | 广州嘉航通信科技有限公司 | 航道变化分析预警方法、***、计算机设备、及存储介质 |
CN115798265A (zh) * | 2022-11-01 | 2023-03-14 | 航科院中宇(北京)新技术发展有限公司 | 一种基于数字孪生技术的数字塔台构建方法及其实现*** |
Also Published As
Publication number | Publication date |
---|---|
CN116541799A (zh) | 2023-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116541799B (zh) | 航道结构数据精度提升方法、***、电子设备及存储介质 | |
CN110515116B (zh) | 一种无人机航空伽马能谱探测和数据处理***及其处理方法 | |
CN107944701B (zh) | 一种飞机着陆过程中冲出跑道风险的检测方法及装置 | |
Iliff et al. | Practical aspects of a maximum likelihood estimation method to extract stability and control derivatives from flight data | |
CN115017578A (zh) | 一种建筑的智能实测实量方法、装置、ugv及存储介质 | |
CN108073865B (zh) | 一种基于卫星资料的飞机尾迹云识别方法 | |
CN114694130A (zh) | 基于深度学习的铁路沿线电线杆及杆号检测方法和装置 | |
CN113311858A (zh) | 一种基于无人机的飞控***及方法 | |
CN112542800A (zh) | 一种输电线路故障识别的方法及其*** | |
CN116774191A (zh) | 一种pos数据和机载激光雷达点云数据的交互质检装置 | |
CN116307822A (zh) | 一种通用型飞行训练评估*** | |
CN114782847B (zh) | 一种基于无人机的矿山产能监测方法及装置 | |
CN103837130B (zh) | 用于机载激光扫描***的数据处理方法及装置 | |
CN114200527B (zh) | 一种基于倾斜摄影的无人机航磁测量方法及*** | |
Behere et al. | A Method for the Parametric Representation of Take-off Time-Series Trajectory Data for Environmental Impact Assessment | |
CN107967714A (zh) | 一种通过无人机数字高程模型自动提取森林郁闭度的方法 | |
CN114035189A (zh) | 一种道路隐患定位方法、装置、电子设备及存储介质 | |
CN107588760B (zh) | 适用于固定翼无人机跟踪道路的定高高精度航迹生成方法 | |
CN111122900B (zh) | 一种磁差修正方法及装置 | |
CN115060457B (zh) | 基于飞行器颠簸探测大气涡旋耗散率的方法、***及设备 | |
CN117685975B (zh) | 用于航空事故调查的飞机着陆轨迹精确计算方法及*** | |
CN116453377B (zh) | 一种对飞机qar数据进行飞行阶段划分的方法 | |
CN117612047B (zh) | 基于ai大模型的电网用无人机巡检图像识别方法 | |
CN115359686B (zh) | 仪表着陆***信号分析诊断方法、装置及电子设备 | |
CN114881543B (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 |