CN106530838B - A kind of method of aircraft flight profiles jump removal - Google Patents
A kind of method of aircraft flight profiles jump removal Download PDFInfo
- Publication number
- CN106530838B CN106530838B CN201610931969.XA CN201610931969A CN106530838B CN 106530838 B CN106530838 B CN 106530838B CN 201610931969 A CN201610931969 A CN 201610931969A CN 106530838 B CN106530838 B CN 106530838B
- Authority
- CN
- China
- Prior art keywords
- longitude
- latitude
- data
- time
- value
- 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
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0004—Transmission of traffic-related information to or from an aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
- Traffic Control Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a kind of methods of aircraft flight profiles jump removal, the following steps are included: the cubic polynomial fitting algorithm of the on-fixed time interval of one corresponding types of building, finds out the corresponding predicted value of height above sea level, angle, ground velocity, vertical speed and stores；When receiving latitude and longitude value, if longitude and latitude number of elements before current time within the scope of 1 minute is less than 3 in set, without calculating and jumping out；Otherwise current time longitude and latitude was compared with the longitude and latitude at a upper time point, the longitude and latitude at a upper time point respectively: calculates Prediction distance and actual range, the greater among the difference of comparison prediction distance and actual range, Prediction distance and actual range the two obtains data reliability weight；Data reliability weight < 1 is then regarded as trust data and is stored.The present invention is able to solve the problem of positioning longitude and latitude or height jump, keeps real time positioning data more acurrate, and with the flight path of smoother data exhibiting aircraft.
Description
Technical field
The present invention relates to a kind of methods of aircraft flight profiles jump removal, belong to aircraft flight profiles optimization method
Technical field.
Background technique
Grow along with the continuous development of China Aviation cause, it is also constantly prominent the problems such as airspace, flow, flight delay,
Flow control leads to passenger, and complaints are heard everywhere.The appearance of this new technology of ADS-B can be obtained automatically from related airborne equipment and be joined
Number is to information such as position, height, speed, course, the identifiers of other aircrafts or earth station's broadcast aircraft, so that controller is to winged
Machine state is monitored.
As the application popularization rate of ADS-B on aircraft increases, the feasible of position of aircraft monitoring is carried out by ADS-B
Property it is also increasingly stronger；But since ADS-B transmitting range is limited, and clock is difficult to synchronize and causes data between data receiver base station
Randomness, ADS-B issue the problems such as data itself are existing wrong, the track for causing the ADS-B initial data received to construct
There are jump, track missings.Based on the above-mentioned state of the art, it is necessary to which design provides a kind of jump removal of aircraft flight profiles
Method keeps real time positioning data more acurrate to solve the above problems, and with the flight rail of smoother data exhibiting aircraft
Mark.
Summary of the invention
The present invention is exactly in view of the deficienciess of the prior art, providing a kind of side of aircraft flight profiles jump removal
Method is solved the problems, such as that positioning longitude and latitude or height jump, keeps real time positioning data more acurrate, and navigated with smoother data exhibiting
The flight path of pocket.
To solve the above problems, the technical solution used in the present invention is as follows:
A kind of method of aircraft flight profiles jump removal, comprising the following steps:
S1, longitude and latitude is being receivedHeight above sea levelAngleGround velocityVertical speedWhen, it is indicated with t
The timestamp for receiving the data, will be in the set A of metadata record to nearest three minute data of correspondence；
S2, it is receivingWhen value, the on-fixed time interval an of corresponding types is constructed
Cubic polynomial fitting algorithm operand O；O is used to carry out interpolation to each time interval of one-dimensional time series data, from
And achieve the purpose that data are predicted in the short time；
S3, it finds outCorresponding predicted value β:
And it is stored β as actual value；
S4, it is receivingWhen value, if AposIn number of elements before current time within the scope of 1 minute less than 3, then
Without calculating and jumping out；Otherwise, by current time longitude and latitudeThat is the longitude and latitude of current point in time in set ARespectively
With the longitude and latitude at a upper time pointThe longitude and latitude at a upper time pointBe compared, comparison process include with
Lower step:
S41, two functions B and K are defined, wherein B expression is asked between ta and corresponding two longitudes and latitudes of two timestamps of tb
Earth spherical distance, K expression seek the earth spherical surface to fly between ta and the corresponding two aircraft ground velocity of two timestamps of tb
Distance:
S42, the actual range that current location Yu the first two position location are found out according to the function B of definition:
S43, because of the timestamp of velocity vector and position vector and asynchronous, need to find out current location time point
(t3), the predetermined speed at the first two position time point (t1, t2):
S44, the Prediction distance between two time points is found out according to the function K defined:
S45, degree of error is set as a, the difference of Prediction distance and actual range is | ε-δ |, both Prediction distance and actual range
Among the greater be max { ε-δ }, will | ε-δ | be compared evaluation with max { ε-δ }, obtain data reliability weight w:
W=w1+w2；
Obtained w is the latitude and longitude value being currently receivedConfidence level, if < 1, willBe regarded as trust data into
Row storage.
As an improvement of the above technical solution, in the step S45, degree of error a is set as 3.
Compared with prior art, implementation result of the invention is as follows by the present invention:
A kind of method of aircraft flight profiles jump removal of the present invention, is sent out by the observation to historical data
It is existing: the longitude and latitude that ADS-B is issued, there are small probability, jumped on a large scale to a fixed position for a long time the case where, should
Data present regular, primary to a fixed position hopping at interval of two anchor points；Based on above-mentioned jump rule, and
Longitude and latitude, height above sea level, angle, ground velocity, the rule of status values in vertical speed are only issued every time according to ADS-B, and are counted
The method for obtaining aircraft flight profiles jump removal of the invention is calculated, asking for positioning longitude and latitude or height jump is able to solve
Topic, keeps real time positioning data more acurrate, and with the flight path of smoother data exhibiting aircraft.
Specific embodiment
Illustrate the contents of the present invention below in conjunction with specific embodiments.
By the observation to historical data, the longitude and latitude that ADS-B is issued, there are small probability, solid to one for a long time
The case where fixed position jumps on a large scale.
The data present it is regular, primary to a fixed position hopping at interval of two anchor points, and due to
ADS-B only issues a status values every time, such as a value in longitude and latitude, height above sea level, angle, ground velocity, vertical speed.
Based on above-mentioned rule and rule, a kind of method of aircraft flight profiles jump removal provided by the present embodiment,
Comprising the following specific steps
1, longitude and latitude is being receivedHeight above sea levelAngleGround velocityVertical speed
When, t indicates to receive the timestamp of the data, will be in the set A of metadata record to nearest three minute data of correspondence.
2, it is receivingWhen value, the on-fixed time interval an of corresponding types is constructed
Cubic polynomial fitting algorithm operand O；
O is used to carry out interpolation to each time interval of one-dimensional time series data, to reach the pre- of data in the short time
The purpose of survey.
3, it finds outCorresponding predicted value:
And it is stored β as actual value.
4, it is receivingWhen value, if AposIn 1 minute range interior element quantity is less than 3 before current time, then not
It is calculated and is jumped out；
Otherwise, willRespectively withIt is compared, process is as follows:
(1) two functions are defined, B indicates to seek the earth spherical surface between ta and corresponding two longitudes and latitudes of two timestamps of tb
The earth spherical distance to fly between ta and the corresponding two aircraft ground velocity of two timestamps of tb is sought in distance, K expression:
(2) actual range of current location Yu the first two position location is found out:
(3) because of the timestamp of velocity vector and position vector and asynchronous, need to find out current location time point (t3),
The predetermined speed at the first two position time point (t1, t2):
(4) Prediction distance between two time points is found out:
(5) degree of error is set as a, by Prediction distance compared with actual range, by the difference and Prediction distance and reality of two distances
The greater of border among the two is compared evaluation, obtains data reliability weight w:
W=w1+w2,
Since data longitude and latitude jump range is larger, generally hundred kilometers or more, and airline carriers of passengers flying speed 1 minute one
As be no more than 60 kilometers, so the larger value can be set by a, avoid the inexactness bring predicted due to velocity fitting will just
The problem of true longitude and latitude data exclude, a is traditionally arranged to be 3；
Obtained w is the latitude and longitude value being currently receivedConfidence level, if < 1, willBe regarded as trust data into
Row storage.
The foregoing is a detailed description of the present invention in conjunction with specific embodiments, and it cannot be said that the present invention is specifically real
It applies and is only limitted to these explanations.For those skilled in the art to which the present invention belongs, before not departing from present inventive concept
It puts, a number of simple deductions or replacements can also be made, all shall be regarded as belonging to the scope of protection of the invention.
Claims (2)
1. a kind of method of aircraft flight profiles jump removal, characterized in that the following steps are included:
S1, longitude and latitude is being receivedHeight above sea levelAngleGround velocityVertical speedWhen, it indicates to receive with t and be somebody's turn to do
The timestamp of data, will be in the set A of metadata record to nearest three minute data of correspondence；
S2, it is receivingWhen value, the three of the on-fixed time interval an of corresponding types are constructed
Order polynomial fitting algorithm operand O；O is used to carry out interpolation to each time interval of one-dimensional time series data, to reach
The purpose that data are predicted in the short time；
S3, it finds outCorresponding predicted value β:
And it is stored β as actual value；
S4, it is receivingWhen value, if AposIn number of elements before current time within the scope of 1 minute less than 3, then not into
Row is calculated and is jumped out；Otherwise, by current time longitude and latitudeThat is the longitude and latitude of current point in time in set ARespectively with it is upper
The longitude and latitude at one time pointThe longitude and latitude at a upper time pointIt is compared, comparison process includes following step
It is rapid:
S41, two functions B and K are defined, wherein B indicates to seek the ground between ta and corresponding two longitudes and latitudes of two timestamps of tb
The earth spherical distance to fly between ta and the corresponding two aircraft ground velocity of two timestamps of tb is sought in ball spherical distance, K expression:
S42, the actual range that current location Yu the first two position location are found out according to the function B of definition:
S43, because of the timestamp of velocity vector and position vector and asynchronous, need to find out current location time point (t3), preceding
Predetermined speed of two position time points (t1, t2):
S44, the Prediction distance between two time points is found out according to the function K defined:
S45, degree of error is set as a, the difference of Prediction distance and actual range is | ε-δ |, among Prediction distance and actual range the two
The greater be max { ε-δ }, will | ε-δ | be compared evaluation with max { ε-δ }, obtain data reliability weight w:
W=w1+w2；
Obtained w is the latitude and longitude value being currently receivedConfidence level, if < 1, willTrust data is regarded as to be deposited
Storage.
2. a kind of method of aircraft flight profiles jump removal as described in claim 1, characterized in that in the step
In S45, degree of error a is set as 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610931969.XA CN106530838B (en) | 2016-10-25 | 2016-10-25 | A kind of method of aircraft flight profiles jump removal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610931969.XA CN106530838B (en) | 2016-10-25 | 2016-10-25 | A kind of method of aircraft flight profiles jump removal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106530838A CN106530838A (en) | 2017-03-22 |
| CN106530838B true CN106530838B (en) | 2019-04-09 |
Family
ID=58292407
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610931969.XA Active CN106530838B (en) | 2016-10-25 | 2016-10-25 | A kind of method of aircraft flight profiles jump removal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106530838B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110017838A (en) * | 2019-02-13 | 2019-07-16 | 浩亚信息科技有限公司 | Schedule flight trajectory predictions method, equipment, medium and system |
| CN111968413B (en) * | 2020-08-26 | 2021-09-24 | 成都民航空管科技发展有限公司 | Flight plan synchronization method for regional control center and terminal area ATC system |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1020778A2 (en) * | 1999-01-18 | 2000-07-19 | Honeywell Inc. | Flight management system (FMS) |
| WO2004006159A1 (en) * | 2002-07-08 | 2004-01-15 | Proteome Systems Intellectuel Property Pty Ltd | Method and system for picking peaks for mass spectra |
| CN101339036A (en) * | 2008-08-20 | 2009-01-07 | 北京航空航天大学 | Terrain auxiliary navigation method and apparatus |
| CN101572017A (en) * | 2009-06-01 | 2009-11-04 | 民航数据通信有限责任公司 | Method for extracting cruise section of flight path from air craft three dimensional position sequence |
| CN101796546A (en) * | 2007-10-30 | 2010-08-04 | 国立大学法人横滨国立大学 | Interpolation processing method and interpolation processor |
| CN103888975A (en) * | 2012-12-24 | 2014-06-25 | 中国移动通信集团河北有限公司 | Base station latitude and longitude data verification method and system |
| CN104504940A (en) * | 2015-01-07 | 2015-04-08 | 江苏理工学院 | Control method of air traffic control system |
| CN105227400A (en) * | 2015-09-16 | 2016-01-06 | 大连理工大学 | A kind of for time become the adaptive estimation method of ADS-B message time delay |
| CN105270642A (en) * | 2014-07-09 | 2016-01-27 | 霍尼韦尔国际公司 | Systems and methods for displaying degraded intruder traffic data on an aircraft display |
-
2016
- 2016-10-25 CN CN201610931969.XA patent/CN106530838B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1020778A2 (en) * | 1999-01-18 | 2000-07-19 | Honeywell Inc. | Flight management system (FMS) |
| WO2004006159A1 (en) * | 2002-07-08 | 2004-01-15 | Proteome Systems Intellectuel Property Pty Ltd | Method and system for picking peaks for mass spectra |
| CN101796546A (en) * | 2007-10-30 | 2010-08-04 | 国立大学法人横滨国立大学 | Interpolation processing method and interpolation processor |
| CN101339036A (en) * | 2008-08-20 | 2009-01-07 | 北京航空航天大学 | Terrain auxiliary navigation method and apparatus |
| CN101572017A (en) * | 2009-06-01 | 2009-11-04 | 民航数据通信有限责任公司 | Method for extracting cruise section of flight path from air craft three dimensional position sequence |
| CN103888975A (en) * | 2012-12-24 | 2014-06-25 | 中国移动通信集团河北有限公司 | Base station latitude and longitude data verification method and system |
| CN105270642A (en) * | 2014-07-09 | 2016-01-27 | 霍尼韦尔国际公司 | Systems and methods for displaying degraded intruder traffic data on an aircraft display |
| CN104504940A (en) * | 2015-01-07 | 2015-04-08 | 江苏理工学院 | Control method of air traffic control system |
| CN105227400A (en) * | 2015-09-16 | 2016-01-06 | 大连理工大学 | A kind of for time become the adaptive estimation method of ADS-B message time delay |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106530838A (en) | 2017-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105894862B (en) | A kind of air traffic control intelligence command system | |
| US11024181B2 (en) | Systems and methods for generating avionic displays including forecast overpressure event symbology | |
| CN107808552B (en) | Flight behavior detection method and device | |
| CN105931497B (en) | Navigation on-air collision detection method, device and all purpose aircraft | |
| US8560148B2 (en) | Method and apparatus for air traffic trajectory synchronization | |
| US9171473B1 (en) | Method and system for dynamic automated corrections to weather avoidance routes for aircraft in en route airspace | |
| US11410558B2 (en) | Traffic control with reinforcement learning | |
| US10678265B2 (en) | Revised speed advisory for an aircraft during flight based on holding time | |
| CN102069915B (en) | Improved stabilized approach monitor | |
| US6420993B1 (en) | Air traffic control system | |
| CN104077928B (en) | Be applicable to the navigation of irregular flight outside air route dot position information update method | |
| EP2837914A1 (en) | Display systems and methods for providing displays indicating a required time of arrival | |
| US20180061245A1 (en) | Community Noise Management with Aircraft Dynamic Path Variation | |
| CN105303893A (en) | A system and method for providing enhanced HMI navigation | |
| CN103294062A (en) | Method for flying an aircraft along a flight path | |
| CN110473431A (en) | A kind of conflict prediction method under the height limitation of approach control airspace | |
| US9666082B2 (en) | Method and system for guidance of an aircraft | |
| CN112396872B (en) | Airplane yaw judging method and device based on computer flight plan CFP data and storage medium | |
| US9493249B2 (en) | System and method for providing an aircraft noise advisory during departure and arrival of an aircraft | |
| US10043401B2 (en) | Methods and apparatus for detecting airport terminal area congestion | |
| US10147327B2 (en) | Method for integrating a constrained route(s) optimization application into an avionics onboard system with open architecture of client server type | |
| EP3407332B1 (en) | System and method for determining uncertainty in a predicted flight path for an aerial vehicle | |
| CN110400061A (en) | A kind of method, apparatus, controller and the storage medium of flight safety comprehensive assessment | |
| CN107209270A (en) | equipment for monitoring position of aircraft | |
| US9406235B2 (en) | Runway location determination |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 230001 Anhui City, Hefei Province, No. 4 Suzhou Road, Xin building, room 3009 Applicant after: Flying friend Technology Co., Ltd. Address before: 230001 Anhui City, Hefei Province, No. 4 Suzhou Road, Xin building, room 3009 Applicant before: Hefei Feiyou Network Technology Co., Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP02 | Change in the address of a patent holder |
Address after: Room 813, Kechuang building, 860 Wangjiang West Road, hi tech Zone, Hefei City, Anhui Province Patentee after: VARIFLIGHT TECHNOLOGY Co.,Ltd. Address before: 3009, room 4, Xin Xin building, No. 230001, Suzhou Road, Hefei, Anhui Patentee before: VARIFLIGHT TECHNOLOGY Co.,Ltd. |
|
| CP02 | Change in the address of a patent holder |