《GB_T 43997.1-2024 地表温度热红外遥感反演 第1部分：单通道法.docx》由会员分享，可在线阅读，更多相关《GB_T 43997.1-2024 地表温度热红外遥感反演 第1部分：单通道法.docx（21页珍藏版）》请在课桌文档上搜索。
1、ICS35.240.70CCSA77中华人民共和家标准3GB43997.12024地表温度热红外遥感反演第1部分：单通道法Surfacetemperatureretrieva1.fromtherma1.infraredremotesensingdata一PartIzSingIechanne1.method2024-04-25发布2024-11-01实施国家市场监督管理总局*国家标准化管理委员会发布目次前言III引吉IN1范阚12规范性引用文件I3术语和定义14方法原理与适用条件24.1 方法原理24.2 适用条件35反演流程和步骡35.1 操作流程35.2 反演步骤46不确定度分析6附录A（规
2、范性）通道辅光度与温度的关系7A1.通道辐亮度与海度的出换方法7A.2通道辐光度时温度一阶儡导的计算方法7附录B（规范性）多项式法系数确定方法9附录C（规范性）不确定度评定方法10CA地表淑度的标准合成不确定度评估IOC.2反演模鞭分骑自身引起的地表温度不确定度评估10C.3大气参数分埴引起的地表讯度不确定度评估IOC.4地表发时率分僦引起的地表温度不确定度评估13C.5入臃箱亮度分曜引起的地表温度不确定度评估13C.6大气气溶胶分设引起的地上温度不确定度评估14参考文献15本文件按照GBT1.1-2020E标准化工作导W1.第1部分：标准化文件的结构和起草规则3的规定起草.本文件是GB439
3、97E地表温度热红外遥感反演的第1部分。GBzT43997己经发布了以下部分：第1部分：单通道法；第2部分：分裂窗法.请注意本文件的某些内容可能涉及专利，本文件的发布机构不承担识别专利的说任。本文件由中国科学院提出。本文件由全国遥能技术标准化技术委员会(SACyrC327)归F1.本文件起我单位：4国农业科学院农业资源与农业区划研究所、电件供大学、昆明理工大学、苏州中科大启遥感科技有限公司、中国科学院水利部成都山地灾害与环境研究所、河北地质大学、桂林理工大学、中国科学院地理科学与资源研究所、中国科学院大学、北京师范大学、北京大学、中国科学院空天信息创新研究院、大连海事大学、西南大学、山西大学、
4、国家卫星气象中心、中国资源卫星应用中心、中科星图股份有限公司。本文件主要起草人：李召良、段四波、吴弊、唐4曲、周成虎、冷佩、赵伟、尚国排、陈银山、唐荣帆姜小光、脩丁戏、任华忠、刘向阳、马灵玲、钱永刚、邱实、高懋芳、张筱、草志豪、赵恩宇、干文凭、僦缎、黄成、张育羽、之照吉、Em巾硼、欧阳麻、葩福龙、帏启金、蔡文文。GRT439970地表阻度热红外遥超反演3旨在给出利用航天或者航空热红外遥感数据进行地表温度反演的方法，规范流程和步骤,评定反演方法的不确定度.GHT43997拟由四个部分构成.第1部分：单通道法。目的在于明确利用单个热红外通道进行地表温度遥感反演的流程、步骤以及不确定度。第2部分：分
5、裂窗法。目的在于明确利用两个热红外通道进行地表温健遥矮反演的流程、步骤以及不确定度。第3部分：TES分离法,目的在于明确利用三个及三个以上热红外通道迸行地表温I度和发射率协同反演的流程、步骤以及不确定度。-第I部分：日夜法.目的在于明确利用日夜双时相多个热红外通道进行地表温度遥建反演的流程、步界以及不确定度。地表温度热红外遥感反演第1部分：单通道法1范BI本文件描述了利用热红外遥感数据反演地表温度的单通道法原理与适用条件、反演流程和步骤以及不确定度分析.本文件适用于在热红外大气沟口区内采用一个通道数据的他表温度反演，2H范性引用文件本文件，GBTGBTGBT415342022415382022
6、415412022下列文件中的内容通过文中的规莅性引用而构成本文件必不可少的条款。其中，注I期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于地表温度遥感产品真实性检验地衣发射率遥感产品JX实性检龄热红外遥感超本术语GB.T415342022.GB，T415382022.GB门415412022界定的以及下列术i?5和定义适用于本文件.3.1地球表面温度surfacetemperature表征地球表画厚度等于穿透深度（范围0.110倍波长）的表皮的粽合甜度.注：逸盛数据反演得到的地表湘度是在性修器获取的硒温度的伽比消除了大气和发射率彰响后的地
7、农非向温混合像元的等效油以即方向辎射“五度）,包括征地表面温度和海洋表面温度，的位为开尔文（2，【来源：GBT415342022.3,23.2光谱】发射率Ispectni1.1.em1.ssivityDtift比辐射率热辐射体的自身光谱J辐射亮度与同温度同波长下的黑体的光谱相射亮度之比.来源：GBT4153820223.1）33J*D*f1.M%ftAtseosarapectra1.radi三jc星上光阳辐射亮度sate1.1.ite（spettra1.）radia11ce在给定方向上，表面离地光诙相射）亮度经大气作用后被P星传感器接收到的光i朴福射亮度.来源：GBrr415412022,3.
8、261Xi普朗克定律P1.ancks1.aw冷朗克函数P1.ancksfunction黑体S射定律b1.ackbodyradiation1.aw描述给定温度下处于热平衡状态的黑体的光谱招射亮度的波长或波数变化的函数关系的定律。注：随波长纪匕的普朗克定律表达W见公式-(J9Ig3,单位为焦每开尔文(JK);T奥体温度.单位为开尔文(K).Iif1.波数变化的普朗克定律衣达式见公式.bII7.:r1(2)*,-1.式中：Be奥体ff尤渊矶射比度,单位为比短平方米球面度负次方厘米W*nRwcmJH;h普朗克常数,取位为6.62607015XIcH.取位为焦杪(J。C-JtiS.取值为Z9979245
9、X10.单位为米每杪(mXy一波数.取位为标厘米(CmkB玻尔兹及常数.取值为1.380&19X1023,单位为供传开尔文(JKkT-黑体温度，单位为井尔文(K)。来源：GBT415412022.3.10)3.5大气柱总水汽含员tota1.coIuiuiWatervapor地面到大气层顶垂直路径上单位板面积内所含气态水的总质量。注：单位为克林平方匣米但m”3.6近地面空气温度surfaceairIenipera1.iire距地表下桀面高度介于1.5m2m处的空气温度。4方法原理与适用条件41方法原理星于热红外相射佻检模型，单通道法通过抉取的大气参数来消除大气效应的影响，进而实现地衣温度的遥感反
10、演.学洒道法的数学模型按公式(3)衣示:T.(0.)If,(。)一1t(.f)式中1TeM观测天顶角0和观测方位角巾上的地表温度，单位为开尔文(0观测天顶角,单位为度()：观测方位角,电位为度C);BJ一通道i的密朗克定律的逆函数，计律方法按附录A的A.I执行;1.i(0,)观测天顶角。和观测方位角小上通道i的入瞳轴亮度,单位为瓦每平方米球面度微XW(mjsrm)或乩每平方米球面度负一次方厘米W(m*srcmJ”;c,(O,)观测天顶角0和观测方位角上通道i的地表发射率：Ia一通道i的大气下行幅射.单位为瓦每平方米球面度微米IWMm37h砌或瓦每平方米球面度负一次方闻米W(nRsrcnM)J:
11、IXo)一观测天顶角0上通道i的大气透过率：1.(O)一观测天顶角。上通道i的大气上行和射，单位为瓦林平方米球面度微米W(m2srm)或f1.3平方米球面度负一次方厘米IW(m*-srcm-)b注IUQg)=B1.T1.勒刈，其中，B.为通道i的普朗克定律；Tj(Oe)为观测天顶角。和观测方位角小上通道i的XJI也b单位为开尔兑Id4.2采用单通道法反演地表温fS.应满足以下条件：a)特反演像元所对应的传感器瞬时视场向内无云：b)待反演像元沙尘气溶胶的光学厚度小于0.3：c)传法器在8Pm14m波长范围内至少仃一个通道：d)待反演像元的通道地表发射率己知；C)传感器观测时刻前后Ih内的大气温湿
12、晚线已知，或大气柱总水汽含量和近地面空气温度已知,或大气柱总水汽含量已知.55.1地表温度热红外遥媾反演单通道法的操作流程见图10图1地表温度热红外遥感反演单通道法流程图5.2反演步赛1 .2.1反演通道选择如果传感器的热纤.外通道数盘大于1,陶选取大气吸收最小(即大气透过率最大)且地衣发射率较大的通道诳行地表温度反演，5 .2.2大气叁数估算方法选择根据大气信息获取的情况，选择合适的大气参数估算方法。a)当与传感器观泅时刻同步(时间相差前后1h内)的大气讯泡脚线已知时,应优先采用保射传输模拟法,通过将同步大气温湿施线、传感器光诣响应函数、观测角度输入到大气辎射传输模型中，计算得到大气透过率大
13、气上行辐射1“t.-()-6wu()TW”AJ=JTW式中：t；(0观测天顶向D上通遒i的人气透过率；0观测天顶角，单位为度()；a三n.；(0)一一人气透过率与大气柱总水汽含量和近地面空气温度之间的二元二次多项式拟合系数，是观测大顶角0和通道i的函数：T,近地面空气制弘单位为开尔文(K);W大气柱总水汽含S1.单位为克每平方厘米(gcm2);1.(0)-观测天顶角。上通道i的大气上行辐射,单位为瓦每平方米球面度微米(W/(mjsrm)或瓦每平方米球面度负一次方瓯米W(m2srcm-,:b三:(0)大气上行辎射与大气柱总水汽含M和近地面空气温度之间的：元：次多项式拟合系数，是观测天顶角0和通遒
14、i的函数；1.a,一通道i的大气下行辐射，单位为瓦每平方米球面度做米W(mrm)或瓦每平方米球面度负一次方厘米(WXnAsrsP);C皿一一大气下行幅射与大气柱总水汽含量和近地面空气淑度之间的二元二次多项式拟合系数,是通道i的函数.C)当未知同步大气温湿成戏，且仅已知大气柱总水汽含Jft时，采用一元二次多项式法见公式(7)公式(9)进行估算.30)=2;(O)W2+a;(0)W+ao;(0)1.u:(。尸2:(O)W%(Q)W+o0)(8)1.a=Y2B+y1W+Yo;式中；U(O)观测天顶角。上通道i的人气透过率：0天顶角，单位为度C)：g(0)火气透过率与大气柱总水汽含量的一元二次多项式的
15、二次项拟合系数，是观测天顶角0和通道i的函数：w大气柱总水汽含Bt单位为克每平方厘米(WCm4：1：(0)大气透过率与大气柱总水汽含量的一元二次多项式的一次项拟合系数，是观测天顶角0和通道i的函数：ao;(0)大气透过率与大气柱总水汽含碳的一元二次多J式的常数JJj拟合系数.是观测天顶角0和通道i的函数：1.U0)-观测大顶角。上通道i的大气上行辐射,单位为阮每平方米球面度微米W/(mjsrm)或瓦每平方米球面度负一次方J里米W(mXrcm-);B?:。)大气上行辐射与大气柱总水汽含麻的一元二次多项式的二次项拟合系数.是观测天顶角0和通道i的函数：1：(0)一一大气上行辎射与大气柱总水汽含法的
16、元.次多项式的一次项拟合系数，是观测天顶角0和通道i的函数：Bo;(。)大气上行幅射与大气柱总水汽含量的一元二次多项式的常数项拟合系数，是观测天顶角0和通道i的函数；1.d.通遒i的大气下行辐射,单位为瓦林平方米球面度微米W(msrm)或瓦好平方米球面度负一次方显米W(nRsrtmJ):Y2；一一大气下行幅射与大气柱总水汽含址的一元二次多项式的二次项拟合系数.是通道i的函数：Yi大气下行辐射与大气柱总水汽含地的一元二次多项式的一次项拟合系数，是通道i的函数：Yo:-大气下行幅射与大气柱总水汽含量的一元二次多项式的常数项拟合系数是通道i的函数：二元二次多项式法和一元二次多项式法的系数确定应符合冏
17、录B5.2.3大气麴数逐像元估算根据选定的大气参数估算方法,逐像元估算给定观测大顶角条件下的大气参数.U)当同步大气甜湿廓找己知时,根据传呼器观测天顶地的变化范用0,OmxJ.在观测天以角正割俏sec(e)的变化区间U,M0mx)J内，按照等步长(例如0.1)设置观测天顶角，再将大气温湿战跳、传感器光谱响应函数、所设置的观测天顶角输入到大气粕射传输模蟹中，计算得到该大气状况下不同观测天顶角对应的大气透过率、大气上行辐射以及与观测角度无关的大气下行辐射；然后，根据待反演像元的实际观测天顶角，按照观测天顶角的正割伯sec(6对大气透过率、大气上行辐射分别进行线性场值,逐像元分别估算大气透过率、大气
18、上行将射，b)当未知同步大气ifi湿廓鼓，但大气柱总水汽含量和近地面空气温度己知时，根据恃反演像元的观测天顶角，按照观测天顶角的正例值SeC(O)对公式(4)和公式(5)的拟合系数分别进行线性排值，获得像元时应观测角度卜公式(力和公式的拟合系数，结合公式的拟合系数，利用公式公式(6)、大气柱总水汽含量和近地面空气温度，逐像元分别估算大气通过率、大气上/下行粗时.c)当未知同步大气湿湿廓线，且仅已知人气柱总水汽含里时，根据待反演像元的观测天顶角，按照观测天顶角的正割值Scc(O)对公式(7)和公式(8)的拟合系数分别进行线性辅值,获得像元对应观测角度下公式(7)和公式(8)的拟合系数，结合公式(
19、9)的拟合系数，利用公式(7)公式(9)、大气柱总水汽含域，逐像元分别估算大气透过率、大气上/下行辎射,5.2.4地表温度反演将入腌辅光度、地表发射率和5.2.3估算的大气透过率、大气上/下行辎射输入到单通道法的数学模型公式（3）中，逐像元反演得到地表器度，6不定度分析不确定度分析应逐像元分析,主要步骤包括：a）分析影响单通道遥嬷反演地表温度的不确定度来源：b）根据不确定度传播律，计算地表温度的合成标准不确定度，合成标准不确定度计匏方法按附录C的C1.执行：c）估算每个不确定度分做引起的地衣温度不确定度，各分法的不确定度计算方法按C2C6执行.附录A(规范性)通道辐亮度与温度的关系A.I通道锚
20、亮度与温度的转换方法在一定的温度范国内(例如240K330K),温度按照一定的步长(例如0.1K)变化,通过普朗克定律与传感器光i?j响应函数卷积计修得到每个温度对应的通道福亮度，见公式(A.1.):Bj(T)iQ)UB.(T)(AJ)J/.Q)U式中：B:(T)通道i的思体播亮度,单位为瓦每平方米球面度微米M(mrm”:T温度,单位为开尔文(K);,通道i的光谱响应函数的最小波长.舱位为微米(Um)：?一通道i的光谱响应函数的最大波长，单位为微米(m);Ba(T)一温度为T的黑体尤诺福亮度，单位为瓦年平方米球面度微米W(nPs11m3:f;C)-通道i的光谱响的函数.利用公式3)建立通道箱亮
21、度与温度的查找表，然后通过战性插值进行通道锚亮度与温度的相互转换，A.2通道辐亮度对温度一阶偏导的计算方法通道幅亮度对度一阶偏导可利用公式(a.2)进行计算：一心(A2)7f-式中：B(T)通道i的黑体辑亮度B.(T)对温度T的一阶偏3,单位为瓦每平方米球面度微米开尔文Wf(m1srmK);T温度，单位为开尔文(K);一一通道i的光谱响应函数的嫉大波长，单位为微米(Uin);：一一通道i的光谱响应函数的皎小波长，单位为微米(“m);B2(T)温度为T的黑体光i普辅亮度，单位为瓦每平方米球面度微米W(m3srm)J;C2第二辐射常数，单位为1.4的XIOaPmKf1.()通JHi的光谱响应函数.
22、在实除使用中常用下列简单的近似公式(A.3)计算:-HI(T)H(T+7)11f式中：B,(T)通道i的黑体辐亮质B;(T)对温度T的一阶偏导，单位为乩每平方米球面度微米开尔文(W(m2-srmK);T温度，单位为开尔文(K);B1(T)温度为T的通道i的黑体辎亮度,单位为瓦林平方米球面度微米WNsrm),按照A.1的步骤计算；T温度变化出附录B多贩法系数础防法二元二次多项式法公式(力公式(6)和一元二次多项式法公式(7)公式(9)的系数确定操作流程见图B.I.具体步骤如卜Za)大气温湿麻线选取，根据待反演区域的大气特性，从现有的大气温湿麻线数据库或待反演区域的大气温湿除级历史数据中,按照大气
23、柱总水汽含量均匀分布和大气底层空气温度均匀分布的原则，选取一定数6(例如，不少于60条)且涵茂符反演区域晒空大气温湿度变化范掴的大气温湿廊线。b)观测天顶角设留根据传嬷器观测天顶角的变化范国0,OmaX1.,在观测天顶角正割伯SeC(O)的变化区间1.scmax)内,按照等步长(例如QD设四观测天顶角.C)数据模拟。将步骤a)选取的大气温湿廊线和步骤b)设置的观测天顶角以及传覆器光谱响应函数输入到大气辐射传输模受中，运行模型得到大气透过率、大气上厅行箱射。d)大气柱总水汽含故与近地面空气温度获取.根据步骤a)选取的大气温湿廊线计算大气柱总水汽含量,并提取大气阻湿廓线的最底层温度作为近地面空气温
24、度。e)多项式法系数拟合。将步骤C)获得的大气透过率、大气上/下行辐射和步骤获得的大气柱总水汽含G与近地面空气温度输入到公式(4)公式(6),或将步骤C)投御的大气透过率、大气上下行辐射和步骤d)获得的大气柱总水汽含量输入到公式公式(9),进行多项式拟合确定多项式法系数，ff1.H.1多理式法系败流程网附录。(规楚性)不确定度评定方法C.1地表温度的标准合成不确定度评估单通道法反演地表温度的数学模型按照公式(3)执行。根据不确定度传播律，假设地表M度不确定度来源之间不相关，地表温度标准合成不确定度可按公式9I)计匏：uwu=(ug)2(u)2+(u)2+(u)1+(u)2(C.1)式中：UOr
25、a1.地表温度的合成标准不确定度：Ua1g一一反演模型分属自身引起的地衣温度不确定度：Uan1大气多数分员引起的地表温度不随定度：UCm一地衣发射率分最弓I起的地表温度不确定度：urd入瞳辎亮度分房(包括仪器噪声和辎射定标不确定度)引起的地衣温度不确定度：uero一大气气溶胶分砍引起的他去温哎不确定度，C.2反演模型分量自身引起的地表温度不确定度评估在使用公式(3)反演地表制度过程中，反演模型自身在在不确定度，反演模型分成自身引起的地表温度不确定度按公式(C2)计算：式中：Ua1.g反演模型分量引起的不确定度：HHfiifii的黑体辐亮度对地表显度一阶体导的倒数，可利用A.2计算：1.：通道i
26、的入瞳箱亮度：1.u通道i的大气上行相射：E;通遒i的地表发射率：T;ifiifii的大气透过率；,通道i的大气透过率估算模型的不确定度：u一一通道i的大气上行辐射估算模型的不确定度：u1.-一通Uii的人气下行辅射估算模型的不确定度.注：当选择辐射传输模拟法进行大气参敢估帆，u.U1.、UI分别为大，网射传输模型估悯遒i的人气透过率，大气上行箱射，大气下行秘幡J模里不确定度：当选择二元二次多项式法或一元二次多项式法进行大气参数估算时，1.为公式或公式的模型不确定度，U1.为公式或公式的模里不确定度，u为公式或公式的模型彳浦定度.C.3大气参数分引起的地表温度不辞定度评估大气参数分量引起的地衣
27、温度不确定度计算方法如下.a)当选择辐射传输模拟法进行大气参数估算时,根据5.2.2所使用的大气淑湿怫线的不确定度，在地表温度反演过程中，依据先蛤知识，将所使用的大气温湿靡跳各层加入给定的特定分IO布的茶龙,通过蒙特R罗法评定由大气湿湿凯峻分峡中起的最表温度不确定度，见公式(C3);Uam=E(g(x)-g(xo.y)Fx.y1.2(C3)式中：Uam大气参数分量引起的地表淞度不确定度IEU-算术平均值：g单通道法反演地表湿度的数学模型；X-加入不施定度的大气塌湿而绥：y入瞌辑亮度和通道发射率的数据集：X没有加入不确定度的大气温湿睨淡.b)当选择二元二次多项式法进行大气参数估算时,根据公式(3
28、)公式(6),大气参数分量大气柱总水汽含量不确定度(UW)和近地面空气温度不确定度(ur。)引起的地表温度不确定度表示为公式(C.4)一公式(C.6):m-7(m*,其中,1.t-1.,1.9.11礼.“,(1i)J1.,1.rfMJ1.im无一)B,(T.)1.1通道i的入Itt辐亮度：1.w通道i的大气上行储射：E通道i的地衣发射率：t通道i的火气透过率；一大气透过率对大气柱总水汽含量的阶儡导,按照公式(C7)if-W：一大气上.行相射对大气柱总水汽含砥的一阶倡导.按照公式(C.8)计算：一一大气卜行辎射时大气柱总水汽含量的阶偏导,按照公式9)计算：一一大气透过率对近地面空气温度的一阶编导
29、.按照公式(UI0)计完：一一大气上行辐射对近地面空气温度的阶佰导，按照公式(C.II)计犯：大气下行解射时近地面空气湿度的一阶偏导，按照公式(C.12)计算。大气透过率、大气上行储射和大气下.行辎射对大气柱总水汽含Ift的-阶偏导计算见公式(C7一公式(C9):Jr而/-(C.7)(C.8)(C.9)式机anun.:(0)公式(4)的拟合系数.Ta一近地面空气温度：W一大柱总水汽含量：bmm,;(0)公式的拟合系数,Cmn.一公式(6)的拟合系数.是观测天顶角和通道i的函数:是观测天及角和通道i的函数:是通道i的函数,大气透过率、大气上行辐射和大气下行辐射对近地面空气温度的一阶偏导计算见公式
30、(C.10)公式(CJ2):夺-七2i，Tr俨x-w*-t.,.,I(I-a5w+F7(C.13)式中：Uatm1.1.mEIW.IW一一大气参数分址引起的地上温度不确定度:通道i的黑体辐比度时地表温度一阶偏壮的倒数，可利用A.2计算;通道i的入瞳辐亮度：通道i的大气上行辅射;一一通道i的地表发射率：通道i的大气透过率：一大气透过率对大气柱总水汽含量的一阶偏导,按照公式(C.1率计算;-大气上行羯射对大气柱总水汽含Ift的一阶偏导,按限公式915)计驾:大气下行辎射时大气柱总水汽含盘的阶偏导,按照公式(CI6)计宛:g一大气柱总水汽含量不确定度.大气透过率、大气上行辎射和大气下行辐射对大气柱总
31、水汽含量的一阶偏导计算见公(C.I4)-公式(CI6):；w=21.9W-a(0)(C.14)Uc-15)式中：2;(0)-公式(7)的二次项拟合系数.是观测天顶加0和通道i的函数:W一一大气柱总水汽含量：.：(0)一一公式。)的一次项拟合系数,是观测大顶角0和通道i的函数:Pj()公式的二次项拟合系数，是观测天顶角和通道i的函数;Y2i公式(9)的二次项拟合系数，是道道i的函数：Y一一公式(9)的一次项拟合系数，是通道i的函数.C.I地表发射率分量引起的地表温度不确定度评估地我温度反演模型公式(3)是通道i的发射率的显式函数,根据公式3,地衣发射率分后引起的地表温度不确定度表示为公式(CI7
32、:.T/de(C.17)-JHfT.)%就机Ucms胞表发射率分瞅引起的地表温度不确定度：、一通道i的黑体辐亮度对地表温度一阶偏导的倒数,可利用A2计算：D,fa)B1(T.)温度为T.的通道i照体根亮度.单位为瓦每平方米球而度微米W(mrm);1.u一一通ifii的大气下行辎射：巳通道i的她去发射率：U.通igi的地表发射率不确定度.C.5人睁镉亮度分量引起的地表温度不确定度评估根据公式(3,入瞳辐亮度分量引起的地表温度不确定度衣示为公式(CI8):-1G)式中：Urad入瞭福亮度分奴(包括仪器噪向和辐射定标不确定度)引起的地表显度不确定度：n；：r)通道i的黑体福亮度对地表温度一阶偏导的倒
33、数，可利刖A2计算：E;通道i的地表发射率；t通道i的大气透过率：u1.通道i的入瞳辐亮度不确定度.ISC.6大气气溶胶分量引起的地表温度不确定度评估地表温度反演模型公式(3”是大气气溶胶含地的隐式函数，根据公式(3)分别反演考虑和不考虑大气气溶胶类型与含量时的地表湿度，通过蒙特卡罗法评定由大气气溶胶分矍引起的地表温度不确定度，表示为公式(C1.9):ur=E(Ts-Ts)2立(C.I9)式中：gem大气气溶胶分量引起的地表温度不确定度：EH算术平均值：T.考虑实际大气气溶胶类型和含变化时反潮的地衣温度；Tsdcau1.t仅考虑缺省的大气气溶胶类型和含Ift时反演的地表温度,考文献1GRT36
34、299-2018光学遥感辐射传输基本术语(2UJF1059.12012测量不确定度评定与表示3ISOIECGuide98-3:2008Uncertain1.yofmeasurementPart3:GUideioIheexpressionofuncertain1.yin11eaurc11ent(GUM:1995)4Crisioba1.J.JiinGnez-MunozJ.C.SobrinoJ.A.,ea1.Inprovenentsin1.andsurfacetemperatureretrieva1.fromthe1.andsatseriestherma1.bandusingwatervaporan
35、dairtemperature(J).Journa1.ofGeophysica1.Research-Atmosphcrcs.2(X)9,114,D08103.5Jin泾nezMunozJCandSobrinoJAAgenera1.izedsing1.e-channe1.methodforretrieving1.andsurfacetemperaturefromremotesensingdataJJourna1.ofGeophysica1.Rcseanrh-Atmospheres.2(X)3,108(D22).4688.61.iZ.-1.agB.H.WuH.,eta1.Sate1.1.ite-d
36、erived1.andsurfacetemperature:CUfTentstatusandperspectivesJj.RemoteSensingofEnviro11iient.2()13,131.14-37.71.iZ.-1.WuH.WangN.QiuS.eta1.1.andsurfaceemissivityretrieva1.fromsate1.1.itedataIJjJnternationaIJourna1.ofRemoteSensing.2013.34(9-10),3084-3127.(8QinZHXarnicIiAandBer1.inerP.Amono-windowa1.gorithmforretrieving1.andsurfacetemperaturefn)m1.andsatTMdataanditsapp1.icationtotheIsrae1.-EgyptborderregionJ.Internationa1.Journa1.ofRemoteSensing,2(X)1.223719-3746.
