MODIS数据地表温度反演劈窗算法比较_陈瀚阅

表1 MODIS 部分波段及其参数[14]波段 光谱范围 信噪比 主要用途 分辨率 1 620～670nm 128 陆地、云边界 250m 2 841～876nm 201 陆地、云边界 250m 19 915～965nm 250 大气水汽 1000m 31 10.780～11.280μm 0.05 地球表面和 云顶温度1000m 3211.770～12.270μm0.051000m劈窗算法介绍McMillin （1975年）最早提出了劈窗算法，最先是用于海面温度的反演，这种方法是利用2个相邻的热红外窗口大气水汽吸收特性的差异，把海面温度表达成2个热红外窗口亮度温度的线性组合。

Price （1984年）最先把劈窗算法推广到陆面温度的反演，通过引入比辐射率改正项来减小因陆地表面比辐射率变化而引起的误差。

Becker 从理论上证明了用分裂窗技术反演地表温度的可行性，并且第一次从理论上给出了使用分裂窗技术时大气和比辐射率对地表温度反演的影响。

Becker 和Li 根据热辐射传导的地方性特征，提出了著名的局地劈窗算法，已得到了较广泛的应用。

Wan 和Dozier 在Becker 和Li 的研究基础上，于1996年提出了一种广义的地表温度反演劈窗算法。

Sobrino 和Becker 用Lowtran 7对不同的大气、观测角度以及地表参数进行模拟，得出了各参数的表达式。

在这些表达式里，大气和比辐射率的作用是耦合在一起的。

而Sobrino 等则通过某些近似把这2种作用分开了，通过对大气向下热辐射的近似解和对Planck 辐射函数的线性化。

覃志豪等推导了劈窗算法，该算法仅需要2个因素来进行地表温度的演算，即大气透过率和地表比辐射率[15][ 16]。

在众多的劈窗算法中，覃志豪等提出的算法由于需要参数少、计算简单且精度较高，被认为是较好的算法之一。

本文主要针对这一算法进行介绍。

覃志豪[15]等提出的针对MODIS 数据反演地表温度的劈窗算法使用的公式如下:0131232Ts A AT A T =+- （1）其中：Ts 是地表温度，31T 、32T 分别是MODIS 第31、32通道的亮温。

第！〇卷，第1期国土资源遥感V。

'！0，）〇( 2018 年 3 月R EM O TESEN SIN G FO R LA N D&R ESO U R C ES Mar.，2018 doi： 10.6046/gtzyyg.2018. 01.07引用格式：陈瀚阅，朱利，李家国，等.基于L anded数据的2种海表温度反演单窗算法对比——以红沿河核电基地海域为例 [J].国土资源遥感，2018，30( 1)：45 -53. (Chen H Y，Zhu L，Li J G，et a' A comparison of t'w o mono- window algoritlims for retrie­ving sea surface temperature from Landsat^^data in coastal water of Hongan River nuclear power station#J].Remote Sensing for Landand Resources，2018，30( 1) ：45 -53.)基于Landsat8数据的2种海表温度反演单窗算法对比—以红沿河核电基地海域为例陈瀚阅U，3,朱利4!李家国5!范协裕1>2,3(1.福建农林大学资源与环境学院，福州350002$2.福建农林大学土壤生态系统健康与调控福建省高校重点实验室，福州350002;3.福建农林大学福建省土壤环境健康与调控重点实验室，福州350002;4.环 境保护部卫星环境应用中心，北京100094;5.中国科学院遥感与数字地球研究所，北京100101)摘要：以辽宁省红沿河核电站附近海域为研究区，对单窗算法用于Landst8 TIRS数据反演沿海海表温度（s a sur­face temperature，SST)的适用性进行比较分析。

首先，基于大气廓线数据 （thermodynamic initial guess retrieval，TIGR)，针对Landst8 TIRS第10波段修订QK&B算法系数；然后，从星地同步验证和参数敏感性2方面对辐射传 输模型（radiation transfer model，RTM)和QK&B算法进行对比分析。

基于劈窗算法的Landsat8影像地表温度反演吴亮;姚昆【期刊名称】《贵州电力技术》【年(卷),期】2018(021)004【摘要】陆地表面温度(LST)是表征地表能量交换和地面特征的重要指标,目前遥感技术逐渐成为区域和全球尺度上LST反演的一种便捷工具,而采样不同算法及不同影像的热红外遥感LST反演研究层出不穷,其中基于Landsat数据的反演成果尤为突出.文章利用劈窗算法对Landsat 8遥感影像进行地表温度反演,对比探讨了根据经验值与借助MODIS热红外数据两种不同方式的LST反演结果,并进行北京市热红外波段辐射亮度温度比较,针对地表温度分级进行统计,分析了当地地表温度分布趋势.结果表明:劈窗算法下Landsat 8数据的反演温度更接近实际温度,精度较高且优于MODIS产品;北京市地表温度空间分布格局受地物结构与反射率所制约,高温区主要集中分布于中东部,中低温区分布与林地及水体分布结构较为吻合.【总页数】8页(P18-25)【作者】吴亮;姚昆【作者单位】中国电建集团贵州电力设计研究院有限公司,贵州贵阳550002;西昌学院资源与环境学院,四川西昌615000【正文语种】中文【中图分类】K99【相关文献】1.基于劈窗算法的青岛地区地表温度反演研究 [J], 孙乐乐;金宝轩2.基于Landsat－8数据和劈窗算法的地表温度反演及城市热岛效应研究 [J], 宋挺;段峥;刘军志;严飞;黄君;吴蔚3.基于Landsat8影像的地表温度反演算法应用分析 [J], 于梦馨;刘波4.基于Landsat8影像和劈窗算法的济南地区地表温度反演 [J], 王菲5.基于劈窗算法的Landsat8影像地表温度反演 [J], 吴亮;姚昆;;因版权原因，仅展示原文概要，查看原文内容请购买。

基于MODIS数据的地表温度反演
杨鹏;杨世植;夏鹏
【期刊名称】《大气与环境光学学报》
【年(卷),期】2011(6)4
【摘要】利用MODIS影像数据,采用劈窗算法来反演安徽地区的地表温度。

结合Sobrine、覃志豪等提出的NDVI^(TEM)方法和地物监督分类方法,对地表比辐射率进行了估算,将反演结果与NASA的地表温度产品进行比较,平均误差在1 K左右.同时利用卫星过境当天从安徽省高密度自动监测站获取的实时数据对反演结果进行验证,发现反演温度与地面实测数据的曲线走势具有高度的一致性且有较高的相关性,能直观地反映安徽地区地表温度的空间分布。

【总页数】7页(P287-293)
【关键词】地表温度;MODIS;劈窗算法;地表比辐射率
【作者】杨鹏;杨世植;夏鹏
【作者单位】中国科学院安徽光学精密机械研究所中国科学院通用光学定标与表征技术重点实验室
【正文语种】中文
【中图分类】P314
【相关文献】
1.基于MODIS数据的云南省地表温度反演研究 [J], 徐虹;戴丛蕊
2.基于劈窗算法MODIS数据建立地表温度反演模型 [J], 何晓琳;吴笑天
3.基于遗传算法综合Terra/Aqua MODIS热红外数据反演地表组分温度 [J], 孙珂;陈圣波
4.基于MODIS数据的乌昌地区地表温度反演 [J], 石玉;肖继东;张旭
5.基于MODIS数据的黑河流域地表温度反演研究 [J], 张福存;李净;吴立宗;张晓因版权原因，仅展示原文概要，查看原文内容请购买。

第41卷第1期华中师范大学学报(自然科学版)Vol.41No.12007年3月 J OU RNAL OF HUAZHON G NORMAL UNIV ERSIT Y (Nat.Sci.) Mar.2007收稿日期:2006206203.基金项目:国家“863”项目(2006AA12Z128).3通讯联系人.Email :lrd @.文章编号:100021190(2007)0120143205基于MOD IS 数据的湖北省地表温度反演研究许国鹏1,2,李仁东13,刘可群3,张　斌1,2(1.中国科学院测量与地球物理研究所,武汉430077;2.中国科学院研究生院,北京100039;3.武汉区域气候中心,武汉430074)摘　要:利用遥感影像获取地表温度信息,对于监测地表状况和生态环境变化有重要意义.利用MODIS 影像和劈窗算法计算了反演地表温度的关键参数:大气透过率和地表比辐射率,并针对湖北省水域单独提取出水体像元来计算地表比辐射率,估算出湖北省地表温度,与地面同步实测数据比较表明,平均误差为0.51℃,精度比较高.关键词:地表温度;MODIS ;精度;大气透过率;地表比辐射率中图分类号:P23文献标识码:A 地表温度是研究区域地表能量平衡和资源环境变化的重要参数之一[1],不仅可以为森林火灾和工厂热污染排放监测提供直接依据,而且还在以地表温度为基础的相关模型进行区域土壤水分估算、农业旱情监测、城市热岛效应、军事伪装辨别等方面都要重要的实际应用价值.相对星上亮温,地表温度才能作为一个重要的基本参数直接参与相关模型的计算,如“全球环流模型”、地表潜热、显热通量方程、土壤热流方程等[2].由于地表温度在时空上的动态分布变化特性,依靠地面观测站大面积获取地表温度参数、宏观地把握其时空分布规律是很难的.而借助遥感技术大范围、快速地反演地表参数的优势,可利用热红外波段获取地表温度信息,进行资源环境的动态监测和研究.由于陆地表面的复杂性和大气影响的许多因素,地表真实温度反演是一个非常复杂的过程.在对大气热辐射传输和相关参数各种近似、假设的基础上,先后提出了多种地表温度的反演算法,其中劈窗算法是目前为止发展最为成熟的地表温度反演算法[325].本文利用MODIS 数据的劈窗算法反演湖北省地表温度,并将反演结果与地面同步实测数据进行了精度分析.1MODIS 数据的地表温度反演方法1.1劈窗算法劈窗算法最初是根据地表热辐射传导方程,利用AV HRR 大气窗口内热红外第4、5两个相邻通道对大气吸收作用的差异,通过两个通道亮度温度的各种组合来剔除大气的影响,进行大气和地表比辐射率的订正来获取地表温度的[2].AV HRR 的两个热通道10.5～11.3μm 、11.5～12.5μm 与MODIS 第31波段(10.7805～11.280μm )和32波段(11.770～12.270μm )的中心波长基本对应,毛克彪[4]等人据此提出了一种利用MODIS 数据31、32波段估算地表温度的劈窗算法.该劈窗算法的表达式如下:T S =A 0+A 1T 31+A 2T 32,(1)其中:T S 为地表温度,T 31、T 32是MODIS 31/32波段的亮度温度,单位是℃,可用Planck 辐射方程获取;A 0、A 1和A 2是系数,可用式(2)计算:A 0=-64.6036E 1-68.7258E 2-273.16,(2a )A 1=1+A +0.440817E 1,(2b )A 2=-(A +0.473453E 2),(2c )E 1=D 31(1-C 31-D 31)/E 0,(2d )E 2=D 31(1-C 32-D 32)/E 0,(2e )A =D 31/E 0,(2f )E 0=D 32C 31-D 31C 32,(2g )C i =εi τi ,(2h )D i =(1-τi )(1+(1-εi )τi ),(2i )其中,E 1、E 2、E 0、C 31、D 31、C 32、D 32均为中间变量,可迭代消除;τi 为大气透射率;εi 为地表辐射率.算144　华中师范大学学报(自然科学版) 第41卷法的关键是计算大气透射率τi 和地表比辐射率εi .1.2大气透射率的求算大气透射率是地表辐射、反射透过大气到达传感器的能量与地表辐射能、反射能的比值,它与大气状况、高度等因素有关.对于热红外波段,最重要的大气变化是大气温度和水汽的变化.在天气稳定情况下,水汽含量是影响大气透射率的主要因素.Kauf man 、Bo 2Cai Gao [6]利用MODIS 第19和第2波段模拟出了大气水汽含量的表达式:w =[(α-ln τw )/β]2,(3)其中:w 是指大气水汽含量;τw 是大气水汽吸收波段第19波段地面反射率与大气窗口波段第2波段地面反射率的比值;α、β是参数,对于复合性地表,α=0.02,β=0.651.MODIS 31/32夏季中纬度标准大气状况下大气水汽含量和透过率的变化之间呈近似线性关系[7]:τ31=-0.10671w +1.04015,(4a )τ32=-0.12577w +0.99229.(4b )1.3地表比辐射率的求算地标比辐射率是物体与黑体在同温度、同波长下的辐射出射度的比值.在传感器的波段区间及像元大小确定情况下,地表比辐射率主要取决于地表物质的组成和结构.在MODIS 1km 的像元尺度下,像元可以粗略视作由水体、植被和裸土3种类型构成.MODIS 混合像元的地表比辐射率可表示为[7]:εi =P w R w εiw +P v R v εiv +(1-P w -P v )R s εis ,(5)其中,P w 和P v 分别是水面和植被在该像元内的构成比例;εiw 、εiv 和εis 分别是水面、植被和裸土在该波段的辐射率,可在ASTER 提供的常用地物比辐射率光谱库内查得;R w 、R v 和R s 分别是水体、植被和裸土的温度比率,在5～45℃范围内,分别为1.00744、0.99240和0.99565[5].该方程的求算关键在于估算混合像元中的P w 和P v 值.对于水面较大的地区来说,可以利用可见光和红外波段水体反射率一般明显低于其它地物以及水体归一化植被指数NDV I <0的特性,提取纯水体像元,并取P w =1.此时,εi =R w εiw .对于水面可以忽略的陆地来说则构成比较复杂,像元中植被构成比例可以表示为:P v =(N DV I -N DV I s )/(N DV I v -N DV I s ),(6)其中,N DV I v 和N DV I s 表示完全植被和裸土的植被指数;N DV I 表示任意像元的植被指数.地表比辐射率可根据像元NDV I 来求算:①当N DV I >N DV I v 时,像元被看作是完全的植被覆盖,取P v =1,则εi =R v εiv ;②当N DV I s <N DV I <N DV I v 时,εi =P v R v εiv +(1-P v )R s εis ;③当N DV I <N DV I s 时,像元被看作完全裸土,取P v =0,则εi =R s εis .2数据处理本文采用2005年10月10日HDF 格式1B 级别的MODIS 数据.首先进行预处理:①对相邻扫描行之间数据重复的“蝴蝶结”现象,利用ENV I 软件IDL 模块中开发的纠正函数去除;②进行物理定标,将DN 值转化为反射亮度或辐射亮度;③对于有云的区域,鉴于云在多波段的光谱特征互补性,可利用热温度信息和云检测指数进行云的综合检测[8].数据处理过程中的编程、计算均在Matlab 环境中完成,主要步骤如下:①用第1、2波段计算N DV I ,对比反射率提取水体;②根据湖北省秋季植被覆盖情况,取N DV I v =0.70,N DV I s =0.05,估算水面和植被的构成P w 、P v ;③根据式(5)及P w和P v ,估计像元的地表比辐射率ε31、ε32;④用第2、19波段计算大气水分含量w ,并进而根据式(4)估计大气透过率τ31、τ32;⑤用第31、32波段的辐射亮度,根据Planck 方程计算星上亮度温度T 31、T 32;⑥运用劈窗算法式(1)、(2),利用ε31和ε32、τ31和τ32以及T 31和T 32,计算地表温度Ts;⑦利用数据自带的坐标信息做精几何校正,叠加省界掩膜,得到湖北省地表温度反演结果图.具体反演流程如图1.图1　地表温度反演流程图Fig.1　Flow chart of retrieved L ST　第1期许国鹏等:基于MODIS数据的湖北省地表温度反演研究145　3结果分析图2　湖北省地表温度分布图Fig.2　Distribution map of L ST in Hubei Province 图2为MODIS反演结果得到的湖北省地表温度空间分布图.经过云检测,鄂西北竹溪、竹山、鄂西南约30°N沿线部分地区和鄂东南黄梅等地区反演温度明显偏低是由于云的覆盖造成的.统计可知,遥感反演的湖北省地表温度平均值为14.24℃;空间分布差别不大,97.66%的地区温度在13℃到27℃之间,其中17℃、18℃、22℃、23℃、24℃的比例均达到了9%,基本符合湖北省10月上中旬白天的季节特点.3.1反演结果精度评价Terra卫星在中国的过境时间为上午10∶30～12∶00,本文利用湖北省气象局野外观测站上午11时观测的地表温度数据对模拟结果进行精度分析.在75个野外站点记录的数据中,由于竹溪、竹山、松滋、黄梅4个观测站刚好处于云覆盖区域,数据缺乏可比性,因此采用其它71个观测站的数据进行评价,如表1.地表实测温度(T m)和反演结果(T r)的平均误差(5)为[9]5=1n ∑ni=1T ri-T mi.(8)分析表明,71个站点的平均误差为0.51℃,反演精度在0.5℃以内为57.7%,反演精度在1.0℃以内为31.0%,反演精度在1.2℃以内为9.9%,反演精度在1.7℃以内为1.4%.部分站点误差相对较大的原因可能是由于影像像元与地面实测点匹配的不确定性造成的.3.2对真实地表温度的模拟将71个站点分别作为拟合和验证数据,用实测值和反演温度进行拟合获取回归方程,再由反演温度模拟真实地表温度.根据数理统计原理,由计算机随机抽取56个观测站点处的反演温度和实测表1　2005年10月10日MOD IS反演温度与地面同步实测值的比较/℃Tab.1　Comparison of retrieved L ST f rom MODIS and measured value(Oct.10,2005)/℃观测站点实测值(T m)反演温度(T r)误差(5)观测站点实测值(T m)反演温度(T r)误差(5)郧西17.0017.430.43枝江23.4023.670.27郧县20.1019.980.12潜江22.5023.040.54十堰18.6017.43 1.17荆州24.2023.800.40房县17.0017.690.69公安24.0024.310.31丹江口19.8019.340.46应城24.1024.310.21老河口20.0019.60.40孝感24.0023.420.58谷城20.0020.490.49天门22.6021.630.97襄樊22.4022.650.25仙桃23.5022.910.59枣阳22.4022.780.38汉川22.7022.140.56巴东19.0018.200.80汉阳25.0025.070.07秭归21.0020.490.51黄陂23.6023.290.31兴山18.4017.560.84新州23.6023.930.33保康17.7018.580.88武汉22.6022.140.46神农架16.3016.920.62鄂州23.5023.670.17南漳20.5019.34 1.16大冶24.4024.180.22远安20.4019.34 1.06咸丰20.1019.980.12宜城23.2023.160.04宣恩20.1020.230.13荆门21.1020.230.87鹤峰17.9016.920.98钟祥21.8021.760.04来凤22.3020.62 1.68随州20.0020.620.62石首23.9023.930.03广水23.8023.930.13监利25.0025.330.33京山23.8023.420.38洪湖19.0019.600.60安陆24.2024.440.24蒲圻23.5023.930.43云梦24.4024.820.42嘉鱼23.4023.040.36大悟24.1023.670.43崇阳21.4020.23 1.17红安22.5023.040.54通城25.0025.200.20麻城24.1024.560.46咸宁25.4024.820.58利川19.2019.730.53通山23.6023.420.18建始20.1018.96 1.14罗田24.4024.180.22恩施19.6019.220.38英山23.2022.14 1.06宜昌县23.3022.650.65浠水23.4023.040.36五峰18.1016.92 1.18黄石23.5023.930.43当阳20.9021.890.99蕲春23.3022.780.52宜昌23.3022.650.65阳新23.0023.290.29长阳22.6022.650.05武穴22.5022.140.36枝城22.2022.270.07值进行一元线性回归(如图3),拟合结果表明,在95%的置信水平下,两者之间一元线性相关性R2146　华中师范大学学报(自然科学版) 第41卷图3　反演温度与实测值的一元线性回归Fig.3　Unary linear regression between retrievedL ST and measured value为93.7%,结果比较理想.利用该线性拟合方程和另外15个站点处的反演温度模拟真实地表温度,并用实测值来验证模拟精度,如表2.分析发现,从整体上看,模拟值(Ts)比反演值更接近实测值,平均误差由原来的0.53℃减小到0.43℃.这说明,用该一元线性回归方程直接模拟区域内各像元的地表温度,相比反演值,更接近真实地表温度,具有更好的效果.4结论本文利用基于劈窗算法的地表温度反演算法,通过采用MODIS可见光第1波段、近红外第2波段、中红外第19波段提取大气透射率和地表比辐射率,结合热红外波段的亮温信息,反演了湖北省的地表温度.结果表明:反演精度比较高, 71个站点的实测值与反演结果平均误差为0.51℃;可以用实测值和反演温度拟合的一元线性回归方程直接模拟各像元的地表温度,能更接近真实地表温度.表2　一元线性模拟值与反演值的比较/℃Tab.2　Comparison of simulated and retrieved temperature/℃|T r-T m||T s-T m||T r-T m||T s-T m||T r-T m||T s-T m||T r-T m||T s-T m| 1 1.060.7650.460.4590.120.14130.580.59 20.040.126 1.140.82100.580.61140.710.56 30.420.3970.650.55110.860.25150.220.21 40.540.6280.070.2120.530.13平均0.530.43 这种劈窗算法能简单、快速、准确地反演地表温度,避免了对地面气象数据的依赖,达到了遥感技术大面积对地表进行资源环境宏观观测的目的.参考文献:[1]　Li Z,Becker F.Feasibility of land surface temperature andemissivity determination from AV HRR data[J].Remote Sensing of Environment,1993(43):67285.[2]　赵英时.遥感应用分析原理与方法[M].北京:科学出版社,2004.[3]　Qin Zhihao,Karnieli A,Berliner P.A mono2window algo2rit hm for retrieving land surface temperature from Landsat TM data and it s application to t he Israel2Egypt border region [J].Int J Remote Sens,2001,22(18):371923746.[4]　毛克彪,覃志豪,施建成.针对MODIS数据的劈窗算法研究[J].武汉大学学报:信息科学版,2005,30(8):7032708. [5]　Becker F,Li Z.Towards a local split window met hod overland surface.Int J Remote Sens,1990,11(3):3692393. [6]　Y oram J Kauf man,Gao Bo2Cai.Remote Sensing of WaterVapor in t he Near IR from EOS/MODIS[J].IEEE Transac2 tions on Geoscience and Remote Sensing,1992,5(30): 8712884.[7]　毛克彪,覃志豪,王建明,等.针对MODIS数据的大气水汽含量及31和32波段透过率计算[J].国土资源遥感,2005,63(1):26229.[8]　宋小宁.基于植被蒸散法的区域缺水遥感监测方法研究[D].中国科学院博士学位研究生论文,2004.[9]　毛克彪,覃志豪,宫　鹏.劈窗算法L ST精度评价和参数敏感性分析[J].中国矿业大学学报,2005,34(3):3182322. 第1期许国鹏等:基于MODIS数据的湖北省地表温度反演研究147 Land surface temperature retrieved fromMODIS d ata in H ubei ProvinceXU Guopeng1,2,L I Rendong1,L IU Kequn3,ZHAN G Bing1,2(1.Institute of Geodesy and G eophysics,Chinese Academy of Sciences,Wuhan430077;2.Graduate University,Chinese Academy of Sciences,Beijing100049;3.Wuhan Regional Climate Center,Wuhan430074)Abstract:Land surface temperat ure information can be retrieved f rom Remote Sensingimage,which is important and significant to observe macroscopically t he erat h’s surfacestat u and ecological and environmental change.The essay ret rieves t he at mo sp herictransmittance and t he surface emissivity which are required when L ST is calculated u2sing MODIS image by split2window algorit hm.The surface emissivity is calculated bymeans of o btaining t he water p ure pix individually in t hat t he water area in Hubei Prov2ince can not be ignored.Then,t he values of L ST in Hubei Province are derived in t hismet hod.The comparison with the ground synchronization observation data shows that the to2tal average error is0.51℃,therefore,the estimated precision is reasonable quitely.K ey w ords:land surface temperat ure;MODIS;p recision;at mosp heric t ransmittance;surface emissivity (上接第137页)The pattern innovation of B eijing City competitiveness evolvement and upgrade in the contexts of globalizationMA Qingbin,WEN Hui(China Center for Town Reform and Development,Beijing100045;Research Center of Economy Environment,Chinese Academy of Sciences,Beijing100085) Abstract:Based on the description of concept and evaluation method of urban competitiveness,the paper quantitatively analyses Beijing city competitiveness and the evolvement of industrystructure under new industry type through defining the concept of city competitiveness andconstructing new evaluation method,brings forward choosing the best pattern to improving thecity competitiveness of Beijing and gives proposals on correlative policy in the end.It wasshowed with the data that the urban competitiveness of Beijing from1995to2003has beenpromoted during the city evolution process on a global scale.Y et Beijing’s urban competitive2ness still has a pretty big ga Pcompared with other equal level city.The rationality of the envi2ronment policy in the era of globalization can not be judged by its strictness degree,but be ap2praised on whether the environment policy coincide with the level of urban competitiveness.The urban policy2makers ought to either find the reasonable binding site among the urban de2velopment strategy,the industry policy and the environment policy or choose the appropriate de2velopment pattern in order to impulse the city to enhance its competitiveness fleetly,steadily,proportionally and sustainable with the least resources input and the minimal environmental cost.K ey w ords:global city;city competitiveness;pattern。

式的系数也不相同，如下表１：表１不同观测角度与分裂窗算法系数的对应表观测角度ＡｌＡ２Ａ３ＢｌＢ２Ｂ３Ｃ０１．０３５９０．１９７２－０．０９７９０．５７９１．０．２７３９１．９００１．１０．１１１６１０１．０３６７９５Ｏ．１３７４７９６．０．７７９７１０５３．６８９２１７３．１０６１８０３．１０３７５９－６．１２７７６２．２０１．０４４８０４０．１４６２１８６．ｏ．７０３３９６６３．３２９９９０２．６６５４６９２．５１９３０９．８．３６３３７３０１．０４２６５６Ｏ．１６４９７９２＿０．４８６６９８２２．３４５１６５１．７２１３９１２．０４７９３０－８．１４１５５５３５１．０３７６７９Ｏ．１９８５２１７．０．０８３３９０７０．４９５１５０２．０．２２９５７３４０．５０６８８５６．８．０７６９６５４０１．０４０７０．１９８９－０．０７４５０．４９００．０．４３０３１．９６５６．１０．７２５４４２１．０４８１０４０．１４８８１６２－０．６８４５１０６３．３０９７６１２．４９３４１８３．２３６０９１９．７．８２８６６２４４１．０５４２４７０．１３６６４０４．０．８４６６７０１４．０４１５９２３．２４３１３３３．５７４４２５．８．５３０２５４４６１．０６２５０９Ｏ．１１７３４９５．１．０８６４４１５．１９９５９８４．０６６４２８４．６８７８７５－９．１６６２４３４８Ｉ．０６１３４２０．１３１２９８５．０．９４２４９３１４．５１４９１０３．７２９６６９３．７６５７３８．９．１５６９９７５０１．０３９７９４Ｏ．１８４２９２４．０．２５２６２２７１．２６５３５４Ｏ．７１８１８２５１．９１５２４９－６．５３４４４７５２１．０２２４１３０．２１１７１７３Ｏ．１８０８９４６．ｏ．５９６３７３８－２．０８５９９８２．２０１６１２－３．８４０６７５５４１．０２０９１２０Ｏ．２１７６８４４０．２６２９８５ｌ．１．０２４１１ｆ－２．３３８８６８１．６４９２０５．３．７０５１９３５６１．００３５０５０．２３８“ｌＯ．６５１６７３－２．６７１４４－４．８０６８０６２．０１９６１９－１．３７５３７６５８１．００８５８００．２２５８９３６０．４９７３８３６．Ｉ．９８１０７－３．７８４２２８２．７２７３“．１．４４４４０５平均１．０３７３３０．１７７０２９－０．２９６３６１．５３２４５６０．５１９６１３２．５２１３８５－６．８７５２６Ｍ０ＤＩＳ扫描观测宽度达２３３０公里，观测角度达到±５８。

表1 MODIS 部分波段及其参数[14]波段 光谱范围 信噪比 主要用途 分辨率 1 620～670nm 128 陆地、云边界 250m 2 841～876nm 201 陆地、云边界 250m 19 915～965nm 250 大气水汽 1000m 31 10.780～11.280μm 0.05 地球表面和 云顶温度1000m 3211.770～12.270μm0.051000m劈窗算法介绍McMillin （1975年）最早提出了劈窗算法，最先是用于海面温度的反演，这种方法是利用2个相邻的热红外窗口大气水汽吸收特性的差异，把海面温度表达成2个热红外窗口亮度温度的线性组合。

Price （1984年）最先把劈窗算法推广到陆面温度的反演，通过引入比辐射率改正项来减小因陆地表面比辐射率变化而引起的误差。

Becker 从理论上证明了用分裂窗技术反演地表温度的可行性，并且第一次从理论上给出了使用分裂窗技术时大气和比辐射率对地表温度反演的影响。

Becker 和Li 根据热辐射传导的地方性特征，提出了著名的局地劈窗算法，已得到了较广泛的应用。

Wan 和Dozier 在Becker 和Li 的研究基础上，于1996年提出了一种广义的地表温度反演劈窗算法。

Sobrino 和Becker 用Lowtran 7对不同的大气、观测角度以及地表参数进行模拟，得出了各参数的表达式。

在这些表达式里，大气和比辐射率的作用是耦合在一起的。

而Sobrino 等则通过某些近似把这2种作用分开了，通过对大气向下热辐射的近似解和对Planck 辐射函数的线性化。

覃志豪等推导了劈窗算法，该算法仅需要2个因素来进行地表温度的演算，即大气透过率和地表比辐射率[15][ 16]。

在众多的劈窗算法中，覃志豪等提出的算法由于需要参数少、计算简单且精度较高，被认为是较好的算法之一。

本文主要针对这一算法进行介绍。

覃志豪[15]等提出的针对MODIS 数据反演地表温度的劈窗算法使用的公式如下:0131232Ts A A T A T =+- （1）其中：Ts 是地表温度，31T 、32T 分别是MODIS 第31、32通道的亮温。