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

区区域域尺尺度度的的clue-s模模型型模模拟拟 区域尺度的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 选择none Cellsize选择100.即栅格的⼤⼩为100*100. （所有因⼦都要保持⼀致） 为了防⽌出错以及⽅便快捷采⽤批处理进⾏栅格化操作。基本设置和上⾯⼀致。 所有因⼦以及底图。其中pd以逻辑回归的结果确定有⼏个 与回归。 转化成栅格图之后，检查⼀下各因⼦是否和底图保持⼀样的形状和⼤⼩。具体⽅法看边⾓以及中间缺少的部分是否完全重合。 1.2.1.2 cov all⽂件准备 栅格图象进⾏重分类处理。 打开base1 （底图）属性表，可以发现LY FSDM和Value值不对应。需要进⾏保持⼀致。同时记下不同地类栅格的数量 （count）。 3D analyst tools——raster reclassly——reclassify 只需改reclassification。⽅法为LY FSDM减去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距离。为回归⽅程的常量值。 第三⾏顶格。为该⽅程中的变量的个数。 第四⾏及以后。为该⽅程中各变量的系数。并在系数后⾯标注变量的编号。中间空格隔开。下⼀个变量另起⼀⾏。 因⼦编号 因⼦编号 pd1 0 pd2 1 gc 2 sw md 3 dxzs 4 dis citis 5 dis tow ns 6 dis roads 7 XZRK 8 DJGNYCZ 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 所有地类的结果： 0.984 8 -0.434 0 -0.360 1 -0.608 3 -1.839 4 -2.026 6 0.863 7 0.314 9 1 -2.642 4 -0.863 0