使用RSD进行高、中分数据交互计算（Landsat8-MODIS为例）.pdfVIP

使用RSD 进行高、中分数据交互计算（Landsat8-MODIS 为例） 想要同时交互处理的不同来源不同分辨率遥感数据大都有区域特性，比方我关心的区域是华 北平原、朝鲜半岛等等。RSD 进行这类处理前最好根据区域和空间分辨率要求先建立框架， 然后将五花八门的遥感数据加载，就可以处理了。RSD 将这类处理叫RO （Region Oriented ）， 区域朝向或者面向区域的。 但是通常大多数同学习惯的是数据是哪儿就处理哪儿。频繁地将数据打开关闭，加载了再存 出去。从RSD 的角度看这叫DO （Data Oriented ），数据朝向或者面向数据的。 为了看起来不那么突兀，适应大家对数据集打开关闭...的习惯，这里咱们先可着数据来，按 数据选择区域进行处理。 1. 打开一个Landsat8 数据集，使用反射通道30m 数据。 图1 打开的Landsat8 OLI 数据 这里是辽宁中西部的一个区域，打开投影为UTM 51 带。 2. 加载一个MODIS 的L1B 数据集 图2 加载MOD02 1km 数据的菜单 注意这里一定是“添加”而不是“打开”，如果选择了打开则RSD 又重新为你创建了一个新 任务，就不会将两种数据放在一起了。 选择一个MODIS L1B 1km 数据加载（这里用的不是一个时间的，暂时代用），在MODIS 数据 的预览图中，图像插值方式按缺省选择“稀疏像元图像”（也可以选择其它图像合成方式）。 注意这里的选择只是图像成像方式，不影响数据。选稀疏图像速度要快一些。 图3 加载到Landsat8 数据上的MODIS 的稀疏图像 注意到Landsat8 数据周围麻麻沙沙小点子了？这就是MODIS 的1km 数据，由于数据分辨率 不同，这里又没有插值，数据都是稀疏分布的。其实整个区域都存在这样的点子，但是在有 Landsat8 数据的区域看不见，只是在四周才能看见这些点子。 这时在层列表的窗口多了一层数据 （层2 ）。并且这个层的类型是 SWATH 。 图4 新加载的MODIS SWATH 层 SWATH 层是 RSD 为多尺度数据协同处理而特意保留的数据组织格式。它保留栅格数据和矢 量数据两种数据特征，就是说可以当栅格数据处理，也可以当矢量数据处理。例如：我们需 要同时处理1m 空间尺度的高分数据和1km 的MODIS 等中分数据。空间范围是整个东亚地 区，这时我们要是将MODIS 数据按栅格数据重采样为1m，那数据量大的岂不是要爆炸？可 是如果仍按SWATH 像元位置描述这个MODIS 数据集，就不会发生数据膨胀（这些内容容后 介绍）。 在这个例子里，数据尺度差异不是很大 （30m 到1km），数据范围比较小，我们可以将其重 采样为GRID 数据，与Landsat8 数据一起处理。 3. MODIS 数据重采样 可以根据需要选择重采样方法。 图5 对选中的SWATH 层进行最近邻点法重采样 重采样结果如下图。 图6 对MODIS SWATH 层进行最近邻点法重采样的结果 最近邻点法重采样的结果将 MODIS 的每个像元表示为一个个的小格子，虽然看起来粗糙， 但是真实反映出了像元位置的实际通道值，和其它重采样方式比较各有特点。 注意现在我们多了一个层3，MODIS 数据的GRID 层。Landsat8 的数据呢？ 是层1，被挡在 后面了。 当然可以分开观察这两种数据集的图像，也可以融合在一起看一看，大损失这里还是直接设 置半透明吧简单观察一下吧。 图7 Landsat8 和MODIS 半透明显示 也可以点击外面层的小眼睛，关闭外面的层使之不可见，里面层的图像就露出来了。 RSD 的图像和数据集是严格分开的，这儿显示不管怎么折腾后面的数据集是不变的，所以这 个图就是给你直观看看而已。 4. Landsat8 数据和MODIS 数据放在一起算一算 都是GRID 数据时就可以相互之间进行任意计算了，或者任意交换两个数据集的任意通道数 据都是可以的。这里使用 RSD 的数据集运算工具 （DSM）进行一下示范，将 Landsat8 的红 光波段数据与 MODIS 的红光波段数据计算一下平均值。当然也可以组合其它任何需要的运 算。 图8 启动数据集运算工具的菜单命令 点击后出现下面的对话框。 图9 DSM 对话框 在这个对话框中的公式输入窗口（f(x)= ）输入计算公式 : (_L1_B4+_L3_B1)/2 在这个公式中，一对括号表示优先计算，“_L1 ”表示层1，“_B4 ”表示4 通道，这样，括号 内部分表示“层1 的4 通道加上层3 的1 通道