一种土地利用空间布局智能优化的种群初始化方法.pdf

《一种土地利用空间布局智能优化的种群初始化方法.pdf》由会员分享，可在线阅读，更多相关《一种土地利用空间布局智能优化的种群初始化方法.pdf（12页珍藏版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 103577897 A (43)申请公布日 2014.02.12 CN 103577897 A (21)申请号 201310565628.1 (22)申请日 2013.11.14 G06Q 10/04(2012.01) (71)申请人 武汉大学 地址 430072 湖北省武汉市武昌区珞珈山武 汉大学 (72)发明人 赵翔 刘耀林 刘艳芳 刘殿锋 何建华 焦利民 (74)专利代理机构 武汉科皓知识产权代理事务 所 ( 特殊普通合伙 ) 42222 代理人 严彦 (54) 发明名称 一种土地利用空间布局智能优化的种群初始 化方法 (57) 摘要 一种土地利用空间布局智能优。

2、化的种群初始 化方法， 包括将区域规划年的土地利用数量结构 与区域土地利用现状的数量结构做差值， 得到区 域内各类用地的规划方案与现状面积的像元个数 差值， 得到需要调入像元的地类的列表 A 和需要 调出像元的地类的列表 D ； 在个体上生成一个随 机的位置i得到邻域P ； 统计同时在邻域P和列表 A 中同时出现的多种地类得到集合 Y， 获得集合 Y 中像元数量最多的地类为L ； 对邻域P的所有像元 执行地类变更操作 ； 直到列表 A 和 D 均为空， 则初 始化过程结束并输出初始个体。本发明提高了初 始个体的生成质量和效率， 支持土地利用规划的 土地资源空间布局优化方案高效执行， 进而为实 。

3、现土地资源的合理、 高效、 可持续利用提供重要保 障。 (51)Int.Cl. 权利要求书 1 页 说明书 7 页 附图 3 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书1页 说明书7页 附图3页 (10)申请公布号 CN 103577897 A CN 103577897 A 1/1 页 2 1. 一种土地利用空间布局智能优化的种群初始化方法， 其特征在于 ： 生成各初始个体 分别包括以下步骤， 步骤 1， 读取规划区域的土地利用现状图， 将对应的规划区域现状的土地利用类型信息 作为个体的初始基因值 ； 同时读取规划区域的土地利用适宜性评价结果图， 作为个体的基 。

4、因属性值 ； 步骤 2， 将区域规划年的土地利用数量结构与区域土地利用现状的数量结构做差值， 得 到区域内各类用地的规划方案与现状面积的像元个数差值， 得到需要调入像元的地类的列 表 A 和需要调出像元的地类的列表 D ； 步骤 3， 在个体上生成一个随机的位置 i ； 步骤 4， 获得位置 i 的邻域像元集， 并标记为邻域 P ； 步骤5， 若邻域P内只有一种地类的像元， 则回到步骤3重新生成一个随机的位置， 否则 进入步骤 6 ； 步骤 6， 统计同时在邻域 P 和列表 A 中同时出现的地类， 并标记所得地类集为集合 Y ； 步骤 7， 统计集合 Y 中各地类在邻域 P 中的像元个数， 获。

5、得集合 Y 中像元数量最多的地 类并记为 L ； 步骤 8， 对邻域 P 的所有像元执行地类变更操作 ； 步骤9， 若列表A和D均为空， 则初始化过程结束并输出初始个体， 否则返回步骤3取新 的随机的位置 i 进行邻域像元地类变更， 直至列表 A 和 D 均为空。 2. 根据权利要求 1 所述土地利用空间布局智能优化的种群初始化方法， 其特征在于 ： 步骤 8 包括以下子步骤， 步骤 8.1， 设邻域 P 的像元个数为 R， 邻域 P 的像元依次标记为 0、 1、R-1， 令当前像 元 x=0 ； 步骤 8.2， 判断当前像元的地类 Px 是否属于集合 D， 若当前像元的地类 Px 不属于集合。

6、 D， 则执行步骤 8.6 ； 若当前像元的地类 Px 属于集合 D， 即 Px D， 则进入步骤 8.3 ； 步骤 8.3， 获取当前像元对 Y 中所有地类的适宜性分值 S ； 步骤 8.4， 判断是否当前像元对 Y 中所有地类的适宜性分值均为 0， 若当前像元对 Y 中所有地类的适宜性分值均为 0， 则跳转到步骤 8.6 ； 若当前像元对 Y 中的某几个地类的适宜性分值不为 0， 则进入步骤 8.5 ； 步骤 8.5， 列表 D 中 Px 对应的像元差值 DPx 加 1， 若 DPx=0， 则把 Px 从列表 D 中移除， 然 后令 Px=L， 列表 A 中 L 对应的像元差值减去 1， 。

7、若 AL=0， 则把地类 L 从列表 A 中移除 ； 步骤 8.6， 判断是否列表 A 和列表 D 为空， 若列表 A 为空， 列表 D 为空， 则退出步骤 8 并 直接执行步骤 9， 否则进入步骤 8.7 ； 步骤 8.7， 令 x=x+1， 若 x=R 则终止循环， 否则返回步骤 8.2， 对新的当前像元 x 开始处 理。 权 利 要 求 书 CN 103577897 A 2 1/7 页 3 一种土地利用空间布局智能优化的种群初始化方法 技术领域 0001 本发明属于土地利用空间布局优化配置技术领域， 特别是涉及一种土地利用空间 布局智能优化的种群初始化方法。 背景技术 0002 土地利用。

8、的空间布局优化是土地利用规划工作的重要组成部分， 其核心内容是根 据区域土地的自然 （位置、 形状、 土壤条件等） 等属性特征， 在区域土地利用数量的约束下， 科学合理的配置个类型用地在区域空间上的布局。 进而提高区域土地利用整体的适宜性和 空间紧凑度， 使得土地利用的空间布局更有利于生态环境保护和生产经营管理。 因此， 从本 质上看， 土地利用的空间布局优化问题是一类典型的带约束条件的组合优化问题。 0003 土地利用空间布局优化问题的求解离不开多目标优化技术和空间信息技术的支 持。 因此， 由于技术条件的限制， 早期的土地利用空间布局优化主要依靠规划专家的认为主 观操作完成， 存在工作量大。

9、、 精度较低、 主观性过强等缺点。 因此， 随着地理信息技术和多目 标优化计算的快速发展， 一些新方法和新技术被不断应用到土地利用的空间布局优化问题 求解中。其中， 基于生物群智能的多目标优化模型被认为是一类求解多目标优化问题的有 效工具。因此， 在最近十年内， 国内外有关专家开始将生物群智能算法和空间信息技术相 结合， 提出了面向土地利用空间布局优化配置的生物群智能优化技术。例如， 基于遗传算 法的智能优化技术有关文献 :1 王世忠 . 基于 LUPO 模型的德清县土地利用空间布局优化 J. 农业工程学报， 2012， 28(20):230-238.、 基于蚁群优化算法有关文献 :2 高小永。

10、 . 基 于多目标蚁群算法的土地利用优化配置 D. 武汉 ; 武汉大学， 2010、 基于微粒群算法的智 能优化技术有关文献 :3LIU Y L， LIU D F， LIU Y F， et al.Rural land use spatial allocation in the semiarid loess hilly area in China:Using a Particle Swarm Optimization model equipped with multi-objective optimization techniquesJ. Sci China-Earth Sci， 2012， 5。

11、5(7):1166-1177. 的智能优化技术等。 0004 对于土地利用空间布局优化问题， 其领域知识主要包括土地配置中的适宜性知识 和区域土地利用现状知识。其中， 前者主要是指， 区域土地利用的空间配置过程中， 不能将 某类不适宜的用途分配给某个地块 ； 后者是指， 区域土地利用的空间配置需要考虑地块当 前的用途。 此外， 区域土地利用空间布局方案中， 各类用地的数量必须符合规划期规定的数 量。 0005 土地利用的空间布局优化问题从本质上看是一类带约束条件的多目标组合优化 问题。按照多目标优化问题的一般定义， 土地利用的空间布局优化问题的数学描述如下 : 0006 0007 0008 说。

12、 明 书 CN 103577897 A 3 2/7 页 4 0009 0010 式中， S 表示土地利用的最大适宜度优化目标， 即土地利用空间布局应当本着因地 制宜的原则， 尽可能提高区域内土地利用的适宜性， 做到 “地尽其用” 。S 值越大， 则表明优 化方案对应的土地利用空间布局的土地利用适宜度越优， 土地得到了最有效利用。Comp 为 最大紧凑度优化目标， 即土地利用的空间布局应当紧凑， 以便于生产经营管理。 Comp的取值 越小， 则表明空间布局方案中各斑块的形状越优。 0011 式 (1) 中， N 表示规划区内栅格单元的总数， Suiti为第 i 个栅格单元对其当前所 配置的地类的。

13、适宜性分值， 其取值范围为0-100。 土地对某类用途的适宜性通过土地适宜 性评价工作得到， 其适宜性分值越高， 则表明该土地对该类用地越适宜 . 例如， 某块土地对 耕地的适宜性为 50， 对林地的适宜性为 80， 则表明该土地更适宜于作为园地使用 . 0012 式 (2) 中， LSIj为第 j 个斑块的形状指数， 式 (3) 中， P 表示斑块的周长， A 表示斑 块的面积， LSI为相应斑块的形状指数。 其中， 斑块是指土地利用空间布局方案中， 空间上相 邻、 且用地类型相同的栅格像元构成了斑块。斑块的大小、 形状对区域生态系统功能保护、 物质和能量流动产生重要影响。例如， 对于耕地、。

14、 园地、 林地等农业或生态用地的安排应当 尽可能形成集中连片的大斑块， 以方便管理、 节约成本， 同时也有利于农用地保护。 0013 式 (4) 中， Al为区域内规划年土地资源数量结构中规定的第 l 类用地的面积， N 表 示规划区内栅格单元的个数， at为第 t 个栅格单元所代表地块的真实面积。xtl为一个决策 变量， 当栅格单元的地类为 l 时取值为 1， 否则取值为 0。式 (4) 的内涵是土地利用空间布 局优化方案对应的土地利用数量结构必须等于区域规划设定的土地利用数量结构。 0014 空间布局优化问题同时也是一类典型的空间优化问题。在使用遗传算法、 微粒群 算法等生物群智能算法进行。

15、求解时， 必须对优化问题的解进行编码， 将优化问题的解映射 为生物个体 （染色体、 微粒、 抗体等） 。对于土地利用空间布局优化问题， 其编码原理见附图 1。土地利用类型采取整数编码的方式， 例如图 1 中地类名称采用地类编码对应表示， 耕地 为 1， 园地为 2， 林地为 3， 草地为 4， 城镇为 5， 居民点为 6。空间上的二维网格划分及其相应 的地类编码构成了一个土地利用空间布局方案 （即土地资源空间布局方案） ， 其中的编码是 单元格的某一个状态， 有可能是土地利用现状， 有可能是随机确定的， 也有可能是根据适宜 性分值的大小确定的。附图 1 中， Col 表示像元的列号， Row 。

16、表示行号， 取值与行数 M、 列数 N 有关， A1， A2 等表示该地块对某类用地 1、 2的适宜性分值。在常用的生物群智能算法 中， 遗传算法的染色体， 人工免疫系统模型的抗体， 微粒群算法中微粒和蚁群优化算法中的 蚂蚁， 在本发明中统称为个体。 染色体种群、 抗体种群、 微粒种群和蚂蚁种群， 在本发明中统 称为种群。 编码的主要思路是， 将规划区采用一定大小的栅格进行存储， 一个栅格像元表示 具体的某个地块， 映射为个体中的基因或维度。 通常为了便于计算， 各基因或维度上同时存 储了栅格像元的位置、 土地适宜性信息等。 栅格像元构成的二维空间则对应一个生物个体。 0015 相对于传统的多。

17、目标组合优化问题， 土地利用的空间布局优化配置问题的求解难 度更大， 主要表现为决策变量个数众多。对于一个中等规模的县， 通常决策变量 （即栅格像 元） 的个数通常数以万计。遗传算法、 蚁群优化算法、 微粒群算法和人工免疫系统等生物群 智能优化算法虽然较好的在求解多目标优化问题领域具有传统数学规划方法无可比拟的 优势， 然而， 且劣势也十分明显。即， 需要通过算法的反复迭代来完成 Pareto 最优解的搜 索。例如一种基于多目标人工免疫系统的土地资源空间布局优化方法， 实现过程为 ： 说 明 书 CN 103577897 A 4 3/7 页 5 0016 步骤一、 将规划区土地适宜性评价结果图。

18、采用栅格文件进行存储 ； 0017 步骤二、 编码， 将一种土地资源空间布局方案通过二维整数编码， 映射为人工免疫 系统的抗体 ； 包括用规划区内的一个栅格单元代表实际的一个地块， 对应人工抗体中的一 个基因位， 每个基因位包含了对应地块的空间位置行列号和地类的特征信息， 所述特征信 息包括对所有地类的适宜性分值和当前地块配置的用地类型信息， 适宜性分值根据土地适 宜性评价结果图的栅格文件取得 ； 0018 步骤三、 目标函数设置， 包含两个优化目标， 分别为最大适宜度 S 和最大紧凑度 Comp ； 0019 步骤四、 约束条件设置， 包括设定土地资源空间布局方案对应的土地资源数量结 构必须。

19、等于区域规划设定的土地资源数量结构 ； 0020 步骤五、 抗体种群初始化， 设染色体种群规模为 AN， 结合规划期土地资源数量结 构， 采取随机的方式产生 AN 个初始抗体 ； 0021 步骤六、 抗体目标值向量计算， 包括采用式 （1） 和式 （2） 计算各抗体的目标值向 量， 并进行 Pareto 排序， 获得 Pareto 优化解集， 作为初始的记忆抗体种群 ； 0022 步骤七、 克隆， 包括将当前的记忆抗体根据预设的克隆系数 C 复制 C 份， 形成新的 抗体种群 ； 0023 步骤八、 变异， 包括遍历新的抗体种群， 对每个抗体实施变异操作 ； 0024 步骤九、 抗体目标值向量。

20、计算， 采用式 （1） 和式 （2） 计算各抗体的目标值向量。 0025 步骤十、 种群更新， 将抗体新种群和原有的记忆抗体种群进行合并， 并对合并后的 种群重新进行 Pareto 排序， 获得新的 Pareto 解集， 得到新的记忆抗体种群 ； 0026 步骤十一、 终止条件判断， 若算当前迭代次数达到预设的最大迭代次数 G， 则终止 迭代， 并解码得到 Pareto 优化方案。 0027 随着决策变量的增大， 算法寻找到 Pareto 最优解所需要的迭代次数会急剧增加， 这一过程也被称之为 “维度灾难” 。 采用随机生成初始解的方式对于土地利用空间布局优化 问题并不一定合适， 具体表现为 。

21、：（1） 采用随机搜索方式产生的初始解的质量通常较差， 经 过多次迭代后得到的解集质量甚至还不如区域土地利用空间布局的现状 ；（2） 采用随机搜 索方式产生的初始解容易违反土地利用的适宜性约束条件， 例如将某块土地的用途分配为 其不适合的用途 ： 将河流水面配置为林地等。 因此， 必须结合土地利用空间布局优化问题的 基本特点， 将土地利用空间布局的相关领域知识融入土地利用空间布局优化模型的初始解 生成过程， 以提高算法性能和效率。 发明内容 0028 针对已有智能优化模型中种群初始化算子所存在的问题， 本发明提出了一种基于 领域知识指导的种群初始化方法， 以提高基于生物群智能的土地利用空间布局。

22、优化模型产 生的初始个体的质量和效率。 0029 为实现上述目的， 本发明提供一种土地利用空间布局智能优化的种群初始化方 法 , 生成各初始个体分别包括以下步骤 ： 0030 步骤 1， 读取规划区域的土地利用现状图， 将对应的规划区域现状的土地利用类型 信息作为个体的初始基因值 ； 同时读取规划区域的土地利用适宜性评价结果图， 作为个体 说 明 书 CN 103577897 A 5 4/7 页 6 的基因属性值 ； 0031 步骤 2， 将区域规划年的土地利用数量结构与区域土地利用现状的数量结构做差 值， 得到区域内各类用地的规划方案与现状面积的像元个数差值， 得到需要调入像元的地 类的列表。

23、 A 和需要调出像元的地类的列表 D ； 0032 步骤 3， 在个体上生成一个随机的位置 i ； 0033 步骤 4， 获得位置 i 的邻域像元集， 并标记为邻域 P ； 0034 步骤5， 若邻域P内只有一种地类的像元， 则回到步骤3重新生成一个随机的位置， 否则进入步骤 6 ； 0035 步骤 6， 统计同时在邻域 P 和列表 A 中同时出现的地类， 并标记所得地类集为集合 Y ； 0036 步骤 7， 统计集合 Y 中各地类在邻域 P 中的像元个数， 获得集合 Y 中像元数量最多 的地类并记为 L ； 0037 步骤 8， 对邻域 P 的所有像元执行地类变更操作 ； 0038 步骤 9。

24、， 若列表 A 和 D 均为空， 则初始化过程结束并输出初始个体， 否则返回步骤 3 取新的随机的位置 i 进行邻域像元地类变更， 直至列表 A 和 D 均为空。 0039 而且， 步骤 8 包括以下子步骤， 0040 步骤 8.1， 设邻域 P 的像元个数为 R， 邻域 P 的像元依次标记为 0、 1、R-1， 令当 前像元 x=0 ； 0041 步骤 8.2， 判断当前像元的地类 Px 是否属于集合 D， 0042 若当前像元的地类 Px 不属于集合 D， 则执行步骤 8.6 ； 0043 若当前像元的地类 Px 属于集合 D， 则进入步骤 8.3 ； 0044 步骤 8.3， 获取当前像。

25、元对 Y 中所有地类的适宜性分值 S ； 0045 步骤 8.4， 判断是否当前像元对 Y 中所有地类的适宜性分值均为 0， 0046 若当前像元对 Y 中所有地类的适宜性分值均为 0， 则跳转到步骤 8.6 ； 0047 若当前像元对 Y 中的某几个地类的适宜性分值不为 0， 则进入步骤 8.5 ； 0048 步骤 8.5， 列表 D 中 Px 对应的像元差值 DPx 加 1， 若 DPx=0， 则把 Px 从列表 D 中移 除， 然后令 Px=L， 列表 A 中 L 对应的像元差值减去 1， 若 AL=0， 则把地类 L 从列表 A 中移除 ； 0049 步骤 8.6， 判断是否列表 A 。

26、和列表 D 为空， 若列表 A 为空， 列表 D 为空， 则退出步骤 8 并直接执行步骤 9， 否则进入步骤 8.7 ； 0050 步骤8.7， 令x=x+1， 若x=R则终止循环， 否则返回步骤8.2， 对新的当前像元x开始 处理。 0051 相对于目前在实际工程中广泛应用的各种土地利用空间布局优化生物群智能优 化过程的种群初始化方法， 本发明主要的效果有 ：（1） 根据土地利用空间布局优化问题的 特点， 在种群初始化过程中引入领域知识的指导， 避免了违反约束条件解的产生， 有助于简 化后续的变异和交叉算子 ；（2） 充分利用区域土地利用空间布局现状信息， 同时融入部分 随机因素， 提高初始。

27、个体的质量， 减少算法收敛所需要的迭代次数， 极大缩短流程执行时 间， 节约了系统资源。本发明的技术方案支持土地利用规划的土地资源空间布局优化方案 高效执行， 进而为实现土地资源的合理、 高效、 可持续利用提供重要保障。 说 明 书 CN 103577897 A 6 5/7 页 7 附图说明 0052 图 1 为本发明实施例的面向土地利用空间布局优化问题求解的生物群智能编码 原理图。 0053 图 2 为本发明实施例的种群初始化原理图。 0054 图 3 为本发明实施例的栅格像元的邻域空间示意图。 0055 图 4 为本发明实施例的种群初始化主要流程图。 0056 图 5 为本发明实施例的种群。

28、初始化方法与传统随机初始化方法的性能对比图。 具体实施方式 0057 本发明提供的技术方案是一种领域知识指导的土地利用空间布局生物群智能优 化的种群初始化的快速、 高效方法。 本发明不针对特定的生物群智能算法， 仅对各类算法的 初始解的构建技术进行改进， 以提高算法初始解集的质量。 该方法的基本思路是 ： 从土地利 用现状空间布局为基础， 通过比较土地利用数量结构优化方案 （规划的用地结构） 和土地利 用现状数量结构， 设计初始化算子， 确保所有初始解对应的数量结构能完全符合预先设定 的土地利用数量结构优化方案。以下结合附图和实施例进行详细说明。 0058 本发明实施例的领域知识指导的土地利用。

29、空间布局生物群智能优化的种群初始 化方法， 基本原理见附图 2， 主要流程见附图 4， 可由本领域技术人员采用计算机软件技术 实现自动运行流程， 生成各初始个体分别具体实现过程如下 ： 0059 （1） 根据编码方案， 读取规划区域的土地利用现状图， 将其对应的规划区域现状的 土地利用类型信息作为个体的初始基因值 ； 同时读取规划区域的土地利用适宜性评价结果 图， 将其作为个体的基因属性值。如附图 1 所示， 个体中的某个基因位 Col， Row 的土地利 用类型信息来自于该区域的土地利用现状信息， 基因属性信息 A1,A2等来自对应土地类 型的适宜性评价结果。 0060 （2） 将区域规划年。

30、的土地利用数量结构与区域土地利用现状的数量结构做差值， 得到区域内各类用地的规划方案与现状面积差值， 并以像元个数进行表达。 例如， 区域城镇 建设用地的现状面积为1000个像元， 优化方案中的城镇建设用地面积为1500个像元， 则其 面积差值为 500 个像元。依次类推， 得到各类用地面积差值的像元表达。其中， 差值为正值 的地类需要调入相应数量像元土地， 差值为负值的地类则需要调出相应数量像元土地。在 此基础上构建需要调入像元的地类列表 A 和需要调出像元的地类列表 D， 分别记录各类用 地需要调入或需要调出的像元个数。如图 2 中， 耕地、 园地、 林地、 草地、 城镇和居民点的现 状面。

31、积分别为 C1、 C2、 C3、 C4、 C5、 C6， 规划面积分别为 P1、 P2、 P3、 P4、 P5、 P6， 求得相应面积 差值为 a1、 a2、 a3、 a4、 a5、 a6。其中 a1、 a3、 a4 为负值， 需要调出像元的地类列表 D 记录了 耕地、 林地、 草地分别需要调出的数目 a1、 a3、 a4 ； 其中 a2、 a5、 a6 为正值， 需要调入像元的 地类列表 A 记录了园地、 城镇和居民点分别需要调入的数目 a2、 a5、 a6。以列表 A 和 D 为数 量结构约束进行土地利用空间的布局调整， 可得到知识指导的初始化算子。 0061 （3） 在个体上生成一个随机。

32、的位置 i， 其中 i 的取值范围为 1 到 MN。即取规划区 域土地利用现状信息上的任一像元。 0062 （4） 获得随机的位置 i 的邻域像元集， 并标记为邻域 P。其中， 位置 i 的邻域一般 可按照栅格像元的8邻域空间的定义获得。 栅格像元的邻域空间见附图3， 随机的位置i所 说 明 书 CN 103577897 A 7 6/7 页 8 在当前像元周围有 8 个领域像元。 0063 （5） 若邻域 P 内只有一种地类的像元， 则回到步骤 （3） 重新生成一个随机位置。否 则进入 （6） 。 0064 （6） 此时邻域P内地类的种数大于1， 则统计同时在P和列表A中同时出现的地类， 并标。

33、记所得地类集为集合 Y。 0065 （7） 统计集合 Y 中各地类在邻域 P 中的像元个数， 并按照像元个数对 Y 中的地类进 行降序排列， 并获得 Y 中像元数量最多的地类， 记为 L。 0066 （8） 对邻域 P 的所有像元执行地类变更操作。 0067 （9） 若列表 A 和 D 均为空， 则初始化过程结束并输出初始解， 即初始个体。否则返 回步骤 （3） 取新的随机的位置 i 进行邻域像元地类变更， 直至列表 A 和 D 均为空即可完成 初始抗体的生成。 0068 执行上述初始化过程可以得到一个初始个体， 根据土地利用空间布局智能优化的 具体初始化种群需要的个体数目 n， 可以 n 次。

34、执行上述初始化过程， 得到 n 个初始个体。具 体实施时， 在执行第 2、 3、 n 个初始个体的相应初始化过程时， 可以省略步骤 1、 2， 直接利 用执行第 1 个初始个体的相应初始化过程时步骤 1、 2 的结果， 以节约系统资源。 0069 上述步骤 （8） 中， 邻域 P 的像元地类变更操作的基本步骤如下 ： 0070 （8.1） 设邻域 P 的像元个数为 R， 邻域 P 的像元依次标记为 0、 1、R-1， 开始依次 遍历 P 中的每一个像元， 令当前像元 x=0。实施例中 R=8。 0071 （8.2） 判断当前像元的地类 Px 是否属于集合 D， 分别处理 ： 0072 若当前像。

35、元的地类 Px 不属于集合 D， 则执行 （8.6） 。 0073 若当前像元的地类 Px 属于集合 D， 即 Px D， 则进入步骤 （8.3） ； 0074 （8.3） 获取当前像元对 Y 中所有地类的适宜性分值 S。步骤 1 读取规划区域的土 地利用适宜性评价结果图作为个体的基因属性值， 因此可知当前像元对 Y 中所有地类的适 宜性分值 S。当前像元对某类用地的 S 值越大则表明其越适合作为该类用地。当前像元对 某类用地的 S=0 时， 则表明地块不能转换为该类用地。 0075 （8.4） 判断是否当前像元对 Y 中所有地类的适宜性分值均为 0， 0076 若当前像元对 Y 中所有地类的。

36、适宜性分值均为 0， 则跳转到步骤 （8.6） 。 0077 若当前像元对 Y 中的某几个地类的适宜性分值不为 0， 则进入步骤 （8.5） 。 0078 （8.5） 列表 D 中 Px 对应的像元差值 DPx 加 1( 即 DPx=DPx+1)， 若 DPx=0， 则把 Px 从 D 中移除， 然后取步骤 7 所得 Y 在邻域 P 中所占像元数最多的地类 L 作为当前像元的地类， 即 Px=L。列表 A 中 L 对应的像元差值减去 1( 即 AL=AL-1)， 若减去后 AL=0， 则把地类 L 从 A 中移除 ； 0079 （8.6） 判断是否 A 列表和 D 列表为空， 若列表 A 为空。

37、， 列表 D 为空， 则终止循环。根 据本发明设计， 列表 A 和列表 D 将同时达到空的状态。当两者为空时， 退出步骤 （8） ， 并直接 执行步骤 （9） ， 即初始化过程结束并输出初始解。否则进入 （8.7） 。 0080 （8.7） x=x+1， 若 x=R 则说明对邻域 P 的像元 0、 1、 R-1 均已处理完成， 终止循环， 否则返回步骤 （8.2） ， 对新的当前像元X开始处理。 选取某地土地利用空间布局优化问题作 为实施案例， 需要对其耕地、 园地、 林地、 草地、 城镇建设用地和农村居民点用地的空间布局 进行优化。 分别采用本发明设计的领域知识指导的种群初始化算子和传统的随。

38、机初始化算 说 明 书 CN 103577897 A 8 7/7 页 9 子进行种群初始化操作， 得到初始种群的Pareto前端见附图5。 由图可知， 本发明提出的种 群初始化算子所产生的初始解集的质量较随机初始化方法所产生的初始解集有了较大提 高， 紧凑度值和适宜度值都更高。 0081 以上仅是用以说明本发明的具体实施案例而已， 并非用以限定本发明的可实施范 围。熟悉本领域的技术人员在不违背本发明所指示的精神与原理下所完成的一切等效变 形、 替换或修饰， 仍包含在本发明权利要求所限定的范围内。 说 明 书 CN 103577897 A 9 1/3 页 10 图 1 图 2 图 3 说 明 书 附 图 CN 103577897 A 10 2/3 页 11 图 4 说 明 书 附 图 CN 103577897 A 11 3/3 页 12 图 5 说 明 书 附 图 CN 103577897 A 12 。