留言板

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

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

云环境下基于CA-Markov的土地利用变化预测方法

康俊锋 李爽 方雷

康俊锋, 李爽, 方雷. 云环境下基于CA-Markov的土地利用变化预测方法[J]. 武汉大学学报 ● 信息科学版, 2020, 45(7): 1021-1026, 1034. doi: 10.13203/j.whugis20180319
引用本文: 康俊锋, 李爽, 方雷. 云环境下基于CA-Markov的土地利用变化预测方法[J]. 武汉大学学报 ● 信息科学版, 2020, 45(7): 1021-1026, 1034. doi: 10.13203/j.whugis20180319
KANG Junfeng, LI Shuang, FANG Lei. Land Use Change Prediction Method Based on CA-Markov Model Under Cloud Computing Environment[J]. Geomatics and Information Science of Wuhan University, 2020, 45(7): 1021-1026, 1034. doi: 10.13203/j.whugis20180319
Citation: KANG Junfeng, LI Shuang, FANG Lei. Land Use Change Prediction Method Based on CA-Markov Model Under Cloud Computing Environment[J]. Geomatics and Information Science of Wuhan University, 2020, 45(7): 1021-1026, 1034. doi: 10.13203/j.whugis20180319

云环境下基于CA-Markov的土地利用变化预测方法

doi: 10.13203/j.whugis20180319
基金项目: 

国家重点研发计划 2016YFC0803105

国家自然科学基金 41701462

国家留学基金 201808360065

江西省教育厅科学技术研究项 3204704062

详细信息

Land Use Change Prediction Method Based on CA-Markov Model Under Cloud Computing Environment

Funds: 

The National Key Research and Development Program of China 2016YFC0803105

the National Natural Science Foundation of China 41701462

Chinese Scholarship Council Foundation 201808360065

Jiangxi Provincial Department of Education Science and Technology Research Projects 3204704062

More Information
    Author Bio:

    KANG Junfeng, PhD, associate professor, specializes in researches of high-performance GIS algorithms and applications. E-mail: junfeng.kang@jxust.edu.cn

    Corresponding author: FANG Lei, PhD. E-mail: fanglei@fudan.edu.cn
  • 摘要: 传统土地利用变化预测方法通常效率低, 无法满足土地利用变化大数据分析和处理的需求。采用MapReduce编程模型对元胞自动机-马尔可夫模型进行并行化扩展,设计了基于Hadoop的土地利用变化预测方法(land-use change prediction method based on cloud computing,Cloud-CMLP),并选取杭州市进行实验,包括:①不同数据量下Cloud-CMLP核心算法的运行效率实验;②利用Cloud-CMLP方法模拟杭州2013年土地利用变化,并将模拟结果与2013年遥感影像解译结果进行对比分析,验证了预测方法的正确性;③杭州2020年的土地利用变化模拟预测分析, 研究区中心地带受发达交通系统影响, 建设用地面积整体呈快速上升趋势, 且主要来源于农业用地的转换。
  • 图  1  基于MapReduce的CA-Markov并行算法

    Figure  1.  Land Use Change Prediction Algorithm of CA-Markov Based on MapReduce

    图  2  Cloud-Markov与串行Markov效率对比

    Figure  2.  Running Efficiency Between Cloud-Markov and Serial Markov Algorithms

    图  3  Cloud-Markov与串行Markov加速比

    Figure  3.  Speedup Ratio of Cloud-Markov and Serial Markov Algorithms

    图  4  Cloud-CELUC与串行CELUC效率对比

    Figure  4.  Running Efficiency Between Cloud-CELUC and Serial CELUC Algorithms

    图  5  Cloud-CELUC与串行CELUC加速比

    Figure  5.  Speedup Ratio of Cloud-CELUC and Serial CELUC Algorithms

    图  6  2020年研究区土地利用变化预测图

    Figure  6.  Land Use Change Prediction of Research Area in 2020

    表  1  2006—2013年面积转移矩阵/km2

    Table  1.   Area Transfer Matrix in 2006—2013/km2

    地类 农用地 建设用地 自然保留地 总计
    农用地 1 282.95 409.71 95.67 1 788.33
    建设用地 210.94 1 381.69 12.52 1 605.15
    自然保留地 93.39 50.79 4 409.60 4 553.78
    总计 1 587.28 1 842.19 4 517.79 7 947.26
    下载: 导出CSV

    表  2  2006—2013年转移概率矩阵/%

    Table  2.   Transition Probability Matrix in 2006—2013/%

    地类 农用地 建设用地 自然保留地
    农用地 71.74 22.91 5.35
    建设用地 13.14 86.08 0.78
    自然保留地 2.05 1.12 96.83
    下载: 导出CSV

    表  3  自然保留地Kappa系数检验表/km2

    Table  3.   Kappa Coefficient of Natural Reserve/km2

    地类 自然保留地 非自然保留地 总计 准确度/% Kappa系数
    自然保留地 4 221.35 288.47 4 509.82 93.60 0.86
    非自然保留地 296.44 3 431.50 3 727.94 92.05
    总计 4 517.79 3 717.97
    下载: 导出CSV

    表  4  建设用地Kappa系数检验表/km2

    Table  4.   Kappa Coefficient of Construction Land/km2

    地类 建设用地 非建设用地 总计 准确度/% Kappa系数
    建设用地 1 391.41 452.77 1 844.18 75.45 0.68
    非建设用地 450.78 5 942.80 6 393.58 92.95
    总计 1 842.19 6 395.57
    下载: 导出CSV

    表  5  农用地Kappa系数检验表/km2

    Table  5.   Kappa Coefficient of Agricultural Land/km2

    地类 农用地 非农用地 总计 准确度/% Kappa系数
    农用地 1 152.64 427.17 1 579.81 72.96 0. 66
    非农用地 434.65 6 223.30 6 657.95 93.47
    总计 1 587.29 6 650.47
    下载: 导出CSV

    表  6  2013—2020年土地利用面积转移预测矩阵/km2

    Table  6.   Area Transfer Matrix in 2013—2020/km2

    地类 农用地 建设用地 自然保留地 总计
    农用地 1 133.35 361.94 84.52 1 579.81
    建设用地 242.35 1 587.45 14.38 1 844.18
    自然保留地 92.49 50.30 4 367.03 4 509.82
    总计 1 468.19 1 999.69 4 465.93 7 933.81
    下载: 导出CSV
  • [1] 李秀彬.全球环境变化研究的核心领域—土地利用/土地覆被变化的国际研究动向[J].地理学报, 1996, 51(6):558-565

    Li Xiubin. A Review of the International Researches on Land Use/Land Cover Change[J]. Acta Geographica Sinica, 1996, 51(6): 558–565
    [2] Wijesekara G N, Farjad B, Cupta A, et al. A Comprehensive Land-Use/Hydrological Modeling System for Scenario Simulations in the Elbow River Watershed, Alberta, Canada[J]. Environmental Mana-gement, 2014, 53(2): 357-381 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=0d74809ebe15ec588c323350c7e60b7b
    [3] Verburg P H, Overmars K P. Combining Top-Down and Bottom-Up Dynamics in Land Use Modeling: Exploring the Future of Abandoned Farmlands in Europe with the Dyna-CLUE Model[J]. Landscape Ecology, 2009, 24(9): 1 167 doi:  10.1007/s10980-009-9355-7
    [4] Luo G, Yin C, Chen X, et al. Combining System Dynamic Model and CLUE-S Model to Improve Land Use Scenario Analyses at Regional Scale: A Case Study of Sangong Watershed in Xinjiang, China[J]. Ecological Complexity, 2010, 7(2): 198-207
    [5] Clarke K C, Hoppen S, Gaydos L, et al. Modifying Historical Urbanization in the San Franeisco Bay Area[J]. Environment and Planning B: Planning and Design, 1997, 24(2): 247-261 doi:  10.1068/b240247
    [6] Fu X, Wang X H, Yang Y J. Deriving Suitability Factors for CA-Markov Land Use Simulation Model Based on Local Historical Data[J]. Journal of Environmental Management, 2018, 206: 10-19 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=4cbce65eb9c77d2d6955a3772e46fdb6
    [7] Le Q B, Park S J, Vlek P L G, et al. Land-Use Dynamic Simulator (LUDAS): A Multi-Agent System Model for Simulating Spatio-Temporal Dynamics of Coupled Human-Landscape System. I. Structure and Theoretical Specification[J]. Ecological Informatics, 2008, 3(2): 135-153 doi:  10.1016/j.ecoinf.2008.04.003
    [8] Dietzel C, Clarke K C. Toward Optimal Calibration of the SLEUTH Land Use Change Model[J]. Transactions in GIS, 2007, 11(1): 29-45 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1111/j.1467-9671.2007.01031.x
    [9] Hu Z, Lo C P. Modeling Urban Growth in Atlanta Using Logistic Regression[J]. Computers, Environment and Urban Systems, 2007, 31(6): 667-688 doi:  10.1016/j.compenvurbsys.2006.11.001
    [10] Memarian H, Balasundram S K, Talib J B, et al. Validation of CA-Markov for Simulation of Land Use and Cover Change in the Langat Basin, Malaysia[J]. Journal of Geographic Information System, 2012, 4(6): 542-554 doi:  10.4236/jgis.2012.46059
    [11] Liu X H, Andersson C.Assessing the Impact of Temporal Dynamics on Land-Use Change Modeling[J]. Computers, Environment and Urban Systems, 2004, 28:107-124 doi:  10.1016/S0198-9715(02)00045-5
    [12] Keeley J E, Zedler P H. Large, High-Intensity Fire Events in Southern California Shrublands: Debunking the Fine-Grain Age Patch Model[J]. Ecological Applications, 2009, 19(1): 69-94
    [13] Araya Y H. Analysis and Modeling of Urban Land Cover Change in Setúbal and Sesimbra, Portugal[J]. Remote Sensing, 2010, 2(6): 1 549-1 563 doi:  10.3390/rs2061549
    [14] Guan Q F, Clarke K C. A General-Purpose Parallel Raster Processing Programming Library Test Application Using a Geographic Cellula Automata Model[J]. International Journal of Geographical Information Science, 2010, 24 (5):695-722 doi:  10.1080/13658810902984228
    [15] Li X, Zhang X H, Ye H A, et al. Parallel Cellular Automata for Large-Scale Urban Simulation Using Load-Balancing Techniques[J]. International Journal of Geographical Information Science, 2010, 24(6):803-820 doi:  10.1080/13658810903107464
    [16] 李丹, 黎夏, 刘小平, 等. GPU-CA模型及大尺度土地利用变化模拟[J].科学通报, 2012(11): 959-969 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=kxtb201211011

    Li Dan, Li Xia, Liu Xiaoping, et al. GPU-CA Model and Large-scale Land Use Change Simulation[J]. Chinese Science Bulletin, 2012(11):959-969 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=kxtb201211011
    [17] Wang H J, Zhang A Q, Zhang W T, et al.A High-Performance Cellular Automata Model for Urban Simulation Based on Vectorization and Parallel Computing Technology[J]. International Journal of Geographical Information Science, 2018, 32(2): 399-424 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1080/13658816.2017.1390118
    [18] Guan Q F, Shi X, Huang M Q, et al. A Hybrid Parallel Cellular Automata Model for Urban Growth Simulation over GPU/CPU Heterogeneous Architectures[J]. International Journal of Geographical Information Science, 2016, 30 (3): 494-514 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1080/13658816.2015.1039538
    [19] 姚晓闯, 杨建宇, 李林, 等.云环境下海量空间矢量数据并行划分算法[J].武汉大学学报·信息科学版, 2018, 43(7):1 092-1 097 doi:  10.13203/j.whugis20160271

    Yao Xiaochuang, Yang Jianyu, Li Lin, et al. Parallel Algorithm for Partitioning Massive Spatial Vector Data in Cloud Environment[J]. Geomatics and Information Science of Wuhan University, 2018, 43(7):1 092-1 097 doi:  10.13203/j.whugis20160271
    [20] Almeer M H. Cloud Hadoop Map Reduce for Remote Sensing Image Analysis[J]. Journal of Emerging Trends in Computing and Information Science, 2012, 3(4):637-644
    [21] 吴翼平.多准则评价方法的评价[J].决策与决策支持系统, 1994(1):74-80 http://d.old.wanfangdata.com.cn/Periodical/xnjtdxxb200301021

    Wu Yiping. The Evaluation of the Multieriteria Evaluation Methods[J]. Journal of Decision and Decision Support Systems, 1994(1):74-80 http://d.old.wanfangdata.com.cn/Periodical/xnjtdxxb200301021
  • [1] 姚晓闯, 杨建宇, 李林, 叶思菁, 郧文聚, 朱德海.  云环境下海量空间矢量数据并行划分算法 . 武汉大学学报 ● 信息科学版, 2018, 43(7): 1092-1097. doi: 10.13203/j.whugis20160271
    [2] 徐文, 邵俊, 喻文勇, 冯鹏.  陆地观测卫星数据中心:大数据挑战及一种解决方案 . 武汉大学学报 ● 信息科学版, 2017, 42(1): 7-13. doi: 10.13203/j.whugis20150172
    [3] 何亚文, 魏海涛, 杜云艳.  基于云计算的海洋环境数据共享平台设计与实现 . 武汉大学学报 ● 信息科学版, 2016, 41(10): 1307-1312. doi: 10.13203/j.whugis20140466
    [4] 王成, 施闯, 段兵兵, 刘晓莉.  利用位置与云计算建立电离层TEC主动式播发平台 . 武汉大学学报 ● 信息科学版, 2016, 41(3): 304-309. doi: 10.13203/j.whugis20130307
    [5] 段苗苗, 马盈盈, 龚威, 王伦澈.  利用微波辐射计观测计算云衰减 . 武汉大学学报 ● 信息科学版, 2015, 40(12): 1606-1612. doi: 10.13203/j.whugis20130820
    [6] 李德仁, 姚远, 邵振峰.  智慧城市中的大数据 . 武汉大学学报 ● 信息科学版, 2014, 39(6): 631-640. doi: 10.13203/j.whugis20140135
    [7] 黄斌, 彭宇行, 彭小宁.  云计算环境中高效分布式索引技术 . 武汉大学学报 ● 信息科学版, 2014, 39(11): 1375-1381.
    [8] 利用MapReduce进行批量遥感影像瓦片金字塔构建 . 武汉大学学报 ● 信息科学版, 2013, 38(3): 278-282.
    [9] 冯梦喆, 何建华, 汤青慧.  土地利用变化反事实情景模拟与评价——以湖北省嘉鱼县为例 . 武汉大学学报 ● 信息科学版, 2013, 38(9): 1122-1125.
    [10] 陈江平, 韩青, 胡晶, 张鹏林.  顾及小波变换的土地利用变化与经济因子的多尺度相关性分析 . 武汉大学学报 ● 信息科学版, 2013, 38(9): 1118-1121.
    [11] 刘 义, 陈 荦, 景 宁, 熊 伟.  利用MapReduce进行批量遥感影像瓦片金字塔构建 . 武汉大学学报 ● 信息科学版, 2013, 38(3): 278-282.
    [12] 李德仁, 眭海刚, 单杰.  论地理国情监测的技术支撑 . 武汉大学学报 ● 信息科学版, 2012, 37(5): 505-512.
    [13] 李雪, 舒宁, 王琰, 李亮.  利用土地利用状态转移分析的变化检测 . 武汉大学学报 ● 信息科学版, 2011, 36(8): 952-955.
    [14] 农宇, 王坤, 杜清运.  利用多分类Logistic回归进行土地利用变化模拟——以湖北省嘉鱼县为例 . 武汉大学学报 ● 信息科学版, 2011, 36(6): 743-746.
    [15] 徐昔保, 杨桂山, 张建明.  兰州市城市土地利用优化研究 . 武汉大学学报 ● 信息科学版, 2009, 34(7): 878-881.
    [16] 熊华, 刘耀林, 车珊珊, 张俊峰.  基于支持向量机的土地利用变化模拟模型 . 武汉大学学报 ● 信息科学版, 2009, 34(3): 366-369.
    [17] 王汉花, 刘艳芳.  基于MOP-CA整合模型的土地利用优化研究 . 武汉大学学报 ● 信息科学版, 2009, 34(2): 174-177.
    [18] 杨国清, 刘耀林, 吴志峰.  基于CA-Markov模型的土地利用格局变化研究 . 武汉大学学报 ● 信息科学版, 2007, 32(5): 414-418.
    [19] 肖平, 李德仁.  基于人工神经元网络技术的土地利用/覆盖变化探测 . 武汉大学学报 ● 信息科学版, 2002, 27(6): 586-590.
    [20] 贾华, 祝国瑞, 佐藤洋平.  土地利用变化研究中的细胞自动机与灰色局势决策 . 武汉大学学报 ● 信息科学版, 1999, 24(2): 166-169.
  • 加载中
图(6) / 表(6)
计量
  • 文章访问数:  317
  • HTML全文浏览量:  63
  • PDF下载量:  36
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-06-20
  • 刊出日期:  2020-07-30

云环境下基于CA-Markov的土地利用变化预测方法

doi: 10.13203/j.whugis20180319
    基金项目:

    国家重点研发计划 2016YFC0803105

    国家自然科学基金 41701462

    国家留学基金 201808360065

    江西省教育厅科学技术研究项 3204704062

    作者简介:

    康俊锋,博士,副教授,主要从事高性能GIS算法及应用研究。junfeng.kang@jxust.edu.cn

    通讯作者: 方雷,博士。fanglei@fudan.edu.cn
  • 中图分类号: P208

摘要: 传统土地利用变化预测方法通常效率低, 无法满足土地利用变化大数据分析和处理的需求。采用MapReduce编程模型对元胞自动机-马尔可夫模型进行并行化扩展,设计了基于Hadoop的土地利用变化预测方法(land-use change prediction method based on cloud computing,Cloud-CMLP),并选取杭州市进行实验,包括:①不同数据量下Cloud-CMLP核心算法的运行效率实验;②利用Cloud-CMLP方法模拟杭州2013年土地利用变化,并将模拟结果与2013年遥感影像解译结果进行对比分析,验证了预测方法的正确性;③杭州2020年的土地利用变化模拟预测分析, 研究区中心地带受发达交通系统影响, 建设用地面积整体呈快速上升趋势, 且主要来源于农业用地的转换。

English Abstract

康俊锋, 李爽, 方雷. 云环境下基于CA-Markov的土地利用变化预测方法[J]. 武汉大学学报 ● 信息科学版, 2020, 45(7): 1021-1026, 1034. doi: 10.13203/j.whugis20180319
引用本文: 康俊锋, 李爽, 方雷. 云环境下基于CA-Markov的土地利用变化预测方法[J]. 武汉大学学报 ● 信息科学版, 2020, 45(7): 1021-1026, 1034. doi: 10.13203/j.whugis20180319
KANG Junfeng, LI Shuang, FANG Lei. Land Use Change Prediction Method Based on CA-Markov Model Under Cloud Computing Environment[J]. Geomatics and Information Science of Wuhan University, 2020, 45(7): 1021-1026, 1034. doi: 10.13203/j.whugis20180319
Citation: KANG Junfeng, LI Shuang, FANG Lei. Land Use Change Prediction Method Based on CA-Markov Model Under Cloud Computing Environment[J]. Geomatics and Information Science of Wuhan University, 2020, 45(7): 1021-1026, 1034. doi: 10.13203/j.whugis20180319
  • 土地利用变化与人类各类活动相互影响和作用,研究不同时空下的土地利用变化,预测未来的土地利用结构和空间布局,可为区域土地资源的科学利用、区域生态环境的保护与社会经济可持续发展提供科学参考[1-2]

    很多学者提出了各自的土地利用变化模拟与预测模型,如CLUE(conversion of land use and its effects)模型[3]、CLUE-S(CLUE at small region extent)模型[4]、元胞自动机(cellular automata, CA)[5]、马尔可夫(Markov)链[6]、系统动力学[7]、SLEUTH[8]、空间Logistic模型[9],CA-Markov模型等。其中CA-Markov模型结合了Markov模型长期预测的能力和CA模型模拟复杂系统空间变化的能力,已被广泛验证[10-13]。然而土地利用变化模拟预测方法涉及大量的数据和计算,传统的串行算法或者半手工计算方法效率较低,有学者设计了单机并行[14]、MPI(message-passing-inter face)[15]、GPU(graphic processing unit)[16-17]、GPU/CPU(central processing unit)混合架构[18]的并行CA模型,应用于城市变化模拟,并验证了其算法性能。但近年来, 土地利用大数据的应用相比原有一城一地,在空间范围、时间跨度、空间分辨率和时间分辨率上都有巨大飞跃,目前,采用MapReduce对传统算法进行并行化处理以提高运行效率,已成为解决这一问题的主流方法和思路,并在很多方面得到了应用和验证[19-20]

    本文首次采用MapReduce对基于CA-Markov的土地利用变化预测方法(land-use change prediction method based on cloud computing,Cloud-CMLP)进行改进,在兼顾预测正确性的基础上,提高了土地利用变化预测的效率,实现了云环境下的土地利用变化预测。

    • 采用MapReduce编程模型对CA-Markov模型进行并行化处理,设计云环境下基于CA-Markov模型的Cloud-CMLP,其可并行化部分主要包括两个算法:云环境下并行Markov算法(parallel Markov algorithm on cloud computing, Cloud-Markov)及云环境下土地利用变化综合评价算法(land use change comprehensive evaluation algorithm on cloud computing, Cloud-CELUC)。

      图 1所示,利用Cloud-Markov对研究区进行土地利用类型变化总量预测,获得各地类数量变化; 利用Cloud-CELUC来确定每个元胞对应各地类转变方向的评价值,得到所有评价值(即土地利用综合评价值统计总表),再结合第一部分Cloud-Markov的变化总量预测结果,判定所有元胞转换方向,从而实现土地利用变化空间分布预测。

      图  1  基于MapReduce的CA-Markov并行算法

      Figure 1.  Land Use Change Prediction Algorithm of CA-Markov Based on MapReduce

      其中,Cloud-CELUC包括云环境下邻域影响度评价(cloud-based evaluation of neighborhood influence degree, C-ENID)、云环境下多准则评价(cloud-based multi-criteria evaluation, C-MCE)及评价值综合计算3个部分。C-ENID利用MapReduce基于CA模型计算得到邻域影响度评价值,C-MCE利用MapReduce基于MCE计算得到土地利用约束性评价值(constraint evaluation value, CEV)与土地利用适宜性评价值(suitability evaluation value,SEV),评价值综合计算使用前两部分计算得到的约束性评价值与适宜性评价值、邻域影响度评价值计算得到每个元胞转换为每个地类的评价值。

    • 将相同空间位置不同时相的两幅栅格影像进行叠加分析,获取每个栅格的地类转换方向,并统计各转换方向的栅格数量,得到地类转换的面积转移矩阵,并计算得到土地利用面积转移概率矩阵。但通过该概率矩阵只能获得研究初期至末期时间间隔m年的各类土地利用面积,将该面积值除以m,得到每年各类土地利用转换面积,可以得到一年间隔的面积转移矩阵。研究给出的Markov模型的Map函数、Combiner函数和Reduce函数的定义公式为:

      $$ \left\{\begin{array}{l} \operatorname{Map}:(N, (T 1, T 2)) \rightarrow \operatorname{List}\left(C_{m k}, i\right) \\ \text { Combiner: } M\left(C_{m k}, \operatorname{List}(i)\right) \rightarrow \operatorname{List}\left(C_{m k}, s\right) \\ \operatorname{Reduce}: L\left(C_{m k}, \operatorname{List}(s)\right) \rightarrow \operatorname{List}\left(C_{m k}, q\right) \end{array}\right. $$ (1)

      式中,N为输入行的行偏移量; T1表示土地利用初期数据栅格地类; T2表示土地利用末期数据栅格的地类; Cmk表示栅格由m地类转换为k地类; i表示栅格发生Cmk转换的个数; M为节点本地数据有M种地类; s为节点本地该种转换方向的栅格总数; L为地类总数; q为整个研究区发生Cmk转换的栅格总数与转移概率的组合值。

    • 1) 元胞邻域读取处理

      邻域影响度评价值的计算需要邻域元胞的状态数据,常用的元胞邻域有Von Neumann型和Moore型邻域,本文算法采用Moore 3×3型邻域进行设计,每个中心元胞都以4周8个方向的元胞为邻居。以mn列的栅格影像为例,各元胞是按行顺序依次读取,并记录每个元胞的序号Kx(其中K为行顺序号,下标x为列顺序号),因此Kx的邻域元胞为(K-1)x-1、(K-1)x、(K-1)x+1Kx-1Kx+1、(K+1)x-1、(K+1)x、(K+1)x+1

      基于MapReduce计算CA模型时,需要查找各元胞的邻域元胞,因此为更好地减少MapReduce计算过程中节点间的数据交换,元胞读取时是以行为单位转换为一维数组,同时保留了各元胞的行列号,因此可快速地获取任意元胞的邻域元胞。

      2) 多准则评价因子确定

      研究利用MCE计算约束性评价值与适宜性评价值,MCE是可在具有相互冲突、不可共度的有限(无限)方案集中进行选择的方法,其目的在于确定最优决策方案[21]。评价标准分为适宜性因子和约束性因子两种。适宜性因子是将影响因素标准化为可度量的连续值,进而评价区域并进行选择。约束性因子是将分析要素空间化处理为布尔值,标准只有“0或者1”两种选择,没有中间灰色的地带。

      通过层次分析法并结合专家打分,判定不同适宜性因子的权重。另外,无论上述适宜性分析结果如何,有元胞所对应的土地类型都不会发生变化,则转变概率为0。约束性因子为0或1,不受约束性因子影响的其值为1,否则为0,例如坡度大于25°的山体,水域,生态保护区的约束性因子设置为0。

      3) Cloud-CELUC算法

      C-MCE的计算需要每个元胞的位置与各因子数据位置一一对应,按相应的行号将两者合并。获得邻域影响度评价值、约束性评价值与适宜性评价值后,进行综合评价值计算。所有节点计算出所有元胞的地类转换综合评价值,即为土地利用综合评价值统计总表。给出C-ENID与C-MCE的Map函数定义如下:

      $$ \begin{array}{*{20}{c}} {(N, (i, H1, H2, H3 \cdots Hm)) \to ((i, j), L\mathit{1}, }\\ {L2 \cdots Lm)} \end{array} $$ (2)

      式中,N为输入行的行偏移量; i表示影像中行号; H1为需计算的元胞行所有元胞状态值; H2与H3表示H1的上下行元胞状态; H4、H5…Hm等为该元胞行相对应的各种约束性因子数据与适宜性因子数据; (i, j)为元胞在影像中的行列号; Lm为该元胞转换方向与对应转换方向综合评价值CELUC的组合值。

    • 采用多目标土地利用竞争方法实现元胞土地类型转换,该方法是解决同一元胞发生土地利用类型转换方向冲突时的抉择方法。以土地利用综合评价值统计总表为基础,并以土地利用变化面积转移矩阵为约束,进行各元胞转换方向确定。其基本算法如下:

      1) 忽略元胞地类转换方向,计算全体元胞空间中评价值最大值,记录其元胞位置以及地类转换方向。

      2) 该种地类转换的土地面积是否到达Cloud-Markov预测的转换上限。

      3) 若已经达到上限,标记该元胞该地类转换评价值为0,即丢弃该种转换方向,返回步骤1);若未达到上限,将该元胞记入该地类,丢弃该元胞其他类型转换评价值,返回步骤1)。

      重复上述过程,直到全体元胞均发生转换,最终获得土地利用变化分布预测。

    • 以中国科学院地理研究所数据中心(http://www.gscloud.cn/)提供的研究区2006年Landsat TM与2013年Landsat 8分辨率为30 m的影像为基础,利用ENVI 5.3和ArcGIS 10.2对影像进行预处理与解译,获得两期研究区土地利用图。其他数据包括杭州市30 m分辨率数字高程模型、路网、行政区划等矢量数据。

    • 实验用5台机器搭建Hadoop云环境。其软件环境为CentOS 7.1.1503系统,Java版本1.8.0_112,Hadoop版本2.7.3;硬件环境为CPU4核2.4 GHz,内存4 GB。串行算法测试环境与Hadoop云节点相同。

      Markov串行算法与Markov并行算法的运行结果及加速比如图 2图 3所示,单机单次CELUC与Cloud-CELUC运行结果及加速比如图 4图 5所示,其中横坐标表示元胞数量(以n为单位,1 n约为9 000 000个),纵坐标表示运行时间(以s为单位)。

      图  2  Cloud-Markov与串行Markov效率对比

      Figure 2.  Running Efficiency Between Cloud-Markov and Serial Markov Algorithms

      图  3  Cloud-Markov与串行Markov加速比

      Figure 3.  Speedup Ratio of Cloud-Markov and Serial Markov Algorithms

      图  4  Cloud-CELUC与串行CELUC效率对比

      Figure 4.  Running Efficiency Between Cloud-CELUC and Serial CELUC Algorithms

      图  5  Cloud-CELUC与串行CELUC加速比

      Figure 5.  Speedup Ratio of Cloud-CELUC and Serial CELUC Algorithms

      从实验结果可知,Cloud-Markov相对串行Markov模型加速比为3.27时增长趋向于平稳,为3.43时达到最高,Cloud-CELUC相对于串行CELUC为1.77时增长趋向于平稳,为1.86时达到最高。

      通过对实验过程中各节点的运行情况进行分析,Cloud-Markov运行效率提高的原因是Map-Reduce分散了两期数据元胞同位置匹配与统计数量的工作量; Cloud-CELUC在Map阶段运行速度较快,而在对Map阶段运算得到的输出进行Reduce合并时,数据迁移占用较长的运行时间,对运行效率有一定影响。未来将考虑设计数据分布优化策略提高Reduce阶段的合并效率。

    • 本文采用基于2006年与2013年土地利用数据计算面积转移概率矩阵,以2006年作为初始年份,水域设为不变,对研究区2013年土地利用数据进行模拟,面积转移矩阵与面积转移概率矩阵如表 1表 2所示。

      表 1  2006—2013年面积转移矩阵/km2

      Table 1.  Area Transfer Matrix in 2006—2013/km2

      地类 农用地 建设用地 自然保留地 总计
      农用地 1 282.95 409.71 95.67 1 788.33
      建设用地 210.94 1 381.69 12.52 1 605.15
      自然保留地 93.39 50.79 4 409.60 4 553.78
      总计 1 587.28 1 842.19 4 517.79 7 947.26

      表 2  2006—2013年转移概率矩阵/%

      Table 2.  Transition Probability Matrix in 2006—2013/%

      地类 农用地 建设用地 自然保留地
      农用地 71.74 22.91 5.35
      建设用地 13.14 86.08 0.78
      自然保留地 2.05 1.12 96.83

      利用模拟结果对C-MCE的各类因子权重参数进行精度评价校正,以提高模拟精度,实现准确预测。目前常用的精度评价方法有目视比较法、分维度检验法、像元对比检验法以及Kappa系数检验法,其中Kappa系数检验法是最常用的定量检验方法。对比计算2013年土地利用模拟数据与2013年土地利用实际数据,得到Kappa系数计算结果如表 3表 4表 5所示。

      表 3  自然保留地Kappa系数检验表/km2

      Table 3.  Kappa Coefficient of Natural Reserve/km2

      地类 自然保留地 非自然保留地 总计 准确度/% Kappa系数
      自然保留地 4 221.35 288.47 4 509.82 93.60 0.86
      非自然保留地 296.44 3 431.50 3 727.94 92.05
      总计 4 517.79 3 717.97

      表 4  建设用地Kappa系数检验表/km2

      Table 4.  Kappa Coefficient of Construction Land/km2

      地类 建设用地 非建设用地 总计 准确度/% Kappa系数
      建设用地 1 391.41 452.77 1 844.18 75.45 0.68
      非建设用地 450.78 5 942.80 6 393.58 92.95
      总计 1 842.19 6 395.57

      表 5  农用地Kappa系数检验表/km2

      Table 5.  Kappa Coefficient of Agricultural Land/km2

      地类 农用地 非农用地 总计 准确度/% Kappa系数
      农用地 1 152.64 427.17 1 579.81 72.96 0. 66
      非农用地 434.65 6 223.30 6 657.95 93.47
      总计 1 587.29 6 650.47

      自然保留地、建筑用地、农用地对应的Kappa系数分别为0.86、0.68和0.66,说明2013年模拟结果精度较高,可以推测对未来土地利用预测可信度较高。2020年土地利用变化预测实验基于2013年土地利用变化数据与表 6面积转移概率矩阵对2020年土地利用进行预测,预测结果如图 6所示。

      表 6  2013—2020年土地利用面积转移预测矩阵/km2

      Table 6.  Area Transfer Matrix in 2013—2020/km2

      地类 农用地 建设用地 自然保留地 总计
      农用地 1 133.35 361.94 84.52 1 579.81
      建设用地 242.35 1 587.45 14.38 1 844.18
      自然保留地 92.49 50.30 4 367.03 4 509.82
      总计 1 468.19 1 999.69 4 465.93 7 933.81

      图  6  2020年研究区土地利用变化预测图

      Figure 6.  Land Use Change Prediction of Research Area in 2020

      从预测结果分析可知,研究区内建设用地面积整体呈上升趋势且主要来源于农业用地转换,农用地面积呈下降趋势,自然保留地面积变化比例不大。建设用地总体增长面积较大,尤其各县市区中心因受路系与公交系统发达的影响,建设面积增速较快,其中西湖区、拱墅区、下城区、滨江区、萧山区建设用地的增长突出,其余县市区因受地形边界限制,建设用地保持基本稳定。

      研究区西北部与北部农用地与建设用地大片面积交织错落,但远离市县区中心且交通系统欠发达,农用地面积也基本不变。其余地区农用地受到城市影响辐射、路系与公交系统的影响,较大面积转换为建设用地。自然保留地由于多在西部与南部丘陵山区,受约束性因子坡度的影响较大,且交通欠发达,基本保持稳定。

    • 本文利用MapReduce编程模型将CA-Markov模型进行并行化设计,包括云环境下Markov算法以及云环境下土地利用变化综合评价算法。通过选取实验区并搭建实验环境,验证了方法的可靠性、正确性和运行效率,并模拟分析了2020年杭州市的土地利用变化预测结果。实验表明,将云计算技术与CA-Markov模型结合模拟预测土地利用变化从技术上是可行的。

      本文首次将云计算应用于土地利用变化预测,并且详细设计了云环境下的并行土地利用变化预测算法,但依然存在一些不足,主要有:土地利用变化除了自然条件的限制,还受到政治、经济、人口等复杂因素的综合影响,由于时间关系以及数据源的限制,研究主要针对交通、地形、区位因素进行建模,未来研究中将考虑更多的社会经济、政策及人口等因素,如国家规定的耕地红线、生态红线对于城市扩张边界的限制,完善土地利用变化空间格局预测模型。

参考文献 (21)

目录

    /

    返回文章
    返回