区域尺度的clue-s模型模拟

《区域尺度的clue-s模型模拟》由会员分享，可在线阅读，更多相关《区域尺度的clue-s模型模拟（15页珍藏版）》请在金锄头文库上搜索。

1、区域尺度的 clue-s 模型模拟任务：区域尺度 三个区的 clue-s 模型模拟 采用以前数据 精度与市域做对比 传统 logistic回归 1 华容区1.1logistic 回归参考 logistic 回归操作。1.2 底图及因子准备1.2.1 底图准备：底图主要为模拟底图，即模拟区域。将底图以及各因子转化为栅格图。在 arcgis 中进行操作。1.2.1.1 将矢量图转化栅格图Conversion toolsto rasterpolygon to raster（右键 batch） （注意 pd 因子的选择）参数设置为：（下图为单独处理设置）其中，Input features 选择你所做的

2、区的矢量文件。Value field 选择你要栅格化的因子。Output raster dataset 栅格文件存放的位置 保存为 tif 格式即可。 Cell assigement type 选择 maxinum_area 即最大面积作为栅格的属性。 （所有因子都要保持一致）Priority field 选择 noneCellsize 选择 100.即栅格的大小为 100*100. （所有因子都要保持一致）为了防止出错以及方便快捷 采用批处理进行栅格化操作。基本设置和上面一致。所有因子以及底图。其中 pd 以逻辑回归的结果确定有几个参与回归。转化成栅格图之后，检查一下各因子是否和底图保持一样

3、的形状和大小。具体方法看边角以及中间缺少的部分是否完全重合。1.2.1.2 cov_all 文件准备栅格图象进行重分类处理。打开 base1（底图）属性表，可以发现 LYFSDM 和 Value 值不对应。需要进行保持一致。同时记下不同地类栅格的数量（count ） 。3D analyst toolsraster reclasslyreclassify 只需改 reclassification。方法为 LYFSDM 减去 1 作为新的值。依次填到 new values 下面。保存为 tif 格式。即可。重分类之后的结果。1.2.1.3 底图文件准备底图需要所有的模拟栅格值都为 0.同样使用 r

4、eclassify 工具。结果如下：1.2.1.4 2009 年现状图准备。用上面制作的底图 base 来裁剪 2009 年的现状图。得到华容区的 09 年现状图。得到各地的 2009 年的数量。首先将 09 年转化为栅格（polygon to raster） 。然后进行重分类（reclassify） 。最后进行裁剪。Spatial analyst toolsextractionextract by mask裁剪结果：记录下 2009 年各地类的数值。后面需要用。1.2.1.5 因子文件准备。将因子以及底图 base、cov_all 由栅格图转化为 ascii 编码文件（txt） 。Conve

5、rsion tools from rasterraster to ascii （右键 batch）（最后一行 23 不要空 不然会报错）结果如下图1.3clue-s 模拟1.3.1 alloc1.reg 文件准备alloc1.reg 文件主要存放 logistic 回归结果。1.3.1.1 文件结构：第一行为地类编码。例如耕地为 0.林地为 1，建设用地、水体、养殖水面、其他用地依次为 2，3 ，4 ，5，6.第二行一个 tab 距离。为回归方程的常量值。第三行顶格。为该方程中的变量的个数。第四行及以后。为该方程中各变量的系数。并在系数后面标注变量的编号。中间空格隔开。下一个变量另起一行。因子

6、编号因子 编号pd1 0pd2 1gc 2swmd 3dxzs 4_dis_citis 5_dis_towns 6_dis_roads 7_XZRK 8_DJGNYCZ 9地类 0 实例：00.984 8-0.434 0-0.360 1-0.608 3-1.839 4-2.026 60.863 70.654 80.314 9所有地类的结果：00.984 8-0.434 0-0.360 1-0.608 3-1.839 4-2.026 60.863 70.654 80.314 91-2.642 4-0.863 0-0.560 15.049 2-1.518 82-1.396 4-0.423 0-0.

7、199 1-0.676 6-0.244 93-1.365 8-1.311 0-0.936 11.837 3-35.603 2-1.720 49.201 6-2.323 7-1.248 84-2.438 41.291 01.004 11.626 4-0.495 85-2.482 8-0.430 0-0.223 10.820 46.611 6-2.120 5-1.203 7-0.428 8-0.468 91.3.1.2 创建文件新建一个 txt 文档。命名为 alloc1.reg 将上面的表格复制粘贴进去。保存即可。如图：将文件拷贝到 clue-s 软件的安装目录下即可。1.3.2 allow.t

8、xt 文件准备allow.txt 文件比较简单。主要用来限定不同地类之间的是否能相互转换。用转换矩阵来表示。我们假定所有的地类之间都可以相互转换。转换矩阵为：地 类 0 1 2 3 4 50 1 1 1 1 1 11 1 1 1 1 1 12 1 1 1 1 1 13 1 1 1 1 1 14 1 1 1 1 1 15 1 1 1 1 1 1所以，allow.txt 文件设置为：同样，将文件拷贝到 clue-s 软件的安装目录下即可。1.3.3 demand.in1 文件准备Demand.in1 存放各地类不同年份的需求，即数量。其中 2004 年和 2009 年已知。中间年份的需求，通过线性

9、内插得到。首先，得到 2004 年和 2009 年的各地类的数量， （上面步骤中可以获得） 。对 2009 年进行平差处理，使得 2009 年的总数量与 2004 年一致。平差方法采用比例平差。然后进行四舍五入，得到整数个栅格。地 类 2004 2009 2009 2009 20090 25200 22018 22048.09 22048.09 220481 529 982 983.3421 983.3421 9832 5212 7759 7769.604 7769.604 77703 6667 5039 5045.887 5045.887 50464 6596 9029 9041.34 90

10、41.34 90415 1957 1271 1272.737 1272.737 127346161 46098 46161 46161 46161然后线性内插出其他年份的各地类的需求。同样四舍五入成整数个。地 类 2004 2005 2006 2007 2008 20090 25200 24569.6 23939.2 23308.8 22678.4 220481 529 619.8 710.6 801.4 892.2 9832 5212 5723.6 6235.2 6746.8 7258.4 77703 6667 6342.8 6018.6 5694.4 5370.2 50464 6596 7

11、085 7574 8063 8552 90415 1957 1820.2 1683.4 1546.6 1409.8 1273all 46161 46161 46161 46161 46161 46161地 类 2004 2005 2006 2007 2008 20090 25200 24570 23939 23309 22678 220481 529 620 711 801 892 9832 5212 5723 6235 6747 7259 77703 6667 6343 6019 5694 5370 50464 6596 7085 7574 8063 8552 90415 1957 1820

12、 1683 1547 1410 1273all 46161 46161 46161 46161 46161 46161最后新建一个 txt 文档。命名为 demand.in1.第一行为年份：设置为 6.后面几行为上面的表格数据。需要转置后，直接复制过去即可。同样，将文件拷贝到 clue-s 软件的安装目录下即可。1.3.4 因子文件 asc 准备将上面的因子 txt 文件重命名为模型适用的文件。 （记得备份）具体命名如下表：命名 因子sc1gr0.fil.asc pd1sc1gr1.fil.asc pd2sc1gr2.fil.asc gcsc1gr3.fil.asc swmdsc1gr4.fi

13、l.asc dxzssc1gr5.fil.asc _dis_citissc1gr6.fil.asc _dis_townssc1gr7.fil.asc _dis_roadssc1gr8.fil.asc 2004rksc1gr8.1.asc 2005rksc1gr8.2.asc 2006rksc1gr8.3.asc 2007rksc1gr8.4.asc 2008rksc1gr8.5.asc 2009rksc1gr9.fil.asc 2004DJGNYCZsc1gr9.1.asc 2005DJGNYCZsc1gr9.2.asc 2006DJGNYCZsc1gr9.3.asc 2007DJGNYCZs

14、c1gr9.4.asc 2008DJGNYCZsc1gr9.5.asc 2009DJGNYCZregionbase.asc regionbase.txtcov_all.0.asc cov_all.txt命名可以稍有不同，但是编号必须要从 0 开始。将重命名后的文件复制到安装目录。1.3.5 main 文件准备打开 main.1 的文件，按照下面进行设置参数。地类个数 6区域个数 1最大因子个数8（回归方程中最大因子数，从 alloc 文件可以得到）总因子个数 10列数 303行数 259单个栅格面积（公顷） 1X 坐标 38555428.5Y 坐标 3358891.71土地利用类型序号 0 1

15、 2 3 4 5转换弹性系数 0.6 0.8 1 1 0.6 0.4迭代变量系数 0 0.3 1模拟的起始年份 2004 2009动态变化驱动因子数字和编码 2 8 9输出文件选择 3特定区域回归选择 0土地利用初试值 1 5邻域选择计算 0区域特定优先值 0可选迭代变量参数 0.05黄色部分可以根据实际情况进行调整。行数 列数 X 坐标 Y 坐标 在因子文件中可以找到。转换弹性系数。为各地类的稳定系数。地类越稳定，值越大，必须小于等于 1.通过修改这个系数，可以得到不同的模拟结果。迭代变量系数。一般设置为 0 0.3 1 也可以设置为 0 0.35 18.其他设置成上面的值即可。保存一下。完

16、成了所有文件的准备。1.3.6 软件操作。双击打开 clues.exe 单击 regionbase.asc 和 demand.in1.即选择了底图和需求文件。单击 edit input 菜单下的 main paraments 、regression results、change matrix。检查上面准备的文件是否正确。不对的可以进行编辑，然后保存。Mode 中选择第一个。File check 中选择 flie check on 进行文件检查。设置好了之后，点击 run-dynaclues 键进行运行。1.4 模拟结果1.5 结果分析5、错误汇总：1、 error: regression cannot be calculated due to large value in cell 0,200 for land cover 3解决方法：查看因子中是否有未归一化的因子。某个因子文件的值较大。与逻辑回归中使用的不