CN108563733B - 一种基于贝叶斯推理的土地利用数据同化方法 - Google Patents
一种基于贝叶斯推理的土地利用数据同化方法 Download PDFInfo
- Publication number
- CN108563733B CN108563733B CN201810311341.9A CN201810311341A CN108563733B CN 108563733 B CN108563733 B CN 108563733B CN 201810311341 A CN201810311341 A CN 201810311341A CN 108563733 B CN108563733 B CN 108563733B
- Authority
- CN
- China
- Prior art keywords
- land use
- assimilation
- use pattern
- model
- observation
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000004088 simulation Methods 0.000 claims abstract description 17
- 238000011160 research Methods 0.000 claims abstract description 14
- 230000001413 cellular effect Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 238000013316 zoning Methods 0.000 claims description 2
- 230000004927 fusion Effects 0.000 abstract description 6
- 239000002689 soil Substances 0.000 description 3
- 238000013480 data collection Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000007500 overflow downdraw method Methods 0.000 description 2
- 238000013439 planning Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/04—Inference or reasoning models
-
- 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/16—Real estate
- G06Q50/165—Land development
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Software Systems (AREA)
- Economics (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Artificial Intelligence (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (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
本发明涉及一种基于贝叶斯推理的土地利用数据同化方法，该方法的核心是将离散的类型分布变量对应到连续的类型分布密度，并将连续的类型分布密度作为中间变量，进而计算类型变量的后验分布概率。本发明将基于共轭先验的贝叶斯推理思想引入到土地利用动态模型中，通过融合观测数据，对离散的、多维的模型模拟的土地利用类型进行同化，提高土地利用动态模型的模拟预测精度。本发明在模型方面，首次利用多维土地利用/覆被变化动态模型对多种土地利用类型进行同化研究。本发明在方法上，首次引入一种同化策略，对离散的类型分布状态变量进行更新。
Description
摘要
本发明发展了一种基于贝叶斯推理的土地利用数据同化方法，该方法的核心是将离散的类型分布变量对应到连续的类型分布密度，并将连续的类型分布密度作为中间变量，进而计算类型变量的后验分布概率。本发明将基于共轭先验的贝叶斯推理思想引入到土地利用动态模型中，通过融合观测数据，对离散的、多维的模型模拟的土地利用类型进行同化，提高土地利用动态模型的模拟预测精度。
技术领域
本发明涉及地理信息科学技术领域，具体来说是一种基于贝叶斯推理的土地利用数据同化方法。
背景技术
土地利用变化是人类活动对地球表层***影响最直接的表现方式，通过改变地球表层***的下垫面，影响了地球表层***各大圈层的循环过程。因此，准确模拟和预测过去、现在、未来地球表层真实的土地利用变化格局和过程，可为土地利用规划、政策决策提供科学依据，也可服务于区域可持续发展。
模型和观测是土地利用变化研究的两种基本手段。模型有助于深入了解土地利用变化成因，重建过去、预测未来发展趋势、评估环境影响以及支持土地利用规划决策分析。目前，元胞自动机（Cellular Automata，CA）已成为土地利用变化模拟最常用的工具之一。但元胞自动机模型结构、静态参数以及边界条件、初始场引起的误差和不确定性会随着模型的运行不断的传递和积累；而观测能够快速准确地获取各种时空尺度上观测时刻的地表真实土地利用状况。但无法了解无观测时刻的状况，不能提供时间上连续的土地利用变化情况。如何将观测信息和模型有机融合，来量化和减少土地利用模拟和预测中的各种不确定性，获取时空连续的土地利用信息，是目前土地利用变化研究中面临的一重大挑战。
数据同化（Data Assimilation），也称“模型-数据融合”（mode-data fusion），是一种基于模型的动力学框架，通过融合不同来源、不同分辨率的直接和间接观测，在考虑模型和各种观测误差的基础上不断地依靠观测而自动调整模型轨迹、优化模型状态，并可同步优化参数、定量评价和减少不确定性的一种方法。数据同化为融合观测和土地利用模型提供了理论框架，为该领域的研究方法注入了新鲜的血液。
然而，无论国际还是国内，土地利用数据同化研究刚刚起步，理论和方法还存在不完善的地方。不同于水文、气象、和陆地过程模型，土地利用动态模型的输出为离散的类型变量（土地利用类型），缺少能够直接用于同化的状态变量。纵观国内外土地利用数据同化相关研究可知，（1）当前研究中使用的均为二维的土地利用动态模型，没有开展对多维的土地利用动态模型的同化。（2）绝大多数研究是基于粒子滤波和集合卡尔曼滤波算法对连续的模型参数进行优化，缺少对离散的、多维的类型分布状态变量（土地利用类型）的更新。
鉴于此，有必要瞄准国际上土地利用研究中模型-观测融合方法。
发明内容
针对上述的两个主要问题，本发明从方法论的角度提出了一种基于贝叶斯推理的土地利用数据同化策略，为提高土地利用动态模型的模拟预测精度，以及离散型变量的数据同化方法研究提供了较好的解决方案。
土地利用动态模型的输出是一个离散的土地利用类型。基于贝叶斯推理方法对土地利用进行同化的核心思想是将离散的类型分布变量对应到连续的类型分布密度，并将连续的类型分布密度作为中间变量，来计算类型变量的后验分布概率。
具体土地利用数据同化实施步骤如下：
S1：运行土地利用动态模型，得到同化年份的模拟值α，同时提取该同化年份的观测值ｃ。此处的模拟值和观测值均指各土地利用类型对应的数目。
S2：运用基于共轭先验的贝叶斯推理公式，得到同化值，即各观测点中各土地利用类型对应的分布密度值。其中，基于共轭先验的贝叶斯推理公式如下：
式中，是同化后的土地利用类型的分布密度值；是土地利用类型的观测数目；代表模型中模拟的土地利用类型数目；是所有土地利用类型的观测数目；是模拟的所有土地利用类型数目。
S3：根据同化值，对当前的模拟结果，即模拟的土地利用类型进行修正，获得该年份同化后的土地利用类型。
本发明的有益效果是：
本发明利用研究中模型-观测融合方法这一崭新的发展方向，从方法论的角度探索土地利用研究的新途径，在多维土地利用动态模型的基础上，融合观测数据，优化多维离散的类型分布状态变量，以其更准确地给出土地利用的空间和时间分布特征。
本发明在模型方面，首次利用多维土地利用/覆被变化动态模型对多种土地利用类型进行同化研究。
本发明在方法上，首次引入一种同化策略，对离散的类型分布状态变量进行更新。
附图说明
图1为本发明实施土地利用数据同化方法的流程图。
图2为对土地利用动态模型的模拟结果进行修正的流程图。
图3为本发明的同化结果与模拟结果以及2001年和2009年真实土地利用空间分布对比示意图
具体实施方式
下面结合附图和具体实施例对本发明做进一步阐述说明。
实施例1
本发明中的研究对象为位于黑河流域中游地区的张掖绿洲甘州区。甘州区是张掖绿洲的行政中心，总面积约3665平方公里。甘州区是2008年张掖绿洲启动的湿地保护政策的重点实施区域，并且该区域的城镇用地扩张及扩耕在张掖绿洲较为典型。本研究中所采用土地利用动态模型为多维逻辑回归马尔科夫元胞自动机模型（MLRMCA），采用的测试数据为：1992-1999年的甘州区土地利用历史变化数据集作为MLRMCA模型的校准数据集，以确定模型参数；1999年-2009年的甘州区土地利用数据集作为同化的观测及结果验证数据，每2年进行一次同化。以上土地利用数据集由TM/ETM+影像解译而成，采用相同的分类***。所有区域的数据均采用同坐标系、同尺度的栅格数据文件。
图1为本发明的方法流程图，主要包括以下几个步骤：
第1步：利用历史变化数据集对土地利用动态模型进行校准，获得模型的参数值。
第2步：将研究区划分为若干个正方形网格，每个网格包含ｋ×ｋ个像元（元胞）。选择其中的若干个正方形网格作为观测点。
第3步：运行土地利用动态模型，当模拟到观测年份时，得到土地利用类型转换概率，以及模拟值α，同时提取该同化年份的观测值ｃ。此处的模拟值和观测值均指各土地利用类型对应的数目。
第4步：运行基于共轭先验的贝叶斯推理公式，得到同化值。其中，基于共轭先验的贝叶斯推理公式如下：
式中，是同化后的土地利用类型的分布密度值；是土地利用类型的观测数目；代表模型中模拟的土地利用类型数目；是所有土地利用类型的观测数目；是模拟的所有土地利用类型数目。
第5步：根据同化值，对当前的模拟结果，即模拟的土地利用类型进行修正，获得该观测年份同化后的土地利用类型。
修正模拟结果的流程见图2，具体修正步骤如下：
S1：在观测点j，基于同化值和正方形方格的面积，计算得到各土地利用类型的变化数量，其计算公式如下：
S2：依次寻找变化数量的土地利用类型。
S3：将的土地利用类型的转换概率从小到大逐像元减少
个元胞（像元），同时将这些像元的土地利用类型赋值为0。
S4：将所有的土地利用类型按照步骤S3进行重新分配，最终得到观测点j中需重新分配的所有元胞。
S5：在待分配的元胞中，将剩余的土地利用类型，即的土地利用类型，依次增加个元胞。
S6：重复步骤S1-S5，直至完成所有观测点的修正。
第6步：把修正后的结果再输入到模型中，作为初始值，进行一下年份的模拟同化。当模拟到下一观测时刻，重复第3步～第5步，直到模拟同化年份结束。
图3为同化结果与模型模拟的结果对比。该图表明基于本发明的方法更新后的土地利用类型空间分布更加接近于真实情况，能够提高模型的模拟预测精度。
具体的同化与模拟精度结果如下：
（OA为总体精度，Kappa为Kappa系数，FM为性能指数）。
Claims (2)
1.一种基于贝叶斯推理的土地利用数据同化方法，其特征在于：所述同化方法包括以下步骤：
第1步：利用历史变化数据集对土地利用动态模型进行校准，获得模型的参数值；
第2步：将研究区划分为若干个正方形网格，每个网格包含ｋ×ｋ个像元；选择其中的若干个正方形网格作为观测点；
第3步：运行土地利用动态模型，当模拟到观测年份时，得到土地利用类型转换概率，以及模拟值α，同时提取该同化年份的观测值ｃ；
第4步：运行基于共轭先验的贝叶斯推理公式，得到同化值；其中，基于共轭先验的贝叶斯推理公式如下：
式中，是同化后的土地利用类型的分布密度值；是土地利用类型的观测数目；代表模型中模拟的土地利用类型数目；N是所有土地利用类型的观测数目；是模拟的所有土地利用类型数目；
第5步：根据同化值，对当前的模拟结果进行修正；
具体修正步骤如下：
S1：在观测点，基于同化值和正方形方格的面积，计算得到各土地利用类型的变化数量，其计算公式如下：
S2：依次寻找变化数量的土地利用类型；
S3：将的土地利用类型的转换概率从小到大逐像元减少
个元胞，即像元；同时将这些像元的土地利用类型赋值为0；
S4：将所有的土地利用类型按照步骤S3进行重新分配，最终得到观测点中需重新分配的所有元胞；
S5：在待分配的元胞中，将剩余的土地利用类型，即的土地利用类型，依次增加个元胞；
S6：重复步骤S1-S5，直至完成所有观测点的修正；
第6步：把修正后的结果再输入到模型中，作为初始值，进行一下年份的模拟同化；
当模拟到下一观测时刻，重复第3步～第5步，直到模拟同化年份结束。
2.如权利要求1所述的一种基于贝叶斯推理的土地利用数据同化方法，其特征在于：在第3步中的模拟值和观测值均指各土地利用类型对应的数目。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810311341.9A CN108563733B (zh) | 2018-04-09 | 2018-04-09 | 一种基于贝叶斯推理的土地利用数据同化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810311341.9A CN108563733B (zh) | 2018-04-09 | 2018-04-09 | 一种基于贝叶斯推理的土地利用数据同化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108563733A CN108563733A (zh) | 2018-09-21 |
| CN108563733B true CN108563733B (zh) | 2019-07-26 |
Family
ID=63534350
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810311341.9A Expired - Fee Related CN108563733B (zh) | 2018-04-09 | 2018-04-09 | 一种基于贝叶斯推理的土地利用数据同化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108563733B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113505948B (zh) * | 2021-09-13 | 2021-11-19 | 四川师范大学 | 基于贝叶斯网络的未来生态安全格局预测和优化方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102737155A (zh) * | 2011-04-12 | 2012-10-17 | 中国科学院寒区旱区环境与工程研究所 | 基于贝叶斯滤波的通用数据同化方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101814117A (zh) * | 2010-04-14 | 2010-08-25 | 北京师范大学 | 一种多源环境生态信息数据同化方法 |
| CN102034027A (zh) * | 2010-12-16 | 2011-04-27 | 南京大学 | 流域尺度土壤湿度遥感数据同化方法 |
-
2018
- 2018-04-09 CN CN201810311341.9A patent/CN108563733B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102737155A (zh) * | 2011-04-12 | 2012-10-17 | 中国科学院寒区旱区环境与工程研究所 | 基于贝叶斯滤波的通用数据同化方法 |
Non-Patent Citations (2)
| Title |
|---|
| Identifying a land use change cellular automaton by Bayesian data assimilation;Judith A.Verstegen et.al;《Environmental Modelling & Software》;20131220;第121-136页 * |
| 数据同化算法研究现状综述;马建文 等;《地球科学进展》;20120710;第27卷(第7期);第747-754页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108563733A (zh) | 2018-09-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liao et al. | Urban sprawl scenario simulations based on cellular automata and ordered weighted averaging ecological constraints | |
| Wijesekara et al. | Assessing the impact of future land-use changes on hydrological processes in the Elbow River watershed in southern Alberta, Canada | |
| Fan et al. | Temporal and spatial change detecting (1998–2003) and predicting of land use and land cover in Core corridor of Pearl River Delta (China) by using TM and ETM+ images | |
| Liu et al. | Rural land use spatial allocation in the semiarid loess hilly area in China: Using a Particle Swarm Optimization model equipped with multi-objective optimization techniques | |
| Li et al. | Pattern of spatial evolution of rural settlements in the Jizhou District of China during 1962–2030 | |
| Long et al. | Reconstruction of historical arable land use patterns using constrained cellular automata: A case study of Jiangsu, China | |
| CN112884226A (zh) | 基于多智能体算法的国土空间格局模拟规划方法及*** | |
| CN113902580A (zh) | 一种基于随机森林模型的历史耕地分布重建方法 | |
| CN109271683A (zh) | 一种面向日照约束的建筑群自动排布算法 | |
| CN103886220B (zh) | 基于bp网络和基尼系数定权重的用水数据离散化方法 | |
| CN105389469A (zh) | 一种暴雨洪水管理模型参数自动率定方法 | |
| Chen et al. | Simulating the optimal land-use pattern in the farming-pastoral transitional zone of Northern China | |
| CN109858180B (zh) | 一种区域作物耗水空间格局优化设计方法 | |
| CN108733952B (zh) | 一种基于序贯模拟的土壤含水量空间变异性三维表征方法 | |
| Xu et al. | Land-use change modeling with cellular automata using land natural evolution unit | |
| CN112528516A (zh) | 一种耦合土地利用类型和气候变化的流域水环境管理方法 | |
| CN113743659B (zh) | 一种基于成分法和马尔可夫元胞自动机的城市布局预测方法及应用 | |
| CN115758856A (zh) | 一种景观格局及气候变化对流域未来水质影响的研究方法 | |
| CN115952421A (zh) | 耦合生态过程模型与机器学习算法的高精度时空模拟方法 | |
| Wen et al. | Application of a loose coupling model for assessing the impact of land-cover changes on groundwater recharge in the Jinan spring area, China | |
| CN108563733B (zh) | 一种基于贝叶斯推理的土地利用数据同化方法 | |
| CN111859303B (zh) | 基于动态贝叶斯平均的土壤湿度融合方法及*** | |
| Zhao et al. | Examining land-use change trends in yucheng district, Ya’an city, China, using ANN-CA modeling | |
| CN106546958B (zh) | 一种优化的雷达数据同化方法 | |
| Schaldach et al. | Conceptual design and implementation of a model for the integrated simulation of large-scale land-use systems |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200714 Address after: 730000 No.318, Donggang West Road, Chengguan District, Lanzhou City, Gansu Province Patentee after: NORTHWEST INSTITUTE OF ECO-ENVIRONMENT AND RESOURCES, CHINESE ACADEMY OF SCIENCES Address before: 730000 Gansu city of Lanzhou province Donggang West Road No. 320 Patentee before: COLD AND ARID REGIONS ENVIRONMENTAL AND ENGINEERING Research Institute,CHINESE ACADEMY OF SCIENCES |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190726 |