首先當然是要下載Landsat 8衛星影像。打開以下網址
https://eos.com/landviewer/或者https://earthexplorer.usgs.gov/注冊賬戶—>選擇感興趣區域——>選擇傳感器類型如下圖。選擇云量較少的景,保存到EOS storage,然后下載。
下載后，你可以查看有許多波段TIF加上一個元數據，元數據很重要，如下：
接下來我可以使用 Landsat 8 熱波段估算或計算地表溫度（LST）。它只需要通過光柵圖像計算器（ArcMap、ArcGIS Pro、QGIS都可以）應用一組方程即可。主要利用熱波段10和波段4、5計算NVDI。這個例子僅僅計算沒有闡述公式的原理，有些公式有固定的修正值，如需知道原理，請知網搜索文章查看，主要有6個步驟。
1. TOA大氣表觀反射率計算。
TOA (L) = ML * Qcal + ALML = 來自元數據（RADIANCE_MULT_BAND_ x，其中 x 是波段編號）Qcal = 對應于波段10。
AL = 來自元數據的波段特定附加縮放因子（RADIANCE_ADD_BAND_ x，其中 x 是波段編號）
數字代入，得出TOA = 0.0003342 * “波段 10” + 0.1使用 ArcMap 中的柵格計算器工具求解方程。（spatial analysis工具—>地圖代數—>柵格計算器
2.TOA 到亮度溫度BT的轉換
BT = (K 2 / (ln (K 1 / L) + 1)) ? 273.15L = TOAK1 = 來自元數據的特定波段熱轉換常數（K1_CONSTANT_BAND_ x，其中 x 是熱波段編號）。K2 = 來自元數據的特定波段熱轉換常數（K2_CONSTANT_BAND_ x，其中 x 是熱波段編號）。要獲得以攝氏度為單位的結果，需要通過添加絕對零（約 -273.15°C）來調整輻射溫度。BT = (1321.0789 / Ln ((774.8853 / “TOA”) + 1)) – 273.15
3.計算 NDVI
什么是 NDVI（歸一化差異植被指數），不同波段組合的應用。NDVI =（波段 5 – 波段 4）/（波段 5 + 波段 4）NDVI 的計算很重要，因為隨后必須計算與 NDVI高度相關的植被比例 (Pv ) 和與 Pv相關的發射率 (ε) 。NDVI = Float(Band 5 – Band 4) / Float(Band 5 + Band 4)4.計算植被比例 Pv
Pv = Square ((NDVI – NDVI min ) / (NDVI max – NDVI min ))注意：如果NDVImin為負數，那么上面的減號-都改為+號。通常NDVI圖像的最小值和最大值可以直接顯示在圖像中（ArcGIS、QGIS、ENVI、Erdas），否則必須打開柵格的屬性才能得到這些值。有爭議的是植被比例 (Pv) 計算中應當使用土壤和植被 NDVI 值而不是最大和最小 NDVI 值。Pv = Square((“NDVI” – 0.216901) / (0.632267 – 0.216901))
5.計算發射率 ε
ε = 0.004 * Pv + 0.986只需應用光柵計算器中的公式，0.986 的值對應于方程的修正值。6.計算地表溫度
LST = (BT / (1 + (0.00115 * BT / 1.4388) * Ln(ε)))最后應用LST方程得到地表溫度圖。下面就是地表溫度圖，需要注意的是它不等于空氣溫度。
本文鏈接：http://www.aiquka.com/blog/531.html
本文標簽：
