
基于无人机遥感影像的盐碱地信息的精准提取方法
徐存东, 李洪飞, 谷丰佑, 张鹏, 高德坤
基于无人机遥感影像的盐碱地信息的精准提取方法
An Accurate Method to Extract Saline-alkali Land Information Based on Unmanned Aerial Vehicle Remote Sensing Image
为探索在无人机遥感影像下可以对盐碱地信息精确分类的方法。选取甘肃省景泰川电力提灌灌区一期灌区为研究区,通过运用Trimble UX5固定机翼无人机采集研究区遥感数据,结合ENVI软件应用监督分类的五种分类器(平行六面体、最小距离、神经网络、最大似然、支持向量机)执行分类,对分类结果进行精度对比,并将典型区域盐碱地面积进行统计提取。两区域分类结果表明:支持向量机对于无人机遥感盐碱地信息提取法优于其他方法,区域1和区域2分类总精度分别为94.66%、96.55%,Kappa系数分别为0.936 1、0.957 3,盐碱地面积占比分别为43.5%、36.1%。支持向量机分类可为干旱荒漠灌区无人机遥感影像下盐碱地信息精确提取提供方法支撑。
Remote sensing technology has been increasingly used in agriculture. This paper presents an improved support vector machine method to estimate saline-alkali soil from imagery obtained from UAV. Taking the first irrigation district at the Jingtai Electric Power Irrigation District of Gansu Province as an example,this paper uses Trimble refused UX5 fixed-wing uav remote sensing data in the study area, combining ENVI software based on supervised classification of five types of classifier (parallelepiped, minimum distance, neural network and maximum likelihood and support vector machine) to perform classification, accuracy compared with the result of classification, statistics and typical area saline-alkali land area is extracted. The classification results of the two regions show that support vector machine (SVM) is superior to other methods in extracting saline-alkali land information from uav remote sensing. The total accuracy of region 1 and region 2 classification is 94.66% and 96.55%, the Kappa coefficient is 0.936 1 and 0.957 3, and the proportion of saline-alkali land area is 43.5% and 36.1%, respectively. SVM classification can provide method support for accurate extraction of saline-alkali land information from uav remote sensing images in arid desert irrigated areas.
无人机遥感 / 盐碱地 / 监督分类 / 景电灌区 {{custom_keyword}} /
UAV remote sensing / saline-alkali land / supervise classification / irrigation district {{custom_keyword}} /
表1 5种监督分类器说明Tab.1 Description of five supervised classifiers |
分类器 | 说明 | 参数设置 |
---|---|---|
平行六面体 | 根据训练样本的亮度值形成一个n维的平行六面体数据空间,其他像元的光谱值如果落在平行六面体任何一个训练样本所对应的区域,就被划分其对应的类别中 | 标准差阈值:None |
最小距离 | 利用训练样本数据计算出每一类样本的均值向量和标准差向量,然后以均值向量作为该类样本在特征空间的中心位置,计算每个像元到各类样本中心的距离,当该距离最小时,则将该像元归入哪一类 | 标准差阈值:4 最大距离误差:None |
神经网络 | 指用人工智能模拟的网络结构,用许多小的分类器模拟生物的神经元,通过算法表达人脑的识别、记忆以及思考过程最后将其应用于图像分类 | 活化函数:Logistic 训练贡献阈值:0.9权重调节速度:0.2 RMS误差值:0.1 |
最大似然 | 假设每一个样本数据在每一类统计中都呈正态分布,计算图像中像元在某一类训练样本下的似然度,并将该像元归到似然度最大的一类样本区间中 | 似然度阈值:None |
支持向量机(SVM) | 支持向量机分类是以统计学习理论为基础的机器学习方法,可以自动寻求对样本分类有较大区分能力的支持向量,由此构造出分类器,可以将各样本之间的阈值最大化,准确地对样本进行分类 | 核心多项式:RBF 惩罚系数:60 间隔r=0.15 |
表2 遥感影像解译标志Tab.2 Interpretation mark of remote sensing image |
类别 | 遥感影像样本 | 解译特征 |
---|---|---|
重度盐碱地 | ![]() | 影像上呈现亮白色斑块分布,中间几乎无植被生长 |
中度盐碱地 | ![]() | 影像上白色斑块较少,且亮斑面积为30%~50% |
轻度盐碱地 | ![]() | 遥感影像上植被中间有亮白色斑块,且亮斑面积较小,约10% |
未耕地块 | ![]() | 影像上呈现出深褐色 |
已耕地块 | ![]() | 影像上呈现出白色纹理,边界清晰 |
荒地 | ![]() | 影像中呈现先浅黄色,纹理较为复杂 |
树木 | ![]() | 影像呈深绿色与黑色,纹理复杂 |
表3 区域1遥感影像分类精度 (%)Tab.3 Classification accuracy of region 1 remote sensing image |
方法 | 类别 | 重度盐渍区 | 中度盐渍区 | 轻度盐渍区 | 已耕地 | 未耕地 | 树木 | 荒地 |
---|---|---|---|---|---|---|---|---|
平行六面体 (总精度:66.96% Kappa系数:0.624 6) | 生产者精度 | 79.41 | 67.30 | 69.29 | 71.59 | 56.38 | 77.00 | 52.53 |
使用者精度 | 98.33 | 98.07 | 63.44 | 98.93 | 86.85 | 71.00 | 100.00 | |
最小距离 (总精度:75.72% Kappa系数:0.708 8) | 生产者精度 | 99.76 | 56.85 | 94.23 | 97.92 | 40.43 | 77.00 | 54.45 |
使用者精度 | 97.00 | 95.14 | 52.30 | 76.58 | 49.25 | 71.00 | 100.00 | |
神经网络 (总精度:89.66% Kappa系数:0.876 1) | 生产者精度 | 98.97 | 57.43 | 100.00 | 99.76 | 94.44 | 92.96 | 85.49 |
使用者精度 | 74.14 | 99.86 | 86.10 | 93.34 | 93.01 | 100.00 | 100.00 | |
最大似然 (总精度:91.15% Kappa系数:0.894 2) | 生产者精度 | 100.00 | 72.95 | 79.94 | 97.51 | 98.97 | 99.06 | 78.99 |
使用者精度 | 97.38 | 96.38 | 86.66 | 93.30 | 96.88 | 89.41 | 100.00 | |
支持向量机 (总精度:94.66% Kappa系数:0.936 1) | 生产者精度 | 100.00 | 77.84 | 98.32 | 96.98 | 99.28 | 99.53 | 95.11 |
使用者精度 | 83.31 | 96.21 | 98.53 | 97.26 | 96.40 | 92.58 | 100.00 |
表4 区域2遥感影像分类精度 (%)Tab.4 Classification accuracy of region 2 remote sensing image |
方法 | 类别 | 重度盐渍区 | 中度盐渍区 | 轻度盐渍区 | 已耕地 | 未耕地 | 荒地 |
---|---|---|---|---|---|---|---|
平行六面体 (总精度:53.67% Kappa系数:0.462 1) | 生产者精度 | 35.23 | 44.96 | 52.73 | 66.37 | 61.19 | 37.58 |
使用者精度 | 100.00 | 100.00 | 60.53 | 99.13 | 66.13 | 100.00 | |
最小距离 (总精度:74.90% Kappa系数:0.690 0) | 生产者精度 | 80.30 | 97.15 | 56.21 | 96.98 | 56.60 | 63.38 |
使用者精度 | 100.00 | 74.96 | 58.31 | 90.54 | 61.47 | 64.19 | |
神经网络 (总精度:90.16% Kappa系数:0.877 5) | 生产者精度 | 97.16 | 87.28 | 85.02 | 99.51 | 98.51 | 42.99 |
使用者精度 | 100.00 | 100.00 | 96.57 | 96.87 | 72.41 | 91.22 | |
最大似然 (总精度:91.15% Kappa系数:0.894 2) | 生产者精度 | 100.00 | 72.95 | 97.38 | 97.51 | 98.97 | 78.99 |
使用者精度 | 79.94 | 96.38 | 86.66 | 93.30 | 96.88 | 100.00 | |
支持向量机 (总精度:96.55% Kappa系数:0.957 3) | 生产者精度 | 99.81 | 98.90 | 91.52 | 100.00 | 100.00 | 80.57 |
使用者精度 | 100.00 | 90.74 | 99.87 | 99.90 | 90.45 | 100.00 |
表5 区域1各分类结果面积占比Tab.5 Area ratio of each classification result in area 1 |
地物类别 | 平行六面体 | 最小距离 | 神经网络 | 最大似然 | 支持向量机 | |||||
---|---|---|---|---|---|---|---|---|---|---|
面积/m2 | 比例/% | 面积/m2 | 比例/% | 面积/m2 | 比例/% | 面积/m2 | 比例/% | 面积/m2 | 比例/% | |
重度盐渍区 | 11 138 | 3.7 | 17 345 | 5.7 | 24 977 | 8.3 | 27 020 | 8.9 | 27 555 | 9.1 |
中度盐渍区 | 78 894 | 26.1 | 77 218 | 25.6 | 34 932 | 11.6 | 57 716 | 19.1 | 56 714 | 18.8 |
轻度盐渍区 | 29 705 | 9.8 | 75 183 | 24.9 | 60 600 | 20.0 | 54 496 | 18.0 | 46 914 | 15.6 |
已耕地块 | 36 228 | 12.0 | 56 868 | 18.8 | 84 091 | 27.8 | 64 877 | 21.4 | 64 428 | 21.3 |
未耕地块 | 44 593 | 14.8 | 55 589 | 18.4 | 77 505 | 25.7 | 75 465 | 25 | 73 856 | 24.4 |
树木 | 17 549 | 5.8 | 15 798 | 5.2 | 11 440 | 3.8 | 16 534 | 5.5 | 22 643 | 7.5 |
荒地 | 4 620 | 1.5 | 3 720 | 1.2 | 8 176 | 2.7 | 5 614 | 1.9 | 9 611 | 3.2 |
未定义区域 | 79 331 | 26.3 | 336 | 0.1 | 336.69 | 0.1 | 336.69 | 0.1 | 336.69 | 0.1 |
表6 区域2各分类结果面积占比Tab.6 Area ratio of each classification result in area 2 |
地物类别 | 平行六面体 | 最小距离 | 神经网络 | 最大似然 | 支持向量机 | |||||
---|---|---|---|---|---|---|---|---|---|---|
面积/m2 | 比例/% | 面积/m2 | 比例/% | 面积/m2 | 比例/% | 面积/m2 | 比例/% | 面积/m2 | 比例/% | |
重度盐渍区 | 2 685 | 1.7 | 6 290 | 4.0 | 10 515 | 6.6 | 6 863 | 4.4 | 8 736 | 5.6 |
中度盐渍区 | 2 004 | 1.3 | 13 996 | 8.9 | 8 904 | 5.7 | 10 807 | 6.9 | 13 019 | 8.4 |
轻度盐渍区 | 27 905 | 17.9 | 34 205 | 22.0 | 42 014 | 27.0 | 40 025 | 25.7 | 34 471 | 22.1 |
已耕地块 | 29 053 | 18.7 | 33 482 | 21.5 | 36 295 | 23.0 | 44 058 | 28.3 | 41 244 | 26.5 |
未耕地块 | 35 290 | 22.7 | 38 351 | 24.6 | 54 605 | 35.1 | 50 352 | 32.3 | 55 036 | 35.3 |
荒地 | 1 832 | 1.2 | 29 306 | 18.8 | 3 188 | 2.0 | 3 416 | 2.2 | 3 015 | 1.9 |
未定义区域 | 56 992 | 36.6 | 240 | 0.16 | 240 | 0.16 | 240 | 0.16 | 240 | 0.16 |
1 |
汪小钦,王苗苗,王绍强,等.基于可见光波段无人机遥感的植被信息提取[J].农业工程学报,2015,31(5):152-157,159,158.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
2 |
汪沛,罗锡文,周志艳,等.基于微小型无人机的遥感信息获取关键技术综述[J].农业工程学报,2014,30(18):1-12.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
3 |
史洁青,冯仲科,刘金成.基于无人机遥感影像的高精度森林资源调查系统设计与试验[J].农业工程学报,2017,33(11):82-90.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
4 |
戴建国,张国顺,郭鹏,等.基于无人机遥感可见光影像的北疆主要农作物分类方法[J].农业工程学报,2018,34(18):122-129.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
5 |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
6 |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
7 |
张茂林.浅析景电灌区马家磨河河道治理方案优选[J].水利建设与管理,2014,34(12):40-44.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
8 |
李建国,濮励杰,朱明,等.土壤盐渍化研究现状及未来研究热点[J].地理学报,2012,67(9):1 233-1 245.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
9 |
程慧. 景电灌区不同灌溉方式下水盐运移规律研究[D].郑州:华北水利水电大学,2018.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
10 |
徐存东,张锐,王荣荣,等.基于改进支持向量机的盐碱地信息精确提取方法研究[J].灌溉排水学报,2018,37(9):62-68.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
11 |
李浩. 南水北调中线水源区生态环境脆弱性研究[D]. 郑州:华北水利水电大学,2018.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
12 |
史路路,郑柯,唐娉,等.一种面向土地覆盖分类的卷积神经网络模型[J].遥感信息,2019,34(3):34-42.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
13 |
孙坤,鲁铁定.监督分类方法在遥感影像分类处理中的比较[J].江西科学,2017,35(3):367-371,468.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
14 |
杨佳佳,姜琦刚,赵静,等.基于环境减灾卫星高光谱数据的盐碱地等级划分[J].农业工程学报,2011,27(10):118-124.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
15 |
吴晗,林晓龙,李曦嵘,等.面向农业应用的无人机遥感影像地块边界提取[J].计算机应用,2019,39(1):298-304.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
16 |
许迪,王少丽,蔡林根,等.利用NDVI指数识别作物及土壤盐碱分布的应用研究[J].灌溉排水学报,2003(6):5-8,32.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
17 |
李爽,丁圣彦,许叔明.遥感影像分类方法比较研究[J].河南大学学报(自然科学版),2002(2):70-73.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
18 |
闻静.传统遥感影像的分类[J].测绘与空间地理信息,2012,35(10):118-120.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
19 |
郑海朋. 羌塘高原典型地段土地覆被遥感信息提取方法研究[D]. 重庆:西南大学,2017.
{{custom_citation.content}}
{{custom_citation.annotation}}
|
{{custom_ref.label}} |
{{custom_citation.content}}
{{custom_citation.annotation}}
|
/
〈 |
|
〉 |