基于纹理信息与神经网络的遥感影像分类Word文件下载.docx

1、目前, 应用和研究最多的是利用反向传播算法（ BP算法）训练权值的多层前馈神经网络。该网络的学习训练过程由正向传播和反向传播组成, 在正向传播过程中, 输入信息从输入层经隐含层逐层处理, 并传向输出层,若在输出层得不到期望的输出, 则输入反向传播, 将误差信号沿原路返回, 通过修改各层神经元间的权值, 达到误差最小。一般说来, 隐含层数目和隐含层结点数难以确定, 增加隐含层数目和结点数, 可以提高精度, 且有效地减少局部极小的概率, 但是需要更长的学习时间。神经网络由于具有较强的非线性逼近能力及自适应、自学习能力, 因此可以处理难以用数学模型描述的系统。对于一个特定的问题,通常很难判断哪种网络

2、是最有效的, 因为对于网络类型的选择取决于很多因素, 包括问题的复杂程度以及所研究问题的性质、训练样本的多少、网络的结构、权值和偏置值的数目、误差目标、参数取值等。实际应用中, 一般要根据具体问题,对几种网络进行比较, 选择较为合适的算法。原理阐述人工神经网络（Artificial Neural Networks，ANNs），有时也称为计算机神经网络（Computational Neural Networks，CNNs），试图模仿生物神经网络的特点，用以解决一系列信息处理、分析和建模的难题。人工神经网络的主要种类即所谓的前馈网络。一个典型的前馈神经网络由3个或更多的相互连接的节点构成一个输入层

3、、一个或多个隐藏的中间层（通常只有1个）和一个输出层。箭头指示信息流方向，从输入到输出传递信息。输入层和隐藏层之间的连接可以用权矩阵W来描述，行/列项值是正或负的实数权重，或如果不存在联系时为零。同样，隐藏层和输出层的联系也可以被视为一个权矩阵Z，同样由一组权值构成。每种情况下，正权值意味着对源节点或输入节点有关的增强作用，而负权值则对应于抑制作用。一般来说，隐藏层和输出层的所有节点也同样被连接到一个偏置节点上。例如：在多层次感知器中，分别有一个偏置节点与每个前馈层连接，n个输入节点，一个带有m个隐藏节点的隐层，和p个输出节点，权重矩阵W和Z的维数分别为和。输入数据根据其权重矩阵的值来加权并在

4、隐藏层中进行组合（主要是求和）。之后这种加权求和由已知的一种激活函数g来修改。可以采用不同的激活函数，但是空间分析中最常用的是logistic或sigmoid函数。这种过程可以由以下节点层来阐述：g如果我们把输入集表示成一个数据矩阵X，把输出集表示成输出矩阵Y，一个人工神经网络就是一个由X到Y的映射f：X Y。典型的，这种映射是一种非线性加权和事实上非线性是ANN应用于大多数问题的一个本质特征。输入数据被被进行两次修改，第一次是加权求和，第二次是利用激活函数。对下一层应用相同的过程，将隐藏层输出值进行求和，并按照和输入层类似进行可选的修改，从而在输出层产生最终结果：和，最后故如果有一个或多个控

5、制数据集，即给定一个输入集i，并且知道目标输出集或者输出矢量t，则可以用该信息来“训练”网络。该学习过程称为“监督”，因为它表示我们具有解的先验知识，并且在监督训练过程。这里的控制数据可被分为两大组：用来训练神经网络的“训练数据”（先前选取的训练样区ROI）和先前未用于神经网络的用于评估网络性能的“测试数据”（测试用到的样区ROI）。监督学习被广泛应用于模式识别/分类和复杂函数逼近（非线性回归）领域中。假设两个加权矩阵的初始值均被随机地设为0,1中的统一值。当值通过网络传播到输出层时，输入数据集或矢量（或训练数据）i将被转换，给出一个输出值的集合或矢量O时，将会与已知的正确值不匹配，即产生误差

6、。理想状态下，我们想通过选择适当的权重集来最小化这些误差，为了实现这一目标，我们通常采用最小化所有输出节点的平方差总和（或均方根误差RMS），它是一种可以反映出最小二乘回归的方法。通过误差平方和表达式的微分来达到最小化的目的，使它等于零，并用该结果来确定对初始权重要求的调整。通常计算出误差平方和均值的平方根（RMSE）:该式提供了一种解的“代价”的标准化方法，目的是最小化该代价函数，最好是通过之前所说的调整权重来使它的值降为零。调整权重矩阵中的值的过程通常在被称为伴随梯度下降的反向传播过程实现的。反向传播部分指一步一步的在网络上应用一个后向学习规则因此第一步是改变Z矩阵中的值，然后改变W矩阵中

7、的值。然后矩阵中值的正负增量是基于它们对输出结果中的正、负误差有关的贡献。由于激活函数g的引入，使得误差信号不是基于输出值和期望输出值之间的差异，而是基于隐藏层值和反向传播隐藏值间的差异，因此，调整起来更加的复杂。这些反向传播过程相当于一个用来找出均方根局部最小值的梯度下降方法（例如：在ENVI进行神经网络分类时，可以看到plot图中理想状态RMS值应逐渐降低，而非波动曲折）。完成调整各层的权重之后，前向传播过程再次进行计算代价函数，现在的代价函数应该小于或等于之前的代价函数。该过程持续进行，直到代价函数停止降低或已经达到预设的最大迭代次数。在ENVI中使用NeuralNet选项可以应用一个分

8、层的正向（feed-forward）神经元网络分类。该技术在进行监督学习时使用标准的后向传播技术（backpropagation）。用户可以选择所用的隐藏层的数量，也可以在对数和双曲线活化（activation）函数之间选择所需函数。由于调整节点中的权重可以使输出节点活化与所需的输出结果间的差异达到最小化，因此神经元网络技术利用该方法对发生的事件进行学习。在学习过程中，误差在网络中后向传播，同时使用递归法调整权重。也可以使用神经元网络来进行一个非线性分类。技术路线3、实验设备与数据与研究区域硬件：计算机软件:Arcgis 10.2 、 ENVI5.1数据： Landset 8影像、福州市shp

9、文件Landset8数据（成像时间：2013年8月4日2点34分；行列号为119，42；云量覆盖度为1.21%） 传感器：OLI（Operational Land Imager ，陆地成像仪） TIRS （Thermal Infrared Sensor，热红外传感器）OLI陆地成像仪序号波段（um）空间分辨率（m）10.433-0.4533020.450-0.51530.525-0.60040.630-0.68050.845-0.88561.560-1.66072.100-2.30080.500-0.6801591.360-1.390TIRS热红外传感器1010.6-11.21001111.5

10、-12.5地成像仪（OLI）包括9个波段，空间分辨率为30米，其中包括一个15米的全色波段；热红外传感器（TIRS）包括2个热红外波段，空间分辨率为100米。研究区域：福州市为福建省省会城市，位于福建省东部、闽江下游沿岸。介于北纬25152639，东经1180812031之间。东濒台湾海峡，西邻南平市和三明市，南接莆田市，北靠宁德市。福州是属于典型的海洋性亚热带季风气候，全年冬短夏长，温暖湿润，气候适宜。福州市的地貌属于典型的河口盆地，是城市中心的发源地。实验中以福州市的5区（台江区、仓山区、鼓楼区、马尾区和晋安区）为例，不包括琅岐岛。4、实验过程4.1 数据预处理4.1.1 数据裁剪 利用s

11、hp文件对遥感影像进行裁剪 4.1.2 辐射校正（包括辐射定标和大气校正）A辐射定标B大气校正（采用ENVI中的FLAASH大气校正模型）在下图的中输入影像的中心经纬度坐标；输入影像所在区域的平均高程0.089；输入卫星的过境时间；在4种输人大气模型，这里选择中纬度夏天；输入气溶胶模型，这里选择Urban。4.2 数据融合 在ENVI软件中，选择Transform - Image Sharpening - HSV ，在Select Low Spatial Resolution Multi Band Input File中选择band8 多光谱波段，Spatial Subset 为Full Sc

12、ene，Spectral Sunset为3个波段Resampling选择Nearest Neighbor，Output Result 选择保存路径4.3 训练样区的选择4.3.1打开Google Earth影像作为训练样区选择的参照使用Envi的Google Earth Bridge将遥感影像和街区矢量图层导入Google Earth，以便在Google Earth寻找参考影像时更快定位。在ENVI 软件中，选择Spectral - SPEAR Tools - Google Earth Bridge。在Google Earth Bridge窗口中，选择添加的遥感影像HVS1，点击NEXT进入下

13、一步，设置遥感影像显示的参数，选择6、5、4波段进行彩色显示，点击NEXT进入下一步，同样添加矢量图层，如图44，保存路径为Google Earth Bridge.kml。4.3.2 建立兴趣区打开兴趣区模版，选择Overlay - Region of Interest ，打开ROI Tool在ROI Tool窗口中，我们可以对ROI Name和Color进行编辑，双击ROI Name的区域，可以对训练样区进行名字的修改，右击Color可以选择训练样区的颜色，windows选择Zoom窗口，即只能在Zoom中选择训练样区；在#1 Zoom窗口中，我们鼠标左键选择训练样区，右击确定，第二次右击即保存该训练样区。在ROI Tool窗口中，点击New Region建立其新的ROI样区，Goto可以跳转各个训练样区，就可以对其进行修改，在实验中，我们将分成5类，进行监督分类，如：water, forest, city，bar