遥感数据处理流程新概念(一) ——RSD的面向任务

目录

一 创建任务... 2

二 为任务添加数据... 3

三 拼接... 5

四 剪切... 11

五 不落地继续分类?... 13

六 创建/保存任务目录... 13

七 讨论... 14

 

提到了流程人们就会想起一步步的处理步骤,或者进行程序化处理。这里我说程序化处理是指写代码控制的处理过程。通常由于需求的多样性没有办法写一个放之四海而皆准的程序,同时你又不可能要求每一个用户都是高手程序员,我们就先看看在平台上的交互处理的流程吧。

RSD是面向任务的,但是就目前来看,虽然用户数增长很快,但大家还仍然把它仅仅当作一个前期处理的工具来用的。处理得到一些初级产品就导出数据然后离开了。一个主要原因是大家对RSD很多设计功能还不了解(主要怪我还没写说明书)。还有就是新的处理概念马上得到大家的认同也不太容易,毕竟有传统的处理习惯的惯性在那儿。我们不着急慢慢来,逐步介绍。这里先举一个例子,如何不落地将6景GF-1 MSS 8m多光谱数据处理为切分好的区域数据。

 

盘锦那个地方产好吃的大米,水稻长势的监测对农业决策部门很重要。现在我们来处理一下来自该地区插秧前的6个GF1-MSS 的8米多光谱数据(图1)。

图1

这6个数据集是从资源卫星中心下载后未经进一步处理的L1A。

创建任务

一开始就需要创建任务框架是RSD推广很头疼的问题。很多人这样抱怨创建任务:还没等开始呢你就给我出了个难题?其实难是不难,但是没有人愿意还没开始前,上来就输入一些莫名其妙的参数。这也是事实。

后来我想了一个办法,把全国各个省、市、县的框架模板都做好,到时候一打开就可以了。即使没有正好合适的模板,也可以在其它模板的基础上裁裁剪剪得到一个。

在这里下载模板文件https://pan.baidu.com/s/1nvIJekT ,找到里面“上传示例数据\RsdTemplate”目录,然后找到你需要的省市县,把对应的 .tpl文件下载来就可以了。比方我需要辽宁省盘锦市的模板文件,找到“盘锦市-UTM-10m.tpl”

下载到本地,然后在RSD中打开,见图2。

图2

预览图有一个盘锦市的轮廓(也可能没有)。打开后就是一个空白,可以添加一个行政区边界的矢量数据看起来有点目的性,也可以就这样放着,没关系。

这里有一个问题,我们的框架是10m的,但是数据是8m,这样直接用等于把分辨率降低了。我们可以按下述步骤复制保存一个8m的框架,再打开这个8m的就可以了。

选择“框架管理->复制框架”,弹出来一个对话框,见图3。

图3

将图2最下面的像元尺寸改成8,这时框架的宽度和高度都会发生变化,不用管它。在最上面可以更换一个名称保存。点“确定”就可以保存为一个新的.tpl文件。关闭当前的任务,再重新打开这个新的.tpl文件就可以了。

 

为任务添加数据

如果希望以后加载的数据都在这个任务里,请使用“添加”的功能。通常RSD默认的操作是“打开”创建一个新任务,而“添加”是将数据添加到当前任务中。

这里我们在“新建->打开”为任务添加数据(这里的打开是个例外,后面有添加功能)。(图4)

图4

定位到预先准备的6个数据集,点击其中的一个TIFF文件,就会在右侧窗口出现一个预览图,右击预览图,在菜单点击“文件添加到当前任务”。就可以看见预览图的添加结果(图5)。

图5

在图5的左下角可以看见添加数据的预览。然后依次添加全部6个数据集(图6)。

图6

图6左上角红框内每行对应一个打开的数据集的层。现在关闭打开对话框,就可以进行后续的操作了。

本节演示的是加载GF-1 MSS数据,加载其它数据或者需要大气校正融合等,请参阅后续介绍。

拼接

其实这时位置上已经对齐了,不拼接也可以单独一个个的处理,或者拼接后再分幅也可以。这里我们要拼接成一个大图。

这里拼接实际上就是把两幅往一块一合并就可以了,要注意的是每景的边缘可能有异常数据的黑边,需要去掉。还有就是被合并的两景中保留谁的多一些的问题。

现在开始合并1、2层。先关掉3、4、5、6层的显示(点层列表的前头的小眼睛,注意TFS时这个不支持自动刷新)。

图7

在图7中,点击“现有通道数据合并到”,出现图8。

图8

注意:层列表中勾选的是层1,这里选择的是层2,是指1、2层合并,交集部分保留第二层的数据。合并后出来一个新层7,这个新层包含了1、2两层的数据,交集部分保留的是第2层的数据。图9是合并后的结果。

图9

这两层合并太完美了,没有遇到一点问题,有问题后面再讨论吧。

现在,再重复上述合并1、2层的步骤,把层3合并到7里面,生成一个新层8。仔细观察合并后层8的结果(图10),注意这时再第3和第7层交界出现了一条绿线,这条绿线是层7(原来的层2)的边缘。

图10

这样的异常边缘可以通过边缘腐蚀消除(或者后面将的切割)。

删除这个层8不要了,勾选层7,点击下面的菜单命令消除边缘的异常数据。点击一次如果还没有完全消除就再来一次,直至满意。

 

再次重新合并3和7(图11)。这次好了,边缘的异常数据消失了。

图11

由于原来合成的层8已经删除了,这次新合成的层还是层8。

下面继续合成层8和层4,见图12。可见有较大的色差,这是合成图像时拉伸不一致产生的,可以不用管它。

图12

合成后结果见图13,色差基本消失了。原因是新添加进来的数据使用了和原数据一样的合成图像拉伸参数,同时还是同一天的数据。不是人为处理的。

图13

可见色调基本上是一致的,但是仔细看右侧水体上还是有一点色差。这是数据本身的差异,不是图像合成的问题。

继续合成层5、6。

图14

最后的合成结果是层11,现在就是所有的6个数据集都合成在一个层(层11)里面了。

我们把中间结果的层7、8、9、10都删掉。层11上升到层7。注意:这些多余的层并不占用很多的内存和CPU资源,只是占据了比较多的磁盘空间,删不删可以视情况而定。如果不需要了,也可以把1~6全部删除。

RSD拼接合并的速度不用担心,都是几秒的时间。

注意:这里是最理想的拼接状况,实际情况可能比这要复杂的多。进一步的拼接问题见《遥感数据处理流程新概念 (二)》

剪切

继续使用上述拼接处理好的数据。先添加一个矢量数据的行政边界,这里需要面数据(Polygon)来转换成RSD蒙板。

图15

然后勾选这个矢量层,再点击

图16

这时候会发现创建了一个蒙版层。打开蒙版属性的浮动窗口,蒙版层的参数见图17。

图17

蒙板属性里面出现的地名等来自于导入的矢量数据的属性。蒙板的一些其它应用容后介绍。

在这里,注意在主窗口移动鼠标时,在主窗口的状态栏里面会动态显示鼠标所在的蒙版区域,比如:“大洼县”、“盘山县”…

勾选一个栅格数据层,比如上述合并后的层7,让光标移动到所在蒙版区域显示的为“大洼县”的位置,右击鼠标出现图18的弹出窗口。

图18

选择“复制当前蒙板区域”,可以看见创建了一个新层。如果刚才添加的矢量数据占用了层8的位置,这个新层应该是层9。关闭其它层,只显示层9 ,见图19。

图19

可见已经复制了我们感兴趣的大洼县的区域数据。

同样的方法我们再复制一个蒙版区域的数据“盘山县”。

图20

复制后的另一个区域,盘山县。余类推。

 

不落地继续分类?

可以,RSD有基于像元的监督、非监督分类,基于知识特征的面向对象分类和面向的支持向量机分类。参阅

https://blog.csdn.net/gordon3000/article/details/84059254

https://blog.csdn.net/gordon3000/article/details/85534116

https://blog.csdn.net/gordon3000/article/details/85586551

 

创建/保存任务目录

为数据安全和大规模数据存取方便,建议使用“任务目录”。

https://blog.csdn.net/gordon3000/article/details/90371911

任务目录的一个福利是下次加载使用时能实现秒开,很少受数据规模大小的影响。

 

讨论

这种不落地的处理操作目标明确,逻辑清晰。不需要用户频繁从磁盘查找和写出数据。其实,本质上RSD与磁盘也是频繁交换数据的(实际上还要更多),只是RSD把这部分繁杂给隐藏了,不需要用户亲自控制这些中间数据。

这里的例子是一个非常理想的情况,应用实践中还存在很多现实问题。比如融合后的大规模数据集如何处理;如何确定拼接交集的比例、是否需要拼接线;拼接中对不齐怎么办等。这些问题RSD都有解决方案,容后慢慢介绍。

下一篇《遥感数据处理流程新概念(二) ——再谈拼接》 

 https://blog.csdn.net/gordon3000/article/details/93974798

 

软件下载地址:https://pan.baidu.com/s/1T-LBvaD_zVCwJsGf_hCyCg

详情加企鹅群136965427,在这里解答和讨论有关遥感数据处理和RSD平台的有关技术问题。

 

 

发布了74 篇原创文章 · 获赞 27 · 访问量 5万+

猜你喜欢

转载自blog.csdn.net/gordon3000/article/details/93171572
今日推荐