高级检索

基于无人机遥感的煤矸石山植被分类

周涛, 胡振琪, 阮梦颖, 刘曙光, 张驭航

周 涛,胡振琪,阮梦颖,等. 基于无人机遥感的煤矸石山植被分类[J]. 煤炭科学技术,2023,51(5):245−259

. DOI: 10.13199/j.cnki.cst.2021-0899
引用本文:

周 涛,胡振琪,阮梦颖,等. 基于无人机遥感的煤矸石山植被分类[J]. 煤炭科学技术,2023,51(5):245−259

. DOI: 10.13199/j.cnki.cst.2021-0899

ZHOU Tao,HU Zhenqi,RUAN Mengying,et al. Classification of coal gangue pile vegetation based on UAV remote sensing[J]. Coal Science and Technology,2023,51(5):245−259

. DOI: 10.13199/j.cnki.cst.2021-0899
Citation:

ZHOU Tao,HU Zhenqi,RUAN Mengying,et al. Classification of coal gangue pile vegetation based on UAV remote sensing[J]. Coal Science and Technology,2023,51(5):245−259

. DOI: 10.13199/j.cnki.cst.2021-0899

基于无人机遥感的煤矸石山植被分类

基金项目: 

国家重点研发计划资助项目(2019YFC1805003)

详细信息
    作者简介:

    周涛: (1997—),男,四川德阳人,硕士。E-mail: cumtzt@foxmail.com

    通讯作者:

    胡振琪: (1963—),男,安徽五河人,教授,博士生导师。E-mail: huzq1963@163.com

  • 中图分类号: P237;TP751

Classification of coal gangue pile vegetation based on UAV remote sensing

Funds: 

National Key Research and Development Program of China (2019YFC1805003)

  • 摘要:

    植被种类的准确分类是实现煤矸石山植被修复效果评价的基础。利用无人机遥感技术获取不同季节的煤矸石山可见光影像,通过色彩空间转换和纹理滤波充分挖掘可见光影像中丰富的色彩、结构及纹理等特征;然后对传统人工特征选取方法做出改进,该方法可快速、简单、高效地筛选特征信息以获取最优分类特征,并将优选结果与RGB影像融合以获得多特征融合影像;最后,利用3种监督分类模型分别对两期RGB影像及多特征融合影像进行分类处理并对结果进行精度评价及植被动态变化分析。结果表明:基于改进的人工特征选取方法可筛选出不同季节煤矸石山影像的最优分类特征,所选特征不仅能有效反映各类地物的差异性,同时可降低特征信息冗余以提高影像分类精度及效率。支持向量机(Support Vector Machine Classification, SVM)分类方法结合多特征融合影像的分类精度最高,总体分类精度最高可达90.60%,相应Kappa系数为0.8780,较同期RGB影像分别提高了9.74%和0.1265;而最大似然(Maximum Likelihood Classification, MLC)和神经网络(Neural Network, NNC)分类方法精度提高较少,总体分类精度较同期RGB影像最多可分别提高6.95%和3.93%,相应Kappa系数分别提高0.0845和0.0541。同时,基于最优分类结果从植被覆盖度和植被配置模式2个角度对常村煤矸石山植被修复效果进行评价,结果表明:该煤矸石山采用了多种不同的植被配置模式,且秋夏两季的植被覆盖度均高于75%,植被修复的整体效果较好。研究可为基于无人机可见光影像的煤矸石山植被信息识别分类提供参考,同时为煤矸石山植被修复的后期管理、维护等提供意见或建议。

    Abstract:

    The accurate classification of vegetation species is the basis for the evaluation of vegetation restoration effect of coal gangue pile. In this paper, the visible image of coal gangue pile in different seasons was obtained by UAV remote sensing technology. The color space conversion and texture filtering were used to adequately explore the rich features of color, structure and texture in the visible image. Then, the traditional artificial feature selection method was improved, which could quickly, simply and efficiently screen features information to obtain the optimal classification features, and the optimized results were fused with RGB images to obtain multi-feature fusion images. Finally, based on two stages of RGB images and multi-feature fusion images, the vegetation of coal gangue pile was classified by three supervised classification methods, including support vector machine (SVM), maximum likelihood (ML) and neural network (NN). Meanwhile, the accuracy of classification results was evaluated by confusion matrix and the dynamic changes of vegetation were analyzed. The results showed that the improved artificial feature selection method could screen out the optimal classification features of coal gangue pile vegetation in different seasons. The selected classification features can not only effectively reflect the differences of various ground features, but also reduce the redundancy of feature information to improve the accuracy and efficiency of image classification. The classification result based on Support Vector Machine Classification (SVM) combined with multi-feature fusion image had highest classification accuracy, and the overall classification accuracy could reach 90.60%, and the corresponding Kappa coefficient is 0.8780, which was 9.74% and 0.1265 higher than that of RGB image of the same period, respectively. And, the accuracy of MLC and NNC classification methods was less improved. Compared with the RGB images of the same period, the overall classification accuracy could be improved by 6.95% and 3.93%, respectively, and the corresponding Kappa coefficient could be improved by 0.0845 and 0.0541, respectively. At the same time, based on the result of optimal classification, this paper evaluated the vegetation restoration effect of coal gangue pile in Changcun from the perspectives of vegetation coverage and vegetation allocation pattern. The results showed that a variety of different vegetation allocation patterns were adopted by the coal gangue pile, and the vegetation coverage in autumn and summer is higher than 75%. The overall effect of vegetation restoration was better. This study could provide reference for the identification and classification of coal gangue piles vegetation information based on UAV visible light image, and meanwhile provide opinions or suggestions for the later management and maintenance of coal gangue piles vegetation restoration.

  • 煤炭资源的过度开采引发诸多环境问题,严重影响矿区周边地区的生态环境[1]。煤矸石山在自然环境中极易发生自燃、淋溶、扬尘等,对周边大气、水体、土壤等造成严重污染[2-4]。然而,煤矸石山的自然修复过程极其缓慢,因此采用人工修复的方法治理煤矸石山显得尤为必要。目前,煤矸石山的人工修复主要侧重于土壤基质改良、立地条件[5]、整形整地技术、优势树种与种间配置、植物栽植与绿化技术[6-7]等方面的研究,其目的是为了建立一个稳定、高效的人工植被生态系统。其中,植被修复效果的准确评价也成为治理工作中必不可少的环节。但目前植被修复效果评价大多基于实地样方调查等传统方法,不仅耗费大量人力物力及时间成本,而且无法从整体上对修复效果做出准确判断。植被种类的准确分类可在一定程度上反映植被的多样性、覆盖度及空间分布情况等。因此,快速、准确的植被信息识别分类对煤矸石山植被修复效果评价意义重大。

    传统卫星遥感技术存在分辨率较低、时效性较差、影像成本高等问题,并不适用于小区域范围内快速、低成本、高精度的地物识别分类。而无人机遥感技术凭借其高分辨率、强时效性、低成本且操作简单灵活等优势[8]广泛应用于精准农业[9-13]、林业资源调查[14-16]、地质灾害识别与监测[17-18]、矿区环境监测与土地复垦[19-21]等领域,并取得显著成果。

    目前,植被信息的识别分类方法主要分为2类。①植被指数结合阈值分割等传统方法提取植被信息,如汪小钦等[22]、高永刚等[23]利用可见光影像有限的波段信息构造不同植被指数并结合不同阈值分割方法有效提取植被信息。此类研究已较为成熟,且多针对背景与前景目标差别明显的情况,如植被与非植被分类,并未对植被种类做进一步细分。②通过色彩空间转换和纹理滤波等图像增强技术充分挖掘高分辨率影像中潜在的色彩、结构及纹理等特征并进行特征优选,然后结合不同分类方法识别植被信息,如LU等[24]根据不同时期草地的多光谱影像提取不同特征并进行特征组合,然后采用面向对象方法结合随机森林算法实现草种分类,结果表明该方法能精准识别草地中的不同物种。郭鹏等[25]对无人机可见光影像进行色彩空间转换和纹理滤波以构建影像特征集,然后通过直接比较差异系数和变异系数获得最优分类特征并利用传统分类方法实现作物分类,同时还利用植被指数区分不同作物类型,最后比较两种分类方法精度,结果表明基于色彩和纹理特征提取的农作物信息精度更高;同样,韩文霆等[26]也采用该方法筛选适宜提取玉米种植信息的特征并最终有效提取拔节期玉米种植信息。戴建国等[27]通过ReliefF-Pearson特征降维方法对可见光影像的色彩、结构及纹理等特征进行优选,然后结合不同监督分类方法对北疆主要农作物进行分类,结果表明该方法可有效区分农作物类型。此类研究其对象多呈现块状、条带状等规则分布的情况,如大田农作物分类等;而缺乏植被种类较多且空间分布复杂的识别分类研究,如煤矸石山植被等。同时,人工特征优选的过程较为复杂,且无法在优选过程中对所选特征及时进行评价,导致特征信息冗余及分类处理的精度和效率降低。

    笔者针对种类较多且空间分布复杂的植被信息识别分类所存在的问题,结合已有研究成果,以常村煤矸石山为典型区域,对传统人工特征选取方法进行改进以简单、快速、高效地获取最优分类特征,并结合多种机器学习算法对煤矸石山植被进行识别分类研究,以期获得一种基于无人机遥感的煤矸石山植被分类方法,为煤矸石山植被修复效果评价提供参考依据。

    潞安矿区位于沁水盆地东部边缘中段,隶属华北地层区山西地层分区宁武−临汾小区[28],地跨长治、襄垣、屯留、潞城和长子等五市县,包括常村、王庄、漳村、石圪节及五阳等大型煤矿。选择常村煤矿矸石山作为研究区域,其位于长治市屯留县境内(112°58′21″E,36°21′22″N),北靠常村矿区工业广场,西倚二广高速,东临208国道,其概略位置如图1a所示。该煤矸石山于1986年启用,由地面起堆,高约55 m,坡度在21°~43°,总堆量约300万t,东西长260 m,南北宽230 m,总占地面积约60000 m2。截止2009年8月完成全山绿化治理工程,无自燃迹象,植被生长态势较好,植被覆盖度较高,可作为煤矸石山植被修复的典型案例进行分析。

    图  1  研究区域
    Figure  1.  Study area

    无人机影像数据采集以大疆M210 RTK V2型四旋翼无人机作为飞行平台,其水平和垂直定位精度达到±0.1 m;同时,该平台搭载禅思(ZENMUSE) XT2双光热成像镜头,有效像素为1200万。影像采集分两次进行,时间分别为2020年10月21日和2021年6月9日,天气状况均良好,适宜进行无人机飞行试验。为避免后期分类处理结果产生误差,2次影像采集的航飞参数均保持一致,其中航向重叠度和旁向重叠度均设置为85%,飞行高度70 m,飞行速度4.0 m/s,飞行时间约75 min,共采集2期数据,总计5098景影像。

    对获取的无人机影像进行预处理,主要包括影像拼接和影像裁剪。其中,利用Pix4D mapper软件中的流程化处理模块进行影像拼接;使用ENVI5.3软件中的地表真实感兴趣区(Ground Truth ROIs)工具进行影像裁剪。研究对可见光影像各波段中心波长及范围无特殊要求,故无需进行辐射定标;同时,无人机定位精度达到厘米级,且研究对影像地理坐标无特殊要求,预处理影像点位精度足以满足要求,故无需布设地面像控点及地理坐标纠正。经过上述影像预处理,最终获取研究区域2020-10-21和2021-06-09两期正射影像,其空间分辨率约0.03 m,如图1b所示。

    所获取的2期不同季节煤矸石山可见光影像只包含红、绿、蓝3个波段信息,若仅依靠RGB色彩信息难以有效区分煤矸石山植被种类。因此,首先借助图像增强技术充分挖掘可见光影像中潜在的色彩、结构及纹理等特征,由此构建影像分类特征集;然后选取各类地物训练样本并统计样本各项特征的变异系数及差异系数;再对传统人工特征选取方法做出改进以简单、快速、高效地获取不同季节煤矸石山可见光影像的最优分类特征,并将优选结果与同期RGB影像融合以获得两期多特征融合影像;最后利用3种监督分类方法分别对2期RGB影像及多特征融合影像进行分类处理并对结果进行精度评价及植被动态变化分析。同时,基于最优分类结果分别计算不同季节煤矸石山的植被覆盖度,从植被覆盖度和植被配置模式2个角度分析煤矸石山植被修复效果。总体研究技术路线如图2所示。

    图  2  研究技术路线
    Figure  2.  Research technology roadmap

    选取合适的颜色空间是有效分割图像的基础[29]。无人机可见光影像采用RGB颜色空间显示,但红、绿、蓝波段之间存在很强的相关性及随机噪声,若直接用于分类处理难以保证精度[30]。而HSI颜色空间采用色度(Hue)、饱和度(Saturation)、亮度(Intensity)3个分量来表征颜色[31],且各分量之间相互独立,与人感受彩色的方式紧密联系[32],有利于提高分类处理精度。

    基于ENVI5.3环境分别对两期无人机可见光影像进行色彩空间转换以获得HSI色彩空间影像,其中包括色度、饱和度及亮度3个特征波段信息,如图3a所示。

    图  3  图像增强结果
    Figure  3.  Results of image enhancement

    无人机可见光影像光谱分辨率较低,但其空间分辨率较高,结构和纹理信息丰富。因此,可借助纹理滤波充分挖掘影像中丰富的结构、纹理等特征[33],并将其应用于影像分类以提高分类处理精度。

    在ENVI5.3环境下采用基于二阶概率统计的滤波方法对HSI颜色空间影像进行滤波处理以获取色调、饱和度和亮度3个分量的均值、方差、协同性、对比度、相异性、信息熵、二阶矩和相关性共计24项纹理特征。综合考虑影像处理的精度及效率,将滤波窗口设置为7×7,空间相关性矩阵XY的变化值均设置为1,灰度量化级别设置为64。纹理滤波结果如图3b所示,其中RGB通道分别为色度二阶矩、亮度相异性和饱和度方差。

    通过实地样方调查,并综合考虑无人机可见光影像的空间分辨率、地物类型差异性程度、地物空间分布复杂程度以及优势植被种类变化等因素,将2020-10-21期影像地物大致分为紫穗槐、铺地柏、草本、紫藤、乔木、枯萎高羊茅、高羊茅、裸地、建筑物及道路等共计10种类型;将2021-06-09期影像地物大致分为紫穗槐、铺地柏、草本、乔木、紫锦木、金叶女贞、紫花苜蓿、裸地、建筑物及道路等共计10种类型。为使各类地物的色彩、结构及纹理特征的统计结果更具代表性,在选取训练样本时尽量选择各类地物内部纯净像元,避免混合像元对后续分类产生影响。同时,选取足够数量且分布均匀的测试样本用于后续精度评价。训练样本及测试样本选取结果见表1

    表  1  训练样本及测试样本选取结果
    Table  1.  Selection results of training samples and test samples
    影像期数地物类型训练样本/个测试样本/个总计
    2020-10-21草本50150200
    紫穗槐100700800
    铺地柏50200250
    乔木50100150
    高羊茅100600700
    枯萎高羊茅100500600
    紫藤305080
    建筑物102030
    道路50200250
    裸地50150200
    2021-06-09紫锦木20100120
    金叶女贞10250260
    乔木30150180
    草本30100130
    紫花苜蓿150700850
    铺地柏50250300
    紫穗槐80600680
    裸地30250280
    建筑物102030
    道路50200250
    总计105052906340
    下载: 导出CSV 
    | 显示表格

    变异系数可有效反映样本内部像元统计特征值的离散程度,值越小表明该特征越能有效代表此类地物,即变异系数值越小的特征可作为此类地物区别于其他地物的典型特征;而差异系数则可反映各类地物间该特征的差异程度,值越大表明该特征越有利于区分各类地物。基于此,分别统计基于2期影像所选训练样本各项特征的平均值、标准差,并分别计算各类样本27项特征的变异系数及类间差异系数,以便用于后续多特征优选,计算公式如下:

    $$ V = \frac{S}{M} \times 100\% $$ (1)
    $$ {D_{\rm{W}}} = \frac{{{M_1} - {M_2}}}{{{M_2}}} \times 100\% $$ (2)

    式中:V为变异系数;S为标准差; M为均值;DW为类间差异系数;M1为第一类均值;M2为第二类均值。2020-10-21期影像计算结果中,草本各项特征的变异系数及草本与其他地物各项特征的差异系数计算结果见表2

    表  2  草本各项特征的变异系数及草本与其他地物各项特征的差异系数
    Table  2.  Coefficients of variation for each feature of herb and coefficient of difference for each feature among herb and other ground features
    特征指标变异系数 差异系数Dw
    紫穗槐铺地柏乔木紫藤高羊茅枯萎高羊茅建筑道路裸地
    色度17.2548.6511.1420.5361.8550.76368.3134.9146.708.18
    亮度14.327.6013.2110.871.0111.8726.5848.4035.6726.33
    饱和度15.379.0766.0018.6544.3269.8313.7582.1221.497.65
    色度均值10.8348.8310.9319.3561.6353.72404.8534.5747.0210.98
    色度方差15.5095.5977.6881.5097.3580.7989.8096.7361.7499.74
    色度协同性18.7515.6316.9818.601.8811.051.3126.747.251.49
    色度对比度10.9797.9585.5266.6398.2379.0495.5798.5442.0499.88
    色度相异性169.4966.2053.5852.9345.9642.6522.3381.1030.1893.14
    色度信息熵4.6936.5930.8436.5311.6131.6314.4322.9114.9718.61
    色度二阶矩75.28165.1378.31179.0326.86137.8824.020.4217.8328.44
    色度相关性32.921.551.248.977.779.9816.3438.1790.2623.97
    亮度均值19.747.2521.7712.691.9216.4134.8058.4145.4434.70
    亮度方差38.7016.4941.6312.3431.3356.4514.1522.20236.9572.02
    亮度协同性7.4041.460.1726.0617.8840.0933.3771.0669.1040.87
    亮度对比度11.03170.414.3369.8927.07228.06125.88178.10806.49243.96
    亮度相异性42.3373.791.4634.7215.3184.6455.13172.74301.5389.02
    亮度信息熵29.902.945.960.465.988.123.9093.1762.549.96
    亮度二阶矩39.0013.2459.469.2659.4632.1819.1897.9692.8839.80
    亮度相关性23.4385.0078.4280.3341.5581.4480.9881.8781.1480.26
    饱和度均值7.0820.5370.3224.7337.1976.9313.7683.0023.377.66
    饱和度方差54.2482.8862.6460.6082.5430.8777.0799.5220.0237.42
    饱和度协同性142.083.459.242.8046.9820.4317.6038.6416.650.93
    饱和度对比度14.4212.6779.370.5770.10255.4152.0999.38100.461.27
    饱和度相异性54.607.0331.293.1649.79102.7731.0891.0260.110.91
    饱和度信息熵11.0523.5513.8612.7735.8818.7927.7526.8533.239.70
    饱和度二阶矩57.08117.4228.7124.32408.4734.45177.7553.4858.9725.29
    饱和度相关性27.6661.5615.4250.0136.9353.0450.3729.0814.0341.55
    注:表中数据为2021-10-21期影像,差异系数DW的计算结果取绝对值。
    下载: 导出CSV 
    | 显示表格

    分析表2可知,草本各项特征的变异系数值差异明显,其中色度相异性的变异系数值为169.49%,而色度信息熵的变异系数值为4.69%;同时,草本与其他地物各项特征的差异系数值变化范围较大,其中草本与道路的亮度对比度差异系数值为806.49%,而草本与建筑物的色度二阶矩差异系数值为0.42%。上述结果表明并不是所有的特征信息都能突出各类地物之间的差异性,因此需要进一步筛选出有利于分类处理的特征信息。

    样本可分离性作为训练样本的评价指标可用于表征训练样本之间的差异性程度,其理论值在0~2.0,若该值大于1.9,说明样本之间差异大、可分离性好,属于合格样本;若小于1.8,则属于不合格样本,需要重新选择样本;若小于1.0,则考虑将两类样本合并为一类样本。当训练样本的选取无法继续优化,即不能通过改变样本的选择来增强类别间的差异性程度时,可增加样本内部所含特征信息以增强类别间的差异性程度。

    传统人工特征选取方法通过直接比较变异系数和差异系数筛选最优分类特征[25,34],但仅在分类对象类型较少时效果较好;若分类对象类型较多,该方法的选取过程会更加复杂且效率有所降低。基于此,对传统人工特征选取方法做出适当改进,具体的方法流程如下:

    1)基于原始RGB影像计算各类地物训练样本的可分离性值,并将可分离性值小于1.9的类别组合按(0,1.0)、(1.0,1.8)及(1.8,1.9)3个区间进行分组。

    2)针对可分离性值在最小区间内的类别组合,依次选取这些类别组合差异系数最大的特征信息,并将所选特征与RGB影像融合,获得多特征融合影像。

    3)基于多特征融合影像再次计算各类地物训练样本的可分离性值,并将可分离性值小于1.9的类别组合按上述方法进行分组。

    4)针对可分离性值在最小区间内的类别组合,继续依次选取这些类别组合差异系数最大的特征信息,若该特征在第二步中已被选取,则按差异系数值降序依次选取其他特征,并将所选特征与上述多特征融合影像再次融合以获得新的多特征融合影像。

    5)以此类推,不断选取并融合色彩、结构及纹理等特征以获得新的多特征融合影像,然后基于新的多特征融合影像计算训练样本的可分离性值,直到所有的样本可分离性值均大于1.9时,特征优选结束。

    对上述改进方法做如下4点说明:①与直接比较所有样本各项特征信息的变异系数及差异系数的传统人工特征选取方法不同,基于原始影像及多特征融合影像多次计算样本可分离性值可快速确定哪些地物种类无法进行有效区分,特征选取更具针对性且效率有所提高;②该方法将各步骤中样本可分离性值的计算结果按区间进行分组,且每次仅优先针对最小区间内的样本组合选取对应特征,尽可能降低特征信息冗余;③通过不断计算多特征融合影像的样本可分离性值,可对上一步骤中所选取的特征信息及时做出评价,即验证所选取的特征信息是否有助于提高样本可分离性,抑或是造成特征信息冗余以便及时做出修正;④最小区间并非固定不变,而是根据样本可分离性值的计算结果确定。

    根据上述改进的人工特征选取方法,基于2020−10−21期和2021−06−09期影像分别筛选出部分色彩、结构及纹理特征作为2期影像的最优分类特征。其中,2020−10−21期影像的最优分类特征选取结果见表3

    表  3  基于改进的人工特征选取方法的分类特征优选结果
    Table  3.  The results of classification features priority selection based on improved artificial features selection method
    样本可分离性值区间样本类别组合RGB影像特征波段优选结果
    RGB/H-Con/
    S-Var/S-Sec
    RGB/H-Con/S-Var/
    S-Sec /I-Sec
    RGB/H-Con/S-Var/
    S-Sec /I-Sec /H-Var/S-Con
    (0,1.0)乔木/紫穗槐0.506521.828821.862381.94924
    高羊茅/紫穗槐0.699361.905301.908221.95325
    高羊茅/乔木0.867071.841471.887031.94388
    乔木/铺地柏0.880971.865411.936991.97809
    (1.0,1.8)铺地柏/紫穗槐1.321961.932841.987731.99751
    高羊茅/铺地柏1.340701.681571.922281.97568
    草本/紫穗槐1.411741.980511.987141.99943
    裸地/枯萎高羊茅1.528591.917911.933911.98570
    草本/铺地柏1.617591.953591.994181.99755
    乔木/草本1.701911.911271.938761.97885
    高羊茅/枯萎高羊茅1.784801.974961.976401.99099
    道路/建筑物1.797361.999741.999942.00000
    (1.8,1.9)高羊茅/草本1.842791.983141.986071.99772
    枯萎高羊茅/紫穗槐1.847621.930091.940931.97808
    注:H-Hon、S-Var、S-Sec、 I-Sec、 H-Var、S-Con分别代表色度对比度、饱和度方差、饱和度二阶矩、亮度二阶矩、色度方差和饱和度对比度;表中样本可分离性值区间仅适用于基于RGB影像的可分离性值计算结果。
    下载: 导出CSV 
    | 显示表格

    表3所示,按照上述改进的人工特征选取方法流程依次筛选2020−10−21期影像的最优分类特征,具体步骤如下:①第一步:基于RGB影像计算各类地物训练样本的可分离性值,选取值小于1.9的类别组合并将其分为3个区间;②第二步:针对可分离性值位于最小区间(0,1.0)中的4种类别组合,依次选择差异系数值最大的特征,其中乔木/紫穗槐为色度对比度,高羊茅/紫穗槐、高羊茅/乔木均为饱和度二阶矩,乔木/铺地柏为饱和度方差,然后将这3类特征信息与RGB影像融合以获得多特征融合影像;③第三步:基于第二步获取的多特征融合影像继续计算训练样本的可分离性值并将结果按上述方法进行分组,针对可分离性值位于最小区间(1.0,1.8)中的类别组合,依次选择差异系数值最大的特征,其中高羊茅/铺地柏为亮度二阶矩,然后将该特征信息与第二步获得的多特征融合影像再次融合以获得新的多特征融合影像;④第四步:基于第三步所获取的新的多特征融合影像继续计算训练样本的可分离性值并将结果按上述方法进行分组,针对可分离性值位于最小区间(1.8,1.9)中的2种类别组合,依次选择差异系数值最大的特征,其中乔木/紫穗槐和高羊茅/乔木2种类别组合差异系数值最大的特征在第二步中已被选取,则按降序依次选择色度方差及饱和度对比度,然后将这2类特征信息与第三步获得的多特征融合影像再次融合以获得最新的多特征融合影像;⑤第五步:基于第四步所获取的最新的多特征融合影像再次计算训练样本的可分离性值,此时所有样本的可分离性值均大于1.9,表明所有训练样本均为合格样本,特征优选结束。

    根据上述2020-10-21期影像最优分类特征的筛选过程绘制多特征优选过程示意如图4所示。

    图  4  多特征优选过程
    虚线框表示最小区间内的类别组合
    Figure  4.  Process of multi-feature priority selection

    图4可知,基于改进的人工特征选取方法依次选取并融合特征信息后,目标类别组合的样本可分离性值均增大,表明所选取的分类特征可有效增强样本间的可分离性。最终,所有的样本可分离性值均大于1.9,即通过增加样本内部所含特征信息,不合格或较差的训练样本均成为合格样本。

    基于上述改进的人工特征选取方法,最终选择色度对比度、饱和度方差、饱和度二阶矩、亮度二阶矩、色度方差以及饱和度对比度等共计6种特征信息作为2020-10-21期影像的最优分类特征;选择亮度相关性、亮度对比度、色度对比度、亮度方差、饱和度对比度、饱和度二阶矩、色度方差以及饱和度相异性等共计8种特征信息作为2021-06-09期影像的最优分类特征。然后,将上述优选特征分别与相应RGB影像融合以获得2期多特征融合影像,并用于后续分类处理。

    基于ENVI5.3环境,选择3种经典机器学习分类器进行影像分类处理,包括神经网络(Neural Net Classification, NNC)、支持向量机(Support vector machine classification, SVM)和最大似然(Maximum Likelihood Classification, MLC)。同时,借助交叉验证不断优化各分类器参数。

    无人机可见光影像空间分辨率达厘米级,可认为影像中无混合像元。基于此,可计算煤矸石山的植被覆盖度。首先,选择最优分类结果进行类别合并,将所有不同类型的植被合并为植被,其余地物合并为非植被。然后,基于合并结果计算植被覆盖度,即植被像元数与影像像元总数的比值[35],计算公式如下:

    $$ C_{\rm{V}} = \frac{{{N_{{\rm{veg}}}}}}{{{N_{{\rm{total}}}}}} \times 100\% $$ (3)

    式中:CV为植被覆盖度;$ {N}_{{\rm{veg}}} $为植被像元数;$ {N}_{{\rm{total}}} $为影像像元总数。

    基于2期不同季节的煤矸石山可见光影像,分别采用上述3种监督分类方法进行地物分类处理,分类结果如图5所示。

    图  5  RGB影像分类结果
    A~G、a~j—误分类出现的区域
    Figure  5.  Classification results based on RGB image

    图5所示,以研究区域正射影像及实地调查结果为参照,两期RGB影像分类结果中均出现较多错分误差和漏分误差。2020-10-21期影像分类结果中,A、G区域将道路错分为建筑物,D区域将3处建筑物均错分为道路;B、E区域将大部分高羊茅错分为枯萎高羊茅;C区域几乎将所有紫穗槐错分为道路和枯萎高羊茅;F区域将大部分紫穗槐错分为枯萎高羊茅。此外,NNC方法分类结果中道路、枯萎高羊茅和其他地物类型之间出现大量错分现象;SVM方法分类结果中大量紫穗槐未被有效区分。2021-06-09期影像分类结果中,a、c区域将部分紫穗槐错分为铺地柏和紫锦木;d、e区域将部分道路错分为建筑物;b、h区域未能有效识别部分草本及乔木;f、j区域将部分铺地柏错分为紫花苜蓿;g、i区域则将部分紫花苜蓿错分为铺地柏。上述结果分析表明,利用RGB影像进行植被分类会产生较多错分和漏分误差,分类效果并不理想。

    基于利用改进的人工特征选取方法所获得的两期多特征融合影像,分别采用上述3种监督分类方法进行地物分类,分类结果如图6所示。

    图  6  多特征融合影像分类结果
    A~E、a~g—误分类出现的区域;F~G、h~j—对照区域
    Figure  6.  Classification results based on multi-feature fusion image

    图6所示,以研究区域正射影像及实地调查结果为参照,基于多特征融合影像的分类效果总体上优于RGB影像,但也出现了部分错分误差和漏分误差。2020-10-21期影像分类结果中,A区域将部分紫穗槐错分为乔木,C区域将部分高羊茅错分为紫穗槐,而F区域并无明显错分和漏分误差;E区域将大量高羊茅错分为枯萎高羊茅;B、D区域将少量的高羊茅错分为裸地和枯萎高羊茅,而G区域并无明显错分和漏分误差。2021-06-09期影像分类结果中,a、d、e及g区域未能有效识别部分紫穗槐,而h、i及j区域则有效区分了大部分紫穗槐和紫花苜蓿;b区域将部分紫穗槐错分为紫锦木;c区域将部分铺地柏错分为紫花苜蓿;f区域未能有效识别部分草本植被。上述结果分析表明,支持向量机分类算法的分类效果最好,没有明显错分和漏分误差出现。

    以上分析均是以研究区域正射影像及实地调查结果为标准,结合人工目视解译做出的初步判断。为判断各分类结果是否准确,还应进行定量化精度评价。

    为削弱人为因素干扰,客观地对比各分类结果的精度,采用总体分类精度和Kappa系数作为精度评价指标对不同分类结果进行精度评价。具体方法如下:根据不同类型地物的空间分布情况,选择适当数量且分布均匀的地表真实感兴趣区(Ground Truth ROIs)作为测试样本,计算其与各分类结果的混淆矩阵,即可得出总体分类精度和Kappa系数。其中,测试样本选取结果见表1,混淆矩阵计算结果见表4

    表  4  植被分类精度评价结果
    Table  4.  Evaluation results of vegetation classification accuracy
    影像期数影像类型分类方法总体分类精度/%Kappa系数
    2020-10-21RGB影像支持向量机80.860.7515
    最大似然80.060.7474
    神经网络72.180.6424
    多特征融合影像支持向量机90.600.8780
    最大似然81.760.7692
    神经网络76.110.6965
    2021-06-09RGB影像支持向量机84.070.7898
    最大似然73.990.6696
    神经网络80.770.7460
    多特征融合影像支持向量机88.490.8491
    最大似然80.940.7541
    神经网络83.790.7888
    下载: 导出CSV 
    | 显示表格

    分析表4可知,不管采用何种分类方法,基于多特征融合影像的分类结果的总体分类精度和Kappa系数均高于同期RGB影像。2020-10-21期影像分类结果中,多特征融合影像结合支持向量机分类方法的精度最高,总体分类精度为90.60%,Kappa系数为0.8780,较RGB影像分别提高了9.74%和0.1265;最大似然分类方法精度次之,总体分类精度为81.76%,Kappa系数为0.7692,较RGB影像分别提高了1.16%和0.0218;而神经网络分类方法精度最低,总体分类精度为76.11%,Kappa系数为0.6965,较RGB影像分别提高3.93%和 0.0541。2021-06-09期影像分类结果中,同样是多特征融合影像结合支持向量机分类方法精度最高,总体分类精度为88.49%,Kappa系数为0.8491,较RGB影像分别提高了4.42%和0.0593;而神经网络方法次之,总体分类精度为83.79%,Kappa系数为0.7888,较RGB影像分别提高了3.02%和0.0428;最大似然分类方法精度最低,总体分类精度为80.94%,Kappa系数为0.7541,但较RGB影像分别提高了6.95%和0.0845。由于两期影像分属秋夏两季,山体植被种类、植被长势的不同会使得无人机影像中各类植被的纹理、结构、光谱等特征有所区别,进而两期影像的分类特征筛选结果并不一致,最终导致分类结果精度有所偏差,但基于改进的人工特征选取方法的植被提取精度均高于原始影像提取结果。

    上述精度评价结果表明,基于改进的人工特征选取方法所筛选的色彩、结构及纹理等特征可有效提高煤矸石山植被种类识别分类的精度。其中,多特征融合影像结合支持向量机分类方法的识别分类效果最佳。

    煤矸石山作为一种人工构建的特殊生态环境,其立地条件与自然环境有所区别,不同季节优势植被种类也相应不同。因此,基于2020−10−21期(秋季)和2021−06−09期(夏季)影像的最优分类结果统计各类地物面积占比,同时结合实地调查结果分析不同季节常村煤矸石山植被动态变化。

    表5可知:①秋季优势植被种类主要包括高羊茅、紫穗槐和铺地柏,面积分别占55.3%、9.8%和9.7%,而夏季主要包括紫花苜蓿、紫穗槐和铺地柏,面积分别占54.3%、11.6%和9.2%;②由于秋季降水减少、植被枯死,导致裸地占比高于夏季,同时紫穗槐、草本、乔木等植被的占比均低于夏季;③铺地柏作为常绿植被,生长态势较好,其占比并无较大变化。

    表  5  不同季节各类地物面积占比统计结果
    Table  5.  Statistical result of the area proportions of various ground features in different seasons
    地物类型不同影像期数各类地物面积占比/%
    2020-10-212021-06-09
    道路16.516.7
    建筑物0.10.1
    裸地7.12.8
    紫穗槐9.811.6
    草本0.40.7
    乔木0.92.6
    铺地柏9.79.2
    紫藤0.2
    高羊茅42.4
    枯萎高羊茅12.9
    紫花苜蓿54.3
    金叶女贞1.3
    紫锦木0.7
    下载: 导出CSV 
    | 显示表格

    基于上述两期影像的最优分类结果,从植被覆盖度及植被配置模式2个角度对常村煤矸石山植被修复效果进行评价分析。

    根据上述研究结果可知,多特征融合影像结合支持向量机分类方法可高精度识别不同植被信息。因此,分别选择2020−10−21期和2021−06−09期多特征融合影像的SVM方法的分类结果进行类别合并,合并结果中仅包含植被和非植被2种类型,结果如图7所示。

    图  7  植被与非植被像元提取结果
    Figure  7.  Extraction result of vegetation and non-vegetation pixels

    图7a所示,2020-10-21期影像中,煤矸石山西侧及南侧下破位出现了较多的非植被像元,结合研究区域正射影像及实地调查结果可知,该区域的部分植被枯死及植被栽植密度较低导致部分土壤裸露在外。如图7b所示,2021-06-09期影像中,煤矸石山南侧下坡位出现较多非植被像元,结合实地调查结果可知,该区域紫穗槐呈明显条带状分布,且带间并未栽植其他植被,故出现较多裸露土壤。同时,除道路及建筑物外,两期影像中植被栽植区域基本实现了植被全覆盖。基于式(3)分别计算两期影像的植被覆盖度,其中2020-10-21期为76.35%,2021-06-09期为80.46%,夏季矸石山植被覆盖度略高于秋季。上述结果分析表明,常村煤矸石山植被修复的整体效果较好。

    基于上述两期影像的最优分类结果,简要分析常村煤矸石山的植被配置模式,并以此评价植被修复效果。分别选取2020-10-21期和2021-06-09期影像各3处典型植被覆盖区域对植被配置模式进行分析,如图8所示。

    图  8  植被配置模式分析
    Figure  8.  Analysis of vegetation configuration patterns

    图8所示,矸石山体不同坡位采用了不同的植被配置模式进行植被修复且实地情况与分类结果基本相符。其中,A区域采用草、灌型配置模式,主要以栽植紫穗槐为主,并搭配少量高羊茅等草本植物;C、F区域也采用草、灌型配置模式,但主要以栽植铺地柏为主,其他草本植物较少;D区域亦采用草、灌型配置模式,但主要以条带形式栽植紫穗槐,带间种植紫花苜蓿等草本植物。以上4个区域均采用草、灌型配置模式,该模式主要用于固坡或熟化土壤,防止水土流失,增加土壤有机质,提高土壤肥力,为引入乔木打下基础。B、E区域采用纯植草型配置模式,主要栽植高羊茅、紫花苜蓿等草本植被,植被多样性较差,该模式能够起到防止降水冲刷的护坡作用,特别是矸石山边坡常采用植草护坡措施以减少水土流失。

    根据上述分析结果可知,常村煤矸石山采用覆土绿化技术进行植被修复,但不同坡位覆土厚度有所不同,植被配置模式也相应不同。其中,大部分区域受限于覆土厚度,采用草、灌型以及纯植草型配置模式,并取得较好的修复效果。

    综上所述,从植被覆盖度和植被配置模式2个角度分析可知,常村煤矸石山植被修复整体效果较好。同时,验证了利用无人机遥感技术可实现煤矸石山植被种类的快速识别分类,为植被修复效果的快速评价提供参考依据。

    研究成果表明,基于无人机遥感技术可实现煤矸石山植被信息快速、准确的识别分类,但仍存在一定的问题和局限性:

    1)基于改进的人工特征选取方法所筛选的最优分类特征对不同分类算法的精度提升有较大区别,其中SVM分类方法总体分类精度最多可提高9.74%,而NNC和MLC分类方法最多可分别提高 3.93%和6.95%。分析这种差异产生的原因,主要是由于当数据量较小时,SVM分类器性能强于NNC及MLC分类器,会取得较好的分类效果。在后续研究中,可适当增加数据量以验证改进的人工特征选取方法对不同分类算法是否具有良好普适性。

    2)改进的人工特征选取方法是否适用于其他区域的地物识别分类,如大田农作物分类、城市土地利用现状普查、森林资源调查等还有待验证。

    3)当不同类型植被的色彩、结构、纹理等特征相似时,会产生较多错分和漏分误差,后续研究可考虑加入其它特征信息,如植被指数、DSM及多光谱信息等以增强植被种类之间的差异性程度进而获得更好的精细化分类结果。

    4)研究仅获取了煤矸石山的可见光影像,影像数据类型较为单一,且并未在现场实测反映植被长势的其他参数,如叶面积指数、叶绿素含量等,仅从植被覆盖度和植被配置模式2个角度对煤矸石山植被修复效果进行评价,内容较为薄弱。在后续研究中,可采集煤矸石山的多光谱或高光谱影像并现场实测叶面积指数、叶绿素含量等数据以反映矸石山整体植被长势,最终实现对煤矸石山植被修复效果更加准确、全面的评价。

    1)色彩空间转换结合纹理滤波可有效构建煤矸石山地物分类特征信息集。

    2)改进的人工特征选取方法可简单、快速、高效地筛选出最优分类特征,所选特征不仅能够有效反映各类地物的差异性,同时可降低特征信息冗余,从而提高后续影像分类处理的精度及效率。

    3)基于改进的人工特征选取方法所获得的多特征融合影像的分类效果均优于原始RGB影像。其中,SVM分类方法结合优选特征的分类精度最高,总体分类精度最高可达90.60%,相应Kappa系数为0.8780,较同期RGB影像分别提高了9.74%和0.1265;而MLC和NNC分类方法精度提升较少,总体分类精度最多可分别提高6.95%和3.93%,相应Kappa系数分别提高0.0845和0.0541。

    4)常村煤矸石山采用了包括纯植草型和草、灌型等多种植被配置模式进行植被修复,且秋夏两季植被覆盖度均高于75%,整体修复效果较好。

  • 图  1   研究区域

    Figure  1.   Study area

    图  2   研究技术路线

    Figure  2.   Research technology roadmap

    图  3   图像增强结果

    Figure  3.   Results of image enhancement

    图  4   多特征优选过程

    虚线框表示最小区间内的类别组合

    Figure  4.   Process of multi-feature priority selection

    图  5   RGB影像分类结果

    A~G、a~j—误分类出现的区域

    Figure  5.   Classification results based on RGB image

    图  6   多特征融合影像分类结果

    A~E、a~g—误分类出现的区域;F~G、h~j—对照区域

    Figure  6.   Classification results based on multi-feature fusion image

    图  7   植被与非植被像元提取结果

    Figure  7.   Extraction result of vegetation and non-vegetation pixels

    图  8   植被配置模式分析

    Figure  8.   Analysis of vegetation configuration patterns

    表  1   训练样本及测试样本选取结果

    Table  1   Selection results of training samples and test samples

    影像期数地物类型训练样本/个测试样本/个总计
    2020-10-21草本50150200
    紫穗槐100700800
    铺地柏50200250
    乔木50100150
    高羊茅100600700
    枯萎高羊茅100500600
    紫藤305080
    建筑物102030
    道路50200250
    裸地50150200
    2021-06-09紫锦木20100120
    金叶女贞10250260
    乔木30150180
    草本30100130
    紫花苜蓿150700850
    铺地柏50250300
    紫穗槐80600680
    裸地30250280
    建筑物102030
    道路50200250
    总计105052906340
    下载: 导出CSV

    表  2   草本各项特征的变异系数及草本与其他地物各项特征的差异系数

    Table  2   Coefficients of variation for each feature of herb and coefficient of difference for each feature among herb and other ground features

    特征指标变异系数 差异系数Dw
    紫穗槐铺地柏乔木紫藤高羊茅枯萎高羊茅建筑道路裸地
    色度17.2548.6511.1420.5361.8550.76368.3134.9146.708.18
    亮度14.327.6013.2110.871.0111.8726.5848.4035.6726.33
    饱和度15.379.0766.0018.6544.3269.8313.7582.1221.497.65
    色度均值10.8348.8310.9319.3561.6353.72404.8534.5747.0210.98
    色度方差15.5095.5977.6881.5097.3580.7989.8096.7361.7499.74
    色度协同性18.7515.6316.9818.601.8811.051.3126.747.251.49
    色度对比度10.9797.9585.5266.6398.2379.0495.5798.5442.0499.88
    色度相异性169.4966.2053.5852.9345.9642.6522.3381.1030.1893.14
    色度信息熵4.6936.5930.8436.5311.6131.6314.4322.9114.9718.61
    色度二阶矩75.28165.1378.31179.0326.86137.8824.020.4217.8328.44
    色度相关性32.921.551.248.977.779.9816.3438.1790.2623.97
    亮度均值19.747.2521.7712.691.9216.4134.8058.4145.4434.70
    亮度方差38.7016.4941.6312.3431.3356.4514.1522.20236.9572.02
    亮度协同性7.4041.460.1726.0617.8840.0933.3771.0669.1040.87
    亮度对比度11.03170.414.3369.8927.07228.06125.88178.10806.49243.96
    亮度相异性42.3373.791.4634.7215.3184.6455.13172.74301.5389.02
    亮度信息熵29.902.945.960.465.988.123.9093.1762.549.96
    亮度二阶矩39.0013.2459.469.2659.4632.1819.1897.9692.8839.80
    亮度相关性23.4385.0078.4280.3341.5581.4480.9881.8781.1480.26
    饱和度均值7.0820.5370.3224.7337.1976.9313.7683.0023.377.66
    饱和度方差54.2482.8862.6460.6082.5430.8777.0799.5220.0237.42
    饱和度协同性142.083.459.242.8046.9820.4317.6038.6416.650.93
    饱和度对比度14.4212.6779.370.5770.10255.4152.0999.38100.461.27
    饱和度相异性54.607.0331.293.1649.79102.7731.0891.0260.110.91
    饱和度信息熵11.0523.5513.8612.7735.8818.7927.7526.8533.239.70
    饱和度二阶矩57.08117.4228.7124.32408.4734.45177.7553.4858.9725.29
    饱和度相关性27.6661.5615.4250.0136.9353.0450.3729.0814.0341.55
    注:表中数据为2021-10-21期影像,差异系数DW的计算结果取绝对值。
    下载: 导出CSV

    表  3   基于改进的人工特征选取方法的分类特征优选结果

    Table  3   The results of classification features priority selection based on improved artificial features selection method

    样本可分离性值区间样本类别组合RGB影像特征波段优选结果
    RGB/H-Con/
    S-Var/S-Sec
    RGB/H-Con/S-Var/
    S-Sec /I-Sec
    RGB/H-Con/S-Var/
    S-Sec /I-Sec /H-Var/S-Con
    (0,1.0)乔木/紫穗槐0.506521.828821.862381.94924
    高羊茅/紫穗槐0.699361.905301.908221.95325
    高羊茅/乔木0.867071.841471.887031.94388
    乔木/铺地柏0.880971.865411.936991.97809
    (1.0,1.8)铺地柏/紫穗槐1.321961.932841.987731.99751
    高羊茅/铺地柏1.340701.681571.922281.97568
    草本/紫穗槐1.411741.980511.987141.99943
    裸地/枯萎高羊茅1.528591.917911.933911.98570
    草本/铺地柏1.617591.953591.994181.99755
    乔木/草本1.701911.911271.938761.97885
    高羊茅/枯萎高羊茅1.784801.974961.976401.99099
    道路/建筑物1.797361.999741.999942.00000
    (1.8,1.9)高羊茅/草本1.842791.983141.986071.99772
    枯萎高羊茅/紫穗槐1.847621.930091.940931.97808
    注:H-Hon、S-Var、S-Sec、 I-Sec、 H-Var、S-Con分别代表色度对比度、饱和度方差、饱和度二阶矩、亮度二阶矩、色度方差和饱和度对比度;表中样本可分离性值区间仅适用于基于RGB影像的可分离性值计算结果。
    下载: 导出CSV

    表  4   植被分类精度评价结果

    Table  4   Evaluation results of vegetation classification accuracy

    影像期数影像类型分类方法总体分类精度/%Kappa系数
    2020-10-21RGB影像支持向量机80.860.7515
    最大似然80.060.7474
    神经网络72.180.6424
    多特征融合影像支持向量机90.600.8780
    最大似然81.760.7692
    神经网络76.110.6965
    2021-06-09RGB影像支持向量机84.070.7898
    最大似然73.990.6696
    神经网络80.770.7460
    多特征融合影像支持向量机88.490.8491
    最大似然80.940.7541
    神经网络83.790.7888
    下载: 导出CSV

    表  5   不同季节各类地物面积占比统计结果

    Table  5   Statistical result of the area proportions of various ground features in different seasons

    地物类型不同影像期数各类地物面积占比/%
    2020-10-212021-06-09
    道路16.516.7
    建筑物0.10.1
    裸地7.12.8
    紫穗槐9.811.6
    草本0.40.7
    乔木0.92.6
    铺地柏9.79.2
    紫藤0.2
    高羊茅42.4
    枯萎高羊茅12.9
    紫花苜蓿54.3
    金叶女贞1.3
    紫锦木0.7
    下载: 导出CSV
  • [1] 胡振琪,理源源,李根生,等. 碳中和目标下矿区土地复垦与生态修复的机遇与挑战[J]. 煤炭科学技术,2023,51(1):474−483.

    HU Zhenqi,LI Yuanyuan,LI Gensheng,et al. Opportunities and challenges of land reclamation and ecological restoration in mining areas under carbon neutral target[J]. Coal Science and Technology,2023,51(1):474−483.

    [2] 王绍清,张 路,赵云刚,等. 激光诱导法制备煤基石墨烯的探索及其结构特征研究[J]. 煤炭科学技术,2023,51(2):458−465.

    WANG Shaoqing,ZHANG Lu,ZHAO Yungang,et al. Exploration and structural characteristics of laser-induced preparation of coal-based graphene[J]. Coal Science and Technology,2023,51(2):458−465.

    [3] 李 松,万 洁. 煤矸石自燃机理及其防治技术研究[J]. 环境科学与技术,2005(2):82−84,119.

    LI Song,WAN Jie. Coal waste rocks: spontaneous combustion mechanism and its control[J]. Environmental Science & Technology,2005(2):82−84,119.

    [4] 秦波涛,蒋文婕,史全林,等. 矿井粉煤灰基防灭火技术研究进展[J]. 煤炭科学技术,2023,51(1):329−342.

    QIN Botao,JIANG Wenjie,SHI Quanlin,et al. Research progress on fly ash foundation technology to prevent and control spontaneous combustion of coal in mines[J]. Coal Science and Technology,2023,51(1):329−342.

    [5] 杨宇平,许 丽,丰 菲,等. 基于植被重建的宁东矿区煤矸石山立地类型及其特征[J]. 水土保持通报,2020,40(3):202−207.

    YANG Yuping,XU Li,FENG Fei,et al. Site types and characteristics of coal gangue hill based on vegetation reconstruction in Ningdong mining area[J]. Bulletin of Soil and Water Conservation,2020,40(3):202−207.

    [6] 赵方莹,刘 飞,巩 潇. 煤矸石山危害及其植被恢复研究综述[J]. 露天采矿技术,2013(2):77−81.

    ZHAO Fangying,LIU Fei,GONG Xiao. Summarize of coal gangues harm and plant renew research[J]. Opencast Mining Technology,2013(2):77−81.

    [7] 李鹏波,胡振琪,吴 军,等. 煤矸石山的危害及绿化技术的研究与探讨[J]. 矿业研究与开发,2006(4):93−96.

    LI Pengbo,HU Zhenqi,WU Jun,et al. Research and discussion on the hazards and greening technology of coal gangue dump[J]. Mining Research and Development,2006(4):93−96.

    [8] 韩文霆,张立元,张海鑫,等. 基于无人机遥感与面向对象法的田间渠系分布信息提取[J]. 农业机械学报,2017,48(3):205−214.

    HAN Wenting,ZHANG Liyuan,ZHANG Haixin,et al. Extraction method of sublateral canal distribution information based on UAV remote sensing[J]. Transactions of the Chinese Society for Agricultural Machinery,2017,48(3):205−214.

    [9] 吴智超,李长春,冯海宽,等. 基于无人机数码影像的马铃薯覆盖度提取方法[J]. 农业机械学报,2020,51(3):164−170. doi: 10.6041/j.issn.1000-1298.2020.03.019

    WU Zhichao,LI Changchun,FENG Haikuan,et al. Potato coverage extraction method based on digital image[J]. Transactions of the Chinese Society for Agricultural Machinery,2020,51(3):164−170. doi: 10.6041/j.issn.1000-1298.2020.03.019

    [10]

    DEHKORDI R H,JARROUDI M E,KOUADIO L,et al. Monitoring wheat leaf rust and stripe rust in winter wheat using high-resolution UAV-Based Red-Green-Blue Imagery[J]. Remote Sensing,2020,12:3696.

    [11]

    ŁUKASZ J,KONRAD S,WOJCIECH Ostrowski,et al. Evaluation of Rapeseed Winter Crop Damage Using UAV-Based Multispectral Imagery[J]. Remote Sensing,2020,12(16):2618.

    [12] 陈 鹏,冯海宽,李长春,等. 无人机影像光谱和纹理融合信息估算马铃薯叶片叶绿素含量[J]. 农业工程学报,2019,35(11):63−74.

    CHEN Peng,FENG Haikuan,LI Changchun,et al. Estimation of chlorophyll content in potato using fusion of texture and spectral features derived from UAV multispectral image[J]. Transactions of the Chinese Society of Agricultural Engineering,2019,35(11):63−74.

    [13]

    JIN X,LI Z,Atzberger C. Editorial for the special issue “Estimation of Crop Phenotyping Traits using Unmanned Ground Vehicle and Unmanned Aerial Vehicle Imagery”[J]. Remote Sensing,2020,12(6):940. doi: 10.3390/rs12060940

    [14]

    MIDHUN M,CARLOS A S,CARINE K,et al. Individual tree detection from unmanned aerial Vehicle (UAV) derived canopy height model in an open canopy mixed conifer forest[J]. Forests,2017,8(9):340−357. doi: 10.3390/f8090340

    [15] 杨 坤,赵艳玲,张建勇,等. 利用无人机高分辨率影像进行树木高度提取[J]. 北京林业大学学报,2017,39(8):17−23.

    YANG Kun,ZHAO Yanling,ZHANG Jianyong,et al. Tree height extraction using high-resolution imagery acquired from an unmanned aerial vehicle (UAV)[J]. Journal of Beijing Forestry University,2017,39(8):17−23.

    [16]

    SEIFERT E,SEIFERT S,VOGT H,et al. Influence of Drone Altitude, Image Overlap, and Optical Sensor Resolution on Multi-View Reconstruction of Forest Images[J]. Remote Sensing,2019,11(10):1252. doi: 10.3390/rs11101252

    [17] 彭大雷,许 强,董秀军,等. 无人机低空摄影测量在黄土滑坡调查评估中的应用[J]. 地球科学进展,2017,32(3):319−330.

    PENG Dalei,XU Qiang,DONG Xiujun,et al. Application of unmanned aerial vehicles low-altitude photogrammetry in investigation and evaluation of loess landslide[J]. Advance in Earth Science,2017,32(3):319−330.

    [18] 赵星涛,胡 奎,卢晓攀,等. 无人机低空航摄的矿山地质灾害精细探测方法[J]. 测绘科学,2014,39(6):49−52,64.

    ZHAO Xingtao,HU Kui,LU Xiaopan,et al. Precise detection method for mine geological disasters using low-altitude photogrammetry based on unmanned aerial vehicle[J]. Science of Surveying and Mapping,2014,39(6):49−52,64.

    [19]

    HU Xiao,LI Xinju,MIN Xiangyu,et al. Optimal scale extraction of farmland in coal mining areas with high groundwater levels based on visible light images from an unmanned aerial vehicle (UAV)[J]. Earth Science Informatics,2020,13(5):1151−1162.

    [20] 肖 武,陈佳乐,赵艳玲,等. 利用无人机遥感反演高潜水位矿区沉陷地玉米叶绿素含量[J]. 煤炭学报,2019,44(1):295−306.

    XIAO Wu,CHEN Jiale,ZHAO Yanling,et al. Identify maize chlorophyll impacted by coal mining subsidence in high groundwater table area based on UAV remote sensing[J]. Journal of China Coal Society,2019,44(1):295−306.

    [21]

    X HU, X LI. Information extraction of subsided cultivated land in high-groundwater-level coal mines based on unmanned aerial vehicle visible bands[J]. Environmental Earth Sciences, 2019, 78(14).

    [22] 汪小钦,王苗苗,王绍强,等. 基于可见光波段无人机遥感的植被信息提取[J]. 农业工程学报,2015,31(5):152−157, 159, 158.

    WANG Xiaoqin,WANG Miaomiao,WANG Shaoqiang,et al. Extraction of vegetation information from visible unmanned aerial vehicle images[J]. Transactions of the Chinese Society of Agricultural Engineering,2015,31(5):152−157, 159, 158.

    [23] 高永刚,林悦欢,温小乐,等. 基于无人机影像的可见光波段植被信息识别[J]. 农业工程学报,2020,36(3):178−189.

    GAO Yonggang,LIN Yuehuan,WEN Xiaole,et al. Vegetation information recognition in visible band based on UAV images[J]. Transactions of the Chinese Society of Agricultural Engineering,2020,36(3):178−189.

    [24]

    LU B,HE Y. Species classification using Unmanned Aerial Vehicle (UAV)-acquired high spatial resolution imagery in a heterogeneous grassland[J]. ISPRS Journal of Photogrammetry and Remote Sensing,2017,128:73−85. doi: 10.1016/j.isprsjprs.2017.03.011

    [25] 郭 鹏,武法东,戴建国,等. 基于无人机可见光影像的农田作物分类方法比较[J]. 农业工程学报,2017,33(13):112−119. doi: 10.11975/j.issn.1002-6819.2017.13.015

    GUO Peng,WU Fadong,DAI Jianguo,et al. Comparison of farmland crop classification methods based on visible light images of unmanned aerial vehicles[J]. Transactions of the Chinese Society of Agricultural Engineering,2017,33(13):112−119. doi: 10.11975/j.issn.1002-6819.2017.13.015

    [26] 韩文霆,李 广,苑梦婵,等. 基于无人机遥感技术的玉米种植信息提取方法研究[J]. 农业机械学报,2017,48(1):139−147.

    HAN Wenting,LI Guang,YUAN Mengchan,et al. Extraction method of maize planting information based on UAV remote sensing techonology[J]. Transactions of the Chinese Society for Agricultural Machinery,2017,48(1):139−147.

    [27] 戴建国,张国顺,郭 鹏,等. 基于无人机遥感可见光影像的北疆主要农作物分类方法[J]. 农业工程学报,2018,34(18):122−129.

    DAI Jianguo,ZHANG Guoshun,GUO Peng,et al. Classification method of main crops in northern Xinjiang based on UAV visible waveband images[J]. Transactions of the Chinese Society of Agricultural Engineering,2018,34(18):122−129.

    [28] 冯绪兴,倪小明,郝少伟,等. 常村煤矿煤层气井产水/产气曲线类型及其成因分析[J]. 中国矿业,2021,30(1):106−113.

    FENG Xuxing,NI Xiaoming,HAO Shaowei,et al. Curve types about water and gas production and their formation mechanism of coalbed methane Changcun coal mine[J]. China Mining Magazine,2021,30(1):106−113.

    [29] 庞晓敏,闵子建,阚江明. 基于HSI和LAB颜色空间的彩色图像分割[J]. 广西大学学报(自然科学版),2011,36(6):976−980.

    PANG Xiaomin,MIN Zijian,KAN Jiangming. Color image segmentation based on HIS and LAB color space[J]. Journal of Guangxi University (Natural Science Edition),2011,36(6):976−980.

    [30] 韩 洁,郭 擎,李 安. 结合非监督分类和几何—纹理—光谱特征的高分影像道路提取[J]. 中国图象图形学报,2017,22(12):1788−1797.

    HAN Jie,GUO Qing,LI An. Road extraction based on unsupervised classification and geometric-texture-spectral features for high-resolution remote sensing images[J]. Journal of Image and Graphics,2017,22(12):1788−1797.

    [31] 张学习,杨宜民. 彩色图像工程中常用颜色空间及其转换[J]. 计算机工程与设计,2008(5):1210−1212.

    ZHANG Xuexi,YANG Yimin. Common color space and its conversions in color image project[J]. Computer Engineering and Design,2008(5):1210−1212.

    [32] 孙慧贤,张玉华,罗飞路. 基于HSI颜色空间的彩色边缘检测方法研究[J]. 光学技术,2009,35(2):221−224,228.

    SUN Huixian,ZHANG Yuhua,LUO Feilu. Color edge detection based on HIS color space[J]. Optical Technique,2009,35(2):221−224,228.

    [33] 侯群群,王 飞,严 丽. 基于灰度共生矩阵的彩色遥感图像纹理特征提取[J]. 国土资源遥感,2013,25(4):26−32.

    HOU Qunqun,WANG Fei,YAN Li. Extraction of color image texture feature based on gray-level co-occurrence matrix[J]. Remote Sensing for Land and Resources,2013,25(4):26−32.

    [34] 刘 斌,史 云,吴文斌,等. 基于无人机遥感可见光影像的农作物分类[J]. 中国农业资源与区划,2019,40(8):55−63.

    LIU Bin,SHI Yun,WU Wenbin,et al. Crop classification based on UAV remote sensing images[J]. Chinese Journal of Agricultural Resources and Regional Planning,2019,40(8):55−63.

    [35] 刘 峰,刘素红,向 阳. 园地植被覆盖度的无人机遥感监测研究[J]. 农业机械学报,2014,45(11):250−257. doi: 10.6041/j.issn.1000-1298.2014.11.039

    LIU Feng,LIU Suhong,XIANG Yang. Study on monitoring fraction vegetation cover of garden plots by unmanned aerial vehicles[J]. Transactions of the Chinese Society for Agricultural Machinery,2014,45(11):250−257. doi: 10.6041/j.issn.1000-1298.2014.11.039

  • 期刊类型引用(2)

    1. 何继,崔瑞豪,李虎民,王磊,马飞,王培俊. 基于无人机的人工植被分布识别与重建效果评价. 中国煤炭. 2024(11): 142-152 . 百度学术
    2. 陈凯,雷少刚,杨星晨,史运喜,陈树召. 基于UAV-RGB的矿区植物分类及其叶绿素含量时空变化分析. 金属矿山. 2023(12): 227-233 . 百度学术

    其他类型引用(4)

图(8)  /  表(5)
计量
  • 文章访问数:  59
  • HTML全文浏览量:  13
  • PDF下载量:  29
  • 被引次数: 6
出版历程
  • 收稿日期:  2021-08-02
  • 网络出版日期:  2023-05-28
  • 刊出日期:  2023-05-30

目录

/

返回文章
返回