CN116931006A - Troposphere delay value determining method, system, equipment and storage medium - Google Patents
Troposphere delay value determining method, system, equipment and storage medium Download PDFInfo
- Publication number
- CN116931006A CN116931006A CN202310926477.1A CN202310926477A CN116931006A CN 116931006 A CN116931006 A CN 116931006A CN 202310926477 A CN202310926477 A CN 202310926477A CN 116931006 A CN116931006 A CN 116931006A
- Authority
- CN
- China
- Prior art keywords
- value
- delay value
- zenith
- water vapor
- vapor pressure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000005436 troposphere Substances 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012937 correction Methods 0.000 claims abstract description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 92
- 230000008859 change Effects 0.000 claims description 73
- 230000003247 decreasing effect Effects 0.000 claims description 54
- 230000006870 function Effects 0.000 claims description 29
- 238000013507 mapping Methods 0.000 claims description 27
- 238000001035 drying Methods 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 101001072091 Homo sapiens ProSAAS Proteins 0.000 description 1
- 102100036366 ProSAAS Human genes 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012731 temporal analysis Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000000700 time series analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/40—Correcting position, velocity or attitude
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/02—Instruments for indicating weather conditions by measuring two or more variables, e.g. humidity, pressure, temperature, cloud cover or wind speed
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the application discloses a method, a system, equipment and a storage medium for determining a troposphere delay value, which can effectively improve the correction precision of the conventional troposphere experience model by combining a weather observation file and a zenith humidity delay value obtained by the aid of the troposphere experience model.
Description
[ field of technology ]
The application relates to the field of geodetic measurement and satellite navigation, in particular to a tropospheric delay value determination method, a system, equipment and a storage medium.
[ background Art ]
Tropospheric delay generally refers to refraction of electromagnetic wave signals through a non-ionized atmosphere, slowing down and yielding a curved path, which is one of the main sources of error in GNSS (Global Navigation Satellite System, global navigation satellite data) processing. The errors produced by the non-ionized atmosphere cannot be used to eliminate tropospheric delay by the observed value combination method, but rather correction of tropospheric delay is needed by modeling and combining parameter estimation.
In recent years, the existing GNSS data processing strategy selects a single troposphere empirical model for delay correction, and only the annual scale and half annual scale change are considered when the single troposphere empirical model is modeled, so that the deviation of the troposphere delay value is large, and the actual change trend of weather in a short time of a weather observation station is difficult to reflect.
[ application ]
In view of the above, the application provides a method, a system, a device and a storage medium for determining tropospheric delay values, which can effectively improve correction accuracy of an existing tropospheric empirical model, and compared with the existing tropospheric empirical model, the accuracy of a meteorological observation file in the embodiment of the application is high, so that the obtained tropospheric delay values are more accurate, and actual change trend of weather in a short time of a meteorological observation station can be reflected.
The specific technical scheme of the first embodiment of the application is as follows: a method of tropospheric delay value determination, the method comprising: acquiring an air pressure value and a water vapor pressure value of a weather observation station at a preset time point according to a weather observation file; acquiring a longitude value, a latitude value, a geodetic elevation value and a yearly product daily value of the weather observation station; obtaining a specific constant of dry gas, a water vapor pressure decreasing factor, a weighted average temperature and a temperature decreasing rate of the weather observation station in a preset time period and an average temperature and a pixel change coefficient of sea level in the preset time period by using the longitude value, the latitude value, the geodetic altitude value, the annual daily value and a preset troposphere empirical model; the water vapor pressure decreasing factor is the average attenuation rate of the water vapor pressure value; the preset time point is within the preset time period; the meteorological element change coefficient is used for representing the meteorological change amplitude of the meteorological observation station; obtaining a zenith dry delay value of the weather observation station according to the air pressure value, the latitude value and the geodetic altitude value; obtaining a zenith humidity delay value of the weather observation station according to the air pressure value, the water vapor pressure value, the drying gas specific constant, the water vapor pressure decreasing factor, the weighted average temperature, the temperature decreasing rate and the average temperature of the sea level; and obtaining the troposphere delay value of the weather observation station according to the weather element change coefficient, the zenith dry delay value and the zenith wet delay value.
Preferably, the obtaining a specific constant of the dry gas, a water vapor pressure decreasing factor, a weighted average temperature, a temperature decreasing rate of the weather observation station in a preset time period, and an average temperature of the sea level and a pixel change coefficient in the preset time period by using the longitude value, the latitude value, the geodetic altitude, the annual daily value and a preset troposphere empirical model includes: obtaining a specific constant of the dry gas, a weather standard value of the weather observation station in a preset time period, a water vapor pressure decreasing factor, a weighted average temperature, a temperature decreasing rate and an average temperature of sea level by using the longitude value, the latitude value, the geodetic altitude value, the annual product daily value and a preset troposphere empirical model; the weather standard value at least comprises a barometric pressure standard value, a temperature standard value and a water vapor pressure standard value; acquiring an air pressure element value of the weather observation station according to a preset earth surface weather model and the air pressure standard value; the air pressure element value comprises an air pressure annual average value, an air pressure annual period change amplitude and an air pressure half-year period change amplitude; obtaining a temperature element value of the weather observation station according to a preset earth surface weather model and the temperature standard value; the temperature element value comprises a temperature annual average value, a temperature annual period change amplitude and a temperature half-annual period change amplitude; obtaining a water vapor pressure element value of the weather observation station according to a preset earth surface weather model and the water vapor pressure standard value; the water vapor pressure element value comprises a water vapor pressure annual average value, a water vapor pressure annual cycle variation amplitude and a water vapor pressure half-year cycle variation amplitude; and fixing the air pressure element value, the temperature element value and the water vapor pressure element value by using a least square method to obtain a weather element change coefficient of the weather observation station.
Preferably, the zenith dry delay value is obtained by adopting the following formula:
wherein ZHD is the zenith dry delay value, P is the air pressure value,for the latitude value, H is the geodetic high value, k 1 、k 2 And k 3 Are all constant.
Preferably, the zenith wet retardation value is obtained by adopting the following formula:
wherein ZWD is the zenith wet delay value, R d For the specific constant of the drying gas, P is the air pressure value, lambda is the water vapor pressure decreasing factor, alpha is the temperature decreasing rate, e is the water vapor pressure value, T s T is the average temperature of sea level m G is the weighted average temperature m Gravitational acceleration, k 4 、k 5 、k 6 And k 7 Are all constant.
Preferably, the obtaining the tropospheric delay value of the weather observation station according to the weather element change coefficient, the zenith dry delay value and the zenith wet delay value includes: constructing a zenith dry delay value mapping function by using the meteorological element change coefficient, the altitude angle between the meteorological observation station and a target satellite, the geodetic altitude value and the correction coefficient of a preset zenith dry delay value on different components; constructing zenith wet delay value mapping functions by using the meteorological element change coefficients, the altitude angle between the meteorological observation station and a target satellite, the geodetic altitude value and correction coefficients of preset zenith wet delay values on different components; and obtaining the troposphere delay value of the weather observation station by using the zenith dry delay value mapping function, the zenith dry delay value, the zenith wet delay value mapping function and the zenith wet delay value.
Preferably, the zenith dry delay value mapping function is obtained by adopting the following formula:
wherein f h (E) Mapping function a for zenith dry delay value h B is the change coefficient of meteorological elements h And c h E is the altitude angle of the meteorological station and the target satellite, h is a constantFor the ground high value, a ht 、b ht And c ht And (3) a preset correction coefficient for the zenith dry delay value on different components.
Preferably, the tropospheric delay value is obtained using the following formula:
ρ delay =ρ ZHD f h (E)+ρ ZWD f w (E)
wherein ρ is delay For the tropospheric delay value ρ ZHD For the zenith dry delay value, f h (E) Mapping the zenith dry delay value into a function ρ ZWD For the zenith wet retardation value, f w (E is the zenith wet delay value mapping function.
The specific technical scheme of the second embodiment of the application is as follows: a tropospheric delay value determination system, the system comprising: the system comprises a data acquisition module, a prediction module, a zenith dry delay value determination module, a zenith wet delay value determination module and a troposphere delay value determination module; the data acquisition module is used for acquiring an air pressure value and a water vapor pressure value of a weather observation station at a preset time point according to a weather observation file, and acquiring a longitude value, a latitude value, a geodetic altitude value and a yearly product daily value of the weather observation station; the prediction module is used for obtaining a specific constant of dry gas, a water vapor pressure decreasing factor, a weighted average temperature and a temperature decreasing rate of the weather observation station in a preset time period and an average temperature and a pixel change coefficient of sea level in the preset time period by using the longitude value, the latitude value, the earth elevation value, the annual energy daily value and a preset troposphere empirical model; the water vapor pressure decreasing factor is the average attenuation rate of the water vapor pressure value; the preset time point is within the preset time period; the meteorological element change coefficient is used for representing the meteorological change amplitude of the meteorological observation station; the zenith dry delay value determining module is used for obtaining a zenith dry delay value of the weather observation station according to the air pressure value, the latitude value and the geodetic altitude value; the zenith humidity delay value determining module is used for obtaining a zenith humidity delay value of the weather observation station according to the air pressure value, the water vapor pressure value, the drying gas specific constant, the water vapor pressure decreasing factor, the weighted average temperature, the temperature decreasing rate and the average temperature of the sea level; and the tropospheric delay value determining module is used for obtaining the tropospheric delay value of the weather observation station according to the weather element change coefficient, the zenith dry delay value and the zenith wet delay value.
The specific technical scheme of the third embodiment of the application is as follows: a tropospheric delay value determination apparatus comprising a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to perform the steps of the method according to any one of the third embodiments of the present application.
The specific technical scheme of the fourth embodiment of the application is as follows: a computer readable storage medium storing a computer program which, when executed by a processor, causes the processor to perform the steps of the method according to any of the first embodiments of the application.
The implementation of the embodiment of the application has the following beneficial effects:
according to the method, the air pressure value and the water vapor pressure value of a weather observation station at a preset time point are obtained according to weather observation files; acquiring a longitude value, a latitude value, a geodetic elevation value and a yearly product daily value of a weather observation station; obtaining a specific constant of dry gas, a water vapor pressure decreasing factor, a weighted average temperature and a temperature decreasing rate of a weather observation station in a preset time period, and an average temperature of sea level and a pixel element change coefficient in the preset time period by using a longitude value, a latitude value, a geodetic altitude value, a yearup daily value and a preset troposphere empirical model; the water vapor pressure decreasing factor is the average attenuation rate of the water vapor pressure value; the preset time point is within a preset time period; the meteorological element change coefficient is used for representing the meteorological change amplitude of the meteorological observation station; obtaining a zenith dry delay value of the weather observation station according to the air pressure value, the latitude value and the geodetic altitude value; obtaining a zenith humidity delay value of the weather observation station according to the air pressure value, the water vapor pressure value, the specific constant of dry gas, the water vapor pressure decreasing factor, the weighted average temperature, the temperature decreasing rate and the average temperature of the sea level; and obtaining the troposphere delay value of the weather observation station according to the weather element change coefficient, the zenith dry delay value and the zenith wet delay value.
The weather observation file can accurately reflect weather data of the weather observation station in any time period, the zenith dry delay value and the zenith wet delay value obtained by combining the weather observation file and the current layer experience model can effectively improve correction precision of the current troposphere experience model, and compared with the current troposphere experience model, the weather observation file in the embodiment of the application has high precision, so that the obtained troposphere delay value is more accurate, and the actual change trend of weather in a short time of the weather observation station can be reflected.
[ description of the drawings ]
In order to more clearly illustrate the embodiments of the application or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of steps of a method for determining tropospheric delay values;
FIG. 2 is a flowchart showing steps for obtaining the change coefficients of meteorological elements;
FIG. 3 is a flowchart showing steps for obtaining tropospheric delay values;
FIG. 4a is a schematic diagram of air pressure values of different models;
FIG. 4b is a graph showing temperature values of different models;
FIG. 4c is a graph showing water vapor pressure values of different models;
FIG. 5a is a time series of zenith stem delay values;
FIG. 5b is a time series of zenith wet delay values;
FIG. 5c is a tropospheric delay value time series;
FIG. 6 is a schematic diagram of a tropospheric delay value determination system;
601, a data acquisition module; 602. a prediction module; 603. a zenith dry delay value determining module; 604. a zenith wet delay value determining module; 605. and a tropospheric delay value determination module.
[ detailed description ] of the application
The following description of the embodiments of the present application will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present application, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
Referring to fig. 1, a method for determining a tropospheric delay value according to a first embodiment of the present application can reflect an actual trend of weather in a short time at a weather station, the method comprising:
step 101, acquiring an air pressure value and a water vapor pressure value of a weather observation station at a preset time point according to a weather observation file;
102, acquiring a longitude value, a latitude value, a geodetic elevation value and a yearly product daily value of the weather observation station;
step 103, obtaining a specific constant of dry gas, a water vapor pressure decreasing factor, a weighted average temperature and a temperature decreasing rate of the weather observation station in a preset time period and an average temperature and a pixel change coefficient of sea level in the preset time period by using the longitude value, the latitude value, the geodetic altitude value, the annual daily value and a preset troposphere empirical model; the water vapor pressure decreasing factor is the average attenuation rate of the water vapor pressure value; the preset time point is within the preset time period; the meteorological element change coefficient is used for representing the meteorological change amplitude of the meteorological observation station;
104, obtaining a zenith dry delay value of the weather observation station according to the air pressure value, the latitude value and the geodetic altitude value;
step 105, obtaining a zenith humidity delay value of the weather observation station according to the air pressure value, the water vapor pressure value, the drying gas specific constant, the water vapor pressure decreasing factor, the weighted average temperature, the temperature decreasing rate and the average temperature of the sea level;
and 106, obtaining a tropospheric delay value of the weather observation station according to the weather element change coefficient, the zenith dry delay value and the zenith wet delay value.
Specifically, the barometric pressure value P, the temperature T, the relative humidity hr and the water vapor pressure value e of the weather observation station at a preset time point are obtained according to a weather observation file obtained by the global navigation satellite system service, and the water vapor pressure e is calculated by the following formula to be e=hr x exp (-37.2465+0.213266T-0.256908 x 10 -3 )T 2 ) The method comprises the steps of carrying out a first treatment on the surface of the Acquiring a longitude value, a latitude value, a geodetic altitude value h and a yearly product daily value of a weather observation station; obtaining a specific constant R of the dry gas by using a longitude value, a latitude value, a geodetic altitude value, a yearling daily value and a preset tropospheric empirical model GPT3 d The water vapor pressure of the weather observation station decreases by a factor lambda and the weighted average temperature T within a preset time period m Temperature decrease rate a, average temperature T of sea level within a predetermined period of time s And the pixel change coefficient a h The method comprises the steps of carrying out a first treatment on the surface of the GPT3 is an empirical weather model established according to weather re-analysis data of a weather forecast center, is the GPT series model with highest precision, and can be directly used; obtaining a zenith dry delay value ZHD of the weather observation station according to the air pressure value P, the latitude value and the geodetic altitude value h; according to the air pressure value P, the water vapor pressure value e and the specific constant R of the dry gas d A water vapor pressure decreasing factor lambda, a weighted average temperature T m Temperature decrease rate alpha, average temperature T at sea level s Obtaining a zenith wet delay value ZWD of a weather observation station; according to the change coefficient a of the meteorological element h The zenith dry delay value ZHD and zenith wet delay value ZWD obtain tropospheric delay values for the weather observation station.
The weather observation file can accurately reflect weather data of the weather observation station in any time period, the zenith dry delay value and the zenith wet delay value obtained by combining the weather observation file and the current layer experience model can effectively improve correction precision of the current troposphere experience model, and compared with the current troposphere experience model, the weather observation file in the embodiment of the application has high precision, so that the obtained troposphere delay value is more accurate, and the actual change trend of weather in a short time of the weather observation station can be reflected.
In a specific embodiment, referring to fig. 2, in step 103, a specific constant of the dry gas, a water vapor pressure decreasing factor, a weighted average temperature, a temperature decreasing rate of the weather observation station in a preset time period, and an average temperature of the sea level and a pixel element change coefficient in the preset time period are obtained by using the longitude value, the latitude value, the geodetic altitude value, the annual product daily value and a preset tropospheric experience model, including:
step 201, obtaining a specific constant of dry gas, a weather standard value of a weather observation station, a water vapor pressure decreasing factor, a weighted average temperature, a temperature decreasing rate and an average temperature of sea level in a preset time period by using the longitude value, the latitude value, the geodetic altitude value, the annual energy daily value and a preset troposphere empirical model; the weather standard value at least comprises a barometric pressure standard value, a temperature standard value and a water vapor pressure standard value;
step 202, obtaining an air pressure element value of the weather observation station according to a preset earth surface weather model and the air pressure standard value; the air pressure element value comprises an air pressure annual average value, an air pressure annual period change amplitude and an air pressure half-year period change amplitude;
step 203, obtaining a temperature element value of the weather observation station according to a preset earth surface weather model and the temperature standard value; the temperature element value comprises a temperature annual average value, a temperature annual period change amplitude and a temperature half-annual period change amplitude;
step 204, obtaining a water vapor pressure element value of the weather observation station according to a preset earth surface weather model and the water vapor pressure standard value; the water vapor pressure element value comprises a water vapor pressure annual average value, a water vapor pressure annual cycle variation amplitude and a water vapor pressure half-year cycle variation amplitude;
and 205, fixing the air pressure element value, the temperature element value and the water vapor pressure element value by using a least square method to obtain a weather element change coefficient of the weather observation station.
Specifically, the formula of the preset earth surface meteorological model is as follows:
wherein t represents the annual product daily value, and the air pressure standard value is taken as qt to be input into a preset earth surface meteorological model to obtain an air pressure element value, and the air pressure element value comprises an air pressure annual average value A 0 Amplitude of annual change of air pressure A 1 、B 1 And amplitude of half-year cycle change of air pressure A 2 、B 2 The method comprises the steps of carrying out a first treatment on the surface of the Similarly, respectively inputting a temperature standard value and a water vapor pressure standard value into a surface meteorological model to obtain a temperature element value and a water vapor pressure element value; and fixing the air pressure element value, the temperature element value and the water vapor pressure element value by using a least square method to obtain a weather element change coefficient of the weather observation station. The least square method can find the optimal function matching of the data by minimizing the square sum of errors, and the optimal weather change coefficient can be obtained by fixing the air pressure element value, the temperature element value and the water vapor pressure element value by using the least square method.
In a specific embodiment, the zenith dry delay value is obtained using the following formula:
wherein ZHD is the zenith dry delay value, P is the air pressure value,for the latitude value, H is a geodetic high value, and the unit is km, k 1 、k 2 And k 3 Are all constant.
Specifically, a Saastamoinen model is utilized to obtain a zenith dry delay value, k 1 、k 2 And k 3 Is a fixed constant. Specifically, the formula for obtaining the zenith dry delay value is as follows:
in a specific embodiment, the zenith wet delay value is obtained using the following formula:
wherein ZWD is the zenith wet delay value, R d For the specific constant of the drying gas, P is the air pressure value, lambda is the water vapor pressure decreasing factor, alpha is the temperature decreasing rate, e is the water vapor pressure value, T s T is the average temperature of sea level m G is the weighted average temperature m Gravitational acceleration, k 4 、k 5 、k 6 And k 7 Are all constant.
Specifically, A is used&N formula calculates zenith wet delay value, k 4 、k 5 、k 6 And k 7 Can be set empirically, wherein g m =9.80665m/s 2 ,k 4 =77.604.0K/mbar,k 5 =16.5221K/mbar,k 6 =377600.0K 2 /mbar,R d =287.0464。
In a specific embodiment, referring to fig. 3, in step 106, obtaining a tropospheric delay value of the weather observation station according to the weather element change coefficient, the zenith dry delay value and the zenith wet delay value includes:
step 301, constructing a zenith dry delay value mapping function by using the meteorological element change coefficient, the altitude angle between the meteorological observation station and a target satellite, the geodetic altitude value and the correction coefficient of a preset zenith dry delay value on different components;
step 302, constructing a zenith humidity delay value mapping function by using the meteorological element change coefficient, the altitude angle between the meteorological observation station and a target satellite, the geodetic altitude value and correction coefficients of a preset zenith humidity delay value on different components;
and 303, obtaining a tropospheric delay value of the weather observation station by using the zenith dry delay value mapping function, the zenith dry delay value, the zenith wet delay value mapping function and the zenith wet delay value.
Specifically, the weather element change coefficient can better show the weather change characteristic of the area, and the accuracy of the obtained troposphere delay value can be improved by integrating the weather element change coefficient, the zenith dry delay value and the correction coefficient of the zenith wet delay value on different components.
In a specific embodiment, the zenith dry delay value mapping function is obtained by adopting the following formula:
wherein f h (E) Mapping function a for zenith dry delay value h B is the change coefficient of meteorological elements h And c h Is constant, E is the altitude angle of the weather observation station and the target satellite, h is the ground high value, a ht 、b ht And c ht And (3) a preset correction coefficient for the zenith dry delay value on different components.
In a specific embodiment, the tropospheric delay value is obtained using the following formula:
ρ delay =ρ ZHD f h (E+ρ ZWD f w (E)
wherein ρ is delay For the tropospheric delay value ρ ZHD For the zenith dry delay value, f h (E) Mapping the zenith dry delay value into a function ρ ZWD For the zenith wet retardation value, f w (E) Mapping a function for the zenith wet delay value.
Referring to fig. 4a, fig. 4b and fig. 4c, metafile is a weather observation file, which can reflect weather changes of a weather observation station, and the temperature and the water vapor pressure of the metafile have large changes in a short time, for example: 20-30h, the temperature is raised by about 7 ℃; within 40-50h, the vapor pressure is reduced by about 5hPa; whereas the atmospheric pressure varies maximally by about 6hPa in 0-80 h. GPT3 is used as an empirical weather model, and weather elements have no obvious change trend in a short time, so that actual weather element changes are difficult to reflect. VMFsite is used as a numerical meteorological model (ideal model), the time resolution is 6 hours, and compared with GPT3, the VMFsite has better meteorological change trend.
In order to truly and effectively reflect the accuracy of the tropospheric delay value determination method in this embodiment, observation data of WUH2 stations (30.5°n,114.4°e,28.2 m) for 3 days are selected for testing, and five schemes for acquiring tropospheric delay values are designed in total, as shown in table 1:
table 1 tropospheric delay value acquisition scheme
The method comprises the steps that a scheme 1 is to obtain a tropospheric delay value by using an empirical model GPT3, a scheme 2 is to obtain a tropospheric delay value by using a meteorological observation file, a scheme 3 and a scheme 4 are numerical meteorological models, the precision is higher than that of the empirical model GPT3, but the time cost of the numerical model in data processing is higher than that of the empirical model GPT3, and the numerical model is used as a reference and has no contrast significance; scheme 5 is to combine GPT3 model and meteorological observation file to obtain troposphere delay value in this embodiment, wherein SAAS and A & N are existing zenith troposphere delay model, and the function is to calculate troposphere dry and wet delay value respectively; GMF and VMF1 are an empirical mapping function model and a numerical mapping function model, respectively, and function to obtain the mapping amount of the tropospheric delay value in the direction of the inclined path.
TABLE 2 error comparison of different tropospheric delay models
In table 2, RMS is the error, ZHD is the zenith dry delay value, ZHD is the zenith wet delay value, ZTD troposphere delay value, as can be seen from table 2, the troposphere delay value in scheme five is less than that in schemes one and two; wherein, the root mean square of the scheme II and the scheme V ZHD is consistent and is about 3.5mm better than the scheme I; in the calculation of ZWD and ZTD, the second scheme has the worst precision, the first scheme is improved by about 33.6% compared with the second scheme, and the first scheme is improved by 9.3%. The time series analysis of the troposphere delay value errors obtained by different schemes is shown in fig. 5a, 5b and 5c, and the experimental results are obtained by comparing and analyzing as follows:
(1) ZHD the calculation is mainly affected by meteorological elements: ZHD calculated by the meteorological observation file has the best precision, the numerical model is inferior, and the precision difference of the empirical meteorological model GPT3 is the largest;
(2) ZWD calculations are affected by both meteorological elements and models: comparing the scheme I with the scheme V, and improving the water vapor pressure accuracy obtained by using the actually measured water vapor pressure in the weather observation file by 9.3% compared with the water vapor pressure obtained by GPT 3; comparison of the scheme II and the scheme five shows that the precision of the scheme 5 is improved by 33.6%;
(3) Table 2 in combination with fig. 5a, 5b and 5c, it can be seen that when the tropospheric delay value is calculated in scheme 5, the ZTD root mean square value is smaller than ZHD, which is due to the fact that ZHD and ZWD compensate each other in the calculation process, so that the accuracy in scheme 5 is the highest.
Referring to fig. 6, a schematic structural diagram of a tropospheric delay value determining system according to a second embodiment of the present application is provided, the system includes: a data acquisition module 601, a prediction module 602, a zenith dry delay value determination module 603, a zenith wet delay value determination module 604 and a tropospheric delay value determination module 605; the data acquisition module 601 is configured to acquire an air pressure value and a water vapor pressure value of a weather-observing station at a preset time point according to a weather-observing file, and acquire a longitude value, a latitude value, a geodetic altitude value and a yearly-product daily value of the weather-observing station; the prediction module 602 is configured to obtain a specific constant of dry gas, a water vapor pressure decreasing factor, a weighted average temperature, a temperature decreasing rate of a weather observation station in a preset time period, and an average temperature of sea level and a pixel element change coefficient in the preset time period by using a longitude value, a latitude value, a geodetic altitude value, an annual daily value and a preset troposphere empirical model; the water vapor pressure decreasing factor is the average attenuation rate of the water vapor pressure value; the preset time point is within a preset time period; the meteorological element change coefficient is used for representing the meteorological change amplitude of the meteorological observation station; the zenith dry delay value determining module 603 is configured to obtain a zenith dry delay value of the weather observation station according to the barometric pressure value, the latitude value and the geodetic altitude value; the zenith humidity delay value determining module 604 is configured to obtain a zenith humidity delay value of the weather observation station according to the air pressure value, the water vapor pressure value, the specific constant of the dry gas, the water vapor pressure decreasing factor, the weighted average temperature, the temperature decreasing rate and the average temperature of the sea level; the tropospheric delay value determination module 605 is configured to obtain a tropospheric delay value for the weather observation station based on the weather element variation coefficient, the zenith dry delay value, and the zenith wet delay value.
By using the system in the embodiment, the zenith dry delay value and the zenith wet delay value can be obtained by combining the meteorological observation file and the current layer experience model, so that the correction precision of the current troposphere experience model can be effectively improved.
A third embodiment provided by the present application provides a tropospheric delay value determination apparatus comprising a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to perform the steps of the method according to any one of the first embodiments.
A fourth embodiment of the present application provides a computer-readable storage medium storing a computer program, characterized in that the computer program, when executed by a processor, causes the processor to perform the steps of the method according to any of the first embodiments.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The foregoing examples illustrate only a few embodiments of the application, which are described in detail and are not to be construed as limiting the scope of the application. It will be apparent to those skilled in the art that various modifications and improvements can be made without departing from the spirit of the application. Accordingly, the scope of protection of the present application is to be determined by the appended claims.
Claims (10)
1. A method for determining a tropospheric delay value, the method comprising:
acquiring an air pressure value and a water vapor pressure value of a weather observation station at a preset time point according to a weather observation file;
acquiring a longitude value, a latitude value, a geodetic elevation value and a yearly product daily value of the weather observation station;
obtaining a specific constant of dry gas, a water vapor pressure decreasing factor, a weighted average temperature and a temperature decreasing rate of the weather observation station in a preset time period and an average temperature and a pixel change coefficient of sea level in the preset time period by using the longitude value, the latitude value, the geodetic altitude value, the annual daily value and a preset troposphere empirical model; the water vapor pressure decreasing factor is the average attenuation rate of the water vapor pressure value; the preset time point is within the preset time period; the meteorological element change coefficient is used for representing the meteorological change amplitude of the meteorological observation station;
obtaining a zenith dry delay value of the weather observation station according to the air pressure value, the latitude value and the geodetic altitude value;
obtaining a zenith humidity delay value of the weather observation station according to the air pressure value, the water vapor pressure value, the drying gas specific constant, the water vapor pressure decreasing factor, the weighted average temperature, the temperature decreasing rate and the average temperature of the sea level;
and obtaining the troposphere delay value of the weather observation station according to the weather element change coefficient, the zenith dry delay value and the zenith wet delay value.
2. The tropospheric delay value determination method of claim 1, wherein said obtaining a specific constant of dry gas, a water vapor pressure decreasing factor of the weather observation station during a predetermined period, a weighted average temperature, a temperature decreasing rate, and an average temperature of sea level and a pixel change coefficient during the predetermined period using the longitude value, latitude value, altitude value, yearup value, and a predetermined tropospheric empirical model comprises:
obtaining a specific constant of the dry gas, a weather standard value of the weather observation station in a preset time period, a water vapor pressure decreasing factor, a weighted average temperature, a temperature decreasing rate and an average temperature of sea level by using the longitude value, the latitude value, the geodetic altitude value, the annual product daily value and a preset troposphere empirical model; the weather standard value at least comprises a barometric pressure standard value, a temperature standard value and a water vapor pressure standard value;
acquiring an air pressure element value of the weather observation station according to a preset earth surface weather model and the air pressure standard value; the air pressure element value comprises an air pressure annual average value, an air pressure annual period change amplitude and an air pressure half-year period change amplitude;
obtaining a temperature element value of the weather observation station according to a preset earth surface weather model and the temperature standard value; the temperature element value comprises a temperature annual average value, a temperature annual period change amplitude and a temperature half-annual period change amplitude;
obtaining a water vapor pressure element value of the weather observation station according to a preset earth surface weather model and the water vapor pressure standard value; the water vapor pressure element value comprises a water vapor pressure annual average value, a water vapor pressure annual cycle variation amplitude and a water vapor pressure half-year cycle variation amplitude;
and fixing the air pressure element value, the temperature element value and the water vapor pressure element value by using a least square method to obtain a weather element change coefficient of the weather observation station.
3. The tropospheric delay value determination method of claim 1, wherein the zenith dry delay value is obtained using the formula:
wherein ZHD is the zenith dry delay value, P is the air pressure value,for the latitude value, H is the geodetic high value, k 1 、k 2 And k 3 Are all constant.
4. The tropospheric delay value determination method of claim 1, wherein the zenith wet delay value is obtained using the formula:
wherein ZWD is the zenith wet delay value, P d For the specific constant of the drying gas, P is the air pressure value, lambda is the water vapor pressure decreasing factor, alpha is the temperature decreasing rate, e is the water vapor pressure value, T s T is the average temperature of sea level m G is the weighted average temperature m Gravitational acceleration, k 4 、k 5 、k 6 And k 7 Are all constant.
5. The tropospheric delay value determination method of claim 1, wherein the obtaining the tropospheric delay value of the weather observation station from the weather element variation coefficient, the zenith dry delay value, and the zenith wet delay value comprises:
constructing a zenith dry delay value mapping function by using the meteorological element change coefficient, the altitude angle between the meteorological observation station and a target satellite, the geodetic altitude value and the correction coefficient of a preset zenith dry delay value on different components;
constructing zenith wet delay value mapping functions by using the meteorological element change coefficients, the altitude angle between the meteorological observation station and a target satellite, the geodetic altitude value and correction coefficients of preset zenith wet delay values on different components;
and obtaining the troposphere delay value of the weather observation station by using the zenith dry delay value mapping function, the zenith dry delay value, the zenith wet delay value mapping function and the zenith wet delay value.
6. The tropospheric delay value determination method of claim 5, wherein the zenith dry delay value mapping function is obtained using the following equation:
wherein f h (E) Mapping function a for zenith dry delay value h B is the change coefficient of meteorological elements h And c h Is constant, E is the altitude angle of the weather observation station and the target satellite, h is the ground high value, a ht 、b ht And c ht And (3) a preset correction coefficient for the zenith dry delay value on different components.
7. The tropospheric delay value determination method of claim 5, wherein the tropospheric delay value is obtained using the following equation:
ρ delay =ρ ZHD f h (E)+ρ ZWD f w (E)
wherein ρ is delay For the tropospheric delay value ρ ZHD For the zenith dry delay value, f h (E) Mapping the zenith dry delay value into a function ρ ZWD For the zenith wet retardation value, f w (E) Mapping a function for the zenith wet delay value.
8. A tropospheric delay value determination system, the system comprising: the system comprises a data acquisition module, a prediction module, a zenith dry delay value determination module, a zenith wet delay value determination module and a troposphere delay value determination module;
the data acquisition module is used for acquiring an air pressure value and a water vapor pressure value of a weather observation station at a preset time point according to a weather observation file, and acquiring a longitude value, a latitude value, a geodetic altitude value and a yearly product daily value of the weather observation station;
the prediction module is used for obtaining a specific constant of dry gas, a water vapor pressure decreasing factor, a weighted average temperature and a temperature decreasing rate of the weather observation station in a preset time period and an average temperature and a pixel change coefficient of sea level in the preset time period by using the longitude value, the latitude value, the earth elevation value, the annual energy daily value and a preset troposphere empirical model; the water vapor pressure decreasing factor is the average attenuation rate of the water vapor pressure value; the preset time point is within the preset time period; the meteorological element change coefficient is used for representing the meteorological change amplitude of the meteorological observation station;
the zenith dry delay value determining module is used for obtaining a zenith dry delay value of the weather observation station according to the air pressure value, the latitude value and the geodetic altitude value;
the zenith humidity delay value determining module is used for obtaining a zenith humidity delay value of the weather observation station according to the air pressure value, the water vapor pressure value, the drying gas specific constant, the water vapor pressure decreasing factor, the weighted average temperature, the temperature decreasing rate and the average temperature of the sea level;
and the tropospheric delay value determining module is used for obtaining the tropospheric delay value of the weather observation station according to the weather element change coefficient, the zenith dry delay value and the zenith wet delay value.
9. Tropospheric delay value determination apparatus comprising a memory and a processor, characterized in that the memory stores a computer program which, when executed by the processor, causes the processor to perform the steps of the method according to any one of claims 1 to 7.
10. A computer readable storage medium storing a computer program, which when executed by a processor causes the processor to perform the steps of the method according to any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310926477.1A CN116931006A (en) | 2023-07-26 | 2023-07-26 | Troposphere delay value determining method, system, equipment and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310926477.1A CN116931006A (en) | 2023-07-26 | 2023-07-26 | Troposphere delay value determining method, system, equipment and storage medium |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116931006A true CN116931006A (en) | 2023-10-24 |
Family
ID=88392184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310926477.1A Pending CN116931006A (en) | 2023-07-26 | 2023-07-26 | Troposphere delay value determining method, system, equipment and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116931006A (en) |
-
2023
- 2023-07-26 CN CN202310926477.1A patent/CN116931006A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109543353B (en) | Three-dimensional water vapor inversion method, device, equipment and computer readable storage medium | |
CN114518586B (en) | GNSS precise single-point positioning method based on spherical harmonic expansion | |
CN103323888B (en) | Method for eliminating delay errors of troposphere of GNSS atmospheric probing data | |
Ingleby et al. | Quality control of ocean temperature and salinity profiles—Historical and real-time data | |
Chen et al. | A comprehensive evaluation and analysis of the performance of multiple tropospheric models in China region | |
CN104965207B (en) | A kind of acquisition methods of zone convection layer zenith delay | |
US7509838B2 (en) | Correction of humidity measurement results of a radiosonde | |
Liu et al. | An analysis of GPT2/GPT2w+ Saastamoinen models for estimating zenith tropospheric delay over Asian area | |
CN108828690A (en) | A kind of method and device of Retrieval of Atmospheric Water Vapor content | |
CN109117555A (en) | A kind of Bevis model refinement method of sub-region right mean temperature | |
CN106324620A (en) | Tropospheric zenith delay method based not on real-time measurement of surface meteorological data | |
CN108680268A (en) | A kind of Bevis model refinement methods of the sub-region right mean temperature based on sounding data | |
CN109145344A (en) | A kind of experience ZTD model refinement method based on sounding data | |
CN105509770B (en) | A kind of barometer on-line correction method in GNSS and MEMS integrated navigation systems | |
CN113960635B (en) | Troposphere delay correction method considering daily variation | |
Yang et al. | Determination of weighted mean temperature (Tm) lapse rate and assessment of its impact on Tm calculation | |
Yao et al. | A global empirical model for estimating zenith tropospheric delay | |
Zhu et al. | A coalescent grid model of weighted mean temperature for China region based on feedforward neural network algorithm | |
Rózsa | Modelling tropospheric delays using the global surface meteorological parameter model GPT2 | |
CN116794682A (en) | Water vapor real-time inversion method and device, electronic equipment and storage medium | |
CN116931006A (en) | Troposphere delay value determining method, system, equipment and storage medium | |
Yao et al. | Establishment of a regional precipitable water vapor model based on the combination of GNSS and ECMWF data | |
CN116165728A (en) | Method and device for estimating atmospheric precipitation, computer equipment and storage medium | |
CN115755103A (en) | Robust self-adaptive GNSS (Global navigation satellite System) water vapor chromatography method | |
Farah | Accuracy assessment study of UNB3m neutral atmosphere model for global tropospheric delay mitigation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |