CN115810087A - Low-altitude space domain use conflict detection method based on multi-scale space grid - Google Patents

Low-altitude space domain use conflict detection method based on multi-scale space grid Download PDF

Info

Publication number
CN115810087A
CN115810087A CN202310047397.9A CN202310047397A CN115810087A CN 115810087 A CN115810087 A CN 115810087A CN 202310047397 A CN202310047397 A CN 202310047397A CN 115810087 A CN115810087 A CN 115810087A
Authority
CN
China
Prior art keywords
grid
airspace
aircraft
low
altitude
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
Application number
CN202310047397.9A
Other languages
Chinese (zh)
Inventor
张昕
石畅
刘程威
刘玉
杜宏
陶德进
闫嘉明
王刚
胡国文
牛帅
沈驰
刘锡明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CETC 28 Research Institute
Original Assignee
CETC 28 Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CETC 28 Research Institute filed Critical CETC 28 Research Institute
Priority to CN202310047397.9A priority Critical patent/CN115810087A/en
Publication of CN115810087A publication Critical patent/CN115810087A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Traffic Control Systems (AREA)

Abstract

The invention discloses a low-altitude airspace use conflict detection method based on a multi-scale spatial grid, which comprises the following steps of carrying out rasterization modeling and digital coding, and establishing the built-in attribute of the low-altitude airspace grid; aiming at different aircraft physical sizes and safety distances in a given low-altitude airspace, the space-time occupation characteristics of the airspace are digitally described by three-dimensional matrixes corresponding to different grid scales, and multi-scale rasterization representation of a traditional airspace use area is realized; the method comprises the steps of constructing a low-altitude airspace use conflict detection function based on a multi-scale airspace use three-dimensional grid matrix, fully utilizing space-time information expressed by different scales using airspace three-dimensional matrix codes, designing a hierarchical airspace conflict detection flow, and realizing conflict detection of airspaces used by aircrafts of different sizes. The method can simultaneously consider the spatial-temporal information of the airspace occupied by the low-altitude aircraft with different sizes, effectively detect the position and time of the low-altitude airspace use conflict, improve the airspace use conflict detection precision and reduce the false judgment rate.

Description

Low-altitude space domain use conflict detection method based on multi-scale space grid
Technical Field
The invention relates to the field of airspace management and control, in particular to a low-altitude airspace use conflict detection method based on a multi-scale space grid.
Background
The airspace conflict detection is a key component for guaranteeing the safe operation of the airspace in China and orderly carrying out the airspace use plan, and how to quickly and accurately judge the conflict airspace of the airspace use plan is a key problem of the future airspace collaborative planning.
The research on the spatial domain conflict detection starts in the last 40-50 years, various correlation models and algorithms have been proposed by many scholars at home and abroad, geometric floating point calculation is most widely used at present, namely, whether spatial domain conflicts exist is judged by intersecting the edges of the spatial domain required by each spatial domain use plan, although the method can accurately calculate the spatial domain use plan conflicts and the range of the conflict spatial domain, the method has the problems of long calculation time, low efficiency and the like for large-scale spatial domain conflict detection.
There is therefore a need to solve the above problems.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide a low-altitude spatial domain use conflict detection method based on a multi-scale spatial grid, which is efficient and high in precision.
The technical scheme is as follows: in order to achieve the above purpose, the invention discloses a low-altitude space domain use conflict detection method based on a multi-scale space grid, which comprises the following steps:
(1) Performing airspace rasterization modeling and digital coding, namely performing discretization processing on a given continuous low-altitude airspace, dividing the continuous airspace into gapless cube grid models with different scales, wherein a cube grid of each scale can completely describe all the given low-altitude airspaces, and performing three-dimensional matrix type sequential coding on the cube grids of each scale to form a multi-scale low-altitude airspace grid coding system with codes mapped with the cube grids of different scales one by one, so as to form a low-altitude airspace discretization description mode;
(2) The method comprises the steps that a low-altitude airspace is represented by multi-scale rasterization, the size of an aircraft in the low-altitude airspace and a corresponding safety distance of the aircraft are considered, a low-altitude airspace route of the aircraft is represented as cubic grid matrix coded data, grid index coordinates and corresponding grid matrix data respectively describe space and time information of an airspace occupied by the aircraft, and space-time information of the low-altitude airspace occupied by the aircraft is described by three-dimensional grid matrix elements with different scales;
(3) Constructing an airspace use conflict detection function, firstly judging the conflict situation of different aircraft use airspaces on the space based on the three-dimensional grid matrix description of the aircraft use airspace, and further considering the time conflict situation of the aircraft occupied in the airspace aiming at the airspace grid with conflict, thereby determining the space-time conflict range of the different aircraft use airspaces;
(4) Designing an airspace conflict detection flow, starting from a large-scale three-dimensional grid matrix on the basis of a multi-scale airspace use grid matrix, detecting the conflict situation of the different types of aircraft airspaces under the scale, and further detecting the conflict situation of the aircraft use airspace under a small-scale three-dimensional grid in a local conflict area aiming at the aircraft with conflicts.
Wherein, the given continuous low-altitude airspace in the step (1) is the altitude of
Figure SMS_1
A flight zone below, the size being above ground
Figure SMS_2
Left and right continuous low-altitude airspace.
Preferably, the spatial domain rasterization coding of different scales in the step (1) specifically comprises the following steps:
(1.1) constructing a multi-scale space grid model of a given low-altitude airspace, dividing the given low-altitude airspace into cubic grids with different sizes and side lengths, realizing multi-scale gapless cubic grid modeling of a continuous low-altitude airspace, dividing the airspace into three scales, and respectively setting the sizes of the corresponding cubic grids to be three
Figure SMS_5
Figure SMS_7
And
Figure SMS_9
Figure SMS_4
the cubic grid of dimensions may include
Figure SMS_8
An
Figure SMS_10
Cubic grid of dimensions, likewise
Figure SMS_11
The size of the cubic grid can be contained in
Figure SMS_3
An
Figure SMS_6
A cubic grid of dimensions; the original point of Cartesian coordinates is selected as the vertex of the low-altitude airspace, and the cubic grid of the low-altitude airspace divided by three scales can be indexed through three-dimensional coordinates, so that the flight path coordinates of the aircraft represented by a geodetic coordinate system are converted into Cartesian coordinatesLow-altitude airspace coordinates represented by coordinates;
(1.2) coding the low-altitude space domain cube grids divided by different scales in a three-dimensional matrix form, wherein for a given low-altitude space domain, the three-dimensional matrix dimensions corresponding to the low-altitude space domain grids divided by the three scales are respectively
Figure SMS_12
Figure SMS_13
Figure SMS_14
Coding low-altitude space domain grids divided in different scales according to a three-dimensional matrix element index mode, wherein cubic grids in the low-altitude space domain three-dimensional grid division corresponding to each scale can be obtained by corresponding three-dimensional matrix element coordinates
Figure SMS_15
And uniquely determining, and realizing matrix type low-altitude space domain discretization description.
Furthermore, the specific steps of using multi-scale rasterization characterization for the low-altitude airspace in the step (2) are as follows:
(2.1) when the rasterized representation aircraft occupies airspace, the physical sizes of different aircraft need to be considered, and meanwhile, the safety interval of the aircraft during normal flight needs to be considered, so that various aircraft are combined
Figure SMS_19
Radius considered to take into account its physical dimensions and safety interval
Figure SMS_22
The waypoint coordinate point of the aircraft is the center of the circumscribed sphere, and the waypoint represented in the geodetic coordinate system is converted into the waypoint space coordinate represented in the Cartesian coordinate system corresponding to the space grid through coordinate system conversion
Figure SMS_25
And when the aircraft route point is represented in a rasterization manner, the distance between the grid edge and the aircraft route point is set
Figure SMS_18
Is less than
Figure SMS_21
The cubic grids are all regarded as space grids occupied by the air craft routes, wherein
Figure SMS_23
Representing coordinates as
Figure SMS_26
The grid of the space cube of (a),
Figure SMS_16
representing spatial coordinates
Figure SMS_20
The waypoints represented, for a scale of
Figure SMS_24
Figure SMS_27
) Low-altitude airspace grid of (1), grid occupied by course of aircraft
Figure SMS_17
The following inequalities need to be satisfied:
Figure SMS_28
in the formula (I), the compound is shown in the specification,
Figure SMS_30
represent
Figure SMS_34
Coordinate at scale of
Figure SMS_37
The grid of the space cube of (a),
Figure SMS_31
to represent
Figure SMS_33
Cubic grid size at scale, if in course and occupied grid
Figure SMS_35
The space coordinate of the nearest center waypoint is
Figure SMS_38
Then grid
Figure SMS_29
The corresponding time is the aircraft
Figure SMS_32
The waypoint is located at
Figure SMS_36
The time of (d);
(2.2) comparing the flight path data with low-altitude space grid matrixes of different scales based on a low-altitude space aircraft flight path data sequence, namely representing the flight path of the aircraft by using element coordinates and element numerical value sequences of a cubic grid matrix, representing the space grid occupied by the flight path of the aircraft by using the element coordinates of the grid matrix, representing the time of a flight path point of the aircraft to the space grid by using the element numerical value, and aiming at the problem that the flight path data sequence is used by using the space grid matrix
Figure SMS_39
Scale-described low-altitude airspace aircraft represented by cubic grid matrix
Figure SMS_40
Is represented as:
Figure SMS_41
Figure SMS_42
Figure SMS_43
,……
in the formula (I), the compound is shown in the specification,
Figure SMS_46
representing aircraft
Figure SMS_49
In that
Figure SMS_52
The spatial grid coordinate at scale is
Figure SMS_45
Is determined to be the first waypoint of (c),
Figure SMS_48
representing the time corresponding to the waypoint;
Figure SMS_51
representing aircraft
Figure SMS_53
In that
Figure SMS_44
The spatial grid coordinate at scale is
Figure SMS_47
Of the second waypoint of (a) the first waypoint,
Figure SMS_50
representing the time corresponding to the waypoint; by parity of reasoning, the aircraft can be driven
Figure SMS_54
All represented in spatial grid coordinates and their corresponding times.
Further, the collision detection function in the step (3) is specifically constructed by the following steps:
(3.1) for different aircraft
Figure SMS_55
Detecting whether the flight path sequences of different aircrafts contain the same cubic grid coordinate, namely, a certain element coordinate represented by a matrix exists in the flight path sequences of the two aircrafts at the same time
Figure SMS_56
(3.2) the same grid matrix element coordinates exist in different aircraft
Figure SMS_58
Then, whether the time difference value corresponding to the element coordinate of different aircraft route sequences meets the safety time interval or not is further judged
Figure SMS_64
Figure SMS_67
Determining two aircraft routes if the set parameters related to the space grid size are met
Figure SMS_59
And
Figure SMS_61
there is no conflict in the grid space between the grid spaces,
Figure SMS_63
and
Figure SMS_66
are respectively represented by
Figure SMS_57
And
Figure SMS_60
for two marked aircraft routes, collision detection function
Figure SMS_62
Outputting 0, otherwise indicating that two aircraft routes conflict in the grid space, and outputting 1 as a function, wherein the function is expressed as follows:
Figure SMS_65
in the formula (I), the compound is shown in the specification,
Figure SMS_68
is shown in
Figure SMS_69
The grid coordinates in the marked aircraft route are
Figure SMS_70
The time corresponding to the point of the waypoint of,
Figure SMS_71
is shown in
Figure SMS_72
The grid coordinates in the marked aircraft route are
Figure SMS_73
The waypoint of (a) corresponds to the time.
Preferably, the empty domain conflict detection process in the step (4) specifically comprises the following steps:
(4.1) first from the scale
Figure SMS_74
Starting the expressed aircraft route sequence, judging whether the routes of different aircrafts have conflicts under the scale, and if the routes of different aircrafts do not have conflicts, indicating that the corresponding aircrafts do not have conflicts when being used in the current airspace;
(4.2) for on-scale
Figure SMS_75
The space grid with conflict exists in the represented route sequence, and further judgment is carried out according to the scale
Figure SMS_76
Judging whether the flight paths of different aircrafts have conflicts under the scale according to the flight paths of the represented aircrafts in the conflict grid, and if no conflict exists, indicating that no conflict exists when the corresponding aircraft is used in the current airspace;
(4.3) for on-scale
Figure SMS_77
The space grid with conflict exists in the represented route sequence, and further judgment is carried out according to the scale
Figure SMS_78
And judging whether the routes of different aircrafts in the scale conflict with each other according to the routes of the represented aircrafts in the conflict grid, if so, indicating that the corresponding aircraft does not conflict with the current airspace, otherwise, indicating that the corresponding aircraft conflicts with the current airspace.
Has the beneficial effects that: compared with the prior art, the invention has the following remarkable advantages: the method can simultaneously consider the space-time information of the airspace occupied by the low-altitude aircrafts with different sizes, effectively detect the position and time of the low-altitude airspace usage conflict, improve the detection precision of the airspace usage conflict and reduce the misjudgment rate.
Drawings
FIG. 1 is a schematic diagram of a medium and low altitude airspace multi-scale cube grid modeling according to the present invention;
FIG. 2 is a schematic diagram of matrix three-dimensional space grid coding according to the present invention;
FIG. 3 is a first two-dimensional plan view of a first encoded representation of an aircraft flight path according to the present invention;
FIG. 4 is a two-dimensional plane schematic view of a flight path code representation of an aircraft according to the present invention;
FIG. 5 is a schematic two-dimensional plane view of spatial domain usage collision detection at different scales according to the present invention;
FIG. 6 is a flow chart of the low-altitude spatial domain use collision detection based on multi-scale spatial grid modeling according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be further described in detail with reference to the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The invention relates to a low-altitude space domain use conflict detection method based on a multi-scale space grid, which comprises the following steps of:
(1) The spatial domain rasterization modeling and the digital coding are that the given continuous low-altitude spatial domain is discretized, and the continuous spatial domain is divided into gapless cubes with different scalesThe method comprises the steps that a volume grid model is adopted, cubic grids of all scales can completely describe all given low-altitude airspaces, three-dimensional matrix type sequential coding is carried out on the cubic grids of all scales, a multi-scale low-altitude airspace grid coding system with codes mapped with cubic grids of different scales one by one is formed, and a low-altitude airspace discretization description mode is formed; wherein a given continuous low altitude airspace is an altitude of
Figure SMS_79
A flight zone below, the size being above ground
Figure SMS_80
A left and right continuous low-altitude airspace;
the spatial domain rasterization coding of different scales comprises the following specific steps:
(1.1) constructing a multi-scale space grid model of a given low-altitude airspace, as shown in FIG. 1, for the ground size
Figure SMS_83
The method comprises the steps of dividing a given low-altitude airspace into cubic grids with different sizes and side lengths, realizing multi-scale zero-gap cubic grid modeling of a continuous low-altitude airspace, dividing the airspace into three scales, wherein the sizes of the corresponding cubic grids are respectively
Figure SMS_86
Figure SMS_89
And
Figure SMS_82
Figure SMS_85
the size of the cubic grid can be contained in
Figure SMS_88
An
Figure SMS_90
Cubic grid of dimensions, likewise
Figure SMS_81
The cubic grid of dimensions may include
Figure SMS_84
An
Figure SMS_87
A cubic grid of dimensions; selecting a Cartesian coordinate origin as a vertex of a low-altitude airspace, and indexing a low-altitude airspace cube grid divided by three scales through three-dimensional coordinates as shown in FIG. 1, so that the flight path coordinates of an aircraft represented by a geodetic coordinate system can be converted into low-altitude airspace coordinates represented by Cartesian coordinates;
(1.2) coding the cubic grids of the low-altitude airspace divided by different scales in a three-dimensional matrix form, wherein the specific coding mode is as shown in figure 2, and aiming at a given low-altitude airspace, the cubic grids of the low-altitude airspace are coded in order to
Figure SMS_91
Figure SMS_92
Figure SMS_93
The three-dimensional matrix dimensions corresponding to the three-scale low-altitude spatial grid are respectively
Figure SMS_94
Figure SMS_95
Figure SMS_96
Coding low-altitude space domain grids divided in different scales according to a three-dimensional matrix element index mode, wherein cubic grids in the low-altitude space domain three-dimensional grid division corresponding to each scale can be obtained by corresponding three-dimensional matrix element coordinates
Figure SMS_97
The method comprises the steps of determining uniquely, and realizing matrix type low-altitude space domain discretization description;
(2) The method comprises the steps that a low-altitude airspace represents by using multi-scale rasterization, the size of an aircraft in the low-altitude airspace and a corresponding safety distance of the aircraft are comprehensively considered, a low-altitude airspace route of the aircraft is represented as cubic grid matrix coded data, grid index coordinates and corresponding grid matrix data respectively describe space and time information of an airspace occupied by the aircraft, and therefore space-time information of the aircraft occupying the low-altitude airspace is described by three-dimensional grid matrix elements with different scales;
the method comprises the following specific steps of using multi-scale rasterization characterization in a low-altitude airspace:
(2.1) when the rasterized representation aircraft occupies airspace, the physical sizes of different aircraft need to be considered, and meanwhile, the safety interval of the aircraft during normal flight needs to be considered, so that various aircraft are combined
Figure SMS_100
Radius considered to take into account its physical dimensions and safety interval
Figure SMS_103
The envelope sphere of the envelope sphere, the course coordinate point is the sphere center of the envelope sphere, the flight path point of the aircraft expressed by a geodetic coordinate system is converted into an airspace coordinate point expressed by a Cartesian coordinate system in the corresponding low-altitude airspace, and the space grid occupation condition is judged by calculating the distance between the airspace coordinate point and the center point of the nearby grid; converting aircraft waypoints represented in a geodetic coordinate system into waypoint spatial coordinates represented in a Cartesian coordinate system corresponding to a spatial grid by coordinate system conversion
Figure SMS_107
And when the aircraft route point is represented in a rasterization manner, the distance between the grid edge and the aircraft route point is set
Figure SMS_101
Is less than
Figure SMS_105
The cubic grids are all regarded as space grids occupied by the air craft routes, wherein
Figure SMS_108
Representing coordinates as
Figure SMS_110
The grid of the space cube of (a),
Figure SMS_98
representing spatial coordinates
Figure SMS_104
The waypoints represented, for a scale of
Figure SMS_106
Figure SMS_109
) The low-altitude airspace grid and the aircraft
Figure SMS_99
Grid of course occupancy
Figure SMS_102
The following inequalities need to be satisfied:
Figure SMS_111
in the formula (I), the compound is shown in the specification,
Figure SMS_113
to represent
Figure SMS_116
Coordinate under scale of
Figure SMS_119
The grid of the space cube of (a),
Figure SMS_114
to represent
Figure SMS_117
Cubic grid size at scale, if course and occupied grid
Figure SMS_120
The spatial coordinates of the waypoint closest to the center are
Figure SMS_121
Then grid
Figure SMS_112
The corresponding time is the aircraft
Figure SMS_115
The waypoint is located at
Figure SMS_118
The time of (d);
determining the space grid occupancy, as shown in FIG. 3, which is a two-dimensional plan view of a space grid
Figure SMS_123
And
Figure SMS_125
two grids and flight route points
Figure SMS_128
Are respectively at a distance of
Figure SMS_124
And
Figure SMS_127
wherein the distance
Figure SMS_129
Figure SMS_130
Thus, therefore, it is
Figure SMS_122
Occupying a grid for the flight path of the aircraft, and
Figure SMS_126
do not belong to the aircraft track occupancy grid;
determining the time parameter of the flight path of the aircraft stored in the occupancy grid, as shown in FIG. 4, which is a two-dimensional plan view, for example, a space grid
Figure SMS_132
Belonging to aircraft route points at the same time
Figure SMS_135
Figure SMS_140
Figure SMS_133
Occupancy grid of (2), grid
Figure SMS_138
With aircraft waypoints
Figure SMS_141
Figure SMS_143
Figure SMS_131
Are respectively at a distance of
Figure SMS_137
Figure SMS_142
Figure SMS_144
Wherein
Figure SMS_134
Shortest, therefore, grid
Figure SMS_136
The time information expressed in the method is the waypoint
Figure SMS_139
The corresponding time;
(2.2) comparing the air route data with the low-altitude space domain grid matrix of different scales based on the air route data sequence of the low-altitude space domain aircraft, namely, representing the aircraft by using the element coordinates and the element numerical value sequence of the cubic grid matrixThe element coordinate of the grid matrix represents the airspace grid occupied by the aircraft route, the numerical value of the element represents the time of the aircraft route point reaching the airspace grid, and the specific coordinate point and the time represented in the grid are determined by the step (2.1) aiming at the situation that the time represented by the grid is used
Figure SMS_145
Scale-described low-altitude airspace aircraft represented by cubic grid matrix
Figure SMS_146
The course sequence of (a) is represented as:
Figure SMS_147
Figure SMS_148
Figure SMS_149
,……
in the formula (I), the compound is shown in the specification,
Figure SMS_152
representing aircraft
Figure SMS_154
In that
Figure SMS_156
The spatial grid coordinate at scale is
Figure SMS_151
Is determined to be the first waypoint of (c),
Figure SMS_155
representing the time corresponding to the waypoint;
Figure SMS_158
representing aircraft
Figure SMS_160
In that
Figure SMS_150
Spatial grid coordinates at scale
Figure SMS_153
Is the second waypoint of the flight path,
Figure SMS_157
representing the time corresponding to the waypoint; by parity of reasoning, the aircraft can be divided into two parts
Figure SMS_159
The flight path sequence is expressed by space grid coordinates and corresponding time;
(3) Constructing an airspace use conflict detection function, firstly judging the conflict situation of different aircraft use airspaces on the space based on the three-dimensional grid matrix description of the aircraft use airspace, and further considering the time conflict situation of the aircraft occupied in the airspace aiming at the airspace grid with conflict, thereby determining the space-time conflict range of the different aircraft use airspaces;
the conflict detection function is specifically constructed by the following steps:
(3.1) for different aircraft
Figure SMS_161
Detecting whether the flight path sequences of different aircrafts contain the same cubic grid coordinate, namely, a certain element coordinate represented by a matrix exists in the flight path sequences of the two aircrafts at the same time
Figure SMS_162
If the same cube grid coordinate exists, the fact that the corresponding aircraft air route conflicts at the grid is indicated;
(3.2) when the same grid matrix element coordinates exist in different aircraft
Figure SMS_164
Then, whether the time difference value corresponding to the element coordinate of different aircraft route sequences meets the safety time interval or not is further judged
Figure SMS_166
Figure SMS_171
Determining two aircraft routes if a set parameter related to the space grid size is satisfied
Figure SMS_165
And
Figure SMS_168
there is no conflict in the grid space between the grid spaces,
Figure SMS_170
and
Figure SMS_172
are respectively shown in
Figure SMS_163
And
Figure SMS_167
for two marked aircraft routes, collision detection function
Figure SMS_169
Outputting 0, otherwise indicating that two aircraft routes conflict in the grid space, and outputting 1 as a function, wherein the function is expressed as follows:
Figure SMS_173
in the formula (I), the compound is shown in the specification,
Figure SMS_174
is shown in
Figure SMS_175
The grid coordinates in the marked aircraft route are
Figure SMS_176
The time corresponding to the point of the waypoint of,
Figure SMS_177
is shown in
Figure SMS_178
The grid coordinates in the marked aircraft route are
Figure SMS_179
Time corresponding to the waypoint of (1);
(4) Designing an airspace conflict detection process, starting from a large-scale three-dimensional grid matrix on the basis of a multi-scale airspace use grid matrix, detecting conflict situations of different types of aircraft airspaces under the scale, and further detecting the conflict situations of the aircraft use airspace under a small-scale three-dimensional grid in a local conflict area aiming at the aircraft with conflicts, wherein the conflict situations are shown in figure 5;
the airspace conflict detection process specifically comprises the following steps:
(4.1) when the aircraft route conflict is detected, the detection flow is as shown in FIG. 6, firstly, the order of the orders is
Figure SMS_180
Starting the expressed aircraft route sequence, judging whether the routes of different aircrafts have conflicts under the scale, and if the routes of different aircrafts do not have conflicts, indicating that the corresponding aircrafts do not have conflicts when being used in the current airspace;
(4.2) for on-scale
Figure SMS_181
The spatial grid of the expressed flight path sequence with conflict is further judged according to the scale
Figure SMS_182
Judging whether the represented flight paths of the aircrafts in the conflict grid have conflicts or not, if not, indicating that the corresponding aircrafts are used in the current airspace without conflicts;
(4.3) for on-scale
Figure SMS_183
The space grid with conflict exists in the represented route sequence, and further judgment is carried out according to the scale
Figure SMS_184
And judging whether conflicts exist in the routes of different aircrafts under the scale according to the represented routes of the aircrafts in the conflict grid, if no conflict exists, indicating that no conflict exists in the use of the corresponding aircraft in the current airspace, otherwise, indicating that a conflict exists in the use of the corresponding aircraft in the current airspace.

Claims (6)

1. A low-altitude space domain use conflict detection method based on a multi-scale space grid is characterized by comprising the following steps:
(1) Performing airspace rasterization modeling and digital coding, namely performing discretization processing on a given continuous low-altitude airspace, dividing the continuous airspace into gapless cube grid models with different scales, wherein a cube grid of each scale can completely describe all the given low-altitude airspaces, and performing three-dimensional matrix type sequential coding on the cube grids of each scale to form a multi-scale low-altitude airspace grid coding system with codes mapped with the cube grids of different scales one by one, so as to form a low-altitude airspace discretization description mode;
(2) The method comprises the steps that a low-altitude airspace represents by using multi-scale rasterization, the size of an aircraft in the low-altitude airspace and a corresponding safety distance of the aircraft are considered, a low-altitude airspace route of the aircraft is represented as cubic grid matrix coded data, grid index coordinates and corresponding grid matrix data respectively describe space and time information of an airspace occupied by the aircraft, and space-time information of the aircraft occupying the low-altitude airspace is described by three-dimensional grid matrix elements with different scales;
(3) Constructing an airspace use conflict detection function, firstly judging the conflict situation of different aircraft use airspaces on the space based on the three-dimensional grid matrix description of the aircraft use airspace, and further considering the time conflict situation of the aircraft occupied in the airspace aiming at the airspace grid with conflict, thereby determining the space-time conflict range of the different aircraft use airspaces;
(4) Designing an airspace conflict detection process, starting from a large-scale three-dimensional grid matrix on the basis of a multi-scale airspace use grid matrix, detecting the conflict situation of the use of different types of aircraft airspaces under the scale, and further detecting the conflict situation of the use airspace of the aircraft under a small-scale three-dimensional grid in a local conflict region aiming at the aircraft with conflict.
2. The method for detecting the usage of the low-altitude space domain based on the multi-scale space grid according to claim 1, wherein: the given continuous low-altitude airspace in the step (1) is the altitude of
Figure QLYQS_1
Flight area below, sized above ground
Figure QLYQS_2
Left and right continuous low-altitude airspace.
3. The method for detecting the conflict of the use of the low-altitude space based on the multi-scale space grid as claimed in claim 2, wherein: the spatial domain rasterization coding of different scales in the step (1) specifically comprises the following steps:
(1.1) constructing a multi-scale space grid model of a given low-altitude airspace, dividing the given low-altitude airspace into cubic grids with different sizes and side lengths, realizing multi-scale gapless cubic grid modeling of a continuous low-altitude airspace, dividing the airspace into three scales, and respectively setting the sizes of the corresponding cubic grids to be three
Figure QLYQS_4
Figure QLYQS_8
And
Figure QLYQS_10
Figure QLYQS_5
the cubic grid of dimensions may include
Figure QLYQS_6
An
Figure QLYQS_9
Cubic grid of dimensions, likewise
Figure QLYQS_11
The size of the cubic grid can be contained in
Figure QLYQS_3
An
Figure QLYQS_7
A cubic grid of dimensions; selecting a Cartesian coordinate origin as a vertex of a low-altitude airspace, and indexing a low-altitude airspace cube grid divided by three scales through three-dimensional coordinates to convert aircraft route coordinates represented by a geodetic coordinate system into low-altitude airspace coordinates represented by the Cartesian coordinates;
(1.2) coding the low-altitude space domain cube grids divided by different scales in a three-dimensional matrix form, wherein for a given low-altitude space domain, the three-dimensional matrix dimensions corresponding to the low-altitude space domain grids divided by the three scales are respectively
Figure QLYQS_12
Figure QLYQS_13
Figure QLYQS_14
Coding low-altitude space domain grids divided in different scales according to a three-dimensional matrix element index mode, wherein cubic grids in the low-altitude space domain three-dimensional grid division corresponding to each scale can be obtained by corresponding three-dimensional matrix element coordinates
Figure QLYQS_15
And uniquely determining to realize matrix type low-altitude space domain discretization description.
4. The method for detecting the usage of the low-altitude space domain based on the multi-scale space grid according to claim 3, wherein: the specific steps of using multi-scale rasterization characterization for the low-altitude airspace in the step (2) are as follows:
(2.1) when the aircraft occupies airspace in a rasterization representation mode, the physical sizes of different aircraft need to be considered, and meanwhile, the safety interval of the aircraft in normal flight needs to be considered, so that various aircraft can be used
Figure QLYQS_17
Radius considered to take into account its physical dimensions and safety interval
Figure QLYQS_20
The waypoint coordinate point of the aircraft is the center of the circumscribed sphere, and the waypoint represented in the geodetic coordinate system is converted into the waypoint space coordinate represented in the Cartesian coordinate system corresponding to the space grid through coordinate system conversion
Figure QLYQS_23
And when the aircraft route point is represented in a rasterization mode, the distance between the grid edge and the aircraft route point is set
Figure QLYQS_19
Is less than
Figure QLYQS_22
The cubic grids are all regarded as space grids occupied by the air craft routes, wherein
Figure QLYQS_25
Representing coordinates as
Figure QLYQS_27
The grid of the space cube of (a),
Figure QLYQS_16
representing spatial coordinates
Figure QLYQS_21
The waypoints represented, for a scale of
Figure QLYQS_24
Figure QLYQS_26
) Low-altitude airspace grid of (1), grid occupied by course of aircraft
Figure QLYQS_18
The following inequalities need to be satisfied:
Figure QLYQS_28
in the formula (I), the compound is shown in the specification,
Figure QLYQS_30
represent
Figure QLYQS_33
Coordinate under scale of
Figure QLYQS_36
The grid of the space cube of (a),
Figure QLYQS_31
represent
Figure QLYQS_34
Cubic grid size at scale, if course and occupied grid
Figure QLYQS_37
The spatial coordinates of the waypoint closest to the center are
Figure QLYQS_38
Then grid
Figure QLYQS_29
The corresponding time is the aircraft
Figure QLYQS_32
The waypoint is located at
Figure QLYQS_35
The time of (d);
(2.2) comparing the flight path data with low-altitude space domain grid matrixes of different scales based on the airspace flight path data sequence of the low-altitude space domain aircraft, namely, representing the flight path of the aircraft by using element coordinates and element numerical value sequences of a cubic grid matrix, representing the airspace grid occupied by the flight path of the aircraft by using the element coordinates of the grid matrix, representing the time of the flight path point of the aircraft to the airspace grid by using the numerical value of the element, and aiming at the situation that the flight path point of the aircraft reaches the airspace grid by using the numerical value of the element
Figure QLYQS_39
Scale-described low-altitude airspace aircraft represented by cubic grid matrix
Figure QLYQS_40
Is represented as:
Figure QLYQS_41
Figure QLYQS_42
Figure QLYQS_43
,……
in the formula (I), the compound is shown in the specification,
Figure QLYQS_45
representing aircraft
Figure QLYQS_49
In that
Figure QLYQS_52
Spatial grid coordinate at scale of
Figure QLYQS_46
Is determined to be the first waypoint of (c),
Figure QLYQS_47
representing the time corresponding to the waypoint;
Figure QLYQS_50
representing aircraft
Figure QLYQS_53
In that
Figure QLYQS_44
The spatial grid coordinate at scale is
Figure QLYQS_48
Of the second waypoint of (a) is,
Figure QLYQS_51
representing the time corresponding to the waypoint; by parity of reasoning, the aircraft can be divided into two parts
Figure QLYQS_54
All represented in spatial grid coordinates and their corresponding times.
5. The method for detecting the usage of the low-altitude space domain based on the multi-scale space grid as claimed in claim 4, wherein: the collision detection function in the step (3) is specifically constructed by the following steps:
(3.1) for different aircraft
Figure QLYQS_55
Detecting whether the flight path sequences of different aircrafts contain the same cubic grid coordinate, namely, a certain element coordinate represented by a matrix exists in the flight path sequences of the two aircrafts at the same time
Figure QLYQS_56
(3.2) the same grid matrix element coordinates exist in different aircraft
Figure QLYQS_58
Then, whether the time difference value corresponding to the element coordinate of different aircraft route sequences meets the safety time interval or not is further judged
Figure QLYQS_60
Figure QLYQS_63
Determining two aircraft routes if the set parameters related to the space grid size are met
Figure QLYQS_59
And
Figure QLYQS_62
there is no conflict in the grid space between the grid spaces,
Figure QLYQS_65
and
Figure QLYQS_67
are respectively shown in
Figure QLYQS_57
And
Figure QLYQS_61
for two marked aircraft routes, collision detection function
Figure QLYQS_64
Outputting 0, otherwise, indicating that two aircraft routes conflict in the grid space, and outputting 1 by the function, wherein the function is represented as follows:
Figure QLYQS_66
in the formula (I), the compound is shown in the specification,
Figure QLYQS_68
is shown in
Figure QLYQS_69
The grid coordinates in the marked aircraft route are
Figure QLYQS_70
The time corresponding to the point of the waypoint of,
Figure QLYQS_71
is shown in
Figure QLYQS_72
The grid coordinates in the marked aircraft route are
Figure QLYQS_73
The waypoint of (a) corresponds to the time.
6. The method for detecting the usage of the low-altitude space domain based on the multi-scale space grid according to claim 5, wherein: the empty domain conflict detection process in the step (4) specifically comprises the following steps:
(4.1) first from the scale
Figure QLYQS_74
Starting the expressed aircraft route sequence, judging whether the routes of different aircrafts have conflicts under the scale, and if the routes of different aircrafts do not have conflicts, indicating that the corresponding aircrafts do not have conflicts when being used in the current airspace;
(4.2) for on-scale
Figure QLYQS_75
The spatial grid of the expressed flight path sequence with conflict is further judged according to the scale
Figure QLYQS_76
Judging whether the represented flight paths of the aircrafts in the conflict grid have conflicts or not, if not, indicating that the corresponding aircrafts are used in the current airspace without conflicts;
(4.3) needleTo the scale
Figure QLYQS_77
The spatial grid of the expressed flight path sequence with conflict is further judged according to the scale
Figure QLYQS_78
And judging whether conflicts exist in the routes of different aircrafts under the scale according to the represented routes of the aircrafts in the conflict grid, if no conflict exists, indicating that no conflict exists in the use of the corresponding aircraft in the current airspace, otherwise, indicating that a conflict exists in the use of the corresponding aircraft in the current airspace.
CN202310047397.9A 2023-01-31 2023-01-31 Low-altitude space domain use conflict detection method based on multi-scale space grid Pending CN115810087A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310047397.9A CN115810087A (en) 2023-01-31 2023-01-31 Low-altitude space domain use conflict detection method based on multi-scale space grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310047397.9A CN115810087A (en) 2023-01-31 2023-01-31 Low-altitude space domain use conflict detection method based on multi-scale space grid

Publications (1)

Publication Number Publication Date
CN115810087A true CN115810087A (en) 2023-03-17

Family

ID=85487671

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310047397.9A Pending CN115810087A (en) 2023-01-31 2023-01-31 Low-altitude space domain use conflict detection method based on multi-scale space grid

Country Status (1)

Country Link
CN (1) CN115810087A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117935626A (en) * 2024-03-25 2024-04-26 中国电子科技集团公司第二十八研究所 Size-adaptive low-altitude airspace use conflict detection method and system
CN117933675A (en) * 2024-03-25 2024-04-26 中国电子科技集团公司第二十八研究所 Rapid airspace regulating and controlling method and system based on grid divide and conquer
CN117935626B (en) * 2024-03-25 2024-05-28 中国电子科技集团公司第二十八研究所 Size-adaptive low-altitude airspace use conflict detection method and system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103116691A (en) * 2012-12-25 2013-05-22 中国人民解放军军事交通学院 Environment modeling device by means of variable scale grids
CN110782708A (en) * 2019-11-01 2020-02-11 南京智慧航空研究院有限公司 Unmanned aerial vehicle flight network modeling method based on low-altitude airspace limiting conditions
CN111477034A (en) * 2020-03-16 2020-07-31 中国电子科技集团公司第二十八研究所 Large-scale airspace use plan conflict detection and release method based on grid model
WO2021213540A1 (en) * 2020-09-23 2021-10-28 中国民航大学 Three-dimensional safe route planning method for unmanned aerial vehicle

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103116691A (en) * 2012-12-25 2013-05-22 中国人民解放军军事交通学院 Environment modeling device by means of variable scale grids
CN110782708A (en) * 2019-11-01 2020-02-11 南京智慧航空研究院有限公司 Unmanned aerial vehicle flight network modeling method based on low-altitude airspace limiting conditions
CN111477034A (en) * 2020-03-16 2020-07-31 中国电子科技集团公司第二十八研究所 Large-scale airspace use plan conflict detection and release method based on grid model
WO2021213540A1 (en) * 2020-09-23 2021-10-28 中国民航大学 Three-dimensional safe route planning method for unmanned aerial vehicle

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117935626A (en) * 2024-03-25 2024-04-26 中国电子科技集团公司第二十八研究所 Size-adaptive low-altitude airspace use conflict detection method and system
CN117933675A (en) * 2024-03-25 2024-04-26 中国电子科技集团公司第二十八研究所 Rapid airspace regulating and controlling method and system based on grid divide and conquer
CN117935626B (en) * 2024-03-25 2024-05-28 中国电子科技集团公司第二十八研究所 Size-adaptive low-altitude airspace use conflict detection method and system
CN117933675B (en) * 2024-03-25 2024-05-28 中国电子科技集团公司第二十八研究所 Rapid airspace regulating and controlling method and system based on grid divide and conquer

Similar Documents

Publication Publication Date Title
CN111477034B (en) Large-scale airspace use plan conflict detection and release method based on grid model
CN109541634B (en) Path planning method and device and mobile device
CN103247041B (en) A kind of dividing method of the cloud data of the many geometric properties based on local sampling
CN108204814B (en) Unmanned aerial vehicle three-dimensional scene path navigation platform and three-dimensional improved path planning method thereof
CN108805350B (en) Search and rescue range prediction method based on multi-dimensional Monte Carlo theory
CN104317772B (en) A kind of Monte Carlo PARTICLE TRANSPORT FROM fast geometric disposal route based on space lattice segmentation
CN103900573B (en) A kind of underwater research vehicle multiple constraint Route planner based on S57 standard electronic sea chart
CN107228673A (en) Route planner and device
CN115098993A (en) Unmanned aerial vehicle conflict detection method and device for airspace digital grid and storage medium
CN103745459A (en) Detection method of an unstructured point cloud feature point and extraction method thereof
CN114419281B (en) Method for calculating space geometric relation of airspace grid
CN111475597B (en) Non-rigid grid coding, spatial object unique identification and query method and device
CN116069882B (en) Airspace grid diagram generating method
Yang et al. Maritime traffic flow clustering analysis by density based trajectory clustering with noise
CN115661374B (en) Rapid retrieval method based on space division and model voxelization
Spiliopoulos et al. Knowledge extraction from maritime spatiotemporal data: An evaluation of clustering algorithms on Big Data
CN115329220A (en) Low-altitude spatial domain earth subdivision grid data organization and query method and device
CN114282796B (en) Method for calculating collision risk probability by airplane safety envelope based on airspace grid
Wu et al. A Non-rigid hierarchical discrete grid structure and its application to UAVs conflict detection and path planning
CN109191484B (en) A method of the rapidly extracting planar chip from airborne laser radar point cloud
CN115810087A (en) Low-altitude space domain use conflict detection method based on multi-scale space grid
CN113409408A (en) Airspace multilevel grid characterization and conflict detection method based on multi-way tree
CN117935626B (en) Size-adaptive low-altitude airspace use conflict detection method and system
CN117935626A (en) Size-adaptive low-altitude airspace use conflict detection method and system
CN106373189A (en) Space-time body data coding method based on linear hexadecimal tree

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