(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202110693013.1 (22)申请日 2021.06.22 (65)同一申请的已公布的文献号 申请公布号 CN 113469512 A (43)申请公布日 2021.10.01 (73)专利权人 深圳市城市规划设计 研究院有限 公司 地址 518028 广东省深圳市福田区振兴 路3 号建艺大厦10楼 (72)发明人 单樑　李启军　岳隽　唐文魁　 陈君丽　袁富林　 (74)专利代理 机构 广州嘉权专利商标事务所有 限公司 4 4205 专利代理师 黄广龙(51)Int.Cl. G06Q 10/06(2012.01) G06Q 50/26(2012.01) G06F 30/27(2020.01) G06T 17/05(2011.01) (56)对比文件 CN 110334864 A,2019.10.15 CN 110852532 A,2020.02.28 CN 106294289 A,2017.01.04 审查员 耿玲 (54)发明名称 自然保护区划定方法、 装置、 电子设备及存 储介质 (57)摘要 本发明公开了一种自然保护区划定方法、 装 置、 电子设备及存储介质， 涉及生态规划技术领 域， 其中自然保护 区划定方法包括： 获取待分析 土地数据； 从待分析土地数据中提取出土地物种 分布数据， 从待分析土地数据中提取出威胁源数 据； 利用预设生境质量模型、 土地物种分布数据 和威胁源数据计算出生境质量评价数据； 利用预 设智能体演化模型和生境质量评价数据计算出 自然保护区数据； 根据自然保护区数据和生境质 量评价数据生成自然保护区。 上述自然保护区划 定方法， 能够生态适宜性评价结果模拟、 预测、 优 化及显示指示物种在空间聚集过程来划分自然 保护区， 弥补常规元胞自动机建模在对复杂地理 时空过程进行模拟与优化时存在的功能不足。 权利要求书3页 说明书10页 附图5页 CN 113469512 B 2022.08.30 CN 113469512 B 1.自然保护区划定方法， 其特 征在于， 包括： 获取待分析土地数据， 所述待分析土地数据包括： 研究区域的土地利用数据和DEM高程 数据， 所述土地利用数据和所述DE M高程数据的数据类型均为 栅格数据； 从所述待分析土地数据中提取出土地物种分布数据， 从所述待分析土地数据中提取出 威胁源数据包括： 根据所述待分析土地数据拟定威胁因子数据； 获取所述威胁因子数据对应的威胁因子类型； 根据所述 威胁因子类型对所述 威胁因子数据进行 赋值处理， 得到所述 威胁源数据； 所述威胁源数据包括： 生态胁迫因子、 生境类型对胁迫因子敏感度、 胁迫因子最大影响 距离， 所述 威胁源数据可以根据拟定的威胁因子类型确定； 利用预设生境质量模型、 所述土地物种分布数据和所述威胁源数据计算出生境质量评 价数据， 所述预设生 境质量模型为Habitat  Quality模型； 利用预设智能体演化模型和所述生境质量评价数据计算出自然保护区数据； 其中， 所 述生境质量评价数据包括： 生境质量得分； 所述自然保护区数据为形成自然保护区的概率， 所述预设智能体演化模型用于计算t 时刻第k个智能体选择一个栅格单元形成自然保护区 的概率， 具体为： 式中， 表示智能体k于时刻t对栅格单元i的选择概率； allowedk表示智能体k下次 游历可以造访的栅格单元集合； α 为信息素启发因子， 反映智能体在运动过程中受信息素浓 度τi(t)的作用； β 为期望信息启发因子， 用以表征智能体在运动过程中受启发式函数ηi(t) 的作用； 所述启发函数 ηi(t)为： 式中， 为栅格单元i的生境质量得分； 为研究区所有栅格单元的生境质量得分 之和； 信息素浓度 τi(t)的计算公式为： 式中， d(x)是保护区单 元x到中心细胞i的距离； 紧凑度指数根据以下比值 函数⑦、⑧定义： Uutility＝weS+wcP     ⑧权　利　要　求　书 1/3 页 2 CN 113469512 B 2式⑦、⑧中， A是保护区的总面积， L是受保护场景保护区的周长； Uutility表示效益目 标函数， S和P分别为生境质量和空间的紧凑性指数值， we和wc分别为对应的权重， 预计最佳 保护模式所述效益目标函数应在平均总生态适宜性和所述紧凑度指数产生 最高值； 根据所述自然保护区数据和所述 生境质量评价数据生成自然保护区。 2.根据权利要求1所述的方法， 其特 征在于， 所述获取待分析土地数据， 包括： 获取土地利用数据和DE M高程数据； 对所述土地利用数据和所述DE M高程数据进行关联 预处理， 得到所述待分析土地数据。 3.根据权利要求2所述的方法， 其特征在于， 所述从所述待分析土地数据中提取出土地 物种分布数据， 包括： 对所述待分析土地数据进行物种分布提取， 得到所述土地物种分布数据。 4.根据权利要求1所述的方法， 其特征在于， 所述利用预设生境质量模型、 所述土地物 种分布数据和所述 威胁源数据计算出生 境质量评价数据， 包括： 将所述土地物种分布数据和所述威胁源数据代入所述预设生境质量模型中， 计算出生 境退化指数和生 境质量评价指数； 根据所述 生境退化指数和所述 生境质量评价指数 得到所述 生境质量评价数据。 5.根据权利要求1所述的方法， 其特征在于， 所述利用预设智能体演化模型和所述生境 质量评价数据进行计算， 得到自然保护区数据， 包括： 获取生态保护目标参数和效用函数； 根据所述预设智能体演化模型、 所述生态保护目标参数和所述效用函数进行区域优化 计算， 得到所述自然保护区数据。 6.根据权利要求1所述的方法， 其特征在于， 所述根据所述自然保护区数据和所述生境 质量评价数据生成自然保护区， 包括： 根据所述自然保护区数据和所述 生境质量评价数据生成自然保护区划定数据； 根据所述自然保护区划定数据加载 得到所述自然保护区。 7.自然保护区划定装置， 其特 征在于， 包括： 获取模块， 用于获取待分析土地数据， 所述待分析土地数据包括： 研究区域的土地利用 数据和DE M高程数据， 所述土地利用数据和所述DE M高程数据的数据类型均为 栅格数据； 提取模块， 用于从所述待分析土地数据中提取出土地物种分布数据， 从所述待分析土 地数据中提取 出威胁源数据包括： 根据所述待分析土地数据拟定威胁因子数据； 获取所述威胁因子数据对应的威胁因子类型； 根据所述 威胁因子类型对所述 威胁因子数据进行 赋值处理， 得到所述 威胁源数据； 所述威胁源数据包括： 生态胁迫因子、 生境类型对胁迫因子敏感度、 胁迫因子最大影响 距离， 所述 威胁源数据可以根据拟定的威胁因子类型确定； 生境质量评价模块， 用于利用预设生境质量模型、 所述土地物种分布数据和所述威胁 源数据计算出生 境质量评价数据， 所述预设生 境质量模型为Habitat  Quality模型； 智能体演化模块， 用于利用预设智能体演化模型和所述生境质量评价数据计算出自然 保护区数据， 其中， 所述生境质量评价数据包括： 生境质量得分； 所述自然保护区数据为形 成自然保护区的概率， 所述预设智能体演化模型用于计算t时刻第k个智能体选择一个栅格权　利　要　求　书 2/3 页 3 CN 113469512 B 3