arcgis教程-刘小平ppt课件.pptVIP

　五、基于逻辑回归的CA模型－核心代码变量 　 Public data(,) As Int32 ‘2001urban.txt 数据，以列行存储 Public dataFinal(,) As Int32 ‘2005urban.txt 数据，以列行存储 Public tempData(,) As Int32 临时数据 Public PgData(,) As Double PgFile.txt 数据，以列行存储 Public suitableData(,) As Double LandSuitable.txt 数据，以列行存储 Public noDataValue As Int32 无值数据 Public xCor() As Int32 ‘变化元胞的列坐标 Public yCor() As Int32 ‘变化元胞的行坐标 Public upData() As Int32 ‘变化元胞的数据值 Public chgNumber As Int32 变化的点 Public rdm As Random 产生随机数类 Public realUrbanNumber As Int32 实际城市数目 Public simUrbanNumber As Int32 模拟城市数目 　五、基于逻辑回归的CA模型－核心伪代码 　 for 每一行 for 每一列 if data(列,行)=NoData or data(列,行)=1 or data(列,行)=2 then tempdata(列,行)= data(列,行) else 计算领域影响con 计算随机因子影响rdmdata 读取土地适宜性因子suitabledata 读取PgFile.txt中的开发概率Pg 计算总开发概率P=con*rdmdata*suitabledata*Pg if PPthreshold then tempdata(列,行)=1 else tempdata(列,行)=data(列,行) end if 　end for end for 　五、基于逻辑回归的CA模型－核心代码 　 Public Sub stickOne() Ks += 1 ‘Ks为迭代次数 Dim i As Int32, j As Int32 For j = 0 To rows - 1 For i = 0 To cols - 1 　　　　‘如果该元胞值处于无数据状态或者已经是城市或者是水体，则值不变　　　　　　　 If data(i, j) = noDataValue Or data(i, j) = 1 Or data(i, j) = 2 Then tempData(i, j) = data(i, j) 　　　　‘否则，计算该元胞城市开发概率 Else ‘------------第一步，计算领域影响------------------ Dim con As Double = 0 Dim tempI As Int32, tempJ As Int32 tempI = i - 1 If tempI = 0 Then tempJ = j - 1 If tempJ = 0 Then