Integrating Sentinel Active and Passive Remote Sensing Data to Land Cover Classification in a National Park from GEE Platform
-
摘要: 国家公园的土地覆盖分类对于掌握自然资源现状、查明存在的生态安全威胁并快速应对具有基础性数据支撑作用。基于谷歌地球引擎(Google Earth Engine,GEE)平台,结合哨兵(Sentinel)主被动遥感数据及其导出的光谱指数、纹理特征和地形特征,分别采用基于像元的随机森林(random forest,RF)算法和面向对象的简单非迭代聚类(simple non-iterative clustering,SNIC)+RF算法实现了钱江源国家公园异质性景观的土地覆盖(耕地、森林、草地、水体、人造地表和裸地)分类。地面实验表明,在多种输入数据组合中,基于像元和面向对象方法分类获得的最高总体精度分别为92.37%和93.98%。合成孔径雷达(synthetic aperture radar,SAR)数据的纳入能够提高基于像元方法的分类精度,但在面向对象方法中未能体现精度提升效果。通过SNIC+RF算法生成的土地覆盖分类图完整性更好,所需特征数量较少,并且算法能够在GEE环境下快速执行,适用于国家公园管理实践。Abstract:Objectives Land cover classification in national parks plays an important role in understanding the status of natural resources, identifying the existing ecological security threats and responding to them quickly.Methods Two land cover classification methods are developed based on Google Earth Engine (GEE) platform by combining Sentinel active and passive remote sensing data, and spectral indices, textural features and topographic features derived from the data to classify land cover types in the Qianjiangyuan National Park (cropland, forest, grassland, water body, artificial surface and bare land). One used pixel-based random forest (RF) classification algorithm, the other used object-oriented simple non-iterative clustering (SNIC) segmentation in partnership with RF algorithm.Results The ground experimental results show that the highest overall classification accuracies of the pixel-based method and the object-oriented method are 92.37% and 93.98%, respectively. Furthermore, the integration of synthetic aperture radar (SAR) data can substantially improve the classification accuracy when using the pixel-based method, but there is no apparent escalating effect for the object-oriented method.Conclusions Land cover classification map generated by SNIC+RF algorithm in GEE platform is more complete and the algorithm requires fewer features and runs quickly in GEE platform. Thus, this algorithm deserves to be popularized in national park management practices.
-
随着国家海洋战略的实施,水下潜器的导航需求越来越大,对此需要进行高精度、高自主的导航条件作为辅助。现阶段,提供导航定位的方式多样化,如GPS(global positioning system)、GLONASS(global navigation satellite system)、GSNS(Galileo satellite navigation system)、BDS(BeiDou navigation satellite system)、声呐导航、天文导航[1-2]、惯性/测向测距导航系统[3],其定位精度虽然满足水下潜器需求,但由于隐蔽性差、抗干扰能力差等原因,不利于潜艇长时间水下定位。重力辅助惯性导航作为一种不依靠人为因素,借助地球自身特性,即地球重力场信息,完成潜艇水下长航、长航距导航,具有重要的研究意义。重力辅助惯性导航具有自主、无源、高精度、高隐蔽性、抗干扰性的特点,是真正的无源导航。
利用地球独有特性进行无源导航的方式,目前已涉及地形匹配导航、地磁匹配导航[4-6]以及重力匹配导航。由于惯性导航在极区易丢失方向,地磁导航更难实施,因此亟需研究极区多源信息组合导航定位技术[7]。20世纪80年代初,美国洛克希德马丁公司在美国军方资助下研制了重力敏感器系统(gravity sensor system,GSS)。GSS是一个当地水平的稳定平台,平台上安装有一个重力仪和三个重力梯度仪。重力仪是一个垂直安装的高精度加速度计,重力梯度仪输出两组正交的梯度分量,由安装在同一转轮上的4个加速度计组成。GSS用于实时估计垂线偏差,以补偿惯性导航误差。20世纪90年代,该公司在GSS、静电陀螺导航仪(electrostatic gyro navigator,ESGN)、重力基准图和深度探测仪等技术的基础上开发无源重力辅助导航系统(gravity aided inertial navigation system,GAINS),它通过重力匹配获取导航位置坐标,以无源方式限定或修订惯性导航误差,能实现潜艇14 d精确导航[8-9]。国内对重力辅助惯性导航技术的研究开始于20世纪90年代。目前关于重力辅助导航技术及其相关技术的研究包括重力辅助导航系统组成,重力匹配导航算法,卫星测高反演重力异常,重力场对惯性导航系统精度的影响,重力图数据处理,以及重力仪、重力梯度仪的研制等方面[10]。
目前常用的匹配算法有4种,包括地形轮廓匹配算法(terrain contour matching,TERCOM)、迭代最近等值线算法(iterated closest contour point,ICCP)、桑迪亚惯性地形辅助导航算法(Sandia inertia terrain-aided navigation,SITAN)和滤波算法。TERCOM算法主要是使用重力场数据相关技术的匹配算法,其匹配精度高,受初始误差影响小,但实时性不强[11-15]。ICCP算法主要是使用刚性变换不断靠近等值线上最近点并进行配准,该算法有较高的精度,但其应用受到惯导指示位置与载体真实位置之间误差足够小的假设条件的限制,而且随着迭代次数的增加,计算量也不断增大[16]。此外,ICCP算法忽略了参考导航系统在短时间内的误差变化,应用时限制了匹配点数[17]。SITAN算法是基于扩展卡尔曼滤波技术的单点匹配算法,得到导航误差的估计值,然后把误差估值反馈回惯导系统,从而对其导航状态进行修正[18-19];SITAN算法需要获取较为精确的初始误差,对非线性观测模型线性化敏感,线性化精度低会导致滤波发散[20]。在辅助导航中使用较多的滤波算法主要有卡尔曼滤波算法以及粒子滤波算法,卡尔曼滤波算法是将重力场数据的变化作为卡尔曼滤波器的观测量,将惯导运动的误差方程作为卡尔曼滤波器的状态量来建立卡尔曼滤波方程。
发展海洋重力匹配导航技术还存在很多核心问题,如海洋重力基准图的格网分辨率问题,实时重力仪测量值与重力地图快速匹配问题[21]。近年来,重力实时测量与匹配导航算法引起国内学者广泛关注[22-24],如Wang等[25]研制了重力匹配惯性导航系统,并开展了相关的海试工作。不少学者在传统匹配算法的基础上,改进算法的实时性,改善匹配的精度,提高匹配的可实施性[26-29]。
传统重力匹配过程主要考虑实时重力测量值与重力场背景图之间的相关性。国内学者在进行约束改进算法时,提出了对惯导短时间内具有高精度的约束条件,进而增设轨迹位移角度判断方式,通过增加对每一段采样航距的惯导角度和待匹配轨迹对应角度一致的条件,对待匹配轨迹进行了约束,但惯导系统在短时间内的误差虽然很小,但不可忽略。对惯导设备要求苛刻,仅仅单独从每一段采样航距去进行约束,以及对每一次轨迹旋转角度进行绝对约束,不符合实际情况[30]。
本文通过结合运载体的运动速度与航向信息,提出了基于惯性导航特征相似性的约束条件匹配导航算法,用向量的形式表示轨迹的运动轨迹方向,同时计算轨迹在每一采样间隔段的航距以及总匹配航距。利用匹配约束条件对待匹配轨迹进行选择,剔除大量不符合约束条件的随机匹配轨迹,同时考虑惯导设备存在的客观误差因素,提高匹配效率,缩短匹配时间,改善匹配精度。
1 重力匹配算法原理
重力辅助匹配惯性导航系统是利用地球重力场特性,在高分辨率的重力异常图组成的基本特性信息数据库中[31],结合惯性导航系统提供的实时输出位置,以及海洋重力仪提供的实测值,按照一定的匹配算法对惯导位置进行修正,实现高精度水下定位的一种导航技术。系统原理结构图[32]如图 1所示。
概率神经网络匹配算法主要由4个结构层构成,以此为输入层、模式层、求和层以及输出层。
在输入层中,输入样本X经过归一化后传递到模式层中。
在模式层中,需要将输入样本X与权值向量Wx做内积运算:
$$ {Z_X} = {X^{\rm{T}}}{\mathit{\boldsymbol{W}}_X}$$ (1) ZX然后对做非线性化处理,即:
$$ g\left( {{Z_X}} \right) = {\rm{exp}}\left( {\left( {{Z_X} - 1} \right)/{\sigma ^2}} \right)$$ (2) 将式(1)、式(2)经过推导等价于:
$$ g\left( x \right) = {\rm{exp}}[ - \left( {X - {W_X}{)^{\rm{T}}}\left( {X - {W_X}} \right)/2{\sigma ^2}} \right]$$ (3) 式中,σ表示概率神经网络的平滑因子。
在求和层中,将属于同类的非线性函数进行累加,以求得输入样本中属于每一类的概率ηa,一般利用Parzen窗方法来求取概率ηa,即:
$$ {\eta _a} = \frac{1}{{\left( {2{\rm{ \mathsf{ π} }}} \right)\frac{p}{2}{\sigma _p}}}\left( {\frac{1}{{{n_a}}}} \right)\mathop \sum \limits_{x = 1}^{{n_a}} g\left( {{Z_X}} \right)$$ (4) 式中,na表示求和层中属于分类a的神经元的个数;p表示输入样本的维数。
在输出层中,将求和层中利用Parzen窗方法计算而来的所有分类的概率进行比较,并求取属于某个类别的最大概率ηmax,然后以竞争方式使其输出为1。
将所有待匹配轨迹N条分为N个类别,求取实测重力序列与之对应的最大概率匹配序列作为最终匹配轨迹序列[33]。
概率神经网络通过分析实测重力数据与待匹配轨迹数据之间的属性关系,对比实测数据与各类别属性的从属关系程度,选择最大相关序列、作为最佳匹配序列。概率神经网络算法虽能较好实现重力匹配惯性导航,但其在匹配过程中仅依靠了惯导轨迹提供的单点位置信息搜索待匹配点,忽略了惯性导航所提供的连续轨迹属性信息,即航距与航向,这会导致概率神经网络匹配结果只考虑了重力数据的相似性以及顺序性,而忽略了潜艇轨迹本该具有的空间属性。
2 相关约束的重力匹配算法
相关约束的重力匹配算法的核心思想是:首先利用惯性导航系统提供的连续性方向信息构建序列矩阵,对每一条待匹配轨迹的拐点方向矩阵进行判断,排除不符合惯导轨迹方向序列的待匹配轨迹;其次,利用惯导系统提供的相邻采样距离信息和轨迹总运行距离信息,对待匹配轨迹进行约束,排除局部和整体均不符合距离约束条件的轨迹;最后利用概率神经网络算法对符合以上两步约束条件的待匹配轨迹进行分析,选择最终匹配轨迹。该算法充分利用了惯性导航输出数据的短时高精度特性和在长距离匹配背景下的方向、距离整体特征连续性,并且避免了在重力异常图分辨率较低背景下对方向角度出现错误筛选的情况,提高了匹配算法的适用性,体现了重力匹配惯性导航算法中先整体、后局部,惯导为主、重力为辅的匹配思想。
2.1 设定搜索区域
由于惯性导航系统中加速度计以及陀螺仪的固有误差,使得惯导在输出位置时出现漂移,水下潜艇的实际位置应该在以惯导输出位置为中心的一定范围内,通过对惯导数据以及时长的分析来确定对应一定时间的惯导漂移范围。
本文以10个连续惯导输出点作为一组匹配轨迹,根据惯导输出确定其在重力图中的位置,分别以其为中心,以最后一点对应的累计漂移距离λ为1/2矩阵搜索边长,确定最后两点矩阵的搜索区域,即2λ×2λ,第3至8点搜索边长为,第1至2点的搜索边长为λ/2,以此确定一组矩阵的搜索区域[34],如图 2所示。
2.2 寻找等值点
在搜索范围内,根据重力底图,读取范围内对应格网位置的重力异常值gk,并与相应时刻重力仪实测重力值gi进行比较。根据重力仪误差情况,设定等值点范围,若$\left| {{g_i}\left( {{x_i}, {y_j}} \right) - {g_k}\left( {{x_k}, {y_k}} \right)} \right| = {\rm{\Delta }}g \le \delta $,则点k属于对应时刻惯导输出点i的重力等值点。最终在整条轨迹线中,可形成10个点聚集区。
2.3 设定约束条件
重力匹配惯性导航是一项多源综合性导航,目前大多数匹配算法在研究过程中侧重于研究待匹配重力值与实测重力值之间的相关性,忽略了惯性导航提供轨迹数据的精确性,以及对匹配导航位置关系的约束性。约束匹配流程图如图 3所示。
2.3.1 约束轨迹方向
虽然惯性导航系统中的加速度计和陀螺仪在导航中均存在漂移误差,但在机动性强、运动轨迹方向多变情况下的导航,不存在轨迹前进左右方向判别上的错误。因此,在陀螺仪不能完全精确输出轨迹转动角度的情况下,可以更为准确地输出轨迹相对运行前进方向的左右判别,这对于正确判断真实轨迹的运动起到很好的约束作用[35]。
因此,本文对待匹配轨迹的约束条件之一便是:对于一条由10个采样间隔组成的运动轨迹,判别从第3个采样点至第9个采样点处于其对应前两点构成的直线的左右方位,在判断左右方位时,首先需要判断轨迹的航向,即判断潜艇运动方向是向东还是向西,向南还是向北,其中,向东包括东南方向、东北方向以及正东方向,向西包括西北方向、西南方向以及正西方向,向北即为正北方向,向南即为正南方向。式(5)为直线表达公式,由式(6)可知,利用直线上两点坐标,可求得直线斜率。
$$ Y - y = k\left( {X - x} \right) $$ (5) $$ k = \frac{{{y_2} - {y_1}}}{{{x_2} - {x_1}}} $$ (6) 式(5)和式(6)经过推导等价于:
$$ Y = \frac{{{y_2} - {y_1}}}{{{x_2} - {x_1}}}\left( {X - {x_1}} \right) + {y_1} $$ (7) 求取第3点X轴向对应的Y轴坐标:
$$ {{y'}_3} = \frac{{{y_2} - {y_1}}}{{{x_2} - {x_1}}}\left( {{x_3} - {x_1}} \right) + {y_1} $$ (8) 判断其方向及输出:
$$ {{y'}_3} - {y_3} > 0\;\;\;\;\; - 1 $$ (9) $$ {{y'}_3} - {y_3} < 0\;\;\;\;\; + 1 $$ (10) $$ {{y'}_3} - {y_3} = 0\;\;\;\;\;\;0 $$ (11) 对于向东运动,若第3个采样点位于由第1、第2个采样点构成直线的下方,则判断其方向为-1;若位于直线上方,则判断其方向为+1;若位于直线内,则判断其方向为0。最终对惯导轨迹可输出一个由-1、0、+1组成的8×1的矩阵。
若潜艇向西运动,则输出结果与向东运动结果相反:
$$ {{y'}_3} - {y_3} > 0\;\;\;\;\; + 1 $$ (12) $$ {{y'}_3} - {y_3} < 0\;\;\;\;\; - 1 $$ (13) $$ {{y'}_3} - {y_3} = 0\;\;\;\;\;\;0 $$ (14) 若潜艇向正北方向运动,则输出结果与X值相关:
$$ {x_1} = {x_2} > {x_3}\;\;\;\; - 1 $$ (15) $$ {x_1} = {x_2} < {x_3}\;\;\;\; + 1 $$ (16) $$ {x_1} = {x_2} = {x_3}\;\;\;\;0 $$ (17) 若潜艇向正南方向运动,则输出结果与正北方向相反:
$$ {x_1} = {x_2} > {x_3}\;\;\;\; + 1 $$ (18) $$ {x_1} = {x_2} < {x_3}\;\;\;\; - 1 $$ (19) $$ {x_1} = {x_2} = {x_3}\;\;\;\;0 $$ (20) 同时,对由待匹配点任意连线构成的待匹配轨迹进行方向判断,生成一个8×N的矩阵,其中N为待匹配轨迹条数。根据惯导提供的方向矩阵,从待匹配方向矩阵中,筛选出符合惯导矩阵的对应轨迹,作为下一步待匹配轨迹,剔除不符合轨迹。
2.3.2 约束轨迹距离
通过对惯导系统在长航距中产生的误差分析,可以发现,在较短的时间间隔内,惯导输出的距离精度较高,在相邻采样间隔的输出距离上,与对应的真实运动轨迹距离相仿,因此,本文对每条待匹配轨迹的9个采样间隔段进行判断,若
$$ \begin{array}{*{20}{l}} {\;\;\left| {L\left( {A\left( {{x_t}, {y_t}} \right), B\left( {{x_{t + 1}}, {y_{t + 1}}} \right)} \right) - } \right.}\\ {\left. {L\left( {A'\left( {{x_t}, {y_t}} \right), B'\left( {{x_{t + 1}}, {y_{t + 1}}} \right)} \right)} \right| \le \lambda {\rm{}}} \end{array} $$ (21) 则认为两条线段的距离相仿[36-37]。式中,λ根据惯导元器件精度、重力图分辨率以及采样时长设定。
此外,增设约束条件:即对于由10个采样点构成的航迹,总的距离误差不是简单地由9个段误差累计和,而是有增有减,因此,对惯导轨迹的总距离和待匹配轨迹的总距离之差设定一个阈值ξ,ξ的值可根据惯导元器件的精度和采样时间间隔来确定。即:
$$ \begin{array}{l} \left| {\mathop \sum \limits_1^n \left( {\sqrt {{{({x_{n + 1}} - {x_n})}^2} + {{({y_{n + 1}} - {y_n})}^2}} - } \right.} \right.\\ \left. {\mathop \sum \limits_1^n \left. {\sqrt {{{(x_{n + 1}^{\rm{'}} - x_n^{\rm{'}})}^2} + {{(y_{n + 1}^{\rm{'}} - y_n^{\rm{'}})}^2}} } \right)} \right| \le \xi \end{array} $$ (22) 式中,根据惯导元器件精度、重力图分辨率以及轨迹时长设定。这样,不仅仅是从单个采样距离上对待匹配轨迹做出约束,而是从局部和整体两方面约束,可有效剔除符合局部条件、但不符合整体实际轨迹特征的待匹配轨迹。
3 仿真实验分析
本文基于MATLAB对实验进行编程仿真,重力异常图分辨率为,假定水下潜艇以10 m/s的速度航行,匹配算法选取概率神经网络匹配,每3 min进行一次采样,在满足匹配导航时效性的前提下,为能够尽量满足匹配导航精度条件,实验采用每10个采样点作为一次匹配序列,在同等条件下进行有无约束条件对比,实验中α设定为2 776 m,即1.5倍的格网距离,ξ设定为4 627 m,即为2.5倍的格网距离。轨迹所处背景场如图 4所示,重力异常图背景标准差为16.942 mGal,经度方向相关系数为0.769,纬度方向的相关系数为0.547。因约束条件是对概率神经网络待匹配结果进行筛选,并不在最终原有匹配结果轨迹的基础上加以修改,因此该约束算法不受背景场影响,图 4中右下点为起始点。仿真结果如图 5所示。
由图 5可见,在经纬度方向重力异常相关系数均大于0.5的情况下,未增加约束条件的概率神经网络匹配,在匹配位置上会存在一些不符合实际运动方向、航距的情况。基于约束条件的重力匹配辅助导航算法,对轨迹的航向做出修正,真实轨迹与惯导轨迹均在第8条采样段的航向上向东北方向偏移,但未约束的概率神经网络匹配结果将第9段匹配航迹修改成向西南方向航行,明显不符合实际情况,通过约束,可将其进行正确修正。同时,算法通过剔除大量不符轨迹,大大缩短了神经网络匹配的耗时,如表 1所示,将精度由千米级提高到百米级,较好地提高了匹配位置精度。以匹配误差低于一个格网密度作为成功匹配,将成功匹配点占所有匹配点的百分比作为成功匹配率,在东、北方向,匹配率均由80%提升至100%。如表 2所示,方向正确率由62.5%提升至100%,耗时缩短50%左右,表明本文算法有效提高了运算速度,保证了匹配运算的精确性和实时性。
表 1 东、北方向距离误差/mTable 1. Distance Errors of East, North Directions/m方向 匹配条件 误差统计 最大值 最小值 平均值 均方差 北向 未约束神经网络匹配 5 852.439 0 5.502 187 1 059.555 0 1 801.611 0 约束神经网络匹配 1 620.846 4 5.502 187 591.739 2 618.789 4 东向 未约束神经网络匹配 5 156.672 0 8.035 500 1 221.924 0 1 727.069 0 约束神经网络匹配 1 442.420 2 8.035 500 361.244 1 413.000 4 表 2 航行方向结果比对Table 2. Results Comparison of Navigation Directions航向 转向 方向正确率/% 耗时/s 1 2 3 4 5 6 7 8 惯导方向 1 -1 -1 1 -1 1 1 -1 未约束轨迹方向 1 -1 -1 1 -1 -1 -1 1 62.5 0.094 约束轨迹方向 1 -1 -1 1 -1 1 1 -1 100 0.040 4 结语
基于相关约束的重力匹配算法,根据惯导元器件精度、重力图分辨率、重力仪精度的具体情况,对搜索区域进行规划,对等值点进行选取,对待匹配轨迹的方向向量进行筛选,以及相邻采样点距离和总轨迹距离进行约束,有效排除大量不符合真实轨迹特点的干扰轨迹,提高了导航效率和导航精度,增加了匹配算法的鲁棒性,避免了无可匹配点的情况。仿真结果表明,该约束算法不仅对均方差匹配和神经网络匹配的位置结果进行修正,同时对神经网络匹配的运算速度有大幅提升,提高了该匹配算法的实时性,更好地满足水下潜艇的需求。基于相关约束的重力匹配算法为进一步提高水下潜器导航的精确性和实时性提供了一条新的技术途径。
-
表 1 分类特征组合
Table 1 Combinations of Classification Features
实验序号 特征类型 变量总数 1-1 Sentinel-2光谱波段 10 2-1 Sentinel-2光谱波段+光谱指数 19 2-2 Sentinel-2光谱波段+Sentinel-2纹理特征 20 2-3 Sentinel-2光谱波段+地形特征 13 3-1 Sentinel-2光谱波段+光谱指数+Sentinel-2纹理特征 29 3-2 Sentinel-2光谱波段+光谱指数+地形特征 22 3-3 Sentinel-2光谱波段+Sentinel-2纹理特征+地形特征 23 4-1 Sentinel-2光谱波段+光谱指数+Sentinel-2纹理特征+地形特征 32 5-1 Sentinel-2光谱波段+光谱指数+地形特征+ Sentinel-1极化波段+ Sentinel-1纹理特征 32 5-2 Sentinel-2光谱波段+光谱指数+Sentinel-2纹理特征+地形特征+ Sentinel-1极化波段 34 5-3 Sentinel-2光谱波段+光谱指数+Sentinel-2纹理特征+地形特征+Sentinel-1纹理特征 40 6-1 Sentinel-2光谱波段+光谱指数+Sentinel-2纹理特征+地形特征+Sentinel-1极化波段+Sentinel-1纹理特征 42 表 2 不同特征组合的分类精度及排序
Table 2 Classification Accuracy and Ordering of Different Feature Combinations
排序 基于像元方法 面向对象方法 实验序号 总体精度/% Kappa系数 实验序号 总体精度/% Kappa系数 1 5-1 92.37 0.90 2-1 92.37 0.90 2 5-2 91.97 0.90 5-1 92.37 0.90 3 5-3 91.97 0.90 4-1 91.97 0.90 4 6-1 91.97 0.90 3-1 91.57 0.89 5 2-3 91.57 0.89 3-2 91.57 0.89 6 3-1 91.57 0.89 5-3 91.16 0.89 7 3-2 91.57 0.89 5-2 90.76 0.88 8 2-1 91.16 0.89 6-1 90.76 0.88 9 4-1 91.16 0.89 2-3 89.56 0.87 10 3-3 90.36 0.88 1-1 88.76 0.86 11 1-1 89.56 0.87 3-3 88.76 0.86 12 2-2 88.35 0.85 2-2 87.15 0.84 表 3 不同分类实验的特征选择结果
Table 3 Feature Selection Results of Different Classification Experiments
筛选次数 基于像元方法(实验5-1) 面向对象方法(实验2-1) 面向对象方法(实验5-1) 变量数 总体精度/% 变量数 总体精度/% 变量数 总体精度/% 1 32 92.37 19 92.37 32 92.37 2 26 92.37 15 93.98 26 92.77 3 21 91.57 12 91.57 21 90.76 4 17 91.97 10 90.36 17 91.57 5 14 91.16 8 88.76 14 91.97 6 11 90.36 — — 11 91.97 7 9 89.56 — — 9 89.96 表 4 基于像元分类图的混淆矩阵和精度评估结果
Table 4 Confusion Matrix and Classification Accuracy Measurement of the Map Classified by Pixel-Based Method
分类 耕地 森林 草地 水体 人造地表 裸地 总和 用户精度/% 错分误差/% 耕地 31 0 0 0 1 4 36 86.11 13.89 森林 0 88 1 0 0 0 89 98.88 1.12 草地 3 2 46 0 0 0 51 90.20 9.80 水体 0 0 0 20 0 0 20 100.00 0.00 人造地表 1 0 0 3 26 1 31 83.87 16.13 裸地 1 0 1 1 0 19 22 86.36 13.64 总和 36 90 48 24 27 24 249 生产者精度/% 86.11 97.78 95.83 83.33 96.30 79.17 漏分误差/% 13.89 2.22 4.17 16.67 3.70 20.83 总体精度/% 92.37 Kappa 0.90 表 5 面向对象分类图的混淆矩阵和精度评估结果
Table 5 Confusion Matrix and Classification Accuracy Measurement of the Map Classified by Object-Oriented Method
分类 耕地 森林 草地 水体 人造地表 裸地 总和 用户精度/% 错分误差/% 耕地 32 0 0 0 1 2 35 91.43 8.57 森林 0 89 1 0 0 0 90 98.89 1.11 草地 2 1 47 0 0 0 50 94.00 6.00 水体 0 0 0 21 0 0 21 100.00 0.00 人造地表 2 0 0 2 25 2 31 80.65 19.35 裸地 0 0 0 1 1 20 22 90.91 9.09 总和 36 90 48 24 27 24 249 生产者精度/% 88.89 98.89 97.92 87.50 92.59 83.33 漏分误差/% 11.11 1.11 2.08 12.50 7.41 16.67 总体精度/% 93.98 Kappa 0.92 -
[1] He S, Su Y, Wang L, et al. Taking an Ecosystem Services Approach for a New National Park System in China[J]. Resources, Conservation and Recycling, 2018, 137: 136-144. doi: 10.1016/j.resconrec.2018.04.030
[2] Wang Y Q, Lu Z, Sheng Y W, et al. Remote Sensing Applications in Monitoring of Protected Areas[J]. Remote Sensing, 2020, 12(9): 1370. doi: 10.3390/rs12091370
[3] Chen J, Chen J, Liao A, et al. Global Land Cover Mapping at 30 m Resolution: A POK-Based Operational Approach[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2015, 103: 7-27. doi: 10.1016/j.isprsjprs.2014.09.002
[4] Gong P, Liu H, Zhang M N, et al. Stable Classification with Limited Sample: Transferring a 30 m Resolution Sample Set Collected in 2015 to Mapping 10 m Resolution Global Land Cover in 2017[J]. Science Bulletin, 2019, 64(6): 370-373. doi: 10.1016/j.scib.2019.03.002
[5] Fawzi N I, Husna V N, Helms J A. Measuring Deforestation Using Remote Sensing and Its Implication for Conservation in Gunung Palung National Park, West Kalimantan, Indonesia[J]. IOP Conference Series: Earth and Environmental Science, 2018, 149: 012038. doi: 10.1088/1755-1315/149/1/012038
[6] Recanatesi F, Giuliani C, Ripa M. Monitoring Mediterranean Oak Decline in a Peri-Urban Protected Area Using the NDVI and Sentinel-2 Images: The Case Study of Castelporziano State Natural Reserve[J]. Sustainability, 2018, 10(9): 3308. doi: 10.3390/su10093308
[7] Lucas R, Kerchove R V, Otero V, et al. Structural Characterisation of Mangrove Forests Achieved Through Combining Multiple Sources of Remote Sensing Data[J]. Remote Sensing of Environment, 2020, 237: 111543. doi: 10.1016/j.rse.2019.111543
[8] Clerici N, Valbuena C C A, Posada J M. Fusion of Sentinel-1A and Sentinel-2A Data for Land Cover Mapping: A Case Study in the Lower Magdalena Region, Colombia[J]. Journal of Maps, 2017, 13(2): 718-726. doi: 10.1080/17445647.2017.1372316
[9] Tavares P A, Ely S B N, Guimarães U S, et al. Integration of Sentinel-1 and Sentinel-2 for Classification and LULC Mapping in the Urban Area of Belém, Eastern Brazilian Amazon[J]. Sensors, 2019, 19(5): 1140. doi: 10.3390/s19051140
[10] Talukdar S, Singha P, Mahato S, et al. Land-Use Land-Cover Classification by Machine Learning Classifiers for Satellite Observations: A Review[J]. Remote Sensing, 2020, 12(7): 1135. doi: 10.3390/rs12071135
[11] Li C, Gong P, Wang J, et al. The First All-Season Sample Set for Mapping Global Land Cover with Landsat-8 Data[J]. Science Bulletin, 2017, 62(7): 508-515. doi: 10.1016/j.scib.2017.03.011
[12] Breiman L. Random Forests[J]. Machine Learning, 2001, 45(1): 5-32. doi: 10.1023/A:1010933404324
[13] Wendelberger K S, Gann D, Richards J H. Using Bi-Seasonal WorldView-2 Multi-spectral Data and Supervised Random Forest Classification to Map Coastal Plant Communities in Everglades National Park[J]. Sensors, 2018, 18(3): 829. doi: 10.3390/s18030829
[14] Tsai Y, Stow D, Chen H, et al. Mapping Vegetation and Land Use Types in Fanjingshan National Nature Reserve Using Google Earth Engine[J]. Remote Sensing, 2018, 10(6): 927. doi: 10.3390/rs10060927
[15] 郎芹, 牛振国, 洪孝琪, 等. 青藏高原湿地遥感监测与变化分析[J]. 武汉大学学报(信息科学版), 2021, 46(2): 230-237. doi: 10.13203/j.whugis20180277 Lang Qin, Niu Zhenguo, Hong Xiaoqi, et al. Remote Sensing Monitoring and Change Analysis of Wetlands in Qinghai-Tibet Plateau[J]. Geomatics and Information Science of Wuhan University, 2021, 46(2): 230-237. doi: 10.13203/j.whugis20180277
[16] Lü X W, Ming D P, Chen Y Y, et al. Very High Resolution Remote Sensing Image Classification with SEEDS-CNN and Scale Effect Analysis for Superpixel CNN Classification[J]. International Journal of Remote Sensing, 2019, 40(2): 506-531. doi: 10.1080/01431161.2018.1513666
[17] 袁旭. 基于超像素的图像分割方法研究[D]. 武汉: 华中科技大学, 2019. Yuan Xu. Research on Image Segmentation Method Based on Super-Pixel[D]. Wuhan: Huazhong University of Science and Technology, 2019.
[18] Achanta R, Süsstrunk S. Superpixels and Polygons Using Simple Non-iterative Clustering[C]//IEEE Conference on Computer Vision and Pattern Recognition, Honolulu, USA, 2017.
[19] Tu Y, Chen B, Zhang T, et al. Regional Mapping of Essential Urban Land Use Categories in China: A Segmentation-Based Approach[J]. Remote Sensing, 2020, 12(7): 1058. doi: 10.3390/rs12071058
[20] Herrero H, Southworth J, Bunting E. Utilizing Multiple Lines of Evidence to Determine Landscape Degradation Within Protected Area Landscapes: A Case Study of Chobe National Park, Botswana from 1982 to 2011[J]. Remote Sensing, 2016, 8(8): 623. doi: 10.3390/rs8080623
[21] 余建平, 申云逸, 宋小友, 等. 钱江源国家公园体制试点区功能分区对黑麂保护的有效性评估[J]. 生物多样性, 2019, 27(1): 5-12. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDY201901002.htm Yu Jianping, Shen Yunyi, Song Xiaoyou, et al. Evaluating the Effectiveness of Functional Zones for Black Muntjac(Muntiacus Crinifrons) Protection in Qianjiangyuan National Park Pilot Site[J]. Biodiversity Science, 2019, 27(1): 5-12. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDY201901002.htm
[22] 米湘成. 生物多样性监测与研究是国家公园保护的基础[J]. 生物多样性, 2019, 27(1): 1-4. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDY201901001.htm Mi Xiangcheng. Biodiversity Monitoring and Research are Basis of National Park Conservation[J]. Biodiversity Science, 2019, 27(1): 1-4. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDY201901001.htm
[23] Farr T G, Rosen P A, Caro E, et al. The Shuttle Radar Topography Mission[J]. Reviews of Geophysics, 2007, 45(2): RG2004.
[24] Haklay M, Weber P. OpenStreetMap: User-Generated Street Maps[J]. IEEE Pervasive Computing, 2008, 7(4): 12-18. doi: 10.1109/MPRV.2008.80
[25] Immitzer M, Vuolo F, Atzberger C. First Experience with Sentinel-2 Data for Crop and Tree Species Classifications in Central Europe[J]. Remote Sensing, 2016, 8(3): 166. doi: 10.3390/rs8030166
[26] Kaplan G. Sentinel-2 Pan Sharpening: Comparative Analysis[C]//The 2nd International Electronic Conference on Remote Sensing, 2018.
[27] 杨晓涵, 李雨彤, 王宁, 等. 遥感影像融合方法比较研究[J]. 重庆工商大学学报(自然科学版), 2017, 34(2): 90-94. doi: 10.16055/j.issn.1672-058X.2017.0002.018 Yang Xiaohan, Li Yutong, Wang Ning, et al. Comparative Research on Remote Sensing Image Fusion Method[J]. Journal of Chongqing Technology and Business University (Natural Science Edition), 2017, 34(2): 90-94. doi: 10.16055/j.issn.1672-058X.2017.0002.018
[28] Pettorelli N, Vik J O, Mysterud A, et al. Using the Satellite-Derived NDVI to Assess Ecological Responses to Environmental Change[J]. Trends in Ecology & Evolution, 2005, 20(9): 503-510.
[29] Huete A, Didan K, Miura T, et al. Overview of the Radiometric and Biophysical Performance of the MODIS Vegetation Indices[J]. Remote Sensing of Environment, 2002, 83(1/2): 195-213.
[30] Gitelson A A, Merzlyak M N. Remote Sensing of Chlorophyll Concentration in Higher Plant Leaves[J]. Advances in Space Research, 1998, 22(5): 689-692. doi: 10.1016/S0273-1177(97)01133-2
[31] Louhaichi M, Borman M M, Johnson D E. Spatially Located Platform and Aerial Photography for Documentation of Grazing Impacts on Wheat[J]. Geocarto International, 2001, 16(1): 65-70. doi: 10.1080/10106040108542184
[32] Gitelson A A, Stark R, Grits U, et al. Vegetation and Soil Lines in Visible Spectral Space: A Concept and Technique for Remote Estimation of Vegetation Fraction [J]. International Journal of Remote Sensing, 2002, 23(13): 2537-2562. doi: 10.1080/01431160110107806
[33] Gitelson A A, Kaufman Y J, Merzlyak M N. Use of a Green Channel in Remote Sensing of Global Vegetation from EOS-MODIS[J]. Remote Sensing of Environment, 1996, 58(3): 289-298. doi: 10.1016/S0034-4257(96)00072-7
[34] Merzlyak M N, Gitelson A A, Chivkunova O B, et al. Non-destructive Optical Detection of Pigment Changes During Leaf Senescence and Fruit Ripening[J]. Physiologia Plantarum, 1999, 106(1): 135-141. doi: 10.1034/j.1399-3054.1999.106119.x
[35] Qi J, Chehbouni A, Huete A R, et al. A Modified Soil Adjusted Vegetation Index[J]. Remote Sensing of Environment, 1994, 48(2): 119-126. doi: 10.1016/0034-4257(94)90134-1
[36] 徐涵秋. 利用改进的归一化差异水体指数(MNDWI)提取水体信息的研究[J]. 遥感学报, 2005, 9(5): 589-595. https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB200505011.htm Xu Hanqiu. A Study on Information Extraction of Water Body with the Modified Normalized Difference Water Index (MNDWI)[J]. Journal of Remote Sensing, 2005, 9(5): 589-595. https://www.cnki.com.cn/Article/CJFDTOTAL-YGXB200505011.htm
[37] 李鹏, 黎达辉, 李振洪, 等. 黄河三角洲地区GF-3雷达数据与Sentinel-2多光谱数据湿地协同分类研究[J]. 武汉大学学报(信息科学版), 2019, 44(11): 1641-1649. doi: 10.13203/j.whugis20180258 Li Peng, Li Dahui, Li Zhenhong, et al. Wetland Classification Through Integration of GF-3 SAR and Sentinel-2B Multispectral Data over the Yellow River Delta[J]. Geomatics and Information Science of Wuhan University, 2019, 44(11): 1641-1649. doi: 10.13203/j.whugis20180258
[38] Díaz-Uriarte R, Alvarez de Andrés S. Gene Selection and Classification of Microarray Data Using Random Forest[J]. BMC Bioinformatics, 2006, 7: 3.
[39] Nembrini S, König I R, Wright M N. The Revival of the Gini Importance?[J]. Bioinformatics, 2018, 34(21): 3711-3718.
[40] Erinjery J J, Singh M, Kent R. Mapping and Assessment of Vegetation Types in the Tropical Rainforests of the Western Ghats Using Multispectral Sentinel-2 and SAR Sentinel-1 Satellite Imagery[J]. Remote Sensing of Environment, 2018, 216: 345-354. https://www.sciencedirect.com/science/article/pii/S003442571830333X
[41] Persson M, Lindberg E, Reese H. Tree Species Classification with Multi-temporal Sentinel-2 Data[J]. Remote Sensing, 2018, 10(11): 1794.
[42] Hamad R, Balzter H, Kolo K. Multi-criteria Assessment of Land Cover Dynamic Changes in Halgurd Sakran National Park (HSNP), Kurdistan Region of Iraq, Using Remote Sensing and GIS[J]. Land, 2017, 6(1): 18.
[43] Wang M M, He G J, Ishwaran N, et al. Monitoring Vegetation Dynamics in East Rennell Island World Heritage Site Using Multi-sensor and Multi-temporal Remote Sensing Data[J]. International Journal of Digital Earth, 2020, 13(3): 393-409.
[44] 林栋, 秦志远, 童晓冲, 等. 融合光谱及形态学信息的对象级空间特征提取方法[J]. 武汉大学学报(信息科学版), 2018, 43(5): 704-710. doi: 10.13203/j.whugis20150627 Lin Dong, Qin Zhiyuan, Tong Xiaochong, et al. Objected-Based Structural Feature Extraction Method Using Spectral and Morphological Information[J]. Geomatics and Information Science of Wuhan University, 2018, 43(5): 704-710. doi: 10.13203/j.whugis20150627
[45] 王艳琦, 秦福莹, 银山, 等. 基于GF-1卫星数据与面向对象分类的达里诺尔湿地自然保护区土地覆盖信息提取[J]. 中国农学通报, 2019, 35(10): 137-141. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNTB201910025.htm Wang Yanqi, Qin Fuying, Yin Shan, et al. Land Cover Information Extraction of Dalinuoer Wetland Nature Reserve: Based on GF-1 Satellite Data and Object-Oriented Classification[J]. Chinese Agricultural Science Bulletin, 2019, 35(10): 137-141. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNTB201910025.htm
-
期刊类型引用(1)
1. 陈鸿鑫,马天霆,周阳,简彦辰,高犇,戴明露. 基于CNN-GAN数据增强网络的电厂锅炉管道温度压力及健康状态预测. 电子器件. 2023(06): 1593-1600 . 百度学术
其他类型引用(5)