Dinamica_土地利用/覆被变化模拟

1.获取转移矩阵

##################30*30分辨率#######################

1980-2000

Cells per category (initial landscape):
1       8273849
2       13003062
3       39709577
4       10146092
5       30551912
6       660601
7       424942
8       340740
9       5467

Cells per transition:
 From \ To |         1          2          3          4          5          6          7          8          9 
---------------------------------------------------------------------------------------------------------------
         1 |      XXXX       1701       3001        607       6816       7065      26824       3370         -- 
         2 |      2339       XXXX      11601       7815      22004       4497      12916       1092         -- 
         3 |      3153      49721       XXXX      11119      55340        323       2576         86         -- 
         4 |       970      10919       4715       XXXX      31589       1165        544         49         -- 
         5 |     18298      89834      38134     144987       XXXX       5597        347       2290         -- 
         6 |     18826      13337       2857        932       7812       XXXX        220      67827         -- 
         7 |      1959         61          2          5        151          4       XXXX          3         -- 
         8 |      6133      69261         63        131       2802      11179        247       XXXX         -- 
         9 |        --         --         41         --         --         --         --         --       XXXX 

Single Step Transition Matrix:
 From \ To |         1          2          3          4          5          6          7          8          9 
---------------------------------------------------------------------------------------------------------------
         1 |      XXXX  0.0002056  0.0003627  0.0000734  0.0008238  0.0008539  0.0032420  0.0004073         -- 
         2 | 0.0001799       XXXX  0.0008922  0.0006010  0.0016922  0.0003458  0.0009933  0.0000840         -- 
         3 | 0.0000794  0.0012521       XXXX  0.0002800  0.0013936  0.0000081  0.0000649  0.0000022         -- 
         4 | 0.0000956  0.0010762  0.0004647       XXXX  0.0031134  0.0001148  0.0000536  0.0000048         -- 
         5 | 0.0005989  0.0029404  0.0012482  0.0047456       XXXX  0.0001832  0.0000114  0.0000750         -- 
         6 | 0.0284983  0.0201892  0.0043248  0.0014108  0.0118256       XXXX  0.0003330  0.1026747         -- 
         7 | 0.0046100  0.0001435  0.0000047  0.0000118  0.0003553  0.0000094       XXXX  0.0000071         -- 
         8 | 0.0179991  0.2032664  0.0001849  0.0003845  0.0082233  0.0328080  0.0007249       XXXX         -- 
         9 |        --         --  0.0074995         --         --         --         --         --       XXXX 

Multi Step Transition Matrix:
 From \ To |         1          2          3          4          5          6          7          8          9 
---------------------------------------------------------------------------------------------------------------
         1 |      XXXX  0.0000078  0.0000181  0.0000036  0.0000411  0.0000463  0.0001629  0.0000207         -- 
         2 | 0.0000086       XXXX  0.0000447  0.0000300  0.0000850  0.0000188  0.0000499  0.0000037         -- 
         3 | 0.0000039  0.0000627       XXXX  0.0000139  0.0000700  0.0000004  0.0000032  0.0000001         -- 
         4 | 0.0000047  0.0000538  0.0000232       XXXX  0.0001567  0.0000063  0.0000027         --         -- 
         5 | 0.0000300  0.0001474  0.0000627  0.0002389       XXXX  0.0000099  0.0000004  0.0000037         -- 
         6 | 0.0015001  0.0004206  0.0002352  0.0000744  0.0006202       XXXX  0.0000133  0.0064902         -- 
         7 | 0.0002317  0.0000071  0.0000002  0.0000005  0.0000178  0.0000004       XXXX  0.0000003         -- 
         8 | 0.0010096  0.0117324  0.0000003  0.0000161  0.0004536  0.0020711  0.0000341       XXXX         -- 
         9 |        --         --  0.0003769         --         --         --         --         --       XXXX 

"DetermineTransitionMatrix" ran successfully (elapsed 35 s).



2000-2015
Cells per category (initial landscape):
1       8276143
2       13175632
3       39647673
4       10261737
5       30378939
6       578620
7       466431
8       325641
9       5426

Cells per transition:
 From \ To |         1          2          3          4          5          6          7          8          9 
---------------------------------------------------------------------------------------------------------------
         1 |      XXXX       1475       1333       5700      10002      26664     113786       1720       1659 
         2 |       143       XXXX      13782      52308     171112      45239      40793      11402       1020 
         3 |      1073      10232       XXXX      55747       3784      19530      66139        486       1339 
         4 |       910       5189       7001       XXXX       3275       4924      12087         23        293 
         5 |      3323      23478      15149      14550       XXXX       9020      21771        729       1985 
         6 |       420        157         14        133        681       XXXX        320       2812         -- 
         7 |         5          9          6        213          8        364       XXXX         81         21 
         8 |       632         --        185        758        490      60599        794       XXXX         -- 
         9 |        --         --         --         --         --         --         --         --       XXXX 

Single Step Transition Matrix:
 From \ To |         1          2          3          4          5          6          7          8          9 
---------------------------------------------------------------------------------------------------------------
         1 |      XXXX  0.0001782  0.0001611  0.0006887  0.0012085  0.0032218  0.0137487  0.0002078  0.0002005 
         2 | 0.0000109       XXXX  0.0010460  0.0039701  0.0129870  0.0034335  0.0030961  0.0008654  0.0000774 
         3 | 0.0000271  0.0002581       XXXX  0.0014061  0.0000954  0.0004926  0.0016682  0.0000123  0.0000338 
         4 | 0.0000887  0.0005057  0.0006822       XXXX  0.0003191  0.0004798  0.0011779  0.0000022  0.0000286 
         5 | 0.0001094  0.0007728  0.0004987  0.0004790       XXXX  0.0002969  0.0007166  0.0000240  0.0000653 
         6 | 0.0007259  0.0002713  0.0000242  0.0002299  0.0011769       XXXX  0.0005530  0.0048598         -- 
         7 | 0.0000107  0.0000193  0.0000129  0.0004567  0.0000172  0.0007804       XXXX  0.0001737  0.0000450 
         8 | 0.0019408         --  0.0005681  0.0023277  0.0015047  0.1860914  0.0024383       XXXX         -- 
         9 |        --         --         --         --         --         --         --         --       XXXX 

Unable to calculate the Multi-Step Transition Matrix.

"DetermineTransitionMatrix" ran successfully (elapsed 34 s).

2.制作Cube集合数据集

• slope / altitude 提取

ArcGIS教程之DEM应用（坡度坡向、提等高线）-百度经验

• 人口密度，首先使用一个大于研究区的矩形裁剪全球人口密度数据集，用于投影转换（直接使用全球数据集投影转换将花费大量的时间！），然后投影转换（输出分辨率30m）——用制作好的slope栅格mask提取
• d_river 计算欧式距离，注意输出分辨率，以及投影坐标系！

• d_建设用地/d_林地，首先将1980年土地利用重分类，将建设用地以外的土地利用类型变为Nodata。然后计算欧式距离（动态变量不用单独制作，直接在选到初始土地利用中建设用地/林地的距离）
• d_road（公路+铁路）（全国地理信息资源目录服务系统） 地理坐标系转换——shp裁剪shp（ArcToolbox——>分析工具——提取——>clip）——欧氏距离

3.确定证据权重范围

根据转移矩阵，确定x--y类型。

对于动态变量，layer应为“distance_to_X"

4.计算填充每一个范围内的证据权重

5.分析环境因子是否独立

证据权重法唯一假设是其参与模拟土地利用变化的环境因子是相互独立的，因此使用假设检验去检验因子间独立性。

对于第一个检验指标：卡方独立性检验，结果需要与分布临界值表对比，确定是否独立。

Person Chi-Square适用情况是N>40 person ei>1

Cintinuity correction 适用:2X2

Likelihood Ratio 适用N<40

Fisher‘s Exact Test’适用N<40 2X2

对于非2X2做修正使用Cramer’s V检验，且（0,1）越接近1关系越强

6.初步运行土地利用变化模拟

需要逐个设定每一土地利用类型转移的 平均斑块大小/斑块尺寸变异/斑块聚合度，初步随意设定，目的是为了理解下一步（模型验证）

7.使用exponential decay函数验证模拟结果

结果只描述已更改的单元格.

通过在地图查看器中打开它来可视化结果相似地图。使用PseudoColor, Limits to Actual, and Histogram Equalize.

红色和黄色区域显示高到中等的空间拟合，而蓝色表示较差的拟合。

8.另一种验证方法

结果并不好：相似度太低，即使11*11网格搜索，也才40%作用相似性

9.理解调整patch参数，模拟结果更加接近自然（与第6步相比，将参数信息调整，由于上一步校准发现结果不满意！）

10.增加expander参数，继续提升模拟结果！

该步骤的结构化语言也是最终用来预测的语言！