CN112884342A - Water color satellite atmospheric layer top radiation product quality evaluation and cross calibration method - Google Patents
Water color satellite atmospheric layer top radiation product quality evaluation and cross calibration method Download PDFInfo
- Publication number
- CN112884342A CN112884342A CN202110257688.1A CN202110257688A CN112884342A CN 112884342 A CN112884342 A CN 112884342A CN 202110257688 A CN202110257688 A CN 202110257688A CN 112884342 A CN112884342 A CN 112884342A
- Authority
- CN
- China
- Prior art keywords
- radiation
- satellite
- quality
- product
- cross
- 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.)
- Granted
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 177
- 238000000034 method Methods 0.000 title claims abstract description 42
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 28
- 238000013441 quality evaluation Methods 0.000 title claims abstract description 25
- 230000002159 abnormal effect Effects 0.000 claims abstract description 13
- 230000008569 process Effects 0.000 claims abstract description 10
- 238000011156 evaluation Methods 0.000 claims abstract description 8
- 230000003746 surface roughness Effects 0.000 claims description 13
- 238000004458 analytical method Methods 0.000 claims description 11
- 238000004088 simulation Methods 0.000 claims description 9
- 230000001419 dependent effect Effects 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 238000007405 data analysis Methods 0.000 claims description 3
- 238000000611 regression analysis Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 230000005856 abnormality Effects 0.000 abstract description 9
- 238000012544 monitoring process Methods 0.000 abstract description 9
- 230000000694 effects Effects 0.000 abstract description 2
- 230000007246 mechanism Effects 0.000 abstract description 2
- 230000005540 biological transmission Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 4
- 241000282414 Homo sapiens Species 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000005445 natural material Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Educational Administration (AREA)
- Mathematical Optimization (AREA)
- Development Economics (AREA)
- Computational Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Marketing (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a quality evaluation and cross calibration method for a water color satellite atmospheric layer top radiation product, which comprises the steps of establishing a multi-source satellite quasi-synchronous effective pixel selection constraint condition; the radiation signal abnormality caused by natural process is found to have stronger interstellar correlation; designing an error analytical equation of a multi-source satellite radiation product; and combining a dynamic cross radiometric calibration mechanism; the effect of the quality of the radiation product on the water color product is determined. The scheme firstly standardizes the constraint conditions of quasi-synchronous observation, determines the analytical equation of the error of the multi-source satellite radiation product, realizes the dynamic evaluation and cross calibration of the quality of the water body radiation product, overcomes the serious influence of abnormal signals caused by natural processes in the traditional method, and improves the quality of the water body radiation product and the automatic service level; the method is easy to realize, can be used for dynamically monitoring and improving the quality of the water body radiation product, and has wider practical application value and economic value.
Description
Technical Field
The invention belongs to the field of quantitative remote sensing, and particularly relates to a quality evaluation and cross calibration method for a water color satellite atmospheric layer top radiation product.
Background
The coastal lakes, the coastal rivers, the seas and the like are regions with the most potential and vitality for the prosperous development of the human society, are important ecological environment bearing areas and have great economic benefits. As the population is continuously gathered towards the areas such as lakes, rivers, seas and the like in recent years, the water ecology is more and more stressed, and the problems of resources and environment are more and more serious. Strengthening the monitoring and management of water environments in regions such as lakes, rivers, seas and the like has become an important problem to be solved urgently by governments at all levels. The remote sensing data has the advantages of wide coverage area, diversity of spatial scale, continuous time scale, rich spectral information, flexible and convenient observation and the like, and can provide an effective means for dynamically monitoring the resource environmental conditions of lakes, rivers and seas. Since the Landsat satellite was successfully launched in 1972, the remote sensing images collected by thousands of satellites in the life cycle are accumulated by human beings, and the remote sensing images are greatly successful in the aspects of disaster prevention and reduction, forecast and early warning, resource investigation and the like, so that the living environment and the development space of the human beings are protected, and beautiful Chinese construction is assisted.
The water body water-leaving radiation belongs to a typical weak signal, so that the requirement of water environment remote sensing on satellite image quality is relatively high. The dynamic change of the gravity environment and the like is necessary to strengthen the dynamic monitoring of the quality of satellite radiation products so as to guide the development of radiation calibration work. However, most satellites do not have equipment for detecting the quality of the radiation product, and the quality evaluation of the radiation product is usually realized by using image statistics. However, besides the radiation signal abnormality caused by data error, the radiation signal abnormality is caused by natural factors such as atmospheric conditions, observation geometry, sea water characteristics and the like, and the reliability of the traditional method is greatly reduced. In addition, strict quasi-synchronous effective pixel selection constraint conditions are lacked, high-precision radiation products are difficult to obtain through traditional cross calibration, and business application of the water radiation products is seriously influenced.
Disclosure of Invention
The invention provides a quality evaluation and cross calibration method for a water color satellite atmospheric layer top radiation product, aiming at overcoming the defects in the aspects of dynamic monitoring and cross calibration of the quality of the traditional water body radiation product, and dynamically realizing the cross calibration of the satellite radiation product so as to improve the satellite image data quality and better serve the monitoring of water environments such as lakes, rivers, seas and the like.
The invention is realized by adopting the following technical scheme: a quality evaluation and cross calibration method for a water color satellite atmospheric layer top radiation product comprises the following steps:
step S1, establishing a multisource satellite quasi-synchronization effective pixel selection constraint condition, and extracting an effective pixel metadata set:
through radiation transmission simulation and quasi-synchronous satellite data analysis, constraint conditions for judging illumination-observation geometry, hydrological weather, transit time delay, sea surface roughness and bright targets of quasi-synchronous observation are quantitatively established and used for judging the synchronism of multi-source satellite observation data; extracting an effective image metadata set required by quality evaluation and cross calibration of the radiation product from the data of the multi-source satellite radiation product under the constraint condition;
step S2, constructing an error analytic equation: establishing an error analytical equation of a multi-source satellite radiation product by analyzing consistency characteristics of signal abnormity caused by natural changes of atmosphere and marine substances detected by different satellites in a quasi-synchronous manner;
step S3, quality evaluation: extracting error information of the multi-source satellite radiation product according to the effective pixel data set and an error analytic equation of the multi-source satellite radiation product, and performing dynamic evaluation on the quality of the radiation product;
if the error information of the radiation product obtained by the error equation meets the quality requirement, judging that the quality of the radiation product is reasonable; if the quality requirement range is exceeded, executing step S4;
step S4, cross scaling: and taking the effective image metadata set under the constraint condition as input, and carrying out dynamic cross calibration correction on the satellite radiation product.
Further, in step S1, when determining the illumination-observation geometric constraint condition, the illumination-observation geometry is represented by a solar zenith angle, a satellite zenith angle, and a relative azimuth angle, and the maximum value and the minimum value of the difference in the pixel illumination-observation geometry allowed by the quasi-synchronous observation of the multi-source satellite are obtained by using radiation transmission simulation on the premise of satisfying the IOCCG quality target, and are used as the illumination-observation geometric constraint condition for determining the quasi-synchronous observation.
Further, in step S1, when determining the constraint condition of the hydrographic weather, based on the hydrographic weather data, the relationship diagram of the severe hydrographic weather and the interplanetary deviation of the radiation product is drawn by combining the quasi-synchrotron radiation product of the multi-source satellite, where the severe hydrographic weather includes wind speed and atmospheric pressure, and by analyzing the statistical characteristics of the interplanetary deviation of the two satellite radiation products under different severe hydrographic weather conditions, on the premise of meeting the IOCCG quality target, the upper limit of the severe hydrographic weather is determined and used as the constraint condition of the hydrographic weather for distinguishing the quasi-synchrotron observation.
Further, in step S1, when determining the sea surface roughness constraint condition, the maximum spatial variation coefficient of the satellite image is determined as the constraint condition for determining the sea surface roughness of the quasi-synchronous observation on the premise of meeting the IOCCG quality target by analyzing the quantitative relationship between the spatial variation coefficient of the satellite image radiation product and the interstellar bias of the radiation product.
Further, in step S1, when the transit time delay is determined, the maximum transit time delay of the multi-source satellite is determined by analyzing the quantitative relationship between the satellite transit time difference and the interstellar bias of the image radiation product, and the maximum transit time delay is determined as a constraint condition for determining the transit time delay of the quasi-synchronous observation on the premise that the IOCCG quality target is satisfied.
Further, in step S1, when determining the constraint condition of the bright target, the minimum radiation value of all bright targets in the image is determined as the constraint condition by counting the histogram of the radiation signals of the bright targets in the satellite image.
Further, when the error resolution equation is established in step S2, the following method is specifically used:
s21, selecting at least three transit satellites, and acquiring an effective image metadata set required by an error analysis equation based on the constraint conditions in the step S1;
s22, segmenting the satellite image, establishing an N multiplied by N pixel local area sample set, and calculating the deviation between the pixel radiation quantity and the mean value in the local area as a local abnormal signal;
s23, converting the local abnormal signals into interstellar deviation, and eliminating the influence of the abnormal signals caused by natural processes on an error analytic equation;
and S24, establishing an error analysis equation of the multi-source satellite radiation product by calculating the variance of interstellar deviation by using the assumption that the errors of the multi-source satellite radiation product are mutually independent, and carrying out dynamic monitoring on the quality of the radiation product.
Further, in step S4, when the error of a certain satellite radiation product obtained by the error analysis equation is too large, cross calibration is performed, specifically, the radiation data of the reference satellite is used as a dependent variable, the radiation data of the target satellite is used as an independent variable, an empirical function relationship between the dependent variable and the independent variable is established based on regression analysis, and the empirical function relationship is applied to the radiation data of the target satellite, so as to obtain the cross-calibrated radiation product.
Compared with the prior art, the invention has the advantages and positive effects that:
the method comprises the steps of establishing a multisource satellite quasi-synchronization effective pixel selection constraint condition by analyzing the influence of illumination-observation geometry, hydrological weather, transit time delay, sea surface roughness and the like on radiation signal stability; by comparing the correlation characteristics of the quasi-synchronous satellite signal abnormality, the radiation signal abnormality caused by the natural process is determined to have strong interstellar correlation; by analyzing the statistical characteristics of the data errors of the radiation products, an error analytical equation of the multi-source satellite radiation products is designed, and the influence of the quality of the radiation products on the water color products is determined by combining a dynamic cross radiometric calibration mechanism;
the scheme firstly standardizes the constraint conditions of quasi-synchronous observation, determines the analytical equation of the error of the multi-source satellite radiation product, realizes the dynamic evaluation and cross calibration of the quality of the water body radiation product, overcomes the serious influence of abnormal signals caused by natural processes in the traditional method, and improves the quality of the water body radiation product and the automatic service level; the method is easy to realize, can be used for dynamically monitoring and improving the quality of the water body radiation product, and has higher practical application and popularization values.
Drawings
FIG. 1 is a flowchart of the quality evaluation and cross-calibration of a radiation product according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of the interstellar correlation of the anomaly signals caused by natural processes according to the embodiment of the present invention, wherein (a) and (b) represent the spatial anomaly signals of VIIRS and MODIS satellite image radiation products, respectively, and (c) is a histogram of (a) and (b) images;
FIG. 3 is a schematic diagram of the reliability of an error analysis equation for simulation verification of radiation transmission according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of the effect of radiometric calibration errors on an error resolution equation according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of a cross-scaling model according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of the evaluation of the interstellar consistency of the watercolor product after cross-calibration correction according to the embodiment of the invention.
Detailed Description
In order to make the above objects, features and advantages of the present invention more clearly understood, the present invention will be further described with reference to the accompanying drawings and examples. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, however, the present invention may be practiced in other ways than those described herein, and thus, the present invention is not limited to the specific embodiments disclosed below.
In an embodiment, the present invention provides a method for evaluating quality of a top radiation product of an atmospheric layer of a water color satellite and performing cross calibration, as shown in fig. 1, the method includes the following steps:
step S1, establishing a multisource satellite quasi-synchronization effective pixel selection constraint condition, and extracting an effective pixel metadata set:
the method comprises the steps of taking the quality requirement of an IOCCG product as target guidance, quantitatively establishing constraint conditions for judging illumination-observation geometry, hydrological weather, transit time delay, sea surface roughness and bright targets of quasi-synchronous observation through radiation transmission simulation and quasi-synchronous satellite data analysis, and judging the synchronism of multi-source satellite observation data; extracting an effective image metadata set required by quality evaluation and cross calibration of the radiation product from the data of the multi-source satellite radiation product under the constraint condition;
step S2, by analyzing the consistency characteristics of signal abnormity caused by natural change of atmosphere and marine substances detected by different satellites in quasi-synchronization, eliminating the signal abnormity caused by natural change by adopting a method of subtracting total abnormal signals of different satellites, and establishing an error analytic equation of a multi-source satellite radiation product;
step S3, extracting error information of the multi-source satellite radiation product according to the effective pixel data set and the error analytic equation of the multi-source satellite radiation product, and performing dynamic evaluation on the quality of the radiation product;
if the error information of the radiation product obtained by the error equation meets the quality requirement, judging that the quality of the radiation product is reasonable; if the quality requirement range is exceeded, executing step S4;
and step S4, taking the effective image metadata set under the constraint condition as input, and establishing a cross calibration model based on a regression fitting method to realize dynamic cross calibration correction of the satellite radiation product.
In fact, the better the radiation data synchronism, the more beneficial the cross calibration modeling and the correct evaluation of the data quality. However, under the influence of factors such as illumination-observation geometry, hydrological weather, transit time delay, sea surface roughness and the like, great deviation exists between multi-source satellite radiation data, and based on the consideration, a set of scientific and strict constraint conditions is provided in the embodiment, and radiation data with large interstellar deviation is rejected. In this embodiment, when the effective pixel selection constraint condition involved in step S1 is that the maximum or minimum value of the illumination-observation geometry (including the solar zenith angle, the satellite zenith angle and the relative azimuth), the hydrological weather, the transit time delay and the sea surface roughness, and the radiation amount of the bright pixel is satisfied when the IOCCG watercolor satellite data quality requirement is satisfied (for example, the radiation signal error of 443nm band is not more than 2.5%), step S1 is specifically implemented by the following method:
(1) by utilizing radiation transmission simulation, on the premise of meeting the IOCCG quality target, the maximum value and the minimum value of the difference in the pixel illumination-observation geometry allowed by the multisource satellite quasi-synchronous observation are obtained and used as illumination-observation geometry constraint conditions for judging the quasi-synchronous observation:
the illumination-observation geometry is represented by a sun zenith angle, a satellite zenith angle and a relative azimuth angle, the relative azimuth angle and the satellite zenith angle are fixed and unchanged, the sun zenith angle is kept to be dynamically changed within a range of 0-60 degrees, contour maps of radiation errors caused by interstellar deviation of different sun zenith angles and sun zenith angles are drawn, a distribution area of the sun zenith angle and sun zenith angle interstellar deviation corresponding to the radiation product error of <2.5 percent is determined, and the upper, lower, left and right boundaries of a rectangle with the largest built-in area of <2.5 percent are used as the respective minimum value and the maximum value of the sun zenith angle and sun zenith angle interstellar deviation. The maximum allowable dynamic range (maximum and minimum) of the satellite zenith angle and satellite zenith angle deviation and the relative azimuth angle and relative azimuth angle deviation can be obtained in the same way.
Illumination-observation geometric constraint condition A0Specifically, the value is based on the actual situation, VIIRS and FY-3D satellites are taken as columns, and the obtained constraint conditions of illumination-observation geometry selected by the effective pixels are as follows: the solar zenith angle is larger than 19 degrees but smaller than 50 degrees, the solar zenith angle interstellar deviation is not larger than 2.9 degrees, the satellite zenith angle is larger than 20 degrees but smaller than 38 degrees, the satellite zenith angle interstellar deviation is not larger than 4.9 degrees, the relative azimuth angle is larger than 110 degrees but smaller than 150 degrees, and the relative azimuth angle interstellar deviation is not larger than 15 degrees.
(2) Under the same hydrological meteorological condition, discussing the statistical characteristics of satellite image interstellar deviation, establishing the quantitative linkage relationship between the two, determining the upper limit of severe hydrological meteorological under the premise of meeting the IOCCG quality target, and taking the upper limit as a hydrological meteorological constraint condition for judging quasi-synchronous observation: specifically, a relational graph of wind speed, atmospheric pressure and the like and interstellar deviation of radiation products is drawn by combining a multi-source satellite quasi-synchrotron radiation product on the basis of hydrometeological data, and wind speed and atmospheric pressure values corresponding to 2.5% of errors of the radiation products are determined by analyzing statistical characteristics of the interstellar deviation of the two satellite radiation products under different wind speed and atmospheric pressure conditions to serve as hydrometeological constraint conditions.
Hydrological meteorological constraint condition M0The value is specifically determined by the actual situation, VIIRS and FY-3D satellites are taken as columns,the obtained hydrometeorological constraint conditions selected by the effective pixels are as follows: the wind speed is not more than 7.8m/s and the atmospheric pressure is more than 100.5kPa but less than 102.5 kPa.
(3) Analyzing the quantitative relation between the space variation coefficient of the satellite image radiation product and the interstellar deviation of the radiation product, determining the maximum space variation coefficient of the satellite image on the premise of meeting the IOCCG quality target, and taking the maximum space variation coefficient as a constraint condition for judging the sea surface roughness of quasi-synchronous observation: in this embodiment, the spatial variation coefficient is specifically defined as a ratio of a variance to a mean of the radiation amount in a small region, and is used to represent a roughness of the spatial distribution of the radiation product, and a statistical relationship graph of the variation coefficient and the interplanetary deviation is drawn to determine the variation coefficient corresponding to the 2.5% radiation product error, which is used as a maximum value of the sea surface roughness, that is, a constraint upper limit.
And taking VIIRS and FY-3D satellites as columns, and obtaining the sea surface roughness constraint conditions as follows: the coefficient of spatial variation does not exceed 0.1.
Where CV is the spatial coefficient of variation of the radiation in the nxn region, ρ is the radiation signal, and subscript m represents the average.
(4) The method comprises the steps that the space-time variation of radiation signals of bright targets (such as clouds) is large, the radiation signals are not suitable for a data set for cross calibration or data quality evaluation, the radiation signals need to be removed, the minimum radiation values of all the bright targets in an image are determined by counting histograms of the radiation signals of the bright targets of a satellite image, and the minimum radiation values are used as constraint conditions and used for removing the influence of the bright targets on the data quality evaluation and cross calibration results.
And taking VIIRS and FY-3D satellites as columns, and obtaining the bright target constraint conditions as follows: the top reflectivity of the atmosphere layer is not more than 0.18.
(5) Analyzing the quantitative relation between the satellite transit time difference and the interstellar deviation of the image radiation product, determining the maximum transit time delay of the multisource satellite on the premise of meeting the IOCCG quality target, and taking the maximum transit time delay as a transit time delay constraint condition for judging the quasi-synchronous observation: the bigger the deviation of the transit time is, the radiation of two satellitesThe greater the interstellar bias of the product. In order to obtain a stable and reliable cross-calibration and data quality evaluation data set, it is necessary to establish a time constraint condition. By counting interstellar deviation characteristics of the radiation products at different transit times of the satellite, the transit time delay corresponding to 2.5% of errors of the radiation products is used as the maximum time delay of the two satellites in the embodiment. Transit time delay constraint T0The value is specifically based on the actual situation, VIIRS and FY-3D satellites are taken as columns, and the obtained time delay constraint condition is<For 20 minutes.
In addition, in step S2, the local variance of the satellite images not only includes the radiation data error information, but also includes signal anomalies caused by natural space changes of atmospheric and marine components, so that it is difficult for the conventional local variance method to obtain accurate radiation data errors. However, under different satellite quasi-synchronous observation conditions, signal anomalies caused by natural space changes of atmospheric and marine components have better spatial consistency (as shown in fig. 2), and can be used for improving the traditional method for extracting radiation data errors by local variance. In this embodiment, from the statistical characteristics of the data error of the radiation product, considering that the signal abnormality caused by the natural space change of the atmospheric and marine components has better spatial consistency, an error analysis equation is established, and the dynamic evaluation of the quality of the multi-source satellite radiation data is realized, where the step S2 specifically includes the following steps:
s21, selecting at least three transit satellites, and acquiring an effective image metadata set required by an analytic equation by using the constraint conditions in S1;
s22, segmenting the satellite image, establishing an N multiplied by N pixel local area sample set, and calculating the deviation between the pixel radiation quantity and the mean value in the local area as a local abnormal signal;
s23, converting the local abnormal signals into interstellar deviation, and eliminating the influence of the abnormal signals caused by natural processes on an error analytic equation;
and S24, dynamically monitoring the quality of the radiation product by calculating the variance of the interstellar deviation and establishing an error analysis equation of the multi-source satellite radiation product by using the assumption that the errors of the multi-source satellite radiation product are mutually independent.
The specific construction principle of the error analytic equation is further explained as follows:
(1) for industry-accepted sensors with higher radiometric calibration quality (e.g., VIIRS and MODIS satellites), the radiated signal observed by the satellite is expressed as:
ρs(λ)＝ρt(λ)+ε(λ)
in the formula, ρsRadiation signals, p, acquired for the satellitetIs the corresponding true value, ε is the observation error, and λ is the wavelength.
(2) And (3) not considering interstellar deviation caused by spectral response difference, the radiation signals obtained by the quasi-synchronous observation of the multisource satellite are as follows:
(3) the general image radiation signal abnormality mainly comprises signal abnormality caused by radiation product errors and natural space changes of natural compositions of natural substances of the atmosphere and the ocean, and can be expressed by a formula as follows: :
in the formula, delta rho is the total radiation abnormality of the satellite image pixel, and xi is a natural abnormal signal caused by observation geometry, atmospheric conditions, sea surface conditions and the like; the true error of the epsilon radiation product, the subscript represents three different satellites. In the range of 3 x 3 pixel areas, Δ ρ is expressed by the deviation of the pixel radiation signal from the mean value of the radiation signal in the area.
(4) Suppose that natural signal anomalies observed by different satellites have better consistency (i.e., ξ)1＝ξ2＝ξ3) Then, the total anomaly of the three satellite pixel radiation signals is subtracted by two to eliminate, and the interstellar deviation equation can be obtained as follows:
(5) the error of the satellite image radiation product can be represented by the variance of the image. Taking the image as a unit, and taking the variance (delta) of two sides of the formula in the step (4), so as to obtain the analytical equation of the error of the radiation product, wherein the analytical equation is as follows:
by numerically simulating the error generation process of the satellite radiation product by using FY-3D, VIIRS and MODIS satellites as columns, namely artificially adding known data errors into the radiation product data, and comparing the errors obtained by the error analysis equation with the known errors, the accuracy of the error analysis equation can reach 93 percent, and the method is specifically shown in FIG. 3.
In step S3, when the radiometric calibration error of the satellite image is large, the radiometric calibration error is regarded as a part of the radiometric product error and is transmitted to the calculation result of the error resolution equation, so that the error resolution equation established by the present invention is sensitive to the radiometric calibration quality and can be used to guide further radiometric calibration.
Taking FY-3D, VIIRS and MODIS satellite radiation products obtained by radiation transmission simulation as an example, through a numerical simulation method, a radiation scaling error of < + -20% is added to satellite radiation product data, and through comparing errors obtained by an analytic equation with known errors, the satellite radiation product error of < 32% can be increased by the radiation scaling error of < + -20%, as shown in FIG. 4. And (4) taking the effective image metadata set obtained under the constraint condition of effective image element selection in the step (S1) as the input of an error analysis equation, triggering to start cross calibration when the error of a certain satellite radiation product obtained by the error analysis equation exceeds a reasonable range (for example, the error requirement target of 2.5% of IOCCG), otherwise, judging that the radiation product has reasonable quality, and stopping execution.
In step S4, when the error of a certain satellite radiation product obtained by the error resolution equation is too large, the satellite needs to perform cross calibration on the satellite radiation data. And (4) performing cross calibration by taking the effective image element data set obtained under the constraint condition of effective image element selection in the step (S1) and the satellite (such as MODIS or VIIRS) with the industry-recognized high data quality as a reference satellite and the satellite to be subjected to cross calibration as a target satellite.
The specific implementation process comprises the following steps: establishing an empirical function relationship between the dependent variable and the independent variable through regression analysis by taking the radiation data of the reference satellite as the dependent variable and the radiation data of the target satellite as the independent variable; and applying the empirical function relationship to the radiation data of the target satellite to obtain the radiation product after cross calibration. Taking the FY-3D and VIIRS satellite radiation products as an example (VIIRS is the reference satellite and FY-3D is the satellite to be corrected), a stable and reliable cross calibration model can be established with a correlation coefficient of 0.90 or more under the condition of strictly complying with the constraint conditions of pixel selection, as shown in fig. 5. In addition, the cross-scaled FY-3D radiation product is used for producing the backscattering coefficient (b) of 551nm wave bandb(551) Has good consistency with similar products of a VIIRS satellite, and has interstellar deviation not exceeding 24 percent, as shown in figure 6.
In conclusion, the scheme effectively improves the quality of the water body satellite radiation product and the automation level of the water environment remote sensing data processing by constructing the multisource satellite quasi-synchronous effective pixel selection constraint condition, the water body radiation product quality evaluation and the dynamic cross radiometric calibration.
The above description is only a preferred embodiment of the present invention, and not intended to limit the present invention in other forms, and any person skilled in the art may apply the above modifications or changes to the equivalent embodiments with equivalent changes, without departing from the technical spirit of the present invention, and any simple modification, equivalent change and change made to the above embodiments according to the technical spirit of the present invention still belong to the protection scope of the technical spirit of the present invention.
Claims (8)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110257688.1A CN112884342B (en) | 2021-03-10 | 2021-03-10 | Quality evaluation and cross calibration method for water-color satellite atmospheric roof radiation product |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110257688.1A CN112884342B (en) | 2021-03-10 | 2021-03-10 | Quality evaluation and cross calibration method for water-color satellite atmospheric roof radiation product |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112884342A true CN112884342A (en) | 2021-06-01 |
| CN112884342B CN112884342B (en) | 2024-03-12 |
Family
ID=76054023
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110257688.1A Active CN112884342B (en) | 2021-03-10 | 2021-03-10 | Quality evaluation and cross calibration method for water-color satellite atmospheric roof radiation product |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112884342B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113407596A (en) * | 2021-06-23 | 2021-09-17 | 北京华云星地通科技有限公司 | Meteorological remote sensing data precision inspection analysis system based on cross radiometric calibration |
| CN113534203A (en) * | 2021-07-12 | 2021-10-22 | 长光卫星技术有限公司 | On-orbit cross-radiometric calibration method based on AERONET aerosol data |
| CN114383656A (en) * | 2021-12-22 | 2022-04-22 | 国家卫星气象中心(国家空间天气监测预警中心) | A Satellite Data Quality Monitoring Method Based on Bias of Intersection Double Reference Sources |
| CN114579655A (en) * | 2022-05-06 | 2022-06-03 | 国家卫星海洋应用中心 | Automatic planning method and device for ocean water color satellite radiometric calibration |
| CN114971148A (en) * | 2022-03-28 | 2022-08-30 | 陕西九州遥感信息技术有限公司 | Method for selecting pseudo-invariant cross transfer site for marine water color satellite |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107402395A (en) * | 2017-05-17 | 2017-11-28 | 西安交通大学 | A kind of satellite selection method to be navigated for single system and multisystem combinations of satellites |
| CN108919151A (en) * | 2018-04-03 | 2018-11-30 | 西安空间无线电技术研究所 | A kind of microwave radiometer cross polarization error correcting method |
| CN109213964A (en) * | 2018-07-13 | 2019-01-15 | 中南大学 | A kind of satellite AOD product bearing calibration for merging multi-source feature geographic factor |
| CN109655846A (en) * | 2019-01-30 | 2019-04-19 | 清华大学 | A kind of multistation difference post-processing high-precision time synchronization method and system |
| CN109903352A (en) * | 2018-12-24 | 2019-06-18 | 中国科学院遥感与数字地球研究所 | A method for making large-area seamless orthophotos of satellite remote sensing images |
| CN110286354A (en) * | 2019-04-23 | 2019-09-27 | 中国人民解放军63921部队 | The method, apparatus and computer readable storage medium of multi-target detection and differentiation |
| CN111523209A (en) * | 2020-04-09 | 2020-08-11 | 北京空间机电研究所 | A method for land resource satellite calibration orbit planning and benchmark load orbit optimization |
| CN111737850A (en) * | 2020-05-15 | 2020-10-02 | 中国科学院空天信息创新研究院 | A Multi-source Satellite AOD Fusion Method Based on Uncertainty at Pixel Scale |
| CN112200050A (en) * | 2020-09-30 | 2021-01-08 | 同济大学 | Lake water volume space-time change analysis method and device based on multi-source satellite data |
| CN112197865A (en) * | 2020-09-02 | 2021-01-08 | 华中科技大学 | Estimation method and system for observation brightness temperature data error of satellite-borne microwave radiometer |
| AU2020103576A4 (en) * | 2019-12-27 | 2021-02-04 | Wuhan University | Autonomous orbit and attitude determination method of low-orbit satellite based on non-navigation satellite signal |
-
2021
- 2021-03-10 CN CN202110257688.1A patent/CN112884342B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107402395A (en) * | 2017-05-17 | 2017-11-28 | 西安交通大学 | A kind of satellite selection method to be navigated for single system and multisystem combinations of satellites |
| CN108919151A (en) * | 2018-04-03 | 2018-11-30 | 西安空间无线电技术研究所 | A kind of microwave radiometer cross polarization error correcting method |
| CN109213964A (en) * | 2018-07-13 | 2019-01-15 | 中南大学 | A kind of satellite AOD product bearing calibration for merging multi-source feature geographic factor |
| CN109903352A (en) * | 2018-12-24 | 2019-06-18 | 中国科学院遥感与数字地球研究所 | A method for making large-area seamless orthophotos of satellite remote sensing images |
| CN109655846A (en) * | 2019-01-30 | 2019-04-19 | 清华大学 | A kind of multistation difference post-processing high-precision time synchronization method and system |
| CN110286354A (en) * | 2019-04-23 | 2019-09-27 | 中国人民解放军63921部队 | The method, apparatus and computer readable storage medium of multi-target detection and differentiation |
| AU2020103576A4 (en) * | 2019-12-27 | 2021-02-04 | Wuhan University | Autonomous orbit and attitude determination method of low-orbit satellite based on non-navigation satellite signal |
| CN111523209A (en) * | 2020-04-09 | 2020-08-11 | 北京空间机电研究所 | A method for land resource satellite calibration orbit planning and benchmark load orbit optimization |
| CN111737850A (en) * | 2020-05-15 | 2020-10-02 | 中国科学院空天信息创新研究院 | A Multi-source Satellite AOD Fusion Method Based on Uncertainty at Pixel Scale |
| CN112197865A (en) * | 2020-09-02 | 2021-01-08 | 华中科技大学 | Estimation method and system for observation brightness temperature data error of satellite-borne microwave radiometer |
| CN112200050A (en) * | 2020-09-30 | 2021-01-08 | 同济大学 | Lake water volume space-time change analysis method and device based on multi-source satellite data |
Non-Patent Citations (22)
| Title |
|---|
| PENGHAI WU等: "Integrated fusion of multi-scale polar-orbiting and geostationary satellite observations for the mapping of high spatial and temporal resolution land surface temperature", REMOTE SENSING OF ENVIRONMENT, pages 169 - 181 * |
| RUGANG TANG等: "Multi-source high-resolution satellite products in Yangtze Estuary: cross-comparisons and impacts of signal-to-noise ratio and spatial resolution", OPT. EXPRESS, vol. 27, pages 6426 - 6441 * |
| 傅俏燕;闵祥军;李杏朝;沙崇漠;李小英;马国强;潘志强;郭毅;李启明;刘国栋;: "敦煌场地CBERS-02 CCD传感器在轨绝对辐射定标研究", 遥感学报, no. 04, pages 3 - 9 * |
| 吴兴等: "SPARK卫星高光谱数据辐射质量评价", 遥感技术与应用, vol. 33, no. 2, pages 233 - 240 * |
| 崔廷伟等: "北冰洋卫星水色遥感观测能力评价与展望", 中国海洋大学学报, vol. 51, no. 1, pages 125 - 137 * |
| 张志杰;张浩;常玉光;陈正超;: "Landsat系列卫星光学遥感器辐射定标方法综述", 遥感学报, no. 05, pages 21 - 34 * |
| 徐伟进: "国产光学卫星遥感数据 VNIR 静止卫星反射波段辐射性能评估和辐射归一化", 中国优秀硕士学位论文全文数据库基础科学辑（月刊）, no. 01, pages 009 - 128 * |
| 徐文等: ""高分二号"卫星相机影像辐射质量评价", 航天返回与遥感, vol. 36, no. 4, pages 1 - 9 * |
| 李军;朱建华;韩冰;高飞;杨安安;赵屹立;: "VIIRS在中国渤海的遥感反射率产品验证", 海洋技术学报, no. 02, pages 30 - 36 * |
| 杨磊;周兴华;徐全军;柯宝贵;穆博;朱琳;: "卫星高度计定标现状", 遥感学报, no. 03, 25 May 2019 (2019-05-25), pages 38 - 53 * |
| 漆成莉;胡秀清;张里阳;董超华;: "基于高光谱资料对FY-1C/1D气象卫星进行交叉定标", 气象学报, no. 04, 15 August 2012 (2012-08-15), pages 286 - 295 * |
| 牛明慧: "光学遥感仪器月球定标技术研究", 中国博士学位论文全文数据库基础科学辑（月刊）, no. 05, 15 May 2019 (2019-05-15), pages 007 - 3 * |
| 王春琴;张贤国;李嘉巍;黄聪;张效信;荆涛;沈国红;张申毅;曹光伟;孙越强;梁金宝;朱光武;韩英;: "卫星辐射带高能带电粒子观测数据交叉定标-FY3和NOAA应用示例", 中国科学:技术科学, no. 11, 20 November 2013 (2013-11-20), pages 90 - 96 * |
| 程扬等: "GPM_IMERGE卫星遥感降水数据在巢湖流域的精度评价", 水土保持研究, vol. 27, no. 5, pages 188 - 193 * |
| 胡建国, 章传银, 常晓涛: "近海多卫星测高数据联合处理的方法及应用", 测绘通报, no. 01, 25 January 2004 (2004-01-25), pages 1 - 4 * |
| 荣飘;张淳民;刘冬冬;栗彦芬;: "三峡库区CO_2平均柱浓度的卫星和地面观测研究", 三峡生态环境监测, no. 03, pages 12 - 21 * |
| 蒋虎, 黄, 严豪键: "无线电掩星反演大气参数误差分析及应用进展", 地球物理学进展, no. 01, pages 85 - 91 * |
| 郭强;陈博洋;张勇;徐娜;冯绚;杨昌军;王新;韩秀珍;李元;胡秀清;张志清;: "风云二号卫星在轨辐射定标技术进展", 气象科技进展, no. 06, pages 8 - 14 * |
| 韩启金;张学文;乔志远;杨磊;潘志强;刘李;: "高分一号卫星PMS相机多场地宽动态辐射定标", 红外与激光工程, no. 01, pages 133 - 139 * |
| 马艳红: "卫星辐射数据质量分析与研究", 中国优秀硕士学位论文全文数据库工程科技Ⅱ辑（月刊）, no. 02, pages 028 - 260 * |
| 高彩霞;姜小光;马灵玲;霍红元;: "传感器交叉辐射定标综述", 干旱区地理, no. 01, pages 141 - 148 * |
| 黄端;邱玉宝;石利娟;武胜利;刘小生;阮永俭;: "两种星载微波辐射计被动亮温数据的交叉定标", 测绘科学, no. 01, pages 139 - 145 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113407596A (en) * | 2021-06-23 | 2021-09-17 | 北京华云星地通科技有限公司 | Meteorological remote sensing data precision inspection analysis system based on cross radiometric calibration |
| CN113407596B (en) * | 2021-06-23 | 2022-07-15 | 北京华云星地通科技有限公司 | Meteorological remote sensing data precision inspection analysis system based on cross radiometric calibration |
| CN113534203A (en) * | 2021-07-12 | 2021-10-22 | 长光卫星技术有限公司 | On-orbit cross-radiometric calibration method based on AERONET aerosol data |
| CN113534203B (en) * | 2021-07-12 | 2022-08-23 | 长光卫星技术股份有限公司 | On-orbit cross-radiation calibration method based on AERONET aerosol data |
| CN114383656A (en) * | 2021-12-22 | 2022-04-22 | 国家卫星气象中心(国家空间天气监测预警中心) | A Satellite Data Quality Monitoring Method Based on Bias of Intersection Double Reference Sources |
| CN114383656B (en) * | 2021-12-22 | 2023-09-29 | 国家卫星气象中心(国家空间天气监测预警中心) | Satellite data quality monitoring method based on cross point double-reference source deviation |
| CN114971148A (en) * | 2022-03-28 | 2022-08-30 | 陕西九州遥感信息技术有限公司 | Method for selecting pseudo-invariant cross transfer site for marine water color satellite |
| CN114579655A (en) * | 2022-05-06 | 2022-06-03 | 国家卫星海洋应用中心 | Automatic planning method and device for ocean water color satellite radiometric calibration |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112884342B (en) | 2024-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112884342A (en) | Water color satellite atmospheric layer top radiation product quality evaluation and cross calibration method | |
| WO2020207070A1 (en) | Method and system for evaluating shenzhen sea water quality | |
| CN114372707A (en) | High-cold-wetland degradation degree monitoring method based on remote sensing data | |
| CN101980293A (en) | A MTF Detection Method of Hyperspectral Remote Sensing System Based on Edge Image | |
| CN107368617B (en) | Ground-air detection infrared imaging system action distance calculation method based on Lowtran7 atmospheric software | |
| CN113935249B (en) | Upper-layer ocean thermal structure inversion method based on compression and excitation network | |
| CN105303030A (en) | MODIS remote sensing evaluation method for eutrophication lake algae gauss vertical distribution structural parameters | |
| CN114047563B (en) | All-weather assimilation method for infrared hyperspectrum | |
| Yang et al. | Evaluation and mitigation of rain effect on wave direction and period estimation from X-band marine radar images | |
| Li et al. | A method of visibility detection based on the transfer learning | |
| Chen et al. | A simple atmospheric correction algorithm for MODIS in shallow turbid waters: A case study in Taihu Lake | |
| Massari et al. | Error and uncertainty characterization | |
| Peres et al. | Analyses of the positive bias of remotely sensed SST retrievals in the coastal waters of Rio de Janeiro | |
| KR20210018737A (en) | Apparatus and method for calculating optical properties of aerosol | |
| Zhang et al. | Watershed runoff modeling through a multi-time scale approach by multivariate empirical mode decomposition (MEMD) | |
| Treadon et al. | Global and regional moisture analyses at NCEP | |
| Wilson et al. | Surfzone state estimation, with applications to quadcopter-based remote sensing data | |
| Zarei et al. | Evaluating the changes in gavkhuni wetland using modis satellite images in 2000-2016. | |
| Zhang et al. | Assimilation of remote sensing observations into a sediment transport model of China’s largest freshwater lake: spatial and temporal effects | |
| Liu et al. | Identify optimal predictors of statistical downscaling of summer daily precipitation in China from three-dimensional large-scale variables | |
| Liu et al. | Gridded statistical downscaling based on interpolation of parameters and predictor locations for summer daily precipitation in North China | |
| CN109029735B (en) | Method and device for calculating surface emissivity | |
| Chen et al. | Scale correction of two-band ratio of red to near-infrared using imagery histogram approach: A case study on Indian remote sensing satellite in Yellow River estuary | |
| Liu et al. | Image-based recognition and processing system for monitoring water levels in an irrigation and drainage channel | |
| Kornienko | Analysis of errors in estimating changes in water body areas by satellite data: Case study of thermokarst lakes in Yamal Peninsula |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |