img

官方微信

遥感技术与应用, 2020, 35(4): 749-758 doi: 10.11873/j.issn.1004-0323.2020.4.0749

甘肃遥感学会专栏

基于CNN的不同空间分辨率影像土地覆被分类研究

李宏达,, 高小红,, 汤敏

青海师范大学地理科学学院,青海省自然地理与环境过程重点实验室,高原科学与可持续发展 研究院,青藏高原地表过程与生态保育教育部重点实验室,青海 西宁 810008

Land Cover Classification for Different Spatial Resolution Images from CNN

Li Hongda,, Gao Xiaohong,, Tang Min

School of Geographical Sciences, Qinghai Normal University, Qinghai Province Key Laboratory of Physical Geography and Environmental Process, Academy of Plateau Science and Sustainability, MOE Key Laboratory of Tibetan Plateau Land Surface Processes and Ecological Conservationm, Xi'ning 810008, China

通讯作者: 高小红(1963-),女,陕西白水人,教授,博士生导师,主要从事遥感信息提取与土地覆被变化方面的研究。E⁃mail:xiaohonggao226@163.com

收稿日期: 2019-08-27   修回日期: 2020-06-24   网络出版日期: 2020-09-14

基金资助: 青海省科技厅自然科学基金项目.  2016⁃ZJ⁃907

Received: 2019-08-27   Revised: 2020-06-24   Online: 2020-09-14

作者简介 About authors

李宏达(1995-),男,湖北荆门人,硕士研究生,主要从事遥感应用与地理空间数据分析研究E⁃mail:2395789679@qq.com , E-mail:2395789679@qq.com

摘要

基于卷积神经网络(Convolutional Neural Networks, CNN)和5种不同空间分辨率的遥感影像,对西宁市东部一区域开展土地覆被分类研究,旨在探索CNN在不同空间分辨率下进行影像分类的差异性和对不同地物的提取能力。为提高样本的选择效率,引入了窗口滑动方法进行辅助选样。研究表明5种不同空间分辨率影像的总体分类精度均达89%以上,Kappa系数达0.86以上,分类精度较高。在所涉及的分辨率尺度范围内,空间分辨率越高,CNN分类结果越精细,并能保持较高的分类精度,表明CNN更适合高空间分辨率影像分类;但同时影像空间分辨率越高,地物表现出较高的类内变异性和低类间差异性,分类精度有降低的趋势。相比较而言,SPOT 6影像的分类精度最高,同时窗口滑动是一种有效的样本辅助选择方法。研究对今后同类工作具有一定的借鉴意义。

关键词: CNN ; Landsat-8/Sentinel-2A/SPOT-6/GF-2影像 ; 土地覆被分类

Abstract

Based on convolutional neural networks and five different spatial resolution remote sensing images, the land use/land cover classification study was carried out on a small area in the eastern part of Xining City, aiming at exploring the differences of image classification by CNN with different spatial resolutions and CNN’s ability to extract different features. In order to improve the selection efficiency of the samples, a window sliding method was introduced to assist the samples selection. The research shows that the overall classification accuracy of the five different spatial resolution images is above 89%, the Kappa coefficient is above 0.86. The result further shows that within the resolution scale the higher the resolution, the performance of the CNN classification results for the details is better, and can maintain high classification accuracy, indicating that CNN is more suitable for high spatial resolution images; at the same time, the image spatial resolution is too high, the ground objects exhibit high intra-class variability and low inter-class variability, the classification accuracy tends to decrease. In comparison, CNN has the best classification effect on SPOT 6 images in this study, and window sliding is an effective sample-assisted selection method. This research has certain reference significance for similar research in the future.

Keywords: Convolutional Neural Network ; Landsat-8/Sentinel-2A/SPOT-6/GF-2 images ; Land cover classification

PDF (4662KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

李宏达, 高小红, 汤敏. 基于CNN的不同空间分辨率影像土地覆被分类研究. 遥感技术与应用[J], 2020, 35(4): 749-758 doi:10.11873/j.issn.1004-0323.2020.4.0749

Li Hongda, Gao Xiaohong, Tang Min. Land Cover Classification for Different Spatial Resolution Images from CNN. Remote Sensing Technology and Application[J], 2020, 35(4): 749-758 doi:10.11873/j.issn.1004-0323.2020.4.0749

1 引 言

传统参数化遥感数据分类方法(如最大似然法)计算简单、实现方便,但因仅使用少量的影像特征进行分类,同时假设了数据的参数化分布,泛化能力差,精度通常较低。非参数化机器学习方法,如人工神经网络、支持向量机、决策树和随机森林等具有较强的学习能力及泛化能力,在遥感影像分类中得到广泛应用。有研究表明[1-3],除决策树以外,人工神经网络、支持向量机和随机森林在不同的分类场景下,均有可能是最优的影像分类方法。但总体上它们属于浅层学习,很难有效表达复杂的函数变化,对样本缺少足够的适应性[4]。同时,这些机器学习算法也没有或未充分利用影像的上下文结构信息,尤其对于高空间分辨率影像,其波段数量通常较少,光谱信息有限,表现出较高的类内变异性和类间低差异性,仅用少量的波段信息去区分相似的地物无疑增加了分类的工作难度[5]

20世纪80至90年代,在人工神经网络的基础上,兼顾图像上下文结构信息的卷积神经网络(Convolutional Neural Networks, CNN)开始出现,2006年以来随着计算机技术的发展和深度学习理论的提出,CNN被广泛应用于计算机视觉、自然语言处理等方面,并逐步应用到遥感影像分类领域。Masoud等[5]采用7种不同结构的CNN对RapidEye 5 m多光谱影像进行土地覆被分类,总体分类精度最高为96.17%,同时所有CNN的分类结果均优于支持向量机和随机森林方法。孟祥锐等[6]基于CNN对GF-2融合后的1 m多光谱数据进行地物覆被分类,取得了87.58%的总体分类精度和0.86的Kappa系数,验证了CNN对湿地群落精细划分的可行性。 张伟等[7]利用CNN对GF-1 16 m多光谱数据进行了特征提取,并以SVM作为分类器对北京密云水库周边地区开展土地覆被分类,得到了97.99%的总体分类精度和0.98的Kappa系数,表明CNN可以提取更精细、更准确的地表覆盖特征。李亚飞等[4]基于CNN对Landsat-8 30 m的OLI影像开展土地覆被分类研究,取得了97.83% 的总体分类精度和0.97的Kappa系数,在影像分辨率相对较低的情况下仍获得了较高的分类精度。目前,CNN在土地覆被分类方面的研究更多倾向于CNN算法精度本身,而欠缺考虑影像空间分辨率对于CNN分类精度的影响。崔刚等[8]基于深度学习对不同空间分辨率影像的冬小麦识别进行了研究,结果表明随着影像分辨率的提高,冬小麦的识别精度也越高。但该研究为单一地物的提取,影像空间分辨率提高导致地物光谱复杂程度增加对其精度的影响可能不大,这与多分类的情况存在较大差别。

CNN在发展过程中也衍生出了一些新的深度学习方法,在影像分类方面较为出色的如全卷积神经网络(Fully Convolutional Network,FCN)。Maggiori 等[9]基于空间分辨率为1 m的多光谱数据研究表明:相比于CNN,利用FCN对建筑物进行提取会大幅度节省时间,同时提取精度也会有一定上升。杨瑞等[10]基于FCN改进和扩展而来的U-Net对高分二号影像融合后进行了多要素分类,取得了优于支持向量机和面向对象方法的分类精度。同传统机器学习算法相比,CNN不需要复杂的特征工程,拥有更好的泛化能力,分类精度更高;相比于FCN,CNN的计算量较大。但为获得更为普适性的规律,对其他衍生的深度学习方法也形成一定的参考价值,本文选用较为常规的CNN作为研究方法。

就目前文献检索结果来看,CNN对影像进行分类时,影像的空间分辨率多集中在4 m以上,这与CNN对数据上下文结构信息的要求有很大关系。但通常来讲,影像的空间分辨率越高,数据获取的费用及对其进行处理的时间成本也愈高。为明确CNN对不同分辨率影像分类的差异性,选取5种不同空间分辨率的遥感影像数据,分辨率跨度为4~30 m,探索不同空间分辨率对CNN分类精度的影响,为CNN方法下土地覆被分类影像的选择提供合理参考。

2 研究区及数据

2.1 研究区概况

西宁市地处湟水流域中游河谷盆地,面积为497 km2,地理坐标为36°12′27″~37°30′09″ N,100°47′56″~101°56′49″ E,市区平均海拔2 295 m,属于高原大陆性半干旱气候[11],湟水河自西向东流经市区。研究区为西宁市东部的一区域(图1),土地利用类型主要包括建设用地、林地、草地、耕地及水域,且分布无明显规律,能够充分验证CNN的学习能力。

图1

图1   研究区位置图

Fig.1   Location of the study area


2.2 数据及预处理

研究采用的影像数据为GF-2 4 m多光谱数据(免费来源于高分辨率对地观测青海数据与应用中心)、SPOT-6 6 m多光谱数据、Sentinel-2A 10 m多光谱数据(免费来自于欧洲航天局https://scihub.copernicus.eu/dhus/#/home)、Landsat-8 OLI的30 m空间分辨率数据(免费来源于中国地理空间数据云http://www.gscloud.cn/sources/)及其全色波段融合后的15 m空间分辨率数据。其中,GF-2和SPOT-6的多光谱数据含有红、绿、蓝和近红外4个波段,Sentinel-2A的多光谱数据共包含13个波段,红、绿、蓝和近红外4个波段为10 m空间分辨率,研究中仅采用了Sentinel-2A的10 m分辨率数据。为了保持数据波段的一致性,Landsat-8两种空间分辨率的数据也仅采用了红、绿、蓝和近红外4个波段。

表1   卫星影像及其参数

Table 1  Satellite images and their parameters

卫星平台所属国家空间分辨率/m空间参考使用光谱波段影像获取时间
GF-2中国4WGS-1984R、G、B、NIR2015-07-28
SPOT-6法国62016-08-08
Sentinel-2A欧空局102016-07-27
Landsat-8美国15/302016-08-07

注:Sentinel-2A影像的10 m空间分辨率波段为b2(蓝)、b3(绿)、b4(红)和b8(近红外);Landsat⁃8的15 m影像是经全色波段b8与其多光谱数据融合后得到,分类时使用了b2(蓝)、b3(绿)、b4(红)和b5(近红外)4个波段。

新窗口打开| 下载CSV


CNN作为一个深度学习模型,其最初的目的就是避免复杂的特征工程[12],因此CNN可以直接以图像的原始像素作为输入,不需要复杂的数据预处理工作。文中4种多光谱数据的地理参考均为WGS-1984椭球体,其中Sentinel-2A影像数据为Level-1C级数据,是经正射校正和亚像元级几何精校正后的正射影像产品[13],所以仅对SPOT-6和GF-2影像进行了几何校正,并对所有数据进行了大气校正处理。

3 CNN搭建与数据处理

3.1 CNN简介

CNN是多层前馈神经网络的一种变体,与传统的神经网络一样,它也是由可学习的权重和偏置组成,含有较多的隐含层是其区别于人工神经网络的关键[14]。相较于传统神经网络的全连接方式,CNN最大的特点在于局部连接和权值共享。科学家通过研究发现猫的一个视觉神经元只处理一小块区域内的视觉图像,这个区域被称为感受野[15]。具体的,局部连接和权值共享是指每个神经元不需要接收全部像元的信息,只需接受局部像元点的信息输入(感受野),然后将所有的局部信息综合起来便可得到整体信息。这一特点大大降低了训练的参数量及网络的复杂程度,所需训练的参数与影像无关,仅取决于卷积核的大小及数量,减小了训练的难度。同时,权值共享还赋予了CNN对图像平移的容忍性,提高了模型的泛化能力[16]。一个完整的CNN一般由卷积层、池化层、激活函数及全连接层构成。

卷积层由多个卷积核构成,卷积核通俗地理解为一系列的权重阵列[17]。卷积运算是指将卷积核的每个权重与相应的影像像元一一对应,分别求积后相加。一个卷积核负责提取图像的一种特征,为了从原始影像中提取到足够多的特征,一个卷积层往往含有多个卷积核。图2展示了一个卷积核对影像进行卷积运算得到一个特征图的过程,为了不改变卷积运算后图像的大小,在影像边缘部分采取补0的方式进行卷积运算。

图2

图2   卷积运算示意图

Fig.2   The schematic diagram of convolution operation


池化层又称下采样层,一般连接在卷积层之后,它的主要作用是对数据进行降维处理。常用的池化有平均池化和最大池化,本文采用最大池化算法,减小计算量的同时保留数据最显著的特征,增强模型的泛化能力[18]。采用的池化窗口大小为2×2像素,步长为2,进行池化操作前卷积层输出的数据大小为5×5像素,影像边缘部分采用补0的方式进行池化,经池化操作后数据大小为3×3像素(图3)。

图3

图3   最大池化示意图

Fig.3   Maximum pooling schematic diagram


由于卷积为线性操作,无论堆叠多少层,其最终结果都是经线性计算得到。但现实世界中很多问题是线性不可分的,为增强CNN的表达能力,往往会在卷积层之后连接一个激活函数,它是一个非线性函数,用来学习影像的非线性特征,本文采用常见的ReLu函数[19]作为激活函数,其函数图像见图4

图4

图4   ReLu函数图像

Fig.4   Graph of ReLu function


经过若干次的卷积和池化操作后,可以得到一些更高抽象层次的特征图,将所有的二维特征图展开拼接为一个一维向量,作为全连接层的输入。全连接层的每一个神经元都与前一层的所有像元相连接,一方面可以使用全连接层来扩充CNN的容量;另一方面,其主要作用是将学习到的抽象特征映射到样本标记空间,完成分类任务。全连接层常会产生过拟合现象,可以增加一个dropout层,在训练时随机抛弃一部分节点,减轻或者避免过拟合现象。

3.2 CNN的搭建

为实现CNN分类结果的像素级定位,参考FCN经反卷积层实现上采样的思想,研究中根据周围像元确定待分类像元的上下文结构信息,采用三次卷积插值的方法对影像进行上采样后,使用CNN进行分类,对比发现该方法取得的分类结果优于通常采用的patch-based[9]方法。

CNN的搭建基于Google旗下开源的深度学习框架Tensorflow[20]进行,文中搭建的CNN结构共9层,采用了5×5像素和3×3像素两种大小的卷积核,以便充分提取不同尺度的特征,每个卷积层后均跟有一个激活函数。由于输入图像较小,为了避免过多的信息损失,仅在第四层的卷积层和激活函数之后使用了一个大小为2×2像素,步长为2的最大池化层,并在第二个全连接层后使用了一个dropout层减轻过拟合现象,最后使用softmax分类器完成分类任务,网络详细结构如图5所示。

图5

图5   卷积神经网络结构示意图

Fig.5   The schematic diagram of CNN structure


3.3 分类体系及训练样本选取

根据全国遥感监测土地利用覆盖分类体系及城市用地分类与规划建设用地标准,结合研究区实际情况,土地利用共分为8类:林地(包括灌木林和疏林地)、耕地、草地、河流、水库坑塘、富营养化水体、城镇建设用地及工业仓储用地。其中富营养化水体仅在2016年的影像(Landsat-8、Sentinel-2A和SPOT-6)上存在,该水体受藻类影响,整体呈现植被的光谱特征,由于其斑块面积较小,在Landsat-8 30 m空间分辨率影像上未对其进行单独分类。对于影像中由高层建筑所产生的阴影,研究中简单地将其归为城镇建设用地。

训练样本是在参考高空间分辨率影像基础上,基于待分类影像进行选取,考虑到影像空间分辨率及获取时间不同,针对5种数据分别选取了不同的训练样本。CNN要求输入数据大小相同,例如文中选用的为5×5像素,目前常见的选样方式为逐块勾选并标记,受样本大小及数量的影响,该方法的工作量较大。由于CNN具有旋转的不变性,当样本数量不足时,通常采用的方法是进行样本增强,即将已有的样本旋转某个角度后作为新的样本,但这样难免会影响样本的丰富性和代表性。

为提高样本选择效率,采用了GDAL和C++相结合的辅助选样方法。首先对每类地物手动数字化选取一定数量较为纯净的像素块作为样本候选区,并尽量保证其在影像上均匀分布。然后使用一个5×5的窗口对影像进行遍历操作,窗口每前进一次,根据坐标判断该窗口是否落在数字化的某个多边形内,若是便将其裁剪并标记作为一个样本,否则继续遍历操作。使用固定大小的窗口对影像进行遍历,可以有效避免人工选择时的遗漏,既严格保证了样本大小,也能够较为快速地选出足够数量的样本。

由于研究区内不同的地物类型的面积差异显著,手动数字化的样本数量没有固定大小,主要根据地物的具体分布情况确定。朱秀芳等[21]对训练样本与分类方法分类精度响应关系的研究表明,对于同一种分类方法,当训练样本数量达到一定程度时,其分类精度是相对稳定的。考虑到本文的研究区域较小,为了使最终的训练样本数量满足统计要求,减少其对研究结果的影响,手动数字化过程中尽可能勾选较多的样本候选区。在这种情况下,对于空间分辨率较高的影像或斑块面积较大的地物,窗口滑动方法得到的训练样本数量较大,可以采取等间距取样的方法从已选出的样本中抽取一部分作为训练时实际使用的样本,或者在进行窗口滑动时设置合适的步长,以保证样本在空间上分布的相对均匀性。

将窗口滑动方法对每种地物选出的的最大样本数量设定为300个,其中80%的样本作为用于训练,剩余的20%用作交叉验证,对模型的精度进行评估。对于空间分辨率较低的影像或斑块面积较小的地物,选出的样本数量可能达不到300个,这种情况下将所有选出的样本均予采用,并按同样的比例划分训练样本和测试样本。

3.4 CNN分类参数设置

CNN的各项参数设置对训练及最终的分类结果有较大影响,前人的各项研究中参数设置均有较大差异[22],根据研究区的实际情况结合多次实验对比,研究中各项参数设置如下:batchsize设置为15,学习率设置为0.000 1,采用Adam优化器进行优化训练。为了防止训练过程中产生过拟合现象,dropout层的keep_prob在训练时设置为0.6,分类时设置为1。

受影像接收时刻与分辨率不同的影响,地物产生的阴影在GF-2影像中颜色较深,与正常水体的光谱较为类似,CNN无法直接将两者区分开,在CNN分类结束后计算归一化水体指数[23],阴影区域均为负值,选择0为阈值,将两者区分开,分类完成后对分类结果进行了碎斑去除处理。

4 结果和分析

4.1 分类结果

5种不同空间分辨率数据的CNN分类结果如图6所示,影像分辨率越高,分类结果对地物细节的展现越好。由于研究区内地物斑块整体偏小,导致影像分辨率为30 m时,分类结果中较大的斑块连结成片,边界粗糙,过小的斑块无法有效区分,仅少量的小斑块零星分布,加之不同地物间的光谱信息混合严重,出现了一些错分现象,分类结果较为杂乱;影像分辨率上升到15 m后,小的地物斑块增多,但地物分布细节仍未能得到有效展现;影像分辨率为10 m时,CNN分类结果对地物细节的反映出现了较大的提升,呈现较为明显的相间分布,同时地物边缘多为锯齿状;在SPOT-6影像中,大部分地物斑块在分类结果中均得到了很好的展现,边缘清晰平滑,不再有明显的锯齿状;影像空间分辨率进一步提高到4 m后,一方面地物细节表现更加丰富,另一方面分类结果也显得更为破碎。

图6

图6   标准假彩色影像与CNN分类结果

Fig.6   Standard false color images and CNN classification results


4.2 精度评价

分类结果精度评价主要通过总体分类精度、生产者精度、用户精度以及Kappa系数进行评价。验证样本是在野外调查和参考高空间分辨率影像的基础上,基于待分类影像手动数字化选取的感兴趣区,选取验证样本时尽可能覆盖到整个区域。考虑到影像的空间分辨率有较大差异,分别对5种数据选取了不同的验证样本,按空间分辨率从低到高,验证样本所包含的像元数量也依次递增,分别为828、2 636、4 886、10 490和23 287个。

5种不同空间分辨率下,CNN的总体分类精度均高于89.00%(表2),尤其是分辨率高于15 m以后,总体分类精度均高于90.00%,Kappa系数高于0.90,说明CNN对于不同空间分辨率的影像均有着较强的区分能力,能够获得较高的总体分类精度。

表2   分类精度评价

Table2  Classification accuracies evaluation

Landsat-8Lanssat-8(融合后)Sentinel-2ASPOT-6GF-2
生产者精度/%用户精度/%生产者精度/%用户精度/%生产者精度/%用户精度/%生产者精度/%用户精度/%生产者精度/%用户精度/%
林地88.4894.1296.3796.5493.1297.7393.5598.8493.7797.95
耕地89.2978.1393.4287.9894.8398.2193.6891.0596.6091.65
草地78.4387.5990.0887.5786.8278.9394.2480.6093.1775.79
河流76.5696.0884.1199.4589.9496.4097.07100.0090.7298.55
水库坑塘85.71100.0085.1995.8391.03100.0092.86100.0088.6999.62
富营养化水体————57.14100.0096.88100.0098.5391.78————
工业仓储用地87.1093.1094.92100.0099.6093.9899.0194.1796.8984.52
城镇建设用地98.8787.6394.2291.4793.5693.6191.4697.8390.9595.62
总体精度/%89.0192.6492.3193.6592.46
Kappa系数0.860.910.900.920.90

新窗口打开| 下载CSV


伴随着影像空间分辨率从低到高,不同地物的分类精度逐渐趋于稳定,生产者精度和用户精度方差的平均值整体上呈现先降低后升高的趋势,且SPOT-6的平均方差最低,图7以生产者精度为例,展示了这一趋势,不同地物分类精度中,极差最小的为SPOT-6影像。图7也反映出大部分地物随着影像分辨率的提高,生产者精度首先呈现较大幅度上升,但当影像分辨率达到一定程度时,这种上升趋势开始趋于稳定(如林地),甚至开始呈现出下降的趋势(如工业仓储用地及河流),这说明CNN对多光谱影像进行土地覆被分类存在一个最佳的空间分辨率。伴随着影像分辨率的提高,地物间的规律性也会逐渐减弱,出现较高的类内变异性,相对而言,不同地物类型间的差异性也就随之减小,可能会导致CNN的错分误分比例上升。

图7

图7   不同空间分辨率影像下地物生产者精度

Fig.7   Producer accuracies of features under different spatial resolution images


总体来看,SPOT-6的总体分类精度优于GF-2的4 m多光谱数据,且其不同地物精度的方差更小,分类结果更稳定,同时也较好地反映了地物细节,可以认为在所涉及的分辨率范围内,CNN对SPOT-6的6 m多光谱数据分类效果最佳。该结论与其他学者[24-26]的研究成果相吻合,明冬萍等[24-25]研究表明,土地覆被分类的最佳空间分辨率为3~7 m,杨闫君等[26]研究表明多光谱影像土地覆被分类的最佳分辨率为5 m,均与本文的结果较为接近。

图7中反映的信息来看,影像空间分辨率的提升在一定范围内对CNN提取面积较小或宽度较窄的地物提高分类精度有很大帮助,如富营养化水体本身的斑块面积较小,在Landsat 30 m的影像上未对其进行单独分类,与全色波段融合为15 m后仍未能取得较好的分类精度,但当影像分辨率升高到10 m后,其生产者精度出现大幅提高。受此类影响较为明显的还有同样面积较小的工业仓储用地及宽度较窄的河流。与其他地物不同,随着影像分辨率的提高,城镇建设用地的生产者精度一直呈现下降趋势,这是由于其本身的组成成分较为复杂,随着影像分辨率的提高,类内的规律性进一步减弱,增加了正确分类的难度。

4.3 结果分析

结合精度评价发现当影像空间分辨率为15 m时,CNN取得了92.64%的总体分类精度,甚至略高于Sentinel-2A的10 m多光谱数据,但就分类结果中对地物细节的表达来看,Landsat融合后的15 m数据与后者仍有较大差距。在15~10 m空间分辨率范围,CNN开始区分出丰富的地物细节,Sentinel-2A的10 m多光谱数据因为处在该区间,一些较小的地物虽然能够被CNN识别,但由于特征不足够明显,一些地物的分类精度出现了明显下降(如草地、建筑物之间的小块草地开始被识别且精度不高),同时部分地物受影像空间分辨率升高带来的正向影响,分类精度有所提高,总体来看Sentinel-2A影像的分类精度波动幅度较大(图7)。

综合5种不同空间分辨率影像的分类结果,CNN对光谱特征较为单一的水库坑塘分类效果较差,这与传统的机器学习分类方法有较大差别。研究区内水库坑塘的CNN分类结果生产者精度偏低,结合分类混淆矩阵发现水库坑塘比较容易被错分为林地和城镇建设用地。由于CNN主要依靠图像的上下文结构信息进行分类,水库坑塘这类质地较为均一的地物,邻域中像元的光谱差异很小,但靠近边缘部分的像元容易受到岸边其他地物(主要为林地和城镇建设用地)的光谱影响,同时CNN模型使用了最大池化算法和Relu函数,二者均偏向于学习像元值较大的像素特征,容易对光谱吸收特征明显的水体形成抑制作用[7],因此在水库坑塘的边缘部分会出现一定的错分现象。

CNN基于影像的上下文结构信息进行分类,一方面能够有效提高像元识别的准确率,另一方面对于边界部分的混合像元来讲,更容易受到其他地物的影响,造成CNN对地物边缘分类结果不准确。图8中将5种不同空间分辨率的影像进行分类细节展示,道路边缘部分的像元由于同时受到道路和林地的影响,被错分为了草地。同时由于顾及了影像上下文信息,CNN对较小的地物不敏感,在空间分辨率低于6 m的情况下,均未能提取出影像中的完整环路;当影像分辨率达到6 m时,CNN提取出了较为完整的环路轮廓;空间分辨率为4 m的GF-2影像上,CNN虽然完整提取出了该环路,但对于宽度为2~3个像元的道路,CNN仍未能完全识别。

图8

图8   不同空间分辨率影像CNN分类细节对比

Fig.8   Comparison of CNN classification details for different spatial resolution images


5 结 语

研究针对5种不同空间分辨率的影像数据,设计并搭建了CNN模型分类框架,采用了窗口滑动辅助选样方法,能够在仅更改少量参数的情况下,为不同的遥感影像选出足量且满足大小要求的训练样本。就5种数据的分类精度来看,该方法是可行的,为后期CNN在遥感影像分类中的样本选择提供了一种思路。

在不同空间分辨率下,CNN均取得了高于89.00%的总体分类精度和0.86以上的Kappa系数,分类精度较高。但CNN对地物的边缘部分混合像元的分类效果较差。同时对斑块较小的地物不敏感,难以对其进行正确分类,若要提取地物的精细结构,有赖于影像空间分辨率的进一步提高。

在所涉及的空间分辨率范围(4~30 m)内,考虑到总体分类精度、对地物细节信息的反映及不同地物分类精度分布的集中程度等多方面因素,本文结果表明CNN对SPOT-6 6 m多光谱数据进行土地覆被分类的效果最佳。当影像的空间分辨率进一步提高时,地物呈现的规律性逐渐降低,表现出较高类内变异性和低类间差异性,异质性对CNN分类精度的影响超过了影像分辨率提高带来的正向收益,增加了正确分类的困难程度,总体分类精度及Kappa系数反而出现了下降的趋势。

本文仍存在一些不足之处:一是针对采用的窗口滑动辅助选样方法没有量化其对分类结果的影响,仅通过最终的分类结果来间接验证了该方法的可行性;二是CNN分类结果对地物细节的反映程度未能使用一个具体的指标来衡量,仅通过直观感受对其进行了分析;三是在研究过程中未对地物的阴影进行单独处理,而是简单地将其归入了建设用地之中,可能由此产生一些错分误分现象,今后将考虑添加阴影去除的操作,以减小其对分类精度的影响。

参考文献

Tian Fuyou, Wu Bingfang, Zeng Hongwei, et al.

Identifying Soybean Cropped Area with Sentinel-2 Data and Multi-layer Neural Network

[J]. Journal of Geoinformation Science, 2019,21(6):918-927.

[本文引用: 1]

田富有吴炳方曾红伟.

基于多层神经网络与Sentinel-2数据的大豆种植区识别方法

[J].地球信息科学学报,2019,21(6):918-927.

[本文引用: 1]

Thanh H N T, Doan T M, Tomppo E, et al.

Land Use/Land Cover Mapping Using Multitemporal Sentinel-2 Imagery and Four Classification Methods—A Case Study from Dak Nong, Vietnam

[J]. Remote Sensing, 2020,12(9):1367. doi:10.3390/rs12091367.

Gu Xiaotian, Gao Xiaohong, Ma Huijuan, et al.

Comparision of Machine Learning Methods for Land Use/Land Cover Classification in the Complicated Terrain Regions

[J]. Remote Sensing Technology and Application, 2019,34(1):57-67.

[本文引用: 1]

谷晓天高小红马慧娟.

复杂地形区土地利用/土地覆被分类机器学习方法比较研究

[J].遥感技术与应用,2019,34(1):57-67.

[本文引用: 1]

Li Yafei, Dong Hongbin.

Classification of Remote-sensing Image based on Convolutional Neural Network

[J]. CAAI Transactions on Intelligent Systems,2018,13(4): 550556.

[本文引用: 2]

李亚飞董红斌.

基于卷积神经网络的遥感图像分类研究

[J]. 智能系统学报, 2018,13(4):550556.

[本文引用: 2]

Masoud M, Bahram S, Mohammad R, et al.

Very Deep Convolutional Neural Networks for Complex Land Cover Mapping Using Multispectral Remote Sensing Imagery

[J]. Remote Sensing, 2018,10(7):1119. doi:10.3390/rs10071119.

[本文引用: 2]

Meng Xiangrui, Zhang Shuqing,

Zang Shuying. Remote Sensing Classification of Wetland Communities based on Convolutional Neural Networks and High Resolution Images: A Case Study of the Honghe Wetland

[J]. Scientia Geographica Sinica, 2018,38(11):1914-1923.

[本文引用: 1]

孟祥锐张树清臧淑英.

基于卷积神经网络和高分辨率影像的湿地群落遥感分类——以洪河湿地为例

[J].地理科学,2018,38(11):1914-1923.

[本文引用: 1]

Zhang WeiZheng KeTang Pinget al

Land Cover Classification with Features Extracted by Deep Convolutional Neural Network

[J]. Journal of Image and Graphics201722(8):1144-1153.

[本文引用: 2]

张伟郑柯唐娉

深度卷积神经网络特征提取用于地表覆盖分类初探

[J]. 中国图象图形学报,201722(8) : 1144-1153.

[本文引用: 2]

Cui Gang, Wu Jinsheng, Yu Zhen, et al.

Comparison Analysis on Wheat Mapping Using Deep Learning Algorithm from Different Satellite Data Source

[J]. Remote Sensing Technology and Application2019,34(4):748- 755.

[本文引用: 1]

崔刚吴金胜于镇.

深度学习对不同分辨率影像冬小麦识别的适用性研究

[J].遥感技术与应用,2019,34(4):748-755.

[本文引用: 1]

Maggiori E, Tarabalka Y, Charpiat G, et al.

Fully Convolutional Neural Networks for Remote Sensing Image Classification

[C]// 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS). IEEE, 2016.

[本文引用: 2]

Yang RuiQi YuanSu Yang.

U-Net Neural Networks and Its Application in High Resolution Satellite Image Classification

[J]. Remote Sensing Technology and Application2019,34(6):1-8.

[本文引用: 1]

杨瑞祁元苏阳.

深度学习 U-Net方法及其在高分辨卫星影像分类中的应用

[J].遥感技术与应用,2019,34(6):1-8.

[本文引用: 1]

Liu Xuemei, Gao Xiaohong, Jia Wei.

Spatio Temporal Change of Urban Heat Island Effect in Xining from Landsat Image

[J]. Journal of Earth Environment, 2019, 10(2):166-176.

[本文引用: 1]

刘雪梅高小红贾伟.

基于遥感影像的近30年西宁市热岛效应时空变化

[J].地球环境学报,2019,10(2):166-176.

[本文引用: 1]

Wang Can, Wu Xinhui, Li Lianqing, et al.

Convolutional Neural Network in Prediction of Soil Moisture Content

[J]. Spectroscopy and Spectral Analysis, 2018,38(1):36-41.

[本文引用: 1]

王璨武新慧李恋卿.

卷积神经网络用于近红外光谱预测土壤含水率

[J].光谱学与光谱分析,2018, 38(1):36-41.

[本文引用: 1]

Tian Ying, Chen Zhuoqi, Hui Fengming, et al

. ESA Sentinel-2A/B Satellite: Characteristics and Applications

[J]. Journal of Beijing Normal University (Natural Science Edition), 2019,55(1):57-65.

[本文引用: 1]

田颖陈卓奇惠凤鸣.

欧空局哨兵卫星Sentinel-2A/B数据特征及应用前景分析

[J].北京师范大学学报(自然科学版),2019,55(1):57-65.

[本文引用: 1]

Dong Yunya, Zhang Qian.

A Survey of Depth Semantic Feature Extraction of High-resolution Remote Sensing Images based on CNN

[J]. Remote Sensing Technology and Application, 2019,34(1):1-11.

[本文引用: 1]

董蕴雅张倩.

基于CNN的高分遥感影像深度语义特征提取研究综述

[J].遥感技术与应用,2019,34(1):1-11.

[本文引用: 1]

Zhou Feiyan, Jin Linpeng, Dong Jun.

Review of Convolutional Neural Network

[J]. Chinese Journal of Computers, 2017,40(6):1229-1251.

[本文引用: 1]

周飞燕金林鹏董军.

卷积神经网络研究综述

[J].计算机学报, 2017,40(6):1229-1251.

[本文引用: 1]

Huang H, Sun G Y, Zhang X M, et al.

Combined Multiscale Segmentation Convolutional Neural Network for Rapid Damage Mapping from Postearthquake Very High-resolution Images

[J]. Journal of Applied Remote Sensing,2019,13(2): 022007.doi:10.1117/1.JRS.13.022007.

[本文引用: 1]

Alshehhi R, Marpu P R, Woon W L, et al.

Simultaneous Extraction of Roads and Buildings in Remote Sensing Imagery with Convolutional Neural Networks

[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2017, 130:139-149.

[本文引用: 1]

Scherer D, Müller A, Behnke S.

Evaluation of Pooling Operations in Convolutional Architectures for Object Recognition

[C]//International Conference on Artificial Neural Networks. 2010.

[本文引用: 1]

Xu L, Choy C S, Li Y W.

Deep Sparse Rectifier Neural Network for Speech Denoising

[C]//2016 IEEE International Workshop on Acoustic Signal Enhancement (IWAENC). IEEE, 2016.

[本文引用: 1]

Shu Na, Liu Bo, Lin Weiwei, et al.

Survey of Distribute Machine Learning Platforms and Algorithms

[J]. Computer Science, 2019,46(3):9-18.

[本文引用: 1]

舒娜刘波林伟伟.

分布式机器学习平台与算法综述

[J].计算机科学,2019,46(3):9-18.

[本文引用: 1]

Zhu Xiufang, Pan Yaozhong, Zhang Jinshui, et al

. The Effects of Training Samples on the Wheat Planting Area Measure Accuracy in TM Scale(Ⅰ):The Accuracy Response of Different Classifiers to Training Samples

[J]. Journal of Remote Sensing,2007,11(6):826-837.

[本文引用: 1]

朱秀芳潘耀忠张锦水.

训练样本对TM尺度小麦种植面积测量精度影响研究(Ⅰ)——训练样本与分类方法间分类精度响应关系研究

[J].遥感学报,2007,11(6):826-837.

[本文引用: 1]

Nogueira K, Penatti Otávio A B, Santos J A D.

Towards Better Exploiting Convolutional Neural Networks for Remote Sensing Scene Classification

[J]. Pattern Recognition, 2017,61:539-556. doi:10.1016/j.patcog.2016.07.001.

[本文引用: 1]

McFEETERS S.

K. The Use of the Normalized Difference Water Index (NDWI) in the Delineation of Open Water Features

[J]. International Journal of Remote Sensing, 1996,17(7):1425-1432.

[本文引用: 1]

Ming Dongping, Wang Qun, Yang Jianyu.

Spatial Scale of Remote Sensing Image and Selection of Optimal Spatial Resolution

[J]. Journal of Remote Sensing, 2008,12(4): 529-537.

[本文引用: 2]

明冬萍王群杨建宇.

遥感影像空间尺度特性与最佳空间分辨率选择

[J].遥感学报,2008,12(4):529-537.

[本文引用: 2]

Chen J, Deng M, Xiao P F, et al.

Optimal Spatial Scale Choosing for High Resolution Imagery based on Texture Features Frequency Analysis

[J]. Journal of Remote Sensing, 2011,15(3):492-511.

[本文引用: 1]

Yang Yanjun, Tian Qingjiu, Zhan Yulin, et al.

Effects of Spatial Resolution and Texture Features on Multi-spectral Remote Sensing Classification

[J]. Journal of Geo-information Science, 2018,20(1):99-107.

[本文引用: 2]

杨闫君田庆久占玉林.

空间分辨率与纹理特征对多光谱遥感分类的影响

[J].地球信息科学学报,2018,20(1):99-107.

[本文引用: 2]

/