CN107389198B - It is a kind of that window Surface Temperature Retrieval method is split based on radiance - Google Patents
It is a kind of that window Surface Temperature Retrieval method is split based on radiance Download PDFInfo
- Publication number
- CN107389198B CN107389198B CN201710372227.2A CN201710372227A CN107389198B CN 107389198 B CN107389198 B CN 107389198B CN 201710372227 A CN201710372227 A CN 201710372227A CN 107389198 B CN107389198 B CN 107389198B
- Authority
- CN
- China
- Prior art keywords
- radiance
- surface temperature
- formula
- radiation
- follows
- 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 22
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 19
- 230000005855 radiation Effects 0.000 claims abstract description 15
- 238000001514 detection method Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 abstract 1
- 229910052799 carbon Inorganic materials 0.000 abstract 1
- 230000003993 interaction Effects 0.000 abstract 1
- 238000006243 chemical reaction Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010206 sensitivity analysis 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiation Pyrometers (AREA)
Abstract
Surface temperature is the Important Parameters of the physical process of the earth and atmospheric interaction, is widely used in the research such as the hydrology, meteorology, Global carbon balance and whole world change.The present invention passes through the relationship directly established between two Detection Using Thermal Infrared Channel radiances, the solution for being used for radiation transfer equation group instead of simplifying Planck function, introduce an atmospheric parameter (ratio of Downward atmospheric long-wave radiation and atmosphere uplink radiation) simultaneously to improve the precision of simplified radiation transfer equation, the Split window algorithms (RBSWA) based on radiance are proposed, and algorithm is applied to MODIS data Surface Temperature Retrieval.With previously based on the Split window algorithms of brightness temperature compared with, the Surface Temperature Retrieval precision of RBSWA improves about 0.5K.Method And Principle proposed by the present invention is also applied for other IRMSS thermal bands.
Description
Technical field
The present invention relates to a kind of MODIS Surface Temperature Retrieval methods, can apply in agricultural, forestry, meteorology, ecological environment
And the remote sensing such as defense military department.
Background technique
Split window algorithms, also known as binary channels algorithm are thermal infrared remote sensing Surface Temperature Retrieval methods most general at present.Split window
Algorithm corrects atmospheric effect using two Detection Using Thermal Infrared Channel Atmospheric Absorption differences, calculates earth's surface true temperature.Current splits
Window algorithm is all the model established between surface temperature and two Detection Using Thermal Infrared Channel brightness temperatures, be may be collectively referred to as based on brightness temperature
The Split window algorithms of degree.For the model established between surface temperature and two Detection Using Thermal Infrared Channel brightness temperatures, need to simplify Pu Lang
Gram equation, the relationship between planck radiation brightness and brightness temperature is indicated with simple model.Most common simplified Pu Lang
The method of gram equation is Taylor Expansion.Under the conditions of earth's surface specific heat capacity is larger, Water Vapor Content is higher, surface temperature and star
Significant difference between upper brightness temperature, conversion accuracy between radiance and brightness temperature based on Taylor Expansion method compared with
Difference, transformed error are greater than 0.2K.The present invention directly establishes the relationship between two Detection Using Thermal Infrared Channel radiances, general instead of establishing
The bright gram of relationship between radiance and brightness temperature improves the conversion accuracy between two channel radiances, establishes
Relationship between surface temperature and two Detection Using Thermal Infrared Channel radiances proposes the Split window algorithms (letter based on radiance
Claim RBSWA), and it is applied to MODIS data Surface Temperature Retrieval.
Summary of the invention
The invention proposes one to split window Surface Temperature Retrieval method based on radiance, and MODIS data may be implemented
Surface Temperature Retrieval, the method is as follows:
Wherein B31(TS)=A0L31-A1L31+A2(formula 1)
In formula, TSIndicate surface temperature, B31(TS) expression temperature be TSCorresponding 31 wave band of MODIS of black matrix radiation it is bright
Degree, K31,1And K31,2Respectively indicate the thermal infrared radiation constant (can obtain from metadata) of 31 wave band of MODIS, L31And L32Point
It Biao Shi not radiance, A on the star of the wave band of MODIS 31 and 320、A1And A2Window coefficient is split in expression, and calculation formula is as follows:
In formula, k (L31) and b (L31) indicate to convert relevant parameter, C between two Detection Using Thermal Infrared Channel radiances31、C32、
D31And D32Indicate intermediate variable.k(L31) and b (L31) calculation formula is as follows:
In formula, K32,1And K32,2Respectively indicate the thermal infrared radiation constant (can obtain from metadata) of 32 wave band of MODIS.
C31、C32、D31And D32Calculation formula is as follows:
C31=ε31τ31
C32=ε32τ32
In formula, ε31And ε32Respectively indicate the Land surface emissivity of the wave band of MODIS 31 and 32, τ31And τ32It respectively indicates
The atmosphere effective transmissivity of the wave band of MODIS 31 and 32,WithRespectively indicate the atmosphere uplink spoke of the wave band of MODIS 31 and 32
Penetrate the ratio with Downward atmospheric long-wave radiation.τ31、τ32、WithCalculation formula it is as follows:
τ31(θ)=1- (0.088373w+0.04424)/cos (θ)
τ32(θ)=1- (0.096298w+0.118442)/cos (θ)
W indicates Water Vapor Content in formula.Split window algorithms proposed by the present invention based on radiance need 5 input ginsengs
Amount, including L31、L32、ε31、ε32And w.
Specific embodiment
1, radiance on star, radiance L on star are calculated31And L32It can be calibrated and be obtained via radiation by image DN value；
2, Land surface emissivity, Land surface emissivity ε are calculated31And ε32Can with NDVI threshold method [Sobrino, J.A.,
Jiménez-J.C.,Sòria,G.,Romaguera,M.,Guanter,L.,Moreno,J.,Plaza,A.,&Martí
nez,P.(2008).Land surface emissivity retrieval from different VNIR and TIR
Sensors.IEEE Transactions on Geoscience and Remote Sensing, 46,316-327], Huo Zheqi
Its method [Li, Z.-L., Wu, H., Wang, N., Qiu, S., Sobrino, J.A., Wan, Z., Tang, B.-H., &Yan, G.
(2013b).Land surface emissivity retrieval from satellite data.International
Journal of Remote Sensing, 34,3084-3127] it is calculated；
3, Water Vapor Content is calculated, Water Vapor Content w can use ratio method [A new approach for
retrieving precipitable water from ATSR2split-window channel data over land
Area.International Journal of Remote Sensing, 24,5095-5117] it is obtained from near-infrared data inversion
Arrive, can also with based on NDVI split window covariance and Variance ratio method (SWCVR) [Wang, M., He, G., Zhang, Z., Wang,
G.,&Long,T.(2015).NDVI-based split-window algorithm for precipitable water
vapour retrieval from Landsat-8TIRS data over land area.Remote Sensing
Letters, 6,904-913] it is obtained from Thermal Infrared Data inverting；
4, surface temperature is calculated, calculates surface temperature (such as using the Split window algorithms proposed by the present invention based on radiance
Formula 1).
5, arithmetic accuracy is assessed.It is calculated based on the model accuracy of simulated data sets evaluation RBSWA, and with two classical windows of splitting
Method compares (Split window algorithms of the propositions such as Qin, abbreviation SWAQ[Qin,Z.,Dall'Olmo,G.,Karnieli,A.,&
Berliner,P.(2001b).Derivation of split window algorithm and its sensitivity
analysis for retrieving land surface temperature from NOAA‐advanced very high
resolution radiometer data.Journal of Geophysical Research:Atmospheres(1984–
2012),106,22655-22670]；The Split window algorithms that Sobrino et al. is proposed, abbreviation SWAS[Sobrino,J.,&
Raissouni,N.(2000).Toward remote sensing methods for land cover dynamic
monitoring:application to Morocco.International Journal of Remote Sensing,21,
353-366]).The ε that analogue data is concentrated31,ε32,L31,L32, then the input of θ and w as Split window algorithms compares Split window algorithms
The LST that the LST and analogue data that inverting obtains are concentrated.Arithmetic accuracy assessment result such as table 1 shows to calculate with the window of splitting that two classical
Method compares, and RBSWA Surface Temperature Retrieval precision improves about 0.5K.
The different Split window algorithms Surface Temperature Retrieval precision of table 1 compare
Claims (1)
1. a kind of split window Surface Temperature Retrieval algorithm based on radiance, its main feature is that with radiance on star directly as calculation
The input of method；Step are as follows:
The first step calculates radiance on star, radiance L on star31And L32It is calibrated and is obtained via radiation by image DN value；
Second step calculates Land surface emissivity, Land surface emissivity ε31And ε32It is calculated with NDVI threshold method；
Third step calculates Water Vapor Content, and Water Vapor Content w is obtained with ratio method from near-infrared data inversion, or uses base
It is obtained in split window covariance and the Variance ratio method of NDVI from Thermal Infrared Data inverting；
4th step calculates surface temperature, the calculation formula of surface temperature are as follows:
Wherein B31(TS)=A0L31-A1L32+A2；
In formula, TSIndicate surface temperature, B31(TS) expression temperature be TSThe corresponding MODIS31 wave band of black matrix radiance,
K31,1And K31,2Respectively indicate the thermal infrared radiation constant of MODIS31 wave band, L31And L32Respectively indicate MODIS31 and 32 wave bands
Radiance on star, A0、A1And A2Window coefficient is split in expression, and calculation formula is as follows:
In formula, k (L31) and b (L31) indicate to convert relevant parameter, C between two Detection Using Thermal Infrared Channel radiances31、C32、D31With
D32Indicate intermediate variable, k (L31) and b (L31) calculation formula it is as follows:
In formula, K32,1And K32,2Respectively indicate the thermal infrared radiation constant of MODIS32 wave band；
C31、C32、D31And D32Calculation formula it is as follows:
C31=ε31τ31；
C32=ε32τ32；
In formula, ε31And ε32Respectively indicate the Land surface emissivity of MODIS31 and 32 wave bands, τ31And τ32Respectively indicate MODIS31 and
The atmosphere effective transmissivity of 32 wave bands,WithRespectively indicate the atmosphere uplink radiation and downward atmospheric of MODIS31 and 32 wave bands
The ratio of radiation；τ31、τ32、WithCalculation formula it is as follows:
τ31(θ)=1- (0.088373w+0.04424)/cos (θ)
τ32(θ)=1- (0.096298w+0.118442)/cos (θ)
In formula, w indicates Water Vapor Content.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710372227.2A CN107389198B (en) | 2017-05-23 | 2017-05-23 | It is a kind of that window Surface Temperature Retrieval method is split based on radiance |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710372227.2A CN107389198B (en) | 2017-05-23 | 2017-05-23 | It is a kind of that window Surface Temperature Retrieval method is split based on radiance |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107389198A CN107389198A (en) | 2017-11-24 |
| CN107389198B true CN107389198B (en) | 2019-09-20 |
Family
ID=60338425
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710372227.2A Expired - Fee Related CN107389198B (en) | 2017-05-23 | 2017-05-23 | It is a kind of that window Surface Temperature Retrieval method is split based on radiance |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107389198B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109297605B (en) * | 2018-10-09 | 2020-02-11 | 北京大学 | Surface temperature inversion method based on mid-infrared and thermal infrared data |
| CN109596223B (en) * | 2018-12-24 | 2019-10-22 | 中国农业科学院农业资源与农业区划研究所 | A kind of multispectral satellite borne sensor lane device and method of direct inversion surface temperature |
| CN109632106B (en) * | 2019-01-03 | 2020-06-30 | 北京师范大学 | Remote sensing surface temperature product angle effect correction method |
| CN109887615B (en) * | 2019-01-30 | 2020-12-11 | 北京环境特性研究所 | Earth surface temperature periodic daily change simulation method |
| CN114894321B (en) * | 2022-07-13 | 2022-10-21 | 国家卫星气象中心(国家空间天气监测预警中心) | Calibration method of infrared remote sensing instrument, electronic device and computer storage medium |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101629850A (en) * | 2009-08-24 | 2010-01-20 | 中国农业科学院农业资源与农业区划研究所 | Method for inversing land surface temperature from MODIS data |
| US9551619B1 (en) * | 2011-09-23 | 2017-01-24 | Rockwell Collins, Inc. | Terahertz laser |
| CN104748857A (en) * | 2015-03-05 | 2015-07-01 | 北京师范大学 | Method and system for inverting urban surface temperatures |
| CN104897289B (en) * | 2015-06-23 | 2017-08-18 | 中国科学院遥感与数字地球研究所 | A kind of satellite data Surface Temperature Retrieval methods of Landsat 8 |
-
2017
- 2017-05-23 CN CN201710372227.2A patent/CN107389198B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN107389198A (en) | 2017-11-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107389198B (en) | It is a kind of that window Surface Temperature Retrieval method is split based on radiance | |
| Sobrino et al. | Land surface temperature retrieval from MSG1-SEVIRI data | |
| Liang | An optimization algorithm for separating land surface temperature and emissivity from multispectral thermal infrared imagery | |
| Sobrino et al. | A comparative study of land surface emissivity retrieval from NOAA data | |
| Atitar et al. | A split-window algorithm for estimating LST from Meteosat 9 data: Test and comparison with in situ data and MODIS LSTs | |
| Hall et al. | Comparison of satellite-derived and in-situ observations of ice and snow surface temperatures over Greenland | |
| Liu et al. | Spectral calibration of hyperspectral data observed from a hyperspectrometer loaded on an unmanned aerial vehicle platform | |
| CN101629850A (en) | Method for inversing land surface temperature from MODIS data | |
| CN103744069A (en) | Methane profile orthogonal inversion method aiming at AIRS (atmospheric infrared sounder) hyper-spectrum satellite data | |
| CN113776671B (en) | Day and night remote sensing image short wave infrared volcano temperature-emissivity joint inversion method | |
| Kwon et al. | Use of total precipitable water classification of a priori error and quality control in atmospheric temperature and water vapor sounding retrieval | |
| Danodia et al. | Land surface temperature retrieval by radiative transfer equation and single channel algorithms using landsat-8 satellite data | |
| Qi et al. | Calibration and validation of the InfraRed Atmospheric Sounder onboard the FY3B satellite | |
| Mira et al. | Quantifying uncertainties in land surface temperature due to atmospheric correction: Application to Landsat-7 data over a Mediterranean agricultural region | |
| Zhang et al. | An improved 6s code for atmospheric correction based on water vapor content | |
| Xia et al. | Retrieval land surface temperature from visible infrared imager radiometer suite data | |
| Xie et al. | Progress of remote sensing inversion research on surface temperature of lakes and reservoirs | |
| Guo et al. | Near-surface air temperature retrieval from Chinese Geostationary FengYun Meteorological Satellite (FY-2C) data | |
| Jiménez‐Muñoz et al. | Atmospheric water vapour content retrieval from visible and thermal data in the framework of the DAISEX campaigns | |
| Peng et al. | Estimation of evaporative fraction from top-of-atmosphere radiance | |
| Gao et al. | Study of aerosol influence on nighttime land surface temperature retrieval based on two methods | |
| Wu et al. | A Fast Parametric Model of Estimating Atmospheric Parameters for Landsat 8 Thermal Infrared Sensor | |
| Shi et al. | Atmospheric correction of passive microwave brightness temperature on the estimation of snow depth | |
| Al-Baldawi et al. | The Effect of Changes in Land Cover on the Surface Temperature in Hor Al-Al-Dalmaj Technologies Using Remote Sensing and Geographic Information Systems | |
| Chen et al. | Derivation of new split window algorithm for retrieving land surface temperature from FY-3/VIRR data |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190920 Termination date: 20210523 |