124 results about "Zenith" patented technology

The invention relates to an optical satellite absolute radiometric calibration method, device and equipment under ocean low brightness based on an offshore pseudo-invariant field and a medium, and belongs to the technical field of in-orbit absolute radiometric calibration. According to the method, a marine pseudo-invariant field with stable water and atmospheric optical characteristics is determined according to marine observation data of a calibrated optical satellite, the off-water radiance of the marine pseudo-invariant field is calculated according to the chlorophyll a concentration of the marine pseudo-invariant field, the zenith radiance is theoretically calculated, a data source is provided for absolute radiometric calibration of a satellite to be calibrated, and calibration is completed. The method explores the potential of the marine pseudo-invariant field in supplementing the traditional radiometric calibration method, gets rid of the limitation of a fixed and single marine optical buoy MOBY site data source, improves the data volume, the coverage area and the acquisition frequency of effective calibration data, provides guarantee for the continuity of an optical satellite observation task, and has a wide application prospect. And a plurality of calibration sources can be utilized to realize rapid and accurate calibration of the optical satellite, so that the calibration efficiency of the offshore optical satellite is remarkably improved.

The invention relates to the technical field of satellite navigation and positioning, in particular to a zenith troposphere wet delay calculation method based on deep learning, and the basic steps of segmented ZWD multi-source data fusion based on an improved Transform model are as follows: establishing a ZTD model, establishing a TCA-Transform model, and fusing data. The task of establishing the ZTD model is to calculate the atmospheric refractive index according to meteorological data so as to obtain ZTDPP; the TCA-Transform model establishment comprises the following steps of: improving a Transform model by using a domain self-adaptive method; according to the data fusion, the ZWDPP and the ZWDMODIS are fused by using a TCA-Transform model. And finally, a ZWD product with high temporal-spatial resolution and high precision can be obtained. The generalization ability of an original Transform model is improved, so that the troposphere delay correction ability in a specific area is improved, and a zenith troposphere delay (ZWD) product with high resolution and high precision is obtained.

The invention relates to a polarization three-dimensional reconstruction method for solving zenith angle based on mixed reflected light. The method comprises the following steps: obtaining polarization images of a target in four directions of 0 degree, 45 degrees, 90 degrees and 135 degrees; respectively calculating the polarization degrees of the polarization images; obtaining a hybrid model of target surface reflected light, solving a zenith angle and an azimuth angle of the target, and obtaining a normal vector field of the target; based on the normal vector field of the target, the Frankot-Chellappa algorithm is used for integrating the normal vector field to obtain the three-dimensional depth map of the target, the method is suitable for outdoor natural light conditions, especially for three-dimensional reconstruction of offshore and long-distance targets, the limitation of a traditional method in a severe environment is overcome, and through efficient processing of the multi-angle polarization image, the depth map of the target is obtained. Dependence on a single light source and strict experiment conditions is reduced, and higher operation convenience and engineering practical value are achieved.

The invention belongs to the technical field of remote sensing application, and provides a grassland biomass inversion method based on a radiation transfer model and machine learning, and the method comprises the following steps: obtaining and processing remote sensing data, DEM data and actually measured biomass data of a target region; performing terrain correction by adopting an SCS + C method; performing global sensitivity analysis on input parameters in the PROSAIL model, and calibrating the input parameters according to the global sensitivity analysis; introducing topographic factors into the model, and correcting the solar zenith angle and the observation zenith angle; generating simulated canopy reflectivity data by using the model, and constructing a data set; constructing a lookup table LUT by combining the simulated canopy reflectivity data and the above-ground biomass data, and performing machine learning modeling according to the canopy spectral data set; and verifying the inversion precision, and outputting a grassland above-ground biomass spatial distribution map. According to the inversion method, terrain factors are introduced on the basis of a traditional model, and the adaptability and estimation precision of the model to a complex terrain area are improved by combining a machine learning method.

The invention relates to a sun meridian detection and course angle real-time calculation method based on polarized light field vector characteristics, and aims to solve the problems that a traditional sun meridian detection method is easily interfered in low signal-to-noise ratio scenes such as haze and cloud layer shielding, lens distortion causes polarization angle errors, the false detection rate is high, and the calculation accuracy is high. In order to solve the problems that the real-time performance and robustness of course calculation are affected due to the fact that single polarization angle information is difficult to fully reflect the sun direction, polarization parameter calculation is carried out by utilizing a multi-channel polarization image, and the sun meridian is detected by adopting Hough transform through fisheye lens correction, bilateral filtering and candidate region screening, so that the real-time performance and robustness of course calculation are improved. And a multi-point inversion method is combined to accurately obtain a sun vector and a zenith angle and an azimuth angle thereof, and finally a course angle is calculated in real time. The problems of imaging distortion and insufficient real-time processing in the traditional sun positioning technology can be effectively solved, and a powerful navigation positioning technical guarantee is provided for the fields of unmanned aerial vehicles, autonomous navigation, polarized light navigation and the like.

The invention belongs to the crossing field of space geodetics, satellite navigation positioning technology and atmospheric science, and discloses a real-time troposphere delay modeling method for enhancing a global forecasting system. The method comprises the following specific steps: firstly, acquiring zenith fluid static delay and zenith wet delay predicted by earth surface GFS, and constructing a regional grid difference time sequence by taking a fusion result of ERA5 data and a GNSS troposphere delay product as a reference; secondly, fitting the difference to form a GFS enhanced core model; meanwhile, a troposphere delay correction model from the grid height to the target site height is constructed based on ERA5 stratified atmosphere data, time sequence modeling is carried out on model coefficients, and finally refined troposphere delay output at different heights is achieved. The method has the advantages that (1) the real-time troposphere model is high in precision; (2) real-time application can keep robust performance without depending on a GNSS network; (3) expansion and deployment are facilitated; and (4) the real-time performance is strong, and seamless integration precision positioning application can be realized.

The invention discloses a foundation space ionosphere and troposphere integrated joint detection method and device, and belongs to the technical field of satellite navigation, and the method comprises the steps: detecting the vertical total electron content of an ionosphere: extracting the inclined total electron content of the ionosphere at an ionosphere puncture point, and converting the inclined total electron content into the vertical total electron content of the ionosphere; detecting a foundation Beidou ionosphere disturbance index; and detecting an ionosphere phase scintillation index, a troposphere zenith total delay, a troposphere precipitable amount and an ionosphere amplitude scintillation index, and outputting ground-based Beidou space ionosphere and troposphere integrated joint detection parameters. The method effectively improves the reliability, robustness and meteorological perception capability of a navigation positioning system in a complex atmospheric environment, and is especially suitable for special environmental conditions such as severe convection weather, typhoon landing, violent disturbance of an ionized layer and the like.

The invention relates to the field of marine meteorological monitoring, and discloses a GNSS-based offshore troposphere atmospheric parameter profile inversion method, which comprises the following steps of: obtaining actual observation values of total troposphere zenith delay, dry delay and wet delay; obtaining troposphere atmospheric parameter profile statistical characteristics through clustering analysis; data points related to temperature, water vapor pressure and air pressure are constructed by utilizing troposphere atmospheric parameter profile statistical characteristics, and theoretical values of zenith total delay, dry delay and wet delay corresponding to the data points are calculated through forward modeling; and selecting one or more of total troposphere zenith delay, dry delay and wet delay, utilizing an optimization algorithm to invert a troposphere atmospheric parameter profile, performing iteration to obtain an optimal data point, performing data quality control on the temperature, the water vapor pressure, the air pressure and the atmospheric refractive index corresponding to the optimal data point, and outputting the temperature, the water vapor pressure, the air pressure and the atmospheric refractive index. The method disclosed by the invention is low in equipment installation and maintenance cost, is suitable for an offshore observation platform, is reliable in inversion result, and can provide support for offshore troposphere atmosphere research.

The invention provides a troposphere delay correction method and system, electronic equipment and a storage medium, and the method comprises the steps: extracting troposphere delay of a reference station through a non-difference combination observation value method according to the obtained real-time reference station observation data, broadcast ephemeris, ultra-fast orbit and clock correction, carrying out the adjustment of observation values of a plurality of satellites, and carrying out the adjustment of the observation values of a plurality of satellites; the zenith troposphere delay of the base station is obtained; the zenith troposphere level of the base station after normalization processing is interpolated to the position of the moving station through an inverse distance weighted interpolation method, zenith troposphere delay of the moving station is obtained according to a pre-constructed elevation normalization model, and the elevation normalization model is constructed based on ERA5 horizontal atmospheric layer meteorological data. According to the method, the adaptability of the model to complex terrains can be improved, so that the model precision is improved, the precision single-point positioning convergence time of the moving station is effectively shortened, and the precision single-point positioning precision is improved.

The invention relates to the technical field of atmospheric parameter inversion models, provides an atmospheric model construction method and system based on virtual temperature correction, and solves the problem of insufficient dynamic scene adaptability. The method comprises the following steps: acquiring a radiation brightness observation value of a satellite thermal infrared band, synchronously recording a satellite observation zenith angle and an azimuth angle, and collecting earth surface coverage type data, dew point temperature and corresponding height information; and uplink and downlink radiation components are eliminated. A normalized vegetation coverage index is calculated based on satellite red light and near-infrared band reflectivity data, a specific emissivity mapping relation is constructed in combination with a land surface coverage type, and land surface specific emissivity space distribution is inverted. And generating vertical stratified atmosphere data based on the dew point temperature and height, and constructing radiation transmission path correction parameters in combination with the observation angle. And finally, establishing an atmospheric radiation transmission model, and outputting a continuous spatialization virtual temperature field. According to the method, the problem of insufficient dynamic scene adaptability in the prior art is solved, and high-precision inversion of the atmosphere deficiency temperature vertical section is realized.

The present invention discloses a high-precision engineering calculation method and device for the back irradiance of a bifacial photovoltaic module, wherein the calculation method comprises the following steps: (1) inputting relevant parameter information of a photovoltaic power station; (2) simplifying the array into a two-dimensional diagram and calculating the sun position, including the zenith angle and azimuth angle; (3) calculating the ground shadow situation according to the module geometric information and the sun position; (4) calculating the sky viewing factor between rows that can be seen by each ground segment according to the array geometric parameters; calculating the incident angle according to the inclination angle, the module azimuth angle, the solar azimuth angle, and the zenith angle; using the Perez model to decompose the input horizontal direct irradiance and horizontal diffuse irradiance into direct, isotropic diffuse, circumsolar diffuse part, horizon scattering part, and ground reflection part; using the Perez model to calculate the inclined surface diffuse irradiance; (5) calculating the back irradiance after all the above steps are completed. This method simplifies the three-dimensional scene into two dimensions, and simplifies the View Factor method, back shadow analysis, sky view factor algorithm, accurately considers the irradiance received by the ground, and uses the angle-by-angle incidence loss to calculate the irradiance received by the back. It has both high accuracy and high efficiency.

The invention provides a GNSS global ZTD correction method and device considering anisotropy, and the method comprises the steps: dividing the whole world into regular windows with the same size according to a mesh generation algorithm, and obtaining a correction model training sample through a Panguu-Weather numerical weather forecast system; constructing a zenith troposphere delay correction model taking anisotropy into consideration by adopting a BP neural network, taking an annual day, the height, the latitude and the longitude of the current position, the zenith troposphere delay and the height, the latitude and the longitude of the target position as input layer data, and taking the zenith troposphere delay of the target position as output layer data; based on the zenith troposphere delay correction model considering anisotropy, a user can obtain correction information of zenith troposphere delay only by providing the annual day, the height and longitude and latitude of the current position, zenith troposphere delay and the height and longitude and latitude of the target position and retrieving the corresponding model, and the precision of GNSS meteorological monitoring and environment monitoring is improved.

The present invention discloses a method, device and storage medium for calibrating the magnitude of a full-field instrument. The method is as follows: taking a star map, detecting and clustering the entire field of view to obtain a large number of bright spots; extracting observed stars, performing photometry and positioning on all bright spots to obtain the position, gray scale, uncertainty and zenith angle; matching the observed stars with a star catalog, and after successful matching, recording the following information, including the stellar magnitude, position, gray scale and uncertainty corresponding to the observed stars; dividing the full field of view into regions, counting all the stars that have been successfully matched in all frames within the region, fitting them by weighted least squares, and recording the zero points of the instrument magnitude obtained in different regions; performing polynomial interpolation on the zero points of the instrument magnitude in different regions of the full field of view to obtain the zero point correction at different fields of view; repeating the above steps at different time periods to obtain an accurate correction method for the on-orbit instrument magnitude of the full time period. The present invention can calibrate the on-orbit instrument magnitude of optoelectronic devices.

This invention discloses a polarization orientation method based on the antisymmetry of the atmospheric polarization pattern. First, a UNet convolutional neural network is used to perform binary segmentation on the polarization image of the sky acquired by a polarization camera to obtain a de-obstructed polarization image of the sky. The separated sky is then divided into multiple circular rings centered at the zenith, and the corresponding positive and negative feature points within each ring are obtained. The antisymmetry of the atmospheric polarization pattern is used to fit the solar meridian, ultimately completing the calculation of the heading angle. This invention solves the problem of traditional polarization orientation methods failing due to complex weather conditions, especially large-scale obstructions, and effectively improves the robustness of polarization orientation methods under complex weather conditions.

The present application provides a precise photogrammetric control measurement method, device and system. The method includes: obtaining the control point coordinates of a plurality of control points arranged in the surrounding area of the object to be measured; a plurality of measured points are arranged on the object to be measured; by using a first total station instrument installed at a first measurement station, obtaining the first positional relationship between the first total station instrument and the control points, and obtaining the second positional relationship between the first total station instrument and the measured points; the first positional relationship includes the first distance, the first horizontal angle and the first zenith distance between the first total station instrument and the control points, and the second positional relationship includes the second distance, the second horizontal angle and the second zenith distance between the first total station instrument and the measured points; according to the control point coordinates of the control points, the first positional relationship and the second positional relationship, obtaining the measured point coordinates of the measured points, and determining the contour of the object to be measured according to the measured point coordinates. The method of the present application can accurately measure the coordinates of the measured points of the object to be measured.

The application discloses a kind of floating offshore solar radiation measurement data correction method and system, it includes the following steps: obtaining the pitch angle and roll angle of floating offshore solar radiation measurement system;Obtain the unit normal vector of the measuring plane pointing to the sky in the horizontal plane rectangular coordinate system;Calculate the angle between the measuring plane and the horizontal plane;Get horizontal plane solar zenith angle and azimuth angle;Construct the unit vector along the direction of sunlight incidence pointing to the sun;Correct the zenith angle;Calculate the horizontal plane scattered irradiance value;The horizontal plane total irradiance saved by floating offshore solar radiation measurement system is corrected, and the corrected horizontal plane total irradiance is obtained.The measurement data of the floating offshore solar radiation measurement system is corrected by the application, the measurement error caused by the inclination of the measurement platform caused by wind and wave is reduced, and more accurate measurement data is provided for the evaluation of offshore light resources and the planning of offshore photovoltaic power station.

This invention provides a method, apparatus, electronic device, and storage medium for calculating cloud judgment parameters. The method includes: determining the satellite position vector, solar vector, and satellite attitude matrix in the cloud judgment time inertial frame based on the time interval between the current time and the cloud judgment time; calculating the position information of the ground imaging point corresponding to each pixel after the time interval based on the satellite payload installation information, pixel line-of-sight information, and the satellite position vector and satellite attitude matrix in the cloud judgment time inertial frame; and calculating the cloud judgment parameters for each ground imaging point based on the position information of the ground imaging point corresponding to each pixel after the time interval and the solar vector in the cloud judgment time inertial frame. The cloud judgment parameters include at least one of the following: satellite zenith angle, solar zenith angle, geographic latitude and longitude, and relative azimuth angle. This solution can provide cloud judgment parameters for the cloud judgment process of satellite payloads.

The application discloses a kind of sea fog environment multi-coating target multispectral polarization thermal imaging measuring device and method, belong to polarization detection technical field, method includes the initial polarization characteristic information of different coating target acquisition;Get the polarization characteristic information of different incident and detection zenith angle combination, the polarization characteristic information of different incident and detection azimuth angle combination;Collect polarization characteristic information under different temperature;The data collected are classified according to temperature and geometric angle, calculate stokes vector, linear polarization degree and polarization angle, analyze the polarization characteristic difference of coating target under different wave band;Utilize deep neural network to construct SIP-Net fusion network, realize the feature measurement and identification of sea surface coating target.The device includes marine environment simulation unit, lighting unit, detection unit, control and data processing unit, cooperates and is used to realize the method.The application solves the problem that the detection precision of prior art is low under complex sea surface environment, which leads to difficult target identification.

The invention relates to an arbitrary polyhedron spinning space debris space-based observation light variation curve simulation method. The method comprises the following steps: calculating position vectors of debris and an observation satellite in a geocentric inertial system at the current moment according to orbital elements of the debris and the observation satellite; calculating a sun vector and an observation satellite vector under the debris body system according to the attitude of the debris at the current moment; calculating a sun incident zenith angle and a satellite observation zenith angle of each surface of the debris; calculating the actual effective reflection area of each surface of the fragment; calculating the optical reflectivity of each surface of the fragments based on a BRDF scattering model; calculating the total reflection irradiance of each surface of the fragments; calculating to obtain the apparent magnitude of the fragments at the observation satellite according to an apparent magnitude calculation formula; and on the basis of a preset fragment attitude sequence, repeating the steps to calculate and obtain the light variation curve of the fragments. According to the method, the attitude spinning condition of the debris target is fully considered, and a basis can be provided for designing a space-based debris observation satellite remote sensor and inverting the debris attitude by using the light variation curve.

The application discloses an InSAR atmospheric delay phase correction method along a slant range direction, comprising the following steps: S1, acquiring two SAR image subsets covering a research area and acquiring two clean water vapor data; S2, performing differential interference processing on the two SAR image subsets covering the research area to acquire differential interference phases; S3, accurately calculating atmospheric delay phases by using an atmospheric delay phase correction model; and S4, accurately obtaining real deformation of the ground along a radar line-of-sight direction. The application gives the latest weight value of the zenith path delay difference of each ground pixel on the slant range direction of the ground by correcting the sampling mode of the water vapor product along the ground distance direction, finally proposes an accurate and efficient slant range phase total delay calculation method along the ground distance, further optimizes and improves the repeat orbit synthetic aperture radar differential interferometry method, so that the measurement result is more scientific and reasonable, and therefore the application has a wide application prospect in the field of repeat orbit InSAR atmospheric phase correction.

This invention proposes a spaceborne ultra-high-resolution SAR imaging method with real-time tropospheric correction, which meets the requirements of high timeliness and high precision SAR imaging processing without being limited by the imaging area. Specifically, a corrected radiometer capable of covering the SAR observation field of view is used to simultaneously measure atmospheric temperature and water vapor pressure; the measurement data and SAR echo data are simultaneously transmitted to the ground to obtain the wet tropospheric zenith delay; using the total tropospheric delay estimation model, the total tropospheric zenith delay is obtained, and the slant range delay caused by the troposphere is further obtained; the slant range delay is used as a correction term to correct the SAR slant range model; in the ultra-high-resolution SAR imaging processing on the ground, the imaging algorithm adopts the slant range corrected by the real-time tropospheric correction, and finally obtains the ultra-high-resolution SAR image.

The invention discloses an atmospheric background radiation modeling simulation method and device for air optical detection, and belongs to the field of optical simulation. The method comprises the following steps: determining whether a detection scene is a cloud scene or not according to whether the detection scene has a cloud layer or not; if yes, generating a target three-dimensional cloud field by adopting a particle system; according to the initial environment parameters, respectively establishing a first lookup table of the radiance of the clear sky scene about the detector height and the pixel line-of-sight zenith angle, and establishing a second lookup table of the radiance and the transmittance of each segment path of the cloud scene about the detector height and the detection line-of-sight zenith angle; determining air data in all the detected pixel data according to the periodic height and the visual geometric information of the detector; and searching a lookup table of a corresponding scene according to the air data, and determining a radiance value used for generating a radiance image in front of a mouth. According to the invention, complex air optical visual simulation can be realized.

The embodiment of the invention provides a measurement error correction method and device, electronic equipment and a storage medium. The method comprises the following steps: acquiring meteorological data and original data of a global satellite navigation system ground observation station; calculating the atmospheric zenith delay amount through the meteorological data and the original data of the global satellite navigation system ground observation station; and correcting a measurement error caused by atmospheric delay by using the atmospheric zenith retardation. By fusing multi-source data (ERA5 and GNSS), the atmospheric delay error in InSAR measurement is accurately calculated and corrected, so that the precision and reliability of surface deformation monitoring are remarkably improved.

The present invention discloses a tomographic processing method for mining areas based on borehole observation. Utilizing the principle that cosmic ray muons lose intensity after passing through matter, muon detectors are arranged in boreholes to observe muon intensity underground, thereby obtaining information on the distribution of material opacity at various locations underground with respect to zenith angle and azimuth. Through multi-borehole and multi-position depth measurements, a cross-coverage of muon observations of the observation area is formed, thereby establishing a density inversion solution equation, and ultimately solving the accurate three-dimensional density structure under the surface. The muon source used in the present invention is a natural source, and muons do not damage the observed object. The observable depth can reach up to a depth of one kilometer, and the existing boreholes in the mining area are fully utilized to provide a safety detection function, thereby ensuring safe production.

The aim is to suppress the increase in resources and costs that would result from an increase in the number of satellites by reliably observing directional radio waves, which have uncertainty, using a small number of satellites. [Solution] The radio wave monitoring system includes a satellite constellation forming unit that forms a satellite constellation that flies within a pointing change range at an orbital altitude H where directional radio waves may be emitted. The satellite constellation forms a satellite formation that includes a first group of satellites that have the same azimuth angle value and have the same zenith angle difference between adjacent satellites in a tangential plane coordinate system with the monitored object 10 as the origin, at the time when the position coordinates of the monitored object 10 pass through the common range of the observation ranges of each satellite formed when the observation field of view θsat of each satellite reaches the Earth's surface.

The embodiment of the application provides a method and device for calculating data required by a lunar observation task, a specific embodiment of the method comprising: calculating a zenith angle of the earth and the sun respectively at a target time according to a first coordinate and a second coordinate obtained from ephemeris data and longitude and latitude of a target lunar surface position on the lunar surface; determining a first part of the tidal force of the target lunar surface position at the target time according to data for calculating the first part of the tidal force obtained by applying a second-order expansion formula of associated Legendre functions and a third-order expansion formula of associated Legendre functions to the zenith angle of the earth at the target time; determining a second part of the tidal force of the target lunar surface position at the target time according to data for calculating the second part of the tidal force obtained by applying the second-order expansion formula of associated Legendre functions to the zenith angle of the sun at the target time; and calculating the tidal force of the target lunar surface position at the target time according to the first part of the tidal force and the second part of the tidal force.

The embodiment of the invention provides a zenith troposphere delay prediction method, a zenith troposphere delay prediction device and electronic equipment. The method comprises the following steps: acquiring multi-source heterogeneous data; performing alignment on the multi-source heterogeneous data to obtain global grid data; according to the global grid data, zenith statics delay ZHD and zenith wet delay ZWD are modeled to obtain a ZHD long-period prediction model and a ZWD short-time prediction model; according to the global grid data, constructing a correction model which changes along with the elevation in the vertical direction; and according to the ZHD long-period prediction model, the ZWD short-time prediction model and the correction model, performing ZTD prediction on a to-be-measured station. And the prediction accuracy of the ZTD is improved, so that support is provided for improving the stability and precision of GNSS positioning.

The application discloses a star-inertial integrated outfield dynamic calibration method based on a refracting surface constraint mechanism, and relates to the technical field of calibration methods, and specifically discloses the following technical scheme: a relationship model of an observed star position and a zenith direction is established based on a refracting surface constraint mechanism and observation imaging of a star, a first zenith direction measurement model under a star sensor system corresponding to an observation time and an observation position is established based on the relationship model; a second zenith direction measurement model under an inertial navigation system of the zenith direction of a current observation position is established through a filter combination model based on real-time measurement values of angular velocity and angular acceleration of the inertial navigation system; a zenith direction integrated conversion matrix model between the inertial navigation system and the star sensor system is established by simultaneously solving the first zenith direction measurement model under the star sensor coordinate system and the second zenith direction measurement model under the inertial navigation system; and a star-inertial integrated outfield optimization calibration model is established by changing the longitude and latitude of the observation position and the observation direction for multiple times, so that star-inertial integrated outfield calibration is carried out.

The present invention belongs to the technical field of space geodesy, and particularly relates to a method for calculating the zenith tropospheric delay in a region. The method first selects GNSS stations that meet the accuracy requirements for calculating the zenith tropospheric delay in the target region according to the corresponding relationship between the calculation accuracy of the zenith tropospheric delay and the shortest distance of the stations; wherein, the corresponding relationship is obtained by performing simulation calculations on the target region; then, according to the observables of the selected GNSS stations, the zenith tropospheric delay values at each GNSS station are calculated; the zenith tropospheric delays at each GNSS station are extended to other reference surfaces, and after the extension, grid processing is performed to obtain the regional delay values on other reference surfaces; finally, the regional delay values on other reference surfaces are inversely extended to the ground surface to obtain the regional zenith tropospheric delay values. The present invention can obtain accurate calculation results of the zenith tropospheric delay in the target region based on a small number of GNSS stations.

The application relates to the meteorological detection technical field and provides a water vapor density determination method and device and a computer readable storage medium. The method comprises the following steps: obtaining target information of each observation device in a plurality of observation devices; the target information comprises first information and air pressure data of an air pressure sensor; determining troposphere delay of a GNSS chip according to the first information; determining zenith troposphere wet delay corresponding to the observation device according to the troposphere delay of the GNSS chip and the air pressure data; taking the product of the zenith troposphere delay and a troposphere wet delay projection function as ray path troposphere wet delay corresponding to the observation device; and determining water vapor density of a target area according to the ray path troposphere wet delay corresponding to the plurality of observation devices. Through the scheme, the cost of water vapor density determination can be reduced, and the accuracy of water vapor density determination can be improved.