1、基于(jy)TM遙感(yogn)影像的北京地區地表溫度反演杜香霖 2012045001鄧文星(wn xn) 2012045002蔣 偉 2012045004目錄 TOC o 1-3 h z u HYPERLINK l _Toc385484652 1 項目概述 PAGEREF _Toc385484652 h 2 HYPERLINK l _Toc385484653 1.1研究背景 PAGEREF _Toc385484653 h 2 HYPERLINK l _Toc385484654 1.2研究意義 PAGEREF _Toc385484654 h 2 HYPERLINK l _Toc38548465

2、5 1.3研究內容 PAGEREF _Toc385484655 h 2 HYPERLINK l _Toc385484656 1.4研究方法 PAGEREF _Toc385484656 h 2 HYPERLINK l _Toc385484657 2.2研究區概況及數據資料獲取 PAGEREF _Toc385484657 h 4 HYPERLINK l _Toc385484658 2.3遙感數據預處理 PAGEREF _Toc385484658 h 6 HYPERLINK l _Toc385484659 2.3.1輻射定標 PAGEREF _Toc385484659 h 6 HYPERLINK l

3、 _Toc385484660 2.4計算歸一化植被指數NDVI值 PAGEREF _Toc385484660 h 11 HYPERLINK l _Toc385484661 2.5 地表比輻射率計算 PAGEREF _Toc385484661 h 13 HYPERLINK l _Toc385484662 2.6 反演地表溫度 PAGEREF _Toc385484662 h 17 HYPERLINK l _Toc385484663 3. 分析總結 PAGEREF _Toc385484663 h 19 HYPERLINK l _Toc385484664 結果分析: PAGEREF _Toc38548

4、4664 h 191 項目(xingm)概述1.1研究(ynji)背景地球表面溫度是一個(y )重要的水文、氣象參數， 它影響著大氣、海、陸之間的感熱和潛熱交換， 是諸多研究領域不可或缺的基礎資料， 精確定量反演陸面溫度的成果將推動旱災預報和作物缺水研究、農作物產量估算、數值天氣預報、全球氣候變化和全球碳平衡等領域研究的進展。地表溫度在大氣與地表的物質與能量交換， 植被生態，氣候變化等方面有著重要的應用。而遙感信息技術發展以來，以其獨特的優勢在各個領域得到重大應用。利用熱紅外遙感可以在短時間內獲得大范圍、多時相、多分辨率的地表溫度及其溫度場空間分布信息，與傳統的觀測方式相比，具有快速、便捷、觀

5、測范圍大和信息連續等優點。因此，利用熱紅外遙感數據反演地表溫度的應用日趨廣泛， 利用衛星遙感資料進行地表溫度的反演也成為目前遙感定量研究中的重要任務之一。1.2研究意義近年來“城市氣候和環境”問題引起了研究者的廣泛興趣，其中城市熱環境、城市熱島效應及其與環境的相互關系等領域是當前城市氣候與環境研究的中心問題1。地表溫度是很多環境模型的一個重要參數，在大氣與地表的物質與能量交換，天氣預報，氣候變化等方面有著重要的應用。利用熱紅外遙感2可以在短時間內獲得大范圍的地表溫度及其溫度場空間分布信息，與傳統的觀測方式相比，具有速度快、便捷、觀測范圍大和信息連續等優點。本文利用Landsat TM 第6 波

6、段數據反演北京地區地表溫度，分析了北京地區地面溫度分布模式及城市熱島效應，并探討了城市熱島效應與城市土地覆蓋類型及植被指數之間的關系。1.3研究內容利用ENVI軟件，2003年及2009年兩年的Landsat TM 數據的熱紅外數據(第六波段，波長0.4 - 12.5m) 對北京地區的地表溫度夏季分布規律進行了研究，并結合利用遙感數據獲得的土地利用/覆蓋資料和社會經濟統計資料分析了各區溫度差異的原因與變化規律，為北京市的城市建設和土地合理利用規劃提供一定的借鑒。1.4研究(ynji)方法實際操作中，以LandsatTM數據作為主要(zhyo)數據源，首先網上下載符合操作要求的北京地區兩年TM遙

7、感(yogn)影像，在ENVI 5.0下對遙感影像進行圖像增強，輻射定標，大氣矯正，幾何校正，接著在ENVI軟件中測定影像的比輻射率，然后經過一系列的處理，得到地表溫度，并進行結論分析，得到結果。操作部湊如下:第一步對遙感影像進行圖像增強，輻射定標、大氣矯正，幾何校正。第二步，圖像計算NDVI值，運用ENVI軟件對遙感影像測定比輻射率。第三步，計算相同溫度下黑體的輻射亮度值。第四步，反演地表溫度。第五步，將兩幅圖像進行對比，得出結論。2.影像處理過程2.1遙感影像處理軟件簡介ENVI (The Environment for Visualizing Images)是美國ITT Visual I

8、nformation Solutions公司的旗艦產品。ENVI由遙感領域的科學家采用IDL開發的一套功能強大的遙感圖像處理軟件；它是快速、便捷、準確地從地理空間影像中提取信息的首屈一指的軟件解決方案，它提供先進的，人性化的使用工具來方便用戶讀取、準備、 探測、分析和共享影像中的信息。今天，眾多的影像分析師和科學家選擇ENVI來從地理空間影像中提取信息。已經廣泛應用于科研、環境保護、氣象、石油礦產勘探、農業、林業、醫學、國防安全、地球科學、公用設施管理、遙感工程、水利、海洋，測繪勘察和城市與區域規劃等行業。ENVI包含齊全的 HYPERLINK /z/Search.e?sp=S遙感影像&ch=

9、w.search.yjjlink&cid=w.search.yjjlink 遙感影像處理功能：常規處理、 HYPERLINK /z/Search.e?sp=S幾何校正&ch=w.search.yjjlink&cid=w.search.yjjlink 幾何校正、定標、多 HYPERLINK /z/Search.e?sp=S光譜分析&ch=w.search.yjjlink&cid=w.search.yjjlink 光譜分析、高光譜分析、雷達分析、 HYPERLINK /z/Search.e?sp=S地形地貌&ch=w.search.yjjlink&cid=w.search.yjjlink 地形地貌

10、分析、 HYPERLINK /z/Search.e?sp=S矢量&ch=w.search.yjjlink&cid=w.search.yjjlink 矢量應用、神經網絡分析、區域分析、GPS聯接、正射影象圖生成、 HYPERLINK /z/Search.e?sp=S三維圖像&ch=w.search.yjjlink&cid=w.search.yjjlink 三維圖像生成、豐富的可供二次開發調用的函數庫、制圖、數據輸入/輸出等功能組成了 HYPERLINK /z/Search.e?sp=S圖像處理軟件&ch=w.search.yjjlink&cid=w.search.yjjlink 圖像處理軟件中非

11、常全面的系統。ENVI對于要處理的圖像波段數沒有限制，可以處理最先進的衛星格式，如Landsat7、IKONOS、SPOT， RADARSAT， NASA， NOAA， EROS和TERRA，并準備接受未來所有傳感器的信息。主界面圖如圖1。圖1 ENVI主界面(jimin)圖2.2研究區概況(gikung)及數據資料獲取北京地處華北平原北端，位于(wiy)北緯3926- 414，東經 11524-11831 之間。 位于華北平原西北邊緣。毗鄰渤海灣，上靠遼東半島，下臨山東半島。北京與天津相鄰，并與天津一起被河北省環繞。西部是太行山山脈余脈的西山，北部是燕山山脈的軍都山，兩山在南口關溝相交，形成

12、一個向東南展開的半圓形大山彎，人們稱之為“北京彎”，它所圍繞的小平原即為北京小平原。全市土地面積16410 平方公里， 其中平原面積6338 平方公里 ，占38.6%。屬于典型的溫帶大陸性季風氣候，冬季寒冷干燥，夏季高溫多雨。 在該區利用遙感信息進行區域性尺度的地表溫度及時空要素監測非常重要，有實際指導意義。本文(bnwn)選用Landsat TM的影像(yn xin)數據，TM影像(yn xin)包含7個波段，各個波段的特征和用途及遙感數據詳細信息如下所示: TM 影像各波段特征及主要用途波段號波段類型波長空間分辨率功能與主要用途1藍色0.45-0.5230對水體有穿透能力，用來分析土地利用

13、、干燥的土壤、植被特征、編制森林分布圖2綠色0.52-0.6030對水體的穿透能力較強，對植被的反射敏感，位于葉綠素的兩個吸收帶之間，利用這一波段增強鑒別植被的能力；能區分林型、樹種3紅色0.63-0.6930位于葉綠素的吸收區，能增強植被覆蓋與無植被之間的反差，亦能增強同類植被的反差；可根據植被的色調判斷植物的健狀況4近紅外0.76-0.9030集中反映植物的強反射，用于植被類型、生物量和作物長勢的調查，繪制水體邊界和土壤濕度，也可用來增強土壤與農作物和陸地與水域之間的反差5中紅外1.55-1.7530處于水的吸收帶，對含水量反應敏感，可用于土壤濕度、植物含水量調查、水分狀況研究、作物長勢分

14、析6熱紅外10.4-12.5120屬于熱紅外波段，對熱異常敏感。監測與人類活動有關的熱特征，用于熱分布制圖、巖石識別和地質探礦、水體溫度變化制圖7中紅外2.08-2.3530探測高溫輻射源，如監測森里火災、火山活動等，區分巖石類型、地質探礦與制圖 數據源詳細信息 數據日期 分辨率（m） 投影2003-08-25 30 UTM，WGS84 Zone502009-09-12 30 UTM，WGS84 Zone50遙感圖像(t xin)下載過程2.3遙感(yogn)數據預處理由于地球(dqi)自轉、大氣折射、地形起伏、傳感器姿態變化等多種因素的影響，使遙感圖像產生輻射失真和幾何變形，這些失真與畸變會

15、影響影像的質量和應用，必須進行消除。2.3.1輻射定標輻射定標，就是將像素灰度值 DN( Digital Number) 轉換為輻射亮度或者反射率。通過定標，可以消除一些大氣的干擾。本文通過ENVI軟件中自帶的對TM的定標工具，計算輻射亮度或反射率。實現方法： File-Open External File-Landsat-GeoTIFF with MetadataBasicTools-Preprocessing-Calibration Utilitties-Landsat Calibration利用(lyng)ENVI，在輸入(shr)欄中輸入太陽高度角，成像日期等參數(cnsh)。如圖3.

16、定標對話框輻射定標后的圖像如圖4。 影像定標(dn bio)結果2.3.2ENVI 大氣(dq)矯正大氣校正的目的是消除大氣和光照等因素對第五反射的影響，獲得地物反射率、輻射率、地表(dbio)溫度等真實物理模型參數，包括消除大氣中水蒸氣、氧氣、二氧化碳、甲烷和臭氧等對地物反射的影響；消除大氣分子和氣溶膠散射的影響。 第一步：輻射定標第二步：輻射亮度單位轉換 由于輻射定標之后(zhhu)得到的輻射亮度單位為W/(m2msr)， FLAASH要求(yoqi)的輻射亮度單位為W/(cm2nmsr)，兩者相差(xin ch)10倍，并且FLAASH模塊要求BIL的圖像存儲方式，因此需要作出調整。 第

17、三步：輸入FLAASH參數。 FLAASH參數輸入單位調整： Basic Tools-Band Math表達式為：b1/10圖像存儲方式： Basic Tools-Convert Data(BSQ、BIL、BIP) FLAASH模塊： 根據(gnj)圖像元數據輸入參數，高級設置為默認，執行模塊，輸出結果。大氣(dq)校正結果展示2.3.3圖像的幾何(j h)校正通過一系列的數學模型來改正和消除遙感影像成像時因攝影材料變形、物鏡畸變、大氣折光、地球曲率、地球自轉、地形起伏等因素導致的原始圖像上各地物的幾何位置、形狀、尺寸、方位等特征與在參照系統中的表達要求不一致時產生的變形。幾何校正步奏如下所示

18、：2.4計算(j sun)歸一化植被指數NDVI值在假設地表為朗伯體的前提下，用求得的輻射亮度計算表觀反射率。表觀反射率就是指大氣層頂的反射率，輻射定標的(bio de)結果之一，大氣層頂表觀反射率，簡稱表觀反射率，又稱視反射率。分別(fnbi)對2003年以及2009年北京遙感影響做NDVI歸一化植被指數計算，如下所示：2003年北京地區NDVI顯示(xinsh) 2009年北京(bi jn)NDVI植被(zhbi)顯示對比2003年與2009年北京同一地區植被指數：2.5 地表(dbio)比輻射率計算物體(wt)的比輻射率是物體向外輻射電磁波的能力表征。它不僅依賴于地表物體的組成，而且與物

19、體的表面狀態(表面(biomin)粗糙度等)及物理性質(介電常數、含水量等)有關，并隨著所測定的波長和觀測角度等因素有關。在大尺度上對比輻射率精確測量的難度很大，目前只是基于某些假設獲得比輻射率的相對值，本文主要根據可見光和近紅外光譜信息來估計比輻射率。2.5.1植被覆蓋度計算計算植被覆蓋度Fv 采用的是混合像元分解法，將整景影像的地類大致分為水體、植被和建筑，具體的計算公式如下：FV = (NDVI- NDVIS)(NDVIV - NDVIS) (10)其中，NDVI 為歸一化差異植被指數，取NDVIV = 0.70 和NDVIS = 0.00，且有，當某個像元的NDVI 大于0.70 時，

20、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 值，如圖6。圖6 計算(j sun)植被覆蓋度對話框得到(d do)植被蓋度圖像，如圖7。 2003年植被(zhbi)覆蓋度2009年植被(zhbi)覆蓋度2.5.2地表(dbio)比輻射率計算根據前人的研究，將遙感影像分為水體、城鎮和自然(zrn)表面3 種類型。本專題采取以下方法計算研究區地表比輻射

21、率：水體像元的比輻射率賦值為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-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 計算(j sun)比輻射率對話框地表(dbio)比輻射率計算利用(lyng)ENVI 主菜單(ci dn)-Basic Tools-Band Math，在公式輸入欄中輸入：（b2-3.39-0.6*（1-b1）*5.12）/（0.6*b1）b1： 地表比輻射率值；b2：表示熱紅外波段大氣校正后的輻射定標值。得到了溫度為T 的黑體在熱紅外波段的輻射亮度值。如圖