遥感解译方法
为了准确进行遥感地质解译,解译者首先应具备一定的地质、遥感知识;其次应对解译区的地质基础、构造格架、灾害地质、地形地貌和水文情况等要有粗略的了解。
解译的原则应采用由已知到未知、从区域到局部、先易后难、由宏观到微观,从总体到个别,从定性到定量,循序渐进,不断反馈和逐步深化的方法进行工作
操作方法
- 01
1、首先在arggis中,打开影像图和clipshape文件,如图所示:
- 02
再用editor编辑工具
- 03
在clipshape图层上画出分割线,分割出差不多的图版块。(如果分成三块,还要在属性中将其添加“子文件名”属性,命名方式如hy4_1,hy4-2,hy4_3)
- 04
再在endas中将这一块进行裁剪。这个过程的具体步骤已在“遥感影像裁剪流程”一文中写过。将其结果以我个人的所做的为例命名为:hy4_1.img(保存在SubImage文件夹之中)
- 05
用endas进行监督分类1:{在做这个步骤的时候之前要在arggis中将pr外调成果文件夹中的“玉米样放”图层打开,(玉米代号1105)。以裁剪下来的影像或者原影像做底图,观看玉米的颜色是什么,然后再进行监督分类1},如下步骤。 先在viewer窗口中将裁剪后的成果影像图A(hy4_1.img)打开,再点击
- 06
选择
- 07
接着会弹出这个窗口。
- 08
在viewer中点击
- 09
会弹出工具栏如下
- 10
用工具栏中的
- 11
这个工具在viewer选中一定数量的农作物斑块面(与上面的{……}中玉米样色尽量一致)。 在signature ditor窗口中用点击这个进行类添加,一般在数量上依情况而定。 目前做的是玉米解译项目,所以只分为两类,一类为玉米,另一类为非玉米。每一类至少要取六个斑块。然后将所选的同一种斑块进行合并,用的是这个工具进行合并。 如下图只是玉米,非玉米类采集的图版类:
- 12
我们只分两类,所以在合并数据表图中,将signature name改成class1(玉米)和class2(非玉米),并且将value和order数值分别改成1、2如下所示:
- 13
这样就行了,在进行保存(file——save)会弹出这个窗口进行保存,保存到老师设置的文件夹中的sig文件夹中,文件命名为hy4_1.sig。
- 14
监督分类2:接着从
- 15
在input raster file中打开裁剪的影像图(hy4_1.img),在input signature file中的打开监督分类1,保存的sig文件(hy4_1.sig)。在classified中进行保存,保存在class文件夹中,命名为hy4_1_cls.img。在
- 16
中选择第二个。然后点击ok了。 【在做下一步之前要进行评估,所谓的评估就是看你采集的玉米图斑是不是与影像地图的斑块大概一致,如果出现大片或者明显的玉米斑块被非玉米斑块所代替,那监督分类1最好重做一遍,以提高精度。在viewer窗口中打开监督分类1的成果图和裁剪的影像图(最好是原影像图,比较大看的真切),进行滚动观看。
- 17
去小碎点
- 18
会弹出如下窗口。
- 19
在Input file中打开监督分类的结果文件,即class中的文件。在output file中将文件保存在clump文件夹中,命名为hy4_1_clump.img。勾选再确定就ok了,接着在
- 20
会弹出
- 21
在input file中打开上一个步骤的结果文件(hy4_1_clump.img), 在outfile中保存在eliminate文件夹中,命名为hy4_1_eli.img。
- 22
栅格文件转化矢量文件 点击
- 23
会弹出如下窗口。
- 24
从input raster那打开“eliminate文件夹”中的成果文件(hy4_1_eli.img),在output vecter保存文件到vecter文件夹下的cov文件夹中,命名为hy4s1。点击ok就可以了。 到了这里,遥感玉米解译项目在endas软件中的操作就结束了。接下来是在arggis中进行最后解译的步骤。(至于各种文件命名的方式,会因个人情况而不同,下面用到的不一定是)
- 25
在arggis中打开vecter文件夹中的矢量图hy4s1,再打开clipshape文件。如下图所示。
- 26
下一部将矢量文件中的玉米数据导出来,因为在做监督分类的时候就已经将玉米和非玉米数值分别改为1和2,所以只要将含有1的属性部分提取出来就行了。下面以hy4s1为例。 右击再点击open attribute table会弹出如下窗口。
- 27
现在你所看到的值全部是1,那是因为我用的材料已被我做过了。在你刚完成矢量转化的之后,你的文件atttributes of hy4s1窗口中的grid_code属性分别为1和2.有时会出现0。 所以我们要将含有1的部分提取出来,因为我们只做玉米这一部分的。首先要进行选定他。步骤如下。 在窗口的最下面右击一下select by attributes会弹出如下窗口。
- 28
在这个窗口中双击“GRID_CODE”在下面就会出现它,然后点击=,在点击Get Unique values,会出现1、2或者有0的情况,接着双击1,再点击apply就行了。你会发现在atttributes of hy4s1窗口中所有的含有1的都被选中了。 然后右击hy4s1——Date——Export Date,会弹出如下窗口。
- 29
为了区分开要将其保存在特定的文件夹中如vecter文件夹中,并且将其特殊命名hy4_1_cls。再直接点击ok就可以了,之后会弹出如下窗口。
- 30
点击是,就会在左侧的窗口中将你导出的数据添加到你的图层中了。如下:
- 31
然后将hy4_1_cls复制到part文件夹下,并且命名为hy4_1_ipt,再将其在arggis中打开,右击open attribute table,在里面添加CODE,C字段属性,并且在CODE下方的数值全部改成1501(玉米),另一个不负值(为0)
- 32
这一步就是真正的解译步骤,在上一步的得出的数据图中(同时打开影像图进行对比参照,确定哪些是玉米,与得出的数据图斑是否一致),将不是玉米的斑块进行挖洞删除;将是玉米却被列为非玉米的行列的部分要进行融合。如果在一个图斑中出现道路,需要将其画出。 挖洞:假如下面的图斑中有一块不是玉米的要将其挖洞删除如下。
- 33
先Start Editing,点击
- 34
在图斑上将那一块非玉米的画出来。
- 35
再点击Editor中的clip, 会弹出如下窗口。
- 36
(也可以设置快捷键,如按住shift S)直接点击ok,再点击一下Delete键进行删除就可以了。如下图所示:
- 37
融合: 这一步将是玉米部分却显示非玉米的黑洞现象的小图斑,将其融合。下面的小洞就是要处理的对象,将其融合在大图斑中。就以下图为例:
- 38
首先要将其中的黑洞圈选,用的是
- 39
选中之后里面的黑洞会看不见的(有的能看见),但是依然存在。