A ground-based integrated method and device for joint detection of the ionosphere and troposphere.

Abstract

This invention discloses a ground-based integrated joint detection method and device for the ionosphere and troposphere, belonging to the field of satellite navigation technology. It includes: detecting the vertical total electron content of the ionosphere, including: extracting the tilted total electron content of the ionosphere at the ionospheric puncture point and converting it to the vertical total electron content of the ionosphere; detecting the ground-based BeiDou ionospheric disturbance index; detecting the ionospheric phase scintillation index; detecting the tropospheric zenith total delay; detecting tropospheric precipitable water; detecting the ionospheric amplitude scintillation index; and outputting integrated joint detection parameters of the ground-based BeiDou space ionosphere and troposphere. This invention effectively improves the reliability, robustness, and meteorological perception capabilities of navigation and positioning systems in complex atmospheric environments, and is particularly suitable for special environmental conditions such as severe convective weather, typhoon landfall, and severe ionospheric disturbances.

Description

Technical Field

[0001] This invention belongs to the field of satellite navigation technology, and in particular relates to a ground-based integrated method and device for joint detection of the ionosphere and troposphere. Background Technology

[0002] The ionosphere and troposphere, as important components of the Earth's atmospheric system, significantly impact satellite navigation, precise positioning, remote sensing, and weather forecasting. Changes in the distribution of free electrons in the ionosphere directly affect the propagation path and speed of GNSS signals. The total electron content (TEC) of the ionosphere is a core parameter for measuring the ionospheric state. TEC can be retrieved through multi-frequency GNSS carrier phase differential analysis to describe the degree of ionospheric delay. Furthermore, spatial irregularities or small-scale disturbances within the ionosphere can cause rapid changes in the signal propagation path, manifesting as phase scintillation and amplitude scintillation, corresponding to rapid fluctuations in phase and signal strength, respectively. A typical indicator of these phenomena is the phase scintillation index (PI). ) and amplitude flicker index ( These disturbances, when severe, can lead to receiver signal loss, frequent cycle slips, degraded positioning performance, and even interruption of navigation services.

[0003] Correspondingly, tropospheric delay is another key error source in high-precision GNSS positioning and meteorological inversion. The troposphere introduces a delay of approximately 2-3 meters into GNSS signals. This delay can be decomposed into dry delay (ZHD) and wet delay (ZWD), the sum of which is the total zenith delay (ZTD). In particular, the ZWD component is closely related to water vapor in the atmosphere, and precipitable water volume (PWV) can be further derived from ZWD, which has important application value in numerical weather prediction and precipitation monitoring.

[0004] Currently, mainstream technologies often separate the detection processes of the ionosphere and troposphere, establishing different models and processing procedures for each. For example, ionospheric TEC is mostly based on carrier phase difference decomposition, while ROTI is used for ionospheric disturbance monitoring. The index reflects amplitude scintillation intensity; while ZTD and PWV rely more on independent calculations from GNSS weather stations. This "parameter isolation" approach ignores the correlation between ionospheric disturbances and tropospheric water vapor changes in the actual atmospheric system, leading to anomalies or mutual interference among various parameters during severe weather events (such as rainfall, convection, and solar activity). In high PWV environments or heavy rainfall events, the signal-to-noise ratio of the receiver signal decreases significantly due to the attenuation effect of the cloud and rain layer. The index exhibits artificially high values ​​not caused by the ionosphere, which can easily be misinterpreted as strong ionospheric disturbances. Simultaneously, strong ionospheric scintillation events can degrade the quality of observational data, interfere with the PPP solution process, make ZTD estimation unstable, and consequently affect the accuracy of PWV derivation. This cross-layer interference effect is particularly pronounced in tropical regions or during the summer monsoon season, when ionospheric disturbances are frequent and water vapor content is abundant, leading to a significant decrease in the reliability and robustness of existing methods under extreme environments. Summary of the Invention

[0005] To address the aforementioned technical problems, this invention provides an integrated method and apparatus for joint detection of the ionosphere and troposphere in ground-based space, aiming to solve problems such as separate estimation of ionospheric and tropospheric parameters, insufficient information utilization, and poor model correlation in the prior art.

[0006] The specific technical solution is as follows:

[0007] A ground-based method for integrated detection of the ionosphere and troposphere includes the following steps:

[0008] Step S1: Obtain the total electron content of the tilted ionosphere and convert it to the total electron content of the vertical ionosphere;

[0009] Step S2: Calculate the ionospheric disturbance index of each satellite using the total tilted electron content;

[0010] Step S3: Detect the ionospheric phase scintillation index;

[0011] Step S4: Extract the total tropospheric zenith delay using the ionospheric disturbance index and the ionospheric phase scintillation index;

[0012] Step S5: Calculate the dry delay and conversion factor of the troposphere, obtain the wet delay of the troposphere using the total tropospheric zenith delay, and extract the precipitable water in the troposphere by combining the conversion factor.

[0013] Step S6: Calculate the original ionospheric amplitude scintillation information and correct the original ionospheric amplitude scintillation information using precipitable water volume (PWV) to obtain the ionospheric amplitude scintillation index.

[0014] Step S7: Output ground-based BeiDou space-based ionospheric and tropospheric integrated joint detection parameters, including: total vertical electron content of the ionosphere, ionospheric disturbance index, ionospheric phase scintillation index, ionospheric amplitude scintillation index, total zenith delay of the troposphere, and tropospheric precipitable water.

[0015] A ground-based integrated joint detection device for the ionosphere and troposphere includes:

[0016] The ionospheric vertical total electron content detection module acquires the tilted total electron content of the ionosphere and converts it to the ionospheric vertical total electron content.

[0017] The disturbance index detection module calculates the ionospheric disturbance index of each satellite using the total tilted electron content and obtains the ionospheric disturbance index.

[0018] Phase scintillation index detection module, used to calculate the ionospheric phase scintillation index;

[0019] The tropospheric zenith total delay detection module, combined with the ionospheric disturbance index and the ionospheric phase scintillation index, extracts the tropospheric zenith total delay.

[0020] The tropospheric precipitable water detection module calculates the tropospheric dry delay and conversion factor, obtains the tropospheric wet delay using the total tropospheric zenith delay, and extracts the tropospheric precipitable water by combining the conversion factor.

[0021] The ionospheric amplitude scintillation index detection module calculates the original ionospheric amplitude scintillation information and corrects it using precipitable water to obtain the ionospheric amplitude scintillation index.

[0022] The detection parameter output module outputs ground-based BeiDou space-based integrated joint detection parameters of the ionosphere and troposphere. The parameters include: total vertical electron content of the ionosphere, ionospheric disturbance index, ionospheric phase scintillation index, ionospheric amplitude scintillation index, total zenith delay of the troposphere, and tropospheric precipitable water.

[0023] An electronic device includes: one or more processors; and a memory for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors cause the one or more processors to implement the method.

[0024] A computer-readable storage medium having executable instructions stored thereon, which, when executed by a processor, cause the processor to implement the method described thereon.

[0025] The beneficial effects of this invention are as follows:

[0026] This invention achieves unified joint estimation of the total electron content (TEC) of the ionosphere, ionospheric disturbance index, ionospheric amplitude scintillation index, ionospheric phase scintillation index, total tropospheric zenith delay, and tropospheric precipitable water. This invention optimizes the estimation strategy for the total tropospheric zenith delay parameter using real-time estimation results of the ionospheric disturbance index and the ionospheric phase scintillation index. This invention also optimizes the ionospheric amplitude scintillation detection strategy using the estimation results of the tropospheric precipitable water parameter.

[0027] This invention effectively solves the problems of parameter separation, information isolation, and difficulty in identifying error coupling in traditional methods. This method can significantly improve GNSS positioning stability, meteorological parameter inversion accuracy, and space weather monitoring capabilities under complex atmospheric disturbance conditions.

[0028] This invention effectively improves the reliability, robustness, and meteorological perception capabilities of navigation and positioning systems in complex atmospheric environments, and is particularly suitable for special environmental conditions such as severe convective weather, typhoon landfall, and severe ionospheric disturbances. Attached Figure Description

[0029] Figure 1 A flowchart of a ground-based integrated joint detection method for the ionosphere and troposphere;

[0030] Figure 2 The result is the receiver differential code bias estimation.

[0031] Figure 3 The results represent the estimation of the total vertical electron content of the ionosphere.

[0032] Figure 4 The results are the estimation results for the ROTI index of ionospheric disturbance.

[0033] Figure 5 The results are the estimation results of the ionospheric phase scintillation index;

[0034] Figure 6 The results are the ZTD (Zero-Temperature Delay) estimates for the tropospheric zenith.

[0035] Figure 7 The results are the estimates of tropospheric precipitable water (PWV).

[0036] Figure 8 The results are the estimation results of the ionospheric amplitude scintillation index. Detailed Implementation

[0037] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention. Furthermore, the technical features involved in the various embodiments of this invention described below can be combined with each other as long as they do not conflict with each other.

[0038] like Figure 1 As shown, the integrated ground-based space ionosphere and troposphere detection method of the present invention specifically includes the following steps:

[0039] Step S1: Detect the vertical total electron content (VTEC) of the ionosphere, including: extracting the tilted total electron content (STEC) of the ionosphere and converting it to the vertical total electron content (VTEC) of the ionosphere.

[0040] The tilted total electron content (STEC) was extracted using a carrier phase smoothing pseudorange method.

[0041] (1)

[0042] Where, For satellite With receiver At any moment The total electron content of the ionosphere tilted by STEC, For a moment Smooth pseudoranges without geometric combination of observations The ionospheric delay conversion factor, For the first frequency, For the second frequency, For receiver First frequency With the second frequency Receiver-side differential code deviation, For satellite First frequency With the second frequency Satellite-end differential code deviation.

[0043] Since the receiver differential code bias is coupled with ionospheric delay and satellite-end differential code bias, it cannot be directly obtained. Therefore, under controllable conditions, the initial value of the receiver differential code bias will be calibrated, and the optimal processing strategy for online differential code bias calibration will be determined; simultaneously, online real-time calibration of the receiver differential code bias will be carried out. The receiver differential code bias will be extracted using a synchronization estimation method. Correcting satellite differential code bias Then, combined with the ionospheric delay conversion factor, from Get from .

[0044] estimate To mitigate the impact of ionospheric delay, nighttime observation data was used. During this period, nighttime mode was employed, ionospheric activity was calm, delay was minimal, and there were no significant gradient changes. The cutoff satellite elevation angle was set to 20°, and the average of 14 days' results was used as the receiver differential code bias. .

[0045] Using the projection function, STEC can be converted to VTEC (Vertical Total Electron Content):

[0046] (2)

[0047] in, For satellite With receiver At any moment The vertical total electron content (VTEC) of the ionosphere; This is the value of the projection function.

[0048] like Figure 2 As shown, the differential code deviation time series of various GNSS systems obtained based on this method are presented. The results show that the stability of differential code deviation estimation based on this method is better than 0.5 ns.

[0049] Step S2: Detect the ground-based BeiDou ionospheric disturbance index ROTI, including: extracting the ionospheric disturbance index ROTI of each satellite using the tilted total electron content STEC, and obtaining the ground-based BeiDou ionospheric disturbance index ROTI using the elevation angle weighting method;

[0050] The ROTI (Radio Orbital Disturbance Index) is a GNSS ionospheric disturbance monitoring index. It extracts information about ionospheric gradient changes from high-frequency GNSS observations, reflecting the presence of large-scale irregular structures in the ionosphere and characterizing the strength of ionospheric disturbance events. Calculating the ROTI involves first determining the individual satellite ROTI based on the obtained total tilt electron content (STEC) and setting a time window; then, the weights of each satellite are calculated using their elevation angles to determine the ground-based BeiDou ROTI.

[0051] The total tilt electron content (STEC) of each satellite was time-difference processed to calculate the rate of change of STEC. Its unit is TECU / min, and the calculation formula is:

[0052] (3)

[0053] Where, For satellites at certain two time points The STEC difference in the ionosphere, For time difference.

[0054] The standard deviation is calculated by sliding within a fixed time window (e.g., 5 minutes) to obtain the corresponding... Value, that is:

[0055] (4)

[0056] Where, To find the mean function.

[0057] For each visible satellite, based on its elevation angle relative to the receiver... Calculate the weighting coefficients The weighting function is as follows:

[0058] (5)

[0059] Therefore, the ionospheric disturbance index ROTI, weighted by satellite elevation angle, is:

[0060] (6)

[0061] in, The number of satellites.

[0062] Figure 3 The ROTI index sequence obtained based on this method is presented.

[0063] Step S3: Detect the ionospheric phase scintillation index This includes: using ground-based BeiDou carrier phase observation information, removing ionospheric trend term information based on Butterworth filters, and calculating the ionospheric phase scintillation index. ;

[0064] In phase detection, the standard deviation of the carrier phase is used to determine phase flicker. As shown below:

[0065] (7)

[0066] Where, The phase statistics are obtained after processing with a Butterworth filter. The carrier phase observation sequences at various GNSS frequencies are passed through a 6th-order 3dB cutoff frequency of 0.1Hz. The Butterworth high-pass filter removes low-frequency effects at the cutoff frequency, thereby obtaining high-frequency phase statistics and calculating phase flicker information.

[0067] Figure 4 The ionospheric phase scintillation index obtained based on this method is presented. sequence.

[0068] Step S4: Detect the total tropospheric zenith delay (ZTD), including: adjusting the cycle slip detection threshold in the precise single-point positioning (PPP) algorithm by combining the ionospheric disturbance index (ROTI) and the phase scintillation index. The carrier phase observations are weighted, and then the precise single-point positioning (PPP) algorithm is used to extract the tropospheric zenith total delay (ZTD).

[0069] To improve the accuracy and robustness of the tropospheric zenith total delay (ZTD) extraction process, this invention introduces ionospheric sounding data to constrain and optimize the precise point positioning (PPP) algorithm. Specifically, by extracting the ionospheric disturbance index (ROTI) corresponding to each satellite path, strong disturbance periods and paths are identified, thereby dynamically adjusting the cycle slip detection threshold in the PPP algorithm. Simultaneously, considering the impact of ionospheric scintillation on carrier observation stability, the ionospheric phase scintillation index is utilized. Weighted processing was applied to observation data from different epochs and satellites. During the PPP algorithm calculation, the weight of carrier phase observations corresponding to high scintillation indices was reduced to improve anti-interference capability and ensure the stability of parameter estimation results. Based on the above optimization strategy, and combined with auxiliary information such as precise orbit and clock bias provided by the PPP-B2b service broadcast by the BeiDou-3 GEO satellite, the PPP algorithm was used to calculate the tropospheric zenith total delay (ZTD) parameters of the stations. This process not only improved the accuracy of ZTD calculation under ionospheric disturbance environments but also provided a reliable foundation for subsequent derivation of precipitable water volume (PWV).

[0070] In the PPP algorithm, when estimating using carrier phase observations, the observation weights for each satellite are determined by the phase scintillation information detected in step S3, as follows:

[0071] (8)

[0072] Where, For satellite The carrier phase observation weights.

[0073] When performing cycle slip detection based on the geometrically uncombined GF method, the geometrically uncombined difference between adjacent epochs is compared with a fixed threshold to determine whether a cycle slip has occurred. Considering that during periods of ionospheric activity, ionospheric disturbances and scintillation may cause false detections in the fixed-threshold cycle slip detection method, mistaking carrier phase fluctuations caused by ionospheric fluctuations for cycle slips, and consequently causing the PPP algorithm to reconverge, leading to a decrease in the accuracy of ZTD parameter estimation, this method uses the ionospheric disturbance parameter ROTI from step S2 to adaptively adjust the cycle slip detection threshold. The cycle slip detection model is as follows:

[0074] (9)

[0075] Where, The threshold for GF cycle slip detection without geometric combination is given.

[0076] Furthermore, considering the anisotropic nature of zenith tropospheric delay, a horizontal gradient parameter is simultaneously introduced into the precise point positioning (PPP) algorithm to estimate tropospheric parameters. The tropospheric oblique delay of each satellite can be modeled as follows:

[0077] (10)

[0078] Where, The oblique tropospheric delay of observation station r relative to satellite s; and These are the dry and wet delay components of the zenith troposphere, respectively. and These are the projection functions for the dry and wet components, respectively; and These are the satellite elevation angle and azimuth angle, respectively. and These are the horizontal gradient parameters.

[0079] Figure 5 The tropospheric zenith total delay (ZTD) sequence obtained based on this method is presented.

[0080] Step S5: Detect tropospheric precipitable water volume (PWV), including: using temperature, humidity, and pressure observation information obtained from meteorological sensors, calculating the tropospheric dry delay (ZHD) and conversion factor, using the tropospheric zenith total delay (ZTD) to obtain the tropospheric wet delay (ZWD), and combining the conversion factor to extract the tropospheric precipitable water volume (PWV).

[0081] After estimating the tropospheric zenith total delay (ZTD) using the PPP algorithm based on GNSS precise correction products, and extracting the precipitable water volume (PWV) parameter using meteorological observation data, the following conditions are met:

[0082] (11)

[0083] in, The density of liquid water, The constant of water vapor. , Atmospheric refractive index ZHD is the atmospheric weighted average temperature; it can be calculated using the Saastanmoinen model.

[0084] The calculation formula is as follows:

[0085] (12)

[0086] in, The temperature is the surface temperature of the station, obtained from meteorological observation data.

[0087] Step S6: Detect ionospheric amplitude scintillation This includes: using signal-to-noise ratio (SNR) observation information to calculate the amplitude scintillation of the original ionosphere. Information was collected and corrected based on PWV (Precipitable Water Volume) detection data to obtain ionospheric amplitude scintillation parameters. ;

[0088] When ionospheric scintillation occurs, the amplitude of the signal undergoes rapid and random changes. In this case, the ionospheric amplitude scintillation index is commonly used. The amplitude scintillation index is used to describe the extent to which the signal amplitude is affected by ionospheric scintillation. It is typically defined as the ratio of the standard deviation of signal strength to the average signal strength over a given time period, i.e., the standard deviation of the average signal strength after normalization. In this invention, the original ionospheric amplitude scintillation... The calculation is as follows:

[0089] (13)

[0090] In the above formula Indicates the original amplitude flicker index. Signal intensity is expressed as follows:

[0091] (14)

[0092] in, Carrier-to-noise ratio; Noise density; This represents the carrier power of the signal.

[0093] Corrections were made using PWV (Precipitable Water Volume) data from the receiver to obtain ionospheric amplitude scintillation parameters. The details are as follows:

[0094] (15)

[0095] in, The set PWV threshold; This is the scaling factor.

[0096] Step S7: Output the integrated joint detection parameters of the ionosphere and troposphere using the BeiDou-based system. These parameters include: ionospheric vertical total electron content (VTEC), ionospheric disturbance index (ROTI), and ionospheric phase scintillation index. Ionospheric amplitude scintillation index Total zenith delay (ZTD) and tropospheric precipitable water volume (PWV).

[0097] A ground-based integrated joint detection device for the ionosphere and troposphere includes:

[0098] The ionospheric vertical total electron content (VTEC) detection module detects the ionospheric vertical total electron content (VTEC), including: extracting the ionospheric tilt total electron content (STEC) at the ionospheric puncture point and converting it to the ionospheric vertical total electron content (VTEC);

[0099] The ground-based BeiDou ionospheric disturbance index ROTI detection module detects the ground-based BeiDou ionospheric disturbance index ROTI, including: extracting the ionospheric disturbance index ROTI of each satellite using the tilted total electron content STEC, and obtaining the ground-based BeiDou ionospheric disturbance index ROTI using the elevation angle weighting method;

[0100] Ionospheric phase scintillation index Detection module, detecting the ionospheric phase scintillation index This includes: calculating the ionospheric phase scintillation index using ground-based BeiDou carrier phase observation information. ;

[0101] The tropospheric zenith total delay (ZTD) detection module detects the tropospheric zenith total delay (ZTD), including by combining the ionospheric disturbance index (ROTI) and the ionospheric phase scintillation index. Optimize the precise single-point positioning (PPP) algorithm and extract the tropospheric zenith total delay (ZTD).

[0102] The tropospheric precipitable water volume (PWV) detection module detects tropospheric precipitable water volume (PWV), including: calculating the tropospheric dry delay (ZHD) and conversion factor, obtaining the tropospheric wet delay (ZWD) using the tropospheric zenith total delay (ZTD), and extracting the tropospheric precipitable water volume (PWV) by combining the conversion factor.

[0103] Ionospheric amplitude scintillation index Detection module, detecting ionospheric amplitude scintillation index This includes: calculating the original ionospheric amplitude scintillation information. And utilize the precipitable water volume PWV to obtain the amplitude scintillation information of the original ionosphere. Corrections are made to obtain the ionospheric amplitude scintillation index. ;

[0104] The detection parameter output module outputs integrated ground-based BeiDou space-based ionospheric and tropospheric detection parameters, including: ionospheric vertical total electron content (VTEC), ionospheric disturbance index (ROTI), and ionospheric phase scintillation index. Ionospheric amplitude scintillation index Total zenith delay (ZTD) and tropospheric precipitable water volume (PWV).

[0105] An electronic device includes: one or more processors; and a memory for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors cause the one or more processors to implement the method.

[0106] A computer-readable storage medium having executable instructions stored thereon, which, when executed by a processor, cause the processor to implement the method described thereon.

[0107] The specific embodiments described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above are only specific embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A ground-based integrated method for joint detection of the ionosphere and troposphere, characterized in that, Includes the following steps: Step S1: Obtain the total electron content of the tilted ionosphere and convert it to the total electron content of the vertical ionosphere; Step S2: Calculate the ionospheric disturbance index of each satellite using the total tilted electron content, and further obtain the ionospheric disturbance index; Step S3: Detect the ionospheric phase scintillation index; Step S4: Extract the total tropospheric zenith delay using the ionospheric disturbance index and the ionospheric phase scintillation index; Step S5: Calculate the dry delay and conversion factor of the troposphere, obtain the wet delay of the troposphere using the total tropospheric zenith delay, and extract the precipitable water in the troposphere by combining the conversion factor. Step S6: Calculate the original ionospheric amplitude scintillation information and correct the original ionospheric amplitude scintillation information using precipitable water volume (PWV) to obtain the ionospheric amplitude scintillation index. Step S7: Output the integrated joint detection parameters of the ionosphere and troposphere based on BeiDou in space. The parameters include: total vertical electron content of the ionosphere, ionospheric disturbance index, ionospheric phase scintillation index, ionospheric amplitude scintillation index, total zenith delay of the troposphere, and tropospheric precipitable water. In step 1, the satellite differential code bias is obtained based on the differential code bias correction. Combined with the smooth pseudorange geometric combination observation and the ionospheric delay conversion factor, the total ionospheric tilt electron content (STEC) is obtained, and the vertical total ionospheric electron content (VTEC) is calculated using the projection function. In step 2, the total electron content (STEC) of ionospheric tilt for each satellite is processed by time difference, its rate of change is calculated, and the standard deviation is calculated by sliding within a fixed time window. Then, the weighting coefficient is calculated based on the satellite elevation angle, and the ground-based BeiDou ionospheric disturbance index (ROTI) is obtained by weighting the data using the weighting coefficient. In step 4, the ionospheric phase scintillation index is used. The observation weights of each satellite are determined, the cycle slip detection threshold in the single-point positioning algorithm is dynamically adjusted using the ionospheric disturbance index ROTI, and the horizontal gradient parameter is introduced to estimate the total tropospheric zenith delay (ZTD).

2. The integrated joint detection method for the ground-based space ionosphere and troposphere according to claim 1, characterized in that, In step 3, the carrier phase observation sequences of various GNSS frequencies are passed through a Butterworth high-pass filter to obtain high-frequency phase statistics, and then the exoskeleton phase scintillation index is calculated. .

3. The integrated joint detection method for the ground-based space ionosphere and troposphere according to claim 1, characterized in that, In step 5, the total zenith delay (ZTD) of the troposphere, as well as the liquid water density, water vapor gas constant, atmospheric refractive factor, and atmospheric weighted average temperature, are used to extract the tropospheric precipitable water volume (PWV).

4. The integrated joint detection method for the ground-based space ionosphere and troposphere according to claim 1, characterized in that, In step 6, the original ionospheric amplitude scintillation information is calculated by the ratio of the standard deviation of the signal strength to the average signal strength. Then, the precipitable water volume PWV was used to analyze the original ionospheric amplitude scintillation information. Corrections are made to obtain the ionospheric amplitude scintillation index. .

5. The apparatus for an integrated joint detection method of the ground-based space ionosphere and troposphere according to any one of claims 1-4, characterized in that, include: The ionospheric vertical total electron content detection module acquires the tilted total electron content of the ionosphere and converts it to the ionospheric vertical total electron content. The disturbance index detection module extracts the ionospheric disturbance index of each satellite using the total tilted electron content and obtains the ionospheric disturbance index. Phase scintillation index detection module, used to calculate the ionospheric phase scintillation index; The tropospheric zenith total delay detection module, combined with the ionospheric disturbance index and the ionospheric phase scintillation index, extracts the tropospheric zenith total delay. The tropospheric precipitable water detection module calculates the tropospheric dry delay and conversion factor, obtains the tropospheric wet delay using the total tropospheric zenith delay, and extracts the tropospheric precipitable water by combining the conversion factor. The ionospheric amplitude scintillation index detection module calculates the original ionospheric amplitude scintillation information and corrects it using precipitable water to obtain the ionospheric amplitude scintillation index. The detection parameter output module outputs ground-based BeiDou space-based integrated joint detection parameters of the ionosphere and troposphere. The parameters include: total vertical electron content of the ionosphere, ionospheric disturbance index, ionospheric phase scintillation index, ionospheric amplitude scintillation index, total zenith delay of the troposphere, and tropospheric precipitable water.

6. An electronic device, characterized in that, include: One or more processors; A memory for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors cause the one or more processors to implement the method of any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that, It stores executable instructions that, when executed by a processor, cause the processor to implement the method described in any one of claims 1 to 4.

Images