基于无人机影像的喀斯特农耕区地物识别——以桂林市为例
1.
2.
3.
Object recognition of karst farming area based on UAV image: A case study of Guilin
1.
2.
3.
通讯作者: 付波霖(1988-),男,博士,副教授,主要从事RS和GIS的应用研究。Email:fbl2012@126.com。
责任编辑: 李 瑜
收稿日期: 2019-02-1 修回日期: 2019-03-15 网络出版日期: 2020-03-15
基金资助: |
|
Received: 2019-02-1 Revised: 2019-03-15 Online: 2020-03-15
作者简介 About authors
娄佩卿(1995-),男,硕士,主要从事遥感图像智能处理研究。Email:gislou@126.com。 。
为了探究低空无人机遥感技术对喀斯特地貌条件下不同形态农耕区地物类型的识别精度,以桂林市3个200 m×200 m样方的农耕区为研究区,在无人机航拍影像和地面调查数据的支持下,分别将基于像元和面向对象的影像分析技术与支持向量机(support vector machine,SVM)算法相结合,构建不同地貌条件下农耕区地物遥感识别模型,并进行精度对比分析。结果表明,面向对象的SVM分类结果保留了原始地物的大致轮廓,且地块较完整,更为适用于喀斯特地貌条件下的农耕区地物识别,较基于像元的SVM分类方法总体精度高6.54%,Kappa系数高0.135; 基于像元的SVM分类方法适用于地物分布规则的农耕区地物识别,相比面向对象的SVM分类方法总体精度高2.92%,Kappa系数高0.026。
关键词:
In order to explore the recognition accuracy of remote sensing technology of low-altitude UAV for surface features in agricultural areas with different forms under karst landform conditions, the authors chose three agricultural areas (each having an area size of 200 m×200 m) in Guilin City as the research object. Supported by UAV aerial images and ground survey data, the image analysis technology based on pixel and object-oriented was combined with support vector machine (SVM) algorithm, respectively, to build the remote sensing recognition model of agricultural areas under different geomorphological conditions, and the precision was comparatively studied and analyzed. The results show that the object-oriented SVM classification results retain the rough outline of the original ground features, and the plot is relatively complete, and hence this means is more suitable for the recognition of ground features in agricultural areas under karst landform conditions. Compared with the pixel based SVM classification method, the overall accuracy is higher by 6.54% , and the Kappa coefficient is higher by 0.135 . The SVM classification method based on pixel is suitable for feature recognition in agricultural areas with regular feature distribution. Compared with the object-oriented SVM classification method, the overall accuracy is higher by 2.92% and the Kappa coefficient is higher by 0.026 .
Keywords:
本文引用格式
娄佩卿, 陈晓雨, 王疏桐, 付波霖, 黄永怡, 唐廷元, 凌铭.
LOU Peiqing, CHEN Xiaoyu, WANG Shutong, FU Bolin, HUANG Yongyi, TANG Tingyuan, LING Ming.
0 引言
1995年,Vapnik等[12,13]提出支持向量机(support vector machine,SVM)这一新的机器学习算法。该算法具有分类结构简单、训练难度低、分类精度高、收敛速度快等优点[14,15,16,17,18,19]。目前SVM算法主要与基于像元和面向对象的影像分析技术相结合。在基于像元的SVM算法中,Maulik等[20,21]采用基于迭代过程的半监督SVM对城镇SPOT遥感影像进行基于像元的地物分类,结果精度高于标准SVM算法; Yang等[22]提出了一种基于像元的最小二乘SVM的RGB图像分割方法,结合像元及纹理特征,提高了SVM分类器的泛化能力; Mugiraneza等[23]采用基于像元的SVM算法定量评价城市景观变化动态空间格局和可持续城市生态系统的管理。在面向对象的SVM算法中,Li等[24]提出了一种结合栅格分析和矢量分析的面向对象SVM土地覆盖分类方法,对QuickBird影像进行地物分类的精度高于基于像元的分类方法; Malik等[25]首次利用多尺度分割算法结合SVM分类器对Alsat-2A高分影像进行地物分类,总体精度达到94.23%,Kappa系数达到0.93,获得了很好的分类效果; 裴欢等[26]采用基于纹理的面向对象SVM分类方法对河北省石家庄市Landsat8 OLI遥感影像进行地物分类,总体分类精度高达85.67%,为区域土地利用/覆盖信息提取提供了有效的途径。
综上所述,基于SVM算法的遥感影像地物识别研究已有很大进展,但所用遥感影像多以中等尺度空间分辨率的星载光学影像为主,较少使用到低空无人机遥感影像,且多集中应用于城镇区域的地物识别,对于喀斯特地貌条件下典型农耕区的地物高精度识别,国内学者研究较少。而且,由于受制于喀斯特地貌类型,目前农耕区地物识别仍以人工实地勾勒和人工目视解译无人机影像为主,效率低下且费时费力,亟需要一套快速、自动和高精度的影像识别算法。为此,本文以桂林市农耕区域为研究区,选取具有代表性的3种不同地形农耕区的无人机影像,分别将基于像元和面向对象的影像分析技术与SVM机器学习算法进行组合,对农耕区地物类型进行识别,并结合地面调查数据进行精度验证,进而对比了2种影像处理技术在复杂地形条件下的农耕区地物类型识别方面的精度差异,以期为地理国情普查中农耕区地物类型的遥感自动识别提供科学依据及技术参考。
1 研究区概况及数据预处理
1.1 研究区概况及数据源
选取桂林市3个200 m×200 m样方的农耕区为研究区,并分为梯田+林地复合区、梯田区和耕地+居民地复合区3种农耕区类型(表1)。2017年5—6月利用“大疆精灵4 PRO”无人机对3个研究区域进行航摄成像。影像航拍选择光照条件良好的晴天,且无风、无云,拍摄时间为10: 00—15: 00(北京时间),太阳高度角大于45°,阳光近似直射的条件下进行。无人机飞行高度为80~100 m,航向重叠度70%,旁向重叠度80%,光学镜头垂直向下每3 s对地自动对焦拍摄一次,以保证生成研究区空间分辨率为cm级的影像数据(图1)。同时,在研究区布设若干地面控制点,并采用GPS-RTK获取控制点的经纬度坐标用于影像的正射纠正、拼接镶嵌等处理。
表1 不同农耕区及其典型地物类型
Tab.1
农耕区类型 | 典型地物类型 |
---|---|
梯田+林地复合区 | 水体、道路、建筑物、林地、秧田和水田 |
梯田区 | 草地、经济作物、道路、建筑物、林地和果园 |
耕地+居民地复合区 | 经济作物、林地、菜地、果园、草地、水田、秧田、水体、道路和建筑物 |
图1
图1
航摄耕地+居民地复合区的部分原始影像
Fig.1
Some original images of cultivated land and residential land compound area
1.2 数据预处理
分别将3个农耕区航摄的原始影像和地面控制点导入Smart 3D中,进行影像匹配连接与编辑,剔除重叠率不够的影像,联合GPS-RTK获取控制点进行空中三角测量和光束法区域平差,构建地面立体模型并进行绝对定向,生成具备标准地理参考的数字表面模型(digital surface model,DSM),进一步提取研究区的高精度数字正射影像(digital orthophoto maps,DOM)。设置投影坐标系为WGS 84/UTM zone 49N。最后,利用ArcGIS软件对生成的DOM进行镶嵌,获得试验区完整的DSM和DOM,具体参数见表2和图2。由于3个研究区地貌类型差异较大,无人机航摄高度存在差异,故不同地区无人机航射影像空间分辨率也不同。
表2 无人机航摄影像部分参数
Tab.2
农耕区类型 | 图幅大小/像元 | 空间分辨率/m | 飞行高度/m |
---|---|---|---|
梯田+林地复合区 | 5 222×5 063 | 0.069 | 100 |
梯田区 | 8 252×5 801 | 0.042 | 85 |
耕地+居民地复合区 | 9 111×7 454 | 0.044 | 85 |
图2
为保证选取的样本具有代表性,需利用ArcGIS软件中随机点生成的方法在3个研究区域内随机均匀生成150个样本点,采用国家统计局广西调查总队地面调查矢量数据来确定样本点对应的地物类型。将样本点50%作为训练集,50%作为验证集。基于像元和面向对象的分类样本及验证样本坐标一一对应,以此进行精度对比。具体技术路线如图3所示。
图3
2 研究方法
2.1 理论基础
式中: ai为拉格朗日乘子(不为0的支持向量); N为样本数量; b为分类阈值;
式中: C为惩罚系数; w为权重向量; b为分类阈值,均在拉格朗日函数取极小值时求得; K(xi,x)为核函数,其他参数同上。选择不同核函数构成的SVM有所不同,本文选用径向基核函数(radial basis function,RBF)[33],即
2.2 模型构建
2.2.1 基于像元的SVM模型构建
利用基于像元的样本训练集结合无人机DOM进行模型训练,引入RBF对影像进行分类,由于分类影像波段数为3,所以内核函数值为0.333,概率阈值设为0,降低像元未被分类的概率。地物类型面积较大且光谱相近的无人机影像(耕地+居民地复合区)在ENVI软件中再基于灰度共生矩阵提取纹理特征,并将其与DOM进行波段融合,形成多源遥感影像集,再进行算法模型训练。
2.2.2 面向对象的SVM模型构建
表3 影像分割参数
Tab.3
影像名称 | 分割尺度 | 形状因子 | 紧实度 | 波段权重 (R,G,B) |
---|---|---|---|---|
梯田+林地复合区 | 30 | 0.4 | 0.5 | 1,1,1 |
梯田区 | 200 | 0.3 | 0.5 | 1,1,1 |
耕地+居民地复合区 | 300 | 0.3 | 0.5 | 1,1,1 |
图4
图4
梯田+林地复合区及梯田区影像多尺度分割结果
Fig.4
Multi-scale segmentation results of images of terraced +forest area and terraced area
将提取的训练样本结合无人机DOM数据集进行SVM算法分类器的模型训练,采用线性核函数,惩罚系数设为2,分割后结合点类型分类样本训练集对地物进行分类。分割因子中,分割尺度越大,地物分块越大,则影像空间分辨率越低。3幅分类影像中,梯田+林地复合区和耕地+居民地复合区2幅影像分割尺度及地块大小差异最大。
2.3 分类后处理及精度评价
基于像元的SVM分类结果出现了孤立、细小图斑,降低了分类结果的视觉效果,利用ENVI软件中的Majority/Minority,Clump和Sieve等卷积滤波算法将小斑块进行剔除或合并,邻域窗口设置为5×5。
为了进一步对比验证2种影像分析技术在农耕区典型地物识别中的精度差异,利用验证样本对分类结果进行精度验证,分别选取整体精度、Kappa系数、各类别用户精度和制图精度等指标进行定量评价,并利用混淆矩阵进一步探究精度差异的原因。
3 结果与分析
3.1 基于像元和面向对象的影像分类对比分析
基于像元的和面向对象的SVM分类结果对比如图5—7所示。
图5
图5
梯田+林地区分类结果对比
Fig.5
Comparison between the classification results of terrace + forest area and the original image
图6
图6
梯田区分类结果对比
Fig.6
Comparison between the classification results of terrace area and the original image
图7
图7
耕地+居民地复合区分类结果对比
Fig.7
Comparison between the classification results of cultivated land and residential compound area and the original image
3.2 整体分类精度对比分析
基于像元的SVM算法分类结果中,梯田+林地复合区影像基于像元分类的整体精度为91.55%,Kappa系数为0.83; 梯田区影像基于像元的整体分类精度为81.32%,Kappa系数为0.72; 耕地+居民地复合区影像基于像元整体分类精度为69.93%,Kappa系数为0.65。而面向对象的SVM算法分类结果中梯田+林地复合区影像面向对象的分类精度为90.70%,Kappa系数为0.89; 梯田区影像面向对象的分类精度为87.86%,Kappa系数为0.85; 耕地+居民地复合区影像面向对象分类精度为67.01%,Kappa系数为0.62。2种分类方法中均是梯田+林地区影像的总体精度(高于90%)及Kappa系数(高于0.8)最高,耕地+居民地复合区的总体精度(低于70%)及Kappa系数(低于0.7)最低。基于像元及面向对象分类结果中梯田+林地复合区的整体分类精度最高,耕地+居民地复合区的整体分类精度最低。
3.3 混淆矩阵
梯田+林地复合区影像共包含6种地物类别,采用基于像元的SVM算法进行分类时,水体的分类精度最高,无像元被误分; 水田的分类精度最低,像元分别误分到道路(317个)、建筑物(26个)类别,误分率为24.77%。采用面向对象的SVM算法进行分类时建筑物的分类精度最高,仅有1个对象误分到道路类别,误分率为3.57%; 秧田的分类精度最低,对象分别误分到水田(5个)、林地(1个)类别,误分率为16.67%(表4)。
表4 梯田+林地复合区混淆矩阵
Tab.4
分类方法 | 类别 | 水田 | 秧田 | 林地 | 道路 | 建筑物 | 水体 | 总量 | 误分率/% |
---|---|---|---|---|---|---|---|---|---|
基于像元 | 水田 | 1 042 | 0 | 0 | 317 | 26 | 0 | 1 385 | 24.77 |
秧田 | 0 | 930 | 75 | 3 | 0 | 0 | 1 008 | 7.74 | |
林地 | 40 | 439 | 10 032 | 82 | 4 | 77 | 10 674 | 6.01 | |
道路 | 108 | 0 | 0 | 763 | 29 | 0 | 900 | 15.22 | |
建筑物 | 0 | 0 | 0 | 74 | 993 | 3 | 1 070 | 7.20 | |
水体 | 0 | 0 | 0 | 0 | 0 | 73 | 73 | 0 | |
面向对象 | 水田 | 41 | 0 | 0 | 2 | 0 | 0 | 43 | 4.65 |
秧田 | 5 | 30 | 1 | 0 | 0 | 0 | 36 | 16.67 | |
林地 | 0 | 2 | 21 | 0 | 0 | 2 | 25 | 16.00 | |
道路 | 0 | 0 | 0 | 26 | 1 | 0 | 27 | 3.70 | |
建筑物 | 0 | 0 | 0 | 1 | 27 | 0 | 28 | 3.57 | |
水体 | 0 | 0 | 2 | 0 | 0 | 11 | 13 | 15.38 |
梯田区影像共包含6种地物类别,采用基于像元的SVM算法进行分类时,草地的分类精度最高,无像元被误分; 果园的分类精度最低,像元分别误分到经济作物(718个)、草地(343个)、道路(243个)、建筑物(18个)、林地(3个)类别,误分率为30.40%。采用面向对象的SVM算法进行分类时建筑物及草地的分类精度最高,误分率均为0; 而林地的分类精度最低,对象分别误分到果园(7个)、草地(3个)、经济作物(1个)类别,误分率为36.67%(表5)。
表5 梯田区混淆矩阵
Tab.5
分类方法 | 类别 | 果园 | 林地 | 建筑物 | 道路 | 经济作物 | 草地 | 总量 | 误分率/% | |
---|---|---|---|---|---|---|---|---|---|---|
基于像元 | 果园 | 3 034 | 3 | 18 | 243 | 718 | 343 | 4 359 | 30.40 | |
林地 | 1 218 | 8 307 | 3 | 15 | 10 | 427 | 9 980 | 16.76 | ||
建筑物 | 0 | 0 | 1 024 | 56 | 107 | 0 | 1 187 | 13.73 | ||
道路 | 8 | 0 | 3 | 1 045 | 22 | 0 | 1 078 | 3.06 | ||
经济作物 | 98 | 0 | 5 | 3 | 878 | 13 | 997 | 11.94 | ||
草地 | 0 | 0 | 0 | 0 | 0 | 134 | 134 | 0.00 | ||
面向对象 | 果园 | 26 | 2 | 0 | 0 | 1 | 0 | 29 | 10.34 | |
林地 | 7 | 19 | 0 | 0 | 1 | 3 | 30 | 36.67 | ||
建筑物 | 0 | 0 | 26 | 0 | 0 | 0 | 26 | 0 | ||
道路 | 0 | 0 | 0 | 24 | 2 | 0 | 26 | 7.69 | ||
经济作物 | 1 | 0 | 0 | 0 | 15 | 0 | 16 | 6.25 | ||
草地 | 0 | 0 | 0 | 0 | 0 | 13 | 13 | 0 |
耕地+居民地复合区影像共包含10种地物类别,采用基于像元的SVM算法进行分类时,水体的分类精度最高,无像元被误分; 秧田的分类精度最低,像元分别误分到果园(170个)、草地(154个)、菜地(121个)、水体(21个)、建筑物(6个)、水田(1个)类别,误分率为44.44%。采用面向对象的SVM算法进行分类时道路的分类精度最高,仅有2个对象误分到经济作物类别,误分率为6.25%; 草地的分类精度最低,对象分别误分到菜地(7个)、果园(5个)、秧田类别(2个),误分率为46.67%(表6)。
表6 耕地+居民地复合区混淆矩阵
Tab.6
分类方法 | 类别 | 建筑物 | 道路 | 水体 | 秧田 | 水田 | 草地 | 果园 | 菜地 | 林地 | 经济作物 | 总量 | 误分率/% |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
基于像元 | 建筑物 | 518 | 127 | 23 | 0 | 47 | 0 | 0 | 1 | 0 | 1 | 652 | 27.75 |
道路 | 23 | 497 | 0 | 0 | 34 | 0 | 0 | 0 | 0 | 213 | 689 | 35.20 | |
水体 | 0 | 0 | 136 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 217 | 0 | |
秧田 | 6 | 0 | 21 | 670 | 1 | 154 | 170 | 121 | 0 | 63 | 823 | 44.44 | |
水田 | 83 | 62 | 4 | 88 | 604 | 0 | 1 | 0 | 0 | 17 | 686 | 29.69 | |
草地 | 0 | 0 | 2 | 9 | 0 | 76 | 11 | 0 | 0 | 0 | 674 | 22.45 | |
果园 | 1 | 1 | 21 | 56 | 0 | 177 | 492 | 35 | 26 | 0 | 762 | 39.18 | |
菜地 | 5 | 0 | 0 | 0 | 0 | 3 | 0 | 98 | 0 | 0 | 301 | 7.55 | |
林地 | 15 | 2 | 10 | 0 | 0 | 264 | 88 | 46 | 1535 | 0 | 1 561 | 21.68 | |
经济作物 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 102 | 396 | 0.97 | |
面向对象 | 建筑物 | 51 | 1 | 8 | 5 | 10 | 0 | 0 | 0 | 0 | 2 | 652 | 33.77 |
道路 | 0 | 30 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 2 | 39 | 6.25 | |
水体 | 4 | 0 | 20 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 28 | 20.00 | |
秧田 | 0 | 0 | 0 | 13 | 0 | 0 | 6 | 1 | 0 | 1 | 24 | 38.10 | |
水田 | 1 | 4 | 0 | 0 | 16 | 0 | 0 | 0 | 0 | 4 | 36 | 36.00 | |
草地 | 0 | 0 | 0 | 2 | 0 | 16 | 5 | 7 | 0 | 0 | 19 | 46.67 | |
果园 | 0 | 0 | 0 | 3 | 0 | 0 | 17 | 0 | 0 | 0 | 30 | 15.00 | |
菜地 | 0 | 0 | 0 | 0 | 0 | 3 | 2 | 8 | 0 | 1 | 18 | 42.86 | |
林地 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 21 | 0 | 21 | 8.70 | |
经济作物 | 6 | 4 | 0 | 1 | 10 | 0 | 0 | 0 | 0 | 3 | 13 | 87.50 |
4 结论
1)面向对象的SVM分类结果保留了原始地物的大致轮廓,且地块较完整,而基于像元的SVM分类结果出现大量小斑块,地块完整度较低。
2)在地物分布较规则的农耕区(耕地+居民地复合区)基于像元的SVM分类方法整体分类效果较好,相比面向对象的SVM分类方法总体精度高2.92%,Kappa系数高0.026。在地物分布较复杂的农耕区(梯田区)面向对象的SVM分类方法整体分类效果较好,相比基于像元的SVM分类总体精度高6.54%,Kappa系数高0.135。
3)喀斯特地貌条件下的农耕区内主要地物类型为水田、秧田、果园、菜地、经济作物及水体。在地物分布较为复杂的农耕区(梯田区、梯田+林地复合区)面向对象的SVM分类方法对水田(用户精度高20.12%)、果园(用户精度高20.06%)、经济作物(用户精度高5.69%)的分类效果较好; 基于像元的SVM分类方法对秧田(用户精度高8.93%)、水体(用户精度高15.38%)的分类效果较好。在地物分布较为规则的农耕区(耕地+居民地复合区)面向对象的SVM分类方法对秧田(用户精度高6.34%)、果园(用户精度高24.18%)的分类效果较好; 基于像元的SVM分类方法对水田(用户精度高6.31%)、菜地(用户精度高35.31%)、经济作物(用户精度高86.53%)、水体(用户精度高20.00%)的分类效果较好。
由于喀斯特地貌条件下的农耕区普遍为梯田,故面向对象的SVM分类方法应于该地貌区的农耕区地物的识别是适宜的。后续研究可结合多时相分析方法以来进一步提高分类精度。
参考文献
Assessment of a micro-UAV system for microwave tomography Radar imaging
[J].
北斗辅助无人机航摄影像的空中三角测量
[J].
BeiDou-supported aerotriangulation for UAV aerial images
[J].
无人机航摄质量控制与分析
[J].
Quality control and analysis of UAV aerial photography
[J].
轻小型无人机航摄技术现状及发展趋势
[J].
Aerial photogrammetric technology of light small UAV:Status and trend of development
[J].
Radiometric quality assessment of images acquired by UAV’s in various lighting and weather conditions
[J].
移动测量、地理信息更新与城市管理智能化
[J].
Mobile mapping,geographic information update and urban management intelligence
[J].
基于引力搜索算法的植保无人机三维路径规划方法
[J].
3D path planning approach based on gravitational search algorithm for sprayer UAV
[J].
UAV photogrammetry application to the monitoring of rubble mound breakwaters
[J].
Automated co-registration and calibration in SFM photogrammetry for landslide change detection
[J].
Evaluation of riparian condition of Songhua River by integration of remote sensing and field measurements
[J].
Comparison of object-based and pixel-based random forest algorithm for wetland vegetation mapping using high spatial resolution GF-1 and SAR data
[J].
Improving UAV imaging quality by optical sensor fusion:An initial study
[J].
Improving classification accuracy of spectrally similar land covers in the rangeland and plateau areas with a combination of WorldView-2 and UAV images
[J].
基于遥感数据提取建筑物的主动轮廓改进模型
[J].
Improved active contour model to extract buildings based on remotely sensed data
[J].
Spatial and structured SVM for multilabel image classification
[J].
A spectral-spatial SVM-based multi-layer learning algorithm for hyperspectral image classification
[J].
A probabilistic SVM approach for hyperspectral image classification using spectral and texture features
[J].
A self-trained ensemble with semisupervised SVM:An application to pixel classification of remote sensing imagery
[J].
A novel semisupervised SVM for pixel classification of remote sensing imagery
[J].
LS-SVM-based image segmentation using pixel color-texture descriptors
[J].
Urban land cover dynamics and their impact on ecosystem services in Kigali,Rwanda using multi-temporal Landsat data
[J].
Object-oriented classification of high-resolution remote sensing imagery based on an improved colour structure code and a support vector machine
[J].
Object-oriented SVM classifier for ALSAT-2A high spatial resolution imagery:A case study of algiers urban area
[C]//
基于Landsat8 OLI影像纹理特征的面向对象土地利用/覆盖分类
[J].
Object-oriented land use/cover classification based on texture features of Landsat8 OLI image
[J].
A pattern-based definition of urban context using remote sensing and GIS
[J].
Seed growing for interactive image segmentation using SVM classification with geodesic distance
[J].
Accurate ortho-mosaicked six-band multispectral UAV images as affected by mission planning for precision agriculture proposes
[J].
Automatic identification of agricultural terraces through object-oriented analysis of very high resolution DSMs and multispectral imagery obtained from an unmanned aerial vehicle
[J].
Convolutional SVM networks for object detection in UAV imagery
[J].
Unmanned aerial vehicles (UAV) for assessment of qualitative classification of Norway spruce in temperate forest stands
[J].
An introduction interval kernel-based methods applied on support vector machines
[C]//
基于SVM的湖泊咸度等级遥感信息提取方法——以内蒙古巴丹吉林沙漠为例
[J].
Extraction of remote sensing information for lake salinity level based on SVM:A case from Badain Jaran desert in Inner Mongolia
[J].
/
〈 |
|
〉 |
