区域尺度的clues模型模拟.docx
《区域尺度的clues模型模拟.docx》由会员分享，可在线阅读，更多相关《区域尺度的clues模型模拟.docx（19页珍藏版）》请在冰豆网上搜索。
区域尺度的clues模型模拟
区域尺度的clue-s模型模拟
任务：
区域尺度三个区的clue-s模型模拟采用以前数据精度与市域做对比传统logistic回归
1华容区
1.1logistic回归
参考logistic回归操作。
1.2底图及因子准备
1.2.1底图准备：
底图主要为模拟底图，即模拟区域。
将底图以及各因子转化为栅格图。
在arcgis中进行操作。
1.2.1.1将矢量图转化栅格图
Conversiontools——toraster——polygontoraster（右键batch）（注意pd因子的选择）
参数设置为：
（下图为单独处理设置）
其中，Inputfeatures选择你所做的区的矢量文件。
Valuefield选择你要栅格化的因子。
Outputrasterdataset栅格文件存放的位置保存为tif格式即可。
Cellassigementtype选择maxinum_area即最大面积作为栅格的属性。
（所有因子都要保持一致）
Priorityfield选择none
Cellsize选择100.即栅格的大小为100*100.（所有因子都要保持一致）
为了防止出错以及方便快捷采用批处理进行栅格化操作。
基本设置和上面一致。
所有因子以及底图。
其中pd以逻辑回归的结果确定有几个参与回归。
转化成栅格图之后，检查一下各因子是否和底图保持一样的形状和大小。
具体方法看边角以及中间缺少的部分是否完全重合。
1.2.1.2cov_all文件准备
栅格图象进行重分类处理。
打开base1（底图）属性表，可以发现LYFSDM和Value值不对应。
需要进行保持一致。
同时记下不同地类栅格的数量（count）。
3Danalysttools——rasterreclassly——reclassify
只需改reclassification。
方法为LYFSDM减去1作为新的值。
依次填到newvalues下面。
保存为tif格式。
即可。
重分类之后的结果。
1.2.1.3底图文件准备
底图需要所有的模拟栅格值都为0.同样使用reclassify工具。
结果如下：
1.2.1.42009年现状图准备。
用上面制作的底图base来裁剪2009年的现状图。
得到华容区的09年现状图。
得到各地的2009年的数量。
首先将09年转化为栅格（polygontoraster）。
然后进行重分类（reclassify）。
最后进行裁剪。
Spatialanalysttools——extraction——extractbymask
裁剪结果：
记录下2009年各地类的数值。
后面需要用。
1.2.1.5因子文件准备。
将因子以及底图base、cov_all由栅格图转化为ascii编码文件（txt）。
Conversiontools——fromraster——rastertoascii（右键batch）
（最后一行23不要空不然会报错）
结果如下图
1.3clue-s模拟
1.3.1alloc1.reg文件准备
alloc1.reg文件主要存放logistic回归结果。
1.3.1.1文件结构：
第一行为地类编码。
例如耕地为0.林地为1，建设用地、水体、养殖水面、其他用地依次为2，3，4，5，6.
第二行一个tab距离。
为回归方程的常量值。
第三行顶格。
为该方程中的变量的个数。
第四行及以后。
为该方程中各变量的系数。
并在系数后面标注变量的编号。
中间空格隔开。
下一个变量另起一行。
因子编号
因子
编号
pd1
0
pd2
1
gc
2
swmd
3
dxzs
4
_dis_citis
5
_dis_towns
6
_dis_roads
7
_XZRK
8
_DJGNYCZ
9
地类0实例：
0
0.984
8
-0.434
0
-0.360
1
-0.608
3
-1.839
4
-2.026
6
0.863
7
0.654
8
0.314
9
所有地类的结果：
0
0.984
8
-0.434
0
-0.360
1
-0.608
3
-1.839
4
-2.026
6
0.863
7
0.654
8
0.314
9
1
-2.642
4
-0.863
0
-0.560
1
5.049
2
-1.518
8
2
-1.396
4
-0.423
0
-0.199
1
-0.676
6
-0.244
9
3
-1.365
8
-1.311
0
-0.936
1
1.837
3
-35.603
2
-1.720
4
9.201
6
-2.323
7
-1.248
8
4
-2.438
4
1.291
0
1.004
1
1.626
4
-0.495
8
5
-2.482
8
-0.430
0
-0.223
1
0.820
4
6.611
6
-2.120
5
-1.203
7
-0.428
8
-0.468
9
1.3.1.2创建文件
新建一个txt文档。
命名为alloc1.reg
将上面的表格复制粘贴进去。
保存即可。
如图：
将文件拷贝到clue-s软件的安装目录下即可。
1.3.2allow.txt文件准备
allow.txt文件比较简单。
主要用来限定不同地类之间的是否能相互转换。
用转换矩阵来表示。
我们假定所有的地类之间都可以相互转换。
转换矩阵为：
所以，allow.txt文件设置为：
同样，将文件拷贝到clue-s软件的安装目录下即可。
1.3.3demand.in1文件准备
Demand.in1存放各地类不同年份的需求，即数量。
其中2004年和2009年已知。
中间年份的需求，通过线性内插得到。
首先，得到2004年和2009年的各地类的数量，（上面步骤中可以获得）。
对2009年进行平差处理，使得2009年的总数量与2004年一致。
平差方法采用比例平差。
然后进行四舍五入，得到整数个栅格。
然后线性内插出其他年份的各地类的需求。
同样四舍五入成整数个。
最后新建一个txt文档。
命名为demand.in1.
第一行为年份：
设置为6.后面几行为上面的表格数据。
需要转置后，直接复制过去即可。
同样，将文件拷贝到clue-s软件的安装目录下即可。
1.3.4因子文件asc准备
将上面的因子txt文件重命名为模型适用的文件。
（记得备份）
具体命名如下表：
命名
因子
sc1gr0.fil.asc
pd1
sc1gr1.fil.asc
pd2
sc1gr2.fil.asc
gc
sc1gr3.fil.asc
swmd
sc1gr4.fil.asc
dxzs
sc1gr5.fil.asc
_dis_citis
sc1gr6.fil.asc
_dis_towns
sc1gr7.fil.asc
_dis_roads
sc1gr8.fil.asc
2004rk
sc1gr8.1.asc
2005rk
sc1gr8.2.asc
2006rk
sc1gr8.3.asc
2007rk
sc1gr8.4.asc
2008rk
sc1gr8.5.asc
2009rk
sc1gr9.fil.asc
2004DJGNYCZ
sc1gr9.1.asc
2005DJGNYCZ
sc1gr9.2.asc
2006DJGNYCZ
sc1gr9.3.asc
2007DJGNYCZ
sc1gr9.4.asc
2008DJGNYCZ
sc1gr9.5.asc
2009DJGNYCZ
regionbase.asc
regionbase.txt
cov_all.0.asc
cov_all.txt
命名可以稍有不同，但是编号必须要从0开始。
将重命名后的文件复制到安装目录。
1.3.5main文件准备
打开main.1的文件，按照下面进行设置参数。
地类个数
6
区域个数
1
最大因子个数
8（回归方程中最大因子数，从alloc文件可以得到）
总因子个数
10
列数
303
行数
259
单个栅格面积（公顷）
1
X坐标
38555428.5
Y坐标
3358891.71
土地利用类型序号
012345
转换弹性系数
0.60.8110.60.4
迭代变量系数
00.31
模拟的起始年份
20042009
动态变化驱动因子数字和编码
289
输出文件选择
3
特定区域回归选择
0
土地利用初试值
15
邻域选择计算
0
区域特定优先值
0
可选迭代变量参数
0.05
黄色部分可以根据实际情况进行调整。
行数列数X坐标Y坐标在因子文件中可以找到。
转换弹性系数。
为各地类的稳定系数。
地类越稳定，值越大，必须小于等于1.通过修改这个系数，可以得到不同的模拟结果。
迭代变量系数。
一般设置为00.31也可以设置为00.3518.
其他设置成上面的值即可。
保存一下。
完成了所有文件的准备。
1.3.6软件操作。
双击打开clues.exe
单击regionbase.asc和demand.in1.即选择了底图和需求文件。
单击editinput菜单下的mainparaments、regressionresults、changematrix。
检查上面准备的文件是否正确。
不对的可以进行编辑，然后保存。
Mode中选择第一个。
Filecheck中选择fliecheckon进行文件检查。
设置好了之后，点击run-dynaclues键进行运行。
1.4模拟结果
1.5结果分析
5、错误汇总：
1、error:
regressioncannotbecalculatedduetolargevalueincell0,200forlandcover3
解决方法：
查看因子中是否有未归一化的因子。
某个因子文件的值较大。
与逻辑回归中使用的不一致。
2、栅格转化的时候，栅格没有值。
解决方法：
转化的要素文件是否选择了某个值。
升级会员 