CN115112244A - Satellite-borne thermal infrared camera absolute radiometric calibration processing method - Google Patents
Satellite-borne thermal infrared camera absolute radiometric calibration processing method Download PDFInfo
- Publication number
- CN115112244A CN115112244A CN202210425930.6A CN202210425930A CN115112244A CN 115112244 A CN115112244 A CN 115112244A CN 202210425930 A CN202210425930 A CN 202210425930A CN 115112244 A CN115112244 A CN 115112244A
- Authority
- CN
- China
- Prior art keywords
- satellite
- black body
- coefficient
- calibration coefficient
- calibration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 17
- 230000003595 spectral effect Effects 0.000 claims abstract description 64
- 238000012360 testing method Methods 0.000 claims abstract description 60
- 210000001747 pupil Anatomy 0.000 claims abstract description 55
- 238000006243 chemical reaction Methods 0.000 claims abstract description 27
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 15
- 230000005855 radiation Effects 0.000 claims description 29
- 238000000034 method Methods 0.000 claims description 15
- 238000001228 spectrum Methods 0.000 claims description 3
- 241001270131 Agaricus moelleri Species 0.000 claims 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 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
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
Images
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/48—Thermography; Techniques using wholly visual means
-
- 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
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/02—Details
- G01J1/08—Arrangements of light sources specially adapted for photometry standard sources, also using luminescent or radioactive material
-
- 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
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
-
- 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/52—Radiation pyrometry, e.g. infrared or optical thermometry using comparison with reference sources, e.g. disappearing-filament pyrometer
- G01J5/53—Reference sources, e.g. standard lamps; Black bodies
-
- 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
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/02—Details
- G01J1/08—Arrangements of light sources specially adapted for photometry standard sources, also using luminescent or radioactive material
- G01J2001/086—Calibrating drift correction
-
- 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
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J2003/2866—Markers; Calibrating of scan
- G01J2003/2879—Calibrating scan, e.g. Fabry Perot interferometer
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Radiation Pyrometers (AREA)
Abstract
The invention discloses a satellite-borne thermal infrared camera absolute radiometric calibration processing method, which comprises the following steps: obtaining the equivalent spectral radiance of the entrance pupil of the high-temperature black body and the gray level mean value of the image of the high-temperature black body; obtaining the equivalent spectral radiance of the entrance pupil of the low-temperature black body and the gray level mean value of the image of the low-temperature black body; obtaining the slope of the calibration coefficient of the on-satellite black body and the intercept of the calibration coefficient of the on-satellite black body; obtaining the entrance pupil equivalent spectral radiance of the first scene satellite image and the entrance pupil equivalent spectral radiance of the second scene satellite image; obtaining an external field test calibration coefficient slope and an external field test calibration coefficient intercept; obtaining the conversion coefficient of the first and second on-satellite blackbody calibration coefficients to the external field test calibration coefficient; and obtaining the slope of the absolute radiometric calibration coefficient of the thermal infrared spectrum and the intercept of the absolute radiometric calibration coefficient of the thermal infrared spectrum. The invention improves the absolute radiometric calibration precision of the satellite-borne thermal infrared camera.
Description
Technical Field
The invention belongs to the technical field of radiometric calibration, and particularly relates to an absolute radiometric calibration processing method for a satellite-borne thermal infrared camera.
Background
China has rapid aerospace technology development, and a large number of novel high-performance earth observation satellites especially with thermal infrared band detection capability are launched and lifted off. In order to meet the requirements of industries such as environmental protection, China and soil, agriculture, meteorology, disaster reduction and the like on high-precision monitoring of remote sensing of surface temperature, in recent years, the spatial resolution of a thermal infrared spectrum is improved from kilometer scale to ten-meter scale, and the requirement of quantitative application is higher and higher.
High-precision radiometric calibration is the key for ensuring the quantitative application of infrared data, and the common method for absolute radiometric calibration of a satellite-borne thermal infrared camera is to calculate a calibration coefficient by adopting an on-satellite black body, so that the absolute radiometric calibration precision of the infrared camera is insufficient.
Disclosure of Invention
The technical problem solved by the invention is as follows: the method overcomes the defects of the prior art, provides a satellite-borne thermal infrared camera absolute radiometric calibration processing method, and improves the absolute radiometric calibration precision of the satellite-borne thermal infrared camera.
The purpose of the invention is realized by the following technical scheme: a satellite-borne thermal infrared camera absolute radiometric calibration processing method comprises the following steps: obtaining the equivalent spectral radiance of the entrance pupil of the high-temperature black body according to preset calibration data of the high-temperature black body of the satellite-borne thermal infrared camera; obtaining the equivalent spectral radiance of the entrance pupil of the low-temperature black body according to preset calibration data of the low-temperature black body of the satellite-borne thermal infrared camera; obtaining the gray level mean value of the high-temperature black body image according to preset calibration data of the high-temperature black body of the satellite-borne thermal infrared camera
Obtaining the gray level mean value of the low-temperature black body image according to the preset low-temperature black body calibration data of the satellite-borne thermal infrared camera
According to the equivalent spectral radiance L of the entrance pupil of the high-temperature black body eh (T h ) Low temperature black body entrance pupil equivalent spectral radiance L el (T l ) Mean value of gray levels of high-temperature black body image
And the mean value of the gray levels of the low-temperature black body image
Obtaining the slope K of the calibration coefficient of the black body on the satellite s And the intercept C of the calibration coefficient of the on-satellite black body s (ii) a Obtaining a gray mean value of a first scene satellite image according to preset first scene infrared radiation calibration field satellite image data
Obtaining the gray mean value of the second scene satellite image according to the preset second scene infrared radiation calibration field satellite image data
Obtained through a MODTRAN radiation transmission model according to preset outfield test dataEntrance pupil equivalent spectral radiance L of first scene satellite image 1w And the entrance pupil equivalent spectral radiance L of the second scene satellite image 2w (ii) a According to the mean value of the gray scale of the first scene satellite image
Mean value of gray levels of second scene satellite image
Entrance pupil equivalent spectral radiance L of first scene satellite image 1w And the entrance pupil equivalent spectral radiance L of the second scene satellite image 2w Obtaining the slope K of the calibration coefficient of the outfield test w And external field test calibration coefficient intercept C w (ii) a Scaling coefficient slope K according to on-satellite blackbody s The intercept C of the calibration coefficient of the starry black body s External field test calibration coefficient slope K w And external field test calibration coefficient intercept C w Obtaining a conversion coefficient R of the calibration coefficient of the first on-satellite blackbody calibration coefficient to the external field test calibration coefficient 1 And a conversion coefficient R of a second on-satellite blackbody calibration coefficient to an external field test calibration coefficient 2 (ii) a Scaling coefficient slope K according to on-board black body s The intercept C of the calibration coefficient of the starry black body s The conversion coefficient R of the calibration coefficient of the first on-satellite blackbody calibration coefficient to the external field test calibration coefficient 1 And a conversion coefficient R of a second on-satellite blackbody calibration coefficient to an external field test calibration coefficient 2 And obtaining the slope K of the absolute radiometric calibration coefficient of the thermal infrared spectrum section and the intercept C of the absolute radiometric calibration coefficient of the thermal infrared spectrum section.
In the above absolute radiometric calibration processing method for the satellite-borne heat infrared camera, the equivalent spectral radiance of the entrance pupil of the high-temperature black body is obtained by the following formula:
wherein L is eh (T h ) Is the equivalent spectral radiance, L (lambda, T) of the entrance pupil of the high temperature black body h ) Is the radiance, T, of the high temperature black body spectrum h Is the absolute temperature of the high temperature black body, and epsilon isEmissivity, λ is the wavelength, λ 1 、λ 2 The starting and cut-off wavelengths of the imager response, respectively, and R (λ) is the relative spectral response of the imager.
In the absolute radiometric calibration processing method for the satellite-borne heat infrared camera, the low-temperature black body entrance pupil equivalent spectral radiance is obtained through the following formula:
wherein L is el (T l ) Is the low temperature black body entrance pupil equivalent spectral radiance, L (lambda, T) l ) Is the low temperature black body spectral radiance, T l Absolute temperature of low temperature black body, epsilon emissivity, lambda wavelength, lambda 1 、λ 2 The starting and cut-off wavelengths of the imager response, respectively, and R (λ) is the relative spectral response of the imager.
In the absolute radiometric calibration processing method of the satellite-carried infrared camera, the gray level mean value of the high-temperature black body image
Obtained by the following formula:
wherein,
is the mean value of the gray levels of the high-temperature black body image, DN h (i) The gray value of the ith pixel of the high-temperature black body, m is the number of the pixels of the high-temperature black body, and i is the serial number of the pixels of the high-temperature black body.
In the absolute radiometric calibration processing method of the satellite-borne heat infrared camera, the gray level mean value of the low-temperature black body image
Obtained by the following formula:
wherein,
is the mean value of gray levels, DN, of the low-temperature black body image l (j) The gray value of the jth pixel of the low-temperature black body; n is the number of the pixels of the low-temperature black body, and j is the serial number of the pixels of the low-temperature black body.
In the absolute radiometric calibration processing method of the satellite-borne heat infrared camera, the slope K of the calibration coefficient of the black body on the satellite s Obtained by the following formula:
wherein, K s The slope of the scaling factor for the on-satellite black body,
is the mean value of the gray levels of the high-temperature black body image,
is the mean value of gray levels of low-temperature black body image, L el (T l ) Is the low temperature black body entrance pupil equivalent spectral radiance, L eh (T h ) Is the equivalent spectral radiance of the entrance pupil of the high temperature black body.
In the absolute radiometric calibration processing method of the satellite-borne heat infrared camera, the calibration coefficient intercept C of the black body on the satellite s Obtained by the following formula:
wherein, C s Is the intercept of the calibration coefficient of the black body on the satellite,
is high temperature blackThe mean value of the gray levels of the volume image,
is the mean value of gray levels of low-temperature black body image, L el (T l ) Is the low temperature black body entrance pupil equivalent spectral radiance, L eh (T h ) Is the equivalent spectral radiance of the entrance pupil of the high temperature black body.
In the absolute radiation calibration processing method of the satellite-carried heat infrared camera, the slope K of the calibration coefficient of the outfield test w And external field test calibration coefficient intercept C w Are all obtained by the following formula:
wherein,
is the mean value of the gray levels of the first scene satellite image,
is the mean value of the gray levels of the second scene satellite image, L 1w Is the entrance pupil equivalent spectral radiance, L, of the first view satellite image 2w Is the entrance pupil equivalent spectral radiance, K, of the second scene satellite image w Scaling coefficient slope, C, for the outfield test w The coefficient intercept is calibrated for the outfield test.
In the absolute radiometric calibration processing method for the satellite-borne heat infrared camera, the conversion coefficient R of the calibration coefficient from the first satellite blackbody calibration coefficient to the outfield test calibration coefficient 1 Comprises the following steps:
conversion coefficient R of calibration coefficient of second on-satellite blackbody calibration coefficient to external field test calibration coefficient 2 Comprises the following steps:
wherein R is 1 Conversion coefficient, R, of calibration coefficient for first on-satellite blackbody calibration coefficient to external field test 2 Conversion coefficient of calibration coefficient from second on-satellite blackbody calibration coefficient to external field test, K s Scaling coefficient slope, C, for on-satellite blackbodies s Scaling the coefficient intercept, K, for a black body on a satellite w Scaling coefficient slope, C, for the outfield test w The coefficient intercept is calibrated for the outfield test.
In the above method for calibrating absolute radiation of satellite-borne heat infrared camera, the slope K of the calibration coefficient of absolute radiation of thermal infrared spectrum is:
K＝R 1 ·K s ；
the absolute radiometric calibration coefficient intercept C of the thermal infrared spectrum is:
wherein K is the slope of the absolute radiometric calibration coefficient of the thermal infrared spectrum, C is the intercept of the absolute radiometric calibration coefficient of the thermal infrared spectrum, R 1 Conversion coefficient, R, of calibration coefficient for first on-satellite blackbody calibration coefficient to external field test 2 Conversion coefficient of calibration coefficient from second on-satellite blackbody calibration coefficient to external field test, K s Scaling coefficient slope, C, for on-satellite blackbodies s And the coefficient intercept is calibrated for the black body on the satellite.
Compared with the prior art, the invention has the following beneficial effects:
the invention simultaneously uses the on-satellite black body and the external field test to carry out the absolute radiation calibration of the infrared camera, and the absolute calibration precision is improved by obtaining the conversion coefficient from the on-satellite black body calibration coefficient to the external field test calibration coefficient.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
fig. 1 is a flowchart of an absolute radiometric calibration processing method of a satellite-borne thermal infrared camera according to an embodiment of the present invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art. It should be noted that the embodiments and features of the embodiments of the present invention may be combined with each other without conflict. The present invention will be described in detail below with reference to the embodiments with reference to the attached drawings.
Fig. 1 is a flowchart of an absolute radiometric calibration processing method of a satellite-borne thermal infrared camera according to an embodiment of the present invention. As shown in fig. 1, the method comprises the steps of:
the method comprises the following steps: preparing data, namely high-temperature blackbody calibration data, low-temperature blackbody calibration data, satellite image data of two-scene infrared radiation calibration fields and external field test data of the satellite-borne thermal infrared camera; the two sets of infrared radiation calibration field satellite image data comprise first set of infrared radiation calibration field satellite image data and second set of infrared radiation calibration field satellite image data.
Step two: obtaining the equivalent spectral radiance L of the entrance pupil of the high-temperature black body according to the calibration data of the high-temperature black body in the step one eh (T h )。
The equivalent spectral radiance of the entrance pupil of the high-temperature black body is obtained by the following formula:
wherein L is eh (T h ) Is the equivalent spectral radiance (W/m) of the entrance pupil of the high-temperature black body 2 ·sr·μm)，L(λ，T h ) Is the spectral radiance (W/m) of the high temperature black body 2 ·sr·μm)，T h Is the absolute temperature (K) of the high temperature black body, epsilon is the emissivity, lambda is the wavelength (mum), lambda is 1 、λ 2 The starting and cut-off wavelengths (μm) of the imager response, respectively, and R (λ) is the relative spectral response of the imager.
Step three: obtaining the equivalent spectral radiance L of the entrance pupil of the low-temperature black body according to the calibration data of the low-temperature black body in the step I el (T l )。
Low-temperature black body entrance pupil equivalent spectral radiance L el (T l ) Obtained by the following formula:
wherein L is el (T l ) Is the low temperature black body entrance pupil equivalent spectral radiance (W/m) 2 ·sr·μm)，L(λ，T l ) Is the low temperature black body spectrum radiance (W/m) 2 ·sr·μm)，T l Is the absolute temperature (K) of the low temperature black body, epsilon is the emissivity, lambda is the wavelength (mum), lambda 1 、λ 2 The starting and cut-off wavelengths (μm) of the imager response, respectively, and R (λ) is the relative spectral response of the imager.
Step four: obtaining the gray level mean value of the high-temperature black body image according to the high-temperature black body calibration data in the step one
Mean value of gray levels of high-temperature black body image
Obtained by the following formula:
wherein,
is the mean value of the gray levels of the high-temperature black body image, DN h (i) The gray value of the ith pixel of the high-temperature black body, m is the number of the pixels of the high-temperature black body, and i is the serial number of the pixels of the high-temperature black body.
Step five: obtaining the gray level mean value of the low-temperature black body image according to the low-temperature black body calibration data in the step one
Mean value of gray levels of low-temperature black body image
Obtained by the following formula:
wherein,
is the mean value of gray levels, DN, of the low-temperature black body image l (j) The gray value of the jth pixel of the low-temperature black body; n is the number of the pixels of the low-temperature black body, and j is the serial number of the pixels of the low-temperature black body.
Step six: according to the high-temperature black body entrance pupil equivalent spectral radiance L obtained in the step two eh (T h ) And the low-temperature black body entrance pupil equivalent spectral radiance L obtained in the step three el (T l ) The mean value of the gray levels of the high-temperature black body image obtained in the fourth step
The gray average value of the low-temperature black body image obtained in the fifth step
The system of equations is formed:
obtaining the calibration coefficient K of the star blackbody s And C s ：
Wherein L is eh (T h ) The equivalent spectral radiance of the entrance pupil of the high-temperature black body; l is el (T l ) The equivalent spectral radiance of the low-temperature black body entrance pupil;
the gray level mean value of the high-temperature black body image is obtained;
the gray level mean value of the low-temperature black body image is obtained; k s Scaling coefficient slope for the on-satellite blackbody; c s And the coefficient intercept is calibrated for the black body on the satellite.
Step six: respectively obtaining image gray level mean values according to satellite image data of two scene infrared radiation calibration fields in the step one
And
obtaining L through MODTRAN radiation transmission model according to external field test data 1w And L 2w Forming an equation set to obtain an external field test absolute radiometric calibration coefficient K w And C w 。
The system of equations is;
wherein L is 1w The equivalent spectral radiance of the entrance pupil of the first scene satellite image; l is 2w The equivalent spectral radiance of the entrance pupil of the second scene satellite image;
the gray level mean value of the first scene satellite image is obtained;
the gray level mean value of the second scene satellite image is obtained; k is w Scaling coefficient slope for external field test; c w The coefficient intercept is calibrated for the outfield test.
Step seven: obtaining the calibration coefficient K of the star black body according to the step five s And C s And the external field test absolute radiation scaling coefficient K obtained in the sixth step w And C w The following equations are formed:
when the values of the on-satellite calibration data DN and the external field calibration data DN of each pixel are the same, the following conversion relation can be established:
obtaining a conversion coefficient R 1 And R 2 ：
Wherein L is ew Calculating the entrance pupil equivalent spectral radiance for the satellite image by using an external field calibration coefficient; l is es Calculating the entrance pupil equivalent spectral radiance for the satellite image by using the on-satellite black body calibration coefficient; k s Scaling coefficient slope for the on-satellite blackbody; c s Scaling the coefficient intercept for the on-satellite blackbody; k w Scaling coefficient slope for an external field test; c w Calibrating the coefficient intercept for the outfield test; r is 1 Converting the first on-satellite blackbody calibration coefficient into an external field test calibration coefficient; r 2 And testing the conversion coefficient of the calibration coefficient from the second on-satellite blackbody calibration coefficient to the external field.
Step eight: obtaining the calibration coefficient K of the star black body according to the step five s And C s Conversion coefficient R obtained in step seven 1 And R 2 Obtaining an absolute radiometric calibration equation of the thermal infrared camera:
wherein, the scaling factor K, C is:
K＝R 1 ·K s ；
wherein L is image Calculating an entrance pupil equivalent spectral radiance for the infrared image;
the gray level mean value of the infrared image is obtained; k is the absolute calibration coefficient slope calculated by the satellite blackbody and external field test; and C is the absolute calibration coefficient intercept calculated by the satellite blackbody and external field test.
The invention simultaneously uses the on-satellite black body and the external field test to carry out the absolute radiation calibration of the infrared camera, thereby improving the absolute calibration precision. The invention achieves the effect of simultaneous calibration by establishing the conversion relation between the on-satellite blackbody calibration and the field calibration.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to limit the present invention, and those skilled in the art can make variations and modifications of the present invention without departing from the spirit and scope of the present invention by using the methods and technical contents disclosed above.
Claims (10)
1. A satellite-borne thermal infrared camera absolute radiation calibration processing method is characterized by comprising the following steps:
obtaining the equivalent spectral radiance of the entrance pupil of the high-temperature black body according to preset calibration data of the high-temperature black body of the satellite-borne thermal infrared camera;
obtaining the equivalent spectral radiance of the entrance pupil of the low-temperature black body according to preset calibration data of the low-temperature black body of the satellite-borne thermal infrared camera;
obtaining a high-temperature black body image gray level mean value according to preset calibration data of a high-temperature black body of the satellite-borne thermal infrared camera;
obtaining a low-temperature black body image gray level mean value according to preset low-temperature black body calibration data of the satellite-borne thermal infrared camera;
obtaining an on-satellite black body calibration coefficient slope and an on-satellite black body calibration coefficient intercept according to the high-temperature black body entrance pupil equivalent spectral radiance, the low-temperature black body entrance pupil equivalent spectral radiance, the high-temperature black body image gray average value and the low-temperature black body image gray average value;
obtaining a gray average value of a first scene satellite image according to preset first scene infrared radiation calibration field satellite image data; obtaining a gray mean value of a second scene satellite image according to preset second scene infrared radiation calibration field satellite image data; obtaining entrance pupil equivalent spectral radiance of a first scene satellite image and entrance pupil equivalent spectral radiance of a second scene satellite image according to preset external field test data; obtaining an external field test calibration coefficient slope and an external field test calibration coefficient intercept according to the gray average value of the first scene satellite image, the gray average value of the second scene satellite image, the entrance pupil equivalent spectral radiance of the first scene satellite image and the entrance pupil equivalent spectral radiance of the second scene satellite image;
obtaining a conversion coefficient of the first on-satellite blackbody calibration coefficient to the external field test calibration coefficient and a conversion coefficient of the second on-satellite blackbody calibration coefficient to the external field test calibration coefficient according to the on-satellite blackbody calibration coefficient slope, the on-satellite blackbody calibration coefficient intercept, the external field test calibration coefficient slope and the external field test calibration coefficient intercept;
and obtaining the absolute radiation calibration coefficient slope of the thermal infrared spectrum section and the absolute radiation calibration coefficient intercept of the thermal infrared spectrum section according to the satellite blackbody calibration coefficient slope, the satellite blackbody calibration coefficient intercept, the conversion coefficient of the first satellite blackbody calibration coefficient to the external field test calibration coefficient and the conversion coefficient of the second satellite blackbody calibration coefficient to the external field test calibration coefficient.
2. The method for calibrating absolute radiation of the satellite-borne thermal infrared camera according to claim 1, characterized in that: the high-temperature black body entrance pupil equivalent spectral radiance is obtained through the following formula:
wherein L is eh (T h ) Is the equivalent spectral radiance, L (lambda, T) of the entrance pupil of the high temperature black body h ) Is the radiance, T, of the high temperature black body spectrum h Absolute temperature of high temperature black body, epsilon emissivity, lambda wavelength, lambda 1 、λ 2 Respectively, the starting wavelength and the cut-off wavelength of the imager response, and R (λ) is the relative spectral response of the imager.
3. The method for processing the absolute radiation calibration of the satellite-borne thermal infrared camera according to claim 1, characterized in that: the low-temperature black body entrance pupil equivalent spectral radiance is obtained through the following formula:
wherein L is el (T l ) Is the low temperature black body entrance pupil equivalent spectral radiance, L (lambda, T) l ) Is the low temperature black body spectral radiance, T l Absolute temperature of low temperature black body, epsilon emissivity, lambda wavelength, lambda 1 、λ 2 The starting and cut-off wavelengths of the imager response, respectively, and R (λ) is the relative spectral response of the imager.
4. The method for processing the absolute radiation calibration of the satellite-borne thermal infrared camera according to claim 1, characterized in that: mean value of gray levels of high-temperature black body image
Obtained by the following formula:
wherein,
is the mean value of the gray levels of the high-temperature black body image, DN h (i) The gray value of the ith pixel of the high-temperature black body, m is the number of the pixels of the high-temperature black body, and i is the serial number of the pixels of the high-temperature black body.
5. The method for processing the absolute radiation calibration of the satellite-borne thermal infrared camera according to claim 1, characterized in that: mean gray level of low-temperature black body image
Obtained by the following formula:
wherein,
is the mean value of gray levels, DN, of the low-temperature black body image l (j) The gray value of the jth pixel of the low-temperature black body; n is the number of the low-temperature black body pixels,j is the pixel serial number of the low-temperature black body.
6. The method for processing the absolute radiation calibration of the satellite-borne thermal infrared camera according to claim 1, characterized in that: calibration coefficient slope K of black body on satellite s Obtained by the following formula:
wherein, K s The slope of the calibration coefficient is the black body on the satellite,
is the average value of the gray scales of the high-temperature black body image,
is the mean value of gray levels of low-temperature black body image, L el (T l ) Is the low temperature black body entrance pupil equivalent spectral radiance, L eh (T h ) Is the equivalent spectral radiance of the entrance pupil of the high temperature black body.
7. The method for processing the absolute radiation calibration of the satellite-borne thermal infrared camera according to claim 1, characterized in that: intercept C of calibration coefficient of black body on satellite s Obtained by the following formula:
wherein, C s Is the intercept of the calibration coefficient of the black body on the satellite,
is the average value of the gray scales of the high-temperature black body image,
is the mean value of gray levels of low-temperature black body image, L el (T l ) Is the low temperature black body entrance pupil equivalent spectral radiance, L eh (T h ) Is the equivalent spectral radiance of the entrance pupil of the high temperature black body.
8. The method for calibrating absolute radiation of the satellite-borne thermal infrared camera according to claim 1, characterized in that: scaling coefficient slope K of outfield test w And the intercept C of calibration coefficient of external field test w Are all obtained by the following formula:
wherein,
is the mean value of the gray levels of the first scene satellite image,
is the mean value of the gray levels of the second scene satellite image, L 1w Is the entrance pupil equivalent spectral radiance, L, of the first view satellite image 2w Is the entrance pupil equivalent spectral radiance, K, of the second scene satellite image w Scaling coefficient slope, C, for the outfield test w The coefficient intercept is calibrated for the outfield test.
9. The method for processing the absolute radiation calibration of the satellite-borne thermal infrared camera according to claim 1, characterized in that: conversion coefficient R of first on-satellite blackbody calibration coefficient to external field test calibration coefficient 1 Comprises the following steps:
conversion coefficient R of calibration coefficient of second on-satellite blackbody calibration coefficient to external field test calibration coefficient 2 Comprises the following steps:
wherein R is 1 Conversion coefficient, R, of calibration coefficient for first on-satellite blackbody calibration coefficient to external field test 2 Conversion coefficient, K, of the calibration coefficient for the second on-satellite blackbody calibration coefficient to the external field test calibration coefficient s Scaling coefficient slope, C, for on-satellite blackbodies s Scaling the coefficient intercept, K, for a black body on a satellite w Scaling coefficient slope, C, for the outfield test w The coefficient intercept is calibrated for the outfield test.
10. The method for calibrating absolute radiation of the satellite-borne thermal infrared camera according to claim 1, characterized in that: the absolute radiometric calibration coefficient slope K of the thermal infrared spectrum is:
K＝R 1 ·K s ；
the absolute radiometric calibration coefficient intercept C of the thermal infrared spectrum is:
wherein K is the slope of the absolute radiometric calibration coefficient of the thermal infrared spectrum, C is the intercept of the absolute radiometric calibration coefficient of the thermal infrared spectrum, R 1 Conversion coefficient, R, of calibration coefficient for first on-satellite blackbody calibration coefficient to external field test 2 Conversion coefficient, K, of the calibration coefficient for the second on-satellite blackbody calibration coefficient to the external field test calibration coefficient s Scaling coefficient slope, C, for on-satellite blackbodies s And the coefficient intercept is calibrated for the black body on the satellite.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210425930.6A CN115112244B (en) | 2022-04-21 | 2022-04-21 | Absolute radiometric calibration processing method for satellite-borne thermal infrared camera |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210425930.6A CN115112244B (en) | 2022-04-21 | 2022-04-21 | Absolute radiometric calibration processing method for satellite-borne thermal infrared camera |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115112244A true CN115112244A (en) | 2022-09-27 |
| CN115112244B CN115112244B (en) | 2024-07-30 |
Family
ID=83325251
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210425930.6A Active CN115112244B (en) | 2022-04-21 | 2022-04-21 | Absolute radiometric calibration processing method for satellite-borne thermal infrared camera |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115112244B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090272888A1 (en) * | 2008-01-25 | 2009-11-05 | Nugent Paul W | Thermal infrared imaging system and associated methods for radiometric calibration |
| CN104406696A (en) * | 2014-11-27 | 2015-03-11 | 国家海洋环境预报中心 | Calibration method and device for external field spectral radiance of hyperspectral imager |
| US20160349113A1 (en) * | 2015-05-28 | 2016-12-01 | Raytheon Company | Characterization of absolute spectral radiance of an unknown ir source |
| CN109269641A (en) * | 2018-08-21 | 2019-01-25 | 中国科学院合肥物质科学研究院 | A kind of multi-sensor cooperation Calibration Method for drawing No.1 satellite for day |
| RU2685548C1 (en) * | 2018-06-13 | 2019-04-22 | Федеральное государственное унитарное предприятие "Всероссийский научно-исследовательский институт метрологии им. Д.И. Менделеева" | Method of measuring spectral coefficient of radiation of body |
| CN110455417A (en) * | 2019-08-21 | 2019-11-15 | 北京环境特性研究所 | For the quantitative measurment error calibration method of infrared optical system stray radiation |
-
2022
- 2022-04-21 CN CN202210425930.6A patent/CN115112244B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090272888A1 (en) * | 2008-01-25 | 2009-11-05 | Nugent Paul W | Thermal infrared imaging system and associated methods for radiometric calibration |
| CN104406696A (en) * | 2014-11-27 | 2015-03-11 | 国家海洋环境预报中心 | Calibration method and device for external field spectral radiance of hyperspectral imager |
| US20160349113A1 (en) * | 2015-05-28 | 2016-12-01 | Raytheon Company | Characterization of absolute spectral radiance of an unknown ir source |
| RU2685548C1 (en) * | 2018-06-13 | 2019-04-22 | Федеральное государственное унитарное предприятие "Всероссийский научно-исследовательский институт метрологии им. Д.И. Менделеева" | Method of measuring spectral coefficient of radiation of body |
| CN109269641A (en) * | 2018-08-21 | 2019-01-25 | 中国科学院合肥物质科学研究院 | A kind of multi-sensor cooperation Calibration Method for drawing No.1 satellite for day |
| CN110455417A (en) * | 2019-08-21 | 2019-11-15 | 北京环境特性研究所 | For the quantitative measurment error calibration method of infrared optical system stray radiation |
Non-Patent Citations (1)
| Title |
|---|
| 胡秀清, 戎志国, 邱康睦, 张玉香, 张广顺, 黄意玢: "利用青海湖水面辐射校正场对FY-1C气象卫星热红外传感器进行绝对辐射定标", 遥感学报, no. 05, 10 October 2002 (2002-10-10) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115112244B (en) | 2024-07-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108562363B (en) | Method for accurately measuring infrared radiation characteristic transient temperature field | |
| Fujisada et al. | Design and preflight performance of ASTER instrument protoflight model | |
| CN111707376B (en) | Surface temperature inversion method for broadband thermal infrared sensor | |
| CN111323129A (en) | Earth surface temperature inversion method based on broadband thermal infrared image | |
| CN109883957B (en) | MODIS image-based apparent reflectivity model construction method, system and calibration method | |
| CN108830846B (en) | High-resolution full-waveband hyperspectral remote sensing image simulation method | |
| CN111044153B (en) | Nonlinear calibration method and device for infrared spectrum of spectrum correlation system | |
| CN115265803A (en) | Method for measuring infrared spectrum radiation characteristics of remote point target | |
| Wang et al. | A land surface temperature retrieval method for UAV broadband thermal imager data | |
| CN113357666B (en) | Neural network-based furnace temperature measuring method, device and equipment | |
| CN110646099A (en) | Method and device for inverting target infrared radiation image based on measured data | |
| Arai et al. | Method for uncertainty evaluation of vicarious calibration of spaceborne visible to near infrared radiometers | |
| Gonzalez-Chavez et al. | Radiometric calibration of digital counts of infrared thermal cameras | |
| CN115112244A (en) | Satellite-borne thermal infrared camera absolute radiometric calibration processing method | |
| JP2018200295A (en) | Thermal imaging system | |
| Wang et al. | Inter-comparison of NPP/CrIS radiances with VIIRS, AIRS, and IASI: a post-launch calibration assessment | |
| Jiang et al. | Intercalibration of SVISSR/FY-2C infrared channels against MODIS/Terra and AIRS/Aqua channels | |
| CN111595458B (en) | Thermal infrared imager radiometric calibration method and device | |
| CN113029977B (en) | Automatic cross radiometric calibration method for wide-field-angle multispectral sensor | |
| CN115222619A (en) | Wall surface emissivity solving method based on thermal infrared image correction brightness of thermal imager | |
| Guoqing et al. | Nonlinear atmospheric correction based on neural network for infrared target radiometry | |
| CN114894737B (en) | Spectral reflectivity reconstruction method based on infrared image | |
| Chen et al. | Monitoring and assessment on radiometric stability of HJ-1A CCD using MODIS data | |
| CN117723157B (en) | Point target infrared radiation measuring method | |
| CN114894737A (en) | Spectral reflectivity reconstruction method based on infrared image |
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 |