用MODIS影像和劈窗算法反演山东半岛的地表温度.pdf

第3 4 卷第1 期 2 0 0 5 年1 月 中国矿业大学学报 J o u r n a lo fC h i n aU n i v e r s i t yo fM i n i n g T e c h n o l o g y V 0 1 ．3 4N o ．1 J a n ．2 0 0 5 文章编号1 0 0 0 1 9 6 4 2 0 0 5 0 1 0 0 4 6 0 5 用M O D I S 影像和劈窗算法 反演山东半岛的地表温度 毛克彪1 ’2 “，覃志豪2 ”，施建成1 1 ．中国科学院遥感所，北京1 0 0 1 0 1 ； 2 ．农业部资源遥感与数字农业重点实验室，北京1 0 0 0 8 1 ； 3 ．南京大学国际地球系统科学研究所，江苏南京2 1 0 0 9 3 摘要简要介绍了劈窗算法，并对P l a n c k 方程进行了线性简化，从而简化了劈窗算法的方程计 算．针对M O D I S 的波段设置特点，先用近红外波段反演大气水汽含量，然后根据L O W T R A N 模 拟大气水汽含量与透过率的关系，这样可以根据反演得到的水汽计算3 1 和3 2 波段的透过率．最 后运用劈窗算法反演山东半岛地区的地表温度 I 。S T ，结果比较合理． 关键词L S T ；M O D I S ；劈窗算法；透过率 中图分类号7 5 3文献标识码A R e t r i e v a lo fL a n dS u r f a c eT e m p e r a t u r eo fS h a n d o n gP e n i n s u l a U s i n gM O D I SI m a g ea n dS p l i tW i n d o wA l g o r i t h m M A OK e b i a 0 1 ’2 ～，Q I NZ h i h a 0 1 ～，S H IJ i a n c h e n 9 1 1 ．I n s t i t u t eo fR e m o t eS e n s i n gA p p l i c a t i o n s ，C h i n e s eA c a d e m yo fS c i e n c e ，B e i j i n g1 0 0 1 0 1 C b i n a ； 2 ．T h eK e yL a b o r a t o r yo fR e m o t eS e n s i n ga n dD i g i t a lA g r i c u l t u r e ，t h eM i n i s t r yo fA g r i c u hu r e ， B e i j i n g10 0 0 8 1 ，C h i n a ； 3 ．I n t e r n a t i o n a lI n s t i t u t ef o rE a r t hS y s t e m ，N a n j i n gU n i v e r s i t y ，N a n i i n g ，J i a n g s u2 1 0 0 9 3 ，C h i n a A b s t r a c t T h es p l i t w i n d o wa l g o r i t h mi sb r i e f l yi n t r o d u c e d ．T h ec o m p u t a t i o no fs p l i t w i n d o w a l g o r i t h mi ss i m p l i f i e d ．A c c o r d i n gt ot h ec h a r a c t e ro fM O D I S ，t h ew a t e rc o n t e n ti Sd e r i v e df r o m t h ei m a g ea n dt h et r a n s m i t t a n c e so f31a n d3 2b a n d sa r ec a l c u l a t e dt h r o u g hb u i l d i n gt h er e l a t i o n s h i p b e t w e e nt h ew a t e rc o n t e n ta n dt h et r a n s m i t t a n c e ．F i n a l l y ，t h el a n ds u r f a c et e m p e r a t u r e I 。S T o f S h a n d o n gp e n i n s u l ai sr e t r i e v e du s i n gt h es p l i tw i n d o wa l g o r i t h m ． K e yw o r d s L S T ；M O D I S ；S p l i t w i n d o wa l g o r i t h m ；T r a n s m i t t a n c e 地面温度是反映土壤一植被一大气系统能量 流动与物质交换以及土地资源和环境管理的重要 参数[ 1 ] ，也是地球表面能量平衡和温室效应的一个 指标，是区域和全球尺度地表物理过程的一个关键 因子．因此，地面温度信息对于气候、水文、生态学 和生地化学等许多领域中是非常必须的[ 2 ] ．但是依 靠地面观测站的观测来大面积地获取地表温度参 数是不现实的．而借助于热红外遥感影像，可以方 便快捷地获得大面积，甚至全球的地温资料，且数 据更新快，成本低廉． M O D I S 是美国国家航空航天局、日本国际贸 易与工业厅和加拿大空间局、多伦多大学共同合作 发射的卫星T E R R A 上的一个中分辨率传感器．它 包含3 6 波段，分别针对陆地、海洋、水汽、气溶胶等 收稿日期2 0 0 4 0 5 0 8 基金项目国家重点基础研究发展规划 9 7 3 项目 2 0 0 1 C B 3 0 9 4 0 4 ；国家自然科学基金项目 4 0 4 7 1 0 9 6 作者简介毛克彪 1 9 7 7 一 ，男，湖南省沅江市人，博士生，主要从事热红外，微波遥感，空间数据挖掘及G I S 应用等方面的研究 万方数据 第1 期毛克彪等用M O D I S 影像和劈窗算法反演山东半岛的地表温度 4 7 来设计的．其中有6 个波段是专门针对陆地温度来 设计的．M O D I S 与T M 不同，T M 适合于小区域、 高精度的反演，而M O D I S 则适用于反映大区域， 甚至全球的温度分布规律． 劈窗算法主要是针对N O A A 卫星A V H R R 的第4 1 0 ．5 ～1 1 ．3p m ，5 1 1 ．5 ～1 2 ．5 弘m 热红 外通道提出来的．从P r a b h a k a r a 等开发的第一劈 窗算法算起，发展已经有2 0 多年了[ 3 ] ．相对单窗算 法，劈窗算法已经比较成熟．这个算法有很多版本， 但都是基于辐射传输方程来设计的，主要区别在于 针对不同的情况，对大气参数的计算求解不一样． 对地观测卫星M O D I S 数据有1 7 个热波段 2 0 ～ 3 6 ，其中6 个波段 2 0 ～2 3 ，3 1 ，3 2 可用于分析地 表热辐射和地表温度的变化[ 4 ] ．本文拟通过对 P l a n c k 方程的辐射强度和温度计算模拟，发现温度 和辐射强度的近似关系，由此建立线性方程，以简 化计算．然后通过模拟计算得到大气透过率与大气 水汽含量的关系，从同一景M O D I S 影象中计算得 到M O D I S 3 1 ／3 2 波段透过率，最后运用劈窗算法 和M O D I S 的3 1 和3 2 通道反演山东半岛的温度 分布，并对反演结果进行分析． 1地表温度反演方法及方程简化 由于环境管理和土地资源研究的深入，通过遥 感方法反演地表温度已成为遥感技术的一个重要 的研究课题．劈窗算法主要是针对A V H R R 的第 4 ，5 两个热红外波段，根据地表热辐射传导方程推 导出来的．到目前为止，已经公开发表的劈窗算法 有1 7 个，这些算法主要区别在于计算形式和计算 系数[ 51 8 ] ．最通用的表达式为 T ，一T 4 A 7 ’l T 5 B ， 1 式中7 I 为地表温度，L 和丁。为A V H R R 通道4 和5 的亮度温度；A 和B 是与大气状态、视角、地 表反射率有关的系数．由于受到大气状态参数和地 表反射率的影响，所以不管各种算法的系数计算表 达式的差别有多大，一般只有2 种方法．一是通过 地面实际测量来校正亮度温度；二是大气辐射传输 模型，如6 S ，L O W T R A N ，M O D T R A N 等，通过标 准大气剖面资料来模拟计算．地表发射率对地表温 度的反演有很大的影响．通常，对大多数地物发射 率一般可以合理地假定为0 ．9 5 ～0 ．9 8 ，但在通道 4 ，5 发射率的区别对地表温度的反演也有影响． 本文考虑大气参数的计算，选择覃志豪在文献 E 1 8 ] 中提出的算法对该地区的地表温度进行反演． 该算法的优点是参数比较简单，而且精度比较高． 计算公式如下 B ， 7 1 i 一e i r i O B ， 丁， [ 1 一“ y i 口 ] [ 1 1 一￡， r 。 口 ] B ， 丁。 ， 2 式中B ， 7 T i 为波段i 的传感器辐射亮度；矗为i 波 段的地表反射率忍为波段i 的透过率；丁。表示大 气平均温度．针对M O D I S 热红外通道3 1 和3 2 ，分 别用下式表示 B 3 l 丁3 1 e 3 1 r 3 l 咿 B 3 1 E [ 1 一龟1 口 ] [ 1 1 一￡3 1 r 3 1 口 ] B 3 1 丁。 ， 3 B 3 2 丁3 2 E 3 z r 3 2 曰 B 3 2 T 5 [ 1 一r 3 2 口 ] [ 1 1 一￡3 2 r 3 2 口 ] B 3 2 L ． 4 为了求解 3 和 4 ，覃志豪在文献[ 1 8 ] 中通过 对P l a n c k 方程的泰勒展开，简化了计算过程．根据 P l a n c k 方程，辐射强度与温度的关系如下 r 2 B a 7 1 一丽菇≯丽- 5 根据式 5 对M O D I S 的第3 1 波段和3 2 波段 做温度 o ～5 0C 之间与辐射强度的散点图，发现 温度与辐射率在区间 0 ～5 0C 之间，近似线性关 系． f 曼 ； 每 ■ 基 邑 ＼ 魁 秘 摈 驿 2 7 3 2 7 8 2 8 32 8 8 2 9 32 9 8 3 0 3 3 0 8 3 1 3 3 1 8 3 2 3 7 y K 图1第3 1 和3 2 波段辐射强度随温度的变化关系 F i g ．1 R e l a t i o n s h i po ft h et e m p e r a t u r ea n dr a d i a n c e f o rM O D I Sb a n d s3 1a n d3 2 对于波段3 1 1 1 ．0 3 肚m ，从图中可以看出，在 0 “ - - 5 0 度之间，辐射强度与温度基本呈线性关系． 线性方程为y 一0 ．1 4 x 一3 1 ．8 0 ，R 2 1 ．0 ． 同理，对波段3 2 1 2 ．0 2u m ，线性方程为y 一 0 ．1 2 z 一2 6 ．8 0 ，R 2 1 ．0 0 ．从图1 和模拟的相关关 系可以知道，在o ～5 0 ‘C 的区间内用线性模拟的温 度与辐射强度的关系精度是比较高的．在实际应用 中，一个地区的温差间隔在同一时刻，同一景影象 的温差一般小于5 0 ‘C ．因此，我们可以根据研究区 的实际情况，对P l a n c k 方程采用分段线性模拟以 提高近似精度． 利用上述2 个近似线性关系，可以将式 3 和 4 化简为 0 ．1 3 8 3 4 6 3 1r 3 1 T ，一0 ．1 3 8 3 4 T 3 1 3 1 ．8 0 1 4 8 8 3 1 r 3 1 一 1 一r 3 1 L 1 1 一￡3 1 ／“ 3 1 J 0 ．1 3 8 3 4 T ，，一3 】．8 0 1 4 8 一3 1 ．8 0 1 4 8 。 6 万方数据 4 8中国矿业大学学报第3 4 卷 0 ．1 1 9 5 2 E 3 2 r 3 2 T ，一0 ．1 1 9 5 2 T 3 2 2 6 ．8 0 5 4 8 3 2 r 3 2 一 1 一f 3 2 [ 1 1 一e 3 2 毛2 ] 0 ．1 1 9 5 2 T 。一2 6 ．8 0 5 4 一2 6 ．8 0 5 4 ． 7 式 6 ， 7 中7 1 扪丁。。可以从影象算出，只有 兀，L ，如白是未知的．L 可以通过2 个方程消去， 而丁，为所求，所以只剩如e 。．要精确求t ，必须知 道实时大气剖面资料，但通常很难获取，所以一般 是采用标准大气进行模拟计算或求经验表达式．本 文采用M O D T R A N 对中纬度地区大气辐射传输 进行模拟计算，求取t ．本文针对M O D I S 的波段设 置特点，先反演大气水汽含量，然后根据大气水汽 含量和透过率的关系，求取％下面介绍透过率的 计算方法． 2 透过率估算 对透过率的计算往往是用M O D T R A N 、 L O W T R A N 、6 S 等软件模拟计算．然而实时大气剖 面资料难以获取，往往用标准大气代替，再用地面 资料修正．但天气状况变化无常，所以用标准大气 模拟往往精度不高．覃志豪等[ 1 如针对T M 的影像 提出了求算透过率的线性模拟方法，就是用 M O D T R A N 来计算大气水汽含量和透过率的关 系，然后通过大气水汽含量来计算透过率．但对于 T M 和N O A A ／A V H R R 来说，大气水汽含量很难 通过影像反演得到，故这一表达式很难得到推广应 用，但对于M O D I S ，水汽含量可以从近红外波段反 演出来，从而使得从遥感影像来反演透过率成为可 能．从M O D I S 影像中反演水汽的推导和反演请参 见文献[ 2 0 3 ．表1 是用L O W T R A N 模拟中纬度夏 季大气水汽含量与透过率的部分对应数据，模拟的 表1大气水汽含量与透过率模拟结果 T a b l e1R e s u l t so fs i m u l a t i o nb e t w e e nt h e t r a n s m i t t a n c ea n dw a t e rc o n t e n t 水蒸气含量／ g c m 一2 夏季模拟 0 ．9 5 13 0 6 6 6 7 0 ．9 4 13 1 3 3 3 3 0 ．9 3 04 2 6 6 6 7 0 ．9 1 84 8 0 ．9 0 53 8 0 ．8 9 11 0 6 6 6 7 0 ．8 7 56 8 0 ．8 5 91 1 33 3 3 0 ．8 4 33 4 0 ．8 2 28 8 0 ．8 0 31 6 6 6 6 7 0 ．7 8 29 4 6 6 6 7 0 ．7 6 17 2 6 6 6 7 0 ．7 3 99 0 6 6 6 7 波段范围适用于M O D I S 3 1 ／3 2 通道． 本文分别对上面两个波长范围的透过率与大 气水汽含量的模拟数据进行模拟进行指数关系模 拟．得到M O D I S 一3 1 ／3 2 ，水汽和透过率的关系为 第3 1 波段y 一5 ．7 2 4 ．6 9 e 4 2 ％s ， R 2 1 ．0 0 第3 2 波段y 一一1 ．2 5 2 ．2 8 e j 知， R 2 0 ．9 9 ． 3 结果分析 本文用到的M O D I S 影像是摄取于2 0 0 3 年6 月中旬，对山东半岛地区而言，正是初夏季节，是植 被生长最旺盛之时．为了更准确的反演地表温度， 先用M O D I S 的第2 和1 9 波段反演水汽含量，如 图2 所示 注意右下角为坏条带 ．从图2 的水汽 分布图可以看出，水汽在沿山东半岛海岸线地区水 汽含量相对较低，呈峰谷状分布，这种分布可能是 受地形和天气的影响比较大的缘故．根据本文模拟 大气水汽含量与透过率的关系表达式，分别计算得 到3 1 和3 2 波段的透过率图，如图3 所示．从图3 的剖面图分析可以知道，3 1 波段和3 2 波段的透过 率相差近0 ．0 5 ．其总体趋势基本相同．对于比辐射 率￡，本文根据A S T E R 提供的比辐射率经验给 出[ 20 | ，然后解方程 5 、 6 ，从而可以解出L ． 图2 水汽含量图 F i g ．2 R e t r i e v a lo fw a t e rc o n t e n t 经统计分析，M O D I S 数据中含有少量的噪声 数据，使得最高气温达4 1 ．1 7 2 ℃，在排除这些少量 的异常点之后，得到温度反演结果如4 所示．从图 4 的剖面分析中可知，3 6 ‘C 以上的区域只是很小的 一部分，其中在大部分在3 0C 以下．影像中包含了 黄海，海上的温度相对较低，从而使平均温度与最 高温度相差大，这是可以理解的． 从图4 还可以看出，最高气温分布在中部的丘 陵地区，环渤海地区的温度分布有些异常，在渤海 海岸沿线，温度高，而两面温度都低，在鲁南也有所 体现，这种现象与我们平常由陆地到海洋白天温度 万方数据 第1 期 毛克彪等用M O D I S 影像和劈窗算法反演山东半岛的地表温度4 9 图3M O D I S 的第3 l 和3 2 波段大气透过率图 F i g ．3 R e s u l t so ft h et r a n s m i t t a n c ee s t i m a t i o nf o rM O D I Sb a n d s3 1a n d3 2 o f l a n ds u r f a c et e m p e r a t u r ea n dg r o u n de m i s s i v i t y u s i n gN O A A A V H R Rd a t a [ J ] ．I n t ．J ．R e m o t e S e n s ．，1 9 9 9 2 0 2 3 6 7 2 3 9 3 ． 图4 反演温度图 F i g ．4 R e t r i e v a lo ft e m p e r a t u r e 逐渐递减的规律不同，这可能跟环渤海地区沿岸城 市和工业区以及土地盐碱化的变化有关．相关验证 需要实地和实测数据，本文不做进一步分析．这种 大面积温度分布规律在高分辨率的T M 影像上体 现不出来，反映大区域的温度分布差异是M O D I S 影像的一个优势． 从图4 可以看出，海陆温度差异非常明显，对 于劈窗算法反演地表温度的精度分析评价，请参见 文献E z o ] ． 4结论 本文简要介绍了劈窗算法，通过对P l a n c k 方 程的辐射强度和温度计算模拟，并对P l a n c k 方程 进行了线性简化，从而简化了劈窗算法的计算过 程．针对M O D I S 波段设置特点，通过建立大气水 汽含量与透过率的关系，然后在反演水汽含量的基 础上，计算第3 1 和3 2 波段的大气透过率．最后运 用劈窗算法反演对了山东半岛的地表温度，结果表 明，简化计算和通过影像反演水汽含量和计算透过 率是可行的． 参考文献 [ 1 ] L iZ ，B e c k e rF ．F e a s i b i l i t yo fl a n ds u r f a c et e m p e r a t u r ea n de m i s s i v i t yd e t e r m i n a t i o nf r o mA V H R R d a t a [ J ] ．R e m o t eS e n s i n go fE n v i r o n m e n t ，1 9 9 3 4 3 [ 2 3 6 7 8 5 ． Q i nZ ，K a r n i e l iA ．P r o g r e s si nt h er e m o t es e n s i n g [ 3 ] [ 4 ] [ 5 ] I - 6 - 1 P r a b h a k a r aC ，D a l uG ．K u n d eVG ．E s t i m a t i o no f s e at e m p e r a t u r ef r o mr e m o t es e n s i n gi nt h el 1t O1 2 ”mw i n d o wr e g i o n [ J ] ．J ．G e o p h y s ．R e s ．，19 7 4 7 9 5 0 3 9 5 0 4 4 ． M O D I Sm a n u a l ．h t t p ／／w w w ．i r s a ．a c ．c n ／d o w n l o a d ／ m o d i s - - b r o c h u r e ．p d f P r i c eJC ．E s t i m a t i n gs u r f a c et e m p e r a t u r ef r o ms a t e i l i t et h e r m a li n f r a r e dd a t a As i m p l ef o r m u l a t i o nf o r t h ea t m o s p h e r i ce f f e c t [ J ] ．R e m o t eS e n s ．E n v i r o n ．， 1 9 8 3 1 3 3 5 3 3 6 1 ． P r i c eJC ．L a n ds u r f a c et e m p e r a t u r em e a s u r e m e n t s f r o mt h e s p l i tw i n d o w c h a n n e l so ft h eN O A A7 a d v a n c e dv e r yh i g hr e s o l u t i o nr a d i o m e t e r [ J ] ．J ． G e o p h y s ．R e s ．，1 9 8 4 8 9 7 2 3 1 7 2 3 7 ． [ 7 ] B e c k e rF ，L iZ ．T o w a r d sa l o c a ls p l i tw i n d o wm e t h o d o v e rl a n ds u r f a c e [ J ] ．I n t ．J ．R e m o t eS e n s ．，1 9 9 0 3 3 6 9 3 9 3 ． [ 8 ] H o l b oH ．R ．J ．C ．L u v a l l ．M o d e l i n gs u r f a c et e m p e r a t u r ed i s t r i b u t i o n si nf o r e s tl a n d s c a p e s [ J ] ，R e m o t e S e n s ．E n v i r o n ．，1 9 8 9 2 7 1 1 2 4 ． [ 9 ] C o o p e r ，DI ，A s r a rG ．E v a l u a t i n ga t m o s p h e r e i cc o r r e c t i o nm o d e l sf o rr e t r i e v i n gs u r f a c et e m p e r a t u r e s f r o mt h eA V H R Ro v e rat a l l g r a s sp r a i r i e s [ J ] ． R e m o t eS e n s ．E n v i r o n ．，1 9 8 9 2 7 9 3 1 0 2 ． [ 1 0 ] O t t l eC ，V i d a l M a d j a rD ．E s t i m a t i o no fl a n ds u r f a c e t e m p e r a t u r ew i t hN O A A 一9D a t a [ J ] ．R e m o t eS e n s ． E n v i r o n ．。1 9 9 2 4 2 7 4 1 ． [ 11 ] P r a t aAJ ．L a n ds u r f a c et e m p e r a t u r ed e r i v e df r o m t h ea d v a n c e dv e r yh i g hr e s o l u t i o nr a d i o m e t e ra n d a l o n g t r a c ks c a n n i n gr a d i o m e t e r [ J ] ．J ．G e o p h y s ． R e s ．，1 9 9 3 9 8 1 6 6 8 9 1 6 7 0 2 ． [ 1 2 ] P r a t aAJ ．L a n ds u r f a c et e m p e r a t u r ef r o mt h ea d v a n c e dv e r yh i g hr e s o l u t i o nr a d i o n m e t e ra n dt h e a l o n g t r a c ks c a n n i n gr a d i o m e t e r ．2 ．E x p e r i m e n t a l r e s u l t sa n dv a l i d a t i o no fA V H R Ra l g o r i t h m s [ J ] ．J ． G e o p h y s ．R a s ．，1 9 9 4 9 9 1 3 0 2 5 1 3 0 5 8 ． 万方数据 5 0中国矿业大学学报 第3 4 卷 f 1 3 ] C o l lCV ，C a s e l l e sS o b r i n oA ，V a l o rE ．O nt h ea t m o s p h e r i cd e p e n d e n c eo ft h es p l i t w i n d o we q u a t i o nf o r l a n ds u r f a c et e m p e r a t u r e [ J ] ．I n t ．J ．R e m o t eS e n s ．， 1 9 9 4 1 5 1 0 5 1 2 2 ． [ 1 4 ] F r a n c aGB 。C r a c k n e l lAP ．R e t r i e v a lo fl a n da n d s e as u r f a c et e m p e r a t u r eu s i n gN O A A 一11A V H R R d a t ai nn o r t h e a s t e r nB r a z i l [ J ] ．I n t ．J ．R e m o t e S e n s ．1 9 9 4 i 5 1 6 9 5 1 7 1 2 ． [ 1 5 ] S e g u i nB ，C o u r a u l tD ，G u e r i fM ．S u r f a c et e m p e r a t u r ea n dv a p o rt r a n s p i r a t i o n a p p l i c a t i o no fl o c a l s c a l em e t h o d st or e g i o n a ls c a l e su s i n gs a t e l l i t ed a t a [ J ] ．R e m o t eS e n s ．E n v i r o n ．，1 9 9 4 4 9 2 8 7 2 9 5 ． [ 1 6 ] C h o u d h u r yJB ，D o r m a nJT ，H s uAY ．M o d e l e d a n do b s e r v e dr e l a t i o n sb e t w e e nt h eA V H R Rs p l i t w i n d o wt e m p e r a t u r ed i f f e r e n c ea n da t m o s p h e r i c p e r c e p t i b l ew a t e ro v e rl a n ds u r f a c e s [ J ] ．R e m o t e S e n s ，E n v i r o n ．，1 9 9 5 5 1 2 8 1 2 9 0 ． [ 173 C a l v e tJC ，J u l l i e nJP ．L a n ds u r f a c et e m p e r a t u r e r e t r i e v a li n d r yc o n d i t i o n sf r o mi n f r a r e da n d E 1 8 ] [ 1 9 3 E 2 0 ] m o c r o w a v es a t e l l i t er a d i o m e t r y [ J ] ．R e m o t eS e n s ． R e v ．，1 9 9 6 1 3 2 3 5 2 5 5 ． Q i nZH ，K a r n i e l iA ．D e r i v a t i o no fs p l i tw i n d o wa l g o r i t h ma n di t ss e n s i t i v i t ya n a l y s i sf o rr e t r i e v i n gl a n d s u r f a c et e m p e r a t u r ef r o mN O A A A V H R Rd a t a [ J ] ． J ．G e o p h y s ．R e s ．，2 0 0 1 1 0 6 2 2 6 5 5 - 2 2 6 7 0 ． 覃志豪，L iwJ ，Z h a n gMH ，等．单窗算法的大气 参数估计方法E J ] ．国土资源遥感，2 0 0 3 5 6 3 7 4 3 ． Q i nZH ，L iWJ ，Z h a n gMH ，e ta 1 ．E s t i m a t i o no f t h ee s s e n t i a l a t m o s p h e r ep a r a m e t e r s o fM o n o w i n d o w a l g o r i t h m f o rl a n ds u r f a c et e m p e r a t u r e r e t r i e v a lf r o ml a n d s a tT M 6 [ J ] ．R e m o t eS e n s i n gf o r L a n d8 LR e s o u r c e s ，2 0 0 3 5 6 3 7 4 3 ． 毛克彪．用于M O D l S 数据的地表温度反演方法研 究[ D ] ．南京大学，2 0 0 4 ． M a oKB ．As t u d yo fm e t h o d sf o rl a n ds u r f a c e t e m p e r a t u r er e t r i e v a lf r o mM O D I Sd a t arD ] ． N a n j i n gU n i v e r s i t y ，2 0 0 4 ． 责任编辑李成俊 第十届中国煤炭分选及加工学术研讨会在中国矿业大学召开 2 0 0 4 年1 2 月1 6 ～1 8 日，第十届中国煤炭分选及加工学术研讨会在中国矿业大学举行．参加会议的有来 自全国各地近百家科研院校和企事业单位的1 5 0 多名专家、学者、工程技术和管理人员．此次研讨会由中国 煤炭学会选煤专业委员会主办，中国矿业大学化工学院承办． 会议特别邀请了清华大学金涌院士、太原理工大学谢克昌院士、中国矿业大学陈清如院士等我国矿 业、能源领域的著名专家参加．众多学者、专业人士到会进行了高层学术论坛，从战略角度研讨了煤炭分选 与加工的发展方向和前景．会议分4 个专题，即煤炭重选、浮选及脱水、选煤设计与管理和煤炭加工．多名专 家学者在大会上作了学术报告和主题发言，出版了全国煤炭分选及加工学术研讨会论文集．与会人员一 致认为，煤炭的分选及加工是煤洁净利用的一种有效途径．这次大会是一次展示我国近年来煤炭分选与加 工领域主要成就的学术盛会，同时也使煤炭行业的企业集团同相关的科研院所加强了联系与合作，从而共 同推动煤炭分选与加工工业的创新与发展． 万方数据