A tomographic imaging method for rapid small-scale ionospheric disturbance events

Abstract

The present application relates to a kind of tomography methods for fast small-scale ionospheric disturbance event, based on the disturbance event occurrence time period of target ionosphere in the disturbance event occurrence space above, and with the various navigation satellites to be analyzed of target layout area in the space link of GNSS receiver, according to the GNSS receiver site selection position coordinate range in target time period, analyze the space intersection of the local motion coordinate track of each navigation satellite to be analyzed and the space link between GNSS receiver in disturbance event occurrence space, and with maximum space intersection as target, determine the optimal site selection coordinate of GNSS receiver, then according to the space link communication between GNSS receiver and navigation satellite to be analyzed, complete the tomography of disturbance event in target ionosphere, to improve effective observation time, the contribution degree of GNSS receiver monitoring data in disturbance event occurrence space to ionospheric tomography, solve the problem that the number of receiver is limited in the application of sparse distribution in fast small-scale ionospheric disturbance event tomography, effective monitoring time is short, coverage area is not enough etc., can effectively improve the ability of ionospheric tomography.

Description

TECHNICAL FIELD

[0001] The application relates to a tomography method for rapid small-scale ionospheric disturbance events, and belongs to the technical field of space environment detection. BACKGROUND

[0002] The ionosphere is an ionized region in the near space of the earth with a height of 60-1000 km from the ground, which is ionized by high-energy radiation of the sun and cosmic rays. The ionosphere is in a partially ionized state and is a plasma. When an electromagnetic wave passes through the ionosphere, the electromagnetic wave interacts with the plasma, and effects such as refraction and reflection occur, which are closely related to the electron density parameter of the ionosphere. The ionospheric electron density is one of the most important parameters in the ionosphere, and effective means are needed to monitor the ionospheric electron density state in the upper space of key communication support areas, measurement and control stations, radar stations and the like for a long time and with high coverage, so as to avoid the adverse effects of the ionosphere on satellite communication, space measurement and control and radar detection. At present, the ionospheric electron density monitoring is mainly realized by combining GNSS (Global Navigation Satellite System) receiver monitoring data with tomographic inversion technology. The GNSS receiver receives satellite navigation signals, and the satellite navigation signals are affected by the ionosphere when passing through the ionosphere. The GNSS receiver samples the satellite navigation signals at a high frequency, and the total electron content of the slant path, i.e. the integral quantity of the electron density on the path, can be calculated according to the amplitude and phase changes of the sampled signals. The GNSS receiver receives the data of the satellite-ground link of multiple navigation satellites at the same time, and obtains the total electron content on multiple slant paths, which can realize the monitoring of the ionospheric electron density in the region by combining the tomographic inversion technology.

[0003] However, the accuracy of the electron density data obtained by this monitoring method is seriously dependent on the distribution of the satellite-ground link between the receivers and the navigation satellites in the three-dimensional space of the monitored region. Especially for the monitoring of the ionospheric electron density disturbance of rapid and small-scale events such as ionospheric scintillation, the GNSS monitoring network such as the IGS (International GNSS Service) monitoring station network and the national earthquake monitoring network is too sparse, and the data resolution is usually more than 200 km, which is difficult to meet the monitoring needs of rapid and small-scale ionospheric disturbance events. Especially in the case of limited number of receivers, a corresponding data acquisition method is needed to face the scale of grid division in inversion, and the receiver layout position is selected in a small range near the region that needs to be protected, so as to improve the contribution rate of the receiver monitoring data to the ionospheric tomography in the target event and the target region, and further improve the monitoring efficiency. SUMMARY

[0004] The technical problem solved by the present application is to provide a tomographic imaging method for rapid small-scale ionospheric disturbance events, which is used to realize efficient tomographic imaging of the space of the target ionospheric disturbance events under the condition of a single GNSS receiver, and further effectively improve the monitoring capability of the regional ionospheric electron density.

[0005] In order to solve the above technical problems, the present application adopts the following technical scheme: the present application designs a tomographic imaging method for rapid small-scale ionospheric disturbance events, which performs the following steps A to D according to the target layout area for GNSS receiver site selection and layout, and realizes tomographic imaging of the target ionospheric disturbance events.

[0006] Step A. According to the navigation constellation ephemeris of each navigation satellite, the present application obtains each navigation satellite that can establish a satellite-ground link with the GNSS receiver in the target layout area and passes above the space of the target ionospheric disturbance events during the time period of the occurrence of the target ionospheric disturbance events, as each to-be-analyzed navigation satellite, and obtains the corresponding motion trajectory of each to-be-analyzed navigation satellite, and then enters step B.

[0007] Step B. Based on the site selection range of the GNSS receiver in the target layout area, the present application converts the corresponding motion trajectory of each to-be-analyzed navigation satellite and the site selection of the GNSS receiver in the target layout area into the same coordinate system, obtains the corresponding motion coordinate trajectory of each to-be-analyzed navigation satellite and the site selection coordinate range of the GNSS receiver in the target layout area in the same coordinate system, and then enters step C.

[0008] Step C. Based on the local motion coordinate trajectory of each to-be-analyzed navigation satellite corresponding to the time period of the occurrence of the target ionospheric disturbance events, the present application measures the spatial intersection of the local motion coordinate trajectory of each to-be-analyzed navigation satellite and the satellite-ground link of the GNSS receiver in the target ionospheric disturbance event space under each site selection coordinate in the site selection coordinate range of the GNSS receiver, and determines the optimal site selection coordinate of the GNSS receiver with the goal of maximizing the spatial intersection, and performs the layout of the GNSS receiver in the target layout area according to the optimal site selection coordinate, and then enters step D.

[0009] Step D. Based on the satellite-ground link communication between the GNSS receiver laid in the target layout area and each to-be-analyzed navigation satellite, the present application applies a CT tomographic algorithm to tomographically invert the three-dimensional grid space of the target ionospheric disturbance events, obtains the three-dimensional distribution data of the electron density of the three-dimensional grid space of the disturbance events, and realizes tomographic imaging of the target ionospheric disturbance events.

[0010] As a preferred technical scheme of the present application, the step A includes the following steps A1 to A3.

[0011] Step A1, define the direction of each position on the edge of the target deployment area relative to the center position of the target deployment area as the reference direction, take each position as the starting point, and draw a ray through the edge position on the edge of the top surface of the target ionosphere disturbance event occurrence space opposite the reference direction relative to the center position of the top surface, to obtain a ray with each position on the edge of the target deployment area as the starting point, and then proceed to step A2;

[0012] Step A2, define the area surrounded by all rays and located above the target ionosphere disturbance event occurrence space as the analysis area, and obtain each navigation satellite present in the analysis area during the target ionosphere disturbance event occurrence period according to the navigation constellation ephemeris of each navigation satellite, as each analysis navigation satellite, and then proceed to step A3;

[0013] Step A3, obtain the motion trajectory corresponding to each analysis navigation satellite.

[0014] As a preferred technical solution of the present application: in step B, taking the position directly below the center position of the target ionosphere disturbance event occurrence space corresponding to the height of the target deployment area as the origin o, the two orthogonal directions in the horizontal plane as the x coordinate axis and the y coordinate axis, and the vertical upward direction as the z coordinate axis, an o-xyz coordinate system is constructed, and the motion trajectory corresponding to each analysis navigation satellite and the site selection of the GNSS receiver in the target deployment area are converted to the o-xyz coordinate system to obtain the motion coordinate trajectory corresponding to each analysis navigation satellite and the site selection coordinate range of the GNSS receiver in the target deployment area in the o-xyz coordinate system.

[0015] As a preferred technical solution of the present application: step C includes steps C1-1 to C1-2.

[0016] Step C1-1, based on the disturbance event three-dimensional grid space under the three-dimensional grid division of the target ionosphere disturbance event occurrence space, obtain the local motion coordinate trajectory corresponding to the motion coordinate trajectory of each analysis navigation satellite during the target ionosphere disturbance event occurrence period, construct the motion trajectory surface function between each site selection coordinate in the GNSS receiver site selection coordinate range and the local motion coordinate trajectory, and construct the way grid statistical function of the motion trajectory surface passing through the grid number in the disturbance event three-dimensional grid space, and then obtain the way grid statistical function corresponding to each analysis navigation satellite, and then proceed to step C1-2.

[0017] Step C1-2. Construct a summation function for the path grid statistical function corresponding to each navigation satellite to be analyzed, calculate the optimal address coordinates of the GNSS receiver that maximize the result of the summation function, and deploy the GNSS receiver in the target deployment area according to the optimal address coordinates, and then proceed to step D.

[0018] As a preferred technical solution of the present invention: step C includes the following steps C2-1 to C2-4;

[0019] Step C2-1. Obtain the local motion coordinate trajectory of each navigation satellite to be analyzed within the time period of the disturbance event in the target ionosphere, that is, obtain the local motion coordinate trajectory corresponding to each navigation satellite to be analyzed, and then proceed to step C2-2.

[0020] Step C2-2. Based on the GNSS receiver's location coordinate range, define the coordinates of the GNSS receiver's location in the target deployment area as (Gx, Gy, 0). For each navigation satellite to be analyzed, based on the local motion coordinate trajectory of the navigation satellite to be analyzed, obtain the starting point A (Ax, Ay, Az) and the ending point A' (A'x, A'y, A'z) of the space top surface where the disturbance event occurs in the ionosphere of the target, corresponding to the local motion coordinate trajectory and (Gx, Gy, 0). The plane containing the starting point A (Ax, Ay, Az), the ending point A' (A'x, A'y, A'z), and the GNSS receiver coordinates (Gx, Gy, 0) is used as the analysis plane corresponding to the navigation satellite to be analyzed. Construct the analysis plane function corresponding to the navigation satellite to be analyzed, and then construct the analysis plane function corresponding to each navigation satellite to be analyzed. Then proceed to step C2-3.

[0021] Step C2-3. Based on the three-dimensional grid space of the disturbance event occurrence space in the target ionosphere, for each navigation satellite to be analyzed, based on the analysis plane function corresponding to the navigation satellite to be analyzed, construct the intersection grid statistical function of the number of intersection grids between the analysis plane and the three-dimensional grid space of the disturbance event, and then obtain the intersection grid statistical function corresponding to each navigation satellite to be analyzed, and then proceed to step C2-4.

[0022] Step C2-4. Construct a summation function for the intersection grid statistical function corresponding to each navigation satellite to be analyzed, calculate the optimal address coordinates of the GNSS receiver that maximize the result of the summation function, and deploy the GNSS receiver in the target deployment area according to the optimal address coordinates, and then proceed to step D.

[0023] The tomographic imaging method for rapid small-scale ionospheric disturbance events described in this invention has the following technical advantages compared with existing technologies:

[0024] This invention presents a tomographic imaging method for rapid small-scale ionospheric disturbance events. Based on the time period of a disturbance event in the target ionosphere, it analyzes the spatial intersection of the local motion coordinate trajectories of each satellite and the satellite-to-ground link between the satellite and the GNSS receiver within the target deployment area. With the goal of maximizing this spatial intersection, the optimal location coordinates of the GNSS receiver are determined. Then, based on the communication between the GNSS receiver and the satellite, tomographic imaging of the disturbance event in the target ionosphere is completed. This improves the effective observation time and the contribution of GNSS receiver monitoring data to ionospheric tomography in the disturbance event space. It solves the problems of short effective monitoring time and insufficient coverage area in the application of rapid small-scale ionospheric disturbance event tomography with sparsely deployed stations when the number of receivers is limited, effectively improving the ionospheric tomographic imaging capability. Attached Figure Description

[0025] Figure 1 This is a schematic diagram illustrating the application of the tomographic imaging method for rapid small-scale ionospheric disturbance events designed in this invention. Detailed Implementation

[0026] The specific embodiments of the present invention will be further described in detail below with reference to the accompanying drawings.

[0027] This invention designs a tomographic imaging method for rapid, small-scale ionospheric perturbation events. In practical applications, it is specifically designed for... Figure 1 In the scenario shown, based on the target deployment area used for GNSS receiver site selection and deployment, the following steps A to D are performed to achieve tomographic imaging of disturbance events in the target ionosphere.

[0028] Step A. Based on the ephemeris of each navigation satellite constellation, obtain the navigation satellites that pass above the space where the disturbance event occurred in the target ionosphere during the time period of the disturbance event and can establish a satellite-to-ground link with the GNSS receiver in the target deployment area. These satellites are identified as the navigation satellites to be analyzed. The motion trajectories of each navigation satellite to be analyzed are obtained. Then proceed to Step B.

[0029] In practice, step A above is specifically implemented as follows: steps A1 to A3.

[0030] Step A1. For each continuous position on the edge of the target deployment area, define the direction of that position relative to the center of the target deployment area as the reference direction. Starting from that position, draw a ray through the edge position on the side opposite to the center of the top surface of the space where the disturbance event occurs in the target ionosphere. Obtain rays starting from each continuous position on the edge of the target deployment area, and then proceed to step A2.

[0031] Step A2. Define the region surrounded by all rays and located above the space where the disturbance event occurred in the target ionosphere as the region to be analyzed. Based on the navigation constellation ephemeris of each navigation satellite, obtain each navigation satellite that appeared in the region to be analyzed during the time period of the disturbance event in the target ionosphere, and then proceed to step A3.

[0032] Step A3. Obtain the motion trajectory of each navigation satellite to be analyzed.

[0033] Step B. Design an o-xyz coordinate system with the origin o located directly below the spatial center of the disturbance event in the target ionosphere and at the height of the target deployment area. The two orthogonal directions in the horizontal plane are the x-axis and y-axis, and the vertical direction upward is the z-axis. Then, transform the motion trajectories of each navigation satellite to be analyzed and the location of the GNSS receiver in the target deployment area into the o-xyz coordinate system to obtain the motion coordinate trajectories of each navigation satellite to be analyzed and the location coordinate range of the GNSS receiver in the target deployment area under the o-xyz coordinate system. Then proceed to step C.

[0034] Step C. Based on the local motion coordinate trajectories of each navigation satellite to be analyzed within the time period of the disturbance event in the target ionosphere, measure the spatial intersection of the local motion coordinate trajectories of each navigation satellite to be analyzed and the satellite-to-ground link between the GNSS receiver and the disturbance event in the target ionosphere at each selected location coordinate. With the goal of maximizing spatial intersection, determine the optimal location coordinates of the GNSS receiver, and deploy the GNSS receiver in the target deployment area according to the optimal location coordinates. Then proceed to step D.

[0035] Regarding step C above, two specific implementations are designed for practical application. One implementation specifically executes steps C1-1 to C1-2 as follows.

[0036] Step C1-1. Based on the three-dimensional grid space of the disturbance event in the target ionosphere, which is divided into three-dimensional grids, for each navigation satellite to be analyzed, obtain the local motion coordinate trajectory of the navigation satellite to be analyzed within the time period of the disturbance event in the target ionosphere. Construct the motion trajectory surface function between each site selection coordinate in the GNSS receiver site selection coordinate range and the local motion coordinate trajectory, and construct the path grid statistical function of the number of grids in the three-dimensional grid space of the disturbance event through the motion trajectory surface. Then obtain the path grid statistical function corresponding to each navigation satellite to be analyzed, and then proceed to step C1-2.

[0037] Step C1-2. Construct a summation function for the path grid statistical function corresponding to each navigation satellite to be analyzed, calculate the optimal address coordinates of the GNSS receiver that maximize the result of the summation function, and deploy the GNSS receiver in the target deployment area according to the optimal address coordinates, and then proceed to step D.

[0038] The second embodiment specifically executes the following steps C2-1 to C2-4.

[0039] Step C2-1. Obtain the local motion coordinate trajectory of each navigation satellite to be analyzed within the time period of the disturbance event in the target ionosphere, that is, obtain the local motion coordinate trajectory corresponding to each navigation satellite to be analyzed, and then proceed to step C2-2.

[0040] Step C2-2. Based on the GNSS receiver's location coordinate range, define the coordinates of the GNSS receiver's location in the target deployment area as (Gx, Gy, 0). For each navigation satellite to be analyzed, based on the local motion coordinate trajectory of the navigation satellite to be analyzed, obtain the starting point A (Ax, Ay, Az) and the ending point A' (A'x, A'y, A'z) of the space top surface where the disturbance event occurs in the ionosphere of the target, corresponding to the local motion coordinate trajectory and (Gx, Gy, 0). The plane containing the starting point A (Ax, Ay, Az), the ending point A' (A'x, A'y, A'z), and the GNSS receiver coordinates (Gx, Gy, 0) is used as the analysis plane corresponding to the navigation satellite to be analyzed. Construct the analysis plane function corresponding to the navigation satellite to be analyzed, and then construct the analysis plane function corresponding to each navigation satellite to be analyzed. Then proceed to step C2-3.

[0041] Here is the coordinate range for GNSS receiver location: , , These represent the minimum value in the negative x-axis and the maximum value in the positive and negative x-axis directions of the target deployment area in the o-xyz coordinate system, respectively. These represent the minimum value in the negative direction and the maximum value in the positive and negative directions of the y-axis in the o-xyz coordinate system corresponding to the target deployment area, respectively.

[0042] Step C2-3. Based on the three-dimensional grid space of the disturbance event occurrence space in the target ionosphere, for each navigation satellite to be analyzed, based on the analysis plane function corresponding to the navigation satellite to be analyzed, construct the intersection grid statistical function of the number of intersection grids between the analysis plane and the three-dimensional grid space of the disturbance event, and then obtain the intersection grid statistical function corresponding to each navigation satellite to be analyzed, and then proceed to step C2-4.

[0043] Step C2-4. Construct a summation function for the intersection grid statistical function corresponding to each navigation satellite to be analyzed, calculate the optimal address coordinates of the GNSS receiver that maximize the result of the summation function, and deploy the GNSS receiver in the target deployment area according to the optimal address coordinates, and then proceed to step D.

[0044] In the above embodiments, the space where disturbance events occur in the target ionosphere is within a vertical range of 60km to 1000km. The size of the three-dimensional grid space for disturbance events in the target ionosphere is determined according to the required analysis accuracy. Furthermore, regarding the spatial intersection of the local motion coordinate trajectories of each navigation satellite to be analyzed and the satellite-to-ground link between the GNSS receiver and the disturbance event space in the target ionosphere, i.e., regarding the three-dimensional grid space for disturbance events, the goal is to maximize the number of covered grids. The more grids covered, the greater the effect of the GNSS receiver on CT tomography of the space where disturbance events occur in the target ionosphere.

[0045] Step D. Based on the satellite-to-ground link communication between the GNSS receivers deployed in the target deployment area and each navigation satellite to be analyzed, the CT tomography algorithm is applied to perform tomographic inversion on the three-dimensional grid space of disturbance events in the target ionosphere to obtain the three-dimensional distribution data of electron density in the three-dimensional grid space of disturbance events, thus realizing tomographic imaging of disturbance events in the target ionosphere.

[0046] The aforementioned technical solution presents a tomographic imaging method for rapid small-scale ionospheric disturbance events. Based on the time period of the disturbance event in the target ionosphere, it analyzes the spatial intersection of the local motion coordinate trajectories of each satellite and the satellite-to-ground link between the satellite and the GNSS receiver within the target deployment area. With the goal of maximizing spatial intersection, it determines the optimal location coordinates of the GNSS receiver. Then, based on the communication between the GNSS receiver and the satellite, it performs tomographic imaging of the disturbance event in the target ionosphere. This improves the effective observation time and the contribution of GNSS receiver monitoring data to ionospheric tomography in the disturbance event space. It solves the problems of short effective monitoring time and insufficient coverage area in the application of rapid small-scale ionospheric disturbance event tomography when the number of receivers is limited, effectively improving the ionospheric tomographic imaging capability.

[0047] The embodiments of the present invention have been described in detail above with reference to the accompanying drawings. However, the present invention is not limited to the above embodiments. Within the scope of knowledge possessed by those skilled in the art, various changes can be made without departing from the spirit of the present invention.

Claims

1. A tomographic imaging method for fast, small-scale ionospheric perturbation events, characterized in that: Based on the target deployment area used for GNSS receiver site selection, perform the following steps A to D to achieve tomographic imaging of disturbance events in the target ionosphere. Step A. Based on the ephemeris of each navigation satellite constellation, obtain each navigation satellite that passes above the space where the disturbance event occurred in the target ionosphere during the time period of the disturbance event and can establish a satellite-to-ground link with the GNSS receiver in the target deployment area. These satellites are selected as the navigation satellites to be analyzed. The motion trajectories of each navigation satellite to be analyzed are obtained. Then proceed to Step B. Step B. Based on the location range of the GNSS receiver in the target deployment area, transform the motion trajectory of each navigation satellite to be analyzed and the location of the GNSS receiver in the target deployment area to the same coordinate system, obtain the motion coordinate trajectory of each navigation satellite to be analyzed and the location coordinate range of the GNSS receiver in the target deployment area under the same coordinate system, and then proceed to step C; Step C. Based on the local motion coordinate trajectories of each navigation satellite to be analyzed within the time period of the disturbance event in the target ionosphere, measure the spatial intersection of the local motion coordinate trajectories of each navigation satellite to be analyzed and the satellite-to-ground link between the GNSS receiver and the GNSS receiver in the space where the disturbance event occurs in the target ionosphere, and determine the optimal location coordinates of the GNSS receiver with the goal of maximizing spatial intersection. Then, deploy the GNSS receiver in the target deployment area according to the optimal location coordinates, and then proceed to step D. Step D. Based on the satellite-to-ground link communication between the GNSS receivers deployed in the target deployment area and each navigation satellite to be analyzed, the CT tomography algorithm is applied to perform tomographic inversion on the three-dimensional grid space of disturbance events in the target ionosphere to obtain the three-dimensional distribution data of electron density in the three-dimensional grid space of disturbance events, thus realizing tomographic imaging of disturbance events in the target ionosphere.

2. The tomographic imaging method for fast small-scale ionospheric perturbation events according to claim 1, characterized in that: Step A includes the following steps A1 to A3; Step A1. For each continuous position on the edge of the target deployment area, define the direction of that position relative to the center of the target deployment area as the reference direction. Starting from that position, draw a ray through the edge position on the side opposite to the center of the top surface of the space where the disturbance event occurs in the target ionosphere to obtain rays starting from each continuous position on the edge of the target deployment area. Then proceed to step A2. Step A2. Define the region surrounded by all rays and located above the space where the disturbance event occurred in the target ionosphere as the region to be analyzed. Based on the navigation constellation ephemeris of each navigation satellite, obtain each navigation satellite that appeared in the region to be analyzed during the time period of the disturbance event in the target ionosphere, and then proceed to step A3. Step A3. Obtain the motion trajectory of each navigation satellite to be analyzed.

3. The tomographic imaging method for fast small-scale ionospheric perturbation events according to claim 1, characterized in that: In step B, the origin o is defined as the position directly below the spatial center of the disturbance event in the target ionosphere, corresponding to the height of the target deployment area. The two orthogonal directions in the horizontal plane are the x-axis and y-axis, respectively, and the vertical upward direction is the z-axis. An o-xyz coordinate system is constructed, and the motion trajectories of each navigation satellite to be analyzed and the location of the GNSS receiver in the target deployment area are transformed into the o-xyz coordinate system. The motion coordinate trajectories of each navigation satellite to be analyzed and the location coordinate range of the GNSS receiver in the target deployment area are obtained in the o-xyz coordinate system.

4. The tomographic imaging method for fast small-scale ionospheric perturbation events according to claim 1, characterized in that: Step C includes the following steps C1-1 to C1-2; Step C1-1. Based on the three-dimensional grid space of the disturbance event in the target ionosphere, which is divided into three-dimensional grids, for each navigation satellite to be analyzed, obtain the local motion coordinate trajectory of the navigation satellite to be analyzed within the time period of the disturbance event in the target ionosphere. Construct the motion trajectory surface function between each site selection coordinate in the GNSS receiver site selection coordinate range and the local motion coordinate trajectory, and construct the path grid statistical function of the number of grids in the three-dimensional grid space of the disturbance event through the motion trajectory surface. Then obtain the path grid statistical function corresponding to each navigation satellite to be analyzed, and then proceed to step C1-2. Step C1-2. Construct a summation function for the path grid statistical function corresponding to each navigation satellite to be analyzed, calculate the optimal address coordinates of the GNSS receiver that maximize the result of the summation function, and deploy the GNSS receiver in the target deployment area according to the optimal address coordinates, and then proceed to step D.

5. The tomographic imaging method for fast small-scale ionospheric perturbation events according to claim 1, characterized in that: Step C includes the following steps C2-1 to C2-4; Step C2-1. Obtain the local motion coordinate trajectory of each navigation satellite to be analyzed within the time period of the disturbance event in the target ionosphere, that is, obtain the local motion coordinate trajectory corresponding to each navigation satellite to be analyzed, and then proceed to step C2-2. Step C2-2. Based on the GNSS receiver's location coordinate range, define the coordinates of the GNSS receiver's location in the target deployment area as (Gx, Gy, 0). For each navigation satellite to be analyzed, based on the local motion coordinate trajectory of the navigation satellite to be analyzed, obtain the starting point A (Ax, Ay, Az) and the ending point A' (A'x, A'y, A'z) of the space top surface where the disturbance event occurs in the ionosphere of the target, corresponding to the local motion coordinate trajectory and (Gx, Gy, 0). The plane containing the starting point A (Ax, Ay, Az), the ending point A' (A'x, A'y, A'z), and the GNSS receiver coordinates (Gx, Gy, 0) is used as the analysis plane corresponding to the navigation satellite to be analyzed. Construct the analysis plane function corresponding to the navigation satellite to be analyzed, and then construct the analysis plane function corresponding to each navigation satellite to be analyzed. Then proceed to step C2-3. Step C2-3. Based on the three-dimensional grid space of the disturbance event occurrence space in the target ionosphere, for each navigation satellite to be analyzed, based on the analysis plane function corresponding to the navigation satellite to be analyzed, construct the intersection grid statistical function of the number of intersection grids between the analysis plane and the three-dimensional grid space of the disturbance event, and then obtain the intersection grid statistical function corresponding to each navigation satellite to be analyzed, and then proceed to step C2-4. Step C2-4. Construct a summation function for the intersection grid statistical function corresponding to each navigation satellite to be analyzed, calculate the optimal address coordinates of the GNSS receiver that maximize the result of the summation function, and deploy the GNSS receiver in the target deployment area according to the optimal address coordinates, and then proceed to step D.

Images