想要同時(shí)交互處理的不同來(lái)源不同分辨率遙感數(shù)據(jù)大都有區(qū)域特性，比方我關(guān)心的區(qū)域是華北平原、朝鮮半島等等。RSD進(jìn)行這類處理前最好根據(jù)區(qū)域和空間分辨率要求先建立框架，然后將五花八門(mén)的遙感數(shù)據(jù)加載，就可以處理了。RSD將這類處理叫RO（Region Oriented），區(qū)域朝向或者面向區(qū)域的。

但是通常大多數(shù)同學(xué)習(xí)慣的是數(shù)據(jù)是哪兒就處理哪兒。頻繁地將數(shù)據(jù)打開(kāi)關(guān)閉，加載了再存出去。從RSD的角度看這叫DO（Data Oriented），數(shù)據(jù)朝向或者面向數(shù)據(jù)的。

為了看起來(lái)不那么突兀，適應(yīng)大家對(duì)數(shù)據(jù)集打開(kāi)關(guān)閉...的習(xí)慣，這里咱們先可著數(shù)據(jù)來(lái)，按數(shù)據(jù)選擇區(qū)域進(jìn)行處理。

1. 打開(kāi)一個(gè)Landsat8數(shù)據(jù)集，使用反射通道30m數(shù)據(jù)。

圖1 打開(kāi)的Landsat8 OLI數(shù)據(jù)

這里是遼寧中西部的一個(gè)區(qū)域，打開(kāi)投影為UTM 51帶。

2. 加載一個(gè)MODIS的L1B數(shù)據(jù)集

圖2 加載MOD02 1km數(shù)據(jù)的菜單

注意這里一定是“添加”而不是“打開(kāi)”，如果選擇了打開(kāi)則RSD又重新為你創(chuàng)建了一個(gè)新任務(wù)，就不會(huì)將兩種數(shù)據(jù)放在一起了。

選擇一個(gè)MODIS L1B1km數(shù)據(jù)加載（這里用的不是一個(gè)時(shí)間的，暫時(shí)代用），在MODIS數(shù)據(jù)的預(yù)覽圖中，圖像插值方式按缺省選擇“稀疏像元圖像”（也可以選擇其它圖像合成方式）。注意這里的選擇只是圖像成像方式，不影響數(shù)據(jù)。選稀疏圖像速度要快一些。

圖3 加載到Landsat8數(shù)據(jù)上的MODIS的稀疏圖像

注意到Landsat8數(shù)據(jù)周?chē)槁樯成承↑c(diǎn)子了？這就是MODIS的1km數(shù)據(jù)，由于數(shù)據(jù)分辨率不同，這里又沒(méi)有插值，數(shù)據(jù)都是稀疏分布的。其實(shí)整個(gè)區(qū)域都存在這樣的點(diǎn)子，但是在有Landsat8數(shù)據(jù)的區(qū)域看不見(jiàn)，只是在四周才能看見(jiàn)這些點(diǎn)子。

這時(shí)在層列表的窗口多了一層數(shù)據(jù)（層2）。并且這個(gè)層的類型是 SWATH。

圖4 新加載的MODIS SWATH 層

SWATH層是RSD為多尺度數(shù)據(jù)協(xié)同處理而特意保留的數(shù)據(jù)組織格式。它保留柵格數(shù)據(jù)和矢量數(shù)據(jù)兩種數(shù)據(jù)特征，就是說(shuō)可以當(dāng)柵格數(shù)據(jù)處理，也可以當(dāng)矢量數(shù)據(jù)處理。例如：我們需要同時(shí)處理1m空間尺度的高分?jǐn)?shù)據(jù)和1km的MODIS等中分?jǐn)?shù)據(jù)。空間范圍是整個(gè)東亞地區(qū)，這時(shí)我們要是將MODIS數(shù)據(jù)按柵格數(shù)據(jù)重采樣為1m，那數(shù)據(jù)量大的豈不是要爆炸？可是如果仍按SWATH像元位置描述這個(gè)MODIS數(shù)據(jù)集，就不會(huì)發(fā)生數(shù)據(jù)膨脹（這些內(nèi)容容后介紹）。

在這個(gè)例子里，數(shù)據(jù)尺度差異不是很大（30m到1km），數(shù)據(jù)范圍比較小，我們可以將其重采樣為GRID數(shù)據(jù)，與Landsat8數(shù)據(jù)一起處理。

3.MODIS數(shù)據(jù)重采樣

可以根據(jù)需要選擇重采樣方法。

圖5 對(duì)選中的SWATH層進(jìn)行最近鄰點(diǎn)法重采樣

重采樣結(jié)果如下圖。

圖6 對(duì)MODIS SWATH層進(jìn)行最近鄰點(diǎn)法重采樣的結(jié)果

最近鄰點(diǎn)法重采樣的結(jié)果將MODIS的每個(gè)像元表示為一個(gè)個(gè)的小格子，雖然看起來(lái)粗糙，但是真實(shí)反映出了像元位置的實(shí)際通道值，和其它重采樣方式比較各有特點(diǎn)。

注意現(xiàn)在我們多了一個(gè)層3，MODIS 數(shù)據(jù)的GRID層。Landsat8的數(shù)據(jù)呢？ 是層1，被擋在后面了。

當(dāng)然可以分開(kāi)觀察這兩種數(shù)據(jù)集的圖像，也可以融合在一起看一看，大損失這里還是直接設(shè)置半透明吧簡(jiǎn)單觀察一下吧。

圖7 Landsat8和MODIS 半透明顯示

也可以點(diǎn)擊外面層的小眼睛，關(guān)閉外面的層使之不可見(jiàn)，里面層的圖像就露出來(lái)了。

RSD的圖像和數(shù)據(jù)集是嚴(yán)格分開(kāi)的，這兒顯示不管怎么折騰后面的數(shù)據(jù)集是不變的，所以這個(gè)圖就是給你直觀看看而已。

4.? Landsat8數(shù)據(jù)和MODIS數(shù)據(jù)放在一起算一算

都是GRID數(shù)據(jù)時(shí)就可以相互之間進(jìn)行任意計(jì)算了，或者任意交換兩個(gè)數(shù)據(jù)集的任意通道數(shù)據(jù)都是可以的。這里使用RSD的數(shù)據(jù)集運(yùn)算工具（DSM）進(jìn)行一下示范，將Landsat8的紅光波段數(shù)據(jù)與MODIS的紅光波段數(shù)據(jù)計(jì)算一下平均值。當(dāng)然也可以組合其它任何需要的運(yùn)算。

圖8 啟動(dòng)數(shù)據(jù)集運(yùn)算工具的菜單命令

點(diǎn)擊后出現(xiàn)下面的對(duì)話框。

圖9 DSM對(duì)話框

在這個(gè)對(duì)話框中的公式輸入窗口（f(x)=）輸入計(jì)算公式 :

(_L1_B4+_L3_B1)/2

在這個(gè)公式中，一對(duì)括號(hào)表示優(yōu)先計(jì)算，“_L1”表示層1，“_B4”表示4通道，這樣，括號(hào)內(nèi)部分表示“層1的4通道加上層3的1通道”。然后括號(hào)外面除2，計(jì)算二者平均值。

注意這種計(jì)算是有效數(shù)據(jù)計(jì)算，遇到無(wú)效數(shù)據(jù)按空值處理。有關(guān)這個(gè)計(jì)算工具的更多信息請(qǐng)參閱相關(guān)的說(shuō)明。

點(diǎn)擊“開(kāi)始計(jì)算”。計(jì)算完畢后，RSD任務(wù)多了一個(gè)新層（層4），這就是計(jì)算結(jié)果。

圖10 DSM計(jì)算Landsat8和MODIS紅光波段平均值

值得注意的是：很多同學(xué)由于使用過(guò)其它軟件的原因，習(xí)慣性地上磁盤(pán)文件上找計(jì)算結(jié)果，還有人問(wèn)我怎么沒(méi)有指定輸出文件。這里統(tǒng)一回答下，RSD所有計(jì)算中間結(jié)果都在任務(wù)中保留，不會(huì)向外存輸出一大堆的中間文件，只有任務(wù)結(jié)束時(shí)才整體保存任務(wù)，或者導(dǎo)出特定的輸出結(jié)果。

RSD任務(wù)中的中間數(shù)據(jù)可以隨時(shí)刪除，例如此例中的層2（MODIS的SWATH），如果你確定不用了，就可以右擊它，在菜單中選擇刪除。

DSM計(jì)算結(jié)果保存為4字節(jié)浮點(diǎn)型數(shù)據(jù)，也可以右擊隨時(shí)更改為需要的數(shù)據(jù)類型。

5. 討論一下

① 不用關(guān)心投影問(wèn)題，框架是什么投影，后面的數(shù)據(jù)自動(dòng)使用該投影。事實(shí)上RSD基本上不需要用戶進(jìn)行投影轉(zhuǎn)換，RSD的投影引擎在后臺(tái)隨時(shí)為你進(jìn)行投影變換。此例中由于開(kāi)始打開(kāi)Landsat8數(shù)據(jù)時(shí)創(chuàng)建的框架是UTM，后面就全部自動(dòng)是UTM了。

② 自動(dòng)的尺度轉(zhuǎn)換和空間位置對(duì)齊。雖然MODIS數(shù)據(jù)范圍和像元尺度與Landsat8完全不同，但是在一致的框架下自動(dòng)對(duì)齊。

③ RSD當(dāng)然不是只將Landsat8和MODIS一起計(jì)算。舉一反三，只要RSD能加載的數(shù)據(jù)交互計(jì)算協(xié)同處理都是暢通無(wú)阻的。

④ 就目前相對(duì)成熟的遙感數(shù)據(jù)處理平臺(tái)來(lái)看，對(duì)異源的、多尺度的遙感數(shù)據(jù)協(xié)同計(jì)算處理，RSD提供了前所未有的便利。該功能對(duì)日益增多的遙感數(shù)據(jù)源是非常重要的，同時(shí)也是RSD獨(dú)有的。

⑤ RSD同時(shí)提供基于SWATH的大范圍中分?jǐn)?shù)據(jù)的高分辨率框架稀疏處理機(jī)制，請(qǐng)關(guān)注。

?

?

在百度云盤(pán)https://pan.baidu.com/s/1T-LBvaD_zVCwJsGf_hCyCg下載新版的SetupRsd2.1.4.rar。如果你是第一次安裝，參考下軟件包里面的安裝說(shuō)明，最好安裝相關(guān)輔助數(shù)據(jù)。

詳情加企鵝群136965427，在這里解答和討論有關(guān)遙感數(shù)據(jù)處理和RSD平臺(tái)的有關(guān)技術(shù)問(wèn)題。

?

?

總結(jié)

以上是生活随笔為你收集整理的使用RSD进行高、中分数据交互计算（以Landsat8-MODIS为例）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。