CN116541799B - Channel structure data precision improving method and system, electronic equipment and storage medium - Google Patents
Channel structure data precision improving method and system, electronic equipment and storage medium 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
Abstract
The application discloses a channel structure data precision improving method, a system, electronic equipment and a storage medium, wherein the method comprises the following steps: s1, acquiring QAR data and extracting flight path data and flight performance data related to the flight of an aircraft; s2, converting longitude and latitude coordinate data of an entrance of a flight landing runway and a position of a course platform into a plane rectangular coordinate system through horizontal-axis mercator projection; s3, carrying out interpolation processing on coordinates of two adjacent sampling points of the longitude and latitude coordinate data of the track; s4, calculating the vertical distance Di from the sampling point to the runway extension line and the distance Di from the vertical intersection point of the sampling point and the runway extension line to the course platform on a plane rectangular coordinate system; s5, carrying out data corresponding comparison and correction fusion processing on the modulation degree difference DDMi and the modulation degree difference DDM extracted by the acquisition module. After the processing and correction are carried out on the track data and the channel structure data, the precision and the accuracy of the channel structure data are greatly improved.
Description
Technical Field
The present application relates to the field of channel structure data processing, and in particular, to a method, a system, an electronic device, and a storage medium for improving channel structure data accuracy.
Background
The instrument landing system is the most important navigation equipment in the last approach landing stage of the airport, the running state of the instrument landing system directly influences the landing and airport running standards of the airplane, and the instrument landing system plays an important role in the running and safety of the airport; the device consists of a course table, a lower sliding table and a range finder table. In recent years, an event that a pilot reflects an unstable course signal in the operation of a large busy airport flight occurs, and an unsafe event such as approaching or flying is stopped due to the unstable course signal, so that huge pressure is brought to airport operators. Once unsafe events such as approaching or flying off are stopped at the airport, airport operators need to trace back and analyze reasons of the unsafe events.
At present, the main modes adopted for checking unsafe retrospective events and reasons are checking equipment faults, checking peripheral interference sources or checking through actual flight verification. The event cause investigation time is long, a great amount of manpower and material resources are spent, meanwhile, the cause investigation efficiency is reduced, and huge risks and threats are brought to the safe operation of the airport. Because the flight QAR data (QAR, english full name Quick Access Recorder) is a quick access recorder, various data in the running process of an aircraft are collected and recorded through airborne equipment, the recording capacity is 128MB, the continuous recording time is 600h, hundreds of data can be collected at the same time, the data information closely related to the running process of the aircraft, such as various parameters, aircraft position information, unit operation and the like, can be analyzed through ground decoding software, and the flight state and condition of the aircraft can be analyzed), and the flight QAR data collecting device has the advantages of rich data content, high-dimensional time sequence characteristic, regular sampling interval and the like, and can be used for acquiring the information of the related data of the course signals received by the airborne equipment in the flight process of the aircraft through analyzing the flight QAR data. Due to various reasons such as flight interference and disturbance, accuracy and data acquisition modes of airborne equipment, performance and accuracy of the OAR recording equipment, the accuracy of flight OAR data is relatively low, so that high-accuracy channel structure data cannot be obtained to influence analysis and judgment of channel structure data results. Based on the method, the precision of channel structure data (the modulation degree difference data is at the core) is improved, the corrected modulation degree difference data with higher precision is conveniently obtained, and the efficiency and the accuracy of tracing unsafe events and searching reasons in the flight process of the aircraft can be effectively improved.
Disclosure of Invention
The application aims to overcome the technical problems pointed out by the background art and provides a channel structure data precision improving method, a system, electronic equipment and a storage medium.
The aim of the application is achieved by the following technical scheme:
a channel structure data precision improving method comprises the following steps:
s1, acquiring QAR data, analyzing, decoding and filtering, and extracting flight path data and flight performance data related to the flight of an aircraft, wherein the flight path data comprise flight path longitude and latitude coordinate data, and the flight performance data comprise the ground speed, airspeed, horizontal acceleration and true heading of the aircraft:
s2, collecting longitude and latitude of entrance of flight landing runwayThe coordinate data and the longitude and latitude coordinate data of the position of the course platform are converted into a plane rectangular coordinate system by the horizontal axis mercator projection, the edge line of the sector of the course is drawn on the plane rectangular coordinate system, and the longitude and latitude coordinate data of the entrance of the landing runway of the flight is recorded as (x) 0 ，y 0 ) Longitude and latitude coordinate data of the position of the course table is recorded as (x) L ，y L ) The longitude and latitude coordinate data of the track is recorded as (x) j ，y j )，i＝1，2，3，……，N；
S3, carrying out interpolation processing on coordinates of two adjacent sampling points of the longitude and latitude coordinate data of the track:
the method for obtaining the new sampling point by interpolation through a double integration method according to the flight performance data comprises the following steps:
wherein (x' j ，y′ j ) Representing that the coordinate of the jth sampling point is interpolated by a double integration method to obtain a new sampling point coordinate, (x) a ，y a ) Representing the coordinates of a known point A, wherein the known point A is a known sampling point or the position of an entry or a course platform of a flight landing runway; v a Represents the ground speed of point A, a j Represents the horizontal acceleration, θ, of the jth sampling point j Representing the true heading of the jth sampling point, t j Time t representing the j-th sampling point j+1 The time of the j+1th sampling point;
inserting N in turn between all adjacent two sampling point coordinates 0 Combining the sampling points into new track longitude and latitude coordinate data, and simultaneously converting the new sampling points into a plane rectangular coordinate system through the horizontal axis mercator projection;
s4, in planeCalculating the vertical distance d from each sampling point to the runway extension line on a rectangular coordinate system i And the distance D from the vertical intersection point of each sampling point and the runway extension line to the course platform i The calculation formula is as follows:
wherein k is 0 Representing the slope of the runway extension line, c 0 Representing the intercept of a straight line segment from the entrance of a flight landing runway to the location of a course on the runway extension; (X) i ，Y i ) Representing the coordinates of a sampling point i in the new track longitude and latitude coordinate data; if the sampling point i is positioned at the left side of the runway extension line, the vertical distance d i Positive and noted as d' i ，d′ i ＝d i The method comprises the steps of carrying out a first treatment on the surface of the If the sampling point i is positioned on the right side of the runway extension line, the vertical distance d i Negative and noted as d' i ，d′ i ＝-d i ；
Wherein k is i Represents the slope of the straight line where the vertical intersection point of the sampling point i and the runway extension line is located, c i Representing the intercept of a straight line segment between a sampling point i and a vertical intersection point of a runway extension line;
the vertical distance d of the sampling point i i Distance D from i Corresponding to the association storage, calculating and obtaining the modulation degree difference DDM of the sampling point i according to the following formula i ：
Wherein W is i And/2, the distance from the vertical intersection point of the sampling point i and the runway extension line to the boundary line of the runway sector;
thereby calculating and obtaining all sampling points along with the distance D i Modulation degree difference DDM of (a) i Thereby obtaining the modulation degree difference DDM of channel structure data i A data collection;
s5, extracting a modulation degree difference DDM from OAR data, setting a threshold F, and carrying out data comparison and correction fusion processing according to the following method;
s51, if modulation degree is poor DDM i Modulation degree difference DDM in data set i The data positive-negative difference exists between the extracted modulation degree difference DDM and the data negative-positive difference, and the modulation degree difference DDM is reserved i And deleting the extracted modulation degree difference DDM;
s52, if the modulation degree is poor DDM i Modulation degree difference DDM in data set i The difference value between the extracted modulation degree difference DDM and the extracted modulation degree difference DDM is smaller than or equal to a threshold value F, and the extracted modulation degree difference DDM is adopted to replace the corresponding modulation degree difference DDM i ；
S53, if modulation degree is poor DDM i Modulation degree difference DDM in data set i The difference value between the extracted modulation degree difference DDM and the extracted modulation degree difference DDM is larger than a threshold value F, and the reserved modulation degree difference DDM is adopted i And deleting the extracted modulation degree difference DDM;
thereby obtaining the corrected modulation degree difference DDM i And (5) data collection.
In order to better realize the channel structure data precision improving method of the application, the sampling time interval of the track data in the QAR data is T 0 The method comprises the steps of carrying out a first treatment on the surface of the In step S3, if the interval time between two adjacent sampling point coordinates is less than or equal to T 0 N is inserted between two adjacent sampling point coordinates in sequence according to the flight performance data by a double integration method 0 Sampling points; if the interval time between the coordinates of two adjacent sampling points is greater than T 0 Sequentially inserting at least 2n between two adjacent sampling point coordinates according to the flight performance data by a double integration method 0 And sampling points.
The channel structure data precision improving method of the application adopts the preferable technical scheme that: in step S4, the slope k of the runway extension line 0 The method is calculated by the following formula:
on the extension line of the runwayIntercept c of straight line segment from entrance of flight landing runway to position of course 0 The method is calculated by the following formula:
the channel structure data precision improving method of the application adopts the preferable technical scheme that: slope k of straight line where vertical intersection point of sampling point i and runway extension line is located i The method is calculated by the following formula:
intercept c of straight line segment between sampling point i and runway extension line vertical intersection point i The method is calculated by the following formula:
a channel structure data accuracy improvement system, comprising:
the acquisition module is used for acquiring OAR data, analyzing, decoding and filtering, extracting flight path data and flight performance data related to the flight of the aircraft, wherein the flight path data comprise flight path longitude and latitude coordinate data, and the flight performance data comprise the ground speed, airspeed, horizontal acceleration and true heading of the aircraft: the acquisition module is also used for extracting a modulation degree difference DDM from the OAR data and acquiring longitude and latitude coordinate data of an entrance of a flight landing runway and longitude and latitude coordinate data of a position of a course platform;
the computing processing system comprises an interpolation computing module, a modulation degree difference computing module and a modulation degree difference correction fusion module, wherein a plane rectangular coordinate system is built in the computing processing system, longitude and latitude coordinate data of an entrance of a flight landing runway, longitude and latitude coordinate data of a position of a course platform are converted into the plane rectangular coordinate system through a horizontal axis mercator projection, and an edge line of a channel sector is drawn on the plane rectangular coordinate system; interpolationThe calculation module carries out interpolation processing on coordinates of two adjacent sampling points of the track longitude and latitude coordinate data, and N is sequentially inserted between the coordinates of all the two adjacent sampling points 0 Combining the sampling points into new track longitude and latitude coordinate data, and simultaneously converting the new sampling points into a plane rectangular coordinate system through the horizontal axis mercator projection; the modulation degree difference calculation module calculates all sampling points along with the distance D i Modulation degree difference DDM of (a) i And obtains the modulation degree difference DDM of channel structure data i A data collection; the modulation degree difference correction fusion module sets a threshold F and sets a modulation degree difference DDM i Data corresponding comparison is carried out on the data and the modulation degree difference DDM extracted by the acquisition module, and fusion processing is corrected, so that the corrected modulation degree difference DDM is obtained i A data collection;
an output module for outputting the corrected modulation degree difference DDM i Data.
An electronic device, comprising: at least one processor; and a memory communicatively coupled to the at least one processor; the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor, so that the at least one processor executes the steps of the channel structure data precision improving method.
A storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the channel structure data accuracy improvement method of the present application.
Compared with the prior art, the application has the following advantages:
(1) The application greatly improves the precision and accuracy of the channel structure data by extracting the data related to the course platform recorded by the receiving end of the flight at the near-flight stage and processing and correcting the channel structure data, and can obtain the high-precision modulation degree difference DDM i The data provides powerful accurate data support for airport related personnel in the event investigation process, and is convenient for realizing signal evaluation and event effective investigation of an instrument landing system.
(2) The application converts the track longitude and latitude data into a plane rectangular coordinate system, creatively carries out interpolation processing on the sampling points of the track data through double integration, obtains the vertical distance between the sampling points (or called track points) and the track center line and the equidistant distance between the course platform and the intersection point of two straight lines, has simple method and accurate result, can quickly realize correction processing on the track structure data, and improves the precision of the track structure data.
(3) The application can directly analyze the channel structure data of the course platform through the flight OAR data, saves the time cost of manual signal measurement and improves the efficiency of checking unsafe events; the accuracy of the channel structure data is improved and the accuracy is high by correcting the data related to the channel structure data of the course platform in the QAR data; the method provides accurate data support for checking and analyzing the instability reasons of the course table signals, provides technical support for safe operation of the airport aircraft, and has the advantages of high data production efficiency, high quality, reliable data and the like.
Drawings
FIG. 1 is a flow chart of a method for improving the precision of channel structure data according to the present application;
FIG. 2 is a flow chart of the method for comparing and correcting the fusion process of data in step S5 of the present application;
fig. 3 is a schematic diagram of a relationship of a rectangular planar coordinate system taking a sampling point i as an example in the embodiment;
FIG. 4 is a schematic diagram showing the comparison of the channel structure data before and after correction in the embodiment;
FIG. 5 is a functional block diagram of a channel structure data accuracy enhancement system of the present application.
Detailed Description
The application is further illustrated by the following examples:
examples
As shown in fig. 1 to 4, a channel structure data precision improving method includes:
s1, acquiring OAR data (the acquired OAR data comprises flight number, flight date, longitude and latitude, flight altitude, horizontal acceleration, ground speed, airspeed, true heading and magnetic heading)Related data such as a DDM value of a course platform, a working frequency of the course platform, a flight landing runway number and the like; because the flight OAR data is original code data, the OAR data is converted into readable data after being analyzed and decoded; the sampling time interval of track data in OAR data is T 0 The method comprises the steps of carrying out a first treatment on the surface of the The QAR data used in this embodiment has a sampling time interval of 1 second), and then analyzing, decoding, and filtering, extracting flight path data and flight performance data related to the flight of the aircraft, where the flight path data includes flight path longitude and latitude coordinate data, and the flight performance data includes an aircraft ground speed, an airspeed, a horizontal acceleration, and a true heading:
s2, collecting longitude and latitude coordinate data of an entrance of a flight landing runway and longitude and latitude coordinate data of a position of a course platform, converting the longitude and latitude coordinate data of the entrance of the flight landing runway, the longitude and latitude coordinate data of the position of the course platform and the longitude and latitude coordinate data of the position of the course platform into a plane rectangular coordinate system through a horizontal axis ink-card-based projection, and drawing an edge line of a channel sector (an area limited by a track of points extending from two sides of the channel line to 0.155 (150 mu A) DDM at 105m at two sides of the entrance of the channel is 0.155 (150 mu A), for example a) s For the runway sector angle, the values of the runway sector angle can be calculated according to the distance from the runway platform to the runway end, the runway length and the widths of the two sides of the runway entrance, and the longitude and latitude coordinate data of the landing runway entrance of the flight is recorded as (x) 0 ，y 0 ) See the example of fig. 3; longitude and latitude coordinate data of the position of the course table is recorded as (x) L ，y L ) See the example of fig. 3; the longitude and latitude coordinate data of the track is recorded as (x) j ，y j )，i＝1，2，3，……，N。
And S3, carrying out interpolation processing on coordinates of two adjacent sampling points of the longitude and latitude coordinate data of the track.
The method for obtaining the new sampling point by interpolation through a double integration method according to the flight performance data comprises the following steps:
wherein (x' j ，y′ j ) Representing that the coordinate of the jth sampling point is interpolated by a double integration method to obtain a new sampling point coordinate, (x) a ，y a ) Representing the coordinates of a known point A, wherein the known point A is a known sampling point or the position of an entry or a course platform of a flight landing runway; v a Represents the ground speed of point A, a j Represents the horizontal acceleration, θ, of the jth sampling point j Representing the true heading of the jth sampling point, t j Time t representing the j-th sampling point j+1 The time of the j+1th sampling point is indicated.
Inserting N in turn between all adjacent two sampling point coordinates 0 And combining the sampling points into new track longitude and latitude coordinate data, and simultaneously converting the new sampling points into a plane rectangular coordinate system through the horizontal axis mercator projection. Normally, track data collected in QAR data will be at sampling time intervals T 0 Collecting track sampling points, namely, the interval time of coordinates of any two sampling points in track data is equal to T 0 . Because the original data can have abnormal conditions such as dislocation of partial data fields, information deletion, parameter value exceeding a theoretical value range, occurrence of logic-free jump and the like in the decoding and analyzing process, the abnormal data are required to be identified and deleted by combining the parameter data related to the aircraft state in a period of time near the time point where the abnormal data are located.
If the interval time between two adjacent sampling point coordinates is less than or equal to T 0 N is inserted between two adjacent sampling point coordinates in sequence according to the flight performance data by a double integration method 0 The number of sampling points (i.e. N 0 ＝n 0 The two adjacent sampling points can be two adjacent sampling points in the track data, a new sampling point can be inserted between one sampling point and the adjacent sampling points in the track data, and a new sampling point can be inserted between the two adjacent sampling points. If the interval time between the coordinates of two adjacent sampling points is greater than T 0 Then, the two products are passed between the coordinates of two adjacent sampling points according to the flight performance dataSequentially inserting at least 2n by a dividing method 0 The number of sampling points (i.e. N 0 ＝2n 0 ). The above mentioned sampling points or insertion sampling points, i.e. track points in the track data.
S4, as shown in FIG. 3, calculating the vertical distance d from each sampling point to the runway extension line on a plane rectangular coordinate system i (see FIG. 3 for an example of the distance between sampling point i and course, and the distance D between the vertical intersection of each sampling point and the course extension and the course i (see fig. 3 for an example of distances for sampling point i) the calculation formula is as follows:
wherein k is 0 Representing the slope of the runway extension line, c 0 Representing the intercept of a straight line segment from the entrance of a flight landing runway to the location of a course on the runway extension; (X) i ，Y i ) Representing the coordinates of a sampling point i in the new track longitude and latitude coordinate data; if the sampling point i is positioned at the left side of the runway extension line, the vertical distance d i Positive and noted as d' i ，d′ i ＝d i The method comprises the steps of carrying out a first treatment on the surface of the If the sampling point i is positioned on the right side of the runway extension line, the vertical distance d i Negative and noted as d' i ，d′ i ＝-d i 。
Slope k of runway extension line 0 The method is calculated by the following formula:
intercept c of straight line segment from entrance of flight landing runway to position of course platform on runway extension line 0 The method is calculated by the following formula:
wherein k is i Represents the slope of the straight line where the vertical intersection point of the sampling point i and the runway extension line is located, c i Representing the intercept of the straight line segment between the vertical intersection point of the sampling point i and the runway extension line. The DDM of the runway extension line closest to the runway center in the horizontal plane is zero point set, namely the extension line of the runway center line.
Slope k of straight line where vertical intersection point of sampling point i and runway extension line is located i The method is calculated by the following formula:
intercept c of straight line segment between sampling point i and runway extension line vertical intersection point i The method is calculated by the following formula:
the vertical distance d of the sampling point i i Distance D from i Corresponding to the association storage, calculating and obtaining the modulation degree difference DDM of the sampling point i according to the following formula i ：
Wherein W is i And/2, the distance from the vertical intersection point of the sampling point i and the runway extension line to the boundary line of the runway sector;
thereby calculating and obtaining all sampling points along with the distance D i Modulation degree difference DDM of (a) i Thereby obtaining the modulation degree difference DDM of channel structure data i A data collection;
s5, extracting a modulation degree difference DDM (namely the modulation degree difference DDM extracted by OAR data, specifically, the difference of the modulation degree percentages of two audio signals of 150Hz and 90Hz to radio frequency) from the OAR data, setting a threshold value F, and carrying out data comparison and correction fusion processing according to the following method (see figure 2);
s51, if modulation degree is poor DDM i Modulation degree difference DDM in data set i The data positive-negative difference exists between the extracted modulation degree difference DDM and the data negative-positive difference, and the modulation degree difference DDM is reserved i And deleting the extracted modulation degree difference DDM;
s52, if the modulation degree is poor DDM i Modulation degree difference DDM in data set i The difference value between the extracted modulation degree difference DDM and the extracted modulation degree difference DDM is smaller than or equal to a threshold value F, and the extracted modulation degree difference DDM is adopted to replace the corresponding modulation degree difference DDM i ；
S53, if modulation degree is poor DDM i Modulation degree difference DDM in data set i The difference value between the extracted modulation degree difference DDM and the extracted modulation degree difference DDM is larger than a threshold value F, and the reserved modulation degree difference DDM is adopted i And deleting the extracted modulation degree difference DDM;
thereby obtaining the corrected modulation degree difference DDM i And (5) data collection. Taking course station channel structure data of I-type operation of an airport and OAR data of landing flights of the airport as examples, the corrected modulation degree difference DDM i The DDM values directly read from the data and OAR data are compared, the comparison curve is shown in FIG. 4, and the corrected modulation degree difference DDM can be seen from the graph i The accuracy of the data (namely the channel structure data) is obviously improved.
As shown in fig. 5, a channel structure data precision improving system includes:
the acquisition module is used for acquiring OAR data, analyzing, decoding and filtering, extracting flight path data and flight performance data related to the flight of the aircraft, wherein the flight path data comprise flight path longitude and latitude coordinate data, and the flight performance data comprise the ground speed, airspeed, horizontal acceleration and true heading of the aircraft: the acquisition module is also used for extracting a modulation degree difference DDM from the OAR data and acquiring longitude and latitude coordinate data of an entrance of a flight landing runway and longitude and latitude coordinate data of a position of a course platform;
the calculation processing system comprises an interpolation calculation module, a modulation degree difference calculation module and a modulation degree difference correction fusion module, wherein a plane rectangular coordinate system is built in the calculation processing system, and longitude and latitude coordinates of an entrance of a flight landing runway are obtainedConverting the data, longitude and latitude coordinate data of the position of the course platform and longitude and latitude coordinate data of the track into a plane rectangular coordinate system through the horizontal axis ink-card support projection, and drawing an edge line of the sector of the track on the plane rectangular coordinate system; the interpolation calculation module carries out interpolation processing on coordinates of two adjacent sampling points of the track longitude and latitude coordinate data, and N is sequentially inserted between the coordinates of all the two adjacent sampling points 0 Combining the sampling points into new track longitude and latitude coordinate data, and simultaneously converting the new sampling points into a plane rectangular coordinate system through the horizontal axis mercator projection; the modulation degree difference calculation module calculates all sampling points along with the distance D i Modulation degree difference DDM of (a) i And obtains the modulation degree difference DDM of channel structure data i A data collection; the modulation degree difference correction fusion module sets a threshold F and sets a modulation degree difference DDM i Data corresponding comparison is carried out on the data and the modulation degree difference DDM extracted by the acquisition module, and fusion processing is corrected, so that the corrected modulation degree difference DDM is obtained i A data collection;
an output module for outputting the corrected modulation degree difference DDM i Data.
An electronic device, comprising: at least one processor; and a memory communicatively coupled to the at least one processor; the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor, so that the at least one processor executes the steps of the channel structure data precision improving method.
A storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the channel structure data accuracy improvement method of the present application.
The foregoing description of the preferred embodiments of the application is not intended to be limiting, but rather is intended to cover all modifications, equivalents, and alternatives falling within the spirit and principles of the application.
Claims (7)
1. A channel structure data precision improving method is characterized in that: the method comprises the following steps:
s1, acquiring QAR data, analyzing, decoding and filtering, and extracting flight path data and flight performance data related to the flight of an aircraft, wherein the flight path data comprise flight path longitude and latitude coordinate data, and the flight performance data comprise the ground speed, airspeed, horizontal acceleration and true heading of the aircraft:
s2, acquiring longitude and latitude coordinate data of an entrance of a flight landing runway and longitude and latitude coordinate data of a position of a course platform, converting the longitude and latitude coordinate data of the entrance of the flight landing runway, the longitude and latitude coordinate data of the position of the course platform and the longitude and latitude coordinate data of the course platform into a plane rectangular coordinate system through a horizontal axis mercator projection, drawing an edge line of a sector of the course on the plane rectangular coordinate system, and recording the longitude and latitude coordinate data of the entrance of the flight landing runway as (x) 0 ，y 0 ) Longitude and latitude coordinate data of the position of the course table is recorded as (x) L ，y L ) The longitude and latitude coordinate data of the track is recorded as (x) j ，y j )，i＝1，2，3，……，N；
S3, carrying out interpolation processing on coordinates of two adjacent sampling points of the longitude and latitude coordinate data of the track:
the method for obtaining the new sampling point by interpolation through a double integration method according to the flight performance data comprises the following steps:
wherein (x' j ，y′ j ) Representing that the coordinate of the jth sampling point is interpolated by a double integration method to obtain a new sampling point coordinate, (x) a ，y a ) Representing the coordinates of a known point A, wherein the known point A is a known sampling point or the position of an entry or a course platform of a flight landing runway; v a Represents the ground speed of point A, a j Represents the horizontal acceleration, θ, of the jth sampling point j Representing the true heading of the jth sample point,t j time t representing the j-th sampling point j+1 The time of the j+1th sampling point;
inserting N in turn between all adjacent two sampling point coordinates 0 Combining the sampling points into new track longitude and latitude coordinate data, and simultaneously converting the new sampling points into a plane rectangular coordinate system through the horizontal axis mercator projection;
s4, calculating the vertical distance d from each sampling point to the runway extension line on a plane rectangular coordinate system i And the distance D from the vertical intersection point of each sampling point and the runway extension line to the course platform i The calculation formula is as follows:
wherein k is 0 Representing the slope of the runway extension line, c 0 Representing the intercept of a straight line segment from the entrance of a flight landing runway to the location of a course on the runway extension; (X) i ，Y i ) Representing the coordinates of a sampling point i in the new track longitude and latitude coordinate data; if the sampling point i is positioned at the left side of the runway extension line, the vertical distance d i Positive and noted as d' i ，d′ i ＝d i The method comprises the steps of carrying out a first treatment on the surface of the If the sampling point i is positioned on the right side of the runway extension line, the vertical distance d i Negative and noted as d' i ，d′ i ＝-d i ；
Wherein k is i Represents the slope of the straight line where the vertical intersection point of the sampling point i and the runway extension line is located, c i Representing the intercept of a straight line segment between a sampling point i and a vertical intersection point of a runway extension line;
the vertical distance d of the sampling point i i Distance D from i Corresponding to the association storage, calculating and obtaining the modulation degree difference DDM of the sampling point i according to the following formula i ：
Wherein W is i And/2, the distance from the vertical intersection point of the sampling point i and the runway extension line to the boundary line of the runway sector;
thereby calculating and obtaining all sampling points along with the distance D i Modulation degree difference DDM of (a) i Thereby obtaining the modulation degree difference DDM of channel structure data i A data collection;
s5, extracting a modulation degree difference DDM from QAR data, setting a threshold value F, and carrying out data comparison and correction fusion processing according to the following method;
s51, if the modulation degree difference DDM in the modulation degree difference DDM data set i The data positive-negative difference exists between the extracted modulation degree difference DDM and the data negative-positive difference, and the modulation degree difference DDM is reserved i And deleting the extracted modulation degree difference DDM;
s52, if the modulation degree is poor DDM i Modulation degree difference DDM in data set i The difference value between the extracted modulation degree difference DDM and the extracted modulation degree difference DDM is smaller than or equal to a threshold value F, and the extracted modulation degree difference DDM is adopted to replace the corresponding modulation degree difference DDM i ；
S53, if modulation degree is poor DDM i Modulation degree difference DDM in data set i The difference value between the extracted modulation degree difference DDM and the extracted modulation degree difference DDM is larger than a threshold value F, and the reserved modulation degree difference DDM is adopted i And deleting the extracted modulation degree difference DDM;
thereby obtaining the corrected modulation degree difference DDM i And (5) data collection.
2. The channel structure data precision improving method according to claim 1, wherein: the sampling time interval of track data in QAR data is T 0 The method comprises the steps of carrying out a first treatment on the surface of the In step S3, if the interval time between two adjacent sampling point coordinates is less than or equal to T 0 N is inserted between two adjacent sampling point coordinates in sequence according to the flight performance data by a double integration method 0 Sampling points; if the interval time between the coordinates of two adjacent sampling points is greater than T 0 Sequentially inserting the two adjacent sampling point coordinates according to the flight performance data by a double integration methodAt least 2n 0 And sampling points.
3. The channel structure data precision improving method according to claim 1, wherein: in step S4, the slope k of the runway extension line 0 The method is calculated by the following formula:
intercept c of straight line segment from entrance of flight landing runway to position of course platform on runway extension line 0 The method is calculated by the following formula:
4. the channel structure data precision improving method according to claim 1, wherein: slope k of straight line where vertical intersection point of sampling point i and runway extension line is located i The method is calculated by the following formula:
intercept c of straight line segment between sampling point i and runway extension line vertical intersection point i The method is calculated by the following formula:
5. a channel structure data precision improving system is characterized in that: comprising the following steps:
the acquisition module is used for acquiring QAR data, analyzing, decoding and filtering, extracting flight path data and flight performance data related to the flight of the aircraft, wherein the flight path data comprises flight path longitude and latitude coordinate data, and the flight performance data comprises the ground speed, airspeed, horizontal acceleration and true heading of the aircraft: the acquisition module is also used for extracting a modulation degree difference DDM from the QAR data and acquiring longitude and latitude coordinate data of an entrance of a flight landing runway and longitude and latitude coordinate data of a position of a course platform;
the computing processing system comprises an interpolation computing module, a modulation degree difference computing module and a modulation degree difference correction fusion module, wherein a plane rectangular coordinate system is built in the computing processing system, longitude and latitude coordinate data of an entrance of a flight landing runway, longitude and latitude coordinate data of a position of a course platform are converted into the plane rectangular coordinate system through a horizontal axis mercator projection, and an edge line of a channel sector is drawn on the plane rectangular coordinate system; the interpolation calculation module carries out interpolation processing on coordinates of two adjacent sampling points of the track longitude and latitude coordinate data, and N is sequentially inserted between the coordinates of all the two adjacent sampling points 0 Combining the sampling points into new track longitude and latitude coordinate data, and simultaneously converting the new sampling points into a plane rectangular coordinate system through the horizontal axis mercator projection; the modulation degree difference calculation module calculates all sampling points along with the distance D i Modulation degree difference DDM of (a) i And obtains the modulation degree difference DDM of channel structure data i Data set, sampling point along with distance D i Modulation degree difference DDM of (a) i The obtaining method comprises the following steps:
calculating the vertical distance d from each sampling point to the runway extension line on a plane rectangular coordinate system i And the distance D from the vertical intersection point of each sampling point and the runway extension line to the course platform i The calculation formula is as follows:
wherein k is 0 Representing the slope of the runway extension line, c 0 Representing the intercept of a straight line segment from the entrance of a flight landing runway to the location of a course on the runway extension; (X) i ，Y i ) Representing the coordinates of a sampling point i in the new track longitude and latitude coordinate data; if the sampling point is i bitsAt the left side of the runway extension line, the vertical distance d i Positive and noted as d' i ，d′ i ＝d i The method comprises the steps of carrying out a first treatment on the surface of the If the sampling point i is positioned on the right side of the runway extension line, the vertical distance d i Negative and noted as d' i ，d′ i ＝-d i ；
Wherein k is i Represents the slope of the straight line where the vertical intersection point of the sampling point i and the runway extension line is located, c i Representing the intercept of a straight line segment between a sampling point i and a vertical intersection point of a runway extension line;
the vertical distance d of the sampling point i i Distance D from i Corresponding to the association storage, calculating and obtaining the modulation degree difference DDM of the sampling point i according to the following formula i ：
Wherein W is i And/2, the distance from the vertical intersection point of the sampling point i and the runway extension line to the boundary line of the runway sector; the modulation degree difference correction fusion module sets a threshold F and sets a modulation degree difference DDM i Data corresponding comparison is carried out on the data and the modulation degree difference DDM extracted by the acquisition module, and fusion processing is corrected, so that the corrected modulation degree difference DDM is obtained i A data collection;
an output module for outputting the corrected modulation degree difference DDM i Data.
6. An electronic device, characterized in that: comprising the following steps: at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to cause the at least one processor to perform the steps of the method of any of claims 1-4.
7. A storage medium having a computer program stored thereon, characterized by: the computer program implementing the steps of the method according to any of claims 1-4 when executed by a processor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310506127.XA CN116541799B (en) | 2023-05-06 | 2023-05-06 | Channel structure data precision improving method and system, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310506127.XA CN116541799B (en) | 2023-05-06 | 2023-05-06 | Channel structure data precision improving method and system, electronic equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116541799A CN116541799A (en) | 2023-08-04 |
| CN116541799B true CN116541799B (en) | 2023-11-07 |
Family
ID=87453756
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310506127.XA Active CN116541799B (en) | 2023-05-06 | 2023-05-06 | Channel structure data precision improving method and system, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116541799B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110456655A (en) * | 2019-06-29 | 2019-11-15 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | Microwave landing analogue simulation system |
| CN114754775A (en) * | 2022-04-07 | 2022-07-15 | 广州嘉航通信科技有限公司 | Channel change analysis early warning method, system, computer equipment and storage medium |
| WO2022257620A1 (en) * | 2021-06-08 | 2022-12-15 | 东航技术应用研发中心有限公司 | Method and system for measuring in-flight landing distance |
| CN115798265A (en) * | 2022-11-01 | 2023-03-14 | 航科院中宇（北京）新技术发展有限公司 | Digital tower construction method based on digital twinning technology and implementation system thereof |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6921339B1 (en) * | 2020-06-11 | 2021-08-18 | 三菱電機株式会社 | Radar device and radar image generation method |
-
2023
- 2023-05-06 CN CN202310506127.XA patent/CN116541799B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110456655A (en) * | 2019-06-29 | 2019-11-15 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | Microwave landing analogue simulation system |
| WO2022257620A1 (en) * | 2021-06-08 | 2022-12-15 | 东航技术应用研发中心有限公司 | Method and system for measuring in-flight landing distance |
| CN114754775A (en) * | 2022-04-07 | 2022-07-15 | 广州嘉航通信科技有限公司 | Channel change analysis early warning method, system, computer equipment and storage medium |
| CN115798265A (en) * | 2022-11-01 | 2023-03-14 | 航科院中宇（北京）新技术发展有限公司 | Digital tower construction method based on digital twinning technology and implementation system thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116541799A (en) | 2023-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110515116B (en) | Unmanned aerial vehicle aviation gamma energy spectrum detection and data processing system and processing method thereof | |
| Iliff et al. | Practical aspects of a maximum likelihood estimation method to extract stability and control derivatives from flight data | |
| CN107944701B (en) | Method and device for detecting runway rushing-out risk in aircraft landing process | |
| CN113205706A (en) | ILS signal quality monitoring method based on flight QAR data | |
| CN113487862A (en) | Green wave speed determination method and device, electronic equipment and storage medium | |
| CN116541799B (en) | Channel structure data precision improving method and system, electronic equipment and storage medium | |
| CN115017578A (en) | Intelligent actual measurement method and device for building, UGV and storage medium | |
| CN108073865B (en) | Aircraft trail cloud identification method based on satellite data | |
| CN112542800A (en) | Method and system for identifying transmission line fault | |
| CN116774191A (en) | POS data and machine-mounted laser radar point cloud data interaction quality inspection device | |
| CN114782847B (en) | Mine productivity monitoring method and device based on unmanned aerial vehicle | |
| CN103837130B (en) | For data processing method and the device of airborne lidar system | |
| CN113190040B (en) | Method and system for line inspection based on unmanned aerial vehicle video and railway BIM | |
| CN114035189A (en) | Road hidden danger positioning method and device, electronic equipment and storage medium | |
| CN113311858A (en) | Flight control system and method based on unmanned aerial vehicle | |
| CN107588760B (en) | Fixed-height high-precision track generation method suitable for fixed-wing unmanned aerial vehicle tracking road | |
| CN111122900B (en) | Magnetic difference correction method and device | |
| CN115060457B (en) | Method, system and equipment for detecting atmospheric vortex dissipation rate based on aircraft bump | |
| Behere et al. | A Method for the Parametric Representation of Take-off Time-Series Trajectory Data for Environmental Impact Assessment | |
| CN116453377B (en) | Method for carrying out flight phase division on airplane QAR data | |
| CN117612047B (en) | Unmanned aerial vehicle inspection image recognition method for power grid based on AI large model | |
| CN115359686B (en) | Instrument landing system signal analysis and diagnosis method and device and electronic equipment | |
| CN114115348B (en) | Switching method and device of inspection device of unmanned aerial vehicle, electronic equipment and storage medium | |
| CN115980676B (en) | Radar signal data analysis system and method based on big data | |
| CN116243342A (en) | GNSS occultation ionosphere profile file quality analysis method and system |
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 |