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

区域尺度的clue-s模型模拟任务：区域尺度 三个区的clue-s模型模拟 采用以前数据 精度与市域做对比 传统logistic回归 1华容区1.1logistic回归参考logistic回归操作。1.2底图及因子准备1.2.1底图准备：底图主要为模拟底图，即模拟区域。将底图以及各因子转化为栅格图。在arcgis中进行操作。1.2.1.1将矢量图转化栅格图Conversion tools——to raster——polygon 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 tools——raster reclassly——reclassify 只需改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 tools——extraction——extract by mask裁剪结果：记录下2009年各地类的数值。后面需要用。1.2.1.5 因子文件准备。将因子以及底图base、cov_all由栅格图转化为ascii编码文件（txt）。Conversion tools—— from raster——raster 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软件的安装目录下