Method, device and use method for confirming observation error of atmospheric bending angle of occultation

Abstract

The application relates to the technical field of occultation observation, and discloses a method for confirming an observation error of an occultation atmospheric bending angle. The method comprises the following steps: obtaining an occultation atmospheric bending angle in occultation data, and first and second meteorological data sets corresponding to an occultation profile in the occultation data and not derived from a current occultation; determining a first atmospheric bending angle according to the first meteorological data set, and determining a second atmospheric bending angle according to the second meteorological data set; obtaining a first bending angle error of the occultation atmospheric bending angle and the first atmospheric bending angle, a second bending angle error of the occultation atmospheric bending angle and the second atmospheric bending angle, and a third bending angle error of the first atmospheric bending angle and the second atmospheric bending angle; and determining the observation error of the occultation atmospheric bending angle through a formula. The method can improve the objectivity of the observation error. The application also discloses a device for confirming the observation error of the occultation atmospheric bending angle and a use method.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of occultation observation, for example to a method and device for confirming an observation error of an atmospheric bending angle of occultation and a use method thereof. BACKGROUND

[0002] Global Navigation Satellite System (GNSS) occultation (RO) data has become an indispensable data source for numerical weather prediction due to its global coverage, high vertical resolution, long-term stability, and immunity to clouds and rain, and plays an important role in improving prediction results.

[0003] In the process of obtaining meteorological data, the atmospheric bending angle of occultation is an important parameter. Further, the observation error of the atmospheric bending angle is an important parameter for evaluating the reliability of the atmospheric bending angle of occultation and is a key parameter for using occultation data. There are currently various models for estimating the observation error of the atmospheric bending angle of occultation. These models are mostly divided into two categories. The first category relies on comparison with reanalysis fields, and the second category relies on the statistics of “observation-background” (O-B) differences. Among them, the first category is highly dependent on the accuracy of the reanalysis data itself, and the reliability of the data in climate-sensitive areas or areas with severe weather changes will be significantly reduced. The second category of O-B differences contains both observation errors and background field errors, which are coupled with each other and difficult to separate, resulting in significant deviation in the estimated observation error. SUMMARY

[0004] To provide a basic understanding of some aspects of the disclosed embodiments, a brief summary is given below. The summary is not an extensive review of the application, is not intended to identify key / critical elements of the application, and is not intended to delineate the scope of the application, but to present some aspects of the application in a simplified form as a prelude to the more detailed description that follows.

[0005] The embodiments of the present application provide a method and device for confirming an observation error of an atmospheric bending angle of occultation and a use method thereof, which confirm the observation error of the atmospheric bending angle of occultation in a new direction, thereby avoiding the interference of reanalysis data errors and background field errors and improving the objectivity of the observation error.

[0006] In some embodiments, the method for confirming the observation error of the atmospheric bending angle of occultation comprises:

[0007] The atmospheric bending angle in occultation data is obtained, and a first meteorological data set and a second meteorological data set corresponding to the occultation profile in the occultation data are obtained, wherein the correlation between the sources of the first meteorological data set and the second meteorological data set is lower than a first preset threshold, and the sources of the two data sets are considered to be irrelevant; the correlation between the source of the atmospheric bending angle and the first meteorological data set and the second meteorological data set is lower than a second preset threshold, so that the source of the atmospheric bending angle is irrelevant to the first meteorological data set and the second meteorological data set;

[0008] A first atmospheric bending angle is determined according to the first meteorological data set, and a second atmospheric bending angle is determined according to the second meteorological data set.

[0009] A first bending angle error between the atmospheric bending angle and the first atmospheric bending angle, a second bending angle error between the atmospheric bending angle and the second atmospheric bending angle, and a third bending angle error between the first atmospheric bending angle and the second atmospheric bending angle are obtained.

[0010] The observation error of the atmospheric bending angle in occultation is determined by the following formula:

[0011] ;

[0012] Wherein, is the observation error of the atmospheric bending angle in occultation, is the first bending angle error, is the second bending angle error, is the third bending angle error.

[0013] Optionally, the first atmospheric bending angle is determined according to the first meteorological data set, and the second atmospheric bending angle is determined according to the second meteorological data set, comprising:

[0014] The atmospheric refractive index corresponding to the first meteorological data set and the second meteorological data set is determined according to the current temperature, the current relative humidity and the current atmospheric pressure in the two data sets respectively; wherein the atmospheric refractive index is one-to-one corresponding to the altitude;

[0015] The first atmospheric bending angle and the second atmospheric bending angle corresponding to the two data sets are determined according to the atmospheric refractive index corresponding to the two data sets respectively.

[0016] Optionally, the first atmospheric bending angle and the second atmospheric bending angle corresponding to the two data sets are determined according to the atmospheric refractive index corresponding to the two data sets respectively, comprising: the refractive index profile corresponding to the two data sets is determined according to the refractive index corresponding to the two data sets respectively; based on the Abel integral equation, the first atmospheric bending angle and the second atmospheric bending angle are determined according to the refractive index profile corresponding to the two data sets and the altitude.

[0017] Optionally, the first atmospheric bending angle and the second atmospheric bending angle corresponding to the two data sets are determined according to the atmospheric refractive index corresponding to the two data sets respectively, comprising:

[0018] The height below the preset height is divided into multiple collision height layers. The collision height of each collision height layer is determined based on the reference height, the local radius of curvature, and the atmospheric refractive index and altitude corresponding to the first and second meteorological datasets, respectively. The reference height is the height of the occultation point of the occultation profile data from the reference ellipsoid, and the local radius of curvature is the local radius of curvature of the reference ellipsoid of the last data in the occultation profile data.

[0019] Based on the atmospheric refractive index of the first and second meteorological datasets corresponding to each collision height layer, determine the refractive index gradient of each collision height layer corresponding to the first and second meteorological datasets.

[0020] Based on the atmospheric refractive index, refractive index gradient, and collision height of each collision height layer corresponding to the first and second meteorological datasets, determine the atmospheric curvature angle of the first and second meteorological datasets in each collision height layer;

[0021] The first atmospheric curvature angle is determined based on multiple atmospheric curvature angles corresponding to the first meteorological dataset, and the second atmospheric curvature angle is determined based on multiple atmospheric curvature angles corresponding to the second meteorological dataset.

[0022] Optionally, the collision height of each collision height layer is determined as follows:

[0023] ;

[0024] in, The collision height layer number, For the first The collision height of each collision height layer For the first meteorological dataset or the second meteorological dataset, corresponding to the first The altitude of the collision height layer For reference height, For the local radius of curvature, For the first meteorological dataset or the second meteorological dataset, corresponding to and the The atmospheric refractive index of the collision altitude layer.

[0025] Optionally, the refractive index gradient of each collision height layer is determined based on the atmospheric refractive index of the first and second meteorological datasets corresponding to each collision height layer, including:

[0026] ;

[0027] in, The collision height layer number. a refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, a refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, a refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, a refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, a refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, a refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, a collision height of the first collision height layer, a collision height of the first collision height layer. a collision height of the first collision height layer. a collision height of the first collision height layer.

[0028] Optionally, determining the atmospheric bending angle of the first meteorological data set and the second meteorological data set in each collision height layer according to the atmospheric refractive index, the refractive index gradient and the collision height of each collision height layer corresponding to the first meteorological data set and the second meteorological data set includes:

[0029] ;

[0030] wherein, the collision height layer sequence number, the atmospheric bending angle of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, the atmospheric bending angle of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, the collision height of the first collision height layer, the collision height of the first collision height layer, the total number of collision height layers, the refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, the refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, the refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, the refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, the refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, the refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set, the refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set. the refractive index gradient of the first collision height layer corresponding to the first meteorological data set or the second meteorological data set.

[0031] Optionally, determining the atmospheric refractive index of the first meteorological data set and the second meteorological data set according to the current temperature, the current relative humidity and the current atmospheric pressure in the first meteorological data set and the second meteorological data set respectively includes:

[0032] ;

[0033] ;

[0034] ;

[0035] ;

[0036] ;

[0037] ;

[0038] wherein, is the ice surface saturation vapor pressure, is the water surface saturation vapor pressure, is the current temperature, is the lower limit temperature of pure water phase of water, is the upper limit temperature of pure ice phase of water, is and a mixing weight, is the mixed saturation vapor pressure, is the actual vapor pressure, is the current relative humidity, is the current atmospheric pressure.

[0039] Optionally, the first meteorological data set and the second meteorological data set are obtained in the following manner:

[0040] In the time dimension, the linear weighted interpolation method is used to obtain the occultation observation time based on the data of two time points before and after the occultation observation time;

[0041] In the horizontal spatial dimension, the numerical values of four grid points around the occultation tangent point are selected with the tangent point longitude and latitude coordinates as the center, and the double linear interpolation algorithm is used to calculate the position of the tangent point;

[0042] In the vertical spatial dimension, the logarithmic pressure coordinate conversion is used, and the cubic spline interpolation method is used to obtain the occultation observation height layer;

[0043] The first meteorological data and the second meteorological data corresponding to the occultation observation time, the tangent point position, and the occultation height layer are obtained at two unrelated meteorological data sources;

[0044] According to the occultation profile, a plurality of first meteorological data and second meteorological data are obtained, the plurality of first meteorological data are combined into a first meteorological data set, and the plurality of second meteorological data are combined into a second meteorological data set.

[0045] Optionally, the occulting star can receive signals of navigation satellites under different types of navigation systems; the occulting star atmospheric bending angle for obtaining occulting star data comprises: obtaining the type of navigation system corresponding to the source of the occulting star data; determining the occulting star atmospheric bending angle corresponding to the type of navigation system.

[0046] In some embodiments, the confirmation device of the occulting star atmospheric bending angle observation error comprises:

[0047] The first obtaining module is configured to obtain a bending angle of atmosphere in the occultation data, and a first meteorological data set and a second meteorological data set corresponding to the occultation profile in the occultation data, wherein the first meteorological data set and the second meteorological data set are not correlated with each other, and the correlation between the first meteorological data set and the second meteorological data set is lower than a first preset threshold; the bending angle of atmosphere is not correlated with the first meteorological data set and the second meteorological data set, and the correlation between the bending angle of atmosphere and the first meteorological data set and the second meteorological data set is lower than a second preset threshold.

[0048] The first determining module is configured to determine a first bending angle of atmosphere according to the first meteorological data set, and determine a second bending angle of atmosphere according to the second meteorological data set.

[0049] The second obtaining module is configured to obtain a first bending angle error between the bending angle of atmosphere and the first bending angle of atmosphere, a second bending angle error between the bending angle of atmosphere and the second bending angle of atmosphere, and a third bending angle error between the first bending angle of atmosphere and the second bending angle of atmosphere.

[0050] The second determining module is configured to determine the observation error of the bending angle of atmosphere according to the following formula:

[0051] ;

[0052] wherein, the observation error of the bending angle of atmosphere, the first bending angle error, the second bending angle error, the third bending angle error.

[0053] In some embodiments, the method for using the observation error of the bending angle of atmosphere includes:

[0054] obtaining a navigation system type, a region range, and a height range input by a user, and determining a mapping bending angle error according to the fitting function;

[0055] wherein, the mapping bending angle error is determined by: mapping the observation error of the bending angle of atmosphere to the mapping bending angle error according to a preset mapping mode; and the preset mapping mode corresponds to the use mode of the observation error of the bending angle of atmosphere.

[0056] The fitting function is determined in the following manner: obtaining the types of navigation systems capable of receiving the occultation received signal, dividing the earth into a plurality of regional ranges according to the earth latitude, dividing the atmospheric height into a plurality of height ranges according to the height; for each type of navigation system, each regional range, and each height range, determining the occultation atmospheric bending angle observation error by using the confirmation method provided in the foregoing embodiments; mapping each occultation atmospheric bending angle observation error to a mapped bending angle error; according to a preset fitting mode, taking the one-to-one mapped bending angle, the type of navigation system, the regional range, and the height range as the fitting data set, and fitting a fitting function taking the type of navigation system, the regional range, and the height range as the independent variable and the mapped bending angle as the dependent variable.

[0057] In some embodiments, the confirmation device of the occultation atmospheric bending angle observation error comprises a processor and a memory storing program instructions, and the processor is configured to execute the confirmation method of the occultation atmospheric bending angle observation error provided in the foregoing embodiments when executing the program instructions; or execute the use method of the occultation atmospheric bending angle observation error provided in the foregoing embodiments.

[0058] The confirmation method, device, and use method of the occultation atmospheric bending angle observation error provided in the embodiments of the present application can achieve the following technical effects:

[0059] After determining the first atmospheric bending angle according to the first meteorological data set and determining the second atmospheric bending angle according to the second meteorological data set, the error of the first atmospheric bending angle is caused by the source of the first meteorological data set, and the error of the second atmospheric bending angle is caused by the source of the second meteorological data set.

[0060] Since the sources of the first meteorological data set and the second meteorological data set are considered to be irrelevant, the errors of the first atmospheric bending angle and the second atmospheric bending angle are also considered to be irrelevant and independent of each other. At the same time, the occultation atmospheric bending angle also has an error, and the source of the occultation atmospheric bending angle is irrelevant to the first meteorological data set and the second meteorological data set, so the error of the occultation atmospheric bending angle is also irrelevant to the errors of the first atmospheric bending angle and the second atmospheric bending angle, and is independent of each other.

[0061] Therefore, in the process of obtaining the first bending angle error of the occultation atmospheric bending angle and the first atmospheric bending angle, the second bending angle error of the occultation atmospheric bending angle and the second atmospheric bending angle, and the third bending angle error of the first atmospheric bending angle and the second atmospheric bending angle, the part of the error will not be "offset" because the error has a certain correlation, so as to indirectly obtain the occultation atmospheric bending angle observation error through the three bending angle errors.

[0062] On this basis, the observation error of the occultation atmospheric bending angle is calculated by the formula . .

[0063] The confirmation method is not dependent on comparison with reanalysis fields and statistics of "observation-background" (O-B) differences, and avoids interference of reanalysis data errors and background field errors. The confirmation method indirectly evaluates and confirms the observation error of the atmospheric bending angle by taking two meteorological data sets that can be considered as irrelevant as references, and realizes objective and pure estimation of the observation error of the atmospheric bending angle.

[0064] The foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the application. BRIEF DESCRIPTION OF DRAWINGS

[0065] One or more embodiments are illustrated by way of example in the figures that are not intended to be limiting of the application as defined by the claims. The same numbers in different figures identify the same components or features. In the figures:

[0066] Figure 1 is a flowchart of a confirmation method for an observation error of an atmospheric bending angle provided by an embodiment of the application;

[0067] Figure 2 is a flowchart of a method for calculating a first atmospheric bending angle and a second atmospheric bending angle provided by an embodiment of the application;

[0068] Figure 3 is a schematic diagram of a confirmation device for an observation error of an atmospheric bending angle provided by an embodiment of the application;

[0069] Figure 4 is a schematic diagram of another confirmation device for an observation error of an atmospheric bending angle provided by an embodiment of the application. DETAILED DESCRIPTION

[0070] In order to enable a more detailed understanding of the features and technical content of the embodiments of the application, the implementation of the embodiments of the application is described in detail below with reference to the accompanying drawings, which are only used for reference and do not limit the embodiments of the application. In the following technical description, for the purpose of easy explanation, a plurality of details are provided to provide a full understanding of the disclosed embodiments. However, one or more embodiments can still be implemented without these details. In other cases, well-known structures and devices can be simplified to facilitate the drawings.

[0071] The terms "first", "second", etc. in the description and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and do not necessarily indicate a specific order or sequence. It should be understood that the data thus used can be interchanged under appropriate circumstances, so that the embodiments of the present application described herein can be implemented. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion.

[0072] Unless otherwise specified, the term "a plurality of" means two or more.

[0073] In the embodiments of the present application, the character " / " represents a "or" relationship between the objects before and after it. For example, A / B means: A or B.

[0074] The term "and / or" is a description of the association relationship between objects, which means that there can be three relationships. For example, A and / or B means: A or B, or, A and B, three relationships.

[0075] Figure 1 is a flowchart of a method for confirming the observation error of the atmospheric bending angle of occultation provided by the embodiments of the present application.

[0076] In combination Figure 1 As shown, the method for confirming the observation error of the atmospheric bending angle of occultation comprises:

[0077] S101, obtaining the atmospheric bending angle of occultation in the occultation data, and the first meteorological data set and the second meteorological data set corresponding to the occultation profile in the occultation data, which are not derived from the current occultation.

[0078] The term "occultation profile" refers to a series of continuous parameters that vary with height, including but not limited to: atmospheric bending angle, atmospheric refractive index, air temperature, air humidity, air density, etc.

[0079] In a direct understanding, the occultation profile can be regarded as a line approximately perpendicular to the ground, and then according to the different heights on this line, two sets of meteorological data are obtained from other two sources, one of which is the first meteorological data set, and the other is the second meteorological data set.

[0080] If the correlation between the sources of the first meteorological data set and the second meteorological data set is lower than the first preset threshold, the sources of the two data sets are considered to be irrelevant.

[0081] If the correlation between the sources of the first meteorological data set and the second meteorological data set is lower than the first preset threshold, the sources of the two data sets are considered to be irrelevant; otherwise, the sources of the two data sets are considered to be irrelevant.

[0082] The similarity between the first and second meteorological data sets can be determined by analyzing the assimilated data and assimilation methods of the first and second meteorological data sets.

[0083] For example, the Global Forecast System (GFS) uses the Global Data Assimilation System (GDAS) as initialization, and the GDAS adds the following types of observation data to the gridded three-dimensional model space: ground observation data, sounding balloon data, wind profile data, aircraft reports, buoy observations, radar observations, and satellite observations.

[0084] In contrast, the EAR5 (ECMWF Reanalysis v5) is the fifth generation of atmospheric reanalysis data sets of global climate from January 1950 to the present by the European Centre for Medium-Range Weather Forecasts (ECMWF), which is generated based on the 4D-Var data assimilation and model prediction of the ECMWF Integrated Forecasting System (IFS), and the assimilated data is mostly from satellites, but also from some ground and aircraft reports.

[0085] That is, the generation methods and assimilated data of the GFS and the EAR5 are different, and the GFS and the EAR5 can be considered as unrelated, and the GFS and the EAR5 are respectively taken as the first meteorological data set and the second meteorological data set, which satisfies the condition that the correlation between the sources of the first and second meteorological data sets is lower than the first preset threshold.

[0086] The correlation between the source of the atmospheric bending angle of the occultation and the first and second meteorological data sets is lower than the second preset threshold, and the source of the atmospheric bending angle is unrelated to the first and second meteorological data sets.

[0087] If the correlation between the source of the atmospheric bending angle of the occultation and the first and second meteorological data sets is lower than the second preset threshold, it means that the source of the atmospheric bending angle can be considered as unrelated to the first and second meteorological data sets; otherwise, it cannot be considered as unrelated.

[0088] The term "atmospheric bending angle" refers to the angle of deflection of an electromagnetic wave when it passes through the atmosphere, and the atmospheric bending angle of the occultation can be directly calculated according to the received electromagnetic wave signal.

[0089] In the case that the satellite sending the electromagnetic wave signal is different from the satellite involved in the GFS and the EAR5, the occultation atmospheric bending angle calculated directly according to the received electromagnetic wave signal is irrelevant to the GFS and the EAR5, and the relevance of the source of the occultation atmospheric bending angle to the first meteorological data set and the second meteorological data set is of course lower than the second preset threshold.

[0090] If there is a certain overlap between the satellite sending the electromagnetic wave signal of the occultation atmospheric bending angle and the source of the first meteorological data set or the second meteorological data set, there is a certain relevance between the source of the occultation atmospheric bending angle and the first meteorological data set and the second meteorological data set; or, the occultation atmospheric bending angle has been optimized based on other meteorological data, then it is also necessary to judge whether the other meteorological data is relevant to the first meteorological data and the second meteorological data, if relevant, the source of the occultation atmospheric bending angle is relevant to the first meteorological data set and the second meteorological data set.

[0091] In this case, the relevance needs to be carefully evaluated, and it is judged whether the source of the occultation atmospheric bending angle is considered to be irrelevant to the first meteorological data set and the second meteorological data set.

[0092] S102, determining a first atmospheric bending angle according to the first meteorological data set and a second atmospheric bending angle according to the second meteorological data set.

[0093] The first atmospheric bending angle can be obtained by forward modeling from the first meteorological data set, and the second atmospheric bending angle can be obtained by forward modeling from the second meteorological data set.

[0094] For example, determining the first atmospheric bending angle according to the first meteorological data set and the second atmospheric bending angle according to the second meteorological data set can include: determining the atmospheric refractive index corresponding to the first meteorological data set and the second meteorological data set respectively according to the current temperature, the current relative humidity and the current atmospheric pressure in the two meteorological data sets; wherein the atmospheric refractive index corresponds to the altitude one-to-one; determining the first atmospheric bending angle and the second atmospheric bending angle corresponding to the two meteorological data sets respectively according to the atmospheric refractive index corresponding to the two meteorological data sets.

[0095] Wherein, determining the first atmospheric bending angle and the second atmospheric bending angle corresponding to the two meteorological data sets respectively according to the atmospheric refractive index corresponding to the two meteorological data sets can include: determining the refractive index profile corresponding to the first meteorological data set and the second meteorological data set respectively according to the atmospheric refractive index corresponding to the two meteorological data sets; determining the first atmospheric bending angle and the second atmospheric bending angle according to the refractive index profile corresponding to the two meteorological data sets and the altitude based on the Abel integral equation.

[0096] The Abel integral equation is specifically:

[0097] ;

[0098] Wherein, is the atmospheric bending angle, is the collision height, is the atmospheric refractive index, is the distance between the refractive point and the center of the earth.

[0099] S103, obtaining a first bending angle error of the occultation atmospheric bending angle and the first atmospheric bending angle, a second bending angle error of the occultation atmospheric bending angle and the second atmospheric bending angle, and a third bending angle error of the first atmospheric bending angle and the second atmospheric bending angle.

[0100] S104, determining the observation error of the occultation atmospheric bending angle by the following formula:

[0101]

[0102] wherein, is the observation error of the occultation atmospheric bending angle, is the first bending angle error, is the second bending angle error, is the third bending angle error.

[0103] In this embodiment, after the first atmospheric bending angle is determined according to the first meteorological data set and the second atmospheric bending angle is determined according to the second meteorological data set, the error of the first atmospheric bending angle is caused by the source of the first meteorological data set, and the error of the second atmospheric bending angle is caused by the source of the second meteorological data set.

[0104] Since the sources of the first meteorological data set and the second meteorological data set are considered to be irrelevant, the error of the first atmospheric bending angle and the error of the second atmospheric bending angle are also considered to be irrelevant and independent of each other. At the same time, the occultation atmospheric bending angle also has an error, and the source of the occultation atmospheric bending angle is irrelevant to the first meteorological data set and the second meteorological data set, so the error of the occultation atmospheric bending angle is also irrelevant to the error of the first atmospheric bending angle and the error of the second atmospheric bending angle, and independent of each other.

[0105] Therefore, in the process of obtaining the first bending angle error of the occultation atmospheric bending angle and the first atmospheric bending angle, the second bending angle error of the occultation atmospheric bending angle and the second atmospheric bending angle, and the third bending angle error of the first atmospheric bending angle and the second atmospheric bending angle, the error will not be "offset" due to a certain degree of correlation, so as to indirectly obtain the observation error of the occultation atmospheric bending angle through the three bending angle errors.

[0106] On this basis, the observation error of the occultation atmospheric bending angle is calculated by the formula .

[0107] ​​​The confirmation method does not rely on comparison with reanalysis fields and statistics of "observation-background" (O-B) differences, avoids interference of reanalysis data errors and background field errors, and indirectly evaluates and confirms the observation error of the atmospheric bending angle of occultation stars by taking two meteorological data sets that can be considered as irrelevant as references, thereby achieving objective and pure estimation of the observation error of the atmospheric bending angle of occultation stars.

[0108] The following exemplary describes a way of determining the atmospheric refractive index corresponding to the first meteorological data set and the second meteorological data set according to the current temperature, the current relative humidity and the current atmospheric pressure in the two meteorological data sets respectively.

[0109] ;

[0110] ;

[0111] ;

[0112] ;

[0113] ;

[0114] ;

[0115] wherein, is the ice surface saturated water vapor pressure, is the water surface saturated water vapor pressure, is the current temperature, is the lower limit temperature of pure water phase of water, and the value is 273.16K, is the upper limit temperature of pure ice phase of water, and the value is 250.16K, is the mixing weight of and is the mixed saturated water vapor pressure, is the actual water vapor pressure, is the current relative humidity, is the current atmospheric pressure.

[0116] Regarding the process of determining the first atmospheric bending angle and the second atmospheric bending angle corresponding to the first meteorological data set and the second meteorological data set according to the atmospheric refractive index corresponding to the two meteorological data sets respectively, the foregoing embodiment can be calculated by the Abel integral equation, and in addition, the embodiment of the application also exemplarily provides a specific calculation method.

[0117] Figure 2 is the flowchart provided by the embodiment of the application for calculating the first atmospheric bending angle and the second atmospheric bending angle.

[0118] In combination with​Figure 2 As shown, based on the atmospheric refractive indices of the two datasets (the first meteorological dataset and the second meteorological dataset), the corresponding first atmospheric curvature angle and second atmospheric curvature angle are determined, including:

[0119] S201. Divide the height below the preset height into multiple collision height layers.

[0120] The preset altitude can be [60, 100] km, for example, the preset altitude can be 60 km, 70 km, 80 km, 90 km or 100 km.

[0121] The total number of collision height layers can be determined based on specific needs. This can be achieved through the formula... In the process of calculating the observation error of the atmospheric curvature angle of an occultation, if the total number of layers is set too much, that is, too many samples, it is easy to have an "overfitting" phenomenon, which reduces the scope of applicability of the assessed and confirmed observation error of the atmospheric curvature angle of an occultation. If the total number of layers is set too little, that is, too few samples, it is not enough to assess and confirm a relatively accurate observation error of the atmospheric curvature angle of an occultation.

[0122] The total number of collision height layers is consistent with the aforementioned calculated atmospheric refractive index data. In other words, the total number of collision height layers can be determined based on the calculated typical atmospheric refractive index. Typical atmospheric refractive index refers to the atmospheric refractive index at the inflection point of the refractive index profile calculated from a meteorological dataset, and / or, typical atmospheric refractive index refers to the atmospheric refractive index where the difference between two meteorological datasets is large enough to satisfy a preset discrimination condition. The smaller the preset discrimination condition, the more prone it is to overfitting; the larger the preset discrimination condition, the more prone it is to insufficient data. Those skilled in the art determine the preset discrimination condition based on this and the required conditions.

[0123] Each collision height layer corresponds to a collision height. Specifically, the collision height of each collision height layer is determined based on the reference height, the local radius of curvature, and the atmospheric refractive index and altitude corresponding to the first and second meteorological datasets, respectively. The reference height is the height of the occultation point of the occultation profile data from the reference ellipsoid, and the local radius of curvature is the local radius of curvature of the reference ellipsoid of the last data in the occultation profile data.

[0124] S202. Based on the atmospheric refractive index of the first meteorological dataset and the second meteorological dataset corresponding to each collision height layer, determine the refractive index gradient of each collision height layer corresponding to the first meteorological dataset and the second meteorological dataset.

[0125] Specifically, the refractive index gradient of each collision height layer can be determined using the following formula:

[0126] ;

[0127] wherein, is a collision height layer number, is a refractive index gradient of the first or second meteorological data set corresponding to the th collision height layer, is an atmospheric refractive index of the first or second meteorological data set corresponding to the th collision height layer, is an atmospheric refractive index corresponding to the th collision height layer, is a collision height of the th collision height layer, is a collision height of the th collision height layer.

[0128] S203, determining the atmospheric bending angles of the first and second meteorological data sets in each collision height layer according to the atmospheric refractive indices, the refractive index gradients, and the collision heights of each collision height layer corresponding to the first and second meteorological data sets.

[0129] S204, determining a first atmospheric bending angle according to the plurality of atmospheric bending angles corresponding to the first meteorological data set, and determining a second atmospheric bending angle according to the plurality of atmospheric bending angles corresponding to the second meteorological data set.

[0130] The first and second atmospheric bending angles determined by the embodiments of the present application are used to calculate the observation error of the atmospheric bending angle of the occultation according to the formula In this process, it is not necessary to calculate the first and second atmospheric bending angles of each height point on the occultation profile based on the original Abbe. The first and second atmospheric bending angles are calculated in the layered manner, which is not to calculate the first and second atmospheric bending angles of each height point on the occultation profile in each layer, but to calculate a set of first and second atmospheric bending angles corresponding to each collision height layer, and then to calculate the observation error of the atmospheric bending angle of the occultation according to the formula which meets the expected requirements in terms of accuracy and expansibility.

[0131] After the above-mentioned embodiments give the layered calculation idea, those skilled in the art can design specific calculation formulas for each step according to experience.

[0132] In the following, the embodiments of the present application will exemplarily give specific calculation formulas for each step, but do not constitute specific limitations on the specific calculation formulas of the specific steps, for example, other formulas can be obtained from other ideas, or formulas with some differences in specific parameters can be used.

[0133] The following calculation formula given in the embodiments of the present application is exemplarily described:

[0134] Optionally, the collision height of each collision height layer is determined by the following method:

[0135] ;

[0136] wherein, is the collision height layer number, is the collision height of the th collision height layer, is the altitude in the first meteorological data set or the second meteorological data set corresponding to the th collision height layer, is the reference height, is the local radius of curvature, is the atmospheric refractive index in the first meteorological data set or the second meteorological data set corresponding to the th and the th collision height layer.

[0137] Optionally, according to the atmospheric refractive index, the refractive index gradient and the collision height of each collision height layer corresponding to the first meteorological data set and the second meteorological data set, the atmospheric bending angle of each collision height layer corresponding to the first meteorological data set and the second meteorological data set is determined, comprising:

[0138]

[0139] wherein, is the collision height layer number, is the atmospheric bending angle of the th collision height layer corresponding to the first meteorological data set or the second meteorological data set, is the collision height of the th collision height layer, is the total number of collision height layers, is the refractive index gradient of the th collision height layer corresponding to the first meteorological data set or the second meteorological data set, is the collision height of the th collision height layer, is the atmospheric refractive index of the th collision height layer corresponding to the first meteorological data set or the second meteorological data set, is the collision height of the th collision height layer.

[0140] The above formulae can be used in combination or independently, for example, the calculation formula given in the above embodiment can be used in both of the above two steps; or the calculation formula of the above embodiment is used when calculating the collision height of the collision height layer, but other calculation formula is used when calculating the first and second atmospheric bending angles corresponding to each height layer; or other calculation formula is used when calculating the collision height of the collision height layer, but the calculation formula of the above embodiment is used when calculating the first and second atmospheric bending angles corresponding to each height layer.

[0141] The way of obtaining the first and second meteorological data sets will be further exemplarily described below.

[0142] Optionally, the way of obtaining the first and second meteorological data sets comprises:

[0143] In the time dimension, the linear weighted interpolation method is used to obtain the occultation observation time based on the data of two time points before and after the occultation observation time;

[0144] In the horizontal spatial dimension, the numerical values of four grid points around the occultation tangent point are selected, and the tangent point position is calculated by the bilinear interpolation algorithm with the tangent point longitude and latitude coordinates as the center;

[0145] In the vertical spatial dimension, the occultation height layer is obtained by the logarithmic pressure coordinate conversion and the cubic spline interpolation method;

[0146] The first and second meteorological data corresponding to the occultation observation time, the tangent point position and the occultation height layer are obtained at two unrelated meteorological data sources respectively;

[0147] According to the occultation profile, a plurality of first and second meteorological data are obtained, a plurality of first meteorological data are combined into a first meteorological data set, and a plurality of second meteorological data are combined into a second meteorological data set.

[0148] After matching in time, space (vertical and horizontal), it can be ensured that the occultation atmospheric bending angle, the first meteorological data set and the second meteorological data set are matched in time and space.

[0149] Finally, in a specific application scenario, the occultation can receive signals of navigation satellites under different types of navigation systems, in which case, the occultation atmospheric bending angle obtained from the occultation data comprises: obtaining the type of navigation system corresponding to the source of the occultation data; determining the occultation atmospheric bending angle corresponding to the type of navigation system.

[0150] The navigation system related by the embodiments of the present application includes, but is not limited to, Beidou Navigation Satellite System (BDS or COMPASS), Global Positioning System (GPS), Global Navigation Satellite System (GLONASS), Galileo Satellite Navigation system (GALILEO).

[0151] The embodiments of the present application also provide a use method of occultation atmospheric bending angle observation error, which comprises: obtaining a navigation system type, a region range and a height range input by a user, and determining a mapping bending angle error according to a fitting function.

[0152] The mapping bending angle error is determined by: mapping the occultation atmospheric bending angle observation error into the mapping bending angle error according to a preset mapping mode; and the preset mapping mode corresponds to a use mode of the occultation atmospheric bending angle observation error.

[0153] The fitting function is determined by: obtaining a navigation system type capable of receiving a signal of occultation, dividing the earth into a plurality of region ranges according to the latitude of the earth, and dividing the atmospheric height into a plurality of height ranges according to the height; for each navigation system type, each region range and each height range, the confirmation method of the occultation atmospheric bending angle observation error provided in the foregoing embodiments is used to determine the occultation atmospheric bending angle observation error; each occultation atmospheric bending angle observation error is mapped into a mapping bending angle error; and according to a preset fitting mode, the mapping bending angle, the navigation system type, the region range and the height range are taken as a fitting data set, and a fitting function is fitted, in which the navigation system type, the region range and the height range are taken as independent variables, and the mapping bending angle is taken as a dependent variable.

[0154] In some application scenarios, the region range can be divided into two ranges according to -40° and 40°, and the height range can be 0-4km, 4-18km, 18-60km, or 0-8km, 8-12km, 12-60km, the height range can include or not include the range endpoints, and the height ranges are not overlapped and represent the overall height completely.

[0155] The preset mapping mode can include:

[0156]

[0157] The mapping bending angle error is determined by: mapping the occultation atmospheric bending angle observation error into the mapping bending angle error according to a preset mapping mode; and the preset mapping mode corresponds to a use mode of the occultation atmospheric bending angle observation error. to map the bending angle error. to map the bending angle error.

[0158] The fitting manner can include:

[0159]

[0160] or,

[0161]

[0162] wherein, , , , are fitting parameters, is a kind of navigation system, area range, height range.

[0163] In some embodiments, the fitting function obtained through the above mapping and the above fitting manner is assimilated into a Gridpoint Statistical Interpolation (GSI) assimilation system.

[0164] Figure 3 is a schematic diagram of a device for confirming an occultation atmospheric bending angle observation error provided by an embodiment of the present application. The device can be realized in the form of software, hardware or a combination of both.

[0165] As shown in Figure 3 , the device for confirming an occultation atmospheric bending angle observation error includes a first obtaining module 31, a first determining module 32, a second obtaining module 33 and a second determining module 34.

[0166] The first obtaining module 31 is configured to obtain an occultation atmospheric bending angle in occultation data and first and second meteorological data sets corresponding to an occultation profile in the occultation data, the first and second meteorological data sets not being derived from a current occultation. The correlation between the sources of the first and second meteorological data sets is lower than a first preset threshold, and the sources of the two data sets are considered to be irrelevant. The correlation between the source of the occultation atmospheric bending angle and the first and second meteorological data sets is lower than a second preset threshold, and the source of the atmospheric bending angle is irrelevant to the first and second meteorological data sets.

[0167] The first determining module 32 is configured to determine a first atmospheric bending angle according to the first meteorological data set and a second atmospheric bending angle according to the second meteorological data set.

[0168] The second obtaining module 33 is configured to obtain a first bending angle error of the occultation atmospheric bending angle and the first atmospheric bending angle, a second bending angle error of the occultation atmospheric bending angle and the second atmospheric bending angle, and a third bending angle error of the first atmospheric bending angle and the second atmospheric bending angle.

[0169] The second determining module 34 is configured to determine the occultation atmospheric bending angle observation error according to the following formula.

[0170] ;

[0171] Wherein, the occultation atmospheric bending angle observation error, the first bending angle error, the second bending angle error, the third bending angle error.

[0172] Optionally, the first determining module 32 includes a first determining unit and a second determining unit.

[0173] Wherein, the first determining unit is configured to determine the atmospheric refractive index corresponding to the first meteorological data set and the second meteorological data set according to the current temperature, the current relative humidity and the current atmospheric pressure in the two meteorological data sets respectively; wherein the atmospheric refractive index corresponds to the altitude one-to-one; the second determining unit is configured to determine the first atmospheric bending angle and the second atmospheric bending angle corresponding to the two atmospheric refractive indexes respectively according to the two atmospheric refractive indexes.

[0174] Optionally, the second determining unit is specifically configured to determine the refractive index profile corresponding to the first meteorological data set and the second meteorological data set according to the two atmospheric refractive indexes respectively; and determine the first atmospheric bending angle and the second atmospheric bending angle according to the refractive index profile corresponding to the two meteorological data sets and the altitude based on the Abel integral equation.

[0175] Optionally, the second determining unit is specifically configured to:

[0176] divide the height below the preset height into a plurality of collision height layers; wherein the collision height of each collision height layer is determined according to the reference height, the local radius of curvature, and the atmospheric refractive index and the altitude corresponding to the first meteorological data set and the second meteorological data set respectively, the reference height is the height of the occultation point of the occultation profile data from the reference ellipsoid, and the local radius of curvature is the local radius of curvature of the last data reference ellipsoid in the occultation profile data;

[0177] determine the refractive index gradient of each collision height layer corresponding to the first meteorological data set and the second meteorological data set according to the atmospheric refractive index of the first meteorological data set and the second meteorological data set corresponding to each collision height layer;

[0178] determining, according to the atmospheric refractive index, the refractive index gradient and the collision height of each collision height layer corresponding to the first meteorological data set and the second meteorological data set, the atmospheric bending angle of each collision height layer corresponding to the first meteorological data set and the second meteorological data set;

[0179] determining the first atmospheric bending angle according to the plurality of atmospheric bending angles corresponding to the first meteorological data set, and determining the second atmospheric bending angle according to the plurality of atmospheric bending angles corresponding to the second meteorological data set.

[0180] Optionally, the collision height of each collision height layer is determined by:

[0181] ;

[0182] wherein, is the collision height layer number, is the collision height of the i-th collision height layer, is the altitude in the first meteorological data set or the second meteorological data set corresponding to the i-th collision height layer, is the reference height, is the local radius of curvature, is the atmospheric refractive index in the first meteorological data set or the second meteorological data set corresponding to the i-th collision height layer and the j-th collision height layer. Optionally, the refractive index gradient of each collision height layer is determined according to the atmospheric refractive index of the first meteorological data set and the second meteorological data set corresponding to each collision height layer, comprising: ;

[0183] wherein,

[0184] is the collision height layer number, is the refractive index gradient of the i-th collision height layer corresponding to the first meteorological data set or the second meteorological data set, is the atmospheric refractive index of the i-th collision height layer corresponding to the first meteorological data set or the second meteorological data set,

[0185] is the atmospheric refractive index corresponding to the i-th collision height layer, is the collision height of the i-th collision height layer, is the collision height of the i-th collision height layer.

[0186] ​​​​​​​​​Optionally, the atmospheric bending angles of each collision height layer corresponding to the first meteorological data set and the second meteorological data set are determined according to the atmospheric refractive index, the refractive index gradient and the collision height of each collision height layer corresponding to the first meteorological data set and the second meteorological data set, and the determination comprises:

[0187] ;

[0188] wherein, is the collision height layer number, is the atmospheric bending angle of the first meteorological data set or the second meteorological data set corresponding to the i th collision height layer, is the collision height of the i th collision height layer, is the total number of collision height layers, is the refractive index gradient of the first meteorological data set or the second meteorological data set corresponding to the i th collision height layer, is the collision height of the i th collision height layer, is the atmospheric refractive index of the first meteorological data set or the second meteorological data set corresponding to the i th collision height layer, is the collision height of the i th collision height layer. Optionally, the first determination unit is specifically configured to: ;

[0189] ;

[0190] ;

[0191] ;

[0192] ;

[0193] ;

[0194] ;

[0195] ;

[0196] wherein, is the ice surface saturated water vapor pressure, is the water surface saturated water vapor pressure, is the current temperature, is the pure water phase lower limit temperature of water, is the pure ice phase upper limit temperature of water, is the mixing weight of and , and is the mixed saturated water vapor pressure. ​​​​is a current relative humidity, is a current atmospheric pressure. is a current atmospheric pressure.

[0197] Optionally, the first meteorological data set and the second meteorological data set are obtained in the following manner: in the time dimension, the data at two time points before and after the occultation observation time is used to obtain the occultation observation time by linear weighted interpolation; in the horizontal spatial dimension, the values of four grid points around the occultation tangent point are selected, and the tangent point position is calculated by bilinear interpolation; in the vertical spatial dimension, the first meteorological data and the second meteorological data corresponding to the occultation observation height layer are obtained by logarithmic pressure coordinate conversion and cubic spline interpolation; the first meteorological data and the second meteorological data corresponding to the occultation observation time, the tangent point position, and the occultation height layer are obtained from two unrelated meteorological data sources; a plurality of first meteorological data and a plurality of second meteorological data are obtained according to the occultation profile, the plurality of first meteorological data are combined into the first meteorological data set, and the plurality of second meteorological data are combined into the second meteorological data set.

[0198] Optionally, the occultation satellite can receive signals of navigation satellites of different types of navigation systems.

[0199] Optionally, the first obtaining module 31 includes an obtaining unit and a third determining unit; the obtaining unit is configured to obtain the type of the navigation system corresponding to the occultation data source; and the third determining unit is configured to determine the occultation atmospheric bending angle corresponding to the type of the navigation system.

[0200] In some embodiments, the present application provides a device for confirming an occultation atmospheric bending angle observation error, which includes a processor and a memory storing program instructions, and the processor is configured to execute the method for confirming an occultation atmospheric bending angle observation error provided in the foregoing embodiments when executing the program instructions.

[0201] In this case, the device is configured to confirm the occultation atmospheric bending angle observation error.

[0202] In some embodiments, the present application provides a device for confirming an occultation atmospheric bending angle observation error, which includes a processor and a memory storing program instructions, and the processor is configured to execute the method for confirming an occultation atmospheric bending angle observation error provided in the foregoing embodiments when executing the program instructions.

[0203] In this case, the device is configured to confirm the occultation atmospheric bending angle observation error.

[0204] Figure 4 is a schematic diagram of a device for confirming an occultation atmospheric bending angle observation error provided by an embodiment of the present application, which is combined with Figure 4 As shown in the figure, the device for confirming an occultation atmospheric bending angle observation error includes:

[0205] The processor 41 and the memory 42 can also include a communication interface 43 and a bus 44. The processor 41, the communication interface 43, and the memory 42 can communicate with each other through the bus 44. The communication interface 43 can be used for information transmission. The processor 41 can invoke the logic instructions in the memory 42 to execute the method for confirming the atmospheric bending angle observation error of occultation stars or the use method provided in the foregoing embodiments.

[0206] In addition, when the logic instructions in the memory 42 described above are implemented in the form of a software functional unit and sold or used as an independent product, the logic instructions can be stored in a computer-readable storage medium.

[0207] The memory 42, as a computer-readable storage medium, can be used to store software programs, computer-executable programs, such as program instructions / modules corresponding to the method in the embodiments of the present application. The processor 41 executes the functions and data processing by running the software programs, instructions, and modules stored in the memory 42, that is, implements the method in the method embodiments described above.

[0208] The memory 42 can include a program storage area and a data storage area. The program storage area can store an operating system and at least one application required by a function; the data storage area can store data created according to the use of the terminal device, and the like. In addition, the memory 42 can include a high-speed random access memory, and can also include a non-volatile memory.

[0209] The embodiments of the present application provide a computer-readable storage medium, which stores computer-executable instructions, and the computer-executable instructions are configured to execute the method for confirming the atmospheric bending angle observation error of occultation stars or the use method provided in the foregoing embodiments.

[0210] The embodiments of the present application provide a computer program product, which includes a computer program stored on a computer-readable storage medium, and the computer program includes program instructions. When the program instructions are executed by a computer, the computer executes the method for confirming the atmospheric bending angle observation error of occultation stars or the use method provided in the foregoing embodiments.

[0211] The computer-readable storage medium described above can be a transitory computer-readable storage medium or a non-transitory computer-readable storage medium.

[0212] The technical solutions of the embodiments of the present application can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes one or more instructions to cause a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods in the embodiments of the present application. The storage medium mentioned above can be a non-transitory storage medium, including: a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, and various media that can store program codes, or can be a transitory storage medium.

[0213] The above description and drawings sufficiently illustrate the embodiments of the present application to enable one skilled in the art to practice them. Other embodiments can include structural, logical, electrical, process, and other changes. The embodiments are merely representative of the possible variations. Individual components and functions are optional unless explicitly required, and the order of operations can be varied. Portions and features of some embodiments can be included in, or alternative to, portions and features of other embodiments. Also, the words used in this specification are words of description, not limitation, and the phrase "including an" followed by a listing of components does not require that there be no others. As used in this application, the term "comprising" and its variants are used to indicate that the underlined feature is an optional element and the use of the term "comprising" does not mean that other elements cannot be included with the claimed disclosure. Without more limitations, an element defined by the phrase "including a" does not exclude the existence of additional identical elements in the process, method, or device including the element. In this document, each embodiment focuses on the differences from other embodiments, and the same or similar parts between embodiments can be referred to each other. For the methods, products, etc. disclosed in the embodiments, if they correspond to the method part disclosed in the embodiments, the relevant part can be referred to the description of the method part.

[0214] Those skilled in the art can appreciate that the units and algorithm steps of the examples described in combination with the embodiments disclosed herein can be realized in electronic hardware, or a combination of computer software and electronic hardware. Whether the functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the embodiments of the present application. Those skilled in the art can clearly appreciate that, for the convenience and brevity of description, the specific working processes of the above-described system, device and unit can refer to the corresponding processes in the foregoing method embodiments, which will not be repeated here.

[0215] In the embodiments disclosed herein, the disclosed methods, products (including but not limited to devices, equipment, etc.) can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of units can be only a logical function division. In actual implementation, another division mode can be used. For example, a plurality of units or components can be combined or integrated into another system, or some features can be omitted or not executed. In addition, the coupling or direct coupling or communication connection between the units shown or discussed can be indirect coupling or communication connection through some interface, device or unit, which can be electrical, mechanical or other form. The units described as separate components can or can not be physically separate, and the components shown as units can or can not be physical units, i.e. they can be located in one place, or distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to implement the embodiments. In addition, the functional units in the embodiments of the present application can be integrated in one processing unit, or each unit can be physically present separately, or two or more units can be integrated in one unit.

[0216] The flowcharts and block diagrams in the drawings show the possible implementation architecture, function and operation of the system, method and computer program product according to the embodiments of the present application. In this regard, each block in the flowchart or block diagram can represent a module, program segment or part of code containing one or more executable instructions for implementing the specified logical function. In some alternative implementations, the functions noted in the blocks can occur in different order from that shown in the drawings. For example, two consecutive blocks can actually be executed substantially in parallel, and sometimes they can be executed in reverse order, depending on the functions involved. Each block in the block diagram and / or flowchart, and the combination of blocks in the block diagram and / or flowchart, can be realized by a dedicated hardware-based system that performs the specified function or action, or can be realized by a combination of special-purpose hardware and computer instructions.

Claims

1. A method for confirming the observation error of the atmospheric curvature angle of an occultation, characterized in that, include: Obtain the atmospheric curvature angle of the occultation data, and the first and second meteorological datasets corresponding to the occultation profile in the occultation data that are not derived from the current occultation; wherein, the correlation between the sources of the first and second meteorological datasets is lower than a first preset threshold, and the sources of the two datasets are considered unrelated; the correlation between the source of the atmospheric curvature angle of the occultation and the first and second meteorological datasets is lower than a second preset threshold, so that the source of the atmospheric curvature angle is unrelated to the first and second meteorological datasets; A first atmospheric curvature angle is determined based on the first meteorological dataset, and a second atmospheric curvature angle is determined based on the second meteorological dataset. The first curvature angle error between the occultation atmospheric curvature angle and the first atmospheric curvature angle, the second curvature angle error between the occultation atmospheric curvature angle and the second atmospheric curvature angle, and the third curvature angle error between the first atmospheric curvature angle and the second atmospheric curvature angle are obtained. The observation error of the atmospheric curvature angle of an occultation is determined by the following formula: ； in, The observation error of the atmospheric curvature angle of the occultation is... The first bending angle error, This is the second bending angle error. This refers to the error of the third bending angle.

2. The method for confirming the observation error of the atmospheric curvature angle of an occultation according to claim 1, characterized in that, Determining a first atmospheric curvature angle based on the first meteorological dataset and a second atmospheric curvature angle based on the second meteorological dataset includes: The atmospheric refractive index is determined based on the current temperature, current relative humidity, and current atmospheric pressure in the first and second meteorological datasets, respectively; wherein the atmospheric refractive index corresponds one-to-one with altitude. Based on their respective atmospheric refractive indices, the first atmospheric curvature angle and the second atmospheric curvature angle are determined respectively.

3. The method for confirming the observation error of the atmospheric curvature angle of an occultation according to claim 2, characterized in that, Based on their corresponding atmospheric refractive indices, the first atmospheric curvature angle and the second atmospheric curvature angle are determined respectively, including: The refractive index profiles of the two are determined based on their respective refractive indices; the first atmospheric curvature angle and the second atmospheric curvature angle are determined based on the Abel integral equation, the refractive index profiles of the two and the altitude. or, The height below the preset height is divided into multiple collision height layers; wherein, the collision height of each collision height layer is determined based on the reference height, the local radius of curvature, and the atmospheric refractive index and altitude corresponding to the first meteorological dataset and the second meteorological dataset, respectively. The reference height is the height of the occultation point of the occultation profile data from the reference ellipsoid, and the local radius of curvature is the local radius of curvature of the reference ellipsoid of the last data in the occultation profile data. Based on the atmospheric refractive indices of the first and second meteorological datasets corresponding to each collision height layer, determine the refractive index gradient of each collision height layer corresponding to the first and second meteorological datasets. Based on the atmospheric refractive index, refractive index gradient, and collision height of each collision height layer corresponding to the first and second meteorological datasets, determine the atmospheric curvature angle of the first and second meteorological datasets in each collision height layer; The first atmospheric curvature angle is determined based on multiple atmospheric curvature angles corresponding to the first meteorological dataset, and the second atmospheric curvature angle is determined based on multiple atmospheric curvature angles corresponding to the second meteorological dataset.

4. The method for confirming the observation error of the atmospheric curvature angle of an occultation according to claim 3, characterized in that, The collision height of each collision height layer is determined as follows: ； in, The collision height layer number. For the first The collision height of each collision height layer For the first meteorological dataset or the second meteorological dataset corresponding to the first The altitude of the collision height layer The reference height, Let be the local radius of curvature. For the first meteorological dataset or the second meteorological dataset corresponding to and the Atmospheric refractive index of the collision height layer; And / or, Based on the atmospheric refractive indices of the first and second meteorological datasets corresponding to each collision height layer, the refractive index gradient of each collision height layer is determined, including: ； in, The collision height layer number. For the first meteorological dataset or the second meteorological dataset, the corresponding first... The refractive index gradient of the collision height layer For the first meteorological dataset or the second meteorological dataset, the corresponding first... Atmospheric refractive index at the collision altitude layer For the corresponding to the first Atmospheric refractive index at the collision altitude layer For the first The collision height of each collision height layer For the first The collision height of each collision height layer; And / or, Based on the atmospheric refractive index, refractive index gradient, and collision height of each collision height layer corresponding to the first and second meteorological datasets, the atmospheric curvature angle corresponding to the first and second meteorological datasets in each collision height layer is determined, including: ； in, The collision height layer number. For the first meteorological dataset or the second meteorological dataset, the corresponding first... The atmospheric curvature angle of the collision altitude layer For the first The collision height of each collision height layer This represents the total number of collision height layers. For the first meteorological dataset or the second meteorological dataset, the corresponding first... The refractive index gradient of the collision height layer For the first The collision height of each collision height layer For the first meteorological dataset or the second meteorological dataset, the corresponding first... Atmospheric refractive index at the collision altitude layer For the first The collision height of each collision height layer.

5. The method for confirming the observation error of the atmospheric curvature angle of an occultation according to claim 2, characterized in that, Based on the current temperature, current relative humidity, and current atmospheric pressure in the first and second meteorological datasets, determine the corresponding atmospheric refractive indices for each dataset, including: ； ； ； ； ； ； in, This is the saturated water vapor pressure on the ice surface. This is the saturated vapor pressure at the water surface. The current temperature. This is the lower limit temperature of the pure aqueous phase of water. This is the upper limit temperature of the pure ice phase of water. for and Mixed weights, For mixed saturated water vapor pressure, This is the actual water vapor pressure. The current relative humidity. This is the current atmospheric pressure.

6. The method for confirming the observation error of the atmospheric curvature angle of an occultation according to claim 1, characterized in that, The methods for obtaining the first meteorological dataset and the second meteorological dataset include: In the time dimension, the occultation observation time is obtained by using linear weighted interpolation based on the data from two time points before and after the occultation observation time. In the horizontal spatial dimension, taking the latitude and longitude coordinates of the occultation tangent point as the center, the values ​​of four grid points around it are selected, and the position of the tangent point is calculated by bilinear interpolation algorithm. In the vertical spatial dimension, the occultation observation height layer is obtained by logarithmic pressure coordinate transformation and cubic spline interpolation method; From two unrelated meteorological data sources, the first meteorological data and the second meteorological data corresponding to the occultation observation time, the tangency point position, and the occultation observation altitude layer are obtained respectively. Based on the occultation profile, a number of first meteorological data and second meteorological data are obtained. The number of first meteorological data is combined into a first meteorological dataset, and the number of second meteorological data is combined into a second meteorological dataset.

7. The method for confirming the observation error of the atmospheric curvature angle of an occultation according to any one of claims 1 to 6, characterized in that, Occultation can receive signals from navigation satellites under different types of navigation systems; Obtaining the atmospheric curvature angle of the occultation data includes: obtaining the navigation system type corresponding to the source of the occultation data; Determine the atmospheric curvature angle of the occultation corresponding to the type of navigation system.

8. A device for confirming the observation error of the atmospheric curvature angle of an occultation, characterized in that, include: The first acquisition module is used to acquire the atmospheric curvature angle of the occultation data, and the first and second meteorological datasets corresponding to the occultation profile in the occultation data that are not derived from the current occultation; wherein, the correlation between the sources of the first and second meteorological datasets is lower than a first preset threshold, and the sources of the two datasets are considered unrelated; the correlation between the source of the atmospheric curvature angle of the occultation and the first and second meteorological datasets is lower than a second preset threshold, so that the source of the atmospheric curvature angle is unrelated to the first and second meteorological datasets; The first determining module is used to determine a first atmospheric curvature angle based on the first meteorological dataset and to determine a second atmospheric curvature angle based on the second meteorological dataset. The second obtaining module is used to obtain the first curvature angle error between the occultation atmospheric curvature angle and the first atmospheric curvature angle, the second curvature angle error between the occultation atmospheric curvature angle and the second atmospheric curvature angle, and the third curvature angle error between the first atmospheric curvature angle and the second atmospheric curvature angle; The second determining module is used to determine the observation error of the occultation atmospheric curvature angle using the following formula: ； in, The observation error of the atmospheric curvature angle of the occultation is... The first bending angle error, This is the second bending angle error. This refers to the error of the third bending angle.

9. A method for using the observation error of the atmospheric curvature angle of an occultation, characterized in that, include: Obtain the navigation system type, area range, and altitude range input by the user, and determine the mapping curvature angle error based on the fitting function; The mapping curvature angle error is determined as follows: the occultation atmospheric curvature angle observation error is mapped to the mapping curvature angle error according to the preset mapping method; wherein the preset mapping method corresponds to the usage method of the occultation atmospheric curvature angle observation error. The fitting function is determined as follows: The types of navigation systems capable of receiving signals from occultations are obtained; the Earth is divided into multiple regions based on latitude; and the atmospheric altitude is divided into multiple altitude ranges based on altitude. For each navigation system type, each region, and each altitude range, the occultation atmospheric curvature angle observation error is determined using the verification method described in any one of claims 1 to 7. Each occultation atmospheric curvature angle observation error is mapped to a mapped curvature angle error. Based on a preset fitting method, using the one-to-one corresponding mapped curvature angle, navigation system type, region, and altitude range as the fitting dataset, a fitting function is fitted with navigation system type, region, and altitude range as independent variables and mapped curvature angle as the dependent variable.

10. A device for confirming the observation error of the atmospheric curvature angle of an occultation, comprising a processor and a memory storing program instructions, characterized in that, The processor is configured to, when executing the program instructions, perform the method for confirming the observation error of the occultation atmospheric curvature angle as described in any one of claims 1 to 7; or perform the method for using the observation error of the occultation atmospheric curvature angle as described in claim 9.

Images