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

区域尺度的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选择你所做的区的矢量文件。Valuefield选择你要栅格化的因子。Outputrasterdataset栅格文件存放的位置保存为tif格式即可。Cellassigementtype选择maxinum_area即最大面积作为栅格的属性。（所有因子都要保持一致）Priorityfield选择noneCellsize选择100.即栅格的大小为100*100.（所有因子都要保持一致）为了防止出错以及方便快捷采用批处理进行栅格化操作。基本设置和上面一致。所有因子以及底图。其中pd以逻辑回归的结果确定有几个参与回归。转化成栅格图之后，检查一下各因子是否和底图保持一样的形状和大小。具体方法看边角以及中间缺少的部分是否完全重合。1.2.1.2cov_all文件准备栅格图象进行重分类处理。打开base1（底图）属性表，可以发现LYFSDM和Value值不对应。需要进行保持一致。同时记下不同地类栅格的数量（count）。3Danalysttools——rasterreclassly——reclassify只需改reclassification。方法为LYFSDM减去1作为新的值。依次填到newvalues下面。保存为tif格式。即可。重分类之后的结果。1.2.1.3底图文件准备底图需要所有的模拟栅格值都为0.同样使用reclassify工具。结果如下：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——fromraster——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距离。为回归方程的常量值。第三行顶格。为该方程中的变量的个数。第四行及以后。为该方程中各变量的系数。并在系数后面标注变量的编号。中间空格隔开。下一个变量另起一行。因子编号因子编号pd10pd21gc2swmd3dxzs4_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.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文件比较简单。主要用来限定不同地类之间的是否能相互转换。用转换矩阵来表示。我们假定所有的地类之间都可以相互转换。转换矩阵为：所以，allow.txt文件设置为：地类012345011111111111112111111311111141111115111111同样，将文件拷贝到clue-s软件的安装目录下即可。1.3.3demand.in1文件准备Demand.in1存放各地类不同年份的需求，即数量。其中2004年和2009年已知。中间年份的需求，通过线性内插得到。首先，得到2004年和2009年的各地类的数量，（上面步骤中可以获得）。对2009年进行平差处理，使得2009年的总数量与2004年一致。平差方法采用比例平差。然后进行四舍五入，得到整数个栅格。然后线性内插出其他年份的各地类的需求。同样四舍五入成整数个。最后新建一个txt文档。命名为demand.in1.地类200420092009200920090252002201822048.0922048.09220481529982983.3421983.34219832521277597769.6047769.60477703666750395045.8875045.88750464659690299041.349041.3490415195712711272.7371272.73712734616146098461614616146161地类20042005200620072008200902520024569.623939.223308.822678.4220481529619.8710.6801.4892.2983252125723.66235.26746.87258.47770366676342.86018.65694.45370.250464659670857574806385529041519571820.21683.41546.61409.81273all461614616146161461614616146161地类200420052006200720082009025200245702393923309226782204815296207118018929832521257236235674772597770366676343601956945370504646596708575748063855290415195718201683154714101273all461614616146161461614616146161第一行为年份：设置为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.5main文件准备打开main.1的文件，按照下面进行设置参数。地类个数6区域个数1最大因子个数8（回归方程中最大因子数，从alloc文件可以得到）总因子个数10列数303行数259单个栅格面积（公顷）1X坐标38555428.5Y坐标3358891.71土地利用类型序号012345转换弹性系数0.60.8110.60.4迭代变量系数00.31模拟的起始年份20042009动态变化驱动因子数字和编码289输出文件选择3特定区域回归选择0土地利用初试值15邻域选择计算0区域特定优先值0可选迭代变量参数0.05黄色部分可以根据实际情况进行调整。行数列数X坐标Y坐标在因子文件中可以找到。转换弹性系数。为各地类的稳定系数。地类越稳定，值越大，必须小于等于1.通过修改这个系数，可以得到不同的模拟结果。迭代变量系数。一般设置为00.31也可以设置为00.3518.其他设置成上面的值即可。保存一下。完成了所有文件的准备。1.3.6软件操作。双击打开clues.exe单击regionbase.asc和demand.in1.即选择了底图和需求文件。单击editinput菜单下的mainparaments、regressionresults、changematrix。检查上面准备的文件是否正确。不对的可以进行编辑，然后保存。Mode中选择第一个。Filecheck中选择fliecheckon进行文件检查。设置好了之后，点击run-dynaclues键进行运行。1.4模拟结果1.5结果分析5、错误汇总：1、error:regressioncannotbecalculatedduetolargevalueincell0,200forlandcover3解决方法：查看因子中是否有未归一化的因子。某个因子文件的值较大。与逻辑回归中使用的不一致。2、栅格转化的时候，栅格没有值。解决方法：转化的要素文件是否选择了某个值。