CN109671003B - 一种集成gcam和ca的全球土地利用和土地覆被序列空间化方法 - Google Patents
一种集成gcam和ca的全球土地利用和土地覆被序列空间化方法 Download PDFInfo
- Publication number
- CN109671003B CN109671003B CN201811569219.8A CN201811569219A CN109671003B CN 109671003 B CN109671003 B CN 109671003B CN 201811569219 A CN201811569219 A CN 201811569219A CN 109671003 B CN109671003 B CN 109671003B
- Authority
- CN
- China
- Prior art keywords
- land
- global
- land utilization
- spatialization
- utilization
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000007704 transition Effects 0.000 claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims abstract description 21
- 230000008859 change Effects 0.000 claims abstract description 20
- 238000013528 artificial neural network Methods 0.000 claims abstract description 15
- 238000011161 development Methods 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims description 20
- 238000004088 simulation Methods 0.000 claims description 11
- 238000011156 evaluation Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 239000002689 soil Substances 0.000 claims description 8
- 210000004027 cell Anatomy 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 claims description 4
- 238000001556 precipitation Methods 0.000 claims description 3
- 230000004913 activation Effects 0.000 claims description 2
- 239000004927 clay Substances 0.000 claims description 2
- 210000002569 neuron Anatomy 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims description 2
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 3
- 238000011160 research Methods 0.000 abstract description 3
- 238000003062 neural network model Methods 0.000 description 4
- 238000012549 training Methods 0.000 description 4
- 238000005192 partition Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000001413 cellular effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000007477 logistic regression Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Development Economics (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Biomedical Technology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Educational Administration (AREA)
- Biophysics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种集成GCAM和CA的全球土地利用和土地覆被序列空间化方法，将GCAM模型与地理CA相结合，实现全球283个子区域、2010年到2100年的全球土地利用序列空间化数据产品，空间分辨率为1km，时间分辨率为5年；综合神经网络初始概率、邻域影响、各子区域的地类转变权重矩阵和规划约束条件，得到各子区域地类的综合转变概率，对各子区域的土地利用和土地覆被进行空间化。本发明主要应用于未来不同发展情景下的全球土地利用变化与覆被建模，具备实现生产高分辨率全球LUCC产品的能力，推动了全球尺度下高分辨率土地利用模型构建的研究。
Description
技术领域
本发明涉及元胞自动机土地利用建模领域，尤其是一种集成GCAM和CA的全球土地利用和土地覆被序列空间化方法。
背景技术
全球土地利用和土地覆盖变化(LUCC)对于理解人类活动与全球环境变化(如生物多样性和碳循环)之间的相互作用至关重要。已有的研究表明，构建LUCC模型已经成为一种揭示全球LUCC变化的驱动机制、研究过去以及未来年份的全球LUCC时空变化的重要工具，全球尺度的LUCC模型已经成为针对全球性问题进行综合评价的关键部分。尽管在许多情况下需要全球尺度的土地利用与土地变化建模，但由于其复杂性，目前很少有全球尺度的LUCC模型存在。
目前，已有的针对全球LUCC构建的模型主要有全球环境综合评估模型(IMAGE)，空间分辨率为0.5°(赤道约55公里)；LUSs模型和CLUMondo模型，空间分辨率为5弧分(赤道约10公里)。这些模型空间分辨率粗糙，很大程度地损失了土地利用和土地覆被的空间细节，对全球不同区域的土地利用模式产生了严重的扭曲。
元胞自动机(CA)以其简单而小规模的组织单元，可以获取其转换规则而用作模拟大尺度复杂的空间现象，这一特性适用于高空间分辨率的土地利用和土地覆盖的模拟。近几十年来，CA已经被用于模拟各种城市与区域尺度的LUCC变化。例如，区域尺度的GEOMOD模型、中国尺度的DynaCLUE模型、以及中国-全球尺度的Landshift模型等。CA还可以用作模拟全球评估模型预测下的地类需求约束的全球尺度的LUCC变化。例如，FLUS整个了CA与IMAGE模型来模拟全球尺度的LUCC变化结果，其空间分辨率达到了1km。空间分辨率的极大提高，大大增强了土地利用与土地变化的模拟效果。
然而，已有的全球规模LUCC模型，忽略了土地利用和土地覆盖的空间异质性。土地利用和土地覆盖的空间异质性在全球影响评估中起着重要作用。例如，一个农业生态区域内部可能具有相同的气候带，湿度条件，土壤类型和地貌，但由于该区域又属于不同国家，而具有不同的土地利用和土地覆盖空间特征。同样，一个国家的邻域范围内具有相同的地缘政治和社会经济条件，但由于气候、水分、土壤或地貌的变化，其土地利用和土地覆盖仍具有显着的空间异质性。在全球评估模型中，GCAM模型较好地反映出了土地利用和土地覆盖的空间异质性，GCAM模型在温度带，湿度与政治经济分区的基础上，综合18个农业生态区(AEZ)与32个经济政治区，形成全球283个子区域，来预测这283个子区域的土地利用和土地覆盖需求，综合神经网络初始概率、邻域影响、各子区域的地类转变权重矩阵和规划约束条件，得到各子区域地类的综合转变概率，对各子区域的土地利用和土地覆被进行序列空间化。本发明通过整合GCAM和CA来构建全球尺度的土地利用序列空间化GCAM-CA模型。
发明内容
本发明所要解决的技术问题在于，提供一种集成GCAM和CA的全球土地利用和土地覆被序列空间化方法，具备实现生产高分辨率全球LUCC产品的能力，推动了全球尺度下高分辨率土地利用模型构建的研究。
为解决上述技术问题，本发明提供一种集成GCAM和CA的全球土地利用和土地覆被序列空间化方法，包括如下步骤：
(1)综合考虑农业生态区划和地缘政治经济区划，将全球陆地划分为283个子区域，获取283个子区域不同地类的土地利用需求；
(2)结合全球土地利用变化的驱动因素，获取12个空间变量，使用神经网络获取初始地类转变概率，分区域获取不同土地利用类型之间转换规则；
(3)计算全球283个子区域的地类转变权重矩阵；
(4)计算全球283个子区域综合转变概率；
(5)基于(4)步骤的综合转变概率进行空间化；
(6)针对全球未来三种发展情景，采取5年的步长，对全球283个区域的空间化结果进行拼接，得到全球土地利用空间化序列结果；
(7)针对2010年的模拟结果进行精度评价与模拟评价，验证模型有效性与实用性。
优选的，步骤(2)中，12个空间变量包括：2000年、2010年的两期GLC30全球土地利用数据、全球GDP数据、人口数据、到主要城市点的距离、到道路的距离、到铁路的距离、DEM数据、坡度数据、土壤含沙比例、土壤含泥比例、土壤含黏土比例、pH数据、温度与降水12类栅格数据，空间分辨率均为1km。
优选的，步骤(2)中，使用神经网络获取初始地类转变概率，分区域获取不同土地利用类型之间转换规则具体为：针对每一个子区域，采用多层感知机，预测十类用地类型变化的概率，输入层为每一个子区域2000年的12个空间变量和土地利用类，隐藏层包含15个神经元，输出层为该子区域2010年的十类土地利用类型，激活函数为logistics函数，学习率为0.01。
优选的，步骤(3)中，计算全球283个子区域的地类转变权重矩阵具体为：地类转换权重矩阵的值的大小表示当前土地利用类型转换为特定的目标土地利用类型的概率大小，在土地利用空间化过程中，283个世界区域的地类转换权重矩阵各不同，所有类型土地利用的转换权重矩阵随时间变化，特定其余的权重矩阵每隔5年进行一次计算一次。
优选的，步骤(4)中，计算全球283个子区域综合转变概率具体为：综合考虑由ANN计算的初始概率，邻域效应，转换权重矩阵和条件约束，计算特定区域某一网格将被特定土地利用类型覆盖的组合概率，计算方法见公式(1)：
其中，
是t时刻元胞i从原始地类转变为目标地类k的综合概率；
是t时刻元胞i从原始地类转变为目标地类k的初始概率，由神经网络计算得到；
是邻域影响，这里采用7×7摩尔邻域，在每一次迭代过程中动态计算；conc→k是原始地类c转变为目标地类k的地类转换权重，Pcons是地类转变过程中的限制条件。
优选的，步骤(5)中，基于(4)步骤的综合转变概率进行空间化具体为：针对每一个世界子区域，根据GCAM提供的未来不同发展情景下的面积变化预测序列，综合神经网络提供的初始土地类型转变概率、随模拟年份动态变化的土地类型转变权重矩阵，以及邻域对未来的土地利用影响，进行不同情景下的未来全球土地利用与土地覆被空间化。
优选的，步骤(6)中，对全球283个区域的空间化结果进行拼接，得到全球土地利用空间化序列结果具体为：2015年至2100年的共18期全球土地利用与土地覆被空间化结果图，时间分辨率为5年，空间分辨率为1km。
本发明的有益效果为：本发明主要应用于未来不同发展情景下的全球土地利用变化与覆被建模，具备实现生产高分辨率全球LUCC产品的能力，推动了全球尺度下高分辨率土地利用模型构建的研究。
附图说明
图1为本发明的基本框架结构示意图。
图2为本发明GCAM-CA模型空间化基本流程示意图。
图3为本发明GCAM-CA模型中神经网络模型结构示意图。
具体实施方式
一种集成GCAM和CA的全球土地利用和土地覆被序列空间化方法，包括如下步骤：
(1)综合考虑农业生态区划和地缘政治经济区划，将全球陆地划分为283个子区域，获取283个子区域不同地类的土地利用需求；
(2)结合全球土地利用变化的驱动因素，获取12个空间变量，使用神经网络获取初始地类转变概率，分区域获取不同土地利用类型之间转换规则；
(3)计算全球283个子区域的地类转变权重矩阵；
(4)计算全球283个子区域综合转变概率；
(5)基于(4)步骤的综合转变概率进行空间化；
(6)针对全球未来三种发展情景，采取5年的步长，对全球283个区域的空间化结果进行拼接，得到全球土地利用空间化序列结果；
(7)针对2010年的模拟结果进行精度评价与模拟评价，验证模型有效性与实用性。
本发明提供一种集成GCAM和CA的全球土地利用和土地覆被序列空间化方法，本方法包含两个部分，如图1所示。第一部分为全球283个区域各地类面积变化序列预测，第二部分为全球283个区域土地利用序列空间化，此部分流程图如图2所示。
一种集成GCAM和CA的全球土地利用和土地覆被序列空间化方法基本框架，如图1所示，包括两部分：
(1)获取不同情景下的283个世界子区域的土地利用变化面积预测的过程。
(11)GCAM模型输出283个全球区域分区，综合了32个经济政治区与18个农业生态区，这一分区结果考虑了社会经济发展因素和温度降水等自然因素。
(12)针对GLC30模型、GCAM模型、IMAGE模型、FLUS模型这四个模型土地利用分类体系的区别，提出一套地类合并规则，将后三类模型的地类统一为GLC30数据的地类标准，以进行产品输出与精度评价。
(13)通过GCAM模型的农业土地利用模块(AgLU)，根据气候，社会经济条件，历史土地利用和技术进步等驱动因素，以5年为间隔输出283个世界子区域的各地类面积变化序列数据预测。
(14)将GCAM模型的输出全球283个区域各地类面积变化结果进行校正和转换，得到适用于GCAM-CA模型的面积变化趋势，作为GCAM-CA模型的面积变化控制模块的输入。
(15)根据GCAM模型动态地输出序列GDP和人口数据，利用多元逻辑回归，依据GDP与人口数据修正城市用地需求预测，得到城市用地面积序列预测。
(21)进行不同情景下、全球283个区域土地利用序列空间化的过程。
(21)对全球283个子区域12个空间变量和两期GLC30土地利用数据进行采集与处理，空间分辨率均重采样1km。
(22)充分考虑283个区域地类组成的差异性特点，提出一种新的按原有地类比例进行采样点分层抽取的方法，获取用于神经网络训练的数据集。
(23)分别针对283个区域训练神经网络模型，并综合神经网络初始概率、邻域影响、各子区域的地类转变权重矩阵和规划约束条件，得到各子区域地类的综合转变概率，获取283个不同区域的CA空间化规则。
(24)根据各子区域地类的综合转变概率，在GCAM预测的三种未来发展情景序列面积约束下，对全球各子区域的土地利用和土地覆盖进行空间化，空间分辨率为1km，时间分辨率为5年。
(25)对不同情景下的283个世界子区域的各地类空间化结果进行投影转换、拼接、存储，从而得到全球土地利用空间化序列产品数据。
第一部分是全球283个区域各地类面积变化序列预测，具体实施步骤包含以下内容：
GCAM模型输出包含不同政策影响情景下未来土地利用和土地变化的面积数据，该数据是每一个区域单元的统计数据。由于GCAM、GLC30地类分类体系差异，需结合真实的全球30米分辨率的数据，确定2010-2100每隔5年的地类面积变化量，从而落实在全球空间上，以进行空间化处理与分析。另外，与本模型对比的IMAGE和FLUS模型的地类体系也存在差异，因此提出一种不同模型之间地类合并的判定标准，如表1所示。
表1 GLC30、GCAM、IMAGE和FLUS之间的地类合并规则
GCAM模型生成序列土地利用需求预测的核心在于农业土地利用模块。GCAM模型下的农业土地利用模块决定了农产品的供需，该模型与气候政策链接，共同决定了生物能源的供需以及生产生物能源土地面积的变化。农业土地利用模型包括一个需求市场和一个供应市场，两个市场之间通过价格、税收和供给来达到市场均衡；需求端有人口和经济决定，供应端依据气候，技术和生物能源的需求来决定；供应端与需求端通过国际市场这一手段进行贸易。基于这样的市场机制，该模块模拟输出未来各种土地利用类型的数量变化、所有农产品价格、产量以及生物能源的可供应量等。基于此，GCAM模型下的农业土地利用模块根据气候，社会经济条件，历史土地利用和技术进步等驱动因素，以5年为间隔预测283个世界子区域的土地利用需求。
第二部分是对全球283个区域分别进行土地利用和土地覆被序列空间化，拼接形成不同情景下、全球土地利用序列产品。具体实施步骤包括以下内容：
(1)全球283个子区域12个空间变量和两期土地利用数据的采集与处理；
2000年和2010年的GLC30全球土地利用原始数据有10种地类，空间率为30米。在本模型中，最初的30米分辨率重采样为1km，以创建CA模型的初始土地利用数据。两期真实土地利用数据用于校准GCAM-CA模型。同时，全球土地利用和土地覆盖变化还受到人为和自然驱动因素的影响，在本方法中，12个空间变量(见表2)被选定为土地利用变化的驱动因子。在模型中，所有空间数据和变量的空间分辨率重采样至1km，投影方式为World_Goode_Homolosine_Land。
表2GCAM-CA模型中用于获取转换规则的空间变量
(2)神经网络训练；
针对283个世界区域、10种土地利用类型，使用分层随机抽样方法提取训练样本，训练样本占整个全球数据集的百分之五，共有6,792,283个样本，用于训练和验证识别土地利用和土地覆盖类型的神经网络。采用三层神经网络模型(如图3所示)，输入层为每一个子区域2000年的12个空间变量和土地利用类，输出层为该子区域2010年的十类土地利用类型。使用训练好的283个对应的神经网络模型计算每个子区域各地类转变的初始概率。
(3)地类转换权重矩阵计算；
不同土地类型的转换权重矩阵表明两种土地利用类型之间转换的可能性，也是影响土地利用空间化的另一个重要因素。地类转换权重矩阵的值的大小表示当前土地利用类型转换为特定的目标土地利用类型的概率大小。在土地利用空间化过程中，283个世界区域的地类转换权重矩阵各不同。所有类型土地利用的转换权重矩阵随时间变化，特定其余的权重矩阵每隔5年进行一次计算一次，如表3所示，显示了编号为223的子区域2000年到2010年地类转变成本。
表3 2000年到2010年地类转变成本矩阵(223区域)
(4)综合转变概率计算
综合考虑由ANN计算的初始概率，邻域效应，转换权重矩阵和条件约束，计算特定区域某一网格将被特定土地利用类型覆盖的组合概率，计算方法见公式(1)：
其中，
是t时刻元胞i从原始地类转变为目标地类k的综合概率；
是t时刻元胞i从原始地类转变为目标地类k的初始概率，由神经网络计算得到；
是邻域影响，这里采用7×7摩尔邻域，在每一次迭代过程中动态计算；conc→k是原始地类c转变为目标地类k的地类转换权重，Pcons是地类转变过程中的限制条件。
(5)GCAM模型与CA的集成；
GCAM-CA模型中，GCAM模型根据气候，社会经济条件，历史土地利用和技术进步，以5年为间隔预测283个世界子区域的土地利用需求。在迭代过程中，GCAM输出的283个世界区域的土地利用需求预测被用作土地利用的约束条件，推动CA模型中土地利用和土地覆盖类型的空间化过程。同时，使用GCAM模型输出的GDP和人口数据，采用多元回归的方法，校准GCAM模型输出中未改变的城市土地面积的预测需求。2010年至2100年的GDP和人口数据系列也被用作CA模型的输入空间变量。此外，CA模型中的全球土地利用和土地覆盖数据与GCAM模型相结合，确保GCAM模型预测的土地利用面积需求与观测到的土地利用面积之间的一致性。
本发明针对GLC30数据、GCAM模型、IMAGE模型、FLUS模型这四个模型土地利用分类体系的区别，提出一套地类合并规则，将后三类模型的地类统一为GLC30数据的地类标准，以进行产品输出与精度评价。
Claims (2)
1.一种集成GCAM和CA的全球土地利用和土地覆被序列空间化方法，其特征在于，包括如下步骤：
(1)综合考虑农业生态区划和地缘政治经济区划，将全球陆地划分为283个子区域，获取283个子区域不同地类的土地利用需求；
(2)结合全球土地利用变化的驱动因素，获取12个空间变量，使用神经网络获取初始地类转变概率，分区域获取不同土地利用类型之间转换规则；针对每一个子区域，采用多层感知机，预测十类用地类型变化的概率，输入层为每一个子区域2000年的12个空间变量和土地利用类，隐藏层包含15个神经元，输出层为该子区域2010年的十类土地利用类型，激活函数为logistics函数，学习率为0.01；
(3)计算全球283个子区域的地类转变权重矩阵；地类转换权重矩阵的值的大小表示当前土地利用类型转换为特定的目标土地利用类型的概率大小，在土地利用空间化过程中，283个世界区域的地类转换权重矩阵各不同，所有类型土地利用的转换权重矩阵随时间变化，特定其余的权重矩阵每隔5年进行一次计算一次；
(4)计算全球283个子区域综合转变概率；具体为：综合考虑由ANN计算的初始概率，邻域效应，转换权重矩阵和条件约束，计算特定区域某一网格将被特定土地利用类型覆盖的组合概率，计算方法见公式(1)：
其中，
是t时刻元胞i从原始地类转变为目标地类k的综合概率；
是t时刻元胞i从原始地类转变为目标地类k的初始概率，由神经网络计算得到；
是邻域影响，这里采用7×7摩尔邻域，在每一次迭代过程中动态计算；conc→k是原始地类c转变为目标地类k的地类转换权重，Pcons是地类转变过程中的限制条件；
(5)基于(4)步骤的综合转变概率进行空间化；具体为：针对每一个世界子区域，根据GCAM提供的未来不同发展情景下的面积变化预测序列，综合神经网络提供的初始土地类型转变概率、随模拟年份动态变化的土地类型转变权重矩阵，以及邻域对未来的土地利用影响，进行不同情景下的未来全球土地利用与土地覆被空间化；
(6)针对全球未来三种发展情景，采取5年的步长，对全球283个区域的空间化结果进行拼接，得到全球土地利用空间化序列结果；2015年至2100年的共18期全球土地利用与土地覆被空间化结果图，时间分辨率为5年，空间分辨率为1km；
(7)针对2010年的模拟结果进行精度评价与模拟评价，验证模型有效性与实用性。
2.如权利要求1所述的集成GCAM和CA的全球土地利用和土地覆被序列空间化方法，其特征在于，步骤(2)中，12个空间变量包括：2000年、2010年的两期GLC30全球土地利用数据、全球GDP数据、人口数据、到主要城市点的距离、到道路的距离、到铁路的距离、DEM数据、坡度数据、土壤含沙比例、土壤含泥比例、土壤含黏土比例、pH数据、温度与降水12类栅格数据，空间分辨率均为1km。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811569219.8A CN109671003B (zh) | 2018-12-21 | 2018-12-21 | 一种集成gcam和ca的全球土地利用和土地覆被序列空间化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811569219.8A CN109671003B (zh) | 2018-12-21 | 2018-12-21 | 一种集成gcam和ca的全球土地利用和土地覆被序列空间化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109671003A CN109671003A (zh) | 2019-04-23 |
| CN109671003B true CN109671003B (zh) | 2023-04-25 |
Family
ID=66145679
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811569219.8A Active CN109671003B (zh) | 2018-12-21 | 2018-12-21 | 一种集成gcam和ca的全球土地利用和土地覆被序列空间化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109671003B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110263111B (zh) * | 2019-05-31 | 2024-03-08 | 南京林业大学 | 基于先验知识的土地利用/覆被信息时空监测方法 |
| CN110390129A (zh) * | 2019-06-11 | 2019-10-29 | 同济大学 | 基于GeoSOS-FLUS的土地利用政策有效性的定量评价方法 |
| CN112700025A (zh) * | 2019-10-22 | 2021-04-23 | 北京京东振世信息技术有限公司 | 覆盖区域的分配方法及装置、存储介质、电子设备 |
| CN112231952B (zh) * | 2020-09-30 | 2024-05-14 | 东南大学 | 一种基于wrf土地覆被数据置换的风暴潮模拟方法 |
| CN112487724B (zh) * | 2020-12-11 | 2022-09-23 | 重庆邮电大学 | 基于分区和改进cnn-ca模型的城市动态扩张模拟方法 |
| CN112801487A (zh) * | 2021-01-21 | 2021-05-14 | 中国科学院新疆生态与地理研究所 | 土地利用类型转换植物功能类型的方法、终端及存储介质 |
| CN113221443A (zh) * | 2021-04-20 | 2021-08-06 | 深圳市数字城市工程研究中心 | 一种集成SDGs和美丽指标的土地利用建模与模拟方法 |
| CN116934026B (zh) * | 2023-07-18 | 2024-02-27 | 中国科学院地理科学与资源研究所 | 一种ssp-rcp情景下城市建成区土地利用数据集制作方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105447235A (zh) * | 2015-11-12 | 2016-03-30 | 中山大学 | 一种土地未来利用情景动态模拟方法 |
| CN107807909A (zh) * | 2017-10-12 | 2018-03-16 | 中国地质大学（北京） | 一种区域土地用途调控政策模拟与效应分析方法 |
-
2018
- 2018-12-21 CN CN201811569219.8A patent/CN109671003B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105447235A (zh) * | 2015-11-12 | 2016-03-30 | 中山大学 | 一种土地未来利用情景动态模拟方法 |
| CN107807909A (zh) * | 2017-10-12 | 2018-03-16 | 中国地质大学（北京） | 一种区域土地用途调控政策模拟与效应分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| GCAM-CA;曹敏;《OpenGMS》;20180722;全文 * |
| 基于元胞自动机降尺度方法的1km分辨率全球土地利用数据集(2010～2100)(英文);李雪草等;《Science Bulletin》;20161115(第21期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109671003A (zh) | 2019-04-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109671003B (zh) | 一种集成gcam和ca的全球土地利用和土地覆被序列空间化方法 | |
| Huang et al. | Analysis of the future land cover change in Beijing using CA–Markov chain model | |
| Wang et al. | Dynamic simulation of land use change and assessment of carbon storage based on climate change scenarios at the city level: A case study of Bortala, China | |
| Dinda et al. | An integrated simulation approach to the assessment of urban growth pattern and loss in urban green space in Kolkata, India: A GIS-based analysis | |
| Domingo et al. | Effect of zoning plans on urban land-use change: A multi-scenario simulation for supporting sustainable urban growth | |
| CN109166309A (zh) | 一种面向复杂城市交通网络的缺失流量数据恢复方法 | |
| CN109830102A (zh) | 一种面向复杂城市交通网络的短时交通流量预测方法 | |
| Yagoub et al. | Prediction of land cover change using Markov and cellular automata models: case of Al-Ain, UAE, 1992-2030 | |
| Saganeiti et al. | Modeling urban sprinkling with cellular automata | |
| Zhang et al. | Spatial suitability and multi-scenarios for land use: Simulation and policy insights from the production-living-ecological perspective | |
| Khwarahm et al. | Past and future prediction of land cover land use change based on earth observation data by the CA–Markov model: a case study from Duhok governorate, Iraq | |
| CN107194504A (zh) | 土地使用状态的预测方法、装置和系统 | |
| Zhang et al. | Scenario analysis of the relationship among ecosystem service values—A case study of Yinchuan Plain in northwestern China | |
| CN102646164A (zh) | 一种结合空间滤波的土地利用变化建模方法及其系统 | |
| CN112016247A (zh) | 基于土地利用变化的高精度未来水文过程耦合模型及构建方法 | |
| Zhou et al. | A novel spatio-temporal cellular automata model coupling partitioning with CNN-LSTM to urban land change simulation | |
| Setturu et al. | Modeling landscape dynamics of policy interventions in Karnataka State, India | |
| Tian et al. | Multi-scale solar radiation and photovoltaic power forecasting with machine learning algorithms in urban environment: A state-of-the-art review | |
| Cao et al. | Analysis of spatiotemporal changes in cultural heritage protected cities and their influencing factors: Evidence from China | |
| Ajeeb et al. | The prediction of urban growth trends and patterns using spatio-temporal CA-MC model in Seremban Basin | |
| Li et al. | Research on land‐use evolution and ecosystem services value response in mountainous counties based on the SD‐PLUS model | |
| Messina et al. | Dynamic spatial simulation modeling of the population—environment matrix in the Ecuadorian Amazon | |
| Parchianloo et al. | Integrated CA model and remote sensing approach for simulating the future development of a city | |
| Pan et al. | Impact of land use change on regional carbon sink capacity: Evidence from Sanmenxia, China | |
| Deal et al. | A spatially explicit urban simulation model: land use evolution and impact assessment model (LEAM) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |