• 中国中文核心期刊
  • 中国科学引文数据库(CSCD)核心库来源期刊
  • 中国科技论文统计源期刊(CJCR)
  • 第二届国家期刊奖提名奖

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

面向三江保护区植被类型识别的最佳波段组合研究

姜怡 刘华 赵峰 马永康 刘浩栋

引用本文:
Citation:

面向三江保护区植被类型识别的最佳波段组合研究

    通讯作者: 赵峰, zhaofeng@ifrit.ac.cn
  • 中图分类号: X87

Optimal Band Combination for Vegetation Type Recognition in Sanjiang Nature Reserve

    Corresponding author: ZHAO Feng, zhaofeng@ifrit.ac.cn ;
  • CLC number: X87

  • 摘要: 目的 针对高分六号影像数据面向植被类型、特别是湿地植被信息识别的最佳波段组合进行研究,为高分六号数据在湿地资源监测中的进一步应用提供参考。 方法 以三江国家级自然保护区为对象,基于高分六号影像数据,结合当地植被特点,开展保护区植被类型识别中的最佳波段组合研究。从影像信息特征、最佳指数因子以及光谱特征曲线和地物可分性角度进行综合分析,获取基于GF6/WFV数据的三江保护区植被类型识别的最佳波段组合。同时,通过对比实验,利用支持向量机分类方法,对样区的植被类型进行信息提取,证实了结果的可靠性。 结果 波段组合4(R)-5(G)-1(B)识别植被类型精度最高。总体精度为89.15%,Kappa系数为0.8463。 结论 波段组合4(R)-5(G)-1(B)为三江保护区植被类型识别的最佳波段组合。本研究可为东北地区湿地保护区的退耕还湿工程监测提供参考。
  • 图 1  研究区高分六号影像(RGB=3,2,1)

    Figure 1.  Image of GF6 in the study area(RGB=3,2,1)

    图 2  相关系数矩阵

    Figure 2.  Image of correlation coefficient matrix

    图 3  5种植被类型的光谱特征曲线

    Figure 3.  5 Spectral characteristic curve of planting quilt type

    图 4  波段组合分类结果

    Figure 4.  Classification results of band combination

    表 1  高分六号WFV数据

    Table 1.  GF6/WFV

    项目 Project高分六号WFV数据GF6/WFV
    光谱范围 Wavelength range/ um B1(蓝波段):0.45~0.52
    B2(绿波段):0.52~0.59
    B3(红波段):0.63~0.69
    B4(近红外波段)0.76~0.90
    B5(红边波段1):0.69~0.73
    B6(红边波段2):0.73~0.77
    B7(紫波段):0.40~0.45
    B8(黄波段):0.59~0.63
    空间分辨率Spatial resolution/m 16
    幅宽 Width/km 864.2
    下载: 导出CSV

    表 2  波段统计信息

    Table 2.  Statistics of spectral information

    波段
    Band
    中心波长
    Central
    wavelength
    最小值
    Min
    最大值
    Max
    均值
    Mean
    标准偏差
    Standard
    deviation
    信息熵
    Information entropy
    1 0.485 0 5 080 343 147 6.28
    2 0.555 0 5 918 605 246 7.26
    3 0.660 0 5 503 604 286 7.22
    4 0.830 0 6 061 1 837 1 064 9.23
    5 0.710 0 6 731 827 337 7.91
    6 0.750 0 5 871 1 482 849 8.98
    7 0.425 0 3 456 162 87 5.32
    8 0.610 0 4 971 641 237 7.23
    下载: 导出CSV

    表 3  最佳波段指数OIF

    Table 3.  The value of OIF

    序号
    Number
    波段组合
    Band
    combination
    OIF
    OIF
    序号
    Number
    波段组合
    Band
    combination
    OIF
    OIF
    序号
    Number
    波段组合
    Band
    combination
    OIF
    OIF
    序号
    Number
    波段组合
    Band
    combination
    OIF
    OIF
    1 4-6-7 924 15 4-7-8 611 29 2-3-6 530 43 3-5-7 264
    2 3-4-6 902 16 2-4-7 604 30 6-7-8 520 44 1-3-5 264
    3 4-5-6 850 17 3-5-6 592 31 2-6-7 513 45 1-2-5 263
    4 4-6-8 837 18 5-6-7 592 32 2-6-8 507 46 2-3-8 240
    5 1-4-6 834 19 1-3-4 587 33 1-3-6 503 47 1-5-7 223
    6 2-4-6 807 20 2-4-8 585 34 1-6-8 489 48 1-2-8 205
    7 4-5-7 685 21 1-4-8 570 35 1-2-6 481 49 1-2-3 202
    8 3-4-5 675 22 5-6-8 564 36 1-6-7 471 50 2-7-8 201
    9 4-5-8 642 23 1-2-4 562 37 2-3-5 298 51 1-3-8 200
    10 1-4-5 629 24 1-4-7 561 38 3-5-8 298 52 2-3-7 197
    11 2-4-5 627 25 2-5-6 550 39 2-5-8 294 53 3-7-8 194
    12 3-4-7 625 26 1-5-6 548 40 5-7-8 268 54 1-2-7 156
    13 3-4-8 617 27 3-6-8 535 41 2-5-7 267 55 1-7-8 155
    14 2-3-4 612 28 3-6-7 532 42 1-5-8 265 56 1-3-7 148
    下载: 导出CSV

    表 4  J-M距离

    Table 4.  Jeffreys-Matusita Distance

    序号
    Number
    波段
    组合
    Band
    combination
    林地-
    玉米
    Forest-Maize
    林地-
    水稻
    Forest-Paddy
    林地-沼泽
    化草甸
    Forest-Swamp meadow
    林地-
    草地
    Forest-Grassland
    玉米-
    水稻
    Maize-Paddy
    玉米-沼泽
    化草甸
    Maize-Swamp meadow
    玉米-
    草地
    Maize-Grassland
    水稻-沼泽
    化草甸
    Paddy-Swamp meadow
    水稻-
    草地
    Paddy-Grassland
    沼泽化草
    甸-草地
    Swamp meadow-Grassland
    1 3-4-5 2.00 2.00 1.89 2.00 2.00 2.00 2.00 1.88 2.00 1.05
    2 4-5-8 2.00 2.00 1.87 2.00 2.00 2.00 2.00 1.95 2.00 0.90
    3 1-4-5 2.00 2.00 1.99 2.00 2.00 1.99 2.00 1.96 1.99 1.75
    4 2-4-5 2.00 2.00 2.00 2.00 2.00 1.99 2.00 1.93 1.99 1.50
    5 3-4-8 2.00 2.00 1.83 2.00 1.98 1.88 2.00 1.95 1.99 1.11
    6 2-3-4 2.00 2.00 2.00 2.00 1.94 2.00 2.00 1.91 2.00 1.11
    7 1-3-4 2.00 2.00 1.99 2.00 1.99 1.89 1.99 1.96 1.99 1.74
    8 2-4-8 2.00 2.00 2.00 2.00 1.99 2.00 2.00 1.99 2.00 1.47
    9 1-4-8 2.00 2.00 1.98 2.00 2.00 1.88 2.00 1.99 1.98 1.76
    10 1-2-4 2.00 2.00 2.00 2.00 1.99 1.98 2.00 1.92 1.99 1.30
    下载: 导出CSV

    表 5  方法一排名

    Table 5.  Rank of method 1

    波段
    组合
    Band
    combination
    林地-
    玉米
    Forest-
    Maize
    林地-
    水稻
    Forest-
    Paddy
    林地-沼泽
    化草甸
    Forest-Swamp
    meadow
    林地-
    草地
    Forest-
    Grassland
    玉米-
    水稻
    Maize-
    Paddy
    玉米-沼泽
    化草甸
    Maize-Swamp
    meadow
    玉米-
    草地
    Maize-
    Grassland
    水稻-沼泽
    化草甸
    Paddy-Swamp
    meadow
    水稻-
    草地
    Paddy-
    Grassland
    沼泽化草
    甸-草地
    Swamp meadow-
    Grassland
    总分
    Total
    Score
    排名
    Rank
    1-4-5 1 1 2 1 1 2 1 2 2 2 15 1
    2-4-8 1 1 1 1 2 1 1 1 1 5 15 1
    1-4-8 1 1 3 1 1 5 1 1 3 1 18 2
    2-4-5 1 1 1 1 1 2 1 4 2 4 18 2
    1-3-4 1 1 2 1 2 4 2 2 2 3 20 3
    1-2-4 1 1 1 1 2 3 1 5 2 6 23 4
    2-3-4 1 1 1 1 4 1 1 6 1 7 24 5
    4-5-8 1 1 5 1 1 1 1 3 1 9 24 5
    3-4-5 1 1 4 1 1 1 1 7 1 8 26 6
    3-4-8 1 1 6 1 3 5 1 3 2 7 30 7
    下载: 导出CSV

    表 6  波段组合综合排名

    Table 6.  Comprehensive rank of band combinations

    波段组合
    Band
    combination
    方法二得分
    Score of
    method 2
    方法二排名
    Rank of
    method 2
    方法一排名
    Rank of
    method 1
    综合排名
    Overall
    1-4-5 19.67 1 1 1
    1-4-8 19.60 2 2 2
    2-4-8 19.43 4 1 3
    1-3-4 19.54 3 3 4
    2-4-5 19.40 5 2 5
    1-2-4 19.18 6 4 6
    2-3-4 18.96 7 5 7
    3-4-5 18.81 8 6 8
    4-5-8 18.71 10 5 9
    3-4-8 18.75 9 7 10
    下载: 导出CSV

    表 7  3种波段组合分类结果评价

    Table 7.  Classification accurary of 3 kinds of band combination

    类型
    Class
    波段组合 1-4-5
    Band combination 1-4-5
    波段组合 1-4-8
    Band combination 1-4-8
    波段组合 2-3-4
    Band combination 2-3-4
    制图精度
    Mapping accuracy
    用户精度
    User accuracy
    制图精度
    Mapping accuracy
    用户精度
    User accuracy
    制图精度
    Mapping accuracy
    用户精度
    User accuracy
    草地 Grassland 86.24 87.97 94.49 73.94 96.92 57.71
    水稻 Paddy 87.02 99.11 78.55 99.83 78.4 99.4
    林地 Forest 92.57 83.77 88.06 89.78 97 90.57
    玉米 Maize 93.35 84.31 92.86 70.57 93.1 69.63
    沼泽化草甸 Swamp meadow 97.61 77.69 87.23 74.89 61.11 87.37
    总体精度 Overall accuracy/% 89.15 83.75 79.07
    Kappa 系数 Kappa coefficient 0.846 3 0.754 8 0.714 9
    下载: 导出CSV
  • [1]

    Mao Dehua, Luo Ling, Wang Zongming, et al. Conversions between natural wetlands and farmland in China: A multiscale geospatial analysis[J]. The Science of the Total Environment, 2018, 634: 550-560. doi: 10.1016/j.scitotenv.2018.04.009
    [2] 张怀清, 郑 曼, 蒋 娴. 基于ArcGIS Engine的退耕还林决策支持系统研究[J]. 林业科学研究, 2008, 21(Z1):65-68. doi: 10.3321/j.issn:1001-1498.2008.z1.013

    [3] 刘文超, 刘纪远, 匡文慧. 陕北地区退耕还林还草工程土壤保护效应的时空特征[J]. 地理学报, 2019, 74(9):1835-1852. doi: 10.11821/dlxb201909010

    [4]

    Zhu Bingmiao, Zhu Xun, Zhang Ran, et al. Study of multiple land use planning based on the coordinated development of wetland farmland: A case study of Fuyuan City, China[J]. Sustainability, 2019, 11(1): 271. doi: 10.3390/su11010271
    [5] 黄 妮, 刘殿伟, 王宗明, 等. 基于GIS的三江平原退耕还湿空间决策分析[J]. 地理科学, 2009, 29(6):874-879.

    [6] 王雅慧, 陈尔学, 郭 颖, 等. 高分辨率多光谱遥感影像森林类型分类深度U-net优化方法[J]. 林业科学研究, 2020, 33(1):11-18.

    [7] 姚 博, 张怀清, 刘 洋, 等. 面向对象CART决策树方法的湿地遥感分类[J]. 林业科学研究, 2019, 32(5):91-98.

    [8] 戴昌达, 雷莉萍. TM图像的光谱信息特征与最佳波段组合[J]. 环境遥感, 1989(4):282-292.

    [9]

    Wang Chunxing, Ma Nan, Ming Yanfang, et al. Classification of hyperspectral imagery with a 3D convolutional neural network and J-M distance[J]. Advances in Space Research, 2019, 64(4): 886-899. doi: 10.1016/j.asr.2019.05.005
    [10] 郭力娜, 李 帅, 张梦华, 等. 最佳波段组合的城市土地利用类型提取[J]. 测绘科学, 2019, 44(8):161-167.

    [11] 王平平, 田淑芳. 基于 WorldView2 数据的岩性信息增强方法——以内蒙古扎嘎乌苏地区为例[J]. 国土资源遥感, 2016, 28(4):176-184.

    [12] 和晓风. 基于GF-1遥感湿地类型提取研究[D]. 长沙: 中南林业科技大学, 2017.

    [13]

    Xie Fuding, Li Fangfei, Lei Cunkuan, et al. Unsupervised band selection based on artificial bee colony algorithm for hyperspectral image classification[J]. Applied Soft Computing, 2019, 75: 428-440. doi: 10.1016/j.asoc.2018.11.014
    [14]

    Wang Mingwei, Wu Chunming, Wang Lizhe, et al. A feature selection approach for hyperspectral image based on modified ant lion optimizer[J]. Knowledge-Based Systems, 2019, 168: 39-48. doi: 10.1016/j.knosys.2018.12.031
    [15] 魏芳洁. 高光谱图像波段选择方法的研究[D]. 哈尔滨: 哈尔滨工程大学, 2013.

    [16] 马 娜, 胡云锋, 庄大方, 等. 基于最佳波段指数和J-M距离可分性的高光谱数据最佳波段组合选取研究——以环境小卫星高光谱数据在东莞市的应用为例[J]. 遥感技术与应用, 2010, 25(3):358-365. doi: 10.11873/j.issn.1004-0323.2010.3.358

    [17]

    Frampton W J, Dash J, Watmough G, et al. Evaluating the capabilities of Sentinel-2 for quantitative estimation of biophysical variables in vegetation[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2013, 82: 83-92. doi: 10.1016/j.isprsjprs.2013.04.007
    [18]

    Adelabu S, Mutanga O, Adam E. Evaluating the impact of red-edge band from Rapideye image for classifying insect defoliation levels[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2014, 95: 34-41. doi: 10.1016/j.isprsjprs.2014.05.013
    [19] 梁志国, 隋 傲, 于 颖, 等. 高分六号遥感数据在有林地识别中的应用[J]. 东北林业大学学报, 2020, 48(5):35-39. doi: 10.3969/j.issn.1000-5382.2020.05.007

    [20]

    Peng L, Liu K, Cao J, et al. Combining GF-2 and RapidEye satellite data for mapping mangrove species using ensemble machine-learning methods[J]. International Journal of Remote Sensing, 2020, 41(3): 813-838. doi: 10.1080/01431161.2019.1648907
    [21]

    Madonsela S, Cho M A, Mathieu R, et al. Multi-phenology WorldView-2 imagery improves remote sensing of savannah tree species[J]. International Journal of Applied Earth Observation and Geoinformation, 2017, 58: 65-73. doi: 10.1016/j.jag.2017.01.018
    [22] 张沁雨, 李 哲, 夏朝宗, 等. 高分六号遥感卫星新增波段下的树种分类精度分析[J]. 地球信息科学学报, 2019, 21(10):1619-1628. doi: 10.12082/dqxxkx.2019.190116

    [23] 魏梦凡. 基于Sentinel-2A卫星遥感影像的开封市冬小麦种植面积提取技术研究[D]. 开封: 河南大学, 2019.

    [24] 梁 继, 郑镇炜, 夏诗婷, 等. 高分六号红边特征的农作物识别与评估[J]. 遥感学报, 2020, 24(10):1168-1179.

  • [1] 何友均崔国发邹大林郑杰董建生李永波 . 三江源自然保护区主要森林群落物种多样性研究. 林业科学研究, 2007, 20(2): 241-245.
    [2] 胡鸿鞠洪波田昕杨雪清孙志超 . 基于GF-1和TM数据的金河林业局天然林保护成效评估研究. 林业科学研究, 2019, 32(1): 141-146. doi: 10.13275/j.cnki.lykxyj.2019.01.019
    [3] 杨海龙李迪强朵海瑞马 剑 . 梵净山国家级自然保护区植被分布与黔金丝猴生境选择. 林业科学研究, 2010, 23(3): 393-398.
    [4] 武晓莉吕相海李国会张友焱党宏忠周泽福 . 晋西黄土区典型植被类型的土壤水分特征. 林业科学研究, 2014, 27(4): 487-492.
    [5] 仲崇禄 . 木麻黄苗期最佳固氮基因型组合体研究. 林业科学研究, 1993, 6(6): 654-660.
    [6] 陈永富刘华邹文涛张怀清 . 三江源湿地变化驱动因子定量研究. 林业科学研究, 2012, 25(5): 545-550.
    [7] 凌成星鞠洪波张怀清孙华 . 基于Worldview-2八波段影像改进指数的湿地类型分类研究. 林业科学研究, 2014, 27(5): 639-643.
    [8] 白黎娜李增元Fabio MaselliLorenzo BottaiAlberto OrtolaniStefano Romanelli . 应用遥感数据识别意大利沿海松林灾害级别. 林业科学研究, 2001, 14(5): 479-483.
    [9] 张时敏李广武 . 山西庞泉沟自然保护区三种落叶松球果花蝇发生与危害研究*. 林业科学研究, 1994, 7(4): 449-452.
    [10] 胡博鞠洪波刘华郝泷刘海 . 基于证据理论组合多分类规则实现大区域植被遥感分类研究*. 林业科学研究, 2017, 30(2): 194-199. doi: 10.13275/j.cnki.lykxyj.2017.02.002
    [11] 纪平易浩若白黎娜张彦忠 . 一种识别NOAA/AVHRR数据异常高温点方法的研究*. 林业科学研究, 1993, 6(6): 699-702.
    [12] 张广才黄利江于卫平马晖王涵 . 毛乌素沙地南缘植被类型及其演替规律初步研究. 林业科学研究, 2004, 17(s1): 131-136.
    [13] 彭佳龙史小华张汝忠 . 浙江括苍山自然保护区种子植物区系分析. 林业科学研究, 2008, 21(5): 713-718.
    [14] . 山西太宽河自然保护区植物区系研究. 林业科学研究, 2009, 22(3): -.
    [15] . 江西官山自然保护区昆虫区系分析. 林业科学研究, 2009, 22(3): -.
    [16] 宋扬徐正会李春良张宁张力蒋华莫福燕 . 云南南滚河自然保护区蚁科昆虫区系分析. 林业科学研究, 2013, 26(6): 773-780.
    [17] 吴兴亮李泰辉章卫民钟琼芯宋斌郭建荣 . 海南岛坝王岭自然保护区的多孔菌研究. 林业科学研究, 1998, 11(2): 163-168.
    [18] 徐正会蒋兴成陈志强吴定敏 . 高黎贡山自然保护区东坡垂直带蚂蚁群落研究. 林业科学研究, 2001, 14(2): 115-124.
    [19] 徐正会吴定敏陈志强蒋兴成 . 高黎贡山自然保护区东坡水平带蚂蚁群落研究. 林业科学研究, 2001, 14(6): 603-609.
    [20] 李文英贺运春王建明张作刚宋东辉 . 庞泉沟自然保护区虫生真菌资源调查. 林业科学研究, 2003, 16(3): 339-342.
  • 加载中
图(4) / 表(7)
计量
  • 文章访问数:  3637
  • HTML全文浏览量:  2075
  • PDF下载量:  39
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-03-15
  • 录用日期:  2021-06-03
  • 网络出版日期:  2021-11-19
  • 刊出日期:  2022-02-20

面向三江保护区植被类型识别的最佳波段组合研究

    通讯作者: 赵峰, zhaofeng@ifrit.ac.cn
  • 1. 中国林业科学研究院资源信息研究所, 国家林业和草原局林业遥感与信息技术重点实验室, 北京 100091

摘要:  目的 针对高分六号影像数据面向植被类型、特别是湿地植被信息识别的最佳波段组合进行研究,为高分六号数据在湿地资源监测中的进一步应用提供参考。 方法 以三江国家级自然保护区为对象,基于高分六号影像数据,结合当地植被特点,开展保护区植被类型识别中的最佳波段组合研究。从影像信息特征、最佳指数因子以及光谱特征曲线和地物可分性角度进行综合分析,获取基于GF6/WFV数据的三江保护区植被类型识别的最佳波段组合。同时,通过对比实验,利用支持向量机分类方法,对样区的植被类型进行信息提取,证实了结果的可靠性。 结果 波段组合4(R)-5(G)-1(B)识别植被类型精度最高。总体精度为89.15%,Kappa系数为0.8463。 结论 波段组合4(R)-5(G)-1(B)为三江保护区植被类型识别的最佳波段组合。本研究可为东北地区湿地保护区的退耕还湿工程监测提供参考。

English Abstract

  • 随着遥感技术的发展和广泛应用,其在退耕还林(草)[1-3]、退耕还湿[4-5]等生态恢复工程的监测中也发挥了重要作用。除了工程进展和实施效果外,植被类型识别也是工程监测的主要内容。现有的植被类型识别方法多采用不同分类器进行监督分类[6-7],其中样本选择是监督分类的必要步骤。通过研究最佳波段组合方式,增强影像信息的可判读性,从而实现分类样本的快速获取。

    最佳波段组合一直是卫星数据研究的热点问题。对于拥有多波段光谱信息的传感器来说,应用最佳波段组合可细致分析卫星原始数据各波段信息特征以及在各个领域的应用价值[8-9]。如郭力娜等[10]采用最佳指数因子与不同土地利用类型光谱曲线相结合的方法,选取了适用于Landsat8 OLI在土地利用类型提取的最佳波段组合为1-5-7;王平平等[11]在研究WorldView2原始波段最佳波段组合的基础上,增加了不同光谱指数,增强岩性信息提取;和晓风等[12]针对湿地类型,从信息量、信息熵、最佳指数3个定量指标出发,确定了GF1影像最佳波段组合为4-3-2。此外,最佳波段组合也是一种数据选择方法,许多学者采用最佳波段组合方式来减少冗余信息与噪声[13-15],进一步提高分类精度。

    高分六号卫星是我国首颗精准农业观测卫星,其新增的红边波段能够有效反映植被的光谱信息,可为退耕还湿工程监测,特别是在区分湿地保护区内农作物与湿地植被方面提供有力的数据支持。但目前,利用高分六号数据针对湿地保护区植被信息识别上的研究较少。本研究以三江国家级自然保护区为对象,基于高分六号数据,从制图中的视觉评估角度,开展保护区植被类型识别中的最佳波段组合研究,并对结果进行定量评估。以期为高分六号数据在湿地资源监测中的进一步应用提供参考。

    • 三江国家级自然保护区位于黑龙江省佳木斯市抚远市和同江市境内,面积198 089 hm2,是以沼泽湿地为主要保护对象的自然保护区。2000年被国务院批准为国家级自然保护区,2002年被列入《国际重要湿地名录》。保护区内有丰富的湿地资源,土壤肥沃、湿地植被茂密且类型复杂多样。由于历史上多次大规模人口迁入及开荒,原有的大部分湿地被农田侵占,原始湿地资源遭到破坏。为恢复保护区内湿地生态系统功能,保护野生动物与鸟类正常繁衍栖息,自2016年起,保护区开始实施“退耕还湿”工程。为方便区划和管理,当地政府将三江自然保护区分为黑龙江流域保护区(面积106 552 hm2);乌苏里江流域保护区(面积91 537 hm2)。本研究选择黑龙江流域保护区作为研究区域(图1)。

      图  1  研究区高分六号影像(RGB=3,2,1)

      Figure 1.  Image of GF6 in the study area(RGB=3,2,1)

    • “高分六号”卫星搭载了8谱段CMOS探测器,首次增加了能够有效反映植被特有光谱特性的“红边”波段(具体参数见表1)。本实验所用数据下载于陆地观测卫星数据(http://36.112.130.153:7777/#/home),2019年9月12号获取,数据经过辐射定标、大气校正以及正射校正最终形成后续研究的遥感数据。此外,样本数据为2020年8月底到9月初实地采集,由于研究区部分地区自然条件复杂,导致某些区域难以进入,故采用无人机进行数据采集作为补充。同时,结合实地考察经验与google earth影像选取足够数量的样本点,作为后续实验的基础数据。

      表 1  高分六号WFV数据

      Table 1.  GF6/WFV

      项目 Project高分六号WFV数据GF6/WFV
      光谱范围 Wavelength range/ um B1(蓝波段):0.45~0.52
      B2(绿波段):0.52~0.59
      B3(红波段):0.63~0.69
      B4(近红外波段)0.76~0.90
      B5(红边波段1):0.69~0.73
      B6(红边波段2):0.73~0.77
      B7(紫波段):0.40~0.45
      B8(黄波段):0.59~0.63
      空间分辨率Spatial resolution/m 16
      幅宽 Width/km 864.2
    • 从影像信息特征与不同植被类型间可分性两个方面出发,选取标准偏差、信息熵、相关系数以及衍生的波段指数来反映影像信息特征,从整体角度缩小最佳波段组合选取范围;再针对保护区植被类型的光谱特征,利用J-M距离方法计算不同植被类型间的可分性,最终确定研究区植被类型识别的最佳波段组合。

    • 标准偏差和信息熵是衡量波段信息量的常用方法。标准偏差可以反映图像像元值与总体均值之间的离散程度,而信息熵表达的是图像灰度分布的聚集特征所包含的信息量。标准偏差与信息熵越大,则信息量就越丰富。统计研究区WFV传感器8个波段的标准偏差、信息熵,最大最小值、均值结果如表2所示。

      表 2  波段统计信息

      Table 2.  Statistics of spectral information

      波段
      Band
      中心波长
      Central
      wavelength
      最小值
      Min
      最大值
      Max
      均值
      Mean
      标准偏差
      Standard
      deviation
      信息熵
      Information entropy
      1 0.485 0 5 080 343 147 6.28
      2 0.555 0 5 918 605 246 7.26
      3 0.660 0 5 503 604 286 7.22
      4 0.830 0 6 061 1 837 1 064 9.23
      5 0.710 0 6 731 827 337 7.91
      6 0.750 0 5 871 1 482 849 8.98
      7 0.425 0 3 456 162 87 5.32
      8 0.610 0 4 971 641 237 7.23
    • 相关系数是用以反映变量之间相关关系密切程度的统计指标。对于遥感影像来说,相关系数表达的是波段间的信息重合度。相关系数越大,则两个波段间相关性越大,二者间的冗余度越大,信息重复度越高。统计研究区各波段间的相关系数并绘制相关系数矩阵图(图2),图中由蓝到红相关性依次增加。考虑到GF6/WFV数据原始波段设置并未按照波长长度排列,不利于相关性分析,故图2坐标按照中心波长长度重新排序,这样更能表达GF6/WFV数据的各波段之间的相关性。

      图  2  相关系数矩阵

      Figure 2.  Image of correlation coefficient matrix

    • 本研究选择应用较为广泛的OIF指数[16]方法。OIF指数不以单个波段为目标寻找最佳波段组合,而是综合了波段组合内部方差与相关系数,从组合角度评价图像信息量。OIF指数越大,则波段组合信息量越大。其计算公式如下:

      $ OIF = \sum\limits_{i = 1}^3 {{S_i}/} \sum\limits_{j = 1}^3 {\left| {{r_{ij}}} \right|} $

      (1)

      其中: Si 为波段i的标准差,rij为波段i与波段j之间的相关系数。

      从GF6/WFV数据的8个波段中任意选取3个波段组合,共计56种组合方式,计算各波段组合的OIF指数,并按照OIF数值从大到小排序,结果如表3所示。

      表 3  最佳波段指数OIF

      Table 3.  The value of OIF

      序号
      Number
      波段组合
      Band
      combination
      OIF
      OIF
      序号
      Number
      波段组合
      Band
      combination
      OIF
      OIF
      序号
      Number
      波段组合
      Band
      combination
      OIF
      OIF
      序号
      Number
      波段组合
      Band
      combination
      OIF
      OIF
      1 4-6-7 924 15 4-7-8 611 29 2-3-6 530 43 3-5-7 264
      2 3-4-6 902 16 2-4-7 604 30 6-7-8 520 44 1-3-5 264
      3 4-5-6 850 17 3-5-6 592 31 2-6-7 513 45 1-2-5 263
      4 4-6-8 837 18 5-6-7 592 32 2-6-8 507 46 2-3-8 240
      5 1-4-6 834 19 1-3-4 587 33 1-3-6 503 47 1-5-7 223
      6 2-4-6 807 20 2-4-8 585 34 1-6-8 489 48 1-2-8 205
      7 4-5-7 685 21 1-4-8 570 35 1-2-6 481 49 1-2-3 202
      8 3-4-5 675 22 5-6-8 564 36 1-6-7 471 50 2-7-8 201
      9 4-5-8 642 23 1-2-4 562 37 2-3-5 298 51 1-3-8 200
      10 1-4-5 629 24 1-4-7 561 38 3-5-8 298 52 2-3-7 197
      11 2-4-5 627 25 2-5-6 550 39 2-5-8 294 53 3-7-8 194
      12 3-4-7 625 26 1-5-6 548 40 5-7-8 268 54 1-2-7 156
      13 3-4-8 617 27 3-6-8 535 41 2-5-7 267 55 1-7-8 155
      14 2-3-4 612 28 3-6-7 532 42 1-5-8 265 56 1-3-7 148
    • 表2可知,从不同波段的中心波长来说,按照波长从短到长依次是b7>b1>b2>b8>b3>b5>b6>b4,标准偏差从小到大顺序为b7>b1>b8>b2>b3>b5>b6>b4,信息熵差从小到大顺序为b7>b1>b3>b8>b2>b5>b6>b4,中心波长长度几乎与信息量呈正相关。信息量总趋势是随中心波长增长而增大,即波长越长,图像反差越大,所含信息量越丰富。其次,由图2可知,波段间相关系数出现分块现象,即:b1、b2、b3、b5、b7、b8(简称:一组)以及b4、b6(简称:二组)。这两组波段组间相关性相对较低,信息重叠程度小,但组内部各波段相关性较高,信息重叠度较大。在最佳波段选择原则中我们认为,最终结果应尽可能选择波段间相关性小的3个波段。因此,最佳波段组合不宜以组内部波段组成。

      表3可知,OIF指数排名前6名波段组合均存在b4与b6。通过前文信息特征分析可知,b4和b6所含信息量虽大,但二者的相关性极高,不宜同时选取。考虑到波段组合的独立性以及信息量,研究决定从所有段组合中选择含有b4或b6(不同时存在)且同时OIF数值排名位于前20的波段组合作为选择范围,用于进一步分析,即OIF指数排名第7~26的波段组合。

    • 研究区植被类型丰富,是三江平原原始沼泽的核心和缩影。但由于历史上大规模开荒活动,三江国家级自然保护区内存在大量耕地,占地面积较大,作物类型单一。研究区耕地类型以玉米和水稻为主,沼泽湿地以沼泽化草甸为主。根据研究区植被类型特点,本研究选择沼泽化草甸、玉米、水稻、林地以及草地5种植被类型作为研究对象,以前期外业调查数据和谷歌影像作为参考,每种植被类型选择30个样本,计算每种植被类型在8个波段上的响应值,并绘制不同植被类型的光谱特征曲线,如图3所示。

      图  3  5种植被类型的光谱特征曲线

      Figure 3.  5 Spectral characteristic curve of planting quilt type

      根据图3可知,5种植被类型的光谱响应值的总体走势相近,在单波段区分程度各不相同。其中,红边波段b5、b6在植被类型识别中表现优异,尤其是波段5对各植被类型区分的综合能力位于前列;b4和b6波段的玉米、水稻比其他植被类型辐射值更高,有利于区分耕地与其他植被类型;可见光波段(b1、b2、b3、b7、b8)的林地辐射值最低而水稻最高,但是b7波段对植被类型区分能力极差,为进一步缩小最佳波段选择范围,故最佳波段组合中不考虑第7波段。此外,综合考虑到波段间信息量,故最佳波段组合中应包含b4,不包含b6与b7。

    • J-M距离是度量不同类别可分性较为有效的方法,公式如下:

      $ {J_{ij}} = {\left\{ {{{\int\limits_x {\left[ {\sqrt {p\left( {{X \mathord{\left/ {\vphantom {X {{\omega _i}}}} \right. } {{\omega _i}}}} \right)} - \sqrt {p\left( {{X \mathord{\left/ {\vphantom {X {{\omega _j}}}} \right. } {{\omega _j}}}} \right)} } \right]} }^2}dX} \right\}^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} $

      (2)

      式中: p( X/ωi)为条件概率密度,即第i个像元属于第ωi个类别的几率。Jij取值在0~2之间,其中,Jij位于0. 0~1. 0之间时,则样本间不具有光谱可分性;Jij位于1.0~1.9之间时,样本间具有一定的光谱可分性,但存在较大程度的重叠;Jij位于1.9~2.0之间时,样本间具有很好的光谱可分性。

      通过前期基于OIF指数的信息量以及植被类型光谱特征的综合分析,从OIF指数法筛选出的20组波段组合中,排除包含波段6、7的波段组合,共选出10组波段组合,依次计算不同地物类型间的J-M距离(结果保留2位小数)。结果如表4所示。

      表 4  J-M距离

      Table 4.  Jeffreys-Matusita Distance

      序号
      Number
      波段
      组合
      Band
      combination
      林地-
      玉米
      Forest-Maize
      林地-
      水稻
      Forest-Paddy
      林地-沼泽
      化草甸
      Forest-Swamp meadow
      林地-
      草地
      Forest-Grassland
      玉米-
      水稻
      Maize-Paddy
      玉米-沼泽
      化草甸
      Maize-Swamp meadow
      玉米-
      草地
      Maize-Grassland
      水稻-沼泽
      化草甸
      Paddy-Swamp meadow
      水稻-
      草地
      Paddy-Grassland
      沼泽化草
      甸-草地
      Swamp meadow-Grassland
      1 3-4-5 2.00 2.00 1.89 2.00 2.00 2.00 2.00 1.88 2.00 1.05
      2 4-5-8 2.00 2.00 1.87 2.00 2.00 2.00 2.00 1.95 2.00 0.90
      3 1-4-5 2.00 2.00 1.99 2.00 2.00 1.99 2.00 1.96 1.99 1.75
      4 2-4-5 2.00 2.00 2.00 2.00 2.00 1.99 2.00 1.93 1.99 1.50
      5 3-4-8 2.00 2.00 1.83 2.00 1.98 1.88 2.00 1.95 1.99 1.11
      6 2-3-4 2.00 2.00 2.00 2.00 1.94 2.00 2.00 1.91 2.00 1.11
      7 1-3-4 2.00 2.00 1.99 2.00 1.99 1.89 1.99 1.96 1.99 1.74
      8 2-4-8 2.00 2.00 2.00 2.00 1.99 2.00 2.00 1.99 2.00 1.47
      9 1-4-8 2.00 2.00 1.98 2.00 2.00 1.88 2.00 1.99 1.98 1.76
      10 1-2-4 2.00 2.00 2.00 2.00 1.99 1.98 2.00 1.92 1.99 1.30

      基于J-M距离评估植被类型的可分性结果可知,从56个波段组合中初选出的10个波段组合在植被类型识别方面均表现优异,仅在区分沼泽化草甸与草地效果不佳,但在其他植被类型间J-M距离均大于1.83,且绝大多数J-M距离大于1.9,验证了筛选结果的可靠性。其中,在不同波段组合中林地、水稻、玉米和草地4种植被类型间的J-M距离均大于1.94,表明该10个组合在这4种植被类型具有较强区分能力。此外,沼泽化草甸与其他植被类型间J-M距离相对较低,不宜区分,尤其是沼泽化草甸与草地,J-M距离最大为1.76(波段组合1-4-8),最小仅为0.90,表明沼泽化草甸与草地的光谱信息十分相近,容易混淆不易区分。

    • 由于J-M距离方法仅能衡量两种地物类型间的可分性,本研究采用两种方法相结合的综合评估法,进一步确定三江保护区植被类型识别的最佳波段组合。方法一:将各植被类型间的J-M距离按大小排序,累加排名次序,统计总排名作为最终分数,总分越小则说明该波段综合排名越高,植被类型识别能力越强(表5)。方法二:累加不同植被类型两两间J-M距离,记为总分,并按总分大小排序,总分越高则表明波段组合综合识别能力越强。综合两种方法结果如表6所示。

      表 5  方法一排名

      Table 5.  Rank of method 1

      波段
      组合
      Band
      combination
      林地-
      玉米
      Forest-
      Maize
      林地-
      水稻
      Forest-
      Paddy
      林地-沼泽
      化草甸
      Forest-Swamp
      meadow
      林地-
      草地
      Forest-
      Grassland
      玉米-
      水稻
      Maize-
      Paddy
      玉米-沼泽
      化草甸
      Maize-Swamp
      meadow
      玉米-
      草地
      Maize-
      Grassland
      水稻-沼泽
      化草甸
      Paddy-Swamp
      meadow
      水稻-
      草地
      Paddy-
      Grassland
      沼泽化草
      甸-草地
      Swamp meadow-
      Grassland
      总分
      Total
      Score
      排名
      Rank
      1-4-5 1 1 2 1 1 2 1 2 2 2 15 1
      2-4-8 1 1 1 1 2 1 1 1 1 5 15 1
      1-4-8 1 1 3 1 1 5 1 1 3 1 18 2
      2-4-5 1 1 1 1 1 2 1 4 2 4 18 2
      1-3-4 1 1 2 1 2 4 2 2 2 3 20 3
      1-2-4 1 1 1 1 2 3 1 5 2 6 23 4
      2-3-4 1 1 1 1 4 1 1 6 1 7 24 5
      4-5-8 1 1 5 1 1 1 1 3 1 9 24 5
      3-4-5 1 1 4 1 1 1 1 7 1 8 26 6
      3-4-8 1 1 6 1 3 5 1 3 2 7 30 7

      表 6  波段组合综合排名

      Table 6.  Comprehensive rank of band combinations

      波段组合
      Band
      combination
      方法二得分
      Score of
      method 2
      方法二排名
      Rank of
      method 2
      方法一排名
      Rank of
      method 1
      综合排名
      Overall
      1-4-5 19.67 1 1 1
      1-4-8 19.60 2 2 2
      2-4-8 19.43 4 1 3
      1-3-4 19.54 3 3 4
      2-4-5 19.40 5 2 5
      1-2-4 19.18 6 4 6
      2-3-4 18.96 7 5 7
      3-4-5 18.81 8 6 8
      4-5-8 18.71 10 5 9
      3-4-8 18.75 9 7 10

      综合两种方法进行评估的优势在于,如果仅采用方法一进行评估,会产生拉平效应。例如表4中的沼泽化草甸与草地,波段组合1-3-4的J-M距离为1.74,波段组合1-4-5为1.75,数值相差0.01,排名相差1。同样在沼泽化草甸与草地中,波段组合2-4-5的J-M距离为1.50,波段组合1-3-4为1.74,数值相差0.24,但排名仍相差1。如果仅采用方法二进行评估,则又会出现因为某两种植被类型的J-M距离数值过小和过大而过分影响该波段的综合排名现象。故在两种评估方法中同时表现优异的波段组合,才应是我们优选的植被类型识别最佳波段组合。

      根据表6可知,综合排名第一的波段组合为1-4-5。本研究根据最佳波段选择原则,在研究区地物光谱特征分析基础上,结合OIF值排序以及地物类型最大差异的分析,最终确定1-4-5波段为研究区植被类型识别的最佳波段组合。

    • 为了更加科学地判定波段1-4-5为最佳波段组合,本研究选取综合排名第二的波段组合1-4-8作为对比数据,以验证最佳波段组合的优越性;并将常用于植被识别的波段组合2-3-4作为另一组对比数据。采用在小样本条件下表现较好支持向量机方法对3种波段组合分别进行监督分类,获取研究区植被类型分布图。分类结果如图4所示。分类精度验证结果如表7所示。

      图  4  波段组合分类结果

      Figure 4.  Classification results of band combination

      表 7  3种波段组合分类结果评价

      Table 7.  Classification accurary of 3 kinds of band combination

      类型
      Class
      波段组合 1-4-5
      Band combination 1-4-5
      波段组合 1-4-8
      Band combination 1-4-8
      波段组合 2-3-4
      Band combination 2-3-4
      制图精度
      Mapping accuracy
      用户精度
      User accuracy
      制图精度
      Mapping accuracy
      用户精度
      User accuracy
      制图精度
      Mapping accuracy
      用户精度
      User accuracy
      草地 Grassland 86.24 87.97 94.49 73.94 96.92 57.71
      水稻 Paddy 87.02 99.11 78.55 99.83 78.4 99.4
      林地 Forest 92.57 83.77 88.06 89.78 97 90.57
      玉米 Maize 93.35 84.31 92.86 70.57 93.1 69.63
      沼泽化草甸 Swamp meadow 97.61 77.69 87.23 74.89 61.11 87.37
      总体精度 Overall accuracy/% 89.15 83.75 79.07
      Kappa 系数 Kappa coefficient 0.846 3 0.754 8 0.714 9

      表7可知,波段组合1-4-5在各植被类型识别精度均表现优异,总体精度和Kappa系数最高,故进一步证明波段波段组合1-4-5为本研究的最佳波段组合。此外,适合的RGB组合方式也是选择最佳波段组合的关键步骤,通过对比发现,R-4、G-5、B-1为RGB最佳组合色彩显示(图4)。

    • 基于GF6/WFV影像,从影像特征和类间可分性角度出发,对标准偏差、信息熵和相关系数等指标进行分析,通过最佳指数因子方法确定各个波段组合的信息量;再结合典型植被类型光谱特征曲线、J-M距离等确定了研究区典型植被类型识别的最佳波段组合,并得到如下结论:

      (1)对于GF6/WFV数据,从可见光到近红外各波段信息量总体为上升趋势,其中b4波段信息熵与标准偏差最大,b7最小,说明b4包含信息量最多而b7最少;依据波段间相关性可分为两组:b1、b2、b3、b7、b8(一组)和b4、b6(二组)。两组波段组内相关性高,信息冗余度高,组间相关性低,冗余程度低。对于影像特征分析可知,在后续研究中,只有从两组中分别选取适当波段进行组合,才能保证最终波段组合的独立性。

      (2)单一的波段选择方法具有一定局限性,综合方法可以实现方法间优势互补。OIF指数可评估波段组合整体信息含量,初步筛选信息丰富的波段组合;光谱特征曲线从单波段角度评估各类型的可分性;J-M距离从组合角度评价不同波段组合对植被类型的区分能力。综合以上3个角度可获得最终的最佳波段组合。

      (3)经过多指标的综合评估,在56个波段组合中,波段组合4(R)-5(G)-1(B)是三江保护区植被类型识别的最佳波段组合。其中林地、水稻、玉米、草地4种植被类型类间具有较强的可分性,但沼泽化草甸与其他植被类型间J-M距离相对较低,尤其与草地间存在部分混淆情况。

      本研究基于GF6/WFV影像对识别三江保护区植被类型的最佳波段组合进行了研究,研究结果中最佳波段组合1-4-5中包含近红与红边两个波段,前人研究表明,该两个波段在植被生物量[17]、健康评估[18]和树种[19-22]、农作物识别[23-24]中发挥了巨大的作用,有助于提高植被信息的提取精度。故本研究提出的最佳波段组合在植被相关的目视解译和分类等方面具有可靠性,对于三江平原的植被类型快速识别具有参考价值。GF6号卫星作为我国自主研发的传感器,其新增加的波段包含紫波段、黄波段以及两个红边波段4个特殊波段,本研究从原始波段组合角度分析了植被类型识别方面的最佳波段组合,在今后的研究中可以从指数进一步挖掘新增波段在植被方面的应用。

参考文献 (24)

目录

    /

    返回文章
    返回