基于時(shí)空濾波方法的

多源遙感數據融合算法

神州信息

劉晶 李佳輝

1.

前言

高時(shí)空分辨率遙感數據在監測土地變化、作物生長(cháng)及物候參數反演等方面具有重要的作用，遙感數據的時(shí)空精細度成為了制約其在各領(lǐng)域應用的重要因素。 受技術(shù)限制的影響，傳感器在獲取地面數據時(shí)通常難以同時(shí)兼顧到數據的空間分辨率和時(shí)間分辨率。對于高時(shí)間分辨率的數據可以得到豐富的地表細節信息，但由于其重訪(fǎng)周期長(cháng)，再加上地面氣象因素( 云、雨和雪等) 對傳感器成像的影響，造成了單個(gè)傳感器有效數據重訪(fǎng)周期的延長(cháng)，使得監測具有很大程度上的“時(shí)空數據缺失”，導致監測研究中基礎觀(guān)測數據出現“空窗”，在研究的關(guān)鍵期不能滿(mǎn)足連續動(dòng)態(tài)跟蹤監測的要求。而對于低分辨遙感影像數，雖然具有較高的時(shí)間分辨率，但是其空間分辨率太低，對于我們監測較小面積種植作物區域具有很大影響。

為滿(mǎn)足地表信息遙感動(dòng)態(tài)監測同時(shí)對具有高空間分辨率和高時(shí)間分辨率特征遙感數據的需求，一些學(xué)者提出了一種能夠綜合高空間分辨率遙感數據的空間分辨率特征和高時(shí)間分辨率遙感數據的時(shí)間分辨率特征的技術(shù)，即遙感數據時(shí)空融合技術(shù)。

2.

多源時(shí)空數據融合技術(shù)

針對遙感數據時(shí)空融合方法，國內外學(xué)者進(jìn)行了大量研究。根據算法原理的不同，遙感數據時(shí)空融合方法可以分為2類(lèi)，即基于變換的模型和基于像元重構的模型。國內外大部分遙感數據時(shí)空融合模型是基于像元分解重構技術(shù)的，其基本思想是通過(guò)一定的規則選擇目標像元周邊的像元參與目標像元的重構。常用的時(shí)空融合方法有STARFM（Spatial and Temporal Adaptive Reflectance Fusion Model）和ESTARFM（Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model）兩種方法。

STARFM方法主要用于融合Landsat和MODIS影像, 生成與MODIS影像相應時(shí)間的模擬Landsat影像, 能較為準確地預測出高時(shí)空分辨率影像的反射率數據。但STARFM算法的預測結果受Landsat數據影響, 若Landsat影像中未記錄地表覆蓋在較短時(shí)間內的劇烈變化, 那么預測的影像也不會(huì )表征該情況，不能監測地表突變的情況。ESTARFM算法, 主要使用兩期對應觀(guān)測時(shí)間的高空間分辨率、低時(shí)間分辨率影像和低空間分辨率、高時(shí)間分辨率影像, 以及在預測日期獲取的低空間分辨率影像獲取對應時(shí)間的模擬高空間分辨率影像，該算法考慮了像元的異質(zhì)性,引入純凈像元與混合像元在一段時(shí)間內反射率變化的轉換系數, 提高了模型在異質(zhì)性較高區域的反射率融合精度, 在小區域實(shí)驗中有較好的效果，大區域效果還有待驗證。這兩種算法對于變化情況較大且未記錄在基準影像的區域預測精度較低，且算法精度受基準影像的數據質(zhì)量影響嚴重，如需構建大范圍的無(wú)云數據，需要保證每一個(gè)區域都能找到無(wú)云覆蓋的參考數據，這是非常困難的，特別是在大量云霧覆蓋區域，一年內都極少有無(wú)云數據，具有一定的局限性，不能滿(mǎn)足我們對全國范圍內作物進(jìn)行動(dòng)態(tài)監測的需求。

對于作物監測，我們主要需要作物關(guān)鍵物候期的時(shí)序性的高分辨率遙感數據，那么為了滿(mǎn)足我們的監測需求，具備全國作物動(dòng)態(tài)監測的能力，我們自研了SSTARFM（Self-tuning spatiotemporal adaptive reflectance fusion model）的時(shí)空遙感數據融合算法，該算法基于權重濾波與多源數據歸一化方法，根據不同尺度影像數據的時(shí)間空間特征，構建全時(shí)序、全空間范圍內的無(wú)云長(cháng)勢產(chǎn)品，完成了多源異構數據的統一,在保證高精度的條件下輸出高時(shí)間高分辨率的監測產(chǎn)品。

3.

SSTARFM多源數據時(shí)空融合方法

目前的時(shí)空融合算法大都是針對2種數據源的遙感數據進(jìn)行融合，沒(méi)有充分利用現有多源數據的優(yōu)勢，SSTARFM方法則基于目前我們常用的Sentinel2、Landsat8/9和MODIS 3種不同傳感器的數據進(jìn)行融合，提高了我們的監測精度。

SSTARFM算法流程

該算法的特點(diǎn)是在填補空缺值的基礎上，最大程度的保留了準確的原始數據。為了實(shí)現這種效果，需要對整個(gè)數據集經(jīng)歷兩次時(shí)空濾波處理。時(shí)空濾波處理首先通過(guò)空間插值的方式填補單幅影像的空缺值，但這樣填補出的空缺值只考慮的空間上的連續性，和真實(shí)值差異較大。因此需要使用時(shí)間濾波，以插值出的空間值為參考計算出在時(shí)間上的合理值。這種先空間后時(shí)間的方法既可以去除數據中的異常值的點(diǎn)，也可以在時(shí)空尺度上為空缺數據計算出一個(gè)相對合理的值。與STARFM算法不同的是，我們在空間維度基礎上增加了時(shí)間尺度進(jìn)行補值，增加了數據合理性和準確度，而且這種方法在長(cháng)時(shí)間序列的數據重構上，隨著(zhù)可參考的時(shí)間樣本增加，構建出的時(shí)序數據也會(huì )更加符合實(shí)際情況，精度更高。

本方法在實(shí)際應用的過(guò)程中，對于MODIS低分辨率數據，首先進(jìn)行一次時(shí)空濾波處理，構建出一個(gè)相對完整的數據產(chǎn)品，其次采用卷積和圖像補值的方法，構建Landsat和MODIS數據融合結果，然后再進(jìn)行時(shí)空濾波處理，保證了融合數據在時(shí)間和空間上的合理性，最后將Landsat8-MODIS融合結果和Sentinel2數據構建回歸模型，保證兩數據源的協(xié)同統一。

SSTARFM多源時(shí)空數據融合流程，注：示例數據為淄博市北部農田區域2022年03月06日數據

其中，雙立方（三次）卷積插值是我們進(jìn)行低分辨率數據降尺度的一個(gè)主要方式。雙三次插值又叫雙立方插值，用于在圖像中“插值”（Interpolating）或增加“像素”（Pixel）數量/密度的一種方法。其主要原理為：雙立方插值計算涉及到16個(gè)像素點(diǎn)，其中(i’, j’)表示待計算像素點(diǎn)在源圖像中的包含小數部分的像素坐標，dx表示X方向的小數坐標，dy表示Y方向的小數坐標，通過(guò)圖像16個(gè)像素點(diǎn)權重卷積之和即可計算得到新的像素值。具體如下圖所示：

雙立方（三次）卷積插值方法

基于SSTARFM算法，我們已經(jīng)具備生產(chǎn)8天一期10m分辨率的NDVI數據的能力，實(shí)現對作物區域進(jìn)行8天1次的生長(cháng)監測，提高了農情監測的時(shí)間顆粒度，可以做到對異常情況的及時(shí)反映，保障了農民的生產(chǎn)和財產(chǎn)安全，為政府部門(mén)的快速決策提供數據支撐。同時(shí)，SSTARFM時(shí)空融合算法解決了監測研究中基礎觀(guān)測數據出現“空窗”的問(wèn)題，能夠為農業(yè)上作物面積監測、長(cháng)勢、產(chǎn)量、苗情、土壤墑情、作物收割、成熟度、氣象災害、蟲(chóng)害等監測模型的構建提供強有力地數據支撐，為農田的精準監測提供保障。

神州信息遙感農情監測模型

4.

結語(yǔ)

在當前單一傳感器數據不能滿(mǎn)足各領(lǐng)域動(dòng)態(tài)監測對于高時(shí)空分辨率遙感數據需求的情況下，如何融合遙感數據綜合利用多源數據的特性和優(yōu)勢，具有重要的理論意義和應用價(jià)值。SSTARFM多源數據融合方法是一種具有良好性能的遙感時(shí)空融合算法，可以解決當前光學(xué)遙感數據空間分辨率和時(shí)間分辨率上相互制約的問(wèn)題，其發(fā)展及應用前景十分廣闊。