中国地质大学遥感图像处理上机实习报告.docx

1、中国地质大学遥感图像处理上机实习报告遥感图像处理课程实习报告学生姓名：王蜀越 班 学 号： 学 号： 指导教师： 王红平、许凯 中国地质大学信息工程学院2017年 7月 1日实习一：影像融合1.1 【实习目的】1.掌握遥感影像融合的基本原理、内容和要点；2.掌握在ENVI中遥感影像融合的操作方法、步骤；3.进一步熟悉ENVI软件的操作。1.2 【实习步骤】（1）实习流程图：（2）参数设置及中间结果：无1.3 【实习过程】（1）影像融合的内容图像融合（Image Fusion）是指将多源信道所采集到的关于同一目标的图像数据经过图像处理和计算机技术等，最大限度的提取各自信道中的有利信息，最后综合成

2、高质量的图像，以提高图像信息的利用率、改善计算机解译精度和可靠性、提升原始图像的空间分辨率和光谱分辨率，利于监测。高效的图像融合方法可以根据需要综合处理多源通道的信息，从而有效地提高了图像信息的利用率、系统对目标探测识别地可靠性及系统的自动化程度。其目的是将单一传感器的多波段信息或不同类传感器所提供的信息加以综合，消除多传感器信息之间可能存在的冗余和矛盾，以增强影像中信息透明度，改善解译的精度、可靠性以及使用率，以形成对目标的清晰、完整、准确的信息描述。（2）图像融合过程 （3）几种常见的融合算法IHS融合算法在计算机内定量处理色彩时，通常采用RGB表色系统，但在视觉上定性地描述色彩时，采用I

3、HS显色系统更直观些。IHS显色系统采用色调（Hue）、饱和度（Saturation）、亮度（Intensity）表示颜色。IHS融合又称为“I分量替换法”，即用全色图像代替多光谱图像的I分量，并与其H和S分量相结合进行IHS逆变换获得融合图像。IHS图像融合的基本步骤：（1）对两幅图像进行几何配准，并将多光谱图像重采样与全色图像分辨率设置一致；（2）对多光谱图像进行IHS变换，将图像变换之IHS空间；（3）对全色图像和多光谱图像的I分量进行直方图匹配；（4）用直方图配准后的全色图像代替多光谱图像的I分量；（5）对I分量替换后的多光谱图像进行IHS逆变换至RGB空间，得到融合图像。通过变换、替

4、代、逆变换获得的融合图像既具有全色高分辨率的特点，又保持了多光谱图像的色调和饱和度。虽然融合后图像清晰度提高了，但光谱信息损失严重，即产生颜色失真。如果融合结果图像应用以光谱分析为主就不太适合选择该方法。PCA变换法 PCA（Principal Component Analysis）变换又称为K-L变换（主成分变换）。它是将N个波段的低分辨率图像进行KL变换，将单波段的高分辨率图像经过灰度拉伸，使灰度的均值与方差和KL变换第一分量图像一致；然后以拉伸过的高分辨率图像代替第一分量图像，经过KL逆变换到原空间。PCA变化法的基本步骤：（1）对多光谱图像进行主成分分析，用相关矩阵求特征值和特征性向量

5、，然后求各主分量；（2）将空间配准后的单波段高分辨率图像与第一主分量图像做直方图配准；（3）用配准后的单波段高分辨率图像代替第一主分量；并将与其余主分量做逆变换得到融合图像。PCA变换法最大的优点是可以应用任意数目的波段，克服了IHS变换只能用于三个波段的缺陷。经过融合后的图像包括了两幅原始图像高空间分辨率和高光谱分辨率特征，保留了原图像的高频信息。利用ENVI进行图像融合 ENVI中提供的融合方法有HSV变换和Brovey变换，这两种方法要求数据具有地理参考或者具有相同的尺寸，RGB输入波段必须为无符号8-bit数据或者从打开的彩色Display中选择。图1 ENVI提供的图像融合办法本次实

6、习我选择了Brovey变换和Gram-Schmidt变换下面是使用Brovey变换操作过程：1选择File-open,打开TM-30m.img和blrd_sp.img文件。2选择Toolbox-Image Sharping-Gram-Color Normalized（Brovery）,在Select Input RGB中选择波段列表中的波段，我分别选择了TM-30m的123波段和321波段来做对比。波段选择后点击OK。3在High Resolution Input File中选择blrd_sp.img文件4设置输出路径，点击OK。图2 Brovey变换图3 Brovey变换细节图使用Gram-

7、Schmidt Pan Sharping融合1选择File-open,打开TM-30m.img和blrd_sp.img文件。2选择Toolbox-Image Sharping-Gram-Schmidt Pan Sharping,Low Spatial Resolution选择TM-30m,High Spatial Resolution选择blrd_sp.img。3设置Pan Sharping Parameters面板，参数设为默认值，设置输出路径，点击OK得到融合后图像。图4 Schmidt Pan Sharping变换实习二：几何校正2 3 3.1 【实习内容】1） 巩固影像几何校正的基本原

8、理。2） 熟练掌握影像几何校正处理流程和操作步骤。3） 掌握几何校正精度评估原理和方法。3.2 【实习步骤】（1）实习流程图：（2）数据待校正影像3.3 【实习过程】（1）几何校正的内容遥感成像的时候，由于飞行器的姿态、高度、速度以及地球自转等因素的影响，造成图相对于地面目标发生几何畸变，这种畸变表现为像元相对于地面目标的实际位置发生挤压、扭曲、拉伸和偏移等，针对几何畸变进行的误差校正就叫几何校正。遥感影像图的几何校正目前有3种方案 ，即系统校正、利用控制点校正以及混合校正。遥感数据接收后，首先由接收部门进行校正，这种校正叫系统校正 (又叫几何粗校正)，即把遥感传感器的校准数据 、传感器的位置

9、 、卫星姿态等测量值代人理论校正公式进行几何 畸变校正；而用户拿到这种产品后，由于使用目的不同或投影及比例尺不同，仍旧需要做进一步 的几何校正，这就需要对其进行几何精校正即利用地面控制点GCP(GroundContro-Point）对因其他因素引起的遥感图像几何畸变进行纠正。混合校正则是由一般地面站提供的遥感CCT已经完成了第一阶段的几何粗校正，用户所要完成的仅仅是对图像做进一步的几何精校正。 （2）几何校正的过程误差太大（3）利用ENVI进行几何校正本次实习采取的是ENVI里面的image to image的方式进行几何校正。1打开并显示图像文件打开ENVI软件，将几何校正前与几何校正后的两

10、幅图通过不同的窗口显示出来。#1为几何校正前的图像，#2为几何校正后的图像。可将图像通过link display简单地连接起来。2启动几何校正模块 点击工具栏中Map，找到registration，选择select GCPs:Image to Image。在弹出的对话框中选择display #2作为base image，也就是基准图像，而display #1作为warp image待校正的图像。选好后点击OK，进入采集地面控制点。图5 选择参考图像和待校正图像3采集地面控制点 在两个Display中移动方框位置，寻找明显的地物特征点作为输入GCP；图6 采取控制点在ground control

11、 points selection上，单击add point按钮，将当前找到的点收集；图7 选取控制点用同样的方法继续寻找其余的点，当选择控制点的数量达到3时，RMS被自动计算。Ground control points selection上的predict按钮可用来在基准图像显示窗口上面定位一个特征点，单击predict按钮，校正图像显示窗口上会自动预测区域，适当调整一下位置，单击add point按钮，将当前找到的点收集。随着控制点数量的增多，预测的的精度越来越精确。 选择option 里的auto predict，打开自动预测功能，这是在基准图像显示窗口上面定位一个特征点时，校正图像显示

12、窗口上会自动预测； 当选择一定数量的控制点之后（至少3个），可利用自动找点功能； 在ground control points selection上，选择options里的automatically generate points，选择一个匹配波段，单击OK； 在automatic tie point method parameter对话框中，选择可将tie点的数量（number of tie points）设为60；其他默认参数，点击OK；选择ground control points selection中的show list按钮，可看到方才选取的所有控制点的相关信息列表；图8 检查GCP精

13、度 可通过选择列表中option里的order point by error由高到低排列RMS值来查看选择点的信息。对于RMS过高的点，可以将其直接删除，或是通过在zoom窗口中进行选择点位置的微调来使RMS值降低；4校正结果的输出在ground control points selection中选择option里的warp file（as image tomap），选择校正文件（TM文件）；在校正参数对话框中，默认投影参数和像元大小与基准图像一致，投影参数不变，点击OK；图9 几何校正后的影像图10 校正前和校正后实习三：影像分类（一）4 4.1 【实习内容】1） 巩固监督分类的基本原理；2

14、） 熟练掌握监督分类的流程和操作步骤；3） 掌握监督分类的影像统计方法。4.2 【实习步骤】4.3 【实习过程】（1）监督分类内容监督分类 （supervised classification）又称训练场地法，是以建立统计识别函数为理论基础，依据典型样本训练方法进行分类的技术。即根据已知训练区提供的样本，通过选择特征参数，求出特征参数作为决策规则，建立判别函数以对各待分类影像进行的图像分类，是模式识别的一种方法。要求训练区域具有典型性和代表性。判别准则若满足分类精度要求，则此准则成立；反之，需重新建立分类的决策规则，直至满足分类精度要求为止。（2）监督分类的过程选择训练样本计算判别函数分类后处

15、理精度评定特征判别（3）用ENVI进行监督分类ENVI提供了多种监督分类的算法：平行六面体法、最大似然法、最小距离法、马氏距离法等。本次实习采用了最大似然法，步骤如下：1选取roi：在影像窗口选择overlay-region of interest.-在弹出的窗口中选择zoom，在zoom窗口中绘制roi；2检查roi：在ROI Tool窗口中选择option-computer ROI Separability.-选择“裸土+水田”文件-选择所有的ROI进行可分离度计算，点击OK；图11-1 选择roi图11-2 所选roi属性图12 计算可分离度3如果可分离度大于1.8，则认为roi选取合格

16、，否则需要重新选择roi。将选择好的roi保存到指定路径。4选择Classification -supervised - maximuum likelihood -选择待分类影像 -选择所有的roi，设置参数，点击OK。图13 参数设置5修改错分区域：在display窗口中选择overlay-classification-选择分类后的影像；在Interactive Class Tool面板中，选择Edit Mode: Polygon Add to Class在Interactive Class Tool面板中，鼠标左键单击Unclassified前面的方型色块，让Unclassified类别处

17、于激活状态；选择一个编辑窗口：Image，在Image窗口中绘制多边形，多边形以内的像元全部归于Unclassified一类；选择要合并的类并保存。图14 最大似然法分类后结果6影像统计： 选择Classification - Post Classification -class statistics-选择分类后的图像和原始图像；设置统计内容，勾选Basic State,Histograms以及Covariance,点击OK。图15 设置统计内容图16 分类信息直方图实习四：影像分类（二）5 5.1 【实习内容】1） 巩固非监督分类的基本原理；2） 熟练掌握非监督分类的流程和操作步骤；3） 掌

18、握非监督分类的精度评定方法。5.2 【实习步骤】5.3 【实习过程】（1）非监督分类的内容非监督分类是以不同影像地物在特征空间中类别特征的差别为依据的一种无先验（已知）类别标准的图像分类，是以集群为理论基础，通过计算机对图像进行集聚统计分析的方法。根据待分类样本特征参数的统计特征，建立决策规则来进行分类。而不需事先知道类别特征。把各样本的空间分布按其相似性分割或合并成一群集，每一群集代表的地物类别，需经实地调查或与已知类型的地物加以比较才能确定。是模式识别的一种方法。一般算法有：回归分析、趋势分析、等混合距离法、集群分析、主成分分析和图形识别等。（2）非监督分类的过程选择分类器设定参数分类后处

19、理精度评定（3）用ENVI进行非监督分类ENVI提供了两种非监督分类的算法：ISODATA和K-Measns。本次实习我选择的是K-Measns分类，下面是具体步骤：1选择File-Open Image File,打开“裸土+水田.tif”2选择Classification-Unsupervised-K-Measns3选择“裸土+水田.tif”图 选择待分类文件4设置参数，最大迭代次数为10次，一般非监督分类的分类数量比最终分类数量多2到3倍为宜，这样有利于提高分类精度。图 参数设置图 分类后的影像5将分类后的图像对照原始图像再次进行合并类别：打开Class color table，对应原始图

20、像分别对各类命名；选择Classification-Post Classification-Combine Classification-选择k-means分类后的图像；选择合并类别，点击OK，设置输出路径。 图 合并类别图 合并之后的图像6精度评定：本次选择的是用混淆矩阵来做精度评定，由于没有标准分类的图像，只能采取ROI的方式进行精度评定。选择Classification - Post Classification - Confusion Matrix - Using Ground Truth ROI；选择合并后的影像和ROI文件，如果各类名相同ENVI会自动匹配，如果不同则需要手动匹配。

21、匹配好后设置输出路径，点击OK。图 匹配ROI和分类结果图 混淆矩阵结果心得与体会 本次实习中遇到不少问题，让我认识到了自己的不足。 一开始本来想用ENVI5.1来完成本次实习，但发现自己对ENVI5.1的工具箱并不熟悉，十分艰难的做完图像融合之后，还是选择了使用比较熟悉的ENVI Classic进行操作。几何校正中遇到的最大的问题就是控制点的选取，最开始做的时候由于基准图像分辨率较低，在GCP的选取上有着较大的困难，从而导致选取的RMS值总是偏大。在上网查了资料后发现控制点应该选择具有比较固定位置的地方，比如马路的交叉，屋顶等。于是删除了之前选的河流入口处的点，选择了一些马路交叉处的点，这样

22、整体的RMS值就下降了。但是由于分辨率不够，个别的RMS值还是偏高。监督分类中遇到的困难就是分类错误的类的处理。一开始我的分类精度很低，在不断的查资料和尝试中发现是因为分类后没有对分类错误的类进行处理。然后对着资料自己操作了一遍之后，将错误的分类删除或者合并到其他类里，最后我的精度也有所提高。非监督分类最后的精度在70%左右，效果一般，由于没有高精度分类结果，导致影像精度不高的原因应该是用来评估的ROI的选取有问题，重新选择了一次ROI并根据可分离度对ROI进行了修改，分类精度有所提高，此外，非监督分类的精度也与参数值有关，迭代次数，合并类别都会影响到最后的精度。本次实习也在不断的尝试中结束了，通过本次实习，我更加深入的了解了遥感图像处理的过程，加深了对于遥感图像处理的理解。当然，在不断的尝试中，我更加的熟悉了遥感图像处理软件的操作，基本能独自完成遥感图像处理中的大部分过程。对于遥感这个学科有了更进一步的理解和认识。最后，对于本次实习中给予帮助的老师和同学们，表达我真挚的感谢！