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

区域尺度的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 选择你所做的区的矢量文件。Value field 选择你要栅格化的因子。Output raster dataset 栅格文件存放的位置 保存为tif格式即可。Cell assigement type 选择maxinum_area 即最大面积作为栅格的属性。（所有因子都要保持一致）Priority field 选择noneCellsize 选择100.即栅格的大小为100*100. （所有因子都要保持一致） 为了防止出错以及方便快捷 采用批处理进行栅格化操作。基本设置和上面一致。所有因子以及底图。其中pd以逻辑回归的结果确定有几个参与回归。转化成栅格图之后，检查一下各因子是否和底图保持一样的形状和大小。具体方法看边角以及中间缺少的部分是否完全重合。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.同样使用reclassify工具。结果如下：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）。Conversion 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距离。为回归方程的常量值。第三行顶格。为该方程中的变量的个数。第四行及以后。为该方程中各变量的系数。并在系数后面标注变量的编号。中间空格隔开。下一个变量另起一行。因子编号因子编号pd10pd21gc2swmd3dxzs4_dis_citis5_dis_towns6_dis_roads7_XZRK8_DJGNYCZ9地类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.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.txt文件准备allow.txt 文件比较简单。主要用来限定不同地类之间的是否能相互转换。用转换矩阵来表示。我们假定所有的地类之间都可以相互转换。转换矩阵为：所以，allow.txt文件设置为：同样，将文件拷贝到clue-s软件的安装目录下即可。1.3.3 demand.in1文件准备Demand.in1 存放各地类不同年份的需求，即数量。其中2004年和2009年已知。中间年份的需求，通过线性内插得到。首先，得到2004年和2009年的各地类的数量，（上面步骤中可以获得）。对2009年进行平差处理，使得2009年的总数量与2004年一致。平差方法采用比例平差。然后进行四舍五入，得到整数个栅格。然后线性内插出其他年份的各地类的需求。同样四舍五入成整数个。最后新建一个txt文档。命名为demand.in1.第一行为年份：设置为6.后面几行为上面的表格数据。需要转置后，直接复制过去即可。同样，将文件拷贝到clue-s软件的安装目录下即可。1.3.4 因子文件asc准备将上面的因子txt文件重命名为模型适用的文件。（记得备份）具体命名如下表：命名因子sc1gr0.fil.ascpd1sc1gr1.fil.ascpd2sc1gr2.fil.ascgcsc1gr3.fil.ascswmdsc1gr4.fil.ascdxzssc1gr5.fil.asc_dis_citissc1gr6.fil.asc_dis_townssc1gr7.fil.asc_dis_roadssc1gr8.fil.asc2004rksc1gr8.1.asc2005rksc1gr8.2.asc2006rksc1gr8.3.asc2007rksc1gr8.4.asc2008rksc1gr8.5.asc2009rksc1gr9.fil.asc2004DJGNYCZsc1gr9.1.asc2005DJGNYCZsc1gr9.2.asc2006DJGNYCZsc1gr9.3.asc2007DJGNYCZsc1gr9.4.asc2008DJGNYCZsc1gr9.5.asc2009DJGNYCZregionbase.ascregionbase.txtcov_all.0.asccov_all.txt命名可以稍有不同，但是编号必须要从0开始。将重命名后的文件复制到安装目录。1.3.5 main文件准备打开main.1的文件，按照下面进行设置参数。地类个数6区域个数1最大因子个数8（回归方程中最大因子数，从alloc文件可以得到）总因子个数10列数303行数259单个栅格面积（公顷）1X坐标38555428.5Y坐标3358891.71土地利用类型序号0 1 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.其他设置成上面的值即可。保存一下。完成了所有文件的准备。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.