(19)中华 人民共和国 国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202110845572.X (22)申请日 2021.07.26 (65)同一申请的已公布的文献号 申请公布号 CN 113536688 A (43)申请公布日 2021.10.2 2 (73)专利权人 武汉大学 地址 430072 湖北省武汉市武昌区珞珈山 武汉大学 (72)发明人 沙宗尧　方虎盛　邱岱　 (74)专利代理 机构 武汉科皓知识产权代理事务 所(特殊普通 合伙) 42222 代理人 罗飞 (51)Int.Cl. G06F 30/27(2020.01) G06N 20/10(2019.01)(56)对比文件 CN 110458048 A,2019.1 1.15 CN 103413293 A,2013.1 1.27 CN 108345897 A,2018.07.31 WO 201720 5385 A1,2017.1 1.30 审查员 张驰 (54)发明名称 一种基于多尺度邻域分析的景观破碎度模 型的构建方法 (57)摘要 本发明设计了一种基于多尺度邻域分析的 景观破碎Frag模型， 提供了Frag的构建方法及 Frag变化分析方法， 包括： 获取研究区的土地分 类覆盖栅格数据产品； 根据土地覆盖分类， 逐像 素计算给定目标类型标记的同类最邻近距离； 设 计尺度相关的景观破碎指数模型Frag， 并使用邻 域分析的方法计算统计Frag； 基于时间序列的 Frag计算结果， 使用机器 学习方法进行Frag时空 变化与相关联的驱动因子分析和尺度影 响分析。 利用本发明方法得到的景观破碎度时空分布特 征， 可以从尺度和结构上， 精确表示目标类型景 观的破碎分布情况， 更好地反映实际景观布局的 差异和动态， 并可 以进一步根据具体问题， 选择 适合的机器学习方法， 实现景观变化驱动力因子 分析， 具有客观、 灵活、 通用、 简便的特点。 权利要求书2页 说明书10页 附图3页 CN 113536688 B 2022.04.15 CN 113536688 B 1.一种基于多尺度邻域分析的景观破碎度模型的构建方法， 其特 征在于， 包括： S1： 从遥感数据或调查数据中获取研究区土地类型覆盖时空数据， 其中， 获取的研究区 土地类型覆盖时空数据为 栅格格式； S2： 以栅格像素为基本单位， 对研究区土地类型覆盖时空数据中感兴趣的土地覆盖类 型C进行二 值化； S3： 采用邻域分析方法计算所有值为1的像素与同类型的其它像素的最邻近距离， 其 中， 值为1的像素表示覆盖类型为C的像素； S4： 在给定的任意尺度上， 根据最邻近距离， 确定用以表征目标类型C分布距离差异的 第一指标以及用以表征景观分布受尺度和密度差异影响的第二指标， 并基于第一指标和 第 二指标， 构建景观破碎度模型。 2.如权利要求1所述的景观破碎度模型的构建方法， 其特 征在于， 步骤S2包括： 将研究区土地类型覆盖时空数据中覆盖类型为C的像素值置为1， 而其它类型置为0 。 3.如权利要求2所述的景观破碎度模型的构建方法， 其特 征在于， 步骤S3包括： 循环过滤二值化后的像素中像素值为1的像素Pi， 在给定以Pi为中心、 边长为奇数w个像 素的最大预定义正方形窗口W范围内， 寻找其与同类型像素Pj的间距， 像素Pi与同类型像素 Pj距离Disi， j的计算公式为： 式中， Δy表示Pi与像素Pj之间在纵向的间距， Δx表示Pi与像素Pj之间在横向的间距， 最大预定义正方形窗口为邻域窗口； 根据计算得到的像素Pi与同类型像素Pj距离， 得到最邻近距离： NDi＝Min{Disi， j， j＝1， 2， …， m}    (2) 式中， NDi表示最邻近距离， 为Disi， j中的最小值， m为窗口W范围内除Pi外标记为1的像素 总数。 4.如权利要求1所述的景观破碎度模型的构建方法， 其特 征在于， 步骤S4包括： 对某邻域窗口范围的最邻近距离特征， 进行统计分析， 得到邻域窗口内的平均景观破 碎特征， 作为表征类型C在尺度s的景观分布差异的指标， 其中， 平均景观破碎特征的计算 公 式为： AVG_NDs， i＝∑NDs， i/nn    (3) 其中， ∑NDs， i表示对当前以s为尺度、 i为中心点的景观破碎总特征的表示， AVG_NDs， i表 示土地覆盖的类型为C时， 在邻域窗口W内， 同类型分布最短距离的平均值， 即第一指标， 其 中尺度s由窗口W定义的空间大小所表示， NDk表示w内像素值标记为1、 序号为k的像素的最邻近距离， nn为当前s尺度下窗口W内像素值标记为1的像素总数， k为窗 口内循环计算的像序号； 另外， 参照景观破碎总特征， 得到用以表征目标类型C的景观布局， 在尺度和密度差异 上的影响指标A VG_NDs， 即第二指标： AVG_NDs＝∑NDs， i/n    (4) n表示当前尺度s、 邻域窗口W内所有像素的总数目， 所有像素包括 值为1和0的像素； 根据平第一指标和第二指标， 构建景观破碎度模型：权　利　要　求　书 1/2 页 2 CN 113536688 B 2Frag＝AVG_NDs， i‑AVG‑NDs   (5) 其中， Fra g表示景观破碎度模型， 量纲为以像素为基本单位的距离表示。权　利　要　求　书 2/2 页 3 CN 113536688 B 3