Ionosphere product data processing method, device and computer equipment

Abstract

This application relates to a method, apparatus, computer equipment, storage medium, and computer program product for processing product data of ionospheric products. The method includes: acquiring product data through ionospheric products and determining the spherical harmonic function model corresponding to the current epoch; determining the first vertical total electron content of the ionosphere corresponding to the current epoch based on observation data from a global navigation satellite system of the target area; correcting the coefficients of the spherical harmonic function model corresponding to the current epoch based on the first vertical total electron content of the ionosphere and the product data to obtain a corrected spherical harmonic function model; and determining the target ionospheric vertical total electron content based on the corrected spherical harmonic function model. This method can improve the accuracy of ionospheric product data.

Description

Technical Field

[0001] This application relates to the field of global navigation satellite system technology, and in particular to a product data processing method, apparatus, computer equipment, storage medium and computer program product for an ionospheric product. Background Technology

[0002] The ionosphere lies between approximately 60 and 1000 kilometers above the Earth's surface, and it produces effects such as reflection, scattering, absorption, and refraction on radio waves that pass through it.

[0003] In related technologies, the temporal resolution of ionospheric products is relatively small, and the broadcast interval is relatively long, resulting in a large error between the product data broadcast by the ionospheric products and the measured ionospheric data, which in turn leads to low accuracy of the product data of the ionospheric products. Summary of the Invention

[0004] Therefore, it is necessary to provide a product data processing method, apparatus, computer equipment, computer-readable storage medium, and computer program product for ionospheric products to address the above-mentioned technical problems, which can improve the accuracy of product data for ionospheric products.

[0005] Firstly, this application provides a method for processing product data of ionospheric products. The method includes:

[0006] By acquiring product data from ionospheric products, the spherical harmonic function model corresponding to the current epoch can be determined.

[0007] Based on observational data from the Global Navigation Satellite System in the target area, determine the total vertical electron content of the first ionosphere corresponding to the current epoch;

[0008] Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model.

[0009] The vertical total electron content of the target ionosphere is determined based on the modified spherical harmonic function model.

[0010] In one embodiment, the total vertical electron content of the first ionosphere corresponding to the current epoch is determined based on observation data from the Global Navigation Satellite System of the target area, including:

[0011] Acquire observation data from the Global Navigation Satellite System, and select observation data for the target area based on the observation data;

[0012] Based on observational data of the target region, the total electron content of the ionosphere is determined by non-difference and non-combination observation equations;

[0013] By projecting the total electron content of the ionosphere, the vertical total electron content of the first ionosphere corresponding to the current epoch is obtained.

[0014] In one embodiment, the total vertical electron content of the first ionosphere corresponding to the current epoch includes the total vertical electron content of the first ionosphere corresponding to multiple preset times.

[0015] Based on the total vertical electron content of the first ionosphere and the coefficients of the spherical harmonic function model corresponding to the current epoch, the corrected spherical harmonic function model is obtained, including:

[0016] For each preset time, the vertical total electron content of the first ionosphere is calculated based on the vertical total electron content of the first ionosphere at the preset time.

[0017] Based on the spherical harmonic function model corresponding to the current epoch, determine the predicted value of the vertical total electron content of the ionosphere corresponding to the puncture point value;

[0018] Based on the puncture point value and product data, determine the vertical total electron content of the second ionosphere at the preset time.

[0019] Using Kalman filtering, based on the predicted value of the vertical total electron content of the ionosphere and the vertical total electron content of the second ionosphere, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model corresponding to the preset time.

[0020] The modified spherical harmonic function model corresponding to the preset time is used as the spherical harmonic function model corresponding to the current epoch. The process of modifying the coefficients of the spherical harmonic function model corresponding to the current epoch is repeated until the modified spherical harmonic function model corresponding to the preset time satisfies the convergence condition of the Kalman filter, thus obtaining the modified spherical harmonic function model.

[0021] In one embodiment, the total vertical electron content of the second ionosphere at a preset time is determined based on the puncture point value and product data, including:

[0022] By using latitude and longitude interpolation formulas with time, and based on puncture point values ​​and product data, the vertical total electron content of the second ionosphere at a preset time is determined.

[0023] The formula for interpolating latitude, longitude, and time is:

[0024]

[0025] Where E(B,L,t) is the total vertical electron content of the second ionosphere at a preset time t, and B and L are the latitude and longitude corresponding to the puncture point, respectively. i+1 ) is T in product data i+1 The vertical total electron content of the second ionosphere at time E(B,L,T) i ) is T in product data iThe vertical total electron content of the second ionosphere at time T; i <t<T i+1 .

[0026] In one embodiment, Kalman filtering is used to correct the coefficients of the spherical harmonic function model corresponding to the current epoch based on the predicted value of the vertical total electron content of the ionosphere and the second vertical total electron content of the ionosphere, thereby obtaining the corrected spherical harmonic function model corresponding to a preset time, including:

[0027] Based on the predicted value of the vertical total electron content of the ionosphere and the vertical total electron content of the second ionosphere, the vertical total electron content of the candidate ionosphere is determined.

[0028] Based on the matrix corresponding to the Kalman filter and the total vertical electron content of the candidate ionosphere, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model.

[0029] In one embodiment, determining the candidate ionospheric vertical total electron content based on the predicted value of the ionospheric vertical total electron content and the second ionospheric vertical total electron content includes:

[0030] The difference between the vertical total electron content of the second ionosphere and the predicted value of the vertical total electron content of the ionosphere is determined to obtain the vertical total electron content of the candidate ionosphere.

[0031] Secondly, this application also provides a product data processing apparatus for ionospheric products. The apparatus includes:

[0032] By acquiring product data from ionospheric products, the spherical harmonic function model corresponding to the current epoch can be determined.

[0033] Based on observational data from the Global Navigation Satellite System in the target area, determine the total vertical electron content of the first ionosphere corresponding to the current epoch;

[0034] Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model.

[0035] The vertical total electron content of the target ionosphere is determined based on the modified spherical harmonic function model.

[0036] Thirdly, this application also provides a computer device. The computer device includes a memory and a processor, the memory storing a computer program, and the processor executing the computer program to perform the following steps:

[0037] By acquiring product data from ionospheric products, the spherical harmonic function model corresponding to the current epoch can be determined.

[0038] Based on observational data from the Global Navigation Satellite System in the target area, determine the total vertical electron content of the first ionosphere corresponding to the current epoch;

[0039] Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model.

[0040] The vertical total electron content of the target ionosphere is determined based on the modified spherical harmonic function model.

[0041] Fourthly, this application also provides a computer-readable storage medium. The computer-readable storage medium stores a computer program thereon, which, when executed by a processor, performs the following steps:

[0042] By acquiring product data from ionospheric products, the spherical harmonic function model corresponding to the current epoch can be determined.

[0043] Based on observational data from the Global Navigation Satellite System in the target area, determine the total vertical electron content of the first ionosphere corresponding to the current epoch;

[0044] Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model.

[0045] The vertical total electron content of the target ionosphere is determined based on the modified spherical harmonic function model.

[0046] Fifthly, this application also provides a computer program product. The computer program product includes a computer program that, when executed by a processor, performs the following steps:

[0047] By acquiring product data from ionospheric products, the spherical harmonic function model corresponding to the current epoch can be determined.

[0048] Based on observational data from the Global Navigation Satellite System in the target area, determine the total vertical electron content of the first ionosphere corresponding to the current epoch;

[0049] Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model.

[0050] The vertical total electron content of the target ionosphere is determined based on the modified spherical harmonic function model.

[0051] The aforementioned ionospheric product data processing method, apparatus, computer equipment, storage medium, and computer program product acquire product data from the ionospheric product and determine the spherical harmonic function model corresponding to the current epoch. Based on observation data from the Global Navigation Satellite System (GNSS) of the target area, the total vertical electron content of the first ionosphere corresponding to the current epoch is determined. Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain a corrected spherical harmonic function model. Based on the corrected spherical harmonic function model, the total vertical electron content of the target ionosphere is determined. The time accuracy of the first ionospheric vertical electron content determined by the GNSS observation data is relatively high. By combining the product data with the first ionospheric vertical electron content, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected. Through the corrected spherical harmonic function model, the total vertical electron content of the target ionosphere with high time accuracy can be obtained, thus improving the accuracy of the ionospheric product data. Attached Figure Description

[0052] Figure 1 This is an application environment diagram of a product data processing method for ionospheric products in one embodiment.

[0053] Figure 2 This is a flowchart illustrating a product data processing method for ionospheric products in one embodiment;

[0054] Figure 3 This is a flowchart illustrating a product data processing method for ionospheric products in another embodiment;

[0055] Figure 4 This is a structural block diagram of a product data processing device for an ionospheric product in one embodiment.

[0056] Figure 5 This is an internal structural diagram of a computer device in one embodiment. Detailed Implementation

[0057] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0058] The product data processing method for ionospheric products provided in this application embodiment can be applied to, for example... Figure 1 In the application environment shown, terminal 102 communicates with server 104 via a network. A data storage system can store the data that server 104 needs to process. The data storage system can be integrated onto server 104, or it can be located in the cloud or on another network server.

[0059] Terminal 102 obtains product data from server 104 through ionospheric products. Using this data, it determines the spherical harmonic function model corresponding to the current epoch. Based on observation data from the global navigation satellite system of the target area, terminal 102 determines the total vertical electron content of the first ionosphere corresponding to the current epoch. Based on this first vertical electron content, terminal 102 corrects the coefficients of the spherical harmonic function model corresponding to the current epoch, obtaining the corrected spherical harmonic function model. Finally, based on the corrected spherical harmonic function model, terminal 102 corrects the product data to obtain the target product data.

[0060] The terminal 102 can be, but is not limited to, various personal computers, laptops, smartphones, tablets, IoT devices, and portable wearable devices. IoT devices can include smart speakers, smart TVs, smart air conditioners, and smart in-vehicle systems. Portable wearable devices can include smartwatches, smart bracelets, and head-mounted devices. The server 104 can be implemented using a standalone server or a server cluster consisting of multiple servers.

[0061] In one embodiment, such as Figure 2 As shown, a product data processing method for ionospheric products is provided, which can be applied to... Figure 1 Taking the terminal in the example, the explanation includes the following steps:

[0062] Step S202: Obtain product data through ionospheric products and determine the spherical harmonic function model corresponding to the current epoch.

[0063] Among them, the ionospheric product is the global ionospheric map (GIM).

[0064] Specifically, product data can be ionospheric maps; ionospheric products obtain ionospheric maps published by UPC; UPC is an international Global Navigation Satellite System (GNSS) service that can provide real-time global ionospheric maps. The product data published by UPC is released every 15 minutes, covering a latitude range of 87.5°N to 87.5°S.

[0065] For example, the product data for the ionospheric map is shown in Table 1; the product data in Table 1 includes ionospheric maps corresponding to 97 epochs.

[0066] Table 1

[0067]

[0068] Specifically, the model of the spherical harmonic function is shown in formula (1).

[0069]

[0070] in, It is an nth-order m-th classical Legendre function, β is the geocentric latitude corresponding to the puncture point value, and a nm and b nm These are the coefficients of the spherical harmonic function model, λ s It is the diurnal longitude of the puncture point; λ s The expression formula (2) is shown.

[0071]

[0072] Where, λ S It is the diurnal longitude of the puncture point, λ IPP λ0 represents the longitude of the puncture point, λ0 represents the longitude of the line connecting the Earth's center and the Sun, and t is the number of seconds corresponding to the preset time.

[0073] Step S204: Based on the observation data of the global navigation satellite system of the target area, determine the total vertical electron content of the first ionosphere corresponding to the current epoch.

[0074] The target area can be any location represented by longitude and latitude. Global Navigation Satellite System (GNSS) observation data includes, but is not limited to: GNSS station data, clock bias files, and ephemeris files; the Vertical Total Electron Content (VTEC) is the total number of electrons per unit area column of the ionosphere perpendicular to the station. The VTEC determined using GNSS observation data is denoted as the First Ionospheric VTEC.

[0075] Specifically, observation data from the Global Navigation Satellite System (GNSS) is acquired. This data is then filtered based on the latitude and longitude of the target area to obtain the GNSS observation data for that target area. Based on this data, the total vertical electron content of the first ionosphere is determined using Precise Point Positioning (PPP).

[0076] Step S206: Based on the total vertical electron content of the first ionosphere and the product data, correct the coefficients of the spherical harmonic function model corresponding to the current epoch to obtain the corrected spherical harmonic function model.

[0077] The time accuracy of determining the total vertical electron content of the first ionosphere using observation data from the Global Navigation Satellite System (GNSS) is higher than that of the product data. For example, the time interval for the product data published by UPC is 15 minutes, but based on GNSS observation data, the total vertical electron content of the first ionosphere can be determined at a time interval of 30 seconds.

[0078] Specifically, the total vertical electron content of the first ionosphere at a preset time within the current epoch is determined by the total vertical electron content of the first ionosphere. It should be noted that when the time interval of the product data is 15 minutes, the current epoch corresponds to a 15-minute time period, and the preset time within the current epoch is a moment within that 15-minute time period. Using Kalman filtering, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected based on the total vertical electron content of the first ionosphere at the preset time and the product data, resulting in the corrected spherical harmonic function model.

[0079] Step S208: Determine the total vertical electron content of the target ionosphere based on the modified spherical harmonic function model.

[0080] Specifically, the total vertical electron content of the ionosphere at any time and any latitude and longitude can be determined by using the modified spherical harmonic function model. The target time and the target location in the target region are input into the modified spherical harmonic function model, and the target total vertical electron content of the ionosphere is determined by the modified spherical harmonic function model.

[0081] In the above-mentioned product data processing method for ionospheric products, product data is obtained through ionospheric products to determine the spherical harmonic function model corresponding to the current epoch. Based on the observation data of the Global Navigation Satellite System (GNSS) of the target area, the total vertical electron content of the first ionosphere corresponding to the current epoch is determined. Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model. Based on the corrected spherical harmonic function model, the total vertical electron content of the target ionosphere is determined. The time accuracy of the first vertical electron content of the ionosphere determined by the observation data of the GNSS is relatively high. By combining the product data with the first vertical electron content of the ionosphere, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected. Through the corrected spherical harmonic function model, the total vertical electron content of the target ionosphere with high time accuracy can be obtained, thus improving the accuracy of the product data of ionospheric products.

[0082] In some embodiments, determining the total vertical electron content of the first ionosphere at the current epoch based on observation data from the Global Navigation Satellite System (GNSS) of the target region includes: acquiring observation data from the GNSS; selecting observation data for the target region based on the observation data; determining the total electron content of the ionosphere using a non-difference non-combined observation equation based on the observation data of the target region; and projecting the total electron content of the ionosphere to obtain the total vertical electron content of the first ionosphere at the current epoch.

[0083] Specifically, the non-difference, non-combined PPP observation equation is shown in formula (3).

[0084]

[0085] Wherein, IF represents the deionization layer combination, and L... IF For the ionospheric combined carrier phase observation, P IF For the ionospheric composite code observations, ρ is the observed distance value. It is the clock bias of the receiver; It's the satellite's clock bias. It's the receiver's hardware delay. The satellite's hardware delay, T, is the tropospheric delay, ε P This is the observation noise of the pseudorange, ε L It is the observation noise of the carrier phase, N IF To eliminate ionospheric ambiguity, It is the phase delay of the receiver. It is the phase delay of the satellite.

[0086] By modifying the non-differential, non-combined PPP observation equation, a formula for calculating the initial ionospheric delay can be obtained. Based on the observation data of the global navigation satellite system in the target area and the formula for calculating the initial ionospheric delay, the initial ionospheric delay is calculated. The initial ionospheric delay is then corrected using an ionospheric delay correction formula to obtain the ionospheric delay. Based on the ionospheric delay, the total electron content of the ionosphere is calculated. Finally, a projection transformation is performed on the total electron content of the ionosphere to obtain the vertical total electron content of the first ionosphere.

[0087] The formula for ionospheric delay correction is shown in formula (4).

[0088]

[0089] in, I1 is the initial ionospheric delay, γ2 is the ratio between the square of the pseudorange frequency band and the square of the carrier frequency band, and DCB is the initial ionospheric delay. r It's a receiver bias. It's a deviation in the satellite's trajectory.

[0090] The total electron content of the ionosphere can be obtained by projecting the total electron content of the ionosphere using the projection formula, as shown in formula (5).

[0091]

[0092] Where TEC is the total electron content of the ionosphere, VTEC is the vertical total electron content of the first ionosphere, Z is the elevation angle, R is the Earth's radius, and H is the height of the puncture point above the ground.

[0093] In some embodiments, the total vertical electron content of the first ionosphere includes the total vertical electron content of the first ionosphere corresponding to multiple preset times; based on the total vertical electron content of the first ionosphere and product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model, including: for the total vertical electron content of the first ionosphere corresponding to each preset time, calculating the puncture point value corresponding to the preset time based on the total vertical electron content of the first ionosphere corresponding to the preset time; determining the predicted value of the total vertical electron content of the ionosphere corresponding to the puncture point value based on the spherical harmonic function model corresponding to the current epoch; and based on the puncture point value and product data... Data is used to determine the total vertical electron content of the second ionosphere at a preset time. Using Kalman filtering, based on the predicted total vertical electron content of the ionosphere and the total vertical electron content of the second ionosphere, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model corresponding to the preset time. This corrected spherical harmonic function model is then used as the spherical harmonic function model corresponding to the current epoch, and the process of correcting the coefficients of the spherical harmonic function model corresponding to the current epoch is repeated until the corrected spherical harmonic function model corresponding to the preset time satisfies the convergence condition of the Kalman filter, thus obtaining the corrected spherical harmonic function model.

[0094] Specifically, the multiple preset times are multiple times within the time period corresponding to the current epoch; for example, if the current epoch corresponds to the time period from 00:00 to 00:15, the multiple preset times are multiple times within that time period. In practical applications, based on the total vertical electron content of the first ionosphere corresponding to each of the multiple preset times, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected sequentially according to the order of the multiple preset times from early to late.

[0095] The following example illustrates how to correct the coefficients of the spherical harmonic function model corresponding to the current epoch based on the total vertical electron content of the first ionosphere at a preset time.

[0096] Based on the total vertical electron content of the first ionosphere at the preset time, the puncture point value corresponding to the preset time is determined. The puncture point value is estimated by the spherical harmonic function model corresponding to the current epoch, and the corresponding predicted value of the total vertical electron content of the ionosphere is obtained, as shown in formula (6).

[0097]

[0098] Among them, VTEC IPP,0 It is the predicted value of the vertical total electron content of the ionosphere. It is an nth-order m-th Legendre function, β is the geocentric latitude corresponding to the puncture point value, and a nm and b nm These are the coefficients of the spherical harmonic function model corresponding to the current epoch, λ. s It is the daily longitude of the puncture point.

[0099] Based on the puncture point values ​​and product data, interpolation is used to determine the total vertical electron content of the second ionosphere at a preset time. Kalman filtering is then applied, and based on the predicted and second ionospheric total electron content, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model at the preset time. This corrected spherical harmonic function model at the preset time is used as the spherical harmonic function model for the current epoch, and the process of correcting the coefficients of the spherical harmonic function model at the current epoch is repeated until the corrected spherical harmonic function model at the preset time satisfies the convergence condition of the Kalman filter, thus obtaining the corrected spherical harmonic function model.

[0100] For example, the coefficients of the spherical harmonic function model corresponding to the current epoch are a. nm and b nm The coefficients of the spherical harmonic function model corresponding to the current epoch are corrected based on the total vertical electron content of the first ionosphere at the first preset time. The coefficients of the corrected spherical harmonic function model are a. nm |1 and b nm |1; Based on the total vertical electron content of the first ionosphere at the second preset time, the coefficients of the modified spherical harmonic function model are corrected, and the coefficients of the modified spherical harmonic function model are a. nm |2 and b nm |2; Based on the total vertical electron content of the first ionosphere at the h-th preset time, the coefficients of the modified spherical harmonic function model are corrected, and the coefficients of the modified spherical harmonic function model are a. nm | h and b nm | h Assume the coefficient is a nm | h and b nm | h The corrected spherical harmonic function model already satisfies the convergence condition of the Kalman filter, so the convergence condition coefficient of the Kalman filter is b. nm | h and b nm | h The spherical harmonic function model is used as the modified spherical harmonic function model.

[0101] In some embodiments, the total vertical electron content of the second ionosphere at a preset time is determined based on the puncture point value and product data, including: determining the total vertical electron content of the second ionosphere at a preset time based on the puncture point value and product data using a latitude and longitude interpolation formula; the latitude and longitude interpolation formula is shown in formula (7).

[0102]

[0103] Where E(B,L,t) is the total vertical electron content of the second ionosphere at a preset time t, and B and L are the latitude and longitude corresponding to the puncture point, respectively. i+1 ) is T in product data i+1 The vertical total electron content of the second ionosphere at time E(B,L,T) i ) is T in product data i The vertical total electron content of the second ionosphere at time T; i <t<T i+1 .

[0104] In some embodiments, Kalman filtering is employed to correct the coefficients of the spherical harmonic function model corresponding to the current epoch based on the predicted value of the ionospheric vertical total electron content and the second ionospheric vertical total electron content, thereby obtaining the corrected spherical harmonic function model corresponding to a preset time. This includes: determining the candidate ionospheric vertical total electron content based on the predicted value of the ionospheric vertical total electron content and the second ionospheric vertical total electron content; and correcting the coefficients of the spherical harmonic function model corresponding to the current epoch based on the matrix corresponding to the Kalman filter and the candidate ionospheric vertical total electron content, thereby obtaining the corrected spherical harmonic function model.

[0105] Specifically, based on the predicted value of the vertical total electron content of the ionosphere and the vertical total electron content of the second ionosphere, the vertical total electron content of the candidate ionosphere is determined. The coefficients of the spherical harmonic function model corresponding to the current epoch are corrected by the matrix corresponding to the Kalman filter and the vertical total electron content of the candidate ionosphere, as shown in formula (8).

[0106]

[0107] Among them, a nm |1 and b nm |1 represents the coefficients of the modified spherical harmonic function model, a nm ,b nm K represents the coefficients of the spherical harmonic function model corresponding to the current epoch. k It is the matrix corresponding to the Kalman filter, V VTEC It is the total vertical electron content of the candidate ionosphere.

[0108] In some embodiments, determining the candidate ionospheric vertical total electron content based on the predicted value of the second ionospheric vertical total electron content includes: determining the difference between the second ionospheric vertical total electron content and the predicted value of the ionospheric vertical total electron content to obtain the candidate ionospheric vertical total electron content.

[0109] Specifically, the total vertical electron content of the candidate ionosphere is determined according to formula (9).

[0110] V VTEC =VTEC GNSS -VTEC IPP,0 (9)

[0111] Among them, V VTEC It is the total vertical electron content of the candidate ionosphere, VTEC IPP,0 It is the predicted value of the vertical total electron content of the ionosphere, VTEC. GNSS It is the total vertical electron content of the second ionosphere.

[0112] In some embodiments, the steps of modifying the coefficients of the spherical harmonic function model corresponding to the current epoch are extended to multiple epochs. For example, after obtaining the spherical harmonic function model corresponding to the current epoch, the spherical harmonic function model corresponding to the current epoch is modified based on the total vertical electron content of the first ionosphere corresponding to the next epoch to obtain the spherical harmonic function model corresponding to the next epoch.

[0113] In some embodiments, such as Figure 3 As shown, the product data processing method for ionospheric products includes:

[0114] S301, obtain product data through ionospheric products to determine the spherical harmonic function model corresponding to the current epoch;

[0115] S302: Acquire observation data from the Global Navigation Satellite System and select observation data for the target area based on the observation data;

[0116] S303, based on observational data of the target region, determines the total electron content of the ionosphere through non-difference and non-combination observational equations;

[0117] S304, Project the total electron content of the ionosphere to obtain the first vertical total electron content of the ionosphere corresponding to the current epoch;

[0118] S305, the total vertical electron content of the first ionosphere includes the total vertical electron content of the first ionosphere corresponding to multiple preset times; for each preset time corresponding to the total vertical electron content of the first ionosphere, the puncture point value corresponding to the preset time is calculated based on the total vertical electron content of the first ionosphere corresponding to the preset time.

[0119] S306, Based on the spherical harmonic function model corresponding to the current epoch, determine the predicted value of the vertical total electron content of the ionosphere corresponding to the puncture point value;

[0120] S307, based on the puncture point value and product data, determine the vertical total electron content of the second ionosphere at the preset time;

[0121] S308, determine the difference between the total vertical electron content of the second ionosphere and the predicted value of the total vertical electron content of the ionosphere, and obtain the candidate total vertical electron content of the ionosphere;

[0122] S309, based on the matrix corresponding to the Kalman filter and the total vertical electron content of the candidate ionosphere, corrects the coefficients of the spherical harmonic function model corresponding to the current epoch, and obtains the corrected spherical harmonic function model;

[0123] S310, take the corrected spherical harmonic function model corresponding to the preset time as the spherical harmonic function model corresponding to the current epoch, and repeat the above process of correcting the coefficients of the spherical harmonic function model corresponding to the current epoch until the corrected spherical harmonic function model corresponding to the preset time satisfies the convergence condition of the Kalman filter, and obtain the corrected spherical harmonic function model.

[0124] S311, the total vertical electron content of the target ionosphere is determined based on the modified spherical harmonic function model.

[0125] It should be understood that although the steps in the flowcharts of the above embodiments are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the above embodiments may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages of other steps.

[0126] In the above-mentioned product data processing method for ionospheric products, product data is obtained through ionospheric products to determine the spherical harmonic function model corresponding to the current epoch. Based on the observation data of the Global Navigation Satellite System (GNSS) of the target area, the total vertical electron content of the first ionosphere corresponding to the current epoch is determined. Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model. Based on the corrected spherical harmonic function model, the total vertical electron content of the target ionosphere is determined. The time accuracy of the first vertical electron content of the ionosphere determined by the observation data of the GNSS is relatively high. By combining the product data with the first vertical electron content of the ionosphere, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected. Through the corrected spherical harmonic function model, the total vertical electron content of the target ionosphere with high time accuracy can be obtained, thus improving the accuracy of the product data of ionospheric products.

[0127] Based on the same inventive concept, this application also provides a product data processing apparatus for ionospheric products to implement the product data processing method for ionospheric products described above. The solution provided by this apparatus is similar to the solution described in the above method. Therefore, the specific limitations in one or more embodiments of the product data processing apparatus for ionospheric products provided below can be found in the limitations of the product data processing method for ionospheric products described above, and will not be repeated here.

[0128] In one embodiment, such as Figure 4 As shown, a product data processing device for ionospheric products is provided, comprising: a spherical harmonic function model construction module, a first ionospheric vertical total electron content determination module, a model correction module, and a target ionospheric vertical total electron content determination module, wherein:

[0129] The spherical harmonic function model construction module is used to obtain product data from ionospheric products and determine the spherical harmonic function model corresponding to the current epoch.

[0130] The module for determining the total vertical electron content of the first ionosphere is used to determine the total vertical electron content of the first ionosphere at the current epoch based on observation data from the global navigation satellite system of the target area.

[0131] The model correction module is used to correct the coefficients of the spherical harmonic function model corresponding to the current epoch based on the total vertical electron content of the first ionosphere and product data, so as to obtain the corrected spherical harmonic function model.

[0132] The target ionosphere vertical total electron content determination module is used to determine the target ionosphere vertical total electron content based on the modified spherical harmonic function model.

[0133] In some embodiments, the first ionosphere vertical total electron content determination module includes:

[0134] The observation data acquisition unit is used to acquire observation data from the global navigation satellite system and select observation data for the target area based on the observation data.

[0135] The total electron content determination unit of the ionosphere is used to determine the total electron content of the ionosphere based on observation data of the target region through non-difference and non-combination observation equations.

[0136] The first ionospheric vertical total electron content determination unit is used to project the total electron content of the ionospheric layer to obtain the first ionospheric vertical total electron content corresponding to the current epoch.

[0137] In some embodiments, the total vertical electron content of the first ionosphere includes the total vertical electron content of the first ionosphere corresponding to multiple preset times; the model correction module includes:

[0138] The puncture point value determination unit is used to calculate the puncture point value corresponding to the preset time based on the total vertical electron content of the first ionosphere at each preset time.

[0139] The unit for determining the predicted value of the vertical total electron content of the ionosphere is used to determine the predicted value of the vertical total electron content of the ionosphere corresponding to the puncture point value based on the spherical harmonic function model corresponding to the current epoch.

[0140] The second ionosphere vertical total electron content determination unit is used to determine the second ionosphere vertical total electron content at a preset time based on the puncture point value and product data.

[0141] The coefficient correction unit is used to correct the coefficients of the spherical harmonic function model corresponding to the current epoch by using Kalman filtering, based on the predicted value of the vertical total electron content of the ionosphere and the vertical total electron content of the second ionosphere, so as to obtain the corrected spherical harmonic function model corresponding to the preset time.

[0142] The iterative unit is used to take the corrected spherical harmonic function model corresponding to the preset time as the spherical harmonic function model corresponding to the current epoch, and repeatedly execute the above process of correcting the coefficients of the spherical harmonic function model corresponding to the current epoch until the corrected spherical harmonic function model corresponding to the preset time satisfies the convergence condition of the Kalman filter, thus obtaining the corrected spherical harmonic function model.

[0143] In some embodiments, the second ionosphere vertical total electron content determination unit is specifically used for:

[0144] By using latitude and longitude interpolation formulas with time, and based on puncture point values ​​and product data, the vertical total electron content of the second ionosphere at a preset time is determined.

[0145] The formula for interpolating latitude, longitude, and time is:

[0146]

[0147] Where E(B,L,t) is the total vertical electron content of the second ionosphere at a preset time t, and B and L are the latitude and longitude corresponding to the puncture point, respectively. i+1 ) is T in product data i+1 The vertical total electron content of the second ionosphere at time E(B,L,T) i ) is T in product data i The vertical total electron content of the second ionosphere at time T; i <t<T i+1 .

[0148] In some embodiments, the coefficient correction unit is specifically used to determine the candidate ionospheric vertical total electron content based on the predicted value of the ionospheric vertical total electron content and the second ionospheric vertical total electron content; and to correct the coefficients of the spherical harmonic function model corresponding to the current epoch based on the matrix corresponding to the Kalman filter and the candidate ionospheric vertical total electron content, so as to obtain the corrected spherical harmonic function model.

[0149] In some embodiments, the coefficient correction unit is further specifically used to determine the difference between the second ionospheric vertical total electron content and the predicted value of the ionospheric vertical total electron content, so as to obtain the candidate ionospheric vertical total electron content.

[0150] Each module in the aforementioned product data processing device for ionospheric products can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device, or stored in the memory of a computer device as software, so that the processor can call and execute the operations corresponding to each module.

[0151] In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as follows: Figure 5As shown, the computer device includes a processor, memory, input / output interface, communication interface, display unit, and input device. The processor, memory, and input / output interface are connected via a system bus, and the communication interface, display unit, and input device are also connected to the system bus via the input / output interface. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The input / output interface is used for exchanging information between the processor and external devices. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, mobile cellular networks, NFC (Near Field Communication), or other technologies. When the computer program is executed by the processor, it implements a product data processing method for an ionospheric product. The display unit is used to form a visually visible image and can be a display screen, projection device, or virtual reality imaging device. The display screen can be an LCD screen or an e-ink screen. The input device of the computer device can be a touch layer covering the display screen, or buttons, trackballs, or touchpads set on the casing of the computer device, or external keyboards, touchpads, or mice, etc.

[0152] Those skilled in the art will understand that Figure 5 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0153] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to perform the following steps:

[0154] By acquiring product data from ionospheric products, the spherical harmonic function model corresponding to the current epoch can be determined.

[0155] Based on observational data from the Global Navigation Satellite System in the target area, determine the total vertical electron content of the first ionosphere corresponding to the current epoch;

[0156] Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model.

[0157] The vertical total electron content of the target ionosphere is determined based on the modified spherical harmonic function model.

[0158] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon, the computer program performing the following steps when executed by a processor:

[0159] By acquiring product data from ionospheric products, the spherical harmonic function model corresponding to the current epoch can be determined.

[0160] Based on observational data from the Global Navigation Satellite System in the target area, determine the total vertical electron content of the first ionosphere corresponding to the current epoch;

[0161] Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model.

[0162] The vertical total electron content of the target ionosphere is determined based on the modified spherical harmonic function model.

[0163] In one embodiment, a computer program product is provided, including a computer program that, when executed by a processor, performs the following steps:

[0164] By acquiring product data from ionospheric products, the spherical harmonic function model corresponding to the current epoch can be determined.

[0165] Based on observational data from the Global Navigation Satellite System in the target area, determine the total vertical electron content of the first ionosphere corresponding to the current epoch;

[0166] Based on the total vertical electron content of the first ionosphere and the product data, the coefficients of the spherical harmonic function model corresponding to the current epoch are corrected to obtain the corrected spherical harmonic function model.

[0167] The vertical total electron content of the target ionosphere is determined based on the modified spherical harmonic function model.

[0168] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data shall comply with the relevant laws, regulations and standards of the relevant countries and regions.

[0169] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0170] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0171] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A product data processing method of ionospheric products, characterized by, The method comprises: obtaining product data through an ionosphere product, determining a spherical harmonic function model corresponding to a current epoch; determining a first ionosphere vertical total electron content corresponding to the current epoch based on observation data of a global navigation satellite system in a target region; correcting coefficients of the spherical harmonic function model corresponding to the current epoch based on the first ionosphere vertical total electron content and the product data, to obtain a corrected spherical harmonic function model; determining a target ionosphere vertical total electron content based on the corrected spherical harmonic function model; the first ionosphere vertical total electron content corresponding to the current epoch comprises first ionosphere vertical total electron contents corresponding to a plurality of preset time points; the correcting of the coefficients of the spherical harmonic function model corresponding to the current epoch based on the first ionosphere vertical total electron content and the product data to obtain the corrected spherical harmonic function model comprises: for each first ionosphere vertical total electron content corresponding to a preset time point, calculating a piercing point value corresponding to the preset time point based on the first ionosphere vertical total electron content corresponding to the preset time point; determining an ionosphere vertical total electron content prediction value corresponding to the piercing point value based on the spherical harmonic function model corresponding to the current epoch; determining a second ionosphere vertical total electron content corresponding to the preset time point based on the piercing point value and the product data; correcting the coefficients of the spherical harmonic function model corresponding to the current epoch based on the ionosphere vertical total electron content prediction value and the second ionosphere vertical total electron content by using Kalman filtering, to obtain a corrected spherical harmonic function model corresponding to the preset time point; taking the corrected spherical harmonic function model corresponding to the preset time point as the spherical harmonic function model corresponding to the current epoch, and repeatedly performing the process of correcting the coefficients of the spherical harmonic function model corresponding to the current epoch until the corrected spherical harmonic function model corresponding to the preset time point satisfies a convergence condition of the Kalman filtering, to obtain a corrected spherical harmonic function model.

2. The method of claim 1, wherein, the determining of the first ionosphere vertical total electron content corresponding to the current epoch based on the observation data of the global navigation satellite system in the target region comprises: obtaining observation data of the global navigation satellite system, and selecting observation data of the target region based on the observation data; determining an ionosphere total electron content through a non-difference non-combined observation equation based on the observation data of the target region; projecting the ionosphere total electron content to obtain the first ionosphere vertical total electron content corresponding to the current epoch.

3. The method of claim 1, wherein, the determining of the second ionosphere vertical total electron content corresponding to the preset time point based on the piercing point value and the product data comprises: determining the second ionosphere vertical total electron content corresponding to the preset time point based on the piercing point value and the product data through a latitude-longitude-time interpolation formula; the latitude-longitude-time interpolation formula is: ； wherein, is the second ionospheric vertical total electron content at the preset time t, B and L are the latitude and longitude corresponding to the piercing point value, respectively, is the second ionospheric vertical total electron content corresponding to the time t in the product data; is the second ionospheric vertical total electron content corresponding to the time t in the product data; is the second ionospheric vertical total electron content corresponding to the time t in the product data; is the second ionospheric vertical total electron content corresponding to the time t in the product data; .

4. The method of claim 1, wherein, the correcting of the coefficients of the spherical harmonic function model corresponding to the current epoch based on the ionosphere vertical total electron content prediction value and the second ionosphere vertical total electron content by using Kalman filtering to obtain the corrected spherical harmonic function model corresponding to the preset time point comprises: determine a candidate ionospheric vertical total electron content based on the ionospheric vertical total electron content prediction value and the second ionospheric vertical total electron content; correct the coefficient of the spherical harmonic function model corresponding to the current epoch based on a matrix corresponding to Kalman filtering and the candidate ionospheric vertical total electron content, to obtain a corrected spherical harmonic function model.

5. The method of claim 4, wherein, The method comprises: determining a candidate ionospheric vertical total electron content based on the ionospheric vertical total electron content prediction value and the second ionospheric vertical total electron content.

6. A product data processing apparatus of an ionospheric product, characterized by, The method comprises: determining a candidate ionospheric vertical total electron content based on the difference between the second ionospheric vertical total electron content and the ionospheric vertical total electron content prediction value. The device comprises: a spherical harmonic function model construction module configured to acquire product data from ionospheric products and determine a spherical harmonic function model corresponding to a current epoch; a first ionospheric vertical total electron content determination module configured to determine a first ionospheric vertical total electron content corresponding to the current epoch based on observation data of a global navigation satellite system in a target region; a model correction module configured to correct the coefficient of the spherical harmonic function model corresponding to the current epoch based on the first ionospheric vertical total electron content and the product data, to obtain a corrected spherical harmonic function model; a target ionospheric vertical total electron content determination module configured to determine a target ionospheric vertical total electron content based on the corrected spherical harmonic function model; The first ionospheric vertical total electron content corresponding to the current epoch comprises a plurality of first ionospheric vertical total electron contents corresponding to a plurality of preset time points; and the model correction module comprises: a puncture point value determination unit configured to, for each first ionospheric vertical total electron content corresponding to a preset time point, calculate a puncture point value corresponding to the preset time point based on the first ionospheric vertical total electron content corresponding to the preset time point; an ionospheric vertical total electron content prediction value determination unit configured to determine an ionospheric vertical total electron content prediction value corresponding to the puncture point value based on the spherical harmonic function model corresponding to the current epoch; a second ionospheric vertical total electron content determination unit configured to determine a second ionospheric vertical total electron content corresponding to the preset time point based on the puncture point value and the product data; 7.A computer device, comprising a memory and a processor, wherein the memory stores a computer program, and the computer device is configured to perform the method according to any one of claims 1-6 when the computer program is executed by the processor. a coefficient correction unit configured to correct the coefficient of the spherical harmonic function model corresponding to the current epoch based on the ionospheric vertical total electron content prediction value and the second ionospheric vertical total electron content by using Kalman filtering, to obtain a corrected spherical harmonic function model corresponding to the preset time point; 8. A computer readable storage medium having stored thereon a computer program, characterized in that, an iteration unit configured to repeat the process of correcting the coefficient of the spherical harmonic function model corresponding to the current epoch until the corrected spherical harmonic function model corresponding to the preset time point meets a convergence condition of the Kalman filtering, to obtain a corrected spherical harmonic function model. The processor executes the computer program to implement the steps of the method in any one of claims 1 to 5. The computer program is executed by the processor to implement the steps of the method in any one of claims 1 to 5.

9. A computer program product comprising a computer program, characterized in that, The computer program, which when executed by the processor, implements the steps of the method of any one of claims 1 to 5.

Images