CN107389198A - 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
- CN107389198A CN107389198A CN201710372227.2A CN201710372227A CN107389198A CN 107389198 A CN107389198 A CN 107389198A CN 201710372227 A CN201710372227 A CN 201710372227A CN 107389198 A CN107389198 A CN 107389198A
- Authority
- CN
- China
- Prior art keywords
- radiance
- surface temperature
- calculated
- retrieval
- window
- 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
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 21
- 230000005855 radiation Effects 0.000 claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000013459 approach Methods 0.000 claims description 2
- 238000000179 transient infrared spectroscopy Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 7
- 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
- 239000012141 concentrate Substances 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000011159 matrix material Substances 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 is by directly establishing the relation between two Detection Using Thermal Infrared Channel radiances, it is used for the solution of radiation transfer equation group instead of simplifying Planck function, an atmospheric parameter (ratio of Downward atmospheric long-wave radiation and air uplink radiation) is introduced 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 Retrievals.With previously based on the Split window algorithms of brightness temperature compared with, RBSWA Surface Temperature Retrieval precision 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 Retrievals method, can apply in agricultural, forestry, meteorology, ecological environment
And the remote sensing such as national defense and military department.
Background technology
Split window algorithms, also known as binary channels algorithm, it is thermal infrared remote sensing Surface Temperature Retrieval method 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, may be collectively referred to as being based on brightness temperature
The Split window algorithms of degree.In order to establish the model between surface temperature and two Detection Using Thermal Infrared Channel brightness temperatures, it is necessary to simplify Pu Lang
Gram equation, the relation between planck radiation brightness and brightness temperature is represented with simple model.The most frequently used simplification Pu Lang
The method of gram equation is Taylor Expansion.Under 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 more than 0.2K.The relation that the present invention is directly established between two Detection Using Thermal Infrared Channel radiances, it is general instead of establishing
The bright gram of relation between radiance and brightness temperature, the conversion accuracy between two passage radiances is improved, is established
Relation between surface temperature and two Detection Using Thermal Infrared Channel radiances, propose a Split window algorithms (letter based on radiance
Claim RBSWA), and it is applied to MODIS data Surface Temperature Retrievals.
The content of the invention
The present invention proposes one and splits window Surface Temperature Retrieval method based on radiance, it is possible to achieve MODIS data
Surface Temperature Retrieval, method are as follows:
Wherein B31(TS)=A0L31-A1L31+A2(formula 1)
In formula, TSRepresent surface temperature, B31(TS) expression temperature is TSBlack matrix corresponding to the wave bands of MODIS 31 radiation it is bright
Degree, K31,1And K31,2The thermal infrared radiation constant (can be obtained from metadata) of the wave bands of MODIS 31, L are represented respectively31And L32Point
Biao Shi not radiance, A on the star of the wave bands 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) represent to convert related parameter, C between two Detection Using Thermal Infrared Channel radiances31、C32、
D31And D32Represent intermediate variable.k(L31) and b (L31) calculation formula is as follows:
In formula, K32,1And K32,2The thermal infrared radiation constant (can be obtained from metadata) of the wave bands of MODIS 32 is represented respectively.
C31、C32、D31And D32Calculation formula is as follows:
C31=ε31τ31
C32=ε32τ32
In formula, ε31And ε32The Land surface emissivity of the wave bands of MODIS 31 and 32, τ are represented respectively31And τ32Represent respectively
The air effective transmissivity of the wave bands of MODIS 31 and 32,WithThe up spoke of air of the wave bands of MODIS 31 and 32 is represented respectively
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 represents 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.
Embodiment
1st, radiance on star is calculated, radiance L on star31And L32It can be calibrated to obtain via radiation by image DN value;
2nd, Land surface emissivity, Land surface emissivity ε are calculated31And ε32Can use NDVI threshold methods [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;
3rd, 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] 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] obtained from Thermal Infrared Data inverting;
4th, surface temperature is calculated, surface temperature is calculated (such as using the Split window algorithms proposed by the present invention based on radiance
Formula 1).
5th, arithmetic accuracy is assessed.Calculated based on simulated data sets evaluation RBSWA model accuracy, and with two classical windows of splitting
Method contrasts (Split window algorithms of the proposition 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, the input of θ and w as Split window algorithms, then compare Split window algorithms
The LST that the LST that inverting obtains concentrates with analogue data.Arithmetic accuracy assessment result such as table 1, show to calculate with the window of splitting that two classical
Method contrasts, and RBSWA Surface Temperature Retrieval precision improves about 0.5K.
The different Split window algorithms Surface Temperature Retrieval ratios of precision of table 1 compared with
Claims (1)
1. one is split window Surface Temperature Retrieval algorithm based on radiance, it is characterized in radiance on star directly as calculation
The input of method;Step is:
The first step, radiance on star is calculated, radiance L on star31And L32It can be calibrated to obtain via radiation by image DN value;
Second step, calculate Land surface emissivity, Land surface emissivity ε31And ε32Can use NDVI threshold methods [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;
3rd step, 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] 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-8 TIRS data over land area.Remote Sensing
Letters, 6,904-913] obtained from Thermal Infrared Data inverting;
4th step, surface temperature is calculated, surface temperature is calculated (such as using the Split window algorithms proposed by the present invention based on radiance
Formula 1 in specification).
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 true CN107389198A (en) | 2017-11-24 |
CN107389198B 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) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109297605A (en) * | 2018-10-09 | 2019-02-01 | 北京大学 | A kind of Surface Temperature Retrieval method with Thermal Infrared Data infrared based in |
CN109596223A (en) * | 2018-12-24 | 2019-04-09 | 中国农业科学院农业资源与农业区划研究所 | A kind of multispectral satellite borne sensor lane device and method of direct inversion surface temperature |
CN109632106A (en) * | 2019-01-03 | 2019-04-16 | 北京师范大学 | A kind of Remote Sensing temperature angles of product effect correction method |
CN109887615A (en) * | 2019-01-30 | 2019-06-14 | 北京环境特性研究所 | Surface temperature period diurnal variation analogy method |
CN114894321A (en) * | 2022-07-13 | 2022-08-12 | 国家卫星气象中心(国家空间天气监测预警中心) | Calibration method of infrared remote sensing instrument, electronic device and computer storage medium |
Citations (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 |
CN104748857A (en) * | 2015-03-05 | 2015-07-01 | 北京师范大学 | Method and system for inverting urban surface temperatures |
CN104897289A (en) * | 2015-06-23 | 2015-09-09 | 中国科学院遥感与数字地球研究所 | Landsat 8 satellite data land surface temperature inversion method |
US9551619B1 (en) * | 2011-09-23 | 2017-01-24 | Rockwell Collins, Inc. | Terahertz laser |
-
2017
- 2017-05-23 CN CN201710372227.2A patent/CN107389198B/en not_active Expired - Fee Related
Patent Citations (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 |
CN104897289A (en) * | 2015-06-23 | 2015-09-09 | 中国科学院遥感与数字地球研究所 | Landsat 8 satellite data land surface temperature inversion method |
Non-Patent Citations (2)
Title |
---|
秦福莹: "热红外遥感地表温度反演方法应用与对比分析研究", 《中国优秀硕士学位论文全文数据库》 * |
陈瀚阅 等: "MODIS数据地表温度反演劈窗算法比较", 《遥感技术与应用》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109297605A (en) * | 2018-10-09 | 2019-02-01 | 北京大学 | A kind of Surface Temperature Retrieval method with Thermal Infrared Data infrared based in |
CN109596223A (en) * | 2018-12-24 | 2019-04-09 | 中国农业科学院农业资源与农业区划研究所 | A kind of multispectral satellite borne sensor lane device and method of direct inversion surface temperature |
CN109632106A (en) * | 2019-01-03 | 2019-04-16 | 北京师范大学 | A kind of Remote Sensing temperature angles of product effect correction method |
CN109632106B (en) * | 2019-01-03 | 2020-06-30 | 北京师范大学 | Remote sensing surface temperature product angle effect correction method |
CN109887615A (en) * | 2019-01-30 | 2019-06-14 | 北京环境特性研究所 | Surface temperature period diurnal variation analogy method |
CN109887615B (en) * | 2019-01-30 | 2020-12-11 | 北京环境特性研究所 | Earth surface temperature periodic daily change simulation method |
CN114894321A (en) * | 2022-07-13 | 2022-08-12 | 国家卫星气象中心(国家空间天气监测预警中心) | Calibration method of infrared remote sensing instrument, electronic device and computer storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN107389198B (en) | 2019-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107389198B (en) | It is a kind of that window Surface Temperature Retrieval method is split based on radiance | |
Schmugge et al. | Recovering surface temperature and emissivity from thermal infrared multispectral data | |
Sobrino et al. | Land surface temperature retrieval from MSG1-SEVIRI data | |
Kerr et al. | Accurate land surface temperature retrieval from AVHRR data with use of an improved split window algorithm | |
Liang | An optimization algorithm for separating land surface temperature and emissivity from multispectral thermal infrared imagery | |
Hall et al. | Comparison of satellite-derived and in-situ observations of ice and snow surface temperatures over Greenland | |
CN108020322B (en) | Airborne thermal infrared hyperspectral remote sensing quantitative detection method for coal field fire area | |
Ingram et al. | Sensitivity of iterative spectrally smooth temperature/emissivity separation to algorithmic assumptions and measurement noise | |
CN103744069A (en) | Methane profile orthogonal inversion method aiming at AIRS (atmospheric infrared sounder) hyper-spectrum satellite data | |
CN104897289B (en) | A kind of satellite data Surface Temperature Retrieval methods of Landsat 8 | |
Liu et al. | Spectral calibration of hyperspectral data observed from a hyperspectrometer loaded on an unmanned aerial vehicle platform | |
CN111707376A (en) | Surface temperature inversion method for broadband thermal infrared sensor | |
Wang et al. | An enhanced single‐channel algorithm for retrieving land surface temperature from Landsat series data | |
CN114838827A (en) | Earth surface temperature inversion channel selection method based on MERSI-II remote sensing data | |
Wang et al. | Estimation of surface longwave radiation over the Tibetan plateau region using MODIS data for cloud-free skies | |
Danodia et al. | Land surface temperature retrieval by radiative transfer equation and single channel algorithms using landsat-8 satellite data | |
Kwon et al. | Use of total precipitable water classification of a priori error and quality control in atmospheric temperature and water vapor sounding retrieval | |
Chekroun et al. | THE APPRAISAL OF NOAA SATELLITES LST-SW ALGORITHMS: NOAA-20 (JPSS-1) PROPOSAL | |
Li et al. | An improved single-channel polar region ice surface temperature retrieval algorithm using landsat-8 data | |
Olioso et al. | Impact of surface emissivity and atmospheric conditions on surface temperatures estimated from top of canopy brightness temperatures derived from Landsat 7 data | |
Qian et al. | Land surface temperature retrieval using nighttime mid-infrared channels data from airborne hyperspectral scanner | |
Xie et al. | Progress of remote sensing inversion research on surface temperature of lakes and reservoirs | |
Xia et al. | Retrieval land surface temperature from visible infrared imager radiometer suite data | |
Peng et al. | Estimation of evaporative fraction from top-of-atmosphere radiance | |
Shi et al. | Atmospheric correction of passive microwave brightness temperature on the estimation of snow depth |
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 |
Granted publication date: 20190920 Termination date: 20210523 |
|
CF01 | Termination of patent right due to non-payment of annual fee |