CN116735008B - Calibration method and device for infrared cross radiation, electronic equipment and medium - Google Patents
Calibration method and device for infrared cross radiation, electronic equipment and medium Download PDFInfo
- Publication number
- CN116735008B CN116735008B CN202311009259.8A CN202311009259A CN116735008B CN 116735008 B CN116735008 B CN 116735008B CN 202311009259 A CN202311009259 A CN 202311009259A CN 116735008 B CN116735008 B CN 116735008B
- Authority
- CN
- China
- Prior art keywords
- load
- target area
- reference load
- observation
- data
- 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 49
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000012545 processing Methods 0.000 claims abstract description 56
- 210000001747 pupil Anatomy 0.000 claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 8
- 238000012937 correction Methods 0.000 claims description 35
- 238000001228 spectrum Methods 0.000 claims description 24
- 230000005540 biological transmission Effects 0.000 claims description 10
- 230000003595 spectral effect Effects 0.000 claims description 9
- 230000003287 optical effect Effects 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 3
- 230000008685 targeting Effects 0.000 claims 1
- 238000004088 simulation Methods 0.000 description 7
- 238000005316 response function Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 240000004282 Grewia occidentalis Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/80—Calibration
- G01J5/802—Calibration by correcting for emissivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/007—Radiation pyrometry, e.g. infrared or optical thermometry for earth observation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/48—Thermography; Techniques using wholly visual means
-
- 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 provides a calibration method, device, electronic equipment and medium of infrared cross radiation, and relates to the technical field of remote sensing. The scaling method comprises the following steps: determining a target area and a reference load, wherein the reference load is a satellite load matched with a thermal infrared spectrum section of a load to be calibrated; acquiring a matched image pair of a target area passing through a reference load and a load to be calibrated, wherein the matched image pair is a thermal infrared image pair of the target area obtained by matching according to preset conditions; performing data processing on the matched image pair to obtain equivalent radiance observation data of the target area at the reference load entrance pupil; correcting the equivalent radiance observation data; calculating to obtain a cross calibration coefficient according to the corrected equivalent radiance observation data; the uncertainty of cross calibration is obtained through calculation according to the cross calibration coefficient, and the problems of few calibration samples, single target type and low cross frequency caused by strict angle matching in the prior art are solved.
Description
Technical Field
The invention relates to the technical field of remote sensing, in particular to a calibration method, device, electronic equipment and medium of infrared cross radiation.
Background
Thermal infrared remote sensing is a main means for acquiring global and regional target infrared radiation characteristics (earth surface emissivity and earth surface temperature), and high-precision load calibration is a key ring for ensuring the accuracy and stability of thermal infrared load acquired data, and the current thermal infrared on-orbit radiation calibration method comprises on-board calibration, field replacement radiation calibration and cross calibration.
The cross calibration is a field-free calibration technology, and the principle is that when an in-orbit satellite remote sensor to be calibrated and a satellite remote sensor with a better calibration result (generally, a satellite remote sensor with an on-board calibration system) are observed in the same target area, the calibration of the satellite remote sensor to be calibrated can be realized by comparing the measured values of the in-orbit satellite remote sensor to be calibrated and the satellite remote sensor with a better calibration result. Compared with the field replacement calibration technology, the method can calibrate the satellite data of multiple remote sensors without establishing a ground correction field. Its advantages are low cost and high-frequency radiation calibration between multiple remote sensors. In recent years, home and abroad scientists use the technology to calibrate various remote sensors. The most widely used space-time point (SNO) cross calibration method in the world at present requires strict space-time spectrum angle matching so as to eliminate the target heat radiation difference caused by the difference of a plurality of observation elements of the space-time spectrum angle. This results in cross-matching points concentrated at the earth's poles, single target types, and observed data satisfying the matching conditions often account for only 0.1% of the total. Therefore, relaxing the constraint conditions in the cross calibration and effectively eliminating the target heat radiation difference caused by the difference of the time-space-spectrum-angle multiple observation factors is a difficulty and a hot spot of the current international calibration field research. However, since the earth emissivity varies with the zenith angle of observation, the satellite observation luminance radiation value varies with the angle of observation. Therefore, the development of the infrared cross radiometric calibration method considering the influence of the directivity of the earth surface emissivity is an urgent need for expanding the existing cross calibration method and effectively improving the cross calibration precision.
Disclosure of Invention
In view of the technical problems, the invention provides a calibration method, a device, electronic equipment and a medium for infrared cross radiation.
A first aspect of an embodiment of the present invention provides a method for calibrating infrared cross radiation, including: determining a target area and a reference load, wherein the reference load is a satellite load matched with a thermal infrared spectrum section of a load to be calibrated; acquiring a matched image pair of a target area passing through a reference load and a load to be calibrated, wherein the matched image pair is a thermal infrared image pair of the target area obtained by matching according to preset conditions; performing data processing on the matched image pair to obtain equivalent radiance observation data of the target area at the reference load entrance pupil; correcting the equivalent radiance observation data; calculating to obtain a cross calibration coefficient according to the corrected equivalent radiance observation data; and calculating the uncertainty of the cross scaling according to the cross scaling coefficient.
In an embodiment of the present invention, performing data processing on a matching image pair to obtain equivalent radiance observation data of a target area at a reference load entrance pupil includes: data processing is carried out on the matched image pair, and a gray scale count value of each pixel point in the matched image pair and a scaling coefficient of a reference load are obtained; calculating an average gray-scale count value of the matched image pair based on the gray-scale count value of each pixel point in the matched image pair; and calculating equivalent radiance observation data of the target area at the reference load entrance pupil based on the average gray scale count value and the scaling coefficient of the reference load.
In one embodiment of the present invention, correcting the equivalent radiance observation data includes: performing spectrum correction processing on the equivalent radiance observation data; and performing angle correction processing on the equivalent radiance observation data subjected to the optical correction processing.
In one embodiment of the present invention, performing spectral correction processing on equivalent radiance observation data includes: obtaining observation geometric information when a reference load and a load to be calibrated pass through an border target area, wherein the observation geometric information comprises an observation zenith angle and an observation azimuth angle and a solar zenith angle and a solar azimuth angle when the reference load and the load to be calibrated pass through the border target area; based on the observation geometric information, simulating by using an atmospheric radiation transmission model to obtain a first data set of equivalent radiance of a channel corresponding to a reference load and a load to be calibrated; calculating to obtain a spectrum matching coefficient based on the first data set of equivalent radiation brightness; and performing spectrum correction processing on the equivalent radiance first data set according to the spectrum matching coefficient.
In an embodiment of the present invention, performing angle correction processing on equivalent radiance observation data after the optical correction processing includes: acquiring the earth surface emissivity of a target area; based on the earth surface emissivity, establishing a relation that the earth surface emissivity changes along with the observation angle; based on the relation of the earth surface emissivity changing along with the observation angle and the observation geometric information, simulating by using an atmospheric radiation transmission model to obtain a second data set of equivalent radiance of a channel corresponding to the reference load and the load to be calibrated; calculating to obtain an angle matching coefficient based on the second data set of the equivalent radiation brightness; and performing angle correction processing on the equivalent radiance second data set according to the angle matching coefficient.
In an embodiment of the present invention, obtaining a matching image of a target area where a reference load and a load to be scaled pass through includes: acquiring thermal infrared images and cloud product data of a target area passing through a reference load and a load to be calibrated; acquiring cloud amount data above a target area; based on cloud product data and cloud amount data, a sunny non-cloud matching image pair of a reference load and a target area of which the load to be calibrated passes through is screened out from the thermal infrared image; and carrying out time matching processing on the cloud-free matching image pair on a sunny day, wherein the time difference between the reference load and the target area of the load to be calibrated passing through is smaller than a preset time period.
In an embodiment of the present invention, data processing is performed on the matched image to obtain equivalent radiance observation data of the target area at the reference load entrance pupil, and the method further includes: and performing space matching processing on the equivalent radiance observation data by adopting an inverse distance weighting algorithm.
A second aspect of the embodiments of the present invention provides an infrared cross radiation scaling apparatus, comprising: the determining module is used for determining a target area and a reference load, wherein the reference load is a satellite load matched with a thermal infrared spectrum band of the load to be calibrated; the acquisition module is used for acquiring a matched image pair of the reference load and the target area passing through the to-be-calibrated load, wherein the matched image is a thermal infrared scanning image of the reference load and the to-be-calibrated load to the target area; the processing module is used for carrying out data processing on the matched image to obtain equivalent radiance observation data of the target area at the reference load entrance pupil; the correction module is used for correcting the equivalent radiance observation data; the first calculation module is used for calculating to obtain a cross calibration coefficient according to the corrected equivalent radiance observation data; and the second calculation module is used for calculating the uncertainty of the cross scaling according to the cross scaling coefficient.
A third aspect of an embodiment of the present invention provides an electronic device, including: one or more processors; and a storage device for storing one or more programs, wherein the one or more programs, when executed by the one or more processors, cause the one or more processors to perform the method of calibrating infrared cross-radiation described above.
A fourth aspect of an embodiment of the present invention provides a computer readable storage medium having stored thereon executable instructions that, when executed by a processor, cause the processor to perform the above-described method of scaling infrared cross-radiation.
The method for calibrating infrared cross radiation provided by the embodiment of the invention has at least the following beneficial effects:
according to the calibration method for infrared cross radiation provided by the embodiment of the invention, through correcting the space-spectrum-angle of equivalent radiation brightness observation data of the target area at the reference load entrance pupil, the influence of the target emissivity directivity on cross calibration is effectively eliminated, the problems of less calibration sample, single target type and low cross frequency caused by strict angle matching in the prior art are broken, and the on-orbit radiation cross calibration precision and frequency of the optical sensor are favorably improved.
According to the infrared cross radiation calibration method provided by the embodiment of the invention, the cross sample matching points are expanded from the two poles of the earth to the middle and low latitude, so that thermal infrared load multi-level cross calibration of widening angle constraint conditions is realized, load wide dynamic sample data are acquired as much as possible, and calibration uncertainty caused by load nonlinear response is effectively reduced.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are only some embodiments of the invention, and other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 schematically shows a flow chart of a method for calibrating infrared cross-radiation according to an embodiment of the invention.
Fig. 2 schematically shows a block diagram of an infrared cross-radiation scaling device according to an embodiment of the invention.
Detailed Description
The present invention will be further described in detail below with reference to specific embodiments and with reference to the accompanying drawings, in order to make the objects, technical solutions and advantages of the present invention more apparent. It will be apparent that the described embodiments are some, but not all, embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention. The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The terms "comprises," "comprising," and/or the like, as used herein, specify the presence of stated features, steps, operations, and/or components, but do not preclude the presence or addition of one or more other features, steps, operations, or components.
In the present invention, unless explicitly specified and limited otherwise, the terms "mounted," "connected," "secured," and the like are to be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally formed; may be mechanically connected, may be electrically connected or may communicate with each other; can be directly connected or indirectly connected through an intermediate medium, and can be communicated with the inside of two elements or the interaction relationship of the two elements. The specific meaning of the above terms in the present invention can be understood by those of ordinary skill in the art according to the specific circumstances.
In the description of the present invention, it should be understood that the terms "longitudinal," "length," "circumferential," "front," "rear," "left," "right," "top," "bottom," "inner," "outer," and the like indicate an orientation or a positional relationship based on that shown in the drawings, merely for convenience in describing the present invention and simplifying the description, and do not indicate or imply that the subsystem or element in question must have a specific orientation, be configured and operated in a specific orientation, and thus should not be construed as limiting the present invention.
Like elements are denoted by like or similar reference numerals throughout the drawings. Conventional structures or constructions will be omitted when they may cause confusion in the understanding of the invention. And the shape, size and position relation of each component in the figure do not reflect the actual size, proportion and actual position relation.
Similarly, in the description of exemplary embodiments of the invention above, various features of the invention are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the various inventive aspects. The description of the terms "one embodiment," "some embodiments," "example," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the present invention, the meaning of "plurality" means at least two, for example, two, three, etc., unless specifically defined otherwise.
Fig. 1 schematically shows a flow chart of a method for calibrating infrared cross-radiation according to an embodiment of the invention.
As shown in fig. 1, the method for calibrating infrared cross radiation provided by the embodiment of the invention may include:
s1, determining a target area and a reference load, wherein the reference load is a satellite load matched with a thermal infrared spectrum band of a load to be calibrated.
Based on thermal infrared cloud-free images and atmospheric analysis data such as a long-time sequence MODIS (moderate resolution imaging spectrometer), landsat8/OLI (Operational Land Imager land imager) and the like, a sliding window method is utilized to carry out spatial filtering on the images in the global scope, and uniform areas with standard deviation of bright temperature superior to 0.3K, area larger than 3km multiplied by 3km and good atmospheric permeability are screened out.
For each selected target area, selecting satellite load with similar or same thermal infrared spectrum band and higher calibration precision as reference load.
S2, a matched image pair of the reference load and the target area of the to-be-scaled load passing border is obtained.
S2, obtaining a matched image pair of a reference load and a target area of a to-be-scaled load passing border comprises the following steps:
s21, acquiring thermal infrared images and cloud product data of the target area where the reference load and the load to be calibrated pass through.
S22, cloud cover data above the target area are acquired.
S23, based on cloud product data and cloud amount data, a clear-day cloud-free matching image pair of a reference load and a target area of which the load to be calibrated passes through is screened out from the thermal infrared image.
And S24, performing time matching processing on the cloud-free matching image pair on a sunny day, wherein the time difference between the reference load and the target area of the to-be-scaled load passing through is less than 10 minutes.
The influence of the change of the infrared radiation characteristics of the targets caused by the observation time difference is reduced as much as possible by screening the matched image pairs which are cloudless on a sunny day and have the time difference of the reference load and the target area of the border crossing of the load to be calibrated smaller than 10 minutes.
And S3, carrying out data processing on the matched image pair to obtain equivalent radiance observation data of the target area at the reference load entrance pupil.
S3, performing data processing on the matched image pair to obtain equivalent radiance observation data of the target area at the reference load entrance pupil, wherein the step of obtaining the equivalent radiance observation data comprises the following steps:
s31, carrying out data processing on the matched image pair to obtain the gray count value of each pixel point in the matched image pair and the scaling coefficient of the reference load.
For a target area, extracting gray count values of each pixel image of the target area from a matched image pair obtained by a reference load and a load to be calibrated, extracting longitude and latitude information of four corners of the image from a data auxiliary file, calculating an observation zenith angle and an observation azimuth angle of a sensor corresponding to a calibration field according to the position of a large-area uniform calibration field in the image, and calculating the observation zenith angle and the azimuth angle of the sun by combining the date and time of a satellite passing through the calibration field; in addition, the scaling coefficient of the corresponding channel of the reference load is extracted from the data auxiliary file.
Because the reference load and the load to be calibrated generally have different spatial resolutions, the two load observation data are spatially aggregated and matched by adopting an inverse distance weighting algorithm before cross calibration, the load data with higher spatial resolution are registered to the load data with lower spatial resolution, and the uncertainty caused by inconsistent spatial sizes is reduced.
The calculation formula of the inverse distance weighting algorithm is as follows:
wherein,a pixel gray count value representing the pair of aggregated lower spatial resolution matching images;DN j representing higher spatial resolution ofjA pixel gray count value of each pixel;1／d j representing higher spatial resolution pixelsjReciprocal of the distance sampling point;DNrepresenting the pixel gray count value.
S32, calculating an average gray scale count value of the matched image pair based on the gray scale count value of each pixel point in the matched image pair.
And carrying out statistical analysis on the gray count value of each channel target area matching image pair of the reference load and the load to be scaled, and calculating the average value of the target area matching image pair gray count value.
S33, calculating equivalent radiance observation data of the target area at the reference load entrance pupil based on the average gray scale count value and the scaling coefficient of the reference load.
Calculating to obtain equivalent radiance observation values of corresponding channels at the reference load entrance pupil of the target region by using a radiometric calibration formula based on the average value of the target region matching image on the gray count value and the calibration coefficient of the reference loadL j 。
Wherein,representing a reference load pathjThe obtained large-area uniform scene count average value,Gain j representing a reference load pathjIs used for the scaling of the gain of (a),Bias j representing a reference load pathjIs provided.
S4, correcting the equivalent radiance observation data.
S4, correcting the equivalent radiance observation data comprises the following steps:
s41, performing spectrum correction processing on the equivalent radiance observation data.
S42, angle correction processing is carried out on the equivalent radiance observation data after the spectrum correction processing.
Wherein, S41, performing the spectrum correction processing on the equivalent radiance observation data includes:
s411, obtaining the observation geometric information when the reference load and the load to be calibrated pass through the border target area, wherein the observation geometric information comprises the observation zenith angle and the observation azimuth angle and the solar zenith angle and the solar azimuth angle when the reference load and the load to be calibrated pass through the border target area.
Four-corner longitude and latitude information of an image is extracted from the data auxiliary file, and according to the positions of large-area uniform calibration sites in the image, the observed zenith angle and the observed azimuth angle of the sensor corresponding to the calibration sites are calculated; and calculating to obtain the observed zenith angle and azimuth angle of the sun by combining the date and time of the satellite passing calibration field.
And S412, based on the observation geometric information, simulating by using an atmospheric radiation transmission model to obtain a first data set of equivalent radiance of the channels corresponding to the reference load and the load to be calibrated.
For a reference load and a load to be calibrated matched image, collecting ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric analysis data with the date spatial resolution of 0.25 degrees and the time resolution of 1 hour of a target area passing through the reference load and the load to be calibrated, extracting an atmospheric temperature wet pressure profile at the nearest moment, selecting a spectrum library or a target emissivity spectrum under an actually measured 0-degree observation angle, setting the surface temperature range to 290-330K, combining information such as an observation geometry (an observation zenith angle and an observation azimuth angle) between the target area and the reference load, a load spectrum response function and the like, setting an input file for driving an atmospheric radiation transmission model MODTRA (atmospheric radiation transmission mode) to operate, and obtaining a first data set of equivalent radiation brightness of a corresponding channel of the load to be calibrated and the reference load through analog calculation under the reference load observation geometry.
S413, calculating to obtain a spectrum matching coefficient based on the first data set of the equivalent radiance.
And calculating a spectrum matching coefficient by a least square linear fitting method based on the first data set of equivalent radiance.
And S414, performing spectrum correction processing on the first data set of the equivalent radiance according to the spectrum matching coefficient.
And multiplying the spectral matching coefficient serving as a multiplicative factor by an equivalent radiance observation value at the entrance pupil of the reference load, and correcting the spectral matching coefficient to a corresponding channel of the load to be calibrated.
Wherein,representing the equivalent radiation brightness value at the entrance pupil of the load to be calibrated obtained through simulation; />Representing the spectrum radiation brightness at the entrance pupil of the load to be calibrated obtained through simulation; />Representing a spectral response function of the load to be scaled; />Representing a reference load spectral response function; />Representing the equivalent radiance value at the reference load entrance pupil obtained by simulation; />Representing spectral radiance at the entrance pupil of the reference load obtained by simulation; />Representing the spectral matching coefficients.
Wherein, S42, the angle correction processing of the equivalent radiance observation data after the spectrum correction processing includes:
s421, the surface emissivity of the target area is obtained.
And downloading and screening the surface temperature emissivity products of the MYD21 in the cloud-free long-time sequence on a sunny day according to each target area, and extracting the surface emissivity of the target area within the range of 0-65 degrees of observation zenith angles.
S422, based on the earth surface emissivity, establishing a relation that the earth surface emissivity changes along with the observation angle.
Carrying out statistical analysis on the earth surface emissivity in the zenith angle range of 0-65 degrees of observation of the target area, calculating the average value of the emissivity of the target area, analyzing the relation of the average value of the emissivity of the target along with the change of the observation angle, establishing the relation of the average value of the emissivity of the target and the change of the observation angle, and constructing the directivity model of the emissivity of the target area.
S423, based on the relation of the earth surface emissivity changing along with the observation angle and the observation geometric information, simulating by using an atmospheric radiation transmission model to obtain a second data set of equivalent radiance of the channels corresponding to the reference load and the load to be calibrated.
Aiming at the image matched with the reference load and the load to be calibrated, based on the extracted target emissivity under the observation angle of 0 degree, the built target emissivity directivity model is utilized to correct the image to the observation geometrical conditions of the reference load and the load to be calibrated. On the basis, combining an atmospheric temperature wet pressure profile, setting the surface temperature range to 290-330K, and obtaining a second data set of equivalent radiance of each channel corresponding to the reference load under the reference load and the load to be calibrated under the observation geometry by using an atmospheric radiation transmission model MODTRA and a reference load spectral response function through simulation calculation under the reference load and the load to be calibrated.
S424, calculating an angle matching coefficient based on the second data set of the equivalent radiance.
And calculating an angle matching coefficient by a least square linear fitting method based on the simulated equivalent radiance second data set of the corresponding channels of the load to be calibrated and the reference load.
Wherein,representing the angle matching coefficient, ++>Representing the top radiance of the atmosphere layer corresponding to the reference load channel under the reference load observation geometry obtained through simulation; />And (5) representing the simulation to obtain the top radiance of the atmosphere layer corresponding to the reference load channel under the observation geometry of the load to be calibrated.
S425, performing angle correction processing on the second data set of the equivalent radiance according to the angle matching coefficient.
And multiplying the angle matching coefficient serving as a multiplicative factor by an equivalent radiance observation value at a reference load entrance pupil, and correcting the difference of the load radiation value caused by the target directivity and the observation geometric difference.
S5, calculating to obtain a cross scaling coefficient according to the corrected equivalent radiance observation data.
After time, space and angle of the observed data of the reference load and the load to be calibrated are matched, the slope and intercept of a fitting straight line are obtained by utilizing an equivalent radiance observed value at the entrance pupil of the reference load and an image pixel counting average value of a target area of the load to be calibrated and adopting a least square linear fitting method, and the slope and intercept are used as the gain and bias of cross calibration, so that the cross calibration coefficient is calculated.
S6, calculating according to the cross scaling coefficient to obtain the uncertainty of the cross scaling.
The method comprises the steps of analyzing uncertainty factors in an infrared cross radiation calibration method considering the influence of the earth surface emissivity directivity, wherein the uncertainty factors comprise uncertainty caused by reference remote sensing load self calibration errors, spectrum matching coefficients, angle matching coefficients and image space matching, analyzing the influence of each factor on a cross calibration result by adopting a control variable method, and synthesizing uncertainty components of the uncertainty factors based on an error transfer theory to obtain the total cross calibration uncertainty.
According to the calibration method for infrared cross radiation provided by the embodiment of the invention, through correcting the space-spectrum-angle of equivalent radiation brightness observation data of the target area at the reference load entrance pupil, the influence of the target emissivity directivity on cross calibration is effectively eliminated, the problems of less calibration sample, single target type and low cross frequency caused by strict angle matching in the prior art are broken, and the on-orbit radiation cross calibration precision and frequency of the optical sensor are favorably improved.
According to the infrared cross radiation calibration method provided by the embodiment of the invention, the cross sample matching points are expanded from the two poles of the earth to the middle and low latitude, so that thermal infrared load multi-level cross calibration of widening angle constraint conditions is realized, load wide dynamic sample data are acquired as much as possible, and calibration uncertainty caused by load nonlinear response is effectively reduced.
Fig. 2 schematically shows a block diagram of an infrared cross-radiation scaling device according to an embodiment of the invention.
As shown in fig. 2, the device for calibrating infrared cross radiation provided by the embodiment of the invention may include: a determining module 201, an obtaining module 202, a processing module 203, a correcting module 204, a first calculating module 205 and a second calculating module 206.
The determining module 201 is configured to determine a target area and a reference load, where the reference load is a satellite load that matches a thermal infrared spectrum band of the load to be scaled.
The obtaining module 202 is configured to obtain a matching image of the reference load and the target area where the load to be calibrated passes through, where the matching image pair is a thermal infrared image pair obtained by matching according to a preset condition.
The processing module 203 is configured to perform data processing on the matching image, so as to obtain equivalent radiance observation data of the target area at the reference load entrance pupil.
The correction module 204 is used for performing correction processing on the equivalent radiance observation data.
The first calculation module 205 is configured to calculate a cross scaling coefficient according to the corrected equivalent radiance observation data.
The second calculation module 206 is configured to calculate the uncertainty of the cross scaling based on the cross scaling coefficients.
According to an embodiment of the present invention, any of the determining module 201, the obtaining module 202, the processing module 203, the correcting module 204, the first calculating module 205, and the second calculating module 206 may be combined in one module to be implemented, or any one of the modules may be split into a plurality of modules, or at least part of the functions of one or more of the modules may be combined with at least part of the functions of other modules to be implemented in one module. According to an embodiment of the invention, at least one of the determining module 201, the obtaining module 202, the processing module 203, the correcting module 204, the first computing module 205 and the second computing module 206 may be implemented at least partly as hardware circuitry, for example as a Field Programmable Gate Array (FPGA), a Programmable Logic Array (PLA), a system on chip, a system on a substrate, a system on a package, an Application Specific Integrated Circuit (ASIC), or as hardware or firmware by any other reasonable way of integrating or packaging the circuitry, or as any one of or a suitable combination of three of software, hardware and firmware. Alternatively, at least one of the preprocessing module, the extraction module, the generation module, the detection module and the measurement module may be at least partially implemented as a computer program module, which when executed may perform the respective functions.
It should be noted that, in the embodiment of the present invention, the infrared cross radiation calibration device corresponds to the infrared cross radiation calibration method in the embodiment of the present invention, and the specific implementation details and the brought technical effects are the same, which is not repeated here.
The embodiment of the invention also provides electronic equipment, which comprises: one or more processors; and a storage device for storing one or more programs, wherein the one or more programs, when executed by the one or more processors, cause the one or more processors to perform the method of calibrating infrared cross-radiation described above.
Embodiments of the present invention also provide a computer readable storage medium having stored thereon executable instructions that, when executed by a processor, cause the processor to perform the above-described method of calibrating infrared cross-radiation.
While the invention has been illustrated and described in detail in the drawings and foregoing description, such illustration and description are to be considered illustrative or exemplary and not restrictive.
Those skilled in the art will appreciate that the features recited in the various embodiments of the invention can be combined and/or combined in a wide variety of ways, even if such combinations or combinations are not explicitly recited in the present invention. In particular, the features recited in the various embodiments of the invention can be combined and/or combined in various ways without departing from the spirit and teachings of the invention. All such combinations and/or combinations fall within the scope of the invention.
While the invention has been shown and described with reference to certain exemplary embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended embodiments and equivalents thereof. Thus, the scope of the invention should not be limited to the embodiments described above, but should be determined not only by the appended embodiments, but also by equivalents of the appended embodiments.
Claims (8)
1. A method for calibrating infrared cross radiation, comprising:
determining a target area and a reference load, wherein the reference load is a satellite load matched with a thermal infrared spectrum band of a load to be calibrated;
acquiring a matched image pair of the target area, which is obtained by matching according to preset conditions, of the reference load and the load to be calibrated;
performing data processing on the matched image pair to obtain equivalent radiance observation data of the target area at the reference load entrance pupil;
performing spectrum correction processing on the equivalent radiance observation data;
performing angle correction processing on the equivalent radiance observation data after the optical correction processing, wherein the angle correction processing on the equivalent radiance observation data after the optical correction processing comprises the following steps: acquiring the earth surface emissivity of the target area; based on the earth surface emissivity, establishing a relation of the earth surface emissivity changing along with an observation angle; based on the relation of the earth surface emissivity changing along with the observation angle and the observation geometric information, simulating by using an atmospheric radiation transmission model to obtain a second data set of equivalent radiance of the channels corresponding to the reference load and the load to be calibrated; calculating to obtain an angle matching coefficient based on the second data set of the equivalent radiance; performing angle correction processing on the second data set of equivalent radiance according to the angle matching coefficient;
calculating to obtain a cross calibration coefficient according to the corrected equivalent radiance observation data;
and calculating to obtain the uncertainty of the cross scaling according to the cross scaling coefficient.
2. The method of calibrating infrared cross-radiation according to claim 1, wherein the data processing the matched image pair to obtain equivalent radiance observation data of the target region at the reference load entrance pupil comprises:
performing data processing on the matched image pair to obtain a gray count value of each pixel point in the matched image pair and a scaling coefficient of the reference load;
calculating an average gray-scale count value of the matched image pair based on the gray-scale count value of each pixel point in the matched image pair;
and calculating equivalent radiance observation data of the target area at the reference load entrance pupil based on the average gray scale count value and the scaling coefficient of the reference load.
3. The method of calibrating infrared cross-radiation of claim 1, wherein the performing a spectral correction process on the equivalent radiance observation data comprises:
obtaining observation geometric information when the reference load and the load to be calibrated pass through the target area, wherein the observation geometric information comprises an observation zenith angle and an observation azimuth angle and a solar zenith angle and a solar azimuth angle when the reference load and the load to be calibrated pass through the target area;
based on the observation geometric information, simulating and obtaining a first data set of equivalent radiation brightness of the channel corresponding to the reference load and the load to be calibrated by using an atmospheric radiation transmission model;
calculating to obtain a spectrum matching coefficient based on the first data set of equivalent radiation brightness;
and carrying out spectrum correction processing on the first data set of equivalent radiance according to the spectrum matching coefficient.
4. The method of calibrating infrared cross-radiation of claim 1, wherein the acquiring a matched image pair of the reference load and the target area across which the load is to be calibrated comprises:
acquiring thermal infrared images and cloud product data of the target area in which the reference load and the load to be calibrated pass;
acquiring cloud amount data above the target area;
screening out a clear-day cloud-free matching image pair of the target area, which is in transit with the reference load and the load to be calibrated, from the thermal infrared image based on the cloud product data and the cloud amount data;
and carrying out time matching processing on the cloud-free matching image pair on a sunny day, wherein the time difference between the reference load and the target area in which the load to be calibrated passes is smaller than a preset time period.
5. The method for scaling infrared cross radiation according to claim 1, wherein said data processing said matching image to obtain equivalent radiance observation data of said target region at said reference load entrance pupil, further comprises:
and performing space matching processing on the equivalent radiance observation data by adopting an inverse distance weighting algorithm.
6. An infrared cross-radiation targeting device, comprising:
the determining module is used for determining a target area and a reference load, wherein the reference load is a satellite load matched with a thermal infrared spectrum band of the load to be calibrated;
the acquisition module is used for acquiring a matched image pair of the target area, which is obtained by matching according to preset conditions, of the reference load and the load to be calibrated;
the processing module is used for carrying out data processing on the matched image pair to obtain equivalent radiance observation data of the target area at the reference load entrance pupil;
the correction module is used for carrying out spectrum correction processing on the equivalent radiance observation data; performing angle correction processing on the equivalent radiance observation data after the optical correction processing, wherein the angle correction processing on the equivalent radiance observation data after the optical correction processing comprises the following steps: acquiring the earth surface emissivity of the target area; based on the earth surface emissivity, establishing a relation of the earth surface emissivity changing along with an observation angle; based on the relation of the earth surface emissivity changing along with the observation angle and the observation geometric information, simulating by using an atmospheric radiation transmission model to obtain a second data set of equivalent radiance of the channels corresponding to the reference load and the load to be calibrated; calculating to obtain an angle matching coefficient based on the second data set of the equivalent radiance; performing angle correction processing on the second data set of equivalent radiance according to the angle matching coefficient;
the first calculation module is used for calculating to obtain a cross calibration coefficient according to the corrected equivalent radiance observation data;
and the second calculation module is used for calculating the uncertainty of the cross scaling according to the cross scaling coefficient.
7. An electronic device, comprising:
one or more processors;
storage means for storing one or more programs,
wherein the one or more programs, when executed by the one or more processors, cause the one or more processors to perform the method of any of claims 1-5.
8. A computer readable storage medium having stored thereon executable instructions which, when executed by a processor, cause the processor to perform the method according to any of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311009259.8A CN116735008B (en) | 2023-08-11 | 2023-08-11 | Calibration method and device for infrared cross radiation, electronic equipment and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311009259.8A CN116735008B (en) | 2023-08-11 | 2023-08-11 | Calibration method and device for infrared cross radiation, electronic equipment and medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116735008A CN116735008A (en) | 2023-09-12 |
| CN116735008B true CN116735008B (en) | 2023-11-14 |
Family
ID=87915389
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311009259.8A Active CN116735008B (en) | 2023-08-11 | 2023-08-11 | Calibration method and device for infrared cross radiation, electronic equipment and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116735008B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102279393A (en) * | 2011-07-15 | 2011-12-14 | 北京航空航天大学 | Cross radiometric calibration method of hyper-spectral sensor based on multi-spectral sensor |
| CN102540166A (en) * | 2011-12-05 | 2012-07-04 | 北京航空航天大学 | Cross radiation calibration method based on optimization algorithm of hyper-spectral sensor |
| CN103728609A (en) * | 2014-01-16 | 2014-04-16 | 中国科学院地理科学与资源研究所 | Intersected radiometric calibration method for satellite-borne multispectral infrared sensor |
| CN104237869A (en) * | 2014-08-29 | 2014-12-24 | 中国资源卫星应用中心 | Multi-spectral load site-free cross calibration method based on hyperspectral loads |
| CN113970376A (en) * | 2021-09-07 | 2022-01-25 | 中国科学院空天信息创新研究院 | Satellite infrared load calibration method based on ocean area reanalysis data |
| CN115265783A (en) * | 2022-08-01 | 2022-11-01 | 国家卫星海洋应用中心 | Multi-platform satellite-borne point instantaneous cross calibration method and device based on hyperspectral data |
-
2023
- 2023-08-11 CN CN202311009259.8A patent/CN116735008B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102279393A (en) * | 2011-07-15 | 2011-12-14 | 北京航空航天大学 | Cross radiometric calibration method of hyper-spectral sensor based on multi-spectral sensor |
| CN102540166A (en) * | 2011-12-05 | 2012-07-04 | 北京航空航天大学 | Cross radiation calibration method based on optimization algorithm of hyper-spectral sensor |
| CN103728609A (en) * | 2014-01-16 | 2014-04-16 | 中国科学院地理科学与资源研究所 | Intersected radiometric calibration method for satellite-borne multispectral infrared sensor |
| CN104237869A (en) * | 2014-08-29 | 2014-12-24 | 中国资源卫星应用中心 | Multi-spectral load site-free cross calibration method based on hyperspectral loads |
| CN113970376A (en) * | 2021-09-07 | 2022-01-25 | 中国科学院空天信息创新研究院 | Satellite infrared load calibration method based on ocean area reanalysis data |
| CN115265783A (en) * | 2022-08-01 | 2022-11-01 | 国家卫星海洋应用中心 | Multi-platform satellite-borne point instantaneous cross calibration method and device based on hyperspectral data |
Non-Patent Citations (5)
| Title |
|---|
| 传感器交叉辐射定标综述;高彩霞;姜小光;马灵玲;霍红元;;干旱区地理(第01期);第139-146页 * |
| 余明.《遥感影像的城市热环境综合信息图谱研究》.测绘出版社,第47页. * |
| 孙建国.《黄土高原土地退化与植被动态的遥感分析》.中国环境出版社,2014,第29页. * |
| 张玉娟等.《遥感数字图像处理》.哈尔滨工程大学出版社,2016,第61页. * |
| 邵振峰.《遥感图像解译》.武汉大学出版社,2022,第152页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116735008A (en) | 2023-09-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102279393B (en) | Cross radiometric calibration method of hyper-spectral sensor based on multi-spectral sensor | |
| Li et al. | An evaluation of the use of atmospheric and BRDF correction to standardize Landsat data | |
| CN102540166B (en) | Cross radiation calibration method based on optimization algorithm of hyper-spectral sensor | |
| CN102901516B (en) | A kind of multispectral image radiation correction method based on absolute radiometric calibration | |
| CN112798013B (en) | Method for verifying on-orbit absolute radiometric calibration result of optical load | |
| CN103438900A (en) | Three-line-array camera image collaborative absolute radiometric calibration and compensation method | |
| CN113177512B (en) | Matching threshold analysis method for cross radiation calibration between stars | |
| CN106887016B (en) | Automatic relative registration method for GF-4 satellite sequence images | |
| Bendig et al. | Solar-induced chlorophyll fluorescence measured from an unmanned aircraft system: sensor etaloning and platform motion correction | |
| CN104977024A (en) | Solar blind ultraviolet remote sensing camera absolute radiometric calibration coefficient in-orbit correction method | |
| Kusnierek et al. | Challenges in using an analog uncooled microbolometer thermal camera to measure crop temperature | |
| CN108226059B (en) | Satellite hyperspectral CO2On-orbit radiation calibration method for detector | |
| CN111257238A (en) | Detection element relative calibration method based on satellite-borne solar diffusion plate | |
| Valorge et al. | Forty years of experience with SPOT in-flight calibration | |
| Bréon et al. | Measuring the directional variations of land surface reflectance from MODIS | |
| CN116735008B (en) | Calibration method and device for infrared cross radiation, electronic equipment and medium | |
| CN113532652A (en) | Infrared remote sensing sensor absolute calibration method based on buoy and atmospheric reanalysis data | |
| CN115096854A (en) | Tower footing hyperspectral remote sensing synchronous atmospheric correction method supported by artificial target | |
| Liu et al. | Atmospheric absorption ratio algorithm for airborne short-wave infrared hyperspectral imagery spectral calibration based on carbon dioxide and water vapor | |
| CN110702228B (en) | Edge radiation correction method for aviation hyperspectral image | |
| Cariou et al. | Automatic georeferencing of airborne pushbroom scanner images with missing ancillary data using mutual information | |
| CN114296061A (en) | Cross calibration method based on multivariate variable detection and different radiation transmission models | |
| Kouyama et al. | One-year lunar calibration result of Hodoyoshi-1, Moon as an ideal target for small satellite radiometric calibration | |
| CN114494377B (en) | Construction method of directional emissivity model of pixel scale | |
| Langheinrich et al. | Smile correction in the ENMAP ground segment processor: a qualitative analysis |
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 |