CN112884342B - Quality evaluation and cross calibration method for water-color satellite atmospheric roof radiation product - Google Patents
Quality evaluation and cross calibration method for water-color satellite atmospheric roof radiation product Download PDFInfo
- Publication number
- CN112884342B CN112884342B CN202110257688.1A CN202110257688A CN112884342B CN 112884342 B CN112884342 B CN 112884342B CN 202110257688 A CN202110257688 A CN 202110257688A CN 112884342 B CN112884342 B CN 112884342B
- Authority
- CN
- China
- Prior art keywords
- radiation
- satellite
- quality
- radiation product
- product
- 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.)
- Active
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 180
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000013441 quality evaluation Methods 0.000 title claims abstract description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 32
- 238000004458 analytical method Methods 0.000 claims abstract description 29
- 230000002159 abnormal effect Effects 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims abstract description 10
- 238000012544 monitoring process Methods 0.000 claims abstract description 9
- 238000011156 evaluation Methods 0.000 claims abstract description 8
- 230000003746 surface roughness Effects 0.000 claims description 13
- 238000010586 diagram Methods 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 8
- 230000001419 dependent effect Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000007405 data analysis Methods 0.000 claims description 3
- 238000000611 regression analysis Methods 0.000 claims description 3
- 239000005445 natural material Substances 0.000 claims description 2
- 238000013139 quantization Methods 0.000 claims 1
- 230000005856 abnormality Effects 0.000 abstract description 5
- 230000007547 defect Effects 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 230000007246 mechanism Effects 0.000 abstract description 2
- 230000008901 benefit Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000001303 quality assessment method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011835 investigation Methods 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
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
Abstract
The invention discloses a quality evaluation and cross calibration method for a water-based satellite atmospheric top radiation product, which comprises the steps of establishing a multisource satellite quasi-synchronization effective pixel selection constraint condition; the radiation signal abnormality caused by the natural process is found to have stronger interstellar correlation; designing an error analysis equation of a multi-source satellite radiation product; and combines a dynamic cross radiometric calibration mechanism; the effect of the radiation product quality on the water color product is determined. The scheme standardizes the constraint condition of quasi-synchronous observation for the first time, 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 radiation product, overcomes the defect that the traditional method is seriously influenced by abnormal signals caused by the natural process, and improves the quality and the automatic service level of the water radiation product; 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-based satellite atmospheric top radiation product.
Background
The river-along, coastal and other areas with the most potential and activity for the prosperous development of human society are important ecological environment bearing areas, and have great economic benefits. As the population continues to gather in recent years in lakes, rivers, coasts, etc., water ecology is increasingly pressurized and resource and environmental problems are more serious. The enhancement of monitoring and management of water environment in areas such as coastal, river-like and coastal areas has become an important problem to be solved by various levels of governments. The remote sensing data has the advantages of wide coverage area, diversity of spatial scales, continuous time scales, abundant spectrum information, flexible and convenient observation and the like, and can provide an effective means for dynamically monitoring the resource environment conditions of lakes, rivers and seas. Since the Landsat satellite is successfully transmitted in 1972, people accumulate remote sensing images acquired in the life cycle of thousands of satellites, and have great success in disaster prevention, disaster reduction, forecast and early warning, resource investigation and other aspects, thereby protecting the living environment and development space of the people and helping to build in China.
The water body water-leaving radiation belongs to typical weak signals, so that the water environment remote sensing has relatively high requirements on satellite image quality. Dynamic changes in gravitational environment, etc., are necessary to enhance the dynamic monitoring of satellite radiation product quality in order to guide the development of radiometric calibration work. However, most satellites do not possess equipment for quality detection of the radiation product, and quality assessment of the radiation product is typically achieved using image statistics. However, besides the data error can cause radiation signal abnormality, natural factors such as atmospheric conditions, observation geometry, seawater characteristics and the like can also cause radiation signal abnormality, so that the reliability of the traditional method is greatly reduced. In addition, due to the lack of strict quasi-synchronous effective pixel selection constraint conditions, high-precision radiation products are difficult to obtain by traditional cross calibration, and the business application of the water radiation products is seriously affected.
Disclosure of Invention
The invention provides a method for evaluating and cross-scaling the quality of a water-based satellite atmospheric top radiation product, which aims to solve the defects in the aspects of dynamic monitoring and cross-scaling of the quality of the traditional water body radiation product, dynamically realizes the cross-scaling of the satellite radiation product, so as to improve the quality of satellite image data and better serve for monitoring 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-based satellite atmospheric 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 data 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 synchronicity of multi-source satellite observation data; under the constraint condition, extracting an effective pixel data set required by radiation product quality evaluation and cross calibration from the multi-source satellite radiation product data;
s2, constructing an error analysis equation: establishing a multisource satellite radiation product error analysis equation by analyzing consistency characteristics of signal anomalies caused by natural changes of atmosphere and marine substances detected by different satellites in quasi-synchronization;
step S3, quality evaluation: extracting error information of the multi-source satellite radiation product according to the effective pixel data set and the multi-source satellite radiation product error analysis equation, and carrying out 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 pixel data set under the constraint condition as input to carry out dynamic cross calibration correction on the satellite radiation product.
Further, in the step S1, when determining the constraint condition of the illumination-observation geometry, the illumination-observation geometry is represented by the zenith angle of the sun, the zenith angle of the satellite and the relative azimuth angle, and the maximum value and the minimum value of the difference in the aspect of the pixel illumination-observation geometry allowed by the quasi-synchronous observation of the multi-source satellite are obtained as the constraint condition of the illumination-observation geometry for judging the quasi-synchronous observation on the premise of meeting the quality target of the IOCCG by using radiation transmission simulation.
Further, in the step S1, when determining the hydrological constraint condition, the relationship diagram of the interplanetary deviation of the severe hydrological weather and the radiation product is drawn by combining the multisource satellite quasi-synchrotron radiation product based on the hydrological weather data, the severe hydrological weather includes wind speed and atmospheric pressure, and the statistical characteristics of the interplanetary deviation of the two satellite radiation products under different severe hydrological weather conditions are analyzed, and the upper limit of the severe hydrological weather is determined as the hydrological constraint condition for judging quasi-synchronous observation on the premise of meeting the IOCCG quality target.
Further, in the step S1, when determining the sea surface roughness constraint condition, the quantized relationship between the space variation coefficient of the satellite image radiation product and the interplanetary deviation of the radiation product is analyzed, and the maximum space 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.
Further, in the step S1, when determining the transit time delay, the maximum transit time delay of the multi-source satellite is determined as a constraint condition for determining the transit time delay of the quasi-synchronous observation by analyzing the quantitative relationship between the satellite transit time difference and the interstellar deviation of the image radiation product on the premise of meeting the IOCCG quality target.
Further, in the step S1, when determining the constraint condition of the bright target, the minimum radiation value of all the bright targets in the image is determined as the constraint condition by counting the histogram of the bright target radiation signal of the satellite image.
Further, the step S2 is specifically performed by the following manner when the error resolution equation is established:
s21, selecting at least three transit satellites, and acquiring an effective pixel data set required by an error analysis equation based on the constraint condition in the step S1;
s22, dividing 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 to be used as a local abnormal signal;
s23, converting the local abnormal signals into interstellar deviations, and eliminating the influence of the abnormal signals caused by the natural process on an error analysis equation;
s24, by using the assumption that the errors of the radiation products of the multi-source satellites are mutually independent, a multi-source satellite radiation product error analysis equation is established by calculating the variance of the interplanetary deviation, and the radiation product quality is dynamically monitored.
Further, in the step S4, when the error of the radiation product of a certain satellite obtained by the error analysis equation is too large, cross calibration is performed, specifically, the radiation data of the reference satellite is taken as a dependent variable, the radiation data of the target satellite is taken as an independent variable, an empirical function relation between the dependent variable and the independent variable is established based on regression analysis, and the empirical function relation is applied to the radiation data of the target satellite, so that the cross calibrated radiation product is obtained.
Compared with the prior art, the invention has the advantages and positive effects that:
according to the method, the influence of illumination-observation geometry, hydrological weather, transit time delay, sea surface roughness and the like on the stability of radiation signals is analyzed, and a multisource satellite quasi-synchronization effective pixel selection constraint condition is established; determining that the radiation signal abnormality caused by the natural process has stronger interstellar correlation by comparing the correlation characteristics of the quasi-synchronous satellite signal abnormality; by analyzing the statistical characteristics of the data errors of the radiation products, a multi-source satellite radiation product error analysis equation is designed, and the influence of the quality of the radiation products on the water color products is determined by combining a dynamic cross radiation calibration mechanism;
the scheme standardizes the constraint condition of quasi-synchronous observation for the first time, 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 radiation product, overcomes the defect that the traditional method is seriously influenced by abnormal signals caused by the natural process, and improves the quality and the automatic service level of the water radiation product; 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 value.
Drawings
FIG. 1 is a flow chart of radiation product quality evaluation and cross-scaling in accordance with an embodiment of the present invention;
FIG. 2 is a schematic diagram showing the inter-satellite correlation of anomaly signals resulting from the natural process according to an embodiment of the present invention, wherein (a) and (b) correspond to the spatial anomaly signals representing VIIRS and MODIS satellite images, respectively, and (c) is a histogram of the images of (a) and (b);
FIG. 3 is a schematic diagram of the radiation transmission simulation verification of the reliability of the error resolution equation according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of the effect of radiometric calibration errors on the error resolution equations 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 an example of the interplanetary consistency evaluation of a cross-scaled corrected water color product.
Detailed Description
In order that the above objects, features and advantages of the invention will be more readily understood, a further description of the invention will be rendered by reference to specific embodiments thereof which are illustrated in the appended drawings. 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 otherwise than as described herein, and therefore the present invention is not limited to the specific embodiments disclosed below.
In an embodiment, the present embodiment proposes a quality evaluation and cross calibration method for a top radiation product of an atmospheric layer of a water-based satellite, as shown in fig. 1, including the following steps:
step S1, establishing a multisource satellite quasi-synchronization effective pixel selection constraint condition, and extracting an effective pixel data set:
taking the IOCCG product quality requirement as a target guide, and quantitatively establishing constraint conditions for judging the illumination-observation geometry, the hydrological weather, the transit time delay, the sea surface roughness and the bright target of quasi-synchronous observation through radiation transmission simulation and quasi-synchronous satellite data analysis, so as to judge the synchronism of multi-source satellite observation data; under the constraint condition, extracting an effective pixel data set required by radiation product quality evaluation and cross calibration from the multi-source satellite radiation product data;
s2, signal anomalies caused by natural changes of the atmosphere and marine substances detected by quasi-synchronization of different satellites are analyzed to obtain consistency characteristics of the signal anomalies, a method of subtracting total anomaly signals of different satellites is adopted to eliminate the signal anomalies caused by the natural changes, and a multi-source satellite radiation product error analysis equation is established;
s3, extracting error information of the multi-source satellite radiation product according to the effective pixel data set and an error analysis equation of the multi-source satellite radiation product, and carrying out 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 S4, taking the effective pixel data 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 advantageous the cross-scale modeling and the correct evaluation of data quality. However, due to the influence of factors such as illumination-observation geometry, hydrological weather, transit time delay, sea surface roughness and the like, great deviation exists among the radiation data of the multi-source satellite, and based on the consideration, the embodiment provides a set of scientific and strict constraint conditions to remove the radiation data with larger interstellar deviation. In this embodiment, the constraint condition of effective pixel selection in step S1 is that when the data quality requirement of the IOCCG water color satellite is satisfied (for example, the radiation signal error of 443nm band is not more than 2.5%), the maximum or minimum value of the illumination-observation geometry (including solar zenith angle, satellite zenith angle and relative azimuth angle), the hydrological weather, the transit time delay, the sea surface roughness and the bright pixel radiation amount is specifically implemented by the following manner:
(1) On the premise of meeting the IOCCG quality target, obtaining the maximum value and the minimum value of the difference in the pixel illumination-observation geometry aspect allowed by the quasi-synchronous observation of the multi-source satellite as the illumination-observation geometry constraint condition for judging the quasi-synchronous observation by utilizing radiation transmission simulation:
the method comprises the steps of representing illumination-observation geometry by using a solar zenith angle, a satellite zenith angle and a relative azimuth angle, fixing the relative azimuth angle and the satellite zenith angle, keeping the solar zenith angle dynamically changed within a range of 0-60 degrees, drawing a contour map of radiation errors caused by different solar zenith angles and solar zenith angle interstar deviations, determining a distribution area of solar zenith angles and solar zenith angle interstar deviations corresponding to the radiation product errors of <2.5%, and taking the upper, lower, left and right boundaries of a rectangle with the largest built-in area of the contour area of <2.5% as the respective minimum value and maximum value of the solar zenith angles and the solar zenith angle interstar deviations. The maximum allowable dynamic range (maximum and minimum) of the satellite zenith angle and the satellite zenith angle deviation and the relative azimuth angle deviation can be obtained by the same method.
Illumination-observation geometry constraint A 0 Taking the actual situation as the reference, taking VIIRS and FY-3D satellites as columns, and obtaining the constraint conditions of illumination-observation geometry of effective pixel selection: 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 condition, discussing the statistical characteristics of the satellite image interstellar deviation, establishing a quantized linkage relation of the satellite image interstellar deviation and the interstellar deviation, and determining the upper limit of the severe hydrological weather on the premise of meeting the IOCCG quality target as the hydrological weather constraint condition for judging quasi-synchronous observation: specifically, based on hydrological meteorological data, a relation diagram of the interplanetary deviation of a radiation product such as wind speed and atmospheric pressure is drawn by combining a multisource satellite quasi-synchrotron radiation product, and the wind speed and atmospheric pressure corresponding to 2.5% of the radiation product error are determined as the hydrological constraint condition by analyzing the statistical characteristics of the interplanetary deviation of two satellite radiation products under different wind speed and atmospheric pressure conditions.
Hydrological constraint M 0 Taking actual conditions as the reference, taking VIIRS and FY-3D satellites as columns, and obtaining the constraint conditions of the effective pixel selection hydrological weather: the wind speed is not more than 7.8m/s, and the atmospheric pressure is more than 100.5kPa but less than 102.5kPa.
(3) The quantitative relation between the space variation coefficient of the satellite image radiation product and the interplanetary deviation of the radiation product is analyzed, and the maximum space variation coefficient of the satellite image is determined on the premise of meeting the IOCCG quality target and is used 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 variance to mean of the radiant quantity in a small area, which is used to characterize the roughness of the spatial distribution of the radiant product, and the variation coefficient corresponding to 2.5% of the radiant product error is determined as the maximum value of the sea surface roughness, i.e. the constraint upper limit, by drawing a statistical relationship diagram of the variation coefficient and the interplanetary deviation.
Taking VIIRS and FY-3D satellites as columns, the obtained sea surface roughness constraint conditions are as follows: the spatial coefficient of variation does not exceed 0.1.
Where CV is the spatial coefficient of variation of radiation in the N region, ρ is the radiation signal, and subscript m represents the average value.
(4) The method is characterized in that the space-time variation of the radiation signals of the bright targets (such as clouds) is large, the data set which is not suitable for cross calibration or data quality evaluation is required to be removed, the minimum radiation value of all the bright targets in the images is determined by counting the histogram of the radiation signals of the bright targets of the satellite images and is used as a constraint condition for removing the influence of the bright targets on the data quality evaluation and the cross calibration result.
Taking VIIRS and FY-3D satellites as columns, the obtained constraint conditions of the bright target are as follows: the atmospheric top reflectivity is no more than 0.18.
(5) The quantitative relation between the time difference of satellite transit and the interstellar deviation of an image radiation product is analyzed, and the maximum transit time delay of the multi-source satellite is determined on the premise of meeting the IOCCG quality target and is used as a transit time delay constraint condition for judging quasi-synchronous observation: the larger the transit time bias, the larger the interplanetary bias for the two satellite radiation products. To obtain a stable, reliable cross-scaling and data quality assessment dataset, it is necessary to establish time constraints. By counting the inter-satellite deviation characteristics of the radiation products under different transit times of the satellites, the embodiment takes the transit time delay corresponding to the 2.5% radiation product error as the maximum time delay of the two satellites. Transit time delay constraint T 0 Taking the actual situation as the reference, taking VIIRS and FY-3D satellites as columns, and obtaining the time delay constraint condition as<20 minutes.
In addition, in step S2, the local variance of the satellite image includes not only the radiation data error information, but also signal anomalies caused by natural spatial variations of the atmospheric and marine components, so that it is difficult to obtain accurate radiation data errors by the conventional local variance method. However, under different satellite quasi-synchronous observation conditions, the signal anomalies caused by natural spatial variation of the atmospheric and ocean 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. The embodiment starts from the statistical characteristics of the radiation product data errors, considers that signal anomalies caused by natural spatial variation of atmospheric and ocean components have better spatial consistency, establishes an error analysis equation, and realizes the dynamic evaluation of the multi-source satellite radiation data quality, and the step S2 specifically comprises the following steps:
s21, selecting at least three transit satellites, and acquiring an effective pixel data set required by an analytical equation by using the constraint conditions in S1;
s22, dividing 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 to be used as a local abnormal signal;
s23, converting the local abnormal signals into interstellar deviations, and eliminating the influence of the abnormal signals caused by the natural process on an error analysis equation;
s24, by using the assumption that the errors of the radiation products of the multi-source satellites are mutually independent, the dynamic monitoring of the quality of the radiation products is carried out by calculating the variance of the interplanetary deviation and establishing an error analysis equation of the radiation products of the multi-source satellites.
The specific construction principle of the error resolution equation is further explained as follows:
(1) For industry-accepted radiometric calibration quality sensors (e.g., VIIRS and MODIS satellites), the radiation signal observed by the satellite is expressed as:
ρ s (λ)＝ρ t (λ)+ε(λ)
wherein ρ is s For the radiation signals obtained by satellites ρ t The true value is the corresponding value, epsilon is the observed error, and lambda is the wavelength.
(2) Irrespective of the interplanetary bias caused by the difference in spectral response, the radiation signal obtained from the quasi-simultaneous observation of the multi-source satellite is as follows:
(3) The total anomalies of the image radiation signals mainly consist of the anomalies of the signals caused by errors of radiation products and natural spatial variations of the components of natural substances of the atmosphere and the ocean, and can be expressed as follows: :
wherein Deltaρ is total anomaly of satellite image pixel radiation, and ζ is natural anomaly signal caused by observation geometry, atmospheric condition, sea surface condition and the like; the true error of the epsilon radiation product, subscript, represents three different satellites. Wherein, within the range of the 3×3 pixel area, Δρ is represented by the deviation of the pixel radiation signal from the average of the radiation signals in the area.
(4) It is assumed 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 eliminated by subtracting two phases, and an interplanetary deviation equation can be obtained as follows:
(5) Errors in the satellite image radiation product can be represented by variances in the images. Taking the image as a unit, taking the variance (delta) from both sides of the formula in the step (4), and obtaining an analytical equation of the radiation product error as follows:
by taking FY-3D, VIIRS and MODIS satellites as columns, the accuracy of an error analysis equation can reach 93% by numerically simulating the satellite radiation product error generation process, namely manually adding known data errors into radiation product data and comparing the errors obtained by the error analysis equation with the known errors, and the accuracy of the error analysis equation is particularly shown in figure 3.
In step S3, when the radiometric calibration error of the satellite image is larger, the radiometric calibration error is regarded as a part of the radiometric product error and is transmitted to the calculation result of the error analysis equation, so that the error analysis equation established by the invention is more sensitive to radiometric calibration quality and can be used for guiding further radiometric calibration.
Taking FY-3D, VIIRS obtained by radiation transmission simulation and MODIS satellite radiation products as examples, adding the radiation calibration error of < + -20% to satellite radiation product data by a numerical simulation method, and comparing the error obtained by comparing an analytical equation with the known error, it is known that the radiation calibration error of < + -20% can be increased by <32% of satellite radiation product error, as shown in figure 4. And (3) taking the effective pixel data set obtained under the effective pixel selection constraint condition 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 the IOCCG), otherwise, judging that the quality of the radiation product is reasonable, and stopping execution.
In step S4, when the error of the radiation product of a certain satellite obtained by the error analysis equation is too large, the satellite needs to perform cross calibration on the radiation data of the satellite. And (3) taking the effective pixel data set obtained under the effective pixel selection constraint condition in the step S1, taking a satellite (such as MODIS or VIIRS) with higher industry-accepted data quality as a reference satellite, taking a satellite to be cross-scaled as a target satellite, and implementing cross scaling.
The specific implementation process is as follows: establishing an empirical function relation between the dependent variable and the independent variable through regression analysis by taking radiation data of a reference satellite as the dependent variable and radiation data of a target satellite as the independent variable; by applying the empirical functional relationship to the radiation data of the target satellite, a cross-scaled radiation product is obtained. Taking FY-3D and VIIRS satellite radiation products as examples (wherein VIIRS is a reference satellite and FY-3D is a satellite to be corrected), a stable and reliable cross calibration model can be established on the premise of strictly adhering to pixel selection constraint conditions, and the correlation coefficient can reach more than 0.90, as shown in figure 5. Furthermore, cross-scaled FY-3D radiation products for producing backscatter coefficients (b) in the 551nm band b (551) Good with similar products of VIIRS satellitesConsistency, with an interplanetary deviation of no more than 24%, as shown in fig. 6.
In conclusion, the scheme effectively improves the quality of the water body satellite radiation product and the water environment remote sensing data processing automation level by constructing multisource satellite quasi-synchronous effective pixel selection constraint conditions, water body radiation product quality evaluation and dynamic cross radiation calibration.
The present invention is not limited to the above-mentioned embodiments, and any equivalent embodiments which can be changed or modified by the technical content disclosed above can be applied to other fields, but any simple modification, equivalent changes and modification made to the above-mentioned embodiments according to the technical substance of the present invention without departing from the technical content of the present invention still belong to the protection scope of the technical solution of the present invention.
Claims (7)
1. The quality evaluation and cross calibration method for the top radiation product of the water-based satellite atmosphere is characterized by comprising the following steps of:
step S1, establishing a multisource satellite quasi-synchronization effective pixel selection constraint condition, and extracting an effective pixel data 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 synchronicity of multi-source satellite observation data; under the constraint condition, extracting an effective pixel data set required by radiation product quality evaluation and cross calibration from the multi-source satellite radiation product data;
s2, constructing an error analysis equation: establishing a multisource satellite radiation product error analysis equation by analyzing consistency characteristics of signal anomalies caused by natural changes of atmosphere and marine substances detected by different satellites in quasi-synchronization;
s21, selecting at least three transit satellites, and acquiring an effective pixel data set required by an error analysis equation based on the constraint condition in the step S1;
s22, dividing 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 to be used as a local abnormal signal;
the total anomalies of the image radiation signals consist of the anomalies of the signals caused by errors of radiation products and natural spatial variations of the components of natural substances of the atmosphere and the ocean, expressed as follows:
wherein Deltaρ is total anomaly of satellite image pixel radiation, subscript represents three different satellites, ζ is natural anomaly signal caused by observation geometry, atmospheric condition and sea surface condition, ε is true error of radiation product, subscript represents three different satellites, λ is wavelength;
s23, converting the local abnormal signals into interstellar deviations, and eliminating the influence of the abnormal signals caused by the natural process on an error analysis equation;
assuming that natural signal anomalies observed by different satellites have consistency, the total anomaly of three satellite pixel radiation signals is eliminated by subtracting two by two, and an interstellar deviation equation is obtained as follows:
s24, by using the assumption that the errors of the radiation products of the multi-source satellites are mutually independent, establishing an error analysis equation of the radiation products of the multi-source satellites by calculating the variance of the interplanetary deviation, and carrying out dynamic monitoring on the quality of the radiation products;
error of radiation product the analytical equation is as follows:
wherein,and->The corresponding unbiased estimated variances of the 1 st, 2 nd and 3 rd satellites are represented respectively, +.>Represents the unbiased estimated variance between two satellites 1 and 2,/and>represents the unbiased estimated variance between the 2 nd and 3 rd satellites, ++>Representing the unbiased estimated variance between the 1 st and 3 rd satellites;
step S3, quality evaluation: extracting error information of the multi-source satellite radiation product according to the effective pixel data set and the multi-source satellite radiation product error analysis equation, and carrying out 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 pixel data set under the constraint condition as input to carry out dynamic cross calibration correction on the satellite radiation product.
2. The method for evaluating and cross-scaling the quality of a water-based satellite atmospheric roof radiation product according to claim 1, wherein: in the step S1, when determining the constraint condition of illumination-observation geometry, the illumination-observation geometry is represented by solar zenith angle, satellite zenith angle and relative azimuth angle, and the maximum value and the minimum value of the difference in pixel illumination-observation geometry allowed by multi-source satellite quasi-synchronous observation are obtained as the constraint condition of illumination-observation geometry for judging quasi-synchronous observation under the premise of meeting the IOCCG quality target by using radiation transmission simulation.
3. The method for evaluating and cross-scaling the quality of a water-based satellite atmospheric roof radiation product according to claim 1, wherein: in the step S1, when the hydrological constraint condition is determined, based on the hydrological data, a relationship diagram of the interplanetary deviation of the severe hydrological weather and the radiation products is drawn by combining the multisource satellite quasi-synchrotron radiation products, and the upper limit of the severe hydrological weather is determined as the hydrological constraint condition for judging quasi-synchronous observation by analyzing the statistical characteristics of the interplanetary deviation of the two satellite radiation products under different severe hydrological weather conditions and satisfying the IOCCG quality target.
4. The method for evaluating and cross-scaling the quality of a water-based satellite atmospheric roof radiation product according to claim 1, wherein: in the step S1, when determining the sea surface roughness constraint condition, the maximum space variation coefficient of the satellite image is determined as the constraint condition for judging the sea surface roughness of quasi-synchronous observation by analyzing the quantization relation between the space variation coefficient of the satellite image radiation product and the interstellar deviation of the radiation product on the premise of meeting the IOCCG quality target.
5. The method for evaluating and cross-scaling the quality of a water-based satellite atmospheric roof radiation product according to claim 1, wherein: in the step S1, when the transit time delay is determined, the maximum transit time delay of the multi-source satellite is determined as a constraint condition for determining the transit time delay of the quasi-synchronous observation by analyzing the quantitative relation between the satellite transit time difference and the interstellar deviation of the image radiation product on the premise of meeting the IOCCG quality target.
6. The method for evaluating and cross-scaling the quality of a water-based satellite atmospheric roof radiation product according to claim 1, wherein: in the step S1, when determining the constraint condition of the bright target, the minimum radiation value of all the bright targets in the image is determined as the constraint condition by counting the histogram of the radiation signals of the bright targets of the satellite image.
7. The method for evaluating and cross-scaling the quality of a water-based satellite atmospheric roof radiation product according to claim 1, wherein: in the step S4, when the error of the radiation product of a certain satellite obtained by the error analysis equation is too large, cross calibration is performed, specifically, the radiation data of the reference satellite is taken as a dependent variable, the radiation data of the target satellite is taken as an independent variable, an empirical function relation between the dependent variable and the independent variable is established based on regression analysis, and the empirical function relation is applied to the radiation data of the target satellite, so that the cross calibrated radiation product is obtained.
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 CN112884342A (en) | 2021-06-01 |
| CN112884342B true 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) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113407596B (en) * | 2021-06-23 | 2022-07-15 | 北京华云星地通科技有限公司 | Meteorological remote sensing data precision inspection analysis system based on cross radiometric calibration |
| CN113534203B (en) * | 2021-07-12 | 2022-08-23 | 长光卫星技术股份有限公司 | On-orbit cross-radiation calibration method based on AERONET aerosol data |
| CN114383656B (en) * | 2021-12-22 | 2023-09-29 | 国家卫星气象中心(国家空间天气监测预警中心) | Satellite data quality monitoring method based on cross point double-reference source deviation |
| CN114579655B (en) * | 2022-05-06 | 2022-07-29 | 国家卫星海洋应用中心 | Automatic planning method and device for ocean water color satellite radiometric calibration |
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 kind of seamless orthography production method in the big region of satellite remote-sensing image |
| 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 | 北京空间机电研究所 | Land resource satellite calibration orbit planning and reference load orbit optimization method |
| CN111737850A (en) * | 2020-05-15 | 2020-10-02 | 中国科学院空天信息创新研究院 | Multisource satellite AOD fusion method based on uncertainty on 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 |
| 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 kind of seamless orthography production method in the big region of satellite remote-sensing image |
| 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 | 北京空间机电研究所 | Land resource satellite calibration orbit planning and reference load orbit optimization method |
| CN111737850A (en) * | 2020-05-15 | 2020-10-02 | 中国科学院空天信息创新研究院 | Multisource satellite AOD fusion method based on uncertainty on 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 |
|---|
| "高分二号"卫星相机影像辐射质量评价;徐文等;航天返回与遥感;第36卷(第4期);1-9 * |
| GPM_IMERGE卫星遥感降水数据在巢湖流域的精度评价;程扬等;水土保持研究;第27卷(第5期);188-193 * |
| Integrated fusion of multi-scale polar-orbiting and geostationary satellite observations for the mapping of high spatial and temporal resolution land surface temperature;Penghai Wu等;Remote Sensing of Environment;169-181 * |
| Landsat系列卫星光学遥感器辐射定标方法综述;张志杰;张浩;常玉光;陈正超;;遥感学报(第05期);21-34 * |
| Multi-source high-resolution satellite products in Yangtze Estuary: cross-comparisons and impacts of signal-to-noise ratio and spatial resolution;RUGANG TANG等;Opt. Express;第27卷;6426-6441 * |
| SPARK卫星高光谱数据辐射质量评价;吴兴等;遥感技术与应用;第33卷(第2期);233-240 * |
| VIIRS在中国渤海的遥感反射率产品验证;李军;朱建华;韩冰;高飞;杨安安;赵屹立;;海洋技术学报(第02期);30-36 * |
| 三峡库区CO_2平均柱浓度的卫星和地面观测研究;荣飘;张淳民;刘冬冬;栗彦芬;;三峡生态环境监测(第03期);12-21 * |
| 两种星载微波辐射计被动亮温数据的交叉定标;黄端;邱玉宝;石利娟;武胜利;刘小生;阮永俭;;测绘科学(第01期);139-145 * |
| 传感器交叉辐射定标综述;高彩霞;姜小光;马灵玲;霍红元;;干旱区地理(第01期);141-148 * |
| 光学遥感仪器月球定标技术研究;牛明慧;中国博士学位论文全文数据库基础科学辑（月刊）;20190515(第05期);A007-3 * |
| 北冰洋卫星水色遥感观测能力评价与展望;崔廷伟等;中国海洋大学学报;第51卷(第1期);125-137 * |
| 卫星辐射带高能带电粒子观测数据交叉定标-FY3和NOAA应用示例;王春琴;张贤国;李嘉巍;黄聪;张效信;荆涛;沈国红;张申毅;曹光伟;孙越强;梁金宝;朱光武;韩英;;中国科学:技术科学;20131120(第11期);90-96 * |
| 卫星辐射数据质量分析与研究;马艳红;中国优秀硕士学位论文全文数据库工程科技Ⅱ辑（月刊）(第02期);C028-260 * |
| 卫星高度计定标现状;杨磊;周兴华;徐全军;柯宝贵;穆博;朱琳;;遥感学报;20190525(第03期);38-53 * |
| 基于高光谱资料对FY-1C/1D气象卫星进行交叉定标;漆成莉;胡秀清;张里阳;董超华;;气象学报;20120815(第04期);286-295 * |
| 徐伟进.静止卫星反射波段辐射性能评估和辐射归一化.中国优秀硕士学位论文全文数据库基础科学辑（月刊）.2020,(第01期),A009-128. * |
| 敦煌场地CBERS-02 CCD传感器在轨绝对辐射定标研究;傅俏燕;闵祥军;李杏朝;沙崇漠;李小英;马国强;潘志强;郭毅;李启明;刘国栋;;遥感学报(第04期);3-9 * |
| 无线电掩星反演大气参数误差分析及应用进展;蒋虎, 黄, 严豪键;地球物理学进展(第01期);85-91 * |
| 近海多卫星测高数据联合处理的方法及应用;胡建国, 章传银, 常晓涛;测绘通报;20040125(第01期);1-4 * |
| 风云二号卫星在轨辐射定标技术进展;郭强;陈博洋;张勇;徐娜;冯绚;杨昌军;王新;韩秀珍;李元;胡秀清;张志清;;气象科技进展(第06期);8-14 * |
| 高分一号卫星PMS相机多场地宽动态辐射定标;韩启金;张学文;乔志远;杨磊;潘志强;刘李;;红外与激光工程(第01期);133-139 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112884342A (en) | 2021-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112884342B (en) | Quality evaluation and cross calibration method for water-color satellite atmospheric roof radiation product | |
| Lymburner et al. | Landsat 8: Providing continuity and increased precision for measuring multi-decadal time series of total suspended matter | |
| Quilfen et al. | The ERS scatterometer wind measurement accuracy: Evidence of seasonal and regional biases | |
| Sherwood et al. | Shear velocity estimates on the inner shelf off Grays Harbor, Washington, USA | |
| Stengel et al. | Assimilation of SEVIRI infrared radiances with HIRLAM 4D‐Var | |
| Laloyaux et al. | Impact of scatterometer surface wind data in the ECMWF coupled assimilation system | |
| Xie et al. | Calculating NDVI for Landsat7-ETM data after atmospheric correction using 6S model: A case study in Zhangye city, China | |
| Vivone et al. | Atmospheric boundary layer height estimation from aerosol lidar: a new approach based on morphological image processing techniques | |
| Stengel et al. | The impact of cloud‐affected IR radiances on forecast accuracy of a limited‐area NWP model | |
| Zhang et al. | Assimilation of passive microwave streamflow signals for improving flood forecasting: A first study in Cubango River Basin, Africa | |
| Chen et al. | A simple atmospheric correction algorithm for MODIS in shallow turbid waters: A case study in Taihu Lake | |
| Zheng et al. | Characterizing drought events occurred in the Yangtze River Basin from 1979 to 2017 by reconstructing water storage anomalies based on GRACE and meteorological data | |
| Treadon et al. | Global and regional moisture analyses at NCEP | |
| Zhuang et al. | The evaluation and downscaling‐calibration of IMERG precipitation products at sub‐daily scales over a metropolitan region | |
| Zhang et al. | Assimilation of remote sensing observations into a sediment transport model of China’s largest freshwater lake: spatial and temporal effects | |
| Peres et al. | Analyses of the positive bias of remotely sensed SST retrievals in the coastal waters of Rio de Janeiro | |
| Wei et al. | On accuracy of SAR wind speed retrieval in coastal area | |
| Cotton et al. | Assimilation of SMOS L‐band wind speeds: impact on Met Office global NWP and tropical cyclone predictions | |
| Liu et al. | Image-based recognition and processing system for monitoring water levels in an irrigation and drainage channel | |
| Chen et al. | A transferable turbidity estimation method for estimating clear-sky solar irradiance | |
| Niedzielski et al. | El Niño's impact on the probability distribution of sea level anomaly fields. | |
| Nahmani et al. | Sensitivity of GPS tropospheric estimates to mesoscale convective systems in West Africa | |
| Gilles et al. | Topography and recent winter rainfall regime change in temperate Western European areas: a case study in the Rhine‐Meuse basin | |
| Lv et al. | Spatial and temporal variability of precipitation in the context of climate change: A case study of the Upper Yellow River Basin, China | |
| Tang et al. | A comparative evaluation of gauge-satellite-based merging products over multiregional complex terrain basin |
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 |