留言板

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

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

面向林地分类的GF-2影像融合算法评价

胡曼 彭道黎

胡曼, 彭道黎. 面向林地分类的GF-2影像融合算法评价[J]. 浙江农林大学学报, 2017, 34(2): 340-348. doi: 10.11833/j.issn.2095-0756.2017.02.019
引用本文: 胡曼, 彭道黎. 面向林地分类的GF-2影像融合算法评价[J]. 浙江农林大学学报, 2017, 34(2): 340-348. doi: 10.11833/j.issn.2095-0756.2017.02.019
HU Man, PENG Daoli. Comparison of fusion algorithms for GF-2 data from extracted forestland information[J]. Journal of Zhejiang A&F University, 2017, 34(2): 340-348. doi: 10.11833/j.issn.2095-0756.2017.02.019
Citation: HU Man, PENG Daoli. Comparison of fusion algorithms for GF-2 data from extracted forestland information[J]. Journal of Zhejiang A&F University, 2017, 34(2): 340-348. doi: 10.11833/j.issn.2095-0756.2017.02.019

面向林地分类的GF-2影像融合算法评价

doi: 10.11833/j.issn.2095-0756.2017.02.019
基金项目: 

国家重点林业工程监测技术示范推广项目 2015-02

详细信息
    作者简介: 胡曼, 从事森林资源监测与评价的研究工作。E-mail:humanamanda@163.com
    通信作者: 彭道黎, 教授, 从事森林资源监测、森林经营管理以及林业信息领域研究。E-mail:dlpeng@bjfu.edu.cn
  • 中图分类号: S757.2

Comparison of fusion algorithms for GF-2 data from extracted forestland information

  • 摘要: 针对林业部门目前常用的遥感影像融合算法,探究适合于林地信息提取的GF-2影像融合算法。以GF-2卫星1 m全色/4 m多光谱分辨率平面影像为数据源,采用HSV变换(hue-saturation-value,颜色空间变换),Brovey变换(彩色标准化变换),PC变换(principle components,主成分变换),HPF变换(high-pass fusion,高通滤波变换),GS变换(gram-schmidt,正交化变换),Pansharp变换(超分辨率贝叶斯变换)等6种常用融合算法,通过目视和定量特征分析对融合效果进行评价,并结合面向对象分类方法对融合后影像进行地类信息提取和分析,探讨6种融合算法对GF-2影像在林区地类提取的适宜性。研究结果表明:基于Brovey和HSV算法的融合结果目视效果良好,清晰度与纹理增强明显;这2种融合算法影像在不同地类信息的提取上各有优势,HSV算法融合结果在不同地类的提取上效果最好,分类总精度可达85.14%,Brovey算法融合结果则在森林类型的提取上具有最高的分类总精度,为75.72%;其余4种融合算法在图像质量及其他地类提取中各有优势,具体融合算法的选取需根据应用目的和影区应用区域的实际情况而定。该研究可为林业部门提高GF-2卫星的适用性及大规模应用提供参考。
  • 图  1  融合影像目视对比图

    Figure  1  Visual contrast on fusion images

    图  2  熵值指标统计图

    Figure  2  Statistical chart of entropy indexes

    图  3  二阶矩指标统计图

    Figure  3  Statistical chart of second moment indexes

    图  4  融合影像分类结果图

    Figure  4  Classification results of fusion images

    表  1  原多光谱影像各波段与融合影像整体影像特征对比表

    Table  1.   Contrast on indexes of each band between multi-spectral image and fusion images

    影像/方法 波段 均值 平均梯度 相关系数 高频信息融人度
    原始影像 R 1 198.725 3.124
    G 1 295.094 3.679
    B 1 203.466 4.255
    NIR 2 006.270 6.043
    PAN 305.893 12.847
    彩色标准变换 (Brovey) G 106.712 13.285 0.105 0.328
    B 86.911 9.939 0.215 0.298
    NIR 112.739 12.815 0.740 0.373
    颜色空间变换 (HSV) G 103.394 6.254 0.261 0.207
    B 83.878 4.257 0.240 0.223
    NIR 108.964 4.961 0.823 0.570
    正交化变换 (GS) R 1 197.940 54.319 0.210 0.475
    G 1 294.331 55.066 0.204 0.474
    B 1 202.743 54.543 0.242 0.479
    NIR 2 007.082 57.072 0.279 0.376
    高通滤波变换 (HPF) R 1 120.665 54.261 0.403 0.287
    G 1 217.212 54.707 0.404 0.293
    B 1 125.426 55.394 0.403 0.297
    NIR 2 011.079 21.692 0.111 0.481
    主成分变换 (PC) R 1 117.781 41.135 0.170 0.357
    G 1 214.369 42.257 0.172 0.356
    B 1 122.488 41.237 0.166 0.361
    NIR 2 011.442 50.553 0.181 0.375
    超分辨率贝叶斯变换(Pansharp) R 1 086.786 53.438 0.103 0.427
    G 1 091.470 53.829 0.358 0.437
    B 1 183.465 55.879 0.358 0.432
    NIR 2 011.589 60.308 0.603 0.501
    说明:相关系数表示融合后波段与原始多光谱对应波段的相关性。
    下载: 导出CSV

    表  2  研究区各地类样本数量表

    Table  2.   Sample's number of surface features on the study area

    地类 分类样本 验证样本
    有林地 80 43
    草地 36 20
    耕地 40 21
    苗圃 20 7
    水域 25 16
    城矿居民用地 38 21
    未利用地 35 20
    总数 274 148
    下载: 导出CSV

    表  3  研究区地类分类精度对比表

    Table  3.   Contrast on classification accuracy of surface features on study area

    地类 精度(生产者精度/使用者精度)/%
    Mul Brovey HSV GS HPF PC Pansharp
    有林地 83.02/71.43 86.05/90.24 97.67/87.50 79.60/66.67 83.72/90.70 75.12/66.37 86.05/67.27
    草地 95.00/86.36 99.61/99.36 95.00/99.80 77.17/79.93 99.16/76.92 99.31/86.96 95.00/99.20
    耕地 33.33/87.5 57.14/99.80 61.90/92.86 29.05/50.00 47.62/99.24 66.67/60.87 57.14/75.00
    苗圃 86.71/54.54 99.03/50.00 85.71/85.71 71.43/45.45 57.14/57.14 71.43/71.43 98.79/50.00
    水域 43.75/98.17 87.50/98.00 93.75/99.62 99.05/94.12 93.75/93.75 93.75/88.23 93.75/93.75
    城矿居民用地 85.71/60.00 90.48/59.38 66.67/58.33 71.43/53.57 80.95/56.67 57.14/35.3 52.38/68.75
    未利用地 50.00/71.43 70.00/93.33 80.00/80.00 35.00/75.00 75.00/88.24 30.62/8.53 40.00/66.67
    总精度/Kappa系数 72.30/0.66 83.10/0.79 85.14/0.82 71.52/0.62 80.40/0.76 63.51/0.56 78.75/0.68
    下载: 导出CSV

    表  4  研究区有林地内分类精度对比表

    Table  4.   Contrast on classification accuracy of woodland on study area

    林分 精度(生产者精度/使用者精度)/%
    Mul Brovey HSV GS HPF PC Pansharp
    阔叶林 73.33/62.86 86.67/76.92 73.44/59.10 90.00/71.05 80.00/55.26 76.67/64.52 83.33/70.37
    针叶林 57.70/60.00 80.77/84.00 67.70/48.39 61.54/69.57 63.85/63.63 54.30/54.17 65.54/57.12
    混交林 22.45/33.33 70.00/46.67 54.30/38.46 32.64/40.00 44.56/50.00 32.64/27.27 40.00/36.36
    总精度/Kappa系数 59.10/0.47 75.72/0.63 66.67/0.51 68.18/0.54 71.10/0.59 64.54/0.51 72.58/0.60
    下载: 导出CSV
  • [1] 梅青. 全国林地"一张图"大数据时代的经典力作[EB/OL]. http://www.greentimes.com/green/news/yaowen/zhxw/content/2013-10/23/content_235974.htm, 2013-10-23.
    [2] 赵英时.遥感应用分析原理与方法[M].北京:科学出版社, 2003:251-262.
    [3] 李霖, 佘梦媛, 罗恒. ZY-3卫星全色与多光谱影像融合方法比较[J].农业工程学报, 2014, 30(16):157-165.

    LI Lin, SHE Mengyuan, LUO Heng. Comparison on fusion algorithms of ZY-3 panchromatic and multi-spectral images[J]. Trans Chin Soc Agric Eng, 2014, 30(16):157-165.
    [4] 杨丽萍, 夏敦胜, 陈发虎. Landsat 7 ETM+全色与多光谱数据融合算法的比较[J].兰州大学学报 (自然科学版), 2007, 43(4):7-11.

    YANG Liping, XIA Dunsheng, CHEN Fahu. On fusion algorithms of Landsat 7 ETM+PAN and multi-spectral images[J]. J Lanzhou Univ Nat Sci, 2007, 43(4):7-11.
    [5] 高志雄. 基于ZY-3卫星的森林分类及郁闭度提取方法研究[D]. 北京: 北京林业大学, 2015.

    GAO Zhixiong. Study on the Forest Classification and Crown Density Extraction Method Based on ZY-3 Image Data[D]. Beijing:Beijing Forestry University, 2015.
    [6] 郑著彬, 李俊, 任静丽. HPF图像融合技术在大理市遥感影像中的运用探讨[J].云南地理环境研究, 2007, 19(6):96-98.

    ZHENG Zhubin, LI Jun, REN Jingli. Discussion for technology of HPF resolution merge using in Dali's remote sensing image[J]. Yunnan Geogr Environ Res, 2007, 19(6):96-98.
    [7] 赵珍梅, 马伟, 王润生. 3种高保真遥感影像融合方法效果评价与分析[J].地址与勘探, 2010, 46(4):705-710.

    ZHAO Zhenmei, MA Wei, WANG Runsheng. Evaluation and analysis of there methods of fusion remote sensing images with high fidelity of information[J]. Geol Explor, 2010, 46(4):705-710.
    [8] 郑中, 祁元, 张金龙.基于光谱角和光谱距离评价指标的遥感影像融合算法比较研究:以QuickBird数据为例[J].遥感技术与应用, 2013, 28(3):437-443.

    ZHENG Zhong, QI Yuan, ZHANG Jinlong. Comparing with different remote sensing image fusion method based on evaluation index of spectral angle and spectral distance:taking Quickbird datas as example[J]. Rem Sens Technol Appl, 2013, 28(3):437-443.
    [9] 郭会敏, 洪运富, 李营, 等.基于高分一号卫星影像的多种融合算法比较[J].地理与地理信息科学, 2015, 31(1):23-26.

    GUO Huimin, HONG Yunfu, LI Ying, et al. Comparison of fusion methods used for GF-1 satellite image[J]. Geogr Geo-Inf Sci, 2015, 31(1):23-26.
    [10] 王海江, 王周龙, 李丽宏, 等.基于高平衡阶多进制多小波包变换的遥感影像融合[J].农业工程学报, 2015, 31(1):178-186.

    WANG Haijiang, WANG Zhoulong, LI Lihong, et al. Remote sensing image fusion based on high-balanced multi-band multiwavelet packet transform[J]. Trans Chin Soc Agric Eng, 2015, 31(1):178-186.
    [11] 王华斌, 李国元, 张本奎, 等.资源三号卫星影像融合算法对比分析[J].测绘科学, 2015, 40(1):47-51.

    WANG Huabin, LI Guoyuan, ZHANG Benkui, et al. Construct and analysis of different fusion algorithms for ZY-3 satellite images[J]. Sci Surv Mapp, 2015, 40(1):47-51.
    [12] 陈春华, 苏逸平, 邹崇尧. WorldView-2遥感影像融合算法实验研究[J].地理空间信息, 2013, 11(6):72-74.

    CHEN Chunhua, SU Yiping, ZOU Chongyao. Study on remote sensing image fusion algorithm[J]. Geospat Inf, 2013, 11(6):72-74.
    [13] 许民. 高分辨率遥感影像融合方法研究及融合效果评价: 以ALOS和SPOT的全色与多光谱影像为例[D]. 兰州: 兰州大学, 2010.

    XU Min. Reseach on High Resolution Remotely Sensed Image Fusion Methods and Evaluation of Fusion Quality:A Case Study of Mutlti-spectral and Panchromatic Data of ALOS and SPOT[D]. Lanzhou:Lanzhou University, 2010
    [14] 黄先德, 周群, 王兴.资源三号卫星全色与多光谱影像融合算法[J].测绘通报, 2015(1):109-114.

    HUANG Xiande, ZHOU Qun, WANG Xing. Fusion of resources satellite-3 remote sensing panchromatic and multispectral Images[J]. Bull Surv Mapp, 2015(1):109-114.
    [15] 郑丽. SPOT多光谱影像与全色影像融合研究[J].佳木斯大学学报 (自然科学版), 2015, 33(6):868-871.

    ZHENG Li. Research on the SPOT multi-spectral image and panchromatic image fusion[J]. J Jiamusi Univ Nat Sci Ed, 2015, 33(6):868-871.
    [16] 任德智, 潘刚, 葛丽雯.基于Titan Image软件的QuickBird影像融合[J].东北林业大学学报, 2013, 41(3):131-134.

    REN Dezhi, PAN Gang, GE Liwen. Image fusion of QuickBird based on Titan Image software[J]. 2013, 41(3):131-134.
    [17] 王乐, 牛雪峰, 魏斌, 等.遥感影像融合质量评价方法研究[J].测绘通报, 2015(2):77-79.

    WNAG Le, NIU Xuefeng, WEI Bin, et al. Study on quality evaluation methods for remotely sensed images fusion[J]. Bull Surv Mapp, 2015(2):77-79.
    [18] 孟京辉, 陆元昌, 刘刚, 等.基于ETM+遥感图像的图像融合试验及评价方法[J].南京林业大学学报 (自然科学版), 2010, 34(1):69-72.

    MENG Jinghui, LU Yuanchang, LIU Gang, et al. Trail of image fusion methods and quality assessment for ETM+ image[J]. J Nanjing For Univ Nat Sci Ed, 2010, 34(1), 69-72.
    [19] 王广杰, 周介铭, 杨存建, 等.基于不同算法的遥感影像融合分析[J].四川师范大学学报 (自然科学版), 2011, 34(2):255-259.

    WANG Guangjie, ZHOU Jieming, YANG Cunjian, et al. Analysis of remote sensing image fusion based on different algorithms[J]. J Sichuan Norm Univ Nat Sci, 2011, 34(2):255-259.
  • [1] 郭瑞霞, 李崇贵, 刘思涵, 马婷, 全青青.  利用多时相特征的落叶松人工林分类 . 浙江农林大学学报, 2020, 37(2): 235-242. doi: 10.11833/j.issn.2095-0756.2020.02.006
    [2] 张瑜, 陈存友, 胡希军.  应用投影寻踪分类技术的森林生态功能评价 . 浙江农林大学学报, 2020, 37(2): 243-250. doi: 10.11833/j.issn.2095-0756.2020.02.007
    [3] 郭文婷, 张晓丽.  基于Sentinel-2时序多特征的植被分类 . 浙江农林大学学报, 2019, 36(5): 849-856. doi: 10.11833/j.issn.2095-0756.2019.05.002
    [4] 李龙伟, 李楠, 陆灯盛.  多时相Sentinel-2影像在浙西北茶园信息提取中的应用 . 浙江农林大学学报, 2019, 36(5): 841-848. doi: 10.11833/j.issn.2095-0756.2019.05.001
    [5] 张兆鹏, 李增元, 田昕.  高分辨率遥感影像林地类型精细识别 . 浙江农林大学学报, 2019, 36(5): 857-867. doi: 10.11833/j.issn.2095-0756.2019.05.003
    [6] 穆亚南, 丁丽霞, 李楠, 陆琳莹, 吴明.  基于面向对象和随机森林模型的杭州湾滨海湿地植被信息提取 . 浙江农林大学学报, 2018, 35(6): 1088-1097. doi: 10.11833/j.issn.2095-0756.2018.06.012
    [7] 王怀警, 谭炳香, 房秀凤, 李世明, 李太兴.  C5.0决策树Hyperion影像森林类型精细分类方法 . 浙江农林大学学报, 2018, 35(4): 724-734. doi: 10.11833/j.issn.2095-0756.2018.04.018
    [8] 周静平, 李存军, 胡海棠, 陶欢, 彭代亮, 谢春春, 葛艳.  森林植被遥感监测影像最佳分辨率选择 . 浙江农林大学学报, 2018, 35(4): 716-723. doi: 10.11833/j.issn.2095-0756.2018.04.017
    [9] 冯云云, 刘丽娟, 陆灯盛, 庞勇.  机载高光谱影像降维方法比较 . 浙江农林大学学报, 2017, 34(5): 765-774. doi: 10.11833/j.issn.2095-0756.2017.05.001
    [10] 向安民, 刘凤伶, 于宝义, 李崇贵.  基于k-NN方法和GF遥感影像的森林蓄积量估测 . 浙江农林大学学报, 2017, 34(3): 406-412. doi: 10.11833/j.issn.2095-0756.2017.03.004
    [11] 梁文海, 刘吉凯, 张伟, 李新伟, 钟仕全.  基于面向对象方法的GF-2影像桉树林信息提取 . 浙江农林大学学报, 2017, 34(4): 721-729. doi: 10.11833/j.issn.2095-0756.2017.04.019
    [12] 董心玉, 范文义, 田甜.  基于面向对象的资源3号遥感影像森林分类研究 . 浙江农林大学学报, 2016, 33(5): 816-825. doi: 10.11833/j.issn.2095-0756.2016.05.013
    [13] 王妮, 彭世揆, 李明诗.  基于树种分类的高分辨率遥感数据纹理特征分析 . 浙江农林大学学报, 2012, 29(2): 210-217. doi: 10.11833/j.issn.2095-0756.2012.02.010
    [14] 郝红, 徐常青, 张新平.  基于非负矩阵分解的航拍图像信息提取 . 浙江农林大学学报, 2012, 29(1): 72-77. doi: 10.11833/j.issn.2095-0756.2012.01.013
    [15] 刘晓双, 黄建文, 鞠洪波.  高空间分辨率遥感的单木树冠自动提取方法与应用 . 浙江农林大学学报, 2010, 27(1): 126-133. doi: 10.11833/j.issn.2095-0756.2010.01.021
    [16] 张艮龙, 冯益明, 贾建华, 武红敢, 李增元.  基于SPOT-5遥感影像的小班区划探讨 . 浙江农林大学学报, 2010, 27(2): 299-303. doi: 10.11833/j.issn.2095-0756.2010.02.023
    [17] 李义平, 徐爱俊, 杨绍钦.  基于GIS的低丘缓坡林地管理系统的关键技术 . 浙江农林大学学报, 2009, 26(4): 554-560.
    [18] 黄初冬, 邵芸, 李静, 柳晶辉, 陈洁琼.  基于回归决策树和ASTER卫星影像的城市森林研究 . 浙江农林大学学报, 2008, 25(2): 240-244.
    [19] 白降丽, 彭道黎, 杨馥宁.  森林资源信息分类及编码体系研究 . 浙江农林大学学报, 2007, 24(3): 326-330.
    [20] 官凤英, 冯仲科, 王小昆, 周彬.  RTD GPS 量测林地面积的模拟试验 . 浙江农林大学学报, 2004, 21(2): 150-153.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2017.02.019

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2017/2/340

图(4) / 表(4)
计量
  • 文章访问数:  2597
  • HTML全文浏览量:  480
  • PDF下载量:  324
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-04-14
  • 修回日期:  2016-06-28
  • 刊出日期:  2017-04-20

面向林地分类的GF-2影像融合算法评价

doi: 10.11833/j.issn.2095-0756.2017.02.019
    基金项目:

    国家重点林业工程监测技术示范推广项目 2015-02

    作者简介:

    胡曼, 从事森林资源监测与评价的研究工作。E-mail:humanamanda@163.com

    通信作者: 彭道黎, 教授, 从事森林资源监测、森林经营管理以及林业信息领域研究。E-mail:dlpeng@bjfu.edu.cn
  • 中图分类号: S757.2

摘要: 针对林业部门目前常用的遥感影像融合算法,探究适合于林地信息提取的GF-2影像融合算法。以GF-2卫星1 m全色/4 m多光谱分辨率平面影像为数据源,采用HSV变换(hue-saturation-value,颜色空间变换),Brovey变换(彩色标准化变换),PC变换(principle components,主成分变换),HPF变换(high-pass fusion,高通滤波变换),GS变换(gram-schmidt,正交化变换),Pansharp变换(超分辨率贝叶斯变换)等6种常用融合算法,通过目视和定量特征分析对融合效果进行评价,并结合面向对象分类方法对融合后影像进行地类信息提取和分析,探讨6种融合算法对GF-2影像在林区地类提取的适宜性。研究结果表明:基于Brovey和HSV算法的融合结果目视效果良好,清晰度与纹理增强明显;这2种融合算法影像在不同地类信息的提取上各有优势,HSV算法融合结果在不同地类的提取上效果最好,分类总精度可达85.14%,Brovey算法融合结果则在森林类型的提取上具有最高的分类总精度,为75.72%;其余4种融合算法在图像质量及其他地类提取中各有优势,具体融合算法的选取需根据应用目的和影区应用区域的实际情况而定。该研究可为林业部门提高GF-2卫星的适用性及大规模应用提供参考。

English Abstract

胡曼, 彭道黎. 面向林地分类的GF-2影像融合算法评价[J]. 浙江农林大学学报, 2017, 34(2): 340-348. doi: 10.11833/j.issn.2095-0756.2017.02.019
引用本文: 胡曼, 彭道黎. 面向林地分类的GF-2影像融合算法评价[J]. 浙江农林大学学报, 2017, 34(2): 340-348. doi: 10.11833/j.issn.2095-0756.2017.02.019
HU Man, PENG Daoli. Comparison of fusion algorithms for GF-2 data from extracted forestland information[J]. Journal of Zhejiang A&F University, 2017, 34(2): 340-348. doi: 10.11833/j.issn.2095-0756.2017.02.019
Citation: HU Man, PENG Daoli. Comparison of fusion algorithms for GF-2 data from extracted forestland information[J]. Journal of Zhejiang A&F University, 2017, 34(2): 340-348. doi: 10.11833/j.issn.2095-0756.2017.02.019
  • GF-2卫星是中国第1颗空间分辨率优于1 m的民用光学遥感卫星,卫星可提供1 m全色/4 m多光谱分辨率平面影像,在森林资源调查、土地利用动态监测、环境保护与监测等领域提供服务支撑,发挥重要作用。目前,中国已经成功绘制全国林地“一张图”,是通过高分辨率遥感影像、全国林地落界数据、二类调查资料、基础地理信息等多源数据集合,以林地界线为核心内容,构建的中国林业有史以来可动态监测、及时决策的林地资源管理系统[1],取得了较好的应用效果,高分辨率遥感影像作为其中重要的一部分,主要用来提取林地及林木参数等信息,为林业数据库的建设提供了重要的基础数据支撑。因此,应用高分辨率国产卫星数据解决实际问题成为林业发展的重要趋势。影像融合是对遥感影像数据高效利用的重要手段,它综合了影像在光谱分辨率和空间分辨率上的优势,既提高了多光谱影像的分辨率,又保留其多光谱特性,从而增强解译的可靠性,提高分类精度[2]。在以往的研究中,遥感影像融合算法主要包括彩色技术、数学运算、图形变换等3种。彩色技术包括HSV变换(hue-saturation-value,颜色空间变换)等;数学运算包括Brovey变换(彩色标准化变换)等;图形变换包括PC变换(principle components,主成分变换)、HPF变换(high-pass fusion,高通滤波变换),GS变换(gram-schmidt,正交化变换)、Pansharp变换(超分辨率贝叶斯变换)等[3]。HSV算法采用全色影像替换多光谱影像的亮度分量,使HSV变换在增强多光谱影像空间细节表现能力的同时,也带来了较大的光谱失真;Brovey变换只能而且必须同时对3个波段进行融合运算,该变换在增强影像的同时保持了原多光谱影像的光谱信息[4];PC变换可以使所有的波段参与融合,可以更好地保留图像的光谱信息[5];HPF变换对于影像纹理和细节处理方面效果较好,但色彩的层次感一般[6];GS变换有明确的数学意义,对待融合的遥感影像无波段数限制,并且具有较高的光谱信息保持度;Pansharp变换对人为主观因素与数据本身质量的依赖性较小,光谱和细节特征保真度高[7]。针对不同的数据源,郑中等[8],郭会敏等[9],王海江等[10],王华斌等[11],陈春华等[12]分别对QuickBird,GF-1,资源一号02C星,资源三号,WorldView-2卫星数据全色和多光谱数据进行影像融合试验,并通过质量评价选择出适宜于影像的最优融合方法。目前研究中对于GF-2影像的融合基本未见报道,其影像数据特征与QuickBird,WorldView-2等卫星数据不完全一致,适合于现有卫星影像的融合方法不一定适合于GF-2卫星影像。因此,针对目前常用的融合方法,探究适合于林地信息提取的GF-2影像融合算法,为林业行业中大范围使用GF-2影像融合应用提供参考具有重要意义。本研究以GF-2高分辨国产卫星数据为数据源,采用6种常用的影像融合算法针对林区样地影像进行融合试验,通过目视分析和定量比较进行影像融合质量评价,并通过面向对象分类方法对不同地类进行提取分析,研究不同影像融合算法在地类信息提取上的适用性,最终选择出适宜于GF-2影像林地信息提取的融合算法,为林地信息提取、林地变化动态监测等应用提供参考。

    • 本研究区位于北京市延庆县,经纬度范围为40°27′~40°32′N,115°44′~115°52′E,研究范围为10 km × 10 km。数据使用GF-2遥感数据,采集时间2015年5月2日,基本无云,影像质量较好。包括有林地、灌木林地、水体、耕地、未利用地等多种地物类型,且各种地物所占比例相当。

      为保证影像融合的质量,全色影像与多光谱影像对应像元需保证较高的一致性,因此,试验在GIS中进行几何校正,并使用ENVI以1:10 000数字高程模型(DEM)及0.5 m航片作为基准影像进行正射校正,校正误差控制在1个像素以内。

    • 本研究采用高分辨影像应用较多的6种融合算法对GF-2全色和多光谱影像进行融合试验:HSV变换[4]、Brovey变换[4, 13]、PC变换[14]、HPF变换[6]、GS变换[15]和Pansharp变换[16]。其中HPF变换由ERDAS IMAGINE 2014软件实现,Pansharp变换基于PCI Geomatica 2013软件平台,其余融合算法均由ENVI 5.0软件实现。由于HSV变换和Brovey变换融合多光谱影像只能由3个波段参与,因此,本研究分别选择BAND2(G),BAND3(R),BAND4(NIR)与全色波段进行融合,其余融合算法均采用四波段多光谱影像和全色影像融合。

    • 为了比较不同方法的融合效果,用2种方法进行评价[16-18]:一是基于视觉效果的定性分析,关注影像的整体亮度、光谱保真度、清晰度和空间纹理特征等,这种方法简单且具有主观性,依赖于评价者的经验;二是基于数理统计的定量特征分析,能减少对融合影像评价的随机性使评价结果更客观,主要从整体影像特征和不同地类纹理特征等方面进行评价分析,整体影像特征主要包括均值、平均梯度、光谱相关系数、高频信息融入度等指标,纹理特征包括熵值、二阶矩等指标。为了更好地判断融合影像的纹理特征,本研究参照参考文献[3],选择水域、建筑用地、耕地、林地等4块不同的区域进行熵值与二阶矩指标分区域分析。本研究中,光谱相关系数和二阶矩指标的计算由ENVI软件实现,其余指标的计算都由Matlab实现。定量分析指标公式[15-19]如下。均值:

      $$ v = \frac{1}{{M \times N}}\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {F\left( {i,j} \right)} } 。 $$ (1)

      式(1)中:MN为图像的行列数,Fi, j)为融合后影像的灰度值。平均梯度:

      $$ \bar G = \frac{1}{{\left( {M - 1} \right) \times \left( {N - 1} \right)}} \times \sum\limits_{i = 1}^{M - 1} {\sum\limits_{j = 1}^{N - 1} {\sqrt {\frac{{{{\left( {{F_{\left( {i,j} \right)}} + {F_{\left( {i + 1,j} \right)}}} \right)}^2} + {{\left( {{F_{\left( {i,j} \right)}} + {F_{\left( {i,j + 1} \right)}}} \right)}^2}}}{2}} } } 。 $$ (2)

      式(2)中:MN为图像的行列数,Fi, j)为融合后影像的灰度值。光谱相关系数:

      $$ C = \frac{{\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {\left\{ {\left| {{R_{\left( {i,j} \right)}} - {v_R}} \right| \times \left| {{F_{\left( {i,j} \right)}} - {v_F}} \right|} \right\}} } }}{{\sqrt {\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {\left\{ {{{\left| {{R_{\left( {i,j} \right)}} - {v_R}} \right|}^2}} \right\}} } \times \sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {\left\{ {{{\left| {{F_{\left( {i,j} \right)}} - {v_F}} \right|}^2}} \right\}} } } }}。 $$ (3)

      式(3)中:Ri, j)为融合前多光谱影像的灰度值,vRvF分别为融合前后2幅多光谱影像的均值。高频信息融入度:

      $$ H = \frac{{\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {\left\{ {\left| {{P_{\left( {i,j} \right)}} - {v_P}} \right| \times \left| {{F_{\left( {i,j} \right)}} - {v_F}} \right|} \right\}} } }}{{\sqrt {\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {\left\{ {{{\left| {{P_{\left( {i,j} \right)}} - {v_P}} \right|}^2}} \right\}} } \times \sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {\left\{ {{{\left| {{F_{\left( {i,j} \right)}} - {v_F}} \right|}^2}} \right\}} } } }}。 $$ (4)

      式(4)中:Pi, j)为重采样到融合后影像像元大小的全色影像的灰度值,vP分别为重采样后全色影像的均值。熵值:

      $$ E = - \sum\limits_{i = 1}^L {{P_i}{{\log }_2}{P_i}} 。 $$ (5)

      式(5)中:Pi为影像中灰度值为i的概率密度,L为图像灰度级数。二阶矩:

      $$ {S_{\rm{M}}} = \sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {G{{\left( {{X_{i,j}}} \right)}^2}} } 。 $$ (6)

      式(6)中:GXij)为灰度共生矩阵元素。

    • 为探讨各方法对林区影像的适用性,研究基于eCognition软件采用面向对象的KNN分类方法,通过多尺度分割,利用对象的灰度信息、空间特征、纹理信息等,进行地类提取分析,对6种融合算法的效果进一步对比,并建立基于地类提取结果与外业实地调查结果面积的混淆矩阵,进行定量评估,评估指标包括总体精度、生产者精度、使用者精度及Kappa系数。为了进一步探究不同融合算法对不同森林类型的提取效果,研究增加有林地地类的分类结果。

    • 从各融合影像中截取1块包含多种地类的子区域(图 1)。6种算法融合后的影像空间分辨率较原始多光谱影像都有了较大的提高,更容易判读。其中HSV和HPF融合影像中建筑物、道路及水体轮廓最为清晰,更利于该地类的提取;Brovey和HSV融合算法的植被区纹理更清晰,层次感和对比度更好。PC融合算法在纹理效果上没有得到明显提高,林区边缘清晰度较低,融合效果较差;GS融合算法则在色调上有较大的改变,目视效果不佳;Pansharp融合影像颗粒感较强,目视效果一般。因此,从目视经验角度认为,对于GF-2卫星影像,使用Brovey和HSV融合算法在林地信息提取上较另外4种算法更为适合。

      图  1  融合影像目视对比图

      Figure 1.  Visual contrast on fusion images

    • 均值反映影像像素的平均灰度值即平均亮度,均值适中则色调自然。除了Brovey和HSV方法,其余4种方法各波段均值与原始影像均值接近,说明Brovey和HSV算法光谱信息损失较为严重,其余融合结果亮度适中(表 1)。平均梯度反映影像的清晰度,其值越大表示影像层次越丰富、对比度越大。在本次试验中6种方法各波段平均梯度在原始影像上都有了较大程度的提升,除了HSV融合结果的B波段外平均梯度增幅均在50%以上,说明其融合结果的清晰度和对比度提高显著。光谱相关系数反映影像光谱信息的保持性,相关系数越高说明融合对于原始影像光谱信息改变越小。HPF融合结果在R,G,B波段上的值最大(>0.4);在NIR波段上,HSV算法的值高达0.8,保持光谱特征上优于其他算法。高频信息融入度反映融合影像波段与全色波段的相关性,其值越大表示融合影像更好地继承了全色波段的空间信息特征。Pansahrp和GS算法在R,G,B波段的值都远高于其他几种算法,NIR波段上HSV和Pansahrp算法的值超过0.5,更好地继承了全色波段的空间信息。另外,植被反射率极高的NIR波段对于林地信息提取具有重要的意义,除HPF算法在NIR波段光谱相关系数低于其他波段,GS算法在NIR波段高频信息融入度低于其他波段,其余各算法在光谱相关系数和高频信息融入度上均优于其他波段,表明NIR波段在继承光谱信息和全色波段的空间信息特征上比其他波段更有优越性。HSV算法下光谱相关系数和高频信息融入度分别达到0.823和0.570,Pansharp算法次之,分别为0.603和0.501。由此推测,这2种影像融合算法在林地信息提取上具有更大的潜力。

      表 1  原多光谱影像各波段与融合影像整体影像特征对比表

      Table 1.  Contrast on indexes of each band between multi-spectral image and fusion images

      影像/方法 波段 均值 平均梯度 相关系数 高频信息融人度
      原始影像 R 1 198.725 3.124
      G 1 295.094 3.679
      B 1 203.466 4.255
      NIR 2 006.270 6.043
      PAN 305.893 12.847
      彩色标准变换 (Brovey) G 106.712 13.285 0.105 0.328
      B 86.911 9.939 0.215 0.298
      NIR 112.739 12.815 0.740 0.373
      颜色空间变换 (HSV) G 103.394 6.254 0.261 0.207
      B 83.878 4.257 0.240 0.223
      NIR 108.964 4.961 0.823 0.570
      正交化变换 (GS) R 1 197.940 54.319 0.210 0.475
      G 1 294.331 55.066 0.204 0.474
      B 1 202.743 54.543 0.242 0.479
      NIR 2 007.082 57.072 0.279 0.376
      高通滤波变换 (HPF) R 1 120.665 54.261 0.403 0.287
      G 1 217.212 54.707 0.404 0.293
      B 1 125.426 55.394 0.403 0.297
      NIR 2 011.079 21.692 0.111 0.481
      主成分变换 (PC) R 1 117.781 41.135 0.170 0.357
      G 1 214.369 42.257 0.172 0.356
      B 1 122.488 41.237 0.166 0.361
      NIR 2 011.442 50.553 0.181 0.375
      超分辨率贝叶斯变换(Pansharp) R 1 086.786 53.438 0.103 0.427
      G 1 091.470 53.829 0.358 0.437
      B 1 183.465 55.879 0.358 0.432
      NIR 2 011.589 60.308 0.603 0.501
      说明:相关系数表示融合后波段与原始多光谱对应波段的相关性。
    • 随机选择水域、建筑用地、耕地、林地等4块区域分析其熵值与二阶矩指标,其中,林地、耕地与建筑用地50 m×50 m样地各30块,水域50 m×50 m样地15块,10 m×10 m样地10块。熵值表示影像所包含的信息量,熵值越大则影像空间细节丰富。分析熵值指标统计图(图 2),影像融合对于林地的熵值增加最为明显,HSV,Pansharp等2种融合算法信息熵分别为2.10,1.82,较原多光谱影像信息熵(1.45)均提高达25%以上,信息量改善良好。二阶矩又叫能量,二阶矩越小说明影像灰度分布均匀、纹理细致。分析二阶矩指标统计图(图 3)发现林地的二阶矩增益较为明显。在林地中,HSV和Pansharp融合影像的二阶矩分别为0.13,0.19,较原始多光谱影像二阶矩(0.31)均有显著改善,其中HSV改善效率超过50%。综合2个纹理特征指标,HSV和Pansharp算法对纹理具有稳定和显著的改善,其中HSV在增益效果上更具优势,融合影像信息量最大,纹理最为细致。

      图  2  熵值指标统计图

      Figure 2.  Statistical chart of entropy indexes

      图  3  二阶矩指标统计图

      Figure 3.  Statistical chart of second moment indexes

    • 为了探究各融合算法下的影像质量及在林地信息提取上的应用效果,本研究利用面向对象多尺度分割的方法,结合光谱、纹理和空间信息等特征值进行KNN分类,分割尺度为220,形状因子0.1,紧实度因子0.5,分类结果部分区域如图 4所示。本研究利用原始多光谱影像,结合2014年森林资源二类调查数据,选择7种地类的典型样本共274个作为分类样本,并实测全球定位系统(GPS)点148个作为验证样本(表 2)。使用eCognition软件建立基于样本的混淆矩阵,分类精度见表 3

      图  4  融合影像分类结果图

      Figure 4.  Classification results of fusion images

      表 2  研究区各地类样本数量表

      Table 2.  Sample's number of surface features on the study area

      地类 分类样本 验证样本
      有林地 80 43
      草地 36 20
      耕地 40 21
      苗圃 20 7
      水域 25 16
      城矿居民用地 38 21
      未利用地 35 20
      总数 274 148

      表 3  研究区地类分类精度对比表

      Table 3.  Contrast on classification accuracy of surface features on study area

      地类 精度(生产者精度/使用者精度)/%
      Mul Brovey HSV GS HPF PC Pansharp
      有林地 83.02/71.43 86.05/90.24 97.67/87.50 79.60/66.67 83.72/90.70 75.12/66.37 86.05/67.27
      草地 95.00/86.36 99.61/99.36 95.00/99.80 77.17/79.93 99.16/76.92 99.31/86.96 95.00/99.20
      耕地 33.33/87.5 57.14/99.80 61.90/92.86 29.05/50.00 47.62/99.24 66.67/60.87 57.14/75.00
      苗圃 86.71/54.54 99.03/50.00 85.71/85.71 71.43/45.45 57.14/57.14 71.43/71.43 98.79/50.00
      水域 43.75/98.17 87.50/98.00 93.75/99.62 99.05/94.12 93.75/93.75 93.75/88.23 93.75/93.75
      城矿居民用地 85.71/60.00 90.48/59.38 66.67/58.33 71.43/53.57 80.95/56.67 57.14/35.3 52.38/68.75
      未利用地 50.00/71.43 70.00/93.33 80.00/80.00 35.00/75.00 75.00/88.24 30.62/8.53 40.00/66.67
      总精度/Kappa系数 72.30/0.66 83.10/0.79 85.14/0.82 71.52/0.62 80.40/0.76 63.51/0.56 78.75/0.68

      从总体精度和Kappa系数指标来看,HSV变换融合结果分类精度最高,总精度为85.14%,Kappa系数0.82,较原始多光谱影像分类精度提高了约13%,其次是Brovey,HPF及Pansharp变换,PC变换融合影像分类精度为63.51%,低于原始多光谱影像分类精度(72.30%),说明对于GF-2影像来说并不是所有融合变换都能增强影像信息提取的能力。

      在对有林地信息提取精度上,HSV和Brovey变换生产者精度分别为97.67%,86.05%,使用者精度分别为87.50%,90.24%,较其他融合变换有明显的优势,GS变换和PC变换不适宜提取林地,生产者角度和使用者精度均较低。

      对草地分类中,除GS变换和HPF变换使用者精度较低外,其余变换均较高,Brovey,HSV和Pansharp变换使用者精度都超过90%,适宜草地的提取。对耕地提取中,除GS变换外,其余5种变换都较原始多光谱影像分类有所提高,其中Brovey变换使用者精度最高,PC变换生产者精度最高,适于耕地提取。苗圃地提取分析中,Brovey和Pansharp变换生产者精度较高都在95%以上,HSV变换使用者精度最高约85.71%。其他3种方法不适宜苗圃地类提取。水域提取分析中,对比其他地类提取,6种融合图像的生产者精度和使用者精度均较高。相比较而言PC变换不适宜提取水体,生产者精度较低。城矿居民用地提取分析中仅Brovey变换在生产者精度上有了显著的提升,适于城矿居民用地的提取。

      为了进一步研究不同融合变换影像对有林地地类中森林类型信息提取效果,本研究利用上一步中有林地地类矢量与影像中NDVI大于0.35且Brightness小于1 500的矢量叠加,重叠区作为有林地范围,进行再次分割和KNN分类,分为阔叶林、针叶林、混交林3个地类。分割尺度为100,形状因子0.1,紧实度因子0.5。分类样本共128个,并实测GPS点阔叶林30个、针叶林30个以及混交林30个作为验证样本,使用eCognition软件建立基于样本的混淆矩阵,分类精度见表 4

      表 4  研究区有林地内分类精度对比表

      Table 4.  Contrast on classification accuracy of woodland on study area

      林分 精度(生产者精度/使用者精度)/%
      Mul Brovey HSV GS HPF PC Pansharp
      阔叶林 73.33/62.86 86.67/76.92 73.44/59.10 90.00/71.05 80.00/55.26 76.67/64.52 83.33/70.37
      针叶林 57.70/60.00 80.77/84.00 67.70/48.39 61.54/69.57 63.85/63.63 54.30/54.17 65.54/57.12
      混交林 22.45/33.33 70.00/46.67 54.30/38.46 32.64/40.00 44.56/50.00 32.64/27.27 40.00/36.36
      总精度/Kappa系数 59.10/0.47 75.72/0.63 66.67/0.51 68.18/0.54 71.10/0.59 64.54/0.51 72.58/0.60

      可见,Brovey和Pansharp变换融合结果在有林地内分类精度最高,总体精度分别达到75.72%和72.58%,Kappa系数分别达到0.63和0.60,其次是HPF和GS变换,PC变换融合结果分类精度最低。阔叶林分类生产者精度普遍高于针叶林,GS,Brovey,Pansharp和HPF变换均达到80%以上,混生林总体分类精度低,仅有Brovey和HPF变换融合结果在生产者精度及使用者精度均能达到50%以上。因此,在森林类型层次的林地信息提取中采用Brovey或Pansharp融合算法更为适宜。

    • 影像融合是为了使不同数据源实现优势互补,从而对影像数据高效利用。针对基于GF-2影像面向林地信息提取的遥感影像融合算法选择问题,本研究选择6种具有代表性的融合算法(HSV变换,Brovey变换,PC变换,HPF变换,GS变换,Pansharp变换),在采用目视与定量特征分析的传统影像质量评价基础上,辅助面向对象分类林地信息效果分析。在融合影像整体质量方面,HSV和Pansharp变换目视效果较好,其中Pansharp变换在在光谱继承和增强信息纹理方面取得了较好的平衡,光谱相关系数及高频信息融入度都较高,对纹理具有稳定和显著的改善。传统的目视及定量指标分析有一定的合理性,但也缺乏针对性,难以对不同层次的森林类型的信息提取作出比较,具有局限性。

      在基于融合影像的面向对象分类中,HSV,Brovey,HPF变换总体分类精度较高,其余3种方法在个别类型中具有明显劣势。对于林地信息提取的最优融合算法需根据应用的具体层次来选择,在有林地层次的信息提取中HSV变换较为适宜,在森林类型信息提取则采用Brovey和Pansharp变换更为有效。因此,综合融合影像整体质量及分类精度,针对GF-2影像的林地信息提取,HSV及Pansharp变换的融合效果优于其他4种算法,可为GF-2影像在林业行业大规模融合提供参考。不同融合方法对不同地类信息的提取各有优势。在进行森林动态监测、地类信息提取、地块变化等不同林业遥感工作中,需考虑具体情况及地类信息比例选择合适的融合方法。

    • 感谢中国林业科学研究院资源信息研究所为本研究提供的GF-2卫星影像数据。

参考文献 (19)

目录

    /

    返回文章
    返回