1、基于TM遙感影像的北京地區地表溫度反演杜香霖 2012045001鄧文星 2012045002蔣 偉 2012045004目錄1 項目概述21.1研究背景21.2研究意義21.3研究內容21.4研究方法22.2研究區概況及數據資料獲取42.3遙感數據預處理62.3.1輻射定標62.4計算歸一化植被指數NDVI值112.5 地表比輻射率計算132.6 反演地表溫度173. 分析總結19結果分析:191 項目概述1.1研究背景地球表面溫度是一個重要的水文、氣象參數， 它影響著大氣、海、陸之間的感熱和潛熱交換， 是諸多研究領域不可或缺的基礎資料， 精確定量反演陸面溫度的成果將推動旱災預報和作物缺水研

2、究、農作物產量估算、數值天氣預報、全球氣候變化和全球碳平衡等領域研究的進展。地表溫度在大氣與地表的物質與能量交換， 植被生態，氣候變化等方面有著重要的應用。而遙感信息技術發展以來，以其獨特的優勢在各個領域得到重大應用。利用熱紅外遙感可以在短時間內獲得大范圍、多時相、多分辨率的地表溫度及其溫度場空間分布信息，與傳統的觀測方式相比，具有快速、便捷、觀測范圍大和信息連續等優點。因此，利用熱紅外遙感數據反演地表溫度的應用日趨廣泛， 利用衛星遙感資料進行地表溫度的反演也成為目前遙感定量研究中的重要任務之一。1.2研究意義近年來“城市氣候和環境”問題引起了研究者的廣泛興趣，其中城市熱環境、城市熱島效應及其

3、與環境的相互關系等領域是當前城市氣候與環境研究的中心問題1。地表溫度是很多環境模型的一個重要參數，在大氣與地表的物質與能量交換，天氣預報，氣候變化等方面有著重要的應用。利用熱紅外遙感2可以在短時間內獲得大范圍的地表溫度及其溫度場空間分布信息，與傳統的觀測方式相比，具有速度快、便捷、觀測范圍大和信息連續等優點。本文利用Landsat TM 第6 波段數據反演北京地區地表溫度，分析了北京地區地面溫度分布模式及城市熱島效應，并探討了城市熱島效應與城市土地覆蓋類型及植被指數之間的關系。1.3研究內容利用ENVI軟件，2003年及2009年兩年的Landsat TM 數據的熱紅外數據(第六波段，波長0.

4、4 - 12.5m) 對北京地區的地表溫度夏季分布規律進行了研究，并結合利用遙感數據獲得的土地利用/覆蓋資料和社會經濟統計資料分析了各區溫度差異的原因與變化規律，為北京市的城市建設和土地合理利用規劃提供一定的借鑒。1.4研究方法實際操作中，以LandsatTM數據作為主要數據源，首先網上下載符合操作要求的北京地區兩年TM遙感影像，在ENVI 5.0下對遙感影像進行圖像增強，輻射定標，大氣矯正，幾何校正，接著在ENVI軟件中測定影像的比輻射率，然后經過一系列的處理，得到地表溫度，并進行結論分析，得到結果。操作部湊如下:第一步對遙感影像進行圖像增強，輻射定標、大氣矯正，幾何校正。第二步，圖像計算N

5、DVI值，運用ENVI軟件對遙感影像測定比輻射率。第三步，計算相同溫度下黑體的輻射亮度值。第四步，反演地表溫度。第五步，將兩幅圖像進行對比，得出結論。2.影像處理過程2.1遙感影像處理軟件簡介ENVI (The Environment for Visualizing Images)是美國ITT Visual Information Solutions公司的旗艦產品。ENVI由遙感領域的科學家采用IDL開發的一套功能強大的遙感圖像處理軟件；它是快速、便捷、準確地從地理空間影像中提取信息的首屈一指的軟件解決方案，它提供先進的，人性化的使用工具來方便用戶讀取、準備、 探測、分析和共享影像中的信息。今

6、天，眾多的影像分析師和科學家選擇ENVI來從地理空間影像中提取信息。已經廣泛應用于科研、環境保護、氣象、石油礦產勘探、農業、林業、醫學、國防安全、地球科學、公用設施管理、遙感工程、水利、海洋，測繪勘察和城市與區域規劃等行業。ENVI包含齊全的遙感影像處理功能：常規處理、幾何校正、定標、多光譜分析、高光譜分析、雷達分析、地形地貌分析、矢量應用、神經網絡分析、區域分析、GPS聯接、正射影象圖生成、三維圖像生成、豐富的可供二次開發調用的函數庫、制圖、數據輸入/輸出等功能組成了圖像處理軟件中非常全面的系統。ENVI對于要處理的圖像波段數沒有限制，可以處理最先進的衛星格式，如Landsat7、IKONO

7、S、SPOT， RADARSAT， NASA， NOAA， EROS和TERRA，并準備接受未來所有傳感器的信息。主界面圖如圖1。圖1 ENVI主界面圖2.2研究區概況及數據資料獲取北京地處華北平原北端，位于北緯39°26'- 41°4'，東經 115°24'-118°31' 之間。 位于華北平原西北邊緣。毗鄰渤海灣，上靠遼東半島，下臨山東半島。北京與天津相鄰，并與天津一起被河北省環繞。西部是太行山山脈余脈的西山，北部是燕山山脈的軍都山，兩山在南口關溝相交，形成一個向東南展開的半圓形大山彎，人們稱之為“北京彎”，它所圍繞的

8、小平原即為北京小平原。全市土地面積16410 平方公里， 其中平原面積6338 平方公里 ，占38.6%。屬于典型的溫帶大陸性季風氣候，冬季寒冷干燥，夏季高溫多雨。 在該區利用遙感信息進行區域性尺度的地表溫度及時空要素監測非常重要，有實際指導意義。本文選用Landsat TM的影像數據，TM影像包含7個波段，各個波段的特征和用途及遙感數據詳細信息如下所示: TM 影像各波段特征及主要用途波段號波段類型波長空間分辨率功能與主要用途1藍色0.45-0.5230對水體有穿透能力，用來分析土地利用、干燥的土壤、植被特征、編制森林分布圖2綠色0.52-0.6030對水體的穿透能力較強，對植被的反射敏感，

9、位于葉綠素的兩個吸收帶之間，利用這一波段增強鑒別植被的能力；能區分林型、樹種3紅色0.63-0.6930位于葉綠素的吸收區，能增強植被覆蓋與無植被之間的反差，亦能增強同類植被的反差；可根據植被的色調判斷植物的健狀況4近紅外0.76-0.9030集中反映植物的強反射，用于植被類型、生物量和作物長勢的調查，繪制水體邊界和土壤濕度，也可用來增強土壤與農作物和陸地與水域之間的反差5中紅外1.55-1.7530處于水的吸收帶，對含水量反應敏感，可用于土壤濕度、植物含水量調查、水分狀況研究、作物長勢分析6熱紅外10.4-12.5120屬于熱紅外波段，對熱異常敏感。監測與人類活動有關的熱特征，用于熱分布制圖

10、、巖石識別和地質探礦、水體溫度變化制圖7中紅外2.08-2.3530探測高溫輻射源，如監測森里火災、火山活動等，區分巖石類型、地質探礦與制圖 數據源詳細信息 數據日期 分辨率（m） 投影2003-08-25 30 UTM，WGS84 Zone502009-09-12 30 UTM，WGS84 Zone50遙感圖像下載過程2.3遙感數據預處理由于地球自轉、大氣折射、地形起伏、傳感器姿態變化等多種因素的影響，使遙感圖像產生輻射失真和幾何變形，這些失真與畸變會影響影像的質量和應用，必須進行消除。2.3.1輻射定標輻射定標，就是將像素灰度值 DN( Digital Number) 轉換為輻射亮度或者反

11、射率。通過定標，可以消除一些大氣的干擾。本文通過ENVI軟件中自帶的對TM的定標工具，計算輻射亮度或反射率。實現方法： File-Open External File-Landsat-GeoTIFF with MetadataBasicTools-Preprocessing-Calibration Utilitties-Landsat Calibration利用ENVI，在輸入欄中輸入太陽高度角，成像日期等參數。如圖3.定標對話框輻射定標后的圖像如圖4。 影像定標結果2.3.2ENVI 大氣矯正大氣校正的目的是消除大氣和光照等因素對第五反射的影響，獲得地物反射率、輻射率、地表溫度等真實物理模型

12、參數，包括消除大氣中水蒸氣、氧氣、二氧化碳、甲烷和臭氧等對地物反射的影響；消除大氣分子和氣溶膠散射的影響。 第一步：輻射定標 第二步：輻射亮度單位轉換 由于輻射定標之后得到的輻射亮度單位為W/(m2·µm·sr)， FLAASH要求的輻射亮度單位為µW/(cm2·nm·sr)，兩者相差10倍，并且FLAASH模塊要求BIL的圖像存儲方式，因此需要作出調整。 第三步：輸入FLAASH參數。 FLAASH參數輸入 單位調整： Basic Tools-Band Math 表達式為：b1/10 圖像存儲方式： Basic Tools-Conv

13、ert Data(BSQ、BIL、BIP) FLAASH模塊： 根據圖像元數據輸入參數，高級設置為默認，執行模塊，輸出結果。大氣校正結果展示2.3.3圖像的幾何校正通過一系列的數學模型來改正和消除遙感影像成像時因攝影材料變形、物鏡畸變、大氣折光、地球曲率、地球自轉、地形起伏等因素導致的原始圖像上各地物的幾何位置、形狀、尺寸、方位等特征與在參照系統中的表達要求不一致時產生的變形。幾何校正步奏如下所示：2.4計算歸一化植被指數NDVI值在假設地表為朗伯體的前提下，用求得的輻射亮度計算表觀反射率。表觀反射率就是指大氣層頂的反射率，輻射定標的結果之一，大氣層頂表觀反射率，簡稱表觀反射率，又稱視反射率。

14、分別對2003年以及2009年北京遙感影響做NDVI歸一化植被指數計算，如下所示：2003年北京地區NDVI顯示 2009年北京NDVI植被顯示對比2003年與2009年北京同一地區植被指數：2.5 地表比輻射率計算物體的比輻射率是物體向外輻射電磁波的能力表征。它不僅依賴于地表物體的組成，而且與物體的表面狀態(表面粗糙度等)及物理性質(介電常數、含水量等)有關，并隨著所測定的波長和觀測角度等因素有關。在大尺度上對比輻射率精確測量的難度很大，目前只是基于某些假設獲得比輻射率的相對值，本文主要根據可見光和近紅外光譜信息來估計比輻射率。2.5.1植被覆蓋度計算計算植被覆蓋度Fv 采用的是混合像元分解

15、法，將整景影像的地類大致分為水體、植被和建筑，具體的計算公式如下：FV = (NDVI- NDVIS)(NDVIV - NDVIS) (10)其中，NDVI 為歸一化差異植被指數，取NDVIV = 0.70 和NDVIS = 0.00，且有，當某個像元的NDVI 大于0.70 時，FV取值為1；當NDVI 小于0.23，FV取值為0。利用ENVI 主菜單->Basic Tools->Band Math，在公式輸入欄中輸入：（b1 gt 0.7）*1+（b1 lt 0）*0+（b1 ge 0 and b1 le 0.7）*（(b1-0.0)/(0.7-0.0)）b1：NDVI 值，如

16、圖6。圖6 計算植被覆蓋度對話框得到植被蓋度圖像，如圖7。 2003年植被覆蓋度2009年植被覆蓋度2.5.2地表比輻射率計算根據前人的研究，將遙感影像分為水體、城鎮和自然表面3 種類型。本專題采取以下方法計算研究區地表比輻射率：水體像元的比輻射率賦值為0.995，自然表面和城鎮像元的比輻射率估算則分別根據下式（3）（4）進行計算：G1= 0.9625 + 0.0614FV - 0.0461FV2 （11） G2 = 0.9589 + 0.086FV - 0.0671FV2 （12）式中，G1和G2 分別代表自然表面像元和城鎮像元的比輻射率。利用ENVI 主菜單->Basic Tools

17、->Band Math，在公式輸入欄中輸入：（b1 le 0）*0.995+（b1 gt 0 and b1 lt 0.7）*(0.9589 + 0.086*b2 - 0.0671*b22)+(b1 ge0.7)*(0.9625 + 0.0614*b2 - 0.0461*b22)，如圖8。b1：NDVI 值；b2：植被覆蓋度值。圖8 計算比輻射率對話框地表比輻射率計算利用ENVI 主菜單->Basic Tools->Band Math，在公式輸入欄中輸入：（b2-3.39-0.6*（1-b1）*5.12）/（0.6*b1）b1： 地表比輻射率值；b2：表示熱紅外波段大氣校正后的

18、輻射定標值。得到了溫度為T 的黑體在熱紅外波段的輻射亮度值。如圖10。 2003年亮溫圖2009亮溫圖2.6 反演地表溫度在獲取溫度為TS 的黑體在熱紅外波段的輻射亮度后，根據普朗克公式的反函數，求得地表真實溫度TS：TS = K2/ln(K1/ B(TS)+ 1) （15）對TM，K1 =666.09W/(m2srm)，K2 =1282.71K。利用ENVI 主菜單->Basic Tools->Band Math，在公式輸入欄中輸入：（1282.71）/alog（666.09/b1 +1）-273b1：溫度為T的黑體在熱紅外波段的輻射亮度值。得到真實的地表溫度值，單位是攝氏度。如圖11。圖11 地表溫度反演對話框結果顯示：3. 分析總結 結果分析:利用 ENVI 彩色制圖工具處理地溫反演