CN112214537A - Track characteristic modeling method and system for track retrieval and electronic device - Google Patents
Track characteristic modeling method and system for track retrieval and electronic device Download PDFInfo
- Publication number
- CN112214537A CN112214537A CN202010975489.XA CN202010975489A CN112214537A CN 112214537 A CN112214537 A CN 112214537A CN 202010975489 A CN202010975489 A CN 202010975489A CN 112214537 A CN112214537 A CN 112214537A
- Authority
- CN
- China
- Prior art keywords
- area
- flight
- track
- region
- earth
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 238000004590 computer program Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000005484 gravity Effects 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 4
- 238000011156 evaluation Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005192 partition Methods 0.000 description 3
- 241000337007 Oceania Species 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000007596 consolidation process Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/248—Presentation of query results
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computational Linguistics (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
Abstract
The invention belongs to the technical field of aircraft track modeling, and aims to solve at least one technical problem that in the prior art, track retrieval has insufficient precision and processing speed cannot meet real-time requirements; then, a multi-stage earth position dictionary is utilized to quickly search similar tracks; because the earth position code adopts a multi-level algorithm, a flight path set meeting the conditions can be preliminarily screened and found by a gradual thinning method, then a thinned regional dictionary is gradually used, and finally the most similar flight path is accurately positioned; the track retrieval efficiency is greatly improved.
Description
Technical Field
The invention relates to the technical field of aircraft trajectory modeling, in particular to a flight path characteristic modeling method and system for flight path retrieval, an electronic device and a nonvolatile storage medium.
Background
The flight path of the aircraft is important data for carrying out target identification, target warning and daily target behavior intention analysis of the aircraft. In most cases, the travel activities of the target are regular each time, so that the target can be identified, warned, analyzed in behavior, supplemented with data and the like by a similar track retrieval technology. However, in practical situations, the acquired real-time flight path is very unstable, and problems of data loss, abnormal data, no cycle of data sampling, inconsistent cycle and the like often occur.
For example, in chinese patent application No. CN201510084399.0, in order to solve the technical problem that the "inertial navigation + GPS" combined navigation system is susceptible to interference, the "inertial navigation + terrain matching" combined navigation system can only be used in an area with specific terrain undulations, a processing method is provided in which a Douglas-Peucker (Douglas-Peucker) algorithm is used to compress the flight path, and then the longest common substring is used for processing. However, these algorithms still exist, either with insufficient accuracy or with insufficient speed to meet real-time requirements.
Therefore, how to design a method that can quickly and accurately retrieve similar tracks under these difficult conditions becomes a technical problem that those skilled in the art need to solve.
Disclosure of Invention
In order to solve at least one technical problem of insufficient precision and processing speed which cannot meet the real-time requirement in the prior art of track retrieval, the invention provides a track characteristic modeling method and system, an electronic device and a nonvolatile storage medium for track retrieval.
The invention provides a track characteristic modeling method for track retrieval, which is characterized by comprising the following steps:
establishing a coordinate system of the earth surface, and determining different earth region division levels according to different precision of the division regions;
respectively establishing earth position dictionary code libraries with different levels of precision according to the difference of the divided region precision; the earth position dictionary code library with each grade precision sets a number for each region in the earth region respectively aiming at the corresponding region precision;
according to the flight path of the aircraft, connecting lines from a flight starting point to a flight terminal point and from points to points of a path region one by one to form a path region chain; and the track area chain comprises: forming corresponding flight tracks by combining a flight starting point, point-by-point of a path region and a flight terminal point in earth position dictionary code libraries based on different grades respectively;
when the flight path of the aircraft needs to be inquired, sequentially searching whether corresponding flight paths exist in the path region chain from a low-precision earth position dictionary code library to a high-precision earth position dictionary code library.
In a preferred implementation manner of the embodiment of the invention, a flight history track area chain is formed based on the flight history of the aircraft; when an aircraft appears in an area to be monitored, according to the flight path of the aircraft, in the flight history path area chain, whether the same or similar flight paths exist in the flight history path area chain is searched through a method of gradually thinning the earth position dictionary code library level.
In a preferred implementation manner of the embodiment of the present invention, determining different earth region division levels according to different accuracy of the divided regions specifically includes:
determining an initialization area parameter of the earth area based on longitudes and latitudes with different sizes;
and merging the initialized areas according to the maximum weighted population number and the maximum area of the areas.
In a further preferred implementation manner of the embodiment of the present invention, determining the initialized area parameter of the earth area includes:
triple for each area
To indicate that the user is not in a normal position,
wherein { (j)1,w1),(j2,w2) Represents the location of the area;
and { (j)1,w1),(j2,w2) The area calculation of the corresponding area comprises:
the average latitude of the area is calculated as:
calculate the area aspect ratio of the area:
calculating the area actual area of the area: s ═ w1-w2|*|j1-j2|*cos((w1+w2)/2)。
In a further preferred implementation manner of the embodiment of the present invention, the merging the initialized areas according to the maximum weighted population number of the areas and the maximum area of the areas specifically includes:
determining world population density according to the population number on a unit land area of a global land surface, and determining the population density corresponding to each region in a global position dictionary code library with different levels of precision;
in the track area chain, the population density of the area where the flight starting point and the flight ending point of each track are positioned is increased by a first preset value, and the population density of each area is increased by a second preset value in the areas where the rest tracks pass through;
calculating an area weighted area (area actual area × (1+ α × log (1+ population density increase predetermined value/area population));
and combining the initialized regions of the earth regions by using a square principle and a gravity center centering algorithm to obtain the combined earth regions.
In a preferred implementation manner of the embodiment of the present invention, the method for establishing the global position dictionary code library includes: setting a number for each region, and then establishing an index by applying a general red-black tree algorithm to form a global position dictionary; any position of the earth region can be directly inquired about the region code through dictionary indexing, and the region information can be directly inquired about through the dictionary according to the region code.
In a preferred implementation manner of the embodiment of the present invention, the track characteristics include a track area chain, a flight time period, a flight duration, and a flight height, and the track block chain includes a flight starting point, a flight ending point, an earth area through which points pass one by one, and an area through which a connection line of two adjacent points passes; the flight time period is used for recording the starting time and the ending time, the flight time length is used for recording the flight time length, and the flight altitude is used for recording the maximum altitude, the maximum altitude and the average altitude; the flight speed is used for recording the maximum speed, the fastest speed and the average speed; the track features are used for executing the same or similar track search, and the track retrieval comprises the steps of carrying out preliminary screening by using the area chain with low precision, and then determining similar tracks by using the area chain with high precision and combining the track area chain, the flight time period, the flight duration and the flight height.
The second aspect of the present invention also provides a track characteristic modeling system for track search, including:
the earth region dividing module is used for creating a coordinate system of the earth surface and determining different earth region dividing grades according to different accuracy of divided regions;
the earth position dictionary code library establishing module is used for respectively establishing earth position dictionary code libraries with different levels of precision according to the difference of the divided region precision; the earth position dictionary code library with each grade precision sets a number for each region in the earth region respectively aiming at the corresponding region precision;
the flight path region chain forming module is used for connecting lines from a flight starting point to a flight terminal point and from points to points of the path region one by one according to the flight path of the aircraft to form a flight path region chain; and the track area chain comprises: forming corresponding flight tracks by combining a flight starting point, point-by-point of a path region and a flight terminal point in earth position dictionary code libraries based on different grades respectively;
when the flight path of the aircraft needs to be inquired, sequentially searching whether corresponding flight paths exist in the path region chain from a low-precision earth position dictionary code library to a high-precision earth position dictionary code library.
The third aspect of the present invention also provides an electronic device, including:
a memory for storing a plurality of data to be transmitted,
a processor, and
a computer program;
wherein the computer program is stored in the memory and configured to be executed by the processor to implement the method for modeling track features for track retrieval as any one of the aspects provided.
The fourth aspect of the present invention also provides a nonvolatile storage medium characterized by having a computer program stored thereon; the computer program is executed by a processor to implement any one of the track feature modeling methods for track retrieval provided by the first aspect.
Therefore, in the technical scheme provided by the application, all positions of the earth are mapped to a certain region code through a multi-stage earth position coding algorithm, and then codes of flight paths of the aircraft passing through the region are concatenated, so that a flight path region chain convenient to identify is formed. In addition, in the process of track identification, as the earth position code adopts a multi-stage algorithm, similar tracks can be searched by a gradual thinning method, and the retrieval efficiency is greatly improved; for example, a coarsest divided regional dictionary is adopted, a flight path set meeting the conditions is found through preliminary screening, then the refined regional dictionary is used step by step, and finally the most similar flight path is accurately positioned.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention may be realized and attained by the structure and/or process particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
Fig. 1 is a flowchart of a track feature modeling method for track search according to an embodiment of the present invention.
Fig. 2 is a flowchart of a method for encoding a global position dictionary according to an embodiment of the present invention.
Fig. 3 is a schematic diagram of a region merging method in an earth location according to an embodiment of the present invention.
Fig. 4 is a schematic diagram of flight path characteristic sample data according to an embodiment of the present invention.
Fig. 5 is a block diagram of a track characteristic modeling apparatus for track search according to an embodiment of the present invention.
Fig. 6 is a block diagram of an electronic device according to an embodiment of the invention.
Detailed Description
The following detailed description of the embodiments of the present invention will be provided with reference to the drawings and examples, so that how to apply the technical means to solve the technical problems and achieve the technical effects can be fully understood and implemented. It should be noted that the detailed description is only for the purpose of making the invention easier and clearer for those skilled in the art, and is not intended to be a limiting explanation of the invention; moreover, as long as there is no conflict, the embodiments and the features of the embodiments of the present invention may be combined with each other, and the technical solutions formed are all within the scope of the present invention.
Additionally, the steps illustrated in the flow charts of the drawings may be performed in a control system such as a set of controller-executable instructions and, although a logical ordering is illustrated in the flow charts, in some cases, the steps illustrated or described may be performed in an order different than that illustrated herein.
The technical scheme of the invention is described in detail by the figures and the specific embodiments as follows:
examples
In order to solve at least one technical problem of insufficient precision and processing speed which cannot meet the real-time requirement in the prior art, the technical scheme provided by the embodiment maps all positions of the earth to a certain region code through a multi-level earth position coding algorithm, then concatenates codes of aircraft tracks passing through the region, and then finds similar tracks by utilizing a preset algorithm.
The flight path retrieval provided by the embodiment includes, but is not limited to, flight path retrieval corresponding to civil aviation, unmanned aerial vehicles and/or reconnaissance aircraft; in addition, in the embodiment, the location where the flight path of the aircraft passes may also be determined by using different signal acquisition modes, for example, signal trackers of GPS, beidou, ADS-B, and the like, in combination with the specific type of the aircraft, and determining the position of the aircraft according to the signal trackers, and then determining the flight path of the aircraft according to the position of the aircraft and the following mentioned flight path feature modeling method. And further, the historical flight path of the aircraft or the allowed flight path of each aircraft can be combined to be set as a flight path area chain, and when the flight path of the new aircraft is monitored, whether the new aircraft corresponds to the flight path or the allowed flight path of each aircraft can be quickly determined; the flight path of the aircraft in a preset airspace is effectively monitored; and the aircraft which does not meet the requirements can be focused, and if the condition of influencing the normal flight of other aircraft or the condition of the aircraft in a flight forbidden zone occurs, the early warning is carried out.
As shown in fig. 1, the present embodiment provides a track feature modeling method for track retrieval, where the track feature modeling method includes:
s110, establishing a coordinate system of the earth surface, and determining different earth region division levels according to different accuracy of the divided regions;
s120, respectively establishing earth position dictionary code libraries with different levels of precision according to different division region precision; the earth position dictionary code library with each grade precision sets a number for each region in the earth region respectively aiming at the corresponding region precision;
s130, connecting lines from a flight starting point to a flight terminal point and from points of an approach area one by one according to a flight path of the aircraft to form a path area chain; and the track area chain comprises: forming corresponding flight tracks by combining a flight starting point, point-by-point of a path region and a flight terminal point in earth position dictionary code libraries based on different grades respectively;
when the flight path of the aircraft needs to be inquired, whether the corresponding flight path exists is respectively searched in the path region chain from the earth position dictionary code library with low precision to the earth position dictionary code library with high precision.
The earth position coding dictionary mentioned in this embodiment is to divide all positions on the earth into a certain area by a reasonable area division means, and the division depends on factors such as population density, data in a track library, airport coordinates and the like.
Taking 4-level coding as an example, requirements can be gathered, and coding modes of other levels can be set. The surface area of the earth is calculated by 5.1 hundred million square kilometers, and 4-level coding is adopted in the embodiment during verification:
(1) first-order coding is intended to divide the earth into no more than 10000 regions (the average region size is about 225km squares on a side).
(2) The two-level encoding is intended to divide the earth into no more than 100000 regions (with an average region size of approximately 71km square on a side).
(3) Three level coding is intended to divide the earth into no more than 1000000 regions (with an average region size of about 23km square on a side).
(4) Four-level coding is intended to divide the earth into no more than 10000000 regions (with an average region size of about 7km squares on a side).
The track characteristic mentioned in this embodiment is that a track is expressed by a set of measurable character strings, numerical values, etc. The preferred track characteristic design satisfies the following two points:
a: if the two tracks are similar, then the features must be similar (feature distance is close to zero);
b: if similar tracks exist in the track library, N tracks are obtained through track characteristic comparison, and the similar tracks are in the N tracks; the more the similar tracks are ranked in the result, the more perfect the track characteristic design is.
In a preferred embodiment of this embodiment, a flight history track area chain is formed based on the flight history of the aircraft; when an aircraft appears in an area to be monitored, according to the flight path of the aircraft, in a flight history path area chain, whether the same or similar flight paths exist in the flight history path area chain is searched through a method of gradually thinning the earth position dictionary code library level.
In a preferred embodiment of this embodiment, determining different earth region division levels according to different precisions of the divided regions specifically includes:
determining an initialization area parameter of the earth area based on longitudes and latitudes with different sizes;
and merging the initialized areas according to the maximum weighted population number and the maximum area of the areas.
In a further preferred implementation manner of the embodiment of the present invention, determining the initialized area parameter of the earth area includes:
triple for each area
To indicate that the user is not in a normal position,
wherein { (j)1,w1),(j2,w2) Represents the location of the area;
and { (j)1,w1),(j2,w2) The area calculation of the corresponding area comprises:
the average latitude of the area is calculated as:
calculate the area aspect ratio of the area:
calculate the regionActual area of the region of (1): s ═ w1-w2|*|j1-j2|*cos((w1+w2)/2)。
In a further preferred implementation manner of the embodiment of the present invention, the merging the initialized areas according to the maximum weighted population number of the areas and the maximum area of the areas specifically includes:
determining world population density according to the population number on a unit land area of a global land surface, and determining the population density corresponding to each region in a global position dictionary code library with different levels of precision;
in the track area chain, the population density of the area where the flight starting point and the flight ending point of each track are positioned is increased by a first preset value, and the population density of each area is increased by a second preset value in the areas where the rest tracks pass through;
calculating an area weighted area (area actual area × (1+ α × log (1+ population density increase predetermined value/area population));
and combining the initialized regions of the earth regions by using a square principle and a gravity center centering algorithm to obtain the combined earth regions.
In a preferred embodiment of this embodiment, the method for establishing the global position dictionary code library includes: setting a number for each region, and then establishing an index by applying a general red-black tree algorithm to form a global position dictionary; any position of the earth region can be directly inquired about the region code through dictionary indexing, and the region information can be directly inquired about through the dictionary according to the region code.
In this embodiment, a specific implementation manner of the global position dictionary coding, as shown in fig. 2, includes:
s210, determining parameters according to the position dictionary in a grading manner:
before the work begins, a coordinate system of the earth surface needs to be established, and since the most famous global satellite positioning systems (GLONASS, GPS, Beidou and the like) all adopt longitude and latitude coordinates, the longitude and latitude coordinates are also adopted in the invention, the longitude range is-180 degrees to-180 degrees, and the latitude range is-90 degrees to-90 degrees. However, since the earth is a sphere, if the length-width ratio and area formula under the classical coordinate system are used, the obtained result is wrong, and therefore, the formula must be corrected in a targeted manner, and the specific correction formula is shown in section 2.3.
In this coordinate system, the first step must determine the parameters, mainly the earth region initialization region parameters and the termination conditions for region merging.
(1) Initializing regional parameters for an earth region
In the case of one-level coding, we divide one region every 0.3 °, in which case 1200 × 600 — 72 ten thousand regions are initialized.
In the case of two-level coding, we divide one region every 0.1 °, in which case we initialize 3600x1800 ═ 648 ten thousand regions.
In the case of three-level coding, we divide one region every 0.03 °, in which case 12000 × 6000 — 0.72 billion regions are initialized.
In the case of four-level coding, we divide one region every 0.01 °, in which case we initialize 36000x18000 — 6.48 billion regions.
(2) Termination condition of region merging (including region maximum weighted population, region maximum area limit)
In the case of first-level coding, the actual area of the maximum region is 25 kilo-square kilometers, and the weighted maximum population is not more than 1 hundred million.
Under the condition of two-level coding, the actual area of the maximum area is 2.5 ten thousand square kilometers, and the weighted maximum population is not more than 1000 ten thousand.
Under the condition of three-level coding, the actual area of the maximum area is 2500 square kilometers, and the weighted maximum population is not more than 100 ten thousand.
In the case of four-level coding, the maximum area actual area is 250 square kilometers, and the weighted maximum population is not more than 10.
S220, initializing earth region division:
for each region
Such a triple is to be represented by,wherein { (j)1,w1),(j2,w2) Denotes the position of the area,
representing the weighted population of the region.
1. Calculation of actual area of region
The method for calculating the area of the regions can intuitively feel that each region near the equator is larger and the area of the regions near the north and south poles is smaller because the earth is spherical. Let the coordinates of a region be { (j)1,w1),(j2,w2) }, then
Average latitude
Aspect ratio of regions
It can be seen that near the pole,
approximately 90 °, aspect ratio approximately 0; in the vicinity near the equator, it is preferable that,
near 0 deg., the warp and weft weighting is near 1, that is, the initialized region near the equator is closest to the square.
Area actual area s ═ w1-w2|*|j1-j2|*cos((w1+w2)/2)
2. Calculation of weighted population
Two factors are mainly considered for calculating the weighted population number of the region, namely the real population number, the real population number is calculated by multiplying the world population density by the area to simulate calculation, and 78707 main city databases in the world are downloaded, wherein the databases comprise longitude and latitude values and population data; another application in conjunction with the present invention fuses historical track data into population size.
Population density of the worldThe degree is the number of people per unit area of land on the global land surface. The usual unit of use is human/km2Or human/hectare. The index of the population density of the world can be used for simply and clearly reflecting the difference of the population distribution among regions. For example, the world average population density was 39 people/km in 19902With Europe exceeding 10.1 persons/km2In oceania only 3 persons/km2. The world population density is a significant advantage of low density areas from the entire earth. One England geography physician measures, without counting two poles ice-covered land, the population density is less than 2 persons/km2The range of (c) accounts for 60% of the land area of the world. From each continent, the high density of the eurasian population has a significantly higher specific gravity than other continents, while the dense regions of the oceania population have a particularly low specific gravity and the european population is most evenly distributed.
By combining the technical scheme provided by the embodiment, the historical track data is fused into the population quantity: for each track, the area where the starting track and the ending track are located, the population density is increased by 5 persons/km2(ii) a The population density of each area is increased by 1/km in areas crossed by other tracks2。
The starting point and the ending point are coordinates of an airport or the like with a high probability, but may be a mobile airport (such as an aircraft carrier) in some cases.
3. Area weighted area calculation:
the area weighted area comprehensively considers the area actual area and the area weighted population number.
For example: area-weighted area (area actual area) (1+ α log (1+ area-weighted population/10000))
In this embodiment, α is 10, that is:
if the population of the area is 1000, the weighted area of the area is increased by about 40 percent;
if the population of the area is 1 ten thousand, the area weighted area is increased by about 3 times;
if the population of the area is 10 thousands of people, the area weighted area is increased by about 10 times;
if the number of the regional population is 100 ten thousand, the weighted area of the region is increased by 20 times;
other numbers are also calculated by the above formula.
S230, combining position areas:
location area consolidation includes but is not limited to the following 3 factors:
1. regarding which regions need to be merged, in the technical solution provided in this embodiment, the merging region selection first solves the first problem, and the regions most needing to be merged have the following characteristics:
A. the smaller the area weighted is, the more the region weighted needs to be merged;
B. the larger the aspect ratio, the more need to be incorporated;
C. the more surrounding mergeable regions, the more merging is required.
With the above 3 items, we can easily conclude that the regions near the north-south pole are merged first, because these regions almost satisfy all the above conditions:
the area is minimum, the population is small, and the aircraft passes through the aircraft, so the weighted population is small, and the final weighted area is also minimum;
the aspect ratio is the largest of all initialization regions;
the nearby areas are similar, and the number of combinable areas is large;
in the technical scheme provided by this embodiment, the 3 conditional quantization indexes of each initialization area are calculated, then pseudo-randomization of the indexes is performed (i.e., random adjustment of ± 1% is performed on the basis of the indexes obtained by final calculation), then the indexes are sorted, an area queue to be merged is generated, and the areas sorted in the front are merged first.
2. With respect to how to merge:
if one region can merge other regions, but if there are multiple merging schemes, how to choose the merging scheme, as shown in fig. 3, the main region can merge the four surrounding regions, which is how to choose?
In the technical solution provided in this embodiment, there are two merging principles:
(1) and a square principle: calculating the aspect ratio of each scheme after combination, wherein the closer to 1, the better;
(2) the center of gravity centering principle: the combined center of gravity (calculated from the weighted population density) for each solution was calculated, with the center of gravity being as centered as possible.
3. Under what circumstances the merging is stopped. Wherein, the condition for terminating the merging is as follows: the termination conditions have been actually determined in the above embodiments; and if the actual area of the area or the population of the area is larger than the threshold, the merging is terminated.
S240, area division evaluation:
the rule of the partition evaluation is simple, that is, the smaller the number of the partitioned areas is, the better the partition evaluation is, and the smaller the standard deviation of the area weighted area is, the better the partition evaluation is.
Wherein the area evaluation value is weighted area mean-alpha weighted area standard deviation
In this example, α is 1.
Due to the fact that the pseudo-random strategy is adopted when the area queues to be merged are generated, the randomness of the merging process in each time can be guaranteed. And after random combination is performed for a plurality of times, selecting the region with the highest score as the basis of the dictionary code of the level.
S250, judging that the division meets the requirement? Determining whether a predetermined requirement is satisfied based on a result of the area division evaluation in S240, and if so, performing S260; otherwise, returning to S220, the process continues to perform the earth region division initialization process for the region.
S260, determining the global position dictionary code: through the embodiment, the earth position dictionary codes corresponding to different precision levels can be obtained; and combining the earth position dictionary codes with different precision levels into an earth position dictionary code library. In each earth position dictionary code, each region is given a number, and then a general red-black tree algorithm is applied to establish an index to form an earth position dictionary; thus, for any position, the codes of the regions can be directly inquired through dictionary indexing; similarly, through the region coding, the region information can also be directly inquired through the dictionary.
In a preferred embodiment of this embodiment, the track characteristics include a track area chain, a flight time period, a flight duration, and a flight height, and the track block chain includes a flight starting point, a flight ending point, an earth area through which points pass one by one, and an area through which a connection line of two adjacent points passes; the flight time period is used for recording the starting time and the ending time, the flight time length is used for recording the flight time length, and the flight height is used for recording the maximum height, the maximum height and the average height; the flight speed is used for recording the maximum speed, the fastest speed and the average speed; the flight path characteristics are used for executing the same or similar flight path searching, and the flight path retrieval comprises the steps of carrying out preliminary screening by using the area chain with low precision, and then determining similar flight paths by using the area chain with high precision and combining the flight path area chain, the flight time period, the flight time length and the flight height.
Specifically, as shown in fig. 4, is a portion of sample data for a flight path, where the first two columns are latitude and longitude coordinates, the third column is height, the fourth column is speed, and the fifth column is data acquisition time.
The track characteristics provided by the embodiment mainly consider factors such as a track area chain, a flight time period, flight duration, flight height and the like, and flight speed, wherein the track area chain is an important consideration.
The track area chain is calculated from a starting point to an end point, an area which is passed by points one by one and an area which is connected with adjacent points and passed by a connecting line. The characteristics of the track area chain are as follows: ordered and non-repeating. The order means that the track areas are connected into an ordered set; non-repeating means that if the aircraft was in zone a at the previous moment, no action was done if the aircraft was still in zone a at that moment, and B was added to the flight path zone chain only when the position of the aircraft changed to zone B at a certain moment. The track area chains are respectively generated according to dictionaries of different levels, each dictionary corresponds to one track area chain, and the track area chain generated by the primary position dictionary is called a primary area chain for short; similarly, the track area chain generated by the secondary position dictionary is called the secondary area chain for short.
And a flight time period, wherein the start time and the end time are mainly recorded.
And (4) the flight time mainly records the flight time.
Flying height, mainly recording maximum height, maximum height and average height.
The flying speed mainly records the maximum speed, the fastest speed and the average speed.
Based on the preferred embodiment, the method for realizing the rapid search of similar tracks comprises the following steps:
firstly, a first-level regional chain is utilized to carry out primary screening, then the first-level regional chain is gradually refined, and finally, a fourth-level regional chain is integrated to be compared with a plurality of indexes such as a flight time period, a flight duration, a flight height and a flight speed until a similar track is found.
Because a primary region chain is adopted firstly, and the primary region chain is roughly divided, the performance is better, and the comparison range can be quickly reduced; and the last comparison is carried out by adopting a four-level regional chain to assist other characteristics, so that the precision is high.
The weighted euclidean distance is used as the similar distance, and is a general algorithm, but the embodiment is not limited thereto. Specifically, the method comprises the following steps:
assuming that tracks P and Q exist, calculating a track similarity distance RPQThe process of (2) is as follows:
RPQ＝α1Dlink+α2Dtime+α3Dcost+α4Dhigh+α5Dspeed
RPQthe smaller the correlation between the two tracks.
Where α 1 to α 5 are weights of respective calculation dimensions, and α 1+ α 2+ α 3+ α 4+ α 5 is 1.
DlinkIndicating regional chain relevance;
Dtimerepresenting a time of flight correlation;
Dcostrepresenting a long-term flight correlation;
Dhighindicating a flying height dependency;
Dspeedand represents the flight speed dependency.
In this embodiment, the area chain correlation considerations include:
1. the degree of overlap of the code sets in the region chain,
2. the region chain is the longest common sub-string,
3. edit distance of region chain.
In this embodiment, the time-of-flight correlation considerations include:
1. the start time of flight (time in minutes and seconds),
2. the end-of-flight time (time in minutes and seconds),
3. the flight is calculated periodically (e.g., once a day every week, once a day every month).
In this embodiment, the correlation consideration regarding the flight duration is the difference in flight duration.
In the present embodiment, the consideration regarding the flying height correlation includes:
1. the average value of the flying height is,
2. the standard deviation of the flying height is obtained,
3. flight height variation trend.
In the present embodiment, the considerations regarding the dependence of the flying speed include:
1. the average value of the flying speeds of the aircraft,
2. the standard deviation of the flying speed is obtained,
3. and (4) flight speed variation trend.
Therefore, in the above technical solution provided in this embodiment, all positions of the earth are mapped to a certain region code by a multi-level earth position coding algorithm, and then codes of aircraft flight path passing regions are concatenated, so as to form a flight path region chain convenient for identification. In addition, in the process of track identification, as the earth position code adopts a multi-stage algorithm, similar tracks can be searched by a gradual thinning method, and the retrieval efficiency is greatly improved; for example, a coarsest divided regional dictionary is adopted, a flight path set meeting the conditions is found through preliminary screening, then the refined regional dictionary is used step by step, and finally the most similar flight path is accurately positioned.
As shown in fig. 5, the present embodiment further provides a track characteristic modeling system 100 for track retrieval, where the track characteristic modeling system 100 for track retrieval includes:
the earth region dividing module 110 is configured to create a coordinate system of the earth surface, and determine different earth region division levels according to different precision of the divided regions;
the earth position dictionary code library establishing module 120 is used for respectively establishing earth position dictionary code libraries with different levels of precision according to the difference of the divided region precision; the earth position dictionary code library with each grade precision sets a number for each region in the earth region respectively aiming at the corresponding region precision;
a track area chain forming module 130, configured to form a track area chain by connecting points, one by one, from a flight starting point to a flight ending point and a path area according to a flight track of the aircraft; and the track area chain comprises: forming corresponding flight tracks by combining a flight starting point, point-by-point of a path region and a flight terminal point in earth position dictionary code libraries based on different grades respectively;
when the flight path of the aircraft needs to be inquired, whether the corresponding flight path exists is respectively searched in the path region chain from the earth position dictionary code library with low precision to the earth position dictionary code library with high precision.
With regard to specific implementation manners of the modules in the track feature modeling system for track retrieval 100, reference may be made to fig. 1 to 4 and descriptions in the above-mentioned track feature modeling method for track retrieval, which are not described herein again.
As shown in fig. 6, the present embodiment further provides an electronic device 200, where the electronic device 200 includes:
the memory(s) 210 are (are),
a processor 220, and
a computer program;
wherein the computer program is stored in the memory and configured to be executed by the processor to implement any of the track feature modeling methods for track retrieval as provided above.
In addition, the present embodiment also provides a nonvolatile storage medium on which a computer program is stored; the computer program is executed by a processor to implement any of the above provided track feature modeling methods for track retrieval.
Those of ordinary skill in the art will understand that: the above-described method according to an embodiment of the present invention may be implemented in hardware, firmware, or as software or computer code storable in a recording medium such as a CD ROM, a RAM, a floppy disk, a hard disk, or a magneto-optical disk, or as computer code originally stored in a remote recording medium or a non-transitory machine-readable medium downloaded through a network and to be stored in a local recording medium, so that the method described herein may be stored in such software processing on a recording medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware such as an ASIC, an FPGA, or an SoC. It will be appreciated that the computer, processor, microprocessor controller or programmable hardware includes memory components (e.g., RAM, ROM, flash memory, etc.) that can store or receive software or computer code that, when accessed and executed by the computer, processor or hardware, implements the processing methods described herein. Further, when a general-purpose computer accesses code for implementing the processes shown herein, execution of the code transforms the general-purpose computer into a special-purpose computer for performing the processes shown herein.
Those of ordinary skill in the art will appreciate that the various illustrative elements and method steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present embodiments.
Finally, it should be understood that the above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention in any way. Those skilled in the art can make many changes and simple substitutions to the technical solution of the present invention without departing from the technical solution of the present invention, and the technical solution of the present invention is protected by the following claims.
Claims (10)
1. A track feature modeling method for track retrieval, comprising:
establishing a coordinate system of the earth surface, and determining different earth region division levels according to different precision of the division regions;
respectively establishing earth position dictionary code libraries with different levels of precision according to the difference of the divided region precision; the earth position dictionary code library with each grade precision sets a number for each region in the earth region respectively aiming at the corresponding region precision;
according to the flight path of the aircraft, connecting lines from a flight starting point to a flight terminal point and from points to points of a path region one by one to form a path region chain; and the track area chain comprises: forming corresponding flight tracks in earth position dictionary code libraries of different levels based on point-by-point combination with a flight starting point and a route region and a flight terminal point;
when the flight path of the aircraft needs to be inquired, sequentially searching whether corresponding flight paths exist in the path region chain from a low-precision earth position dictionary code library to a high-precision earth position dictionary code library.
2. The method of modeling track characteristics according to claim 1, further comprising: forming a flight history track area chain based on the flight history of the aircraft; when an aircraft appears in an area to be monitored, according to the flight path of the aircraft, in the flight history path area chain, whether the same or similar flight paths exist in the flight history path area chain is searched through a method of gradually thinning the earth position dictionary code library level.
3. The track feature modeling method of claim 1, wherein different earth region division levels are determined according to different accuracy of the division regions, and specifically comprises:
determining an initialization area parameter of the earth area based on longitudes and latitudes with different sizes;
and merging the initialized areas according to the maximum weighted population number and the maximum area of the areas.
4. The method of claim 3, wherein determining an earth region initialization zone parameter comprises:
triple for each area
To indicate that the user is not in a normal position,
wherein { (j)1,w1),(j2,w2) Represents the location of the area;
and { (j)1,w1),(j2,w2) The area calculation of the corresponding area comprises:
the average latitude of the area is calculated as:
calculate the area aspect ratio of the area:
calculating the area actual area of the area: s ═ w1-w2|*|j1-j2|*cos((w1+w2)/2)。
5. The track feature modeling method according to claim 3, wherein the initialization area is merged according to the area maximum weighted population number and the area maximum area, and specifically includes:
determining world population density according to the population number on a unit land area of a global land surface, and determining the population density corresponding to each region in a global position dictionary code library with different levels of precision;
in the track area chain, the population density of the area where the flight starting point and the flight ending point of each track are positioned is increased by a first preset value, and the population density of each area is increased by a second preset value in the areas where the rest tracks pass through;
calculating an area weighted area (area actual area × (1+ α × log (1+ population density increase predetermined value/area population));
and combining the initialized regions of the earth regions by using a square principle and a gravity center centering algorithm to obtain the combined earth regions.
6. The track feature modeling method of claim 1, wherein the global position dictionary code library is established in a manner that includes: setting a number for each region, and then establishing an index by applying a general red-black tree algorithm to form a global position dictionary; any position of the earth region can be directly inquired about the region code through dictionary indexing, and the region information can be directly inquired about through the dictionary according to the region code.
7. The track feature modeling method of any of claims 1-6, wherein the track features include track area chain, flight time period, flight duration, flight altitude; the flight path block chain comprises a flight starting point, a flight terminal point, an earth area through which points pass one by one and an area through which a connecting line of two adjacent points passes; the flight time period is used for recording the starting time and the ending time, the flight time length is used for recording the flight time length, and the flight altitude is used for recording the maximum altitude, the maximum altitude and the average altitude; the flight speed is used for recording the maximum speed, the fastest speed and the average speed; the track features are used for executing the same or similar track search, and the track retrieval comprises the steps of carrying out preliminary screening by using the area chain with low precision, and then determining similar tracks by using the area chain with high precision and combining the track area chain, the flight time period, the flight duration and the flight height.
8. A track feature modeling system for track retrieval, comprising:
the earth region dividing module is used for creating a coordinate system of the earth surface and determining different earth region dividing grades according to different accuracy of divided regions;
the earth position dictionary code library establishing module is used for connecting lines from a flight starting point to a flight destination and from points to points of an approach area one by one according to a flight path of an aircraft to form a path area chain; and the track area chain comprises: forming corresponding flight tracks in earth position dictionary code libraries of different levels based on point-by-point combination with a flight starting point and a route region and a flight terminal point;
the flight path region chain forming module is used for connecting lines from a flight starting point to a flight terminal point and from points to points of the path region one by one according to the flight path of the aircraft to form a flight path region chain; and the track area chain comprises: forming corresponding flight tracks by combining a flight starting point, point-by-point of a path region and a flight terminal point in earth position dictionary code libraries based on different grades respectively;
when the flight path of the aircraft needs to be inquired, sequentially searching whether corresponding flight paths exist in the path region chain from a low-precision earth position dictionary code library to a high-precision earth position dictionary code library.
9. An electronic device, comprising:
a memory for storing a plurality of data to be transmitted,
a processor, and
a computer program;
wherein the computer program is stored in the memory and configured to be executed by the processor to implement the method of modeling track features for track retrieval as claimed in any one of claims 1 to 7.
10. A non-volatile storage medium having a computer program stored thereon; the computer program is executed by a processor to implement a method of modeling track features for track retrieval as claimed in any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010975489.XA CN112214537A (en) | 2020-09-16 | 2020-09-16 | Track characteristic modeling method and system for track retrieval and electronic device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010975489.XA CN112214537A (en) | 2020-09-16 | 2020-09-16 | Track characteristic modeling method and system for track retrieval and electronic device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112214537A true CN112214537A (en) | 2021-01-12 |
Family
ID=74049514
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010975489.XA Pending CN112214537A (en) | 2020-09-16 | 2020-09-16 | Track characteristic modeling method and system for track retrieval and electronic device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112214537A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113177037A (en) * | 2021-04-25 | 2021-07-27 | 上海鹰觉科技有限公司 | Flight path storage and retrieval method, system, equipment and computer storage readable medium |
| CN114812568A (en) * | 2022-03-01 | 2022-07-29 | 航天科工智能运筹与信息安全研究院(武汉)有限公司 | Route matching area selection method based on rule scoring |
| CN114898600A (en) * | 2022-07-13 | 2022-08-12 | 安胜（天津）飞行模拟系统有限公司 | UAM four-dimensional track sharing and management method based on block chain technology architecture |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106446194A (en) * | 2016-09-29 | 2017-02-22 | 北京百度网讯科技有限公司 | Determination method and device for acquisition parameters of a geographic region |
| CN106980129A (en) * | 2017-03-10 | 2017-07-25 | 西安交通大学 | A kind of movement locus comparison method based on position encoded map |
| CN110427360A (en) * | 2019-06-28 | 2019-11-08 | 东软集团股份有限公司 | Processing method, processing unit, processing system and the computer program product of track data |
| CN111141294A (en) * | 2019-12-27 | 2020-05-12 | 西安交通大学 | Position coding based motion trajectory comparison method |
| CN111292356A (en) * | 2020-02-20 | 2020-06-16 | 清华大学 | Method and device for matching motion trail with road |
| CN111522892A (en) * | 2019-02-02 | 2020-08-11 | 阿里巴巴集团控股有限公司 | Geographic element retrieval method and device |
-
2020
- 2020-09-16 CN CN202010975489.XA patent/CN112214537A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106446194A (en) * | 2016-09-29 | 2017-02-22 | 北京百度网讯科技有限公司 | Determination method and device for acquisition parameters of a geographic region |
| CN106980129A (en) * | 2017-03-10 | 2017-07-25 | 西安交通大学 | A kind of movement locus comparison method based on position encoded map |
| CN111522892A (en) * | 2019-02-02 | 2020-08-11 | 阿里巴巴集团控股有限公司 | Geographic element retrieval method and device |
| CN110427360A (en) * | 2019-06-28 | 2019-11-08 | 东软集团股份有限公司 | Processing method, processing unit, processing system and the computer program product of track data |
| CN111141294A (en) * | 2019-12-27 | 2020-05-12 | 西安交通大学 | Position coding based motion trajectory comparison method |
| CN111292356A (en) * | 2020-02-20 | 2020-06-16 | 清华大学 | Method and device for matching motion trail with road |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113177037A (en) * | 2021-04-25 | 2021-07-27 | 上海鹰觉科技有限公司 | Flight path storage and retrieval method, system, equipment and computer storage readable medium |
| CN114812568A (en) * | 2022-03-01 | 2022-07-29 | 航天科工智能运筹与信息安全研究院(武汉)有限公司 | Route matching area selection method based on rule scoring |
| CN114812568B (en) * | 2022-03-01 | 2024-05-28 | 航天科工智能运筹与信息安全研究院(武汉)有限公司 | Track matching area selection method based on rule grading |
| CN114898600A (en) * | 2022-07-13 | 2022-08-12 | 安胜（天津）飞行模拟系统有限公司 | UAM four-dimensional track sharing and management method based on block chain technology architecture |
| CN114898600B (en) * | 2022-07-13 | 2022-10-21 | 安胜（天津）飞行模拟系统有限公司 | UAM four-dimensional track sharing and management method based on block chain technology architecture |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112214537A (en) | Track characteristic modeling method and system for track retrieval and electronic device | |
| KR102232138B1 (en) | Neural architecture search method based on knowledge distillation | |
| US11562168B2 (en) | Clustering for K-anonymity in location trajectory data | |
| CN110726418B (en) | Method, device and equipment for determining interest point region and storage medium | |
| CN110954113B (en) | Vehicle pose correction method and device | |
| CN106871932B (en) | Satellite-borne laser on-orbit pointing calibration method based on pyramid search terrain matching | |
| US10546043B1 (en) | Triangulation for K-anonymity in location trajectory data | |
| US6317690B1 (en) | Path planning, terrain avoidance and situation awareness system for general aviation | |
| CN108225341B (en) | Vehicle positioning method | |
| CN108594274B (en) | Combined navigation satellite selection method for adaptive chaotic particle swarm optimization | |
| CN109492580B (en) | Multi-size aerial image positioning method based on neighborhood significance reference of full convolution network | |
| CN112069276B (en) | Address coding method, address coding device, computer equipment and computer readable storage medium | |
| CN115077556B (en) | Unmanned vehicle field operation path planning method based on multi-dimensional map | |
| CN114065364B (en) | Urban engineering planning method and system based on unmanned aerial vehicle remote sensing mapping | |
| CN111291776A (en) | Channel information extraction method based on crowd-sourced trajectory data | |
| CN111256696B (en) | Aircraft autonomous navigation method with multi-feature and multi-level scene matching | |
| US11465782B2 (en) | Systems and methods for autonomous deorbiting of a spacecraft | |
| Jiang et al. | Redundant-coded radial and neighbor star pattern identification algorithm | |
| CN113302897A (en) | User information management method, identity recognition method and device | |
| CN114563003A (en) | Path planning method and device, electronic equipment and storage medium | |
| FR2993973A1 (en) | METHOD OF PROCESSING A FLIGHT PLAN | |
| CN116026352A (en) | Processing method and device of positioning signal, medium and electronic equipment | |
| CN114935773A (en) | Online road matching method, device, medium and equipment based on Beidou positioning | |
| CN112434228B (en) | Method for predicting track position of moving target | |
| CN114664104A (en) | Road network matching method and device |
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 |