卫星遥感地表温度降尺度的光谱归一化指数法

卫星遥感地表温度降尺度的光谱归一化指数法李小军;辛晓洲;江涛;张海龙【摘要】Aimingattheproblemthatthespatialandtemporalresolutionoflandsurfacetemperature(LST)havethecontradictionwitheachother,anewdownscalingmodelwasputforward,basedontheTsHARP(analgorithmforsharpeningthermalimagery)downscalingmethod,thisresearchmakesimprovementsbyselectingthebettercorrelationofspectralindex(normalizeddifferencevegetationindex,NDVI;normalizeddifferencebuild-upindex,NDBI;modifiednormalizeddifferencewaterindex,MNDWI;enhancedbaresoilindex,EBSI)withLST,i.e.,replacestheoriginalNDVIwithnewspectralindexaccordingtothedifferentsurfaceland-covertypes,toassesstheaccuracyofeachdownscalingmethodbasedonqualitativeandquantitativeanalysiswithsynchronousLandsat8TIRSLSTdata.Theresultsshowthatbothmodelscouldeffectivelyenhancethespatialresolutionwhilesimultaneouslypreservingthecharacteristicsandspatialdistributionoftheoriginal1kmMODISLSTimage,andalsoeliminatethe"mosaic"effectintheoriginal1kmimage,bothmodelswereprovedtobeeffectiveandapplicableinourstudyarea;globalscaleanalysisshowsthatthenewmodel(RMSE:1.635°C)isbetterthantheTsHARPmethod(RMSE:2.736C)intermsofthespatialvariabilityandaccuracyoftheresults;thedifferentland-covertypesofdownscalingstatisticalanalysisshowsthattheTsHARPmethodhaspoordownscalingresultsinthelowvegetationcoveragearea,especiallyforthebarelandandbuilding-uparea(|MBE|>3°C),thenewmodelhasobviousadvantagesinthedescriptionofthelowvegetationcoveragearea.Seasonalanalysisshowsthatthedownscalingresultsoftwomodelsinsummerandautumnaresuperiortothoseinspringandwinter,thenewmodeldownscalingresultsarebetterthantheTsHARPmethodinthefourseasons,inwhichthespringandwinterdownscalingimprovementisbetterthansummerandautumn.%针对卫星遥感技术监测地表温度(landsurfacetemperature,LST)存在时空分辨率矛盾这一难题，以TsHARP温度降尺度算法为基础，根据地表覆盖类型的不同,分别选择与LST相关性更好的光谱指数(归一化植被指数,NDVI；归一化建造指数,NDBI；改进的归一化水体指数,MNDWI；增强型裸土指数,EBSI)提出了新的转换模型，并从定性和定量两个角度评价了TsHARP法和新模型的降尺度精度结果表明:两种模型在提高LST空间分辨率的同时又能较好地保持MODISLST影像热特征的空间分布格局，消除了原始1km影像中的马赛克效应，两种模型均能够达到较好的降尺度效果;全局尺度分析表明,不管是在降尺度结果的空间变异性还是精度方面，本文提出的模型(RMSE:1.635C)均要优于TsHARP法(RMSE:2.736C);TsHARP法在水体、裸地和建筑用地这些低植被覆盖区表现出较差的降尺度结果,尤其对于裸地和建筑用地更为明显(|MBE|>3C),新模型提高了低植被覆盖区地物的降尺度精度;不同季节的降尺度结果表明，两种模型都是夏、秋季的降尺度结果优于春、冬季，新模型的降尺度结果四季均好于TsHARP法，其中春、冬季的降尺度精度提升效果要优于夏、秋季.【期刊名称】《测绘学报》【年(卷)，期】2017(046)003【总页数】9页(P353-361)【关键词】MODIS；降尺度;地表温度;TsHARP算法;地表覆盖【作者】李小军;辛晓洲;江涛;张海龙【作者单位】中国科学院遥感与数字地球研究所遥感科学国家重点实验室，北京100101;中国科学院大学，北京100049;中国科学院遥感与数字地球研究所遥感科学国家重点实验室，北京100101;山东科技大学测绘科学与工程学院,山东青岛266590;中国科学院遥感与数字地球研究所遥感科学国家重点实验室，北京100101【正文语种】中文【中图分类】P237地表温度是地-气系统相互作用过程的重要物理参数，在地表蒸散量估算［1］、土壤水分估计［2］及城市热岛［3］等研究领域中应用广泛。但是，受到成像条件的制约，目前常用的卫星热红外传感器均存在时空分辨率之间的矛盾［4］，或时间分辨率较高而空间分辨率较低（如MODIS），或时间分辨率较低而空间分辨率较高（如LandsatTM）。这种矛盾使得LST数据的使用受到极大限制［5］，为获取满足应用需求的LST产品，地表温度降尺度技术应运而生。目前LST降尺度方法有很多，基于像元尺度的LST尺度转换方法主要分为3类：数理统计方法、光谱混合模型法和调制分配方法。数理统计回归法的基础假设是〃关系尺度不变性”，即LST与回归核的统计关系在各个尺度上保持不变。回归核包括反照率、组分权重和光谱指数等［1,5-6］。光谱混合模型的基本原理是线性光谱混合模型，通过关联不同尺度上的LST，进而回归求解低尺度空间（即高分辨）上的LST［7］。基于调制分配的方法是将高尺度空间的LST按照一定的权重分配给子像元，分配因子包括发射率、同一时刻获取的更高分辨率传感器的LST以及二者的组合等［5］。按照LST降尺度方法中涉及的尺度因子个数，降尺度方法分为单一尺度因子法［1,6,8］和多尺度因子法［9-12］。多尺度因子法主要建立LST与多个相关因子间的线性关系，如逐步回归法［9］和主成分分析［10］，或更为复杂的关系，如支持向量机［11］或神经网络［12］。单一尺度因子法因其物理意义明确并且模型简单而得到广泛应用。文献［1］首次提出利用LST与NDVI(normalizeddifferencevegetationindex)的关系来实现LST降尺度的DisTrad(disaggregationprocedureforradiometricsurfacetemperature)算法，通过建立两者间的二次回归模型，成功实现对千米级到百米级LST的降尺度。文献［6］对DisTrad算法进行了改进，通过对LST、植被覆盖度和植被指数进行回归分析的基础上提出了TsHARP(analgorithmforsharpeningthermalimagery)算法，该算法认为LST与植被覆盖度之间的一元线性关系能够达到更好的降尺度效果。文献［8］利用TsHARP法的5种变形讨论了全局、分段回归、局域和土地利用分层的降尺度模型，并认为分层回归降尺度法的精度要高于原始全局方法。但是该方法中涉及的尺度因子均为NDVI，对于下垫面复杂的区域，不同类型地物的LST最佳指示因子并不一定是NDVI指数。上述研究已经取得突破性的进展，但仍存在一些问题需要解决：①多数研究是先将高分辨率的LST数据聚合至低分辨率而后进行降尺度处理，并再利用原始LST对降尺度结果进行验证，原始LST数据在重采样或升尺度过程中会产生新的误差，会对降尺度结果产生大的干扰，因此其结果缺乏说服力。②NDVI虽然被人们最早用于LST降尺度研究，但是对于下垫面复杂的城市地区，不能较好地反映建筑物和裸地的高温区以及水体造成的低温区。针对上述问题，本文以TsHARP降尺度方法为基础，直接对原始MODISLST产品进行降尺度研究，根据地表覆盖类型的不同，选择与LST相关性更好的光谱指数建立二者之间的关系实现降尺度来弥补问题2的不足，由于Landsat8热红外波段的分辨率明显高于MODIS数据，并且文献［13—14］提出的单通道算法进行Landsat8热红外传感器TIRS(thermalinfraredsensor)影像地表温度反演具有很好的精度和敏感性，TIRS和MODIS成像时间也非常接近(本文选择的数据二者相差小于35min)，因此可以把TIRS反演的地表温度作为当日地表温度的参考值，对不同方法的降尺度结果进行验证。研究区包含北京北以及河北部分地区，如图1所示，其地理位置为116°6'15”E一117°45'12”E,40°14'52”N—40°59'19”N,属于暖温带大陆性气候；区内地势平坦，只有西北部和东部有部分山脉，东南部为平原，属于华北平原的西北边缘区。中下部为北京郊区，包括怀柔城区等下垫面为不透水层覆盖地区以及水系比较发达的密云区；周边主要以林地和耕地覆盖，整个研究区地表覆盖类型复杂，空间异质性高，对进行地表温度降尺度研究具有代表性。本文使用的数据主要有MODIS传感器、Landsat8上搭载的多光谱陆地成像仪OLI(operationallandimager)和TIRS3种传感器的数据，数据列表见表1。研究采用对应区域(Path:123;Row:32)的Landsat8影像，成像时间为2014-05-15、2014-8-19、2014-10-06、2014-12-25，成像效果好，平均云量均低于6%，研究区内无云量，利于进行本文的研究。MODIS数据来自于美国NASAEOSDataGateway提供的MODIS产品，包括4个时相覆盖研究区(h26v04)的500m分辨率反射率产品(MOD09GA)、1km分辨率的每日地表温度产品(MOD11A1)和大气水汽产品(MOD05_L2),数据均可在NASA官网(https:/)免费下载。获取的Landsat8OLI_TIRS已经经过几何校正并转至WGS-84/UTM投影下，故只需对其进行辐射定标、大气校正，将OLI多光谱波段DN值转为地表反射率值，使用文献［13］提出的单通道算法来反演地表真实温度，该方法所需参数较少，并且精度和敏感性较高［14-15］。MODIS产品原始的投影为正弦投影，本文使用NASA研发的MODISSwathReprojectionTool(MRTSwath)和MODISReprojectionTool(MRT)分别对MOD05_L2和另外两种MODIS产品将投影转化为WGS-84/UTM，投影带号为50，最后对上述数据进行裁剪，获得研究区数据。LST降尺度转换总是面临更高分辨率上信息不足的问题，因此常通过引入额外信息构造趋势面才能实现降尺度转换［16］。TsHARP降尺度方法就是通过引入NDVI或植被覆盖度来构造趋势面，其基本假设是：LST和植被指数的关系在各个尺度上基本一致，土壤湿度分布的空间变化是最主要影响因素［17］。以MODIS数据为例，TsHARP方法的原理如下：建立1000m尺度上LST与趋势面因子间的关系，如式(1)所示TsHARP算法建立了多种趋势面转换函数，本文选择其中一种模型进行对比研究(式(4))。为解决问题2，文中利用决策树模型将地物分为植被、建筑用地、水体和裸地，在保证最小残差和最大决定系数的前提下，分别选择NDVI、NDBI(normalizeddifferencebuild-upindex)、MNDWI(modifiednormalizeddifferencewaterindex)和EBSI(enhancedbaresoilindex)来建立不同地物光谱指数与LST之间最佳回归关系式(模型2)从而实现降尺度，4种指数的计算方法分别参考文献［18—21］，技术路线如图2所示，两种模型如下：模型1Ts=a+bNDVI+cNDVI2模型2植被Ts=f(NDVI)—Ts=a+bNDVI+cNDVI2建筑用地Ts=f(NDBI)—Ts=a+bNDBI水体裸地对1000mLST到500mLST的降尺度效果可以从定性和定量两个角度进行。定性分析通过趋势面因子对地物指示性特征分析、对计算结果进行目视对比分析。定量分析包括以下指标。决定系数(R2)相对平均偏差(MBE')相对均方根误差(RMSE')降尺度结果定性评价如图3所示，对比图3中两种模型(图3(c)、(d))的降尺度结果与原始1kmMODISLST(图3(a))可以发现，两种降尺度模型均可以捕捉到大部分原始1kmMODISLST中的地表温度信息，即两幅500mLST降尺度影像较好地保持了原始地表温度影像热特征的空间分布格局，且影像中高温区和低温区均与图3(a)中吻合得较好。此外，两种模型的降尺度结果色调更加丰富(尤其是模型2),更加精细地描绘了地表热特征的空间差异性，较好地消除了图3(a)中的马赛克交攵应。升尺度后的Landsat8TIRS地表温度(500m)如图3(b)所示，其与1kmMODISLST的热特征分布基本相同，并且信息更加丰富，影像热特征分布空间变异性更加明显。对比图3(b)发现，模型2(图3(d))的降尺度结果要比模型1(图3(c))信息更加丰富，同种类型地物温度分布的空间差异性更为明显，温度过渡之间没有明显的过渡圈。为了更加充分验证降尺度的效果，下面对两种模型的降尺度结果进行更进一步的定量统计分析。为了对降尺度结果进行定量评价，除了计算上述评价指标(随机选择5000个点)，还对升尺度Landsat8TIRS500mLST和两种模型降尺度的500mMODISLST进行了统计分析，统计量包括均值(Mean)、标准差(SD)，此外对Landsat8TIRSLST和两种降尺度结果还分别进行了散点图分析，以求更加直观地显示降尺度结果的优良，如图4所示。统计结果见表2，由表可知，两种模型降尺度的500mMODISLST均值均接近于Landsat8TIRS升尺度的LST均值(25.865°C),说明两种模型降尺度的结果与实际验证温度之间的总体偏差较小，其中模型1的平均值偏差较大(1.160C)。SD反映了图像像素值的离散程度，对于热红外影像而言，SD反映了地表热特征的空间变异性［17］。两种模型的SD计算结果表明，模型2的SD最大，并且最接近于TIRS500mLST的SD,说明与模型1相比，模型2能够更好地体现不同下垫面地物的温度差异以及相同下垫面覆盖地物的细节信息，这一点从图3中也可以看出。两种模型的R2差异并不明显，其中模型2的R2最大。RMSE和MBE能够反映降尺度影像和验证影像之间的一致性，可以用来指示两种降尺度模型与TIRSLST热特征分布的差异性。模型1的RMSE要大于模型2,说明模型2很好地再现了TIRS500mLST地表热特征的空间分布格局，能够较好地保持地表真实温度的热特征，而模型1较差，从MBE可以看出，两种模型的降尺度结果与TIRSLST的平均偏差均较小(绝对值小于1C),降尺度结果非常接近真实地表温度，并且两种降尺度结果的MBE均小于0,印证了两种模型降尺度的500mMODISLST均值均低于TIRS500mLST的均值，这一点由散点图(图4)也可以得到证实，综合对上述统计量的讨论，两种算法均能够达到较好的降尺度效果，不管是在降尺度结果的空间变异性还是精度方面，模型2均要优于模型1。本文统计了不同地物覆盖类型的降尺度结果(表3),由表3可知，原始全局降尺度模型(模型1)在水体、裸地和建筑用地这些低植被覆盖区表现出较差的降尺度结果，尤其对于裸地和建筑用地更为明显(|MBE|>3C)。两种降尺度模型均表现出对植被、裸地和建筑用地为低估(MBE<0),而水体为高估(MBE>0)，图4也印证了这一点，对整个研究区而言水体表现为低温，低温区位于1:1直线上方，高温区则位于下方，这是由于原始全局降尺度模型在构建趋势面转换函数时剔除了水体，光谱指数二次项的存在会对水体和建筑用地的降尺度结果产生大的干扰［6］。对比模型2和模型1,在植被地区两种模型的差异并不明显，因为在植被地区模型2与模型1使用了相同形式的趋势面转换函数。而对于其他地物类型，新参数的引进和替换均在一定程度上提高了降尺度结果的精度。对于水体地区，模型2的MBE要比模型1提升1.041°C，这是因为原始全局降尺度模型没有考虑水体，而文中改进的模型2，单独建立水体LST与MNDWI之间的关系，较好地避免了其他地物对水体区域的影响。裸地的MBE由原先的-3.203C提升到了-2.3801，建筑用地则由-3.463C提升到了-2.0141，提升了41.842%，这是由于模型2不仅避开了不同类型地物构建趋势面转换函数时产生的影响，而且对于建筑用地，NDBI与城市中心区的LST拟合效果较好，利用其与LST之间的线性关系对于城市地表温度的降尺度效果要优于单独使用NDVI［16］。两种模型降尺度的500mMODISLST与TIRS500mLST差异的空间分布如图5所示，为详细展示不同模型的降尺度差异，文中从研究区选择出一块大小适中的区域，对比图5(a)、图5(b),结合图5(c)可以发现，模型2对水体、建筑用地的降尺度结果有着明显改善，模型1对水体的降尺度结果明显存在高估，而建筑用地表现为明显低估，空间差异分布情况与上述定量统计具有一至攵的结论。不同时相遥感影像反演的LST和4种光谱指数均存在差异，为避免这种差异对降尺度结果评价的影响，本文还分别应用两种模型对不同季节的MODISLST进行了降尺度，同样以对应时期反演的Landsat8TIRSLST作为当期地表温度的参考值对降尺度结果进行验证，验证结果见表4。由表4可知，4个季节两种模型的降尺度结果存在明显差异。总体而言，模型2的降尺度结果要优于模型1,但是不同季节提升效果不同。两种模型均是夏、秋季的降尺度结果优于春、冬季，这是由于NDVI受季节影响比其他3个指数要更加严重，并且NDVI温度降尺度方法用于春、冬季要弱于其他方法［5］。此外，分别从RMSE和MBE的差异幅度来看，模型2相比模型1的降尺度结果在春、冬季的提升效果均要优于夏、秋季，这是因为NDVI温度降尺度方法较适合于夏、秋季，因此分土地利用类型使用不同光谱指数和全局尺度上使用NDVI提升效果并不明显，相反，春、冬季处于弱势地位的NDVI被替换后，降尺度效果有着明显提升，此处得出的结论与文献[5]一致。针对现有研究中地表温度降尺度