CN117113855A - 一种单通道地表温度遥感方法及其系统 - Google Patents
一种单通道地表温度遥感方法及其系统 Download PDFInfo
- Publication number
- CN117113855A CN117113855A CN202311362956.1A CN202311362956A CN117113855A CN 117113855 A CN117113855 A CN 117113855A CN 202311362956 A CN202311362956 A CN 202311362956A CN 117113855 A CN117113855 A CN 117113855A
- Authority
- CN
- China
- Prior art keywords
- channel
- thermal infrared
- remote sensing
- inversion
- surface temperature
- 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 55
- 210000001747 pupil Anatomy 0.000 claims abstract description 68
- 238000000638 solvent extraction Methods 0.000 claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims description 30
- 238000005316 response function Methods 0.000 claims description 22
- 238000004088 simulation Methods 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 7
- 230000005855 radiation Effects 0.000 claims description 7
- 238000005192 partition Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 description 14
- 230000009286 beneficial effect Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 239000002689 soil Substances 0.000 description 3
- 241000132092 Aster Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- GMTYREVWZXJPLF-AFHUBHILSA-N butorphanol D-tartrate Chemical compound OC(=O)[C@@H](O)[C@H](O)C(O)=O.N1([C@@H]2CC3=CC=C(C=C3[C@@]3([C@]2(CCCC3)O)CC1)O)CC1CCC1 GMTYREVWZXJPLF-AFHUBHILSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- -1 vegetation Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明公开了一种单通道地表温度遥感方法及其系统，属于单通道地表温度遥感领域，该方法包括构建单通道遥感反演基础模型；进行反演状况分区，并基于单通道遥感反演基础模型，得到通道辐亮度到亮度温度的查找表、各反演状况区间的单通道遥感反演模型和通用单通道遥感反演模型；根据通道辐亮度到亮度温度的查找表和通用单通道遥感反演模型，得到当前热红外通道的入瞳亮温和地表温度初始值；根据当前热红外通道的入瞳亮温、各反演状况区间的单通道遥感反演模型和地表温度初始值，更新不同反演状况下的地表温度，实现单通道地表温度遥感反演。本发明解决了现有地表温度遥感方法受大气影响，精度不够的问题。
Description
技术领域
本发明属于单通道地表温度遥感领域，尤其涉及一种单通道地表温度遥感方法及其系统。
背景技术
地表温度是表征地气之间的相互作用及能量交换的关键特征物理量，在陆表生态系统研究、全球气候变化、地质勘探等领域具有重要应用。热红外遥感是获取大尺度地表温度的主要方法，然而，由于地表和大气参数的耦合，地表温度的遥感反演存在一定的困难，要从卫星传感器的观测数据中反演得到地表温度，需要消除大气和地表发射率的影响。
对于多光谱遥感数据，最常用的两种反演方法是单通道法和分裂窗法。单通道法仅需要一个热红外通道，但需要已知大气参数和发射率。大气参数的获取通常基于大气再分析数据和MODTRAN等成熟的辐射传输模型，然而，大气再分析数据的获取也较为困难，且在潮湿大气状况下单通道法的精度欠佳。相比之下，分裂窗法基于相邻两个热红外通道的大气吸收差异来消除大气的影响，只需要发射率即可反演得到高精度的地表温度。但是，分裂窗法需要两个热红外通道，在一定程度上限制了其应用。
发明内容
针对现有技术中的上述不足，本发明提供的一种单通道地表温度遥感方法及其系统解决了现有地表温度遥感方法受大气影响，精度不够的问题。
为了达到上述发明目的，本发明采用的技术方案为：一种单通道地表温度遥感方法，包括以下步骤：
S1、构建单通道遥感反演基础模型；
S2、进行反演状况分区，并基于单通道遥感反演基础模型，得到通道辐亮度到亮度温度的查找表、各反演状况区间的单通道遥感反演模型和通用单通道遥感反演模型；
S3、根据通道辐亮度到亮度温度的查找表和通用单通道遥感反演模型，得到当前热红外通道的入瞳亮温和地表温度初始值；
S4、根据当前热红外通道的入瞳亮温、各反演状况区间的单通道遥感反演模型和地表温度初始值，更新不同反演状况下的地表温度，实现单通道地表温度遥感反演。
本发明的有益效果为：本发明使用一个热红外通道的入瞳亮温和两个通道的发射率来反演地表温度。基于本方法，可以规避大气效应的校正过程，直接从单个热红外通道的入瞳观测数据中反演得到地表温度，并提高了地表温度反演精度。
进一步地，所述步骤S1中单通道遥感反演基础模型的表达式为：
其中，为地表温度；/>为相邻的两个热红外通道/>和/>的发射率的平均值；/>为单通道遥感反演基础模型中热红外通道/>的发射率；/>为单通道遥感反演基础模型中热红外通道/>的发射率；/>为相邻的两个热红外通道/>和/>的发射率的差值；/>、/>、/>、/>、/>、和/>均为单通道遥感反演基础模型系数；/>为单通道遥感反演基础模型中热红外通道的入瞳亮温。
上述进一步方案的有益效果为：利用单个通道的入瞳观测数据和两个通道发射率的线性组合消除大气效应的影响，最终构建无需大气效应校正的单通道遥感反演基础模型，避免大气效应的校正过程。
进一步地，所述步骤S2具体为：
S201、获取配套的大气温湿廓线、地表温度和地表发射率波谱值，并根据配套的大气温湿廓线、地表温度和地表发射率波谱值，生成模拟数据；
S202、根据模拟数据，利用MODTRAN辐射传输模型，得到入瞳辐亮度波谱值；
S203、获取热红外通道的通道响应函数、热红外通道/>的通道响应函数和热红外通道/>的通道响应函数；
S204、根据入瞳辐亮度波谱值、热红外通道的通道响应函数、热红外通道/>的通道响应函数、热红外通道/>的通道响应函数和地表发射率波谱值，得到热红外通道/>的通道辐亮度、热红外通道/>的发射率和热红外通道/>的发射率：
其中，为热红外通道/>的通道辐亮度；/>为热红外通道/>的发射率；/>为热红外通道/>的发射率；/>、/>和/>分别为热红外通道/>、热红外通道/>和热红外通道的通道响应函数；/>和/>分别为热红外通道/>的起始波长和终止波长；/>和/>分别为热红外通道/>的起始波长和终止波长；/>和/>分别为热红外通道/>的起始波长和终止波长；/>为波长；/>为/>的微元；/>为入瞳辐亮度波谱值；/>为地表发射率波谱值；
S205、根据模拟数据，建立通道辐亮度到亮度温度的查找表；
S206、根据热红外通道的通道辐亮度，利用通道辐亮度到亮度温度的查找表，得到热红外通道/>的入瞳亮温；
S207、根据热红外通道的入瞳亮温、热红外通道/>的发射率、热红外通道/>的发射率和地表温度，利用最小二乘法，拟合得到通用单通道遥感反演模型系数；
S208、将通用单通道遥感反演模型系数作为单通道遥感反演基础模型系数代入单通道遥感反演基础模型，得到通用单通道遥感反演模型；
S209、进行反演状况分区，得到若干个反演状况区间，并拟合得到各反演状况区间的单通道遥感反演模型系数；
S2010、将各反演状况区间的单通道遥感反演模型系数作为单通道遥感反演基础模型系数，利用单通道遥感反演基础模型，得到各反演状况区间的单通道遥感反演模型。
上述进一步方案的有益效果为：拟合得到的通用单通道遥感反演模型系数不考虑天气和地表条件，为后续地表温度初始值的求解做准备；并且反演状况分区时不同的大气和地表状况允许区间存在部分重叠，抵消了分区选择错误对反演结果的影响。
进一步地，所述步骤S209中反演状况分区的标准为：
地表温度划分为和/>两个区间；
各热红外通道的发射率划分为和/>两个区间；
总水汽含量划分为、/>和/>三个区间；
观测角度划分为、/>、/>和/>四个区间。
上述进一步方案的有益效果为：反演状况分区时不同的大气和地表状况允许区间存在部分重叠，抵消了分区选择错误对反演结果的影响。
进一步地，所述步骤S3具体为：
S301、获取待反演单通道传感器热红外通道的定标系数；
S302、根据待反演单通道传感器热红外通道的定标系数，将传感器记录的计数值转为热红外通道/>的入瞳辐亮度；
S303、根据通道辐亮度到亮度温度的查找表，将热红外通道的入瞳辐亮度转化为入瞳亮温，得到当前热红外通道/>的入瞳亮温；
S304、获取地表覆盖分类；
S305、根据地表覆盖分类，使用分类赋值法，对每类地物在热红外通道和热红外通道/>的发射率进行赋值，得到各类地物热红外通道/>的发射率和各类地物热红外通道/>的发射率；
S306、将当前热红外通道的入瞳亮温作为通用单通道遥感反演模型中热红外通道/>的入瞳亮温、将各类地物热红外通道/>的发射率作为通用单通道遥感反演模型中热红外通道/>的发射率，以及将各类地物热红外通道/>的发射率作为通用单通道遥感反演模型中热红外通道/>的发射率，代入通用单通道遥感反演模型，计算得到地表温度初始值。
上述进一步方案的有益效果为：基于通用单通道遥感反演模型系数，求解得到地表温度初始值，为后续求解不同反演状况下的地表温度做准备，规避了大气效应的校正过程，并提高了地表温度反演精度。
进一步地，所述步骤S4具体为：
S401、获取卫星过境时刻大气的总水汽含量、当前各热红外通道的发射率和当前观测角度；
S402、根据卫星过境时刻大气的总水汽含量、当前各热红外通道的发射率、当前观测角度和地表温度初始值，确定反演状况区间，得到区间结果；
S403、根据区间结果和各反演状况区间的单通道遥感反演模型，得到当前反演状况区间的单通道遥感反演模型；
S404、根据当前热红外通道的入瞳亮温和当前各热红外通道的发射率，利用当前反演状况区间的单通道遥感反演模型，得到当前反演状况下的地表温度，实现单通道地表温度遥感反演。
上述进一步方案的有益效果为：使用一个热红外通道的入瞳亮温和两个热红外通道的发射率来反演地表温度，求解不同反演状况下的地表温度，规避了大气效应的校正过程，直接从单个热红外通道的入瞳观测数据中反演得到地表温度，并提高了地表温度反演精度。
本发明提供了一种单通道地表温度遥感系统，包括模型构建模块、系数率定模块、初值估算模块和地表温度反演模块；
所述模型构建模块，用于构建单通道遥感反演基础模型；
所述系数率定模块，用于获取模拟数据，进行反演状况分区，并根据模拟数据和单通道遥感反演基础模型，利用最小二乘法，得到通道辐亮度到亮度温度的查找表、各反演状况区间的单通道遥感反演模型和通用单通道遥感反演模型；
所述初值估算模块，用于获取待反演单通道传感器热红外通道的定标系数，并根据待反演单通道传感器热红外通道/>的定标系数、通道辐亮度到亮度温度的查找表和通用单通道遥感反演模型，得到当前热红外通道/>的入瞳亮温和地表温度初始值；
所述地表温度反演模块，用于根据当前热红外通道的入瞳亮温、各反演状况区间的单通道遥感反演模型和地表温度初始值，更新不同反演状况下的地表温度，实现单通道地表温度遥感反演。
本发明的有益效果为：本系统利用单通道遥感反演基础模型，通过一个热红外通道的入瞳亮温和两个通道的发射率来反演地表温度，可以规避大气效应的校正过程，并提高了地表温度反演精度。
附图说明
图1为本发明的方法流程图。
图2为本发明的系统结构图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
实施例1
如图1所示，在本发明的一个实施例中，一种单通道地表温度遥感方法，包括以下步骤：
S1、构建单通道遥感反演基础模型；
S2、进行反演状况分区，并基于单通道遥感反演基础模型，得到通道辐亮度到亮度温度的查找表、各反演状况区间的单通道遥感反演模型和通用单通道遥感反演模型；
S3、根据通道辐亮度到亮度温度的查找表和通用单通道遥感反演模型，得到当前热红外通道的入瞳亮温和地表温度初始值；
S4、根据当前热红外通道的入瞳亮温、各反演状况区间的单通道遥感反演模型和地表温度初始值，更新不同反演状况下的地表温度，实现单通道地表温度遥感反演。
本实施例中，本发明提出了一种无需大气效应校正的单通道地表温度遥感方法。基于入瞳辐亮度的通道相关性假设，使用一个热红外通道的入瞳亮温和两个通道的地表发射率来反演地表温度。基于本方法及系统，可以规避大气效应的校正过程，直接从单个热红外通道的入瞳观测数据中反演得到地表温度。针对仅含有单个热红外通道的卫星传感器，基于通道相关性假设，构造单通道地表温度遥感反演模型，以提高该类传感器地表温度反演的精度，推动其高效应用。
所述步骤S1中单通道遥感反演基础模型的表达式为：
其中，为地表温度；/>为相邻的两个热红外通道/>和/>的发射率的平均值；/>为单通道遥感反演基础模型中热红外通道/>的发射率；/>为单通道遥感反演基础模型中热红外通道/>的发射率；/>为相邻的两个热红外通道/>和/>的发射率的差值；/>、/>、/>、/>、/>、和/>均为单通道遥感反演基础模型系数；/>为单通道遥感反演基础模型中热红外通道的入瞳亮温。
本实施例中，对于单通道而言，考虑通道间亮温的相关性，可以将相似区间通道的亮温表示为通道/>的亮温的线性函数，进而可以实现从通道/>至虚拟的通道/>和/>的转换。例如，Landsat 7的第6通道和Landsat 8/9的第10和11通道的亮温数据具有较好的线性关系（Landsat 7和Landsat 8/9均为陆地卫星名），可以用作构建单通道遥感反演基础模型的/>和/>和/>通道。同时，基于热红外辐射传输方程，以及虚拟通道/>和通道/>大气吸收的差异，利用两个虚拟通道的线性组合消除大气效应的影响，最终构建无需大气效应校正的单通道遥感反演模型。
所述步骤S2具体为：
S201、获取配套的大气温湿廓线、地表温度和地表发射率波谱值，并根据配套的大气温湿廓线、地表温度和地表发射率波谱值，生成模拟数据；
S202、根据模拟数据，利用MODTRAN辐射传输模型，得到入瞳辐亮度波谱值；
S203、获取热红外通道的通道响应函数、热红外通道/>的通道响应函数和热红外通道/>的通道响应函数；
S204、根据入瞳辐亮度波谱值、热红外通道的通道响应函数、热红外通道/>的通道响应函数、热红外通道/>的通道响应函数和地表发射率波谱值，得到热红外通道/>的通道辐亮度、热红外通道/>的发射率和热红外通道/>的发射率：
其中，为热红外通道/>的通道辐亮度；/>为热红外通道/>的发射率；/>为热红外通道/>的发射率；/>、/>和/>分别为热红外通道/>、热红外通道/>和热红外通道的通道响应函数；/>和/>分别为热红外通道/>的起始波长和终止波长；/>和/>分别为热红外通道/>的起始波长和终止波长；/>和/>分别为热红外通道/>的起始波长和终止波长；/>为波长；/>为/>的微元；/>为入瞳辐亮度波谱值；/>为地表发射率波谱值；
S205、根据模拟数据，建立通道辐亮度到亮度温度的查找表；
S206、根据热红外通道的通道辐亮度，利用通道辐亮度到亮度温度的查找表，得到热红外通道/>的入瞳亮温；
S207、根据热红外通道的入瞳亮温、热红外通道/>的发射率、热红外通道/>的发射率和地表温度，利用最小二乘法，拟合得到通用单通道遥感反演模型系数；
S208、将通用单通道遥感反演模型系数作为单通道遥感反演基础模型系数代入单通道遥感反演基础模型，得到通用单通道遥感反演模型；
S209、进行反演状况分区，得到若干个反演状况区间，并拟合得到各反演状况区间的单通道遥感反演模型系数；
S2010、将各反演状况区间的单通道遥感反演模型系数作为单通道遥感反演基础模型系数，利用单通道遥感反演基础模型，得到各反演状况区间的单通道遥感反演模型。
本实施例中，基于模拟数据，采用最小二乘法实现无需大气效应校正的单通道遥感反演模型系数的率定。
生成模拟数据需要获取配套的大气温湿廓线、地表温度和地表发射率波谱值。获取典型大气温湿廓线数据，例如从MODTRAN辐射传输模型中挑选晴空的6条标准大气或者从TIGR 2000标准大气廓线数据库中挑选大气总水汽含量均匀分布的100条晴空大气廓线；获取典型地类的地表发射率波谱数据，例如从ASTER光谱库中挑选岩石、植被、土壤、水体和人造地类的地表发射率波谱；设置不同地气状况的下的地表温度，例如根据待选的大气廓线数据，以其大气底层气温为基准，以5K为步长，在大气底层气温±15K的范围内扰动设置所需的地表温度值。
在上述地表温度、地表发射率波谱值和大气温湿廓线获取的基础上，将其输入到MODTRAN辐射传输模型中，获取到入瞳辐亮度波谱值。
获取待反演单通道传感器的通道响应函数。同时将入瞳辐亮度波谱值和地表发射率波谱值转换为通道等效值、/>和/>。
同时根据模拟数据，建立通道辐亮度到亮度温度的查找表，将入瞳通道辐亮度转换为入瞳亮温。
基于上述步骤，可以得到和/>通道的发射率、/>通道的入瞳亮温及地表温度，进而可以拟合得到单通道遥感反演模型系数/>(/>=0,1, 2, …, 6)。
接着，拟合得到不考虑天气和地表条件的通用单通道遥感反演模型系数，例如，对于Landsat 7传感器而言，通用模型的系数为：，/>，/>，/>，，/>，/>。进一步地，针对不同反演状况，分别拟合不同状况下的单通道遥感反演模型系数。
所述步骤S209中反演状况分区的标准为：
地表温度划分为和/>两个区间；
各热红外通道的发射率划分为和/>两个区间；
总水汽含量划分为、/>和/>三个区间；
观测角度划分为、/>、/>和/>四个区间。
本实施例中，反演状况区间可以按地表温度、各热红外通道的发射率、总水汽含量、观测角度进行划分反演状况区间。不同的大气和地表状况允许区间存在部分重叠，以抵消分区选择错误对反演结果的影响。针对不同反演状况区间分别构建对应的模型系数。例如，针对地表温度>290K，和/>通道发射率>0.96，总水汽含量（WVC）小于/>，且观测角度0-30°的反演状况，拟合该条件下的模型系数。例如，对于Landsat 7，部分反演状况的模型系数如表1所示。
表1
所述步骤S3具体为：
S301、获取待反演单通道传感器热红外通道的定标系数；
S302、根据待反演单通道传感器热红外通道的定标系数，将传感器记录的计数值转为热红外通道/>的入瞳辐亮度；
S303、根据通道辐亮度到亮度温度的查找表，将热红外通道的入瞳辐亮度转化为入瞳亮温，得到当前热红外通道/>的入瞳亮温；
S304、获取地表覆盖分类；
S305、根据地表覆盖分类，使用分类赋值法，对每类地物在热红外通道和热红外通道/>的发射率进行赋值，得到各类地物热红外通道/>的发射率和各类地物热红外通道/>的发射率；
S306、将当前热红外通道的入瞳亮温作为通用单通道遥感反演模型中热红外通道/>的入瞳亮温、将各类地物热红外通道/>的发射率作为通用单通道遥感反演模型中热红外通道/>的发射率，以及将各类地物热红外通道/>的发射率作为通用单通道遥感反演模型中热红外通道/>的发射率，代入通用单通道遥感反演模型，计算得到地表温度初始值。
本实施例中，根据待反演的单通道传感器的定标系数，将传感器记录的计数值转变为入瞳通道辐亮度值，再转为入瞳亮温。例如，针对Landsat 7的第6通道而言，根据传感器数据的头文件中的增益和偏置，以及通道辐亮度到亮度温度的查找表，将DN值（DN值为遥感影像像元亮度值）转换为入瞳亮温。
使用分类赋值法进行发射率的确定。获取地表覆盖分类，对每类地物进行通道和通道的发射率的赋值。例如，针对Landsat 8/9的两个热红外通道（10和11通道），植被在10通道的发射率设置为0.986、在11通道的发射率设置为0.990，土壤在10通道的发射率设置为0.965、在11通道的发射率设置为0.975，水体在10通道的发射率设置为0.990，在11通道的发射率设置为0.991。
将单通道的入瞳亮温和双通道/>和/>的发射率代入单通道遥感反演基础模型，使用不考虑天气和地表条件的通用单通道遥感反演模型系数，求解得到地表温度初始值。
所述步骤S4具体为：
S401、获取卫星过境时刻大气的总水汽含量、当前各热红外通道的发射率和当前观测角度；
S402、根据卫星过境时刻大气的总水汽含量、当前各热红外通道的发射率、当前观测角度和地表温度初始值，确定反演状况区间，得到区间结果；
S403、根据区间结果和各反演状况区间的单通道遥感反演模型，得到当前反演状况区间的单通道遥感反演模型；
S404、根据当前热红外通道的入瞳亮温和当前各热红外通道的发射率，利用当前反演状况区间的单通道遥感反演模型，得到当前反演状况下的地表温度，实现单通道地表温度遥感反演。
本实施例中，获取卫星过境时刻的大气总水汽含量信息，根据地表温度初始值、发射率和大气总水汽含量数据，确定反演状况分区，带入对应区间的反演模型，更新不同反演状况下的地表温度，完成地表温度的迭代求解。
实施例2
如图2所示，本发明提供了一种单通道地表温度遥感系统，包括模型构建模块、系数率定模块、初值估算模块和地表温度反演模块；
所述模型构建模块，用于构建单通道遥感反演基础模型；
所述系数率定模块，用于获取模拟数据，进行反演状况分区，并根据模拟数据和单通道遥感反演基础模型，利用最小二乘法，得到通道辐亮度到亮度温度的查找表、各反演状况区间的单通道遥感反演模型和通用单通道遥感反演模型；
所述初值估算模块，用于获取待反演单通道传感器热红外通道的定标系数，并根据待反演单通道传感器热红外通道/>的定标系数、通道辐亮度到亮度温度的查找表和通用单通道遥感反演模型，得到当前热红外通道/>的入瞳亮温和地表温度初始值；
所述地表温度反演模块，用于根据当前热红外通道的入瞳亮温、各反演状况区间的单通道遥感反演模型和地表温度初始值，更新不同反演状况下的地表温度，实现单通道地表温度遥感反演。
本实施例中，模型构建模块，用于确定热红外通道和热红外通道/>及热红外通道的位置及关系，构建单通道遥感反演基础模型。
系数率定模块，用于获取模拟数据，并基于最小二乘法实现无需大气效应校正的单通道遥感反演模型系数的率定。
初值估算模块，用于获取目标传感器的单通道星上亮温，基于分类赋值法计算双通道/>和/>的发射率，进而根据通用单通道遥感反演模型系数估算地表温度初始值。
地表温度反演模块，用于获取大气总水汽含量，根据温度初始值、发射率、大气总水汽含量和观测角度，判断反演状况，进而，利用不同反演状况区间对应的模型系数，实现不同反演状况下的地表温度估算。
Claims (7)
1.一种单通道地表温度遥感方法，其特征在于，包括以下步骤：
S1、构建单通道遥感反演基础模型；
S2、进行反演状况分区，并基于单通道遥感反演基础模型，得到通道辐亮度到亮度温度的查找表、各反演状况区间的单通道遥感反演模型和通用单通道遥感反演模型；
S3、根据通道辐亮度到亮度温度的查找表和通用单通道遥感反演模型，得到当前热红外通道的入瞳亮温和地表温度初始值；
S4、根据当前热红外通道的入瞳亮温、各反演状况区间的单通道遥感反演模型和地表温度初始值，更新不同反演状况下的地表温度，实现单通道地表温度遥感反演。
2.根据权利要求1所述单通道地表温度遥感方法，其特征在于，所述步骤S1中单通道遥感反演基础模型的表达式为：
其中，为地表温度；/>为相邻的两个热红外通道/>和/>的发射率的平均值；/>为单通道遥感反演基础模型中热红外通道/>的发射率；/>为单通道遥感反演基础模型中热红外通道/>的发射率；/>为相邻的两个热红外通道/>和/>的发射率的差值；/>、/>、/>、/>、/>、/>和均为单通道遥感反演基础模型系数；/>为单通道遥感反演基础模型中热红外通道/>的入瞳亮温。
3.根据权利要求1所述单通道地表温度遥感方法，其特征在于，所述步骤S2具体为：
S201、获取配套的大气温湿廓线、地表温度和地表发射率波谱值，并根据配套的大气温湿廓线、地表温度和地表发射率波谱值，生成模拟数据；
S202、根据模拟数据，利用MODTRAN辐射传输模型，得到入瞳辐亮度波谱值；
S203、获取热红外通道的通道响应函数、热红外通道/>的通道响应函数和热红外通道的通道响应函数；
S204、根据入瞳辐亮度波谱值、热红外通道的通道响应函数、热红外通道/>的通道响应函数、热红外通道/>的通道响应函数和地表发射率波谱值，得到热红外通道/>的通道辐亮度、热红外通道/>的发射率和热红外通道/>的发射率：
其中，为热红外通道/>的通道辐亮度；/>为热红外通道/>的发射率；/>为热红外通道的发射率；/>、/>和/>分别为热红外通道/>、热红外通道/>和热红外通道/>的通道响应函数；/>和/>分别为热红外通道/>的起始波长和终止波长；/>和/>分别为热红外通道/>的起始波长和终止波长；/>和/>分别为热红外通道/>的起始波长和终止波长；/>为波长；/>为/>的微元；/>为入瞳辐亮度波谱值；/>为地表发射率波谱值；
S205、根据模拟数据，建立通道辐亮度到亮度温度的查找表；
S206、根据热红外通道的通道辐亮度，利用通道辐亮度到亮度温度的查找表，得到热红外通道/>的入瞳亮温；
S207、根据热红外通道的入瞳亮温、热红外通道/>的发射率、热红外通道/>的发射率和地表温度，利用最小二乘法，拟合得到通用单通道遥感反演模型系数；
S208、将通用单通道遥感反演模型系数作为单通道遥感反演基础模型系数代入单通道遥感反演基础模型，得到通用单通道遥感反演模型；
S209、进行反演状况分区，得到若干个反演状况区间，并拟合得到各反演状况区间的单通道遥感反演模型系数；
S2010、将各反演状况区间的单通道遥感反演模型系数作为单通道遥感反演基础模型系数，利用单通道遥感反演基础模型，得到各反演状况区间的单通道遥感反演模型。
4.根据权利要求3所述单通道地表温度遥感方法，其特征在于，所述步骤S209中反演状况分区的标准为：
地表温度划分为和/>两个区间；
各热红外通道的发射率划分为和/>两个区间；
总水汽含量划分为、/>和/>三个区间；
观测角度划分为、/>、/>和/>四个区间。
5.根据权利要求1所述单通道地表温度遥感方法，其特征在于，所述步骤S3具体为：
S301、获取待反演单通道传感器热红外通道的定标系数；
S302、根据待反演单通道传感器热红外通道的定标系数，将传感器记录的计数值转为热红外通道/>的入瞳辐亮度；
S303、根据通道辐亮度到亮度温度的查找表，将热红外通道的入瞳辐亮度转化为入瞳亮温，得到当前热红外通道/>的入瞳亮温；
S304、获取地表覆盖分类；
S305、根据地表覆盖分类，使用分类赋值法，对每类地物在热红外通道和热红外通道/>的发射率进行赋值，得到各类地物热红外通道/>的发射率和各类地物热红外通道/>的发射率；
S306、将当前热红外通道的入瞳亮温作为通用单通道遥感反演模型中热红外通道/>的入瞳亮温、将各类地物热红外通道/>的发射率作为通用单通道遥感反演模型中热红外通道/>的发射率，以及将各类地物热红外通道/>的发射率作为通用单通道遥感反演模型中热红外通道/>的发射率，代入通用单通道遥感反演模型，计算得到地表温度初始值。
6.根据权利要求4所述单通道地表温度遥感方法，其特征在于，所述步骤S4具体为：
S401、获取卫星过境时刻大气的总水汽含量、当前各热红外通道的发射率和当前观测角度；
S402、根据卫星过境时刻大气的总水汽含量、当前各热红外通道的发射率、当前观测角度和地表温度初始值，确定反演状况区间，得到区间结果；
S403、根据区间结果和各反演状况区间的单通道遥感反演模型，得到当前反演状况区间的单通道遥感反演模型；
S404、根据当前热红外通道的入瞳亮温和当前各热红外通道的发射率，利用当前反演状况区间的单通道遥感反演模型，得到当前反演状况下的地表温度，实现单通道地表温度遥感反演。
7.根据权利要求1-6任一所述单通道地表温度遥感方法的单通道地表温度遥感系统，其特征在于，包括模型构建模块、系数率定模块、初值估算模块和地表温度反演模块；
所述模型构建模块，用于构建单通道遥感反演基础模型；
所述系数率定模块，用于获取模拟数据，进行反演状况分区，并根据模拟数据和单通道遥感反演基础模型，利用最小二乘法，得到通道辐亮度到亮度温度的查找表、各反演状况区间的单通道遥感反演模型和通用单通道遥感反演模型；
所述初值估算模块，用于获取待反演单通道传感器热红外通道的定标系数，并根据待反演单通道传感器热红外通道/>的定标系数、通道辐亮度到亮度温度的查找表和通用单通道遥感反演模型，得到当前热红外通道/>的入瞳亮温和地表温度初始值；
所述地表温度反演模块，用于根据当前热红外通道的入瞳亮温、各反演状况区间的单通道遥感反演模型和地表温度初始值，更新不同反演状况下的地表温度，实现单通道地表温度遥感反演。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311362956.1A CN117113855B (zh) | 2023-10-20 | 2023-10-20 | 一种单通道地表温度遥感方法及其系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311362956.1A CN117113855B (zh) | 2023-10-20 | 2023-10-20 | 一种单通道地表温度遥感方法及其系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117113855A true CN117113855A (zh) | 2023-11-24 |
| CN117113855B CN117113855B (zh) | 2023-12-26 |
Family
ID=88805828
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311362956.1A Active CN117113855B (zh) | 2023-10-20 | 2023-10-20 | 一种单通道地表温度遥感方法及其系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117113855B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102103203A (zh) * | 2011-01-19 | 2011-06-22 | 环境保护部卫星环境应用中心 | 基于环境一号卫星的地表温度单窗反演方法 |
| CN103902839A (zh) * | 2014-04-18 | 2014-07-02 | 山东科技大学 | 卫星遥感产品支持的单通道法地表温度反演方法 |
| CN114781148A (zh) * | 2022-04-14 | 2022-07-22 | 河北地质大学 | 一种热红外遥感云覆盖像元的地表温度反演方法及系统 |
| WO2022161105A1 (zh) * | 2021-01-27 | 2022-08-04 | 河南理工大学 | 一种多时相多等级的城市温度遥感数据分析方法 |
| CN115031854A (zh) * | 2022-06-22 | 2022-09-09 | 中国科学院空天信息创新研究院 | 热红外高光谱地表温度发射率谱反演方法及装置 |
-
2023
- 2023-10-20 CN CN202311362956.1A patent/CN117113855B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102103203A (zh) * | 2011-01-19 | 2011-06-22 | 环境保护部卫星环境应用中心 | 基于环境一号卫星的地表温度单窗反演方法 |
| CN103902839A (zh) * | 2014-04-18 | 2014-07-02 | 山东科技大学 | 卫星遥感产品支持的单通道法地表温度反演方法 |
| WO2022161105A1 (zh) * | 2021-01-27 | 2022-08-04 | 河南理工大学 | 一种多时相多等级的城市温度遥感数据分析方法 |
| CN114781148A (zh) * | 2022-04-14 | 2022-07-22 | 河北地质大学 | 一种热红外遥感云覆盖像元的地表温度反演方法及系统 |
| CN115031854A (zh) * | 2022-06-22 | 2022-09-09 | 中国科学院空天信息创新研究院 | 热红外高光谱地表温度发射率谱反演方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| 张晓等: "地温敏感波段的遥感图像反演算法研究", 《测绘科学》, vol. 40, no. 5, pages 37 - 43 * |
| 樊辉: "基于Landsat TM热红外波段反演地表温度的算法对比分析", 《遥感应用》, pages 36 - 48 * |
| 赵恩宇等: "中红外双通道夜间数据地表温度反演方法", 《红外与毫米波学报》, vol. 33, no. 3, pages 303 - 310 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117113855B (zh) | 2023-12-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liang | An optimization algorithm for separating land surface temperature and emissivity from multispectral thermal infrared imagery | |
| Ye et al. | Land surface temperature estimate from Chinese Gaofen-5 satellite data using split-window algorithm | |
| CN109297605B (zh) | 一种基于中红外与热红外数据的地表温度反演方法 | |
| Francois et al. | Analytical parameterization of canopy directional emissivity and directional radiance in the thermal infrared. Application on the retrieval of soil and foliage temperatures using two directional measurements | |
| Czajkowski et al. | Thermal remote sensing of near surface environmental variables: Application over the Oklahoma Mesonet | |
| Niclos et al. | Accuracy assessment of land surface temperature retrievals from MSG2-SEVIRI data | |
| CN106778516B (zh) | 一种利用中国风云三号卫星遥感数据确定地表温度的方法及装置 | |
| CN1882854B (zh) | 无线电探空仪的湿度测量结果的校正 | |
| CN103438900B (zh) | 三线阵相机影像协同绝对辐射定标和校正方法 | |
| US10064017B2 (en) | Characterizing and using building weather profile parameters | |
| CN101320475B (zh) | 复杂背景条件下红外成像系统作用距离估算方法 | |
| CN110866467A (zh) | 一种航空中红外高光谱数据温度和发射率反演方法 | |
| CN111737912A (zh) | 一种基于深度神经网络的mwhts模拟亮温计算方法 | |
| CN107368617B (zh) | 基于Lowtran7大气软件的地空探测红外成像系统作用距离计算方法 | |
| CN111323129A (zh) | 一种基于宽波段热红外影像的地表温度反演方法 | |
| CN114544003B (zh) | 一种地表温度遥感产品不确定度的逐像元估算方法 | |
| Ye et al. | Cross-calibration of Chinese Gaofen-5 thermal infrared images and its improvement on land surface temperature retrieval | |
| CN102901563B (zh) | 一种同时确定地表窄波段和宽波段比辐射率的方法及装置 | |
| CN117113855B (zh) | 一种单通道地表温度遥感方法及其系统 | |
| CN113532652A (zh) | 一种基于浮标和大气再分析数据的红外遥感传感器绝对定标方法 | |
| CN115859663B (zh) | 一种考虑不确定度的地表温度反演精度评估方法 | |
| CN115859211A (zh) | 一种基于三温不确定度估算模型的地表温度产品融合方法 | |
| CN114623931A (zh) | 一种扣除背景辐射影响的温度与发射率测量方法及装置 | |
| Lukashevich et al. | LST quality evaluation service for heterogeneous Earth observation data | |
| CN207097004U (zh) | 一种人体特征红外消防探测仪 |
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 |