遥感作业—中国地质大学案例.doc

遥感作业 021111班 按作业要求： 第一步、裁剪 23 43 102 37 我的左下角的坐标如截图，15*×15* 用的波段是TM5、4、3 大概制作过程：裁剪在ENVI的5.0版本中可以直接裁剪，然后输出就行。在4.8版本中需要按上课实习的步骤，但要注意，得将所给的坐标，即左下角的坐标换算成左上角和右下角的坐标（加减15*就行了），然后分别输入进去。最后得到如上图的结果，其所使用的波段是TM5、4、3，直接用来解译。（其中图幅的坐标我没有标上去） 第二步、图像增强 对裁剪后的进行线性拉伸，如上图所示，其增强的效果并不是很明显。 第三步、RGB假彩色合成 取TM4、3、2波段赋予红绿蓝色，即标准的假彩色合成方法，以突出水体、城区、山区及线性特征。如图： 另外，选择B4、B5、B3波段依次被赋予红、绿、蓝进行RGB图像合成进行比较,如图： 感觉两者的效果差别不是很大，相对而言还是标准假彩色合成的效果更明显。 这是真彩色的图像： 第四步、铁染异常: 主成分分析提取 用B1、B3、B4、B5四个波段的图像做主成分分析。由于贴的氧化物特征光谱信息集中在TM1--4波段，在TM3波段成高反射，在TM4和TM1波段相对吸收，同时为了避免含羟基矿物的干扰，注意在选取波段时要舍弃7波段。 以铁染信息为例，利用“均值+N倍标准差”作为阈值对异常图像进行密度分割，对研究区铁染异常信息进行提取。本次实验将异常分为三级，标准差系数分别为2、2.5、3倍。 通过上一步分析特征向量后，得出铁染蚀变异常信息提取中1345文件的第4主成分主要包含铁染信息，加载第4主成分的灰度影像。 在主菜单栏中点击Basic Tools，在下拉菜单中单击Statistics，再选择Compute Statistics。选择文件“铁染”，进行均值、方差、最大值、最小值进行统计。 从结果中可以得到第4主成分的均值（Mean）为0，方差（Stdev）为92.3。 3.等级划分。本次实验将异常分为三级，标准差系数分别为2、2.5、3倍。一级异常的下限值为3×92.3=276.9，二级异常2.5×92.3=230.75、三级异常2×92.3=184.6。利用密度分割对异常信息进行提取。 4.在主图像窗口-Overlay-Density Slice，通过修改编辑数据范围和颜色，将铁染异常结果分为3级。 5.选择Density Slice 对话框中的File-Output Range to Class Image，可以将分割结果输出为分类格式显示如下： 6.对分类结果进行滤波处理。选择主菜单中Filter-Convolution and Morphology,打开Convolution and Morphology Tool对话框，点击Convolution菜单选择Median中值滤波，选择内核算子kernel size为5。点击Apply To File选择上一步得到的分类结果图，对铁染蚀变信息进行5*5中值滤波。 7.结果叠加。以灰度显示方式在主显示窗口Image中显示原始数据。点击Overlay-Classfication，选择前面得到的分类结果，可以叠加在原始图像上，中值滤波后结果及与原图叠加后显示结果如下： 最后输出异常图如下： 分析： 主成分分析提取羟基蚀变信息 方法类似，同理可得异常图： 第五节、线环解译 打开第四歩裁减的RGB假彩色图层，对其再进行空间滤波增强处理，主线性拉伸等来突出河流等一些呈线性分布的地形，因而有利于线性和环形解译。得到许多处理的图像，对它们进行对比，进行线环解译。选择的图层是B5、B4、B3三波段RGB假彩色合成图作为解译的底图，这样图像较为清晰。 1.线性构造：该地区线性构造发育明显，在图的西部出现走向为南北向的负地形，如红色所标的位置，预计是断裂，褶皱等等一些构造的线性表现，具体的还的结合其他的资料。 2.环形构造：环形构造主要发育在图幅的东北角，范围不是很大，结合空间滤波分析，在该处发现有一定的Fe染蚀变，与环形构造未知大致相符，解译该处环形构造为侵入的隐伏岩体，当然了，还有可能是一些简单的蚀变带，具体的信息还要结合其他的地质信息，由于我基础有限，分析不出来，老师见谅。 总结： 通过这次系统的遥感实习解译，我们对ENVI这个常用的遥感解译软件的解译方法和步骤有了深刻的了解和学习。结合《遥感导论》的理论基础，以地质知识为前提，在ENVI软件的支持下，我们通过自己的理解，解译出了遥感图像。由于学习还不是很深入，可能解译中存在许