基于ENVI的土壤高光谱影像监督分类方式比较研究.docx
- 文档编号:11894276
- 上传时间:2023-04-08
- 格式:DOCX
- 页数:27
- 大小:1.47MB
基于ENVI的土壤高光谱影像监督分类方式比较研究.docx
《基于ENVI的土壤高光谱影像监督分类方式比较研究.docx》由会员分享,可在线阅读,更多相关《基于ENVI的土壤高光谱影像监督分类方式比较研究.docx(27页珍藏版)》请在冰豆网上搜索。
基于ENVI的土壤高光谱影像监督分类方式比较研究
摘要
基于监督分类方式在遥感影像分类中的普遍应用,介绍了四种ENVI提供的比较经常使用的、算法简单、计算时刻较短的四种分类方式。
对同一土壤剖面影像运用这四种方式进行分类,并对分类结果进行了对照,分析了这四种方式分类精度之间的不同,目的是为了更好地将土壤剖面高光谱遥感影像中的无效数据剔除,为后续的研究提供数据基础。
关键词:
PCA变换;平行六面体法;最小距离法;马氏距离法;最大似然法
Abstract
Widelyusedinremotesensingimageclassificationbasedonsupervisedclassificationmethod,thispaperintroducesfourkindsofENVIprovidesacommon,simplealgorithm,calculationoffourkindsofclassificationmethodsforashorttime.Theclassificationofthesamesoilprofileimagebyusingthesefourmethods,andtheclassificationresultsarecompared,analyzedthedifferencesbetweentheclassificationaccuracyofthefourpurposeistobettertheinvaliddatatoeliminatethesoilprofileofhyperspectralremotesensingimage,providethedatabasisforthesubsequentresearch.
KeyWords:
PCATransform;Parallelpiped;MinimumDistance;MahalanobisDistance;LikelihoodClassification
1.前言
研究问题的由来
在土壤剖面制备进程中,由于土壤的含水量会对土壤的光谱产生阻碍,掩盖其他土壤属性对土壤光谱曲线的阻碍,因此要将搜集的土壤剖面在室内通风晾干,风干时土壤失水产生块状凝聚,有裂痕显现,另外土壤中的大颗粒或凹凸不平的地址会产生阴影,这些裂痕和阴影区域不具有土壤反射光谱特点,在成像光谱仪取得的土壤剖面高光谱影像中属于无效数据,因此需要把影像中的这些区域剔除,有利于后续的高光谱影像的研究分析。
文献综述
高光谱遥感影像具有光谱分辨率高、信息量大、图谱合一的特点,具有传统遥感技术无法比拟的优势。
可是,高光谱数据的高光谱分辨率也带来了大量数据的冗余,在将其运用于各领域之前,必需进行必要的数据处置。
因此,如何有效的提取高光谱数据信息成为当前研究的热点之一(李静,2021)
遥感作为搜集地球数据及其转变信息的重要技术手腕,活着界范围内和我国的许多政府部门、科研单位和公司取得普遍的应用,不同领域遥感影像的应用对遥感影像处置技术提出了不同的要求,影像分类是影像处置的重要环节,因此具有专门大的研究意义,随着各类新理论新方式的接踵涌现,遥感影像存在多种分类方式,通过几种经常使用监督分类方式的比较发觉,每种分类方式都有最适合应用的范围和自身的局限性,没有一种是最普遍最正确的方式,因此必需灵活应用,综合应用多种分类方式,而且与其他影像处置技术结合起来实现最大精度的分类。
(闫琰等,2020)。
非监督分类是指在分类时对需要分类的地物完全没有加入任何已知的信息,而仅仅依托地物的自然特性,非监督分类的实质是聚类分析法,由于在非监督分类进程中没有类别先验知识的阻碍,因此无法判定分类的结果别离代表哪一类实际地物,而且很难保证所有的特点是被分类别最具有判定力的特点。
因此分类精度不够理想,而且关于庞大的高光谱数据来讲,分类速度较慢。
但该方式不需要对分类区域有普遍的了解,而且人为误差的机遇少。
因此能够与监督分类结合利用(李静,2021)。
自20世纪70年代以来,随着数字成像技术和运算机技术的飞速进展,影像处置、分类与信息提取的算法层出不穷。
高光谱影像的分类与多光谱有专门大不同,目前有沿用传统针对多光谱的分类算法,也进展了针对高光谱的算法。
在高光谱遥感分类及信息提取领域各类算法层出不穷,基于光谱特点的分类是高光谱的特色(林娜,2020)。
传统遥感影像分类多是基于统计特点的。
如平行管道、最大似然、最小距离、马氏距离等,原先针对多光谱数据的处置手腕在高光谱的数据分析中多有利用。
如XiuPingJia(2002)就采纳最大似然、主成份变换等对高光谱数据进行了分析。
佘红伟,张艳宁(2020)研究了一种无监督高光谱图像分类算法.利用传统方式处置高光谱影像存在很多问题如运算速度过慢、无法取得足够多的训练样本、维数灾难(Hughes)等问题。
高光谱遥感图像中端元的提取是明白得高光谱数据,继而对数据进行进一步分析的前提条件。
Boardman(1993)进展了凸面几何学分析,Boardman、Kruse、Green(1995)进展了纯象元指数(PPI)端元提取算法。
N-Findr(Winter,1999)算法寻觅一组像元,它们所组成的单纯形体具有最大的体积,是一种全自动的端元选取算法。
最近几年来陈伟,陈伟,余旭初(2020)等利用粒子群优化(PSO)技术,基于凸面几何学理论,设计了一个新的端元提取算法。
李姗姗,田庆(2020)提出了一种高光谱遥感图像的端元递进提取算法。
吴波,张良培,李平湘(2005)依照迭代误差分析思路,结合端元存在的空间信息,自动提掏出端元光谱。
由于高光谱遥感数据具有高维的光谱特点,基于光谱特点的分类及信息提取是一个重要的研究方向,是高光谱遥感数据特有的识别方式。
代表性的有光谱角度填图(SAM)(Kruse,1993),交叉相关光谱匹配技术(CCSM),二值编码匹配(Mazer,1988),波谱特点拟合(SFF)等。
最近几年来也有些学者对此类算法进行了
深切研究,如周源,方圣辉,李德仁(2020)提出了一种新的光谱匹配算法——光谱角灵敏丛林方式。
王毅,张良培,李平湘(2007)提出了一种基于局部自动搜索和光谱匹配技术的监督分类训练样本的纯化方式。
目前许多新的人工智能算法也被引入高光谱图像分类中,如自组织映射、人工神经网络、遗传算法、支持向量机(SupportVectorMachine,SVM)(刘江永,2020)等。
徐宏根、马洪超、李德仁(2007)将自组织映射(SOM)神经网络与混合像元分解相结合对高光谱影像进行分类研究,取得了较好的分类成效。
官莉,刘旸,张雪慧(2020)将人工神经网络算法应用在红外高光谱资料反演大气温度廓线中,取得了专门好的成效。
ZhuoL(2020)将遗传算法应用到高光谱图像波段选
择和图像分类中。
冯静、舒宁(2020)将改良型遗传算法和支持向量机结合进行波段选择和图像分类取得了中意的成效。
吴波、张良培、李平湘(2006)基于支持向量回归进行高光谱混合像元分解。
骆剑承(2002)等将SVM应用于遥感影像空间特点提取与分类,取得了比较高的分类精度。
牛鹏(2020)结合高光谱数据特性,改良了SVM,将其应用到高光谱遥感图像分类中。
研究目的
为了将土壤剖面高光谱影像的裂痕和阴影区域剔除,需要将土壤剖面上的裂痕和阴影与正常突然分成不同的类别,本文采纳几种不同的监督分类法对土壤剖面高光谱影像分类,并比较不同分类方式的精度不同,为后续的研究提供数据基础。
2.材料与方式
材料
本文利用的是Headwall成像光谱仪(Hyperspec™NIR,900-1700nm)取得的风干土壤剖面的高光谱影像,而且是通过几何校正和裁剪等预处置以后的影像,在后面的实验进程中称为原始影像。
高光谱传感器即成像光谱仪,能够在电磁波谱的紫外、可见光、近红外和中红外区域,以数十至数百个持续且细分的光谱波段对目标区域同时成像,由于成像光谱系统取得的持续波段宽度一样在10nm之内,因此高光谱数据能以足够的光谱分辨率区分出那些具有诊断性光谱特点的地表物质,定量地分析地球表层生物、物理、化学进程与参数(董连凤,2007)。
由于土壤反射光谱特点综合反映了土壤各类理化性质,因此高光谱遥感技术在土壤剖面的研究中有专门大的应用价值(张美琴,2021)。
方式
为了隔离噪声和减少数据集的维数,需要对原始影像进行PCA变换,主成份分析(PCA)通过利用PrincipalComponents选项生成互不相关的输出波段,主成份(PC)波段是原始波谱波段的线性合成,它们之间是互不相关的。
能够计算输出主成份波段(与输入的波谱波段数相同)。
第一主成份包括最大的数据方差百分比,第二主成份包括第二大的方差,以此类推,最后的主成份波段由于包括很小的方差(大多数由原始波谱的噪声引发),因此显示为噪声。
遥感影像的分类方式能够分为非监督分类和监督分类,非监督分类不需要人工搜集地物样本点数据,多是通过聚类的方式来自动分类,要紧有isodata,k均值等。
监督分类是需要学习训练的分类方式,如最大似然分类,人工神经网络分类,即是需要事前为每类地物在遥感影像上搜集样本数据,以后通过学习训练进程才来分类;整体来讲,监督分类的成效要优于非监督分类。
本文选用ENVI软件提供的算法简单、计算时刻较短的,在遥感影像监督分类中经常使用的四种分类方式:
平行六面体法、最小距离法、马氏距离法、最大似然法。
监督分类进程
PCA变换
对原始影像xl1n进行主成份变换,具体方式是在ENVI主菜单上选择Transform>PrincipalComponents>ForwardPCRotation>ComputeNewStatisticsandRotate,输入影像选择xl1n,取得的结果存为pc,然后在ENVI主菜单项选择择BasicTools>ResizeData(Spatial/spectral),选取信息量最大的前三个波段,存为pcResize。
选择训练样本
打开影像,File→Open Image File,选择文件,在可用波段列表中,选择波段48,默许选择GrayScale单项选择按钮,再点击LoadBand按钮,打开一幅灰度影像,如以下图:
图1原始影像
Figure1Theoriginalimage
那个波段的影像纹理比较清楚,能够看出土壤剖面上有很多的裂痕和阴影,在图上点击鼠标右键,在快捷菜单上选择ZProfile(Spectrum)...,打开对话框,点击灰度图上的像素点,能够看到对话框中显示出它的光谱曲线,选择不同位置的像素点查看光谱曲线能够发觉,裂痕和阴影处的光谱曲线和正常土壤的有专门大不同,如以下图:
图2正常土壤
Figure2normalsoil
图3正常土壤的波谱剖面曲线图
Figure3Spectralprofilegraphofnormalsoil
图4裂缝和阴影
Figure4cracksandshadows
图5裂缝和阴影的波谱剖面曲线图
Figure5Spectralprofilegraphofcracksandshadows
通过对照图可发觉,裂痕或阴影中的像素的光谱曲线形状与正常土壤的不同专门大,而且在波长1400nm左右没有水分吸收谷。
利用感爱好区(ROI)工具来选择训练样本,在主影像菜单栏当选择Overlay>Region of Interest,打开ROI Tool对话框,在对话框菜单上选择ROTType>Point,Window后面的单项选择按钮当选择Zoom,然后新建两个感爱好区,默许第一个颜色属性为Red,第二个为Green,在区域名称之前点击,选中第一个区域,然后在影像上的裂痕或阴影处寻觅颜色较深的像素点,在Zoom窗口点击,添加点,选的点要尽可能均匀散布,所有的符合特点的区域都要选择,然后在ROI Tool对话框当选中第二个区域,在影像上颜色比较均匀的区域选择颜色较浅的点,取得的结果如图6,最后把选好的ROI保留为xl1n,将ROITool对话框关闭。
图6ROI的选取
Figure6ROIselection
用不同的分类方式进行监督分类
平行六面体(Parallelpiped)法
依照训练样本的亮度值形成一个n维的平行六面体数据空间,其他像元的光谱值若是落在平行六面体任何一个训练样本所对应的区域,就被划分其对应的类别中。
平行六面体的尺度是由标准差阈值所确信的,而该标准差阈值那么是依照所选类的均值求出。
在ENVI主菜单项选择择Classification>Supervised>Parallelepiped,打开ClassificationInputFile对话框,选择pcResize作为输入影像,点击确信,再弹出的对话框上点击SelectAllItems,然后选择输出文件的位置和名称,其它参数默许,如图,最后计算出的分类结果影像如图8。
图7分类参数选择
Figure7Classificationparameters
最小距离(MinimumDistance)法
利用训练样本数据计算出每一类的均值向量和标准差向量,然后以均值向量作为该类在特点空间中的中心位置,计算输入影像中每一个像元到各类中心的距离,到哪一类中心的距离最小,该像元就归入到哪一类。
在ENVI主菜单上选择Classification>Supervised>MinimumDistance,打开ClassificationInputFile对话框,选择pcResize作为输入影像,点击确信,再弹出的对话框上点击SelectAllItems,然后选择输出文件的位置和名称,其它参数默许,最后计算出的分类结果影像如图9。
马氏距离(MahalanobisDistance)法
计算输入影像到各训练样本的马氏距离(一种有效的计算两个未知样本集的相似度的方式),最终统计马氏距离最小的,即为此类别。
在ENVI主菜单上选择Classification>Supervised>MahalanobisDistance,打开ClassificationInputFile对话框,选择pcResize作为输入影像,点击确信,再弹出的对话框上点击SelectAllItems,然后选择输出文件的位置和名称,其它参数默许,如图,最后计算出的分类结果影像如图10。
最大似然(LikelihoodClassification)法
最大似然法分类假定每一个波段每一类统计呈均匀散布,并计算给定像元属于一特定类别的可能性。
除非选择一个可能性阈值,所有像元都将参与分类。
每一个像元被归到可能性最大的那一类里。
在ENVI主菜单上选择Classification>Supervised>MaximumLikelihood,打开ClassificationInputFile对话框,选择pcResize作为输入影像,点击确信,再弹出的对话框上点击SelectAllItems,然后选择输出文件的位置和名称,其它参数默许,如图,最后计算出的分类结果影像如图11。
图8平行六面体法
Figure8Parallelpiped
图9最小距离法
Figure9MinimumDistance
图10马氏距离法
Figure10MahalanobisDistance
图11最大似然法
Figure11LikelihoodClassification
图12原始影像
Figure12Theoriginalimage
3.结果与分析
本文利用ENVI软件中的监督分类模块对通过PCA变换和裁剪的同一影像别离用平行六面体法、最小距离法、马氏距离法、最大似然法四种分类方式进行监督分类,同时这四种分类都采纳了相同的训练样本,有利于比较不同分类方式的不同。
四种分类结果简单目视判定
四种分类结果如图8,9,10,11,图12是原始影像。
通过对照图8,9,10,11,12,平行六面体法的分类结果把影像的大部份像素分类到裂痕和阴影里面,和原始影像明显不符合,为了更好地比较这四种分类结果,能够在4幅分类结果影像和xl1n原始影像打开的情形下,利用LinkDisplays...工具将5幅影像联系起来,使缩放窗口显示的区域相同,移动主窗口和转动窗口,结合原始影像,对照四种分类结果的细节不同,图13为缩放窗口中某个区域的分类结果细节对照图,上面一排从左到右别离为平行六面体法、最小距离法、马氏距离法的分类结果细节图,下面一排左侧的图是最大似然法的分类结果细节图,右边的图是原始影像细节图,对照能够发觉,平行六面体法不仅把大部份影像分类到了裂痕和阴影类中,还有一些像素点没有被分类,分类成效很差,另外三种分类方式看上去和原始影像都比较符合,其中不同不大,只是有几个像素的分类不同,对照最小距离法和马氏距离法的分类结果,最大似然法将较多的像素分类到了裂痕和阴影类中。
图13某个区域的分类结果细节对比图
Figure13Classificationresultsdetailscontrastfigureofaregion
用混淆矩阵进行分类精度评判
为了定量的评判四种分类结果的不同,利用confusionmatrix工具能够把分类结果的精度显示在一个混淆矩阵里(用于比较分类结果和地表真实信息)。
地面真实影像能够是另一幅分类影像,或是依照地面真实测量生成的影像(也能够是感爱好区)。
以平行六面体法分类结果为例,计算混淆矩阵,另外三种分类结果的计算方式类似。
计算分类结果混淆矩阵
将分类结果影像打开,在影像主窗口选择Tools>RegionofInterest>RestoreSavedROIFile,在打开的消息框中点击确信。
在ENVI主菜单上选择Classification>PostClassification>ConfusionMatrix>UsingGroundTruthROIs,在ClassificationInputFile对话框当选择分类结果影像Parallelepiped,点击确信,打开如图14所示的对话框:
图14匹配类参数
Figure14Matchingparameters
在matchclassesparameters对话框中的两个列表当选择所要匹配的名称,点击addcombination按钮,把地表真实感爱好区与最终分类结果相匹配。
类别之间的归并将显示在对话框底部列表中,若是地表真实影像感爱好区中的类别与分类影像中的类别名称相同,他们将自动匹配。
点击OK,在打开的对话框中利用默许值,点击确信,生成份类结果的精度表,如图15.
图15平行六面体法分类结果的精度图
Figure15TheaccuracyofclassificationresultsofParallelpiped
其它三种分类结果的精度图如图16,17,18.
图16最小距离法分类结果的精度图
Figure16Theaccuracyofclassificationresultsofminimumdistance
图17马氏距离法分类结果的精度图
Figure17TheaccuracyofclassificationresultsofMahalanobisDistance
图18最大似然法分类结果的精度图
Figure18TheaccuracyofclassificationresultsofLikelihoodClassification
混淆矩阵中的几项评判指标
整体分类精度(OverallAccuracy)
整体分类精度:
等于被正确分类的像元总和除以总像元数。
地表真实影像或地表真实感爱好区限定了像元的真实分类。
被正确分类的像元沿着混淆矩阵的对角线散布,它显示出被分类到正确地表真实分类中的像元数。
像元总数等于所有地表真实分类中的像元总和。
如图15平行六面体法分类结果的整体分类精度为%。
Kappa系数(KappaCoefficient)
是另外一种计算分类精度的方式,通过把所有地表真实分类中的像元总数乘以混淆矩阵对角线的和,再减去某一类中地表真实像元总数与该类中被分类像元总数之积对所有类别求和的结果,再除以总像元数的平方减去某一类中地表真实像元总数与该类中被分类像元总数之积对所有类别求和的结果取得的。
如图15平行六面体法分类结果的Kappa系数为.
混淆矩阵(像元数)
混淆矩阵通过将每一个地表真实像元的位置和分类相较较计算的。
混淆矩阵的每一列代表了一个地表真实分类,每一列中的数值等于地表真实像元在分类影像中对应于相应类别的数量。
例如,在图15中groundtruth(pixels)表中,正常土壤的地面真实栏中显示该类地表真实值为479个像元,分类能对其中的213个像元进行正确分类,有4个像元未参加分类,有262个像元被归为裂痕和阴影。
混淆矩阵(百分比):
groundtruth(percent)表显示了每一个地表真实分类中相应的分类结果所占的百分比,该百分比值等于每一个地表真实栏里的像元数除以一个给定的地表真实类中的像元总数。
错分误差(Commission)
是指被分为用户感爱好的类事实上属于另一类的像元。
错分误差显示在混淆矩阵的行里。
例如,在图15中,裂痕和阴影类共有地表真实像元515个,其中253个分类是正确的。
262个是其他像元被误以为裂痕和阴影的(262是混淆矩阵中裂痕和阴影一行其他类的和)。
分类不正确的像元数与参与分类的地表真实像元总数的比例确实是错分误差:
262/515=%。
漏分误差(Omission)
本属于地表真实分类,但没有被分类器分到相应类别中的像元数。
漏分误差显示在混淆矩阵的列中。
例如,在图15中,裂痕和阴影共有257个真实像元,253个正确分类,4个是被误分为其他类型的(裂痕和阴影一列里其他类的和)。
被误分类的像元占某一类别的地表真实分类的像元总数的比例:
4/257=%。
制图精度()
是指分类器将整个影像的像元正确分为A类的像元数(对角线值)与A类真实参考总数(混淆矩阵中A类列的总和)的比率。
例如,在图15中,groundtruth(pixels)表中,正常土壤共有479个地表真实像元,其中213个正确分类,因此制图精度是213/479=%。
用户精度()
是指正确分到A类的像元总数(对角线值)与分类器将整个影像的像元分为A类的像元总数(混淆矩阵中A类行的总和)比率。
例如,在图15中,groundtruth(pixels)表中,分类器将515个像元归到裂痕和阴影一类中,但只有253个像元是正确归类的,因此用户精度是:
253/515=%。
四种分类方式精度比较
如表1,本文选取了几个评判标来比较不同的分类方式的精度,平行六面体分类方式的整体分类精度较低,Kappa系数很低,而且裂痕和阴影的错分误差达到50%以上,同时正常土壤的漏分误差达到50%以上,说明有一半以上的正常土壤被分到裂痕和阴影类,分类成效较差。
其他三种分类方式的整体分类
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 ENVI 土壤 光谱 影像 监督 分类 方式 比较 研究