96 results about "Zenith" patented technology

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 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.

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 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 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 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 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.

The present invention discloses an angle effect correction method for geostationary meteorological satellite thermal infrared surface temperature data, which belongs to the field of thermal infrared remote sensing technology. The present invention first obtains geostationary satellite thermal infrared remote sensing surface temperature data of the target area and converts it to the target coordinate system; then calculates the angle parameters of each pixel based on the satellite sub-satellite data; performs system error elimination processing on the nighttime surface temperature data of two satellites screened based on the satellite zenith angle; then uses the nighttime surface temperature data to simplify the kernel-driven model and solve the emissivity kernel parameters; solves the solar nuclear coefficient and hotspot height based on the daytime surface temperature data at multiple consecutive moments, and combines the solved parameters to obtain vertical observations of each pixel. The present invention can quickly and accurately solve the parameters of the kernel-driven model and then accurately solve the vertically observed surface temperature, overcoming the defects of the traditional method such as large data demand, low calculation efficiency, and low accuracy, and can serve a large number of thermal infrared remote sensing surface temperature application scenarios.

The invention provides a dam GNSS deformation monitoring method and system based on a water vapor three-dimensional model, and relates to the technical field of dam safety monitoring, and the method comprises the steps: obtaining GNSS observation data of a plurality of reference stations in a monitoring region, meteorological station data, and remote sensing data containing regional terrain elevation information; then calculating a zenith total delay sequence based on GNSS observation data, and separating target zenith wet delay; then calculating water vapor content horizontal spatial distribution based on target zenith moisture delay, obtaining vertical spatial distribution in combination with regional terrain elevation information, and constructing a water vapor three-dimensional model; then, water vapor delay correction is conducted on GNSS observation data through the model, and a displacement time sequence is obtained through time sequence analysis; and finally, calculating a long-term deformation rate and identifying an abnormal deformation event, and generating a visual analysis result and a standardized report after quality control. The accuracy and reliability of dam vertical deformation monitoring are improved.

This application discloses a new generation BDSBAS time synchronization method and system based on a layer-by-layer nearest neighbor non-repetition strategy, involving satellite communication and data processing technologies. The method includes: obtaining the number of BeiDou navigation satellites observable by the main monitoring station k at the current time; using the least squares method to obtain an estimate and accuracy of the clock bias of the main monitoring station; and based on the estimate and accuracy of the clock bias of the main monitoring station, using the common-view satellite method to synchronize the clocks of other ground monitoring stations. i Time synchronization is performed between the ground-based main monitoring station k and the observable low-Earth orbit (LEO) satellites in the zenith direction using the common-view satellite method. This achieves time synchronization by calculating the time deviation between the ground-based main monitoring station k and the observable LEO satellites in the zenith direction. A layer-by-layer nearest neighbor (NNN) time synchronization strategy is employed to synchronize the LEO satellites with their respective nearest neighbors. This application solves the problem of ground-to-ground time synchronization applicable to the next-generation BDSBAS, achieving synchronization between the clocks of the entire LEO constellation and the clocks of the ground-based main monitoring station.

The invention relates to the technical field of GNSS positioning, in particular to a zenith troposphere delay estimation method and device and a medium, and the method comprises the following steps: constructing a zenith dry delay calculation module for calculating a GNSS site; selecting a combination of a correlation model and a regression model which are most suitable for the current GNSS data set; determining a spatial distribution scheme of an optimal GNSS station; carrying out zenith troposphere delay spatial interpolation modeling and obtaining an interpolation result by adopting the selected combination of the correlation model and the regression model and a spatial distribution scheme, and constructing a regional zenith wet delay interpolation model; and acquiring a troposphere delay estimated value by using the regional zenith wet delay interpolation model. According to the method, the elevation influence of zenith troposphere delay space interpolation modeling is removed, and the GNSS positioning precision is improved.

This application discloses a system, method, and apparatus for frequency discrimination curve calibration. Frequency locking technology is employed to stably lock the output laser frequency of the seed laser onto the absorption peak of the atomic / molecular absorption cell, effectively solving the problem of seed laser frequency drift. Simultaneously, two acousto-optic modulators are used to perform precise extracavity passive frequency modulation of the seed laser, shifting the laser frequency from the atomic / molecular absorption peak position to the preset zero-velocity frequency point and the edge frequency point of the linear operating region of the optical frequency discriminator. During calibration, the telescope is vertically pointed towards the zenith, providing atmospheric echoes with zero Doppler shift as a reference. At this time, the optical frequency discriminator measures and records the transmittance at the preset zero-velocity frequency point at the center of the linear operating region and at the edge frequency point. Based on the linear assumption of the frequency discrimination curve, a new frequency discrimination curve can be derived.

The present invention relates to a method, device, equipment and medium for absolute radiation calibration of optical satellites under low ocean brightness based on a pseudo-invariant field at sea, and belongs to the technical field of on-orbit absolute radiation calibration. The method determines a pseudo-invariant field at sea with stable water and atmospheric optical characteristics from calibrated optical satellite ocean observation data, calculates its water-leaving radiance through the chlorophyll a concentration of the pseudo-invariant field at sea, and theoretically calculates the zenith radiance, providing a data source for absolute radiation calibration of the satellite to be calibrated, thereby completing the calibration. The present invention explores the potential of pseudo-invariant fields at sea in supplementing traditional radiation calibration methods, breaks away from the data source limitations of fixed, single ocean optical buoy MOBY sites, increases the data volume, coverage area and acquisition frequency of effective calibration data, provides a guarantee for the continuity of optical satellite observation missions, and can use multiple calibration sources to achieve rapid and accurate calibration of optical satellites, significantly improving the efficiency of marine optical satellite calibration.

The invention discloses a cloud scene spectral radiance simulation method based on space-based observation and a related device, and relates to the technical field of spectral radiance simulation, and the method comprises the steps: obtaining simulation parameters of a cloud-containing layer, atmosphere, a position and a wave band; calculating a cloud layer optical characteristic coefficient by using a Mie scattering theory, and solving a cloud layer refraction coefficient in combination with a Kramers-Kronig relationship; based on the refraction coefficients of the atmosphere and the cloud layer, the cloud layer surface reflectivity is obtained through the Fresnel reflection law; and solving a radiation transfer equation by adopting a two-flow approximation method, enabling the non-uniform atmosphere to be equivalent to a plurality of uniform sub-layers, introducing reflectivity into the cloud sub-layer to represent the reflection contribution of the cloud sub-layer, and finally obtaining the spectral radiance at the zenith. According to the method, parameters are accurately set to adapt to various scenes, the calculation precision of optical parameters and reflectivity is ensured by means of a physical mechanism, the calculation burden is remarkably reduced through two-stream approximation, efficiency and precision are considered, and the problems that an existing method is time-consuming in calculation and large in burden on complex clouds, especially clouds with high optical thickness are effectively solved; and the spectral radiance of the cloud scene can be accurately simulated.

The application discloses a troposphere delay correction method fusing GNSS and numerical weather prediction products, comprising the following steps: when receiving a delay calculation request of a target point at a target time, determining a grid unit where the target point is located; inputting the coordinates and elevation of the target point into a fusion model constructed in the grid unit to output predicted values of zenith dry delay and zenith wet delay; for the grid unit with a GNSS site, the construction of the fusion model comprises the following steps: obtaining first delay data of the GNSS site at the target time through time extrapolation based on a sliding autoregressive model, synchronously obtaining second delay data of a numerical weather prediction product at the same target time, and aligning the first delay data and the second delay data with the target time; adopting a double dynamic weighting strategy of distance and accuracy to fuse the time-aligned first delay data and the second delay data on zenith dry delay and zenith wet delay respectively, and outputting predicted values of the target point at the target time.

The present disclosure provides an underwater point source polarized light detection device and a heading angle measurement method, which are applied to the field of polarization navigation technologies. The polarized light detection device includes a base, a rotating detection platform provided on the base, and three sets of polarized light detection modules provided on top of the rotating detection platform. The three sets of polarized light detection modules are of the same structure, each set of detection module consists of a pair of polarization analyzing channels with transmission directions perpendicular to each other, the other two sets of detection modules are positioned such that polarized light detection directions thereof are bilaterally symmetric in the same plane about a zenith polarized light detection direction and that an included angle between each of the polarized light detection directions and the zenith polarized light detection direction does not exceed 5°, and the rotating detection platform is configured to be rotatable for 360° about an axis of the zenith polarized light detection direction. The heading angle measurement method using the detection device has advantages of good adaptability to underwater environments and requiring a reduced number and types of polarization analyzing channels.

This paper proposes a water body identification and extraction method and system based on bilateral red band GF6-WFV data, including: preprocessing GF6 satellite WFV image data to obtain atmospherically corrected multi-angle observation data for each satellite band; calculating the reflectivity of each satellite band, normalized to the same zenith angle after atmospheric correction, to obtain the angularly normalized reflectivity of each band; superimposing bands with high and low reflectivity effects and then normalizing them to eliminate the noise impact of areas including shadows and dark surfaces; calculating a targeted high-resolution water index and the normalized difference vegetation index (NDVI). This paper reveals the differences in sensitivity of different bands in water body monitoring and can effectively distinguish water bodies from other land cover types in complex natural and urban environments, making the water body identification process more accurate.

Disclosed is a method for joint estimation of stellar atmospheric refraction and star tracker attitude, including: capturing, by a star tracker, an image of stars, and recognizing the image to obtain a matching relationship between an observed star and a reference star; transforming the reference star to a terrestrial reference frame based on time and position of observation to obtain a true zenith distance of the reference star; calculating an estimated stellar atmospheric refraction based on a simplified atmospheric refraction model, the true zenith distance, and an initial atmospheric refractive coefficient, and subjecting the reference star to stellar atmospheric refraction compensation; calculating the star tracker attitude, and re-projecting the observed star to the terrestrial reference frame to calculate the observed stellar atmospheric refraction and the stellar atmospheric refraction error; if the observed stellar atmospheric refraction is misaligned with the estimated stellar atmospheric refraction, adjusting the atmospheric refractive coefficient based on the stellar atmospheric refraction error, and then recompensating the reference star and calculating the attitude, till the observed stellar atmospheric refraction is aligned with the estimated stellar atmospheric refraction, whereby joint estimation results of the stellar atmospheric refraction and the star tracker attitude are obtained. The disclosure realizes real-time, autonomous estimation and cancelation of stellar atmospheric refraction without being limited by external sensors or empirical formulae.