CN105160631B - A kind of method for seeking radiant correction coefficient - Google Patents
A kind of method for seeking radiant correction coefficient Download PDFInfo
- Publication number
- CN105160631B CN105160631B CN201510381649.7A CN201510381649A CN105160631B CN 105160631 B CN105160631 B CN 105160631B CN 201510381649 A CN201510381649 A CN 201510381649A CN 105160631 B CN105160631 B CN 105160631B
- Authority
- CN
- China
- Prior art keywords
- pixel
- wave band
- values
- radiant correction
- ref
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000009021 linear effect Effects 0.000 claims abstract description 27
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 claims abstract description 25
- 238000002310 reflectometry Methods 0.000 claims description 38
- 230000005855 radiation Effects 0.000 claims description 14
- 239000004568 cement Substances 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000003331 infrared imaging Methods 0.000 claims description 3
- 238000012417 linear regression Methods 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 241001269238 Data Species 0.000 claims 1
- 238000009738 saturating Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 4
- 230000007812 deficiency Effects 0.000 abstract description 2
- 238000001228 spectrum Methods 0.000 abstract description 2
- 230000003595 spectral effect Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 210000001747 pupil Anatomy 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Landscapes
- Image Processing (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
A kind of method for seeking radiant correction coefficient belongs to infrared image processing technology field.First the pixel to match with ground atural object scaling point is found out on the original image, and record the position of pixel;Then radiant correction is carried out to original image, then atmospheric correction is carried out to the image after radiant correction with laboratory radiant correction coefficient;By ground atural object scaling point and it is above-mentioned it is processed after image on matching pixel radiance value composition data pair, again by the data of different atural objects to carrying out linear fit, obtain the slope and intercept of the wave band, the slope of point value is corresponded in last and laboratory radiant correction coefficient and intercept carries out that new slope and intercept, i.e., new radiant correction coefficient is calculated.The method can solve the deficiency of laboratory Absolute Radiometric Calibration Coefficients, with the image after the radiant correction coefficient processing and Absolute Radiometric Calibration Coefficients, treated that image is compared, the former can save the time asked used in radiant correction coefficient compared with can be close to the actual spectrum curve of atural object.
Description
Technical field
The present invention relates to a kind of methods for seeking radiant correction coefficient in infrared technique, belong to infrared image processing technology neck
Domain.
Background technology
The atural object radiation information that thermal infrared sensor is obtained is other than the interference by atmospheric effect, and also there is a systems
Row systematic error, such as the variation of recording noise, reference temperature and detector error.It is accurate in order to be obtained from sensing data
Radiation information, sensor must establish quantitative relationship by radiant correction between output valve and the radiance value of incidence.
Therefore it before image data is acquired, needs in laboratory to the input value of sensor and integrating sphere spectral radiance value into rower
It is fixed, linear relationship coefficient between the two, as absolute radiation correction coefficient is obtained.
When dealing with remote sensing images, original image is carried out to need to provide accurate radiant correction during radiant correction
Coefficient, the data obtained in this way just can be closer to the spectral charactersiticss of targets.Ren Jianwei, Wan Zhi, Li Xiansheng, appoint the paper for building Yue Fabiao
《The radiation transmission characteristic of space optical remote sensor and bearing calibration》(optical precision engineering, 2007,15 (2):1188) it is provided in
Radiant correction coefficient seeks method:Represent that the target Equivalent under certain observation condition is that entrance pupil spoke brightness is expressed as figure with Taylor series
As the function of gray scale:
L=A'0+A′1*DN+A'2*DN2+A′3*DN3+...+A'n*DNn
Wherein A'0,A′1,A'2...A'nIt is radiant correction coefficient, once the calibration coefficient is determined, so that it may with according to this
A functional relation and the digital output DN values of camera are finally inversed by the radiance L at camera entrance pupil, finally restore the spoke of target
Penetrate luminance picture.This radiant correction coefficient asks method comparatively laborious, and without being combined with practical spectral charactersitics of targets curve,
Therefore its accuracy is poor.
Invention content
In order to overcome defect and deficiency of the existing technology, the present invention proposes a kind of side for seeking radiant correction coefficient
Method, this method are based on Absolute Radiometric Calibration Coefficients, with reference to practical atural object radiance value, acquire improved radiant correction coefficient
Method.
Technical scheme is as follows:
A kind of method for seeking radiant correction coefficient carries out integrating sphere calibration data data reading, analysis by computer
And calculation processing, based on Absolute Radiometric Calibration Coefficients, with reference to practical atural object radiance value, acquire improved radiant correction system
Number, this method step are as follows:
(1) the DN values of pixel and integrating sphere radiance value, product are read in by wave band to integrating sphere calibration data by computer
Each pixel DN values of bulb separation calibration data are denoted as DN0(i, j), integrating sphere radiance value are denoted as L0(j), wherein i is pixel
Number, j be wave band number, DN0Represent detector pixel Digital output numerical value, L0For integrating sphere spoke brightness value, due to infrared imaging system
System is linear, so the Digital output number of the pixel of detector is obtained using the linear absolute radiometric calibration formula in laboratory
Value:
DN0(i, j)=k0(i,j)*L0(j),
Wherein k0(i, j) is denoted as in the radiant correction coefficient that wave band is j, pixel is i;
(2) raw image data and radiant correction coefficient data are read by computer, obtains each of raw image data
The pixel DN values DN of wave band1(i, j), by DN1(i, j) substitutes into radiant correction formula, acquires every after original image corrects via radiation
The radiance value of a pixel, is denoted as L1(i, j), radiant correction formula are:
(3) the various radiation energies that remote sensing is utilized are intended to that interaction-or scattering occurs with earth atmosphere or absorb,
And make energy attenuation, and spatial distribution is made to change, the attenuation of air be to the light of different wave length it is selective, because
And air is different to the image of different-waveband, therefore the processing of these atmospheric effects is eliminated to the image of radiant correction,
Reflectivity for Growing Season is obtained in selecting experience linear approach, with ASD spectrometers respectively geodetic object be water body and cement dyke DN values and ground
Table reflects Ref, and water body is shown as dark target here, and cement dyke is shown as bright target, and the DN values and reflectivity Ref measured are by wave band
Storage, the DN values and reflectivity Ref of dark target are DNl(j) and Refl(j), the DN values of bright target and reflectivity Ref are DNh(j)
And Refh(j), the DN values for making each wave band by linear regression technique are converted to reflectivity, and dark target and bright is read in by computer
The DN values and reflectivity Ref of dark target and bright target substitution DN values are converted to reflectivity by the DN values and reflectivity of target respectively
In the linear representation of Ref, the Product-factor A (j) of each wave band and plus-minus item B (j) is calculated, by Product-factor A (j) and is added
Deduction item B (j) is stored by wave band, and here, the linear representation that DN values are converted to reflectivity Ref is:
Ref(iASD, j) and=A (j) * DN (iASD, j) and+B (j), wherein DN (iASD, j) and represent j-th of wave band, i-th of atural object
The ASD spectrometer Digital output numerical value of pixel, Ref (iASD, j) and represent corresponding j-th of wave band, i-th of atural object pixel reflectivity
Value, A (j) is Product-factor, related with atmospheric transmittance and instrument gain, and B (j) is plus-minus item, with dark current and air journey spoke
It penetrates related;
(4) the DN values DN of raw image data is read in by computer1In (i, j) and step (3) product of atmospheric correction because
Sub- A (j) and plus-minus item B (j), by wave band by DN1(i, j), A (j) and B (j) substitute into the linear list that DN values are converted to reflectivity Ref
Up in formula, the reflectivity Ref (i, j) of each pixel of raw image data is obtained;
(5) the radiance value L of image data after radiant correction in step (2) is read in by computer1(i, j) and step
(4) the reflectivity Ref (i, j) of each pixel of raw image data in, by wave band by L1(i, j) and the mutually multiplied L of Ref (i, j)2(i,
J), L2(i, j) is the radiance value that pixel after atmospheric correction is carried out to image data after radiant correction in step (2);
(6) radiance value of ground atural object scaling point is measured by ASD spectrometers, if its radiance value measured isHere n represents different atural object;
(7) the atural object calibration dot position information file of ASD records is opened, is opened simultaneously through step (5) treated image,
The pixel consistent with geographical location information is found out on the image, and records the radiance value of the pixel
(10) k is read by wave band2(j), k1(j) and d1(j), note K (j)=k1(j)*k2(j), D (j)=d1(j), it is obtained
The K (j) and D (j) of each wave band are new radiant correction coefficient by wave band storage K (j) and D (j), the K (j) acquired and D (j),
The radiation that j band image DN value DN (j) can be converted to j band images by conversion formula L (j)=K (j) * DN (j)+D (j) is bright
Angle value L (j).
The DN values refer to Digital output numerical value.
Radiant correction is carried out to original image the beneficial effects of the invention are as follows obtained radiant correction coefficient, after radiant correction
Image close to atural object actual spectrum curve.
Specific embodiment
With reference to embodiment, the invention will be further described, but not limited to this.
Embodiment:
The embodiment of the present invention is as follows, a kind of method for seeking radiant correction coefficient, by computer to integrating sphere calibration data
Data reading, analysis and calculation processing are carried out, based on Absolute Radiometric Calibration Coefficients, with reference to practical atural object radiance value, is acquired
Improved radiant correction coefficient, this method step are as follows:
(1) the DN values of pixel and integrating sphere radiance value, product are read in by wave band to integrating sphere calibration data by computer
Each pixel DN values of bulb separation calibration data are denoted as DN0(i, j), integrating sphere radiance value are denoted as L0(j), wherein i is pixel
Number, j be wave band number, DN0Represent detector pixel Digital output numerical value, L0For integrating sphere spoke brightness value, due to infrared imaging system
System is linear, so the Digital output number of the pixel of detector is obtained using the linear absolute radiometric calibration formula in laboratory
Value:
DN0(i, j)=k0(i,j)*L0(j),
Wherein k0(i, j) is denoted as in the radiant correction coefficient that wave band is j, pixel is i;
(2) raw image data and radiant correction coefficient data are read by computer, obtains each of raw image data
The pixel DN values DN of wave band1(i, j), by DN1(i, j) substitutes into radiant correction formula, acquires every after original image corrects via radiation
The radiance value of a pixel, is denoted as L1(i, j), radiant correction formula are:
(3) the various radiation energies that remote sensing is utilized are intended to that interaction-or scattering occurs with earth atmosphere or absorb,
And make energy attenuation, and spatial distribution is made to change, the attenuation of air be to the light of different wave length it is selective, because
And air is different to the image of different-waveband, therefore the processing of these atmospheric effects is eliminated to the image of radiant correction,
Reflectivity for Growing Season is obtained in selecting experience linear approach, with ASD spectrometers respectively geodetic object be water body and cement dyke DN values and ground
Table reflects Ref, and water body is shown as dark target here, and cement dyke is shown as bright target, and the DN values and reflectivity Ref measured are by wave band
Storage, the DN values and reflectivity Ref of dark target are DNl(j) and Refl(j), the DN values of bright target and reflectivity Ref are DNh(j)
And Refh(j), the DN values for making each wave band by linear regression technique are converted to reflectivity, and dark target and bright is read in by computer
The DN values and reflectivity Ref of dark target and bright target substitution DN values are converted to reflectivity by the DN values and reflectivity of target respectively
In the linear representation of Ref, the Product-factor A (j) of each wave band and plus-minus item B (j) is calculated, by Product-factor A (j) and is added
Deduction item B (j) is stored by wave band, and here, the linear representation that DN values are converted to reflectivity Ref is:
Ref(iASD, j) and=A (j) * DN (iASD, j) and+B (j), wherein DN (iASD, j) and represent j-th of wave band, i-th of atural object
The ASD spectrometer Digital output numerical value of pixel, Ref (iASD, j) and represent corresponding j-th of wave band, i-th of atural object pixel reflectivity
Value, A (j) is Product-factor, related with atmospheric transmittance and instrument gain, and B (j) is plus-minus item, with dark current and air journey spoke
It penetrates related;
(4) the DN values DN of raw image data is read in by computer1In (i, j) and step (3) product of atmospheric correction because
Sub- A (j) and plus-minus item B (j), by wave band by DN1(i, j), A (j) and B (j) substitute into the linear list that DN values are converted to reflectivity Ref
Up in formula, the reflectivity Ref (i, j) of each pixel of raw image data is obtained;
(5) the radiance value L of image data after radiant correction in step (2) is read in by computer1(i, j) and step
(4) the reflectivity Ref (i, j) of each pixel of raw image data in, by wave band by L1(i, j) and the mutually multiplied L of Ref (i, j)2(i,
J), L2(i, j) is the radiance value that pixel after atmospheric correction is carried out to image data after radiant correction in step (2);
(6) radiance value of ground atural object scaling point is measured by ASD spectrometers, if its radiance value measured isHere n represents different atural object;
(7) the atural object calibration dot position information file of ASD records is opened, is opened simultaneously through step (5) treated image,
The pixel consistent with geographical location information is found out on the image, and records the radiance value of the pixel
(10) k is read by wave band2(j), k1(j) and d1(j), note K (j)=k1(j)*k2(j), D (j)=d1(j), it is obtained
The K (j) and D (j) of each wave band are new radiant correction coefficient by wave band storage K (j) and D (j), the K (j) acquired and D (j),
The radiation that j band image DN value DN (j) can be converted to j band images by conversion formula L (j)=K (j) * DN (j)+D (j) is bright
Angle value L (j).
Claims (1)
1. a kind of method for seeking radiant correction coefficient, integrating sphere calibration data are carried out by computer data reading, analysis and
Calculation processing based on Absolute Radiometric Calibration Coefficients, with reference to practical atural object radiance value, acquires improved radiant correction system
Number, this method step are as follows:
(1) the DN values of pixel and integrating sphere radiance value, integrating sphere are read in by wave band to integrating sphere calibration data by computer
Each pixel DN values of calibration data are denoted as DN0(i, j), integrating sphere radiance value are denoted as L0(j), wherein i is pixel number, and j is
Wave band number, DN0Represent detector pixel Digital output numerical value, L0For integrating sphere spoke brightness value, since infrared imaging system is line
Property, so the Digital output numerical value of the pixel of detector is obtained using the linear absolute radiometric calibration formula in laboratory:
DN0(i, j)=k0(i,j)*L0(j),
Wherein k0(i, j) is denoted as in the radiant correction coefficient that wave band is j, pixel is i;
(2) raw image data and radiant correction coefficient data are read by computer, obtains each wave band of raw image data
Pixel DN values DN1(i, j), by DN1(i, j) substitutes into radiant correction formula, acquires each picture after original image corrects via radiation
The radiance value of member, is denoted as L1(i, j), radiant correction formula are:
(3) the various radiation energies that remote sensing is utilized are intended to that interaction-or scattering occurs with earth atmosphere or absorb, and make
Energy attenuation, and spatial distribution is made to change, the attenuation of air is selective to the light of different wave length, thus big
Gas is different to the image of different-waveband, therefore the processing of these atmospheric effects is eliminated to the image of radiant correction, is selected
Reflectivity for Growing Season is obtained in experience linear approach, and with ASD spectrometers, geodetic object is that the DN values and earth's surface of water body and cement dyke are anti-respectively
Ref is penetrated, water body is shown as dark target here, and cement dyke is shown as bright target, and the DN values and reflectivity Ref measured are deposited by wave band
Storage, the DN values and reflectivity Ref of dark target are DNl(j) and Refl(j), the DN values of bright target and reflectivity Ref are DNh(j) and
Refh(j), the DN values for making each wave band by linear regression technique are converted to reflectivity, and dark target and bright mesh are read in by computer
The DN values and reflectivity Ref of dark target and bright target substitution DN values are converted to reflectivity Ref by target DN values and reflectivity respectively
Linear representation in, calculate each wave band Product-factor A (j) and plus-minus item B (j), by Product-factor A (j) and plus-minus
Item B (j) is stored by wave band, and here, the linear representation that DN values are converted to reflectivity Ref is:Ref(iASD, j) and=A (j) * DN
(iASD, j) and+B (j), wherein DN (iASD, j) and represent the ASD spectrometer Digital output numbers of j-th of wave band, i-th atural object pixel
Value, Ref (iASD, j) and represent corresponding j-th of wave band, i-th of atural object pixel reflectance value, A (j) is Product-factor, saturating with air
It is related with instrument gain to cross rate, B (j) is plus-minus item, related with dark current and atmospheric path radiation;
(4) the DN values DN of raw image data is read in by computer1The Product-factor A (j) of atmospheric correction in (i, j) and step (3)
With plus-minus item B (j), by wave band by DN1(i, j), A (j) and B (j) substitute into the linear representation that DN values are converted to reflectivity Ref
In, the reflectivity Ref (i, j) of each pixel of raw image data is obtained;
(5) the radiance value L of image data after radiant correction in step (2) is read in by computer1(i, j) and step (4) Central Plains
The reflectivity Ref (i, j) of each pixel of beginning image data, by wave band by L1(i, j) and the mutually multiplied L of Ref (i, j)2(i, j), L2
(i, j) is the radiance value that pixel after atmospheric correction is carried out to image data after radiant correction in step (2);
(6) radiance value of ground atural object scaling point is measured by ASD spectrometers, if its radiance value measured is
Here n represents different atural object;
(7) the atural object calibration dot position information file of ASD records is opened, opens simultaneously through step (5) treated image, is scheming
The pixel consistent with geographical location information is found out on picture, and records the radiance value of the pixel
(8) identical wave band j, the radiance value of ground culture pointWith the radiance value of the pixel in correspondence imageForm a pair of of discrete point, the radiance value of Different Ground culture pointWith the radiation of the pixel in correspondence image
Brightness valueOne group of discrete point pair is formed, all discrete points datas are read by computer, it is discrete to every group by wave band j
Point pairLeast-squares algorithm linear fitting is carried out, obtains full wave linear fit slope k1(j) and line
Property fitting intercept d1(j), by linear fit slope k1(j) and linear fit intercept d1(j) wave band number storage is pressed, wherein, Linear Quasi
Closing expression formula is:
(9) the pixel i found on step (7) image is recorded in, by DN of all pixels number for i in integrating sphere calibration data
Value DN0(i, j) is preserved by wave band to txt file, then the DN of all pixels of each wave band to finding0(i, j) averages
It arrivesIt willWith integrating sphere radiance value L0(j) it carries out the fitting of two point Linears and acquires k2(j), k2(j) herein
For the radiant correction coefficient after the mean value of all pixel DN values and integrating sphere radiance value linear fit, wherein linear fit table
It is up to formula:
(10) k is read by wave band2(j), k1(j) and d1(j), note K (j)=k1(j)*k2(j), D (j)=d1(j), it is obtained each
The K (j) and D (j) of wave band are new radiant correction coefficient by wave band storage K (j) and D (j), the K (j) acquired and D (j), by turning
J band image DN value DN (j) can be converted to the radiance value L of j band images by changing formula L (j)=K (j) * DN (j)+D (j)
(j)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510381649.7A CN105160631B (en) | 2015-07-02 | 2015-07-02 | A kind of method for seeking radiant correction coefficient |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510381649.7A CN105160631B (en) | 2015-07-02 | 2015-07-02 | A kind of method for seeking radiant correction coefficient |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105160631A CN105160631A (en) | 2015-12-16 |
CN105160631B true CN105160631B (en) | 2018-06-29 |
Family
ID=54801474
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510381649.7A Expired - Fee Related CN105160631B (en) | 2015-07-02 | 2015-07-02 | A kind of method for seeking radiant correction coefficient |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105160631B (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106485683B (en) * | 2016-10-20 | 2019-04-16 | 中国科学院上海技术物理研究所启东光电遥感中心 | A kind of image adaptive non-uniform correction method based on scene |
CN106600646B (en) * | 2016-11-25 | 2019-05-24 | 北京空间机电研究所 | A kind of in-orbit image conformity bearing calibration of infrared sweeping camera |
CN107656289A (en) * | 2017-08-23 | 2018-02-02 | 中国科学院光电研究院 | Spaceborne optics load absolute radiation calibration method and system based on ground spoke brightness |
CN107860478B (en) * | 2017-11-17 | 2019-11-22 | 北京长峰科威光电技术有限公司 | A kind of restorative procedure of infrared imaging system cold emission blackspot |
CN108090883B (en) * | 2018-01-04 | 2020-05-05 | 中煤航测遥感集团有限公司 | Hyperspectral image preprocessing method and device and electronic equipment |
CN108174127A (en) * | 2018-01-30 | 2018-06-15 | 中国科学院长春光学精密机械与物理研究所 | Relative radiometric correction methods of the face battle array CMOS under global shutter working method |
CN109035150B (en) * | 2018-04-24 | 2019-05-10 | 湖南省有色地质勘查研究院 | A kind of remote sensing image data defogging method based on middle-low alloy steels |
CN114528709B (en) * | 2022-02-24 | 2024-08-13 | 山东省水利综合事业服务中心 | Remote sensing reflectivity correction method based on illumination intensity |
CN114965306B (en) * | 2022-05-27 | 2024-02-20 | 淮阴师范学院 | Calibration method of optical sensor for measuring reflectivity |
CN115615938B (en) * | 2022-12-14 | 2023-03-28 | 天津中科谱光信息技术有限公司 | Water quality analysis method and device based on reflection spectrum and electronic equipment |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101442608A (en) * | 2008-12-31 | 2009-05-27 | 中国资源卫星应用中心 | Method for improving relative radiation correction of CCD camera |
CN102324098A (en) * | 2011-08-23 | 2012-01-18 | 中国资源卫星应用中心 | Relative radiation calibration method in combination with laboratory calibration and even landscape statistics |
CN104580944A (en) * | 2013-10-10 | 2015-04-29 | 中国科学院光电研究院 | Method for carrying out relative radiometric calibration on CCD (Charge-Coupled Device) image |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7912180B2 (en) * | 2009-02-19 | 2011-03-22 | Kabushiki Kaisha Toshiba | Scattered radiation correction method and scattered radiation correction apparatus |
-
2015
- 2015-07-02 CN CN201510381649.7A patent/CN105160631B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101442608A (en) * | 2008-12-31 | 2009-05-27 | 中国资源卫星应用中心 | Method for improving relative radiation correction of CCD camera |
CN102324098A (en) * | 2011-08-23 | 2012-01-18 | 中国资源卫星应用中心 | Relative radiation calibration method in combination with laboratory calibration and even landscape statistics |
CN104580944A (en) * | 2013-10-10 | 2015-04-29 | 中国科学院光电研究院 | Method for carrying out relative radiometric calibration on CCD (Charge-Coupled Device) image |
Non-Patent Citations (4)
Title |
---|
Vicarious radiometric calibration of a multispectral sensor from an aerial trike applied to precision agriculture;Herrero-Huerta M et al.;《Computers & Electronics in Agriculture》;20141231;全文 * |
基于实验室定标和均匀景统计的相对辐射定标方法;曾湧 等;《航天返回与遥感》;20120831;第33卷(第4期);全文 * |
机载推帚式高光谱图像预处理技术研究与实现;吕桂龙;《中国优秀硕士学位论文全文数据库 信息科技辑》;20141015(第10期);全文 * |
高光谱成像仪的高精度实验室辐射定标及校正方法;马庆军;《硅谷》;20121231(第19期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105160631A (en) | 2015-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105160631B (en) | A kind of method for seeking radiant correction coefficient | |
CN108662991B (en) | Land block scale winter wheat leaf area index estimation method based on remote sensing satellite data | |
Duan et al. | Inversion of the PROSAIL model to estimate leaf area index of maize, potato, and sunflower fields from unmanned aerial vehicle hyperspectral data | |
Fensholt | Earth observation of vegetation status in the Sahelian and Sudanian West Africa: comparison of Terra MODIS and NOAA AVHRR satellite data | |
Sun et al. | Aerosol optical depth retrieval by HJ-1/CCD supported by MODIS surface reflectance data | |
Du et al. | Detecting the effects of climate change on canopy phenology in coniferous forests in semi-arid mountain regions of China | |
Cheng et al. | Generating the 30-m land surface temperature product over continental China and USA from Landsat 5/7/8 data | |
CN102901516A (en) | Multispectral image radiation correction method based on absolute radiometric calibration | |
CN111323129A (en) | Earth surface temperature inversion method based on broadband thermal infrared image | |
CN116519557B (en) | Aerosol optical thickness inversion method | |
CN109632106B (en) | Remote sensing surface temperature product angle effect correction method | |
Wan et al. | MODIS land surface temperature and emissivity | |
Nazeer et al. | Combining landsat TM/ETM+ and HJ-1 A/B CCD sensors for monitoring coastal water quality in Hong Kong | |
KR20160073845A (en) | Device and method for calibration of high resolution electro optical satellite | |
Kusnierek et al. | Challenges in using an analog uncooled microbolometer thermal camera to measure crop temperature | |
Wen et al. | Parametrized BRDF for atmospheric and topographic correction and albedo estimation in Jiangxi rugged terrain, China | |
Sèze et al. | Evaluation against CALIPSO lidar observations of the multi‐geostationary cloud cover and type dataset assembled in the framework of the Megha‐Tropiques mission | |
Li et al. | A quality control procedure for FY-3A MWTS measurements with emphasis on cloud detection using VIRR cloud fraction | |
CN116977873A (en) | Exposed sand decision tree extraction method based on high-resolution satellite remote sensing | |
Park et al. | Development of sea surface temperature retrieval algorithms for geostationary satellite data (Himawari-8/AHI) | |
Kabsch et al. | Initial results of the land surface temperature (LST) validation with the Evora, Portugal ground‐truth station measurements | |
Liu et al. | Fractional coverage rather than green chromatic coordinate is a robust indicator to track grassland phenology using smartphone photography | |
Borel et al. | Data processing and temperature-emissivity separation for tower-based imaging Fourier transform spectrometer data | |
Schmidt et al. | A method for operational calibration of AVHRR reflective time series data | |
Gascon et al. | Using multi‐directional high‐resolution imagery from POLDER sensor to retrieve leaf area index |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180629 |
|
CF01 | Termination of patent right due to non-payment of annual fee |