基于成對hj-1b影像時空信息的地表溫度反演方法
【專利摘要】本發(fā)明公開了一種基于成對HJ-1B影像時空信息的地表溫度反演方法。該方法是針對現有的單窗（或單通道）算法在應用中所存在的參數獲取困難而提出的一種新策略。實現該方法的步驟為：S1：獲取研究區(qū)兩個成對時相的HJ-1B影像，并進行相應的預處理；S2：對預處理后的影像以對稱的方式進行邊界延拓處理；S3：識別相似像元及其數目，確定相應的反演方法；S4：依據S3的判定結果，構建輻射傳輸方程組，并求解地表輻射值；S5：循環(huán)進行S3~S4,得的整個研究區(qū)在兩個時相的地表輻射值，進行異常值識別和處理；S6：由處理后的地表輻射值反演地表溫度。本發(fā)明能在相關參數無法準確獲取的情況下準確反演地表溫度，將有效促進HJ-1B熱紅外數據的定量應用。
【專利說明】基于成對HJ-1B影像時空信息的地表溫度反演方法

【技術領域】
[0001] 本發(fā)明涉及熱紅外遙感【技術領域】，更具體的是在地表比輻射率和大氣參數均未知 的情況下，利用兩個成對時相HJ-IB影像的時空信息，從單一熱紅外波段的輻射亮度直接 定量反演地表溫度的方法。

【背景技術】
[0002] 地表溫度是影響生態(tài)系統的關鍵因子之一，它與許多地表過程有關。熱紅外衛(wèi)星 遙感技術被認為是現階段獲取區(qū)域和全球尺度地表溫度信息的一個可行且有效的手段。為 了實現從熱紅外衛(wèi)星遙感影像準確提取地表溫度信息，國內外研究人員提出了一系列的反 演方法，如單窗(或單通道)算法、劈窗算法、溫度/比輻射率分離算法、白天/晝夜法。與其 它方法相比，單窗(或單通道）算法的優(yōu)勢在于它能從單一熱紅外波段影像反演地表溫度。 因此，在實際應用中，對于Landsat TM/ETM+，CBERS，HJ-IB等僅具有單一熱紅外波段的傳 感器，單窗(或單通道)算法是唯一的選擇。由輻射傳輸方程可知，為了從單一熱紅外波段影 像準確地反演地表溫度，需要同時對大氣影響和地表比輻射率進行校正。已有的研究通過 大量的模擬實驗和簡化假設提出了相應的經驗算法模型及其改進形式（Qin和Karnieli， 2001 ; Jimenez-Munoz 和 Sobrino, 2003;張勇等，2006;段四波等，2008 ;Jimenez_Munoz 等，2009;Cristobal 等，2009)。
[0003] 然而,這些經驗的算法均以地表比輻射率和大氣參數(如大氣平均溫度、大氣水分 含量）已知為前提。目前，對地表比輻射率的估算主要采用地表覆蓋分類賦值法（Synder 等，1998)和植被指數法（Valor 和 Caselles，1996 ;覃志豪等，2004 ;Sobrino 等，2008)。 基于地表覆蓋分類的方法需要預先知道不同地類的比輻射率，并認為每個類別的比輻射率 是固定不變的，這顯然是不合理的；盡管基于植被指數的方法考慮了像元的組分差異，但依 然需要預先確定代表性組分的比輻射率，如土壤，植被等。另外，受準確大氣參數獲取條件 的限制，很多研究多假定整個研究區(qū)域的大氣狀態(tài)是均一分布的，在地表溫度反演過程中 不考慮大氣的影響，如Artis和Carnahan (1982)在早期所提出的一個單窗算法就沒有考 慮大氣的影響，因此利用該算法不能得到真實的地表溫度。由于大氣狀態(tài)往往存在一定程 度的空間異質性，假定大氣狀況均質，不進行大氣校正，或以某一點代替整個區(qū)域必將給遙 感地表溫度反演帶來誤差（陳峰等，2010)。理論上，可以從MODIS的大氣產品（如MOD 05) 提取現有經驗算法所需的大氣參數，但是該產品自身的不確定性可能會在有些情況下給最 終的地表溫度結果帶來不利影響?？梢姡捎谠跍蚀_獲取相關參數方面存在的困難，將在一 定程度上限制現有的單窗(或單通道）算法實施，從而也將影響HJ-IB單一熱紅外波段數據 的定量應用。
[0004] 本發(fā)明針對HJ-IB星的高回訪頻率的特點，利用兩個成對時相HJ-IB影像的時空 信息，基于兩個合理的假設，提出一種直接從單一熱紅外波段影像反演地表溫度的方法，而 無需過度依賴準確的地表比輻射率和大氣參數信息。因此，該新方法可以視為對現有的單 窗(或單通道）算法的有益補充。


【發(fā)明內容】

[0005] 要解決的技術問題本發(fā)明主要解決的技術問題是：在無法準確獲取研究區(qū)地表 比輻射率和大氣參數(如大氣水分含量，大氣平均作用溫度）的情況下，即無法有效實施現 有的單窗(或單通道）算法時，實現從HJ-IB影像準確反演地表溫度。
[0006] 技術方案為了解決上述技術問題，本發(fā)明提出了一種基于成對HJ-IB影像時空 信息的地表溫度反演方法，包括以下步驟： Sl獲取研究區(qū)兩個時相的HJ-IB影像(含C⑶多光譜和IRS熱紅外影像)，結合官方公 布資料和影像頭文件信息，相應地，計算C⑶影像的各波段的星上反射率值和IRS熱紅外 影像輻射亮度值；并對HJ-IB C⑶反射率影像和熱紅外波段輻射影像進行地理校正等預處 理； S2對預處理后的影像以對稱的方式進行邊界延拓處理； S3依據每個目標像元所在的滑動窗口范圍內各像元的光譜相似度指數，識別相似像 元及其數目，確定相應的反演方法； S4依據S3的結果，相應地提取相似或準相似像元的輻射值，構建一個簡化的輻射傳 輸方程組，通過一個帶約束的非線性優(yōu)化模型實現對該方程組的求解，得到目標像元在兩 個成對時相的地表輻射值； S5對延拓后的影像循環(huán)進行S3~S4，得的整個研究區(qū)兩個時相的地表輻射值，并對其 進行異常值識別和處理； S6利用普朗克反方程或其簡化模型從處理后的地表輻射值反演地表溫度，并同時保 存反演方法的標識信息。
[0007] 有益效果本發(fā)明可以在地表比輻射率和大氣參數均無法準確獲取的情況下，僅 利用兩個成對時相HJ-IB影像的時空信息實現對地表溫度的準確反演，克服了現有的方法 對輔助信息過于依賴的問題，進而能在一定程度上促進單一熱紅外波段影像的定量應用。

【專利附圖】

【附圖說明】
[0008] 圖1是一種基于成對時相HJ-IB影像時空信息的地表溫度反演方法示意圖。
[0009] 圖2是基于本發(fā)明的兩個時相HJ-IB地表溫度反演結果圖。

【具體實施方式】
[0010] 以下結合附圖和具體實施例，對本發(fā)明提供的基于成對HJ-IB影像時空信息的地 表溫度反演方法進行詳細描述。以下實施例僅用于說明本發(fā)明，而不是用來限制本發(fā)明的 范圍。如圖1所示，本發(fā)明基于成對HJ-IB影像時空信息的地表溫度反演方法，主要包括以 下步驟： Sl獲取研究區(qū)影像，并進行輻射校正、地理校正、重采樣和裁剪預處理 在本實施例中，獲取時相相近的兩景HJ-IB IRS影像(獲取時相分別為2011年12月 18日和2011年12月22 (path/row :452/91)，后文分別記為時相T18和T22)，并以珠江三 角洲的主要區(qū)域為實驗區(qū)。由于數據質量的影響，僅有與T18相同時相的一景HJ-IB CCD2 多光譜影像（path/row: 454/90)是可用的(均來自http://www.cresda.cn)。因此，本實 施例中假定兩個時相的地表狀態(tài)相同，則可用T18時相的CCD影像充當T22時相的影像，而 不對最終的結果產生影響。一景Landsat5 TM影像（獲取時相為2009年1月2日，path/ row: 122/44)，作為地理校正的參考影像。將HJ-IB CXD影像和IRS熱紅外波段的DN值轉 換為輻射值，并進一步將C⑶輻射值影像轉換為星上反射率值。以Landsat5 TM為參考影 像，在ENVI軟件平臺上，通過人工目視解譯和對比判別，在Landsat5 TM影像、以上各HJ-IB CCD星上反射率和IRS熱紅外輻射影像選取15個均勻分布的特征同名點，并采用二次多項 式模型對HJ-IB CXD星上反射率和IRS熱紅外輻射影像進行地理校正。由于HJ-IB C⑶影 像和IRS熱紅外影像存在空間分辨率差異，這里將HJ-IB C⑶星上反射率（30mX30m)結果 重采樣到300mX 300m，并裁剪確定最終的研究區(qū)影像； S2對預處理后的影像進行邊界延拓處理 為了得到研究區(qū)完整的反演結果，采用對稱的方式對Sl中預處理后的CCD星上反射率 影像和IRS熱紅外輻射亮度值影像均進行邊界延拓處理，在本實施例中延拓寬度取值為3 個像元； S3識別相似像元及其數目，確定相應的反演方法 依據每個目標像元所在的滑動窗口范圍(本實施例中定義為7X7像元的范圍）內各像 元的光譜相似度指數，識別相似像元及其數目，確定相應的反演方法。本發(fā)明中采用RMSD 指數（Chen等，2011)來度量某個像元與目標像元的光譜相似度，RMSD指數通過下式計算：

【權利要求】
1. 一種基于成對HJ-1B影像時空信息的地表溫度反演方法，其特征在于包括以下步 驟： S1:獲取研究區(qū)兩個成對時相的HJ-1B影像(含(XD多光譜和IRS熱紅外影像)，結合 官方公布資料和影像頭文件信息，相應地，計算CCD影像的各波段的星上反射率值和IRS熱 紅外影像輻射亮度值；并對HJ-1B C⑶反射率影像和熱紅外波段輻射影像進行地理校正等 預處理； S2:對預處理后的影像以對稱的方式進行邊界延拓處理； S3:依據每個目標像元所在的滑動窗口范圍內各像元的光譜相似度指數，識別相似像 元及其數目，確定相應的反演方法； S4:依據S3的結果，相應地提取相似或準相似像元的輻射值，構建一個簡化的輻射傳 輸方程組，通過一個帶約束的非線性優(yōu)化模型實現對方程組的求解，得到目標像元在兩個 時相的地表福射估計值； S5:對延拓后的影像循環(huán)進行S3飛4，得的整個研究區(qū)兩個時相的地表輻射值，并對 其進行異常值識別和處理； S6:利用普朗克反方程或其簡化模型從處理后的地表輻射值反演地表溫度，并同時保 存與反演方法相關的標識信息。
2. 如權利要求1所述的基于成對HJ-1B影像時空信息的地表溫度反演方法，其特征 在于利用了兩個成對HJ-1B影像的時空信息，輸入影像包括：HJ-1B C⑶多光譜星上反射率 影像和HJ-1B IRS熱紅外輻射亮度值影像，所述步驟S1中需分別將C⑶多光譜和IRS熱 紅外的原始DN值影像分別轉換為星上反射率和輻射亮度，并進行地理校正、C⑶影像重采 樣和研究區(qū)裁剪處理；由于目前HJ-1B衛(wèi)星在地理定位方面存在的偏差問題，因此采用與 HJ-1B影像相近時相的Landsat TM/ ETM+作為參考影像，進一步對HJ-1B C⑶星上反射率 影像和IRS熱紅外波段輻射亮度值影像進行地理校正；為了解決HJ-1B C⑶影像和IRS熱 紅外影像在空間分辨率方面存在的差異，這里將HJ-1B C⑶星上反射率（30mX 30m)結果重 采樣到300mX 300m，與熱紅外影像同空間分辨率，并裁剪影像確定最終的研究區(qū)影像。
3. 如權利要求1所述的基于成對HJ-1B影像時空信息的地表溫度反演方法，其特征 在于，所述步驟S2中，為了得到研究區(qū)完整的反演結果，采用對稱的方式對S1中預處理后 的C⑶星上反射率影像和IRS熱紅外輻射亮度值影像均進行邊界延拓處理，延拓寬度為像 元大小(為反演方法中所采用滑動窗口的大小，如3, 5, 7,9,…;鑒于HJ-1B IRS熱紅外影 像的300mX 300m的空間分辨率，結合我們以往的研究，我們建議取7個像元，這樣延拓寬度 為3個像元)。
4. 如權利要求1所述的基于成對時相HJ-1B影像時空信息的地表溫度反演方法，其 特征在于，所述步驟S3中，定義了滑動窗口范圍內的相似或準相似像元，并認為它們具有 相同的光譜和熱輻射特性，進而在每個像元所在的滑動窗口內選取合理的像元參與對目標 像元地表輻射值的優(yōu)化求解；具體而言，采用了 RMSD指數來度量某個像元與目標像元的相 似程度，通過與閾值比較來確定該像元是否為相似像元。
5. 如權利要求1所述的基于成對HJ-1B影像時空信息的地表溫度反演方法，其特征 在于，所述步驟S3中，統計滑動窗口內目標像元的相似像元數N，并依據N確定相應的求解 策略。
6. 如權利要求1所述的基于成對時相HJ-1B影像時空信息的地表溫度反演方法，其 特征在于，本發(fā)明中的兩個合理假設：（a)在正常天氣狀況下，短時間內在同一空間位置， 地表比輻射率不變；（b)小空間范圍內（本發(fā)明中用滑動窗口來表示）的大氣條件相同或差 異可以忽略；并且某一目標像元，在一定空間范圍內存在N個與其類別相同（或相似）的像 元（即相似像元）；在具體的實施中，本發(fā)明的特征還在于，所述步驟S4中，利用目標像元和 所選取的相似或準相似像元的輻射亮度值信息，構建一個簡化的輻射傳輸方程組，通過建 立一個帶約束的非線性優(yōu)化模型實現對該方程組的求解，得到目標像元在兩個成對時相的 地表輻射值。
7. 如權利要求1所述的基于成對時相HJ-1B影像時空信息的地表溫度反演方法，其 特征在于，所述步驟S5中，依據滑動窗口內相似像元(或準相似像元）具有相同的地表比輻 射率，且受到的大氣影響相同，則該空間范圍內各相似像元間的輻射亮度觀測量差異僅與 地表溫度的差異有關，基于這一原則可判定，對某目標像元的地表輻射的估計是否異常。
【文檔編號】G01J5/00GK104236716SQ201310238628
【公開日】2014年12月24日 申請日期:2013年6月17日 優(yōu)先權日:2013年6月17日
【發(fā)明者】陳 峰, 趙小鋒, 全元 申請人:中國科學院城市環(huán)境研究所