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

《区域尺度的clue-s模型模拟》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、区域尺度的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选择你所做的区的矢量文件

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

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

4、下：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——fromraste

5、r——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距离。为回归方程的常量值。第三行顶格。为该方程中的变量的个数。第四行及以后。为该方程中各变量的系数。并在系数后面标注变量的编号。中间空格隔开。下一个变量另起一行。因子编号因子编号pd10pd21gc2swmd3

6、dxzs4_dis_citis5_dis_towns6_dis_roads7_XZRK8_DJGNYCZ9地类0实例：00.9848-0.4340-0.3601-0.6083-1.8394-2.02660.86370.65480.3149所有地类的结果：00.9848-0.4340-0.3601-0.6083-1.8394-2.02660.86370.65480.31491-2.6424-0.8630-0.56015.0492-1.51882-1.3964-0.4230-0.1991-0.6766-0.24493-1.3658-1.3110-0.93611

7、.8373-35.6032-1.72049.2016-2.3237-1.24884-2.43841.29101.00411.6264-0.49585-2.4828-0.4300-0.22310.82046.6116-2.1205-1.2037-0.4288-0.46891.3.1.2创建文件新建一个txt文档。命名为alloc1.reg将上面的表格复制粘贴进去。保存即可。如图：将文件拷贝到clue-s软件的安装目录下即可。1.3.2allow.txt文件准备allow.txt文件比较简单。主要用来限定不同地类之间的是否能相互转换。用转换矩阵来表示。我们假定

8、所有的地类之间都可以相互转换。转换矩阵为：所以，allow.txt