基于纹理信息与神经网络的遥感影像分类.docx
《基于纹理信息与神经网络的遥感影像分类.docx》由会员分享,可在线阅读,更多相关《基于纹理信息与神经网络的遥感影像分类.docx(29页珍藏版)》请在冰豆网上搜索。
基于纹理信息与神经网络的遥感影像分类
基于纹理信息提取与神经网络分类方法的遥感影像分类
基于纹理信息提取与神经网络分类方法的遥感影像分类
1实验目的
1、通过深入研究遥感影像纹理信息的提取原理与方法,针对研究区影像特征,掌握纹理信息对提高影像分类精度的重要作用。
2、了解ENVI中基于神经网络的影像监督分类方法原理,采用纹理信息与神经网络分类相结合的方式,对研究区影像进行分类,与无纹理信息的分类结果进行对比与评价。
3、进一步熟悉影像神经网络分类的原理和方法、训练区选取、分类后处理方法等操作,并根据研究区实际情况进行分类后评价。
2实验原理
基本概念
神经网络指用计算机模拟人脑的结构,用许多小的处理单元模拟生物的神经元,用算法实现人脑的识别、记忆、思考过程应用于图像分类。
神经网络算法是用计算机模拟人类学习的过程,建立输入和输出数据之间联系的方法。
神经网络分类器在遥感影像分类领域得到了普遍的关注。
国内外学者发展了多种形式的神经网络模型和算法,如反向传播网络、模糊神经网络、多层感知网络、Kohonen自组织特征映射网络、Hybrid学习向量分层网络等。
网络的输入和输出节点之间通过隐含层相连,节点之间通过权重连接,因而这种方法可以将多种数据,如纹理信息、地形信息等,方便有效地融合到遥感影像的分类过程中,增强了分类能力。
神经网络是非线性系统,可以在特征空间构造出分类界面比较复杂的子空间,因此对非线性可分的特征子空间尤为有效。
但是,神经网络也存在一定的缺点,如初始权重选择的困难、收敛速度慢、对输入数据的预处理要求高等,对遥感影像分类结果有重要影响。
目前,应用和研究最多的是利用反向传播算法(BP算法)训练权值的多层前馈神经网络。
该网络的学习训练过程由正向传播和反向传播组成,在正向传播过程中,输入信息从输入层经隐含层逐层处理,并传向输出层,若在输出层得不到期望的输出,则输入反向传播,将误差信号沿原路返回,通过修改各层神经元间的权值,达到误差最小。
一般说来,隐含层数目和隐含层结点数难以确定,增加隐含层数目和结点数,可以提高精度,且有效地减少局部极小的概率,但是需要更长的学习时间。
神经网络由于具有较强的非线性逼近能力及自适应、自学习能力,因此可以处理难以用数学模型描述的系统。
对于一个特定的问题,通常很难判断哪种网络是最有效的,因为对于网络类型的选择取决于很多因素,包括问题的复杂程度以及所研究问题的性质、训练样本的多少、网络的结构、权值和偏置值的数目、误差目标、参数取值等。
实际应用中,一般要根据具体问题,对几种网络进行比较,选择较为合适的算法。
原理阐述
人工神经网络(ArtificialNeuralNetworks,ANNs),有时也称为计算机神经网络(ComputationalNeuralNetworks,CNNs),试图模仿生物神经网络的特点,用以解决一系列信息处理、分析和建模的难题。
人工神经网络的主要种类即所谓的前馈网络。
一个典型的前馈神经网络由3个或更多的相互连接的节点构成——一个输入层、一个或多个隐藏的中间层(通常只有1个)和一个输出层。
箭头指示信息流方向,从输入到输出传递信息。
输入层和隐藏层之间的连接可以用权矩阵W来描述,行/列项值
是正或负的实数权重,或如果不存在联系时为零。
同样,隐藏层和输出层的联系也可以被视为一个权矩阵Z,同样由一组权值
构成。
每种情况下,正权值意味着对源节点或输入节点有关的增强作用,而负权值则对应于抑制作用。
一般来说,隐藏层和输出层的所有节点也同样被连接到一个偏置节点上。
例如:
在多层次感知器中,分别有一个偏置节点与每个前馈层连接,n个输入节点,一个带有m个隐藏节点的隐层,和p个输出节点,权重矩阵W和Z的维数分别为
和
。
输入数据根据其
权重矩阵的值来加权并在隐藏层中进行组合(主要是求和)。
之后这种加权求和由已知的一种激活函数g来修改。
可以采用不同的激活函数,但是空间分析中最常用的是logistic或sigmoid函数。
这种过程可以由以下节点层来阐述:
g
如果我们把输入集表示成一个数据矩阵X,把输出集表示成输出矩阵Y,一个人工神经网络就是一个由X到Y的映射f:
XY。
典型的,这种映射是一种非线性加权和——事实上非线性是ANN应用于大多数问题的一个本质特征。
输入数据被
被进行两次修改,第一次是加权求和,第二次是利用激活函数。
对下一层应用相同的过程,将隐藏层输出值
进行求和,并按照和输入层类似进行可选的修改,从而在输出层产生最终结果
:
和
,最后
故
如果有一个或多个控制数据集,即给定一个输入集i,并且知道目标输出集或者输出矢量t,则可以用该信息来“训练”网络。
该学习过程称为“监督”,因为它表示我们具有解的先验知识,并且在监督训练过程。
这里的控制数据可被分为两大组:
用来训练神经网络的“训练数据”(先前选取的训练样区ROI)和先前未用于神经网络的用于评估网络性能的“测试数据”(测试用到的样区ROI)。
监督学习被广泛应用于模式识别/分类和复杂函数逼近(非线性回归)领域中。
例如:
假设两个加权矩阵的初始值均被随机地设为[0,1]中的统一值。
当值通过网络传播到输出层时,输入数据集或矢量(或训练数据)i将被转换,给出一个输出值的集合或矢量O时,将会与已知的正确值不匹配,即产生误差。
理想状态下,我们想通过选择适当的权重集来最小化这些误差,为了实现这一目标,我们通常采用最小化所有输出节点的平方差总和(或均方根误差RMS),它是一种可以反映出最小二乘回归的方法。
通过误差平方和表达式的微分来达到最小化的目的,使它等于零,并用该结果来确定对初始权重要求的调整。
通常计算出误差平方和均值的平方根(RMSE):
该式提供了一种解的“代价”的标准化方法,目的是最小化该代价函数,最好是通过之前所说的调整权重来使它的值降为零。
调整权重矩阵中的值的过程通常在被称为伴随梯度下降的反向传播过程实现的。
反向传播部分指一步一步的在网络上应用一个后向学习规则——因此第一步是改变Z矩阵中的值,然后改变W矩阵中的值。
然后
矩阵中值的正负增量是基于它们对输出结果中的正、负误差有关的贡献。
由于激活函数g的引入,使得误差信号不是基于输出值和期望输出值之间的差异,而是基于隐藏层值和反向传播隐藏值间的差异,因此,调整起来更加的复杂。
这些反向传播过程相当于一个用来找出均方根局部最小值的梯度下降方法(例如:
在ENVI进行神经网络分类时,可以看到plot图中理想状态RMS值应逐渐降低,而非波动曲折)。
完成调整各层的权重之后,前向传播过程再次进行计算代价函数,现在的代价函数应该小于或等于之前的代价函数。
该过程持续进行,直到代价函数停止降低或已经达到预设的最大迭代次数。
在ENVI中使用NeuralNet选项可以应用一个分层的正向(feed-forward)神经元网络分类。
该技术在进行监督学习时使用标准的后向传播技术(backpropagation)。
用户可以选择所用的隐藏层的数量,也可以在对数和双曲线活化(activation)函数之间选择所需函数。
由于调整节点中的权重可以使输出节点活化与所需的输出结果间的差异达到最小化,因此神经元网络技术利用该方法对发生的事件进行学习。
在学习过程中,误差在网络中后向传播,同时使用递归法调整权重。
也可以使用神经元网络来进行一个非线性分类。
技术路线
3、实验设备与数据与研究区域
硬件:
计算机
软件:
Arcgis10.2、ENVI5.1
数据:
Landset8影像、福州市shp文件
Landset8数据(成像时间:
2013年8月4日2点34分;行列号为119,42;云量覆盖度为1.21%)
传感器:
OLI(OperationalLandImager,陆地成像仪)
TIRS(ThermalInfraredSensor,热红外传感器)
OLI陆地成像仪
序号
波段(um)
空间分辨率(m)
1
0.433-0.453
30
2
0.450-0.515
30
3
0.525-0.600
30
4
0.630-0.680
30
5
0.845-0.885
30
6
1.560-1.660
30
7
2.100-2.300
30
8
0.500-0.680
15
9
1.360-1.390
30
TIRS热红外传感器
10
10.6-11.2
100
11
11.5-12.5
100
地成像仪(OLI)包括9个波段,空间分辨率为30米,其中包括一个15米的全色波段;热红外传感器(TIRS)包括2个热红外波段,空间分辨率为100米。
研究区域:
福州市为福建省省会城市,位于福建省东部、闽江下游沿岸。
介于北纬25°15′~26°39′,东经118°08′~120°31′之间。
东濒台湾海峡,西邻南平市和三明市,南接莆田市,北靠宁德市。
福州是属于典型的海洋性亚热带季风气候,全年冬短夏长,温暖湿润,气候适宜。
福州市的地貌属于典型的河口盆地,是城市中心的发源地。
实验中以福州市的5区(台江区、仓山区、鼓楼区、马尾区和晋安区)为例,不包括琅岐岛。
4、实验过程
4.1数据预处理
4.1.1数据裁剪
利用shp文件对遥感影像进行裁剪
4.1.2辐射校正(包括辐射定标和大气校正)
A.辐射定标
B.大气校正(采用ENVI中的FLAASH大气校正模型)
在下图的中输入影像的中心经纬度坐标;输入影像所在区域的平均高程0.089;输入卫星的过境时间;在4种输人大气模型,这里选择中纬度夏天;输入气溶胶模型,这里选择Urban。
4.2数据融合
在ENVI软件中,选择Transform->ImageSharpening->HSV,在SelectLowSpatialResolutionMultiBandInputFile中选择band8多光谱波段,SpatialSubset为FullScene,SpectralSunset为3个波段
Resampling选择NearestNeighbor,OutputResult选择保存路径
4.3训练样区的选择
4.3.1打开GoogleEarth影像作为训练样区选择的参照
使用Envi的GoogleEarthBridge将遥感影像和街区矢量图层导入GoogleEarth,以便在GoogleEarth寻找参考影像时更快定位。
在ENVI软件中,选择Spectral->SPEARTools->GoogleEarthBridge。
在GoogleEarthBridge窗口中,选择添加的遥感影像HVS1,点击NEXT进入下一步,设置遥感影像显示的参数,选择6、5、4波段进行彩色显示,点击NEXT进入下一步,同样添加矢量图层,如图44,保存路径为GoogleEarthBridge.kml。
4.3.2建立兴趣区
打开兴趣区模版,选择Overlay->RegionofInterest,打开ROITool
在ROITool窗口中,我们可以对ROIName和Color进行编辑,双击ROIName的区域,可以对训练样区进行名字的修改,右击Color可以选择训练样区的颜色,windows选择Zoom窗口,即只能在Zoom中选择训练样区;在#1Zoom窗口中,我们鼠标左键选择训练样区,右击确定,第二次右击即保存该训练样区。
在ROITool窗口中,点击NewRegion建立其新的ROI样区,Goto可以跳转各个训练样区,就可以对其进行修改,在实验中,我们将分成5类,进行