(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210879888.5 (22)申请日 2022.07.25 (71)申请人 广东无限阵列科技有限公司 地址 510000 广东省广州市越秀区德 政北 路401-409号 4-9楼968房 (72)发明人 伍容容　 (74)专利代理 机构 广州三环 专利商标代理有限 公司 44202 专利代理师 吴松滨 (51)Int.Cl. G06Q 10/04(2012.01) G06Q 50/16(2012.01) G06Q 50/26(2012.01) G06V 20/13(2022.01) G06V 20/10(2022.01)G06V 10/77(2022.01) (54)发明名称 基于数字孪生的生产建设活动扰动变化预 测方法及系统 (57)摘要 本发明涉及遥感地理信息技术领域， 尤其涉 及一种基于数字孪生的生产建设活动扰动变化 预测方法及系统， 包括： 从预设的空间数据库中 读取生产建设活动扰动图斑矢量； 获取图斑矢量 在其外包矩形缓冲区内的卫星遥感影像； 根据所 述卫星遥感影像， 建立基于数字孪生技术的图斑 变化预测数字孪生模型； 基于所述图斑变化预测 数字孪生模 型， 预测生产建设活动扰动的持续扰 动概率。 本发明根据数字孪生技术的思路， 建立 图斑变化预测数字孪生模型， 并计算生产建设活 动扰动图斑的持续扰动概率， 基于持续扰动概率 进行警报， 从而提高遥感监管工作中对生产建设 项目/非生产建 设项目判别的准确率。 权利要求书3页 说明书9页 附图4页 CN 115222134 A 2022.10.21 CN 115222134 A 1.一种基于数字 孪生的生产建 设活动扰动变化预测方法， 其特 征在于， 所述方法包括： 从预设的空间数据库中读取生产建 设活动扰动图斑矢量； 获取所述生产建设活动扰动图斑矢量在其外包矩形缓冲区内的卫星遥感影 像； 根据所述 卫星遥感影 像， 建立基于数字 孪生技术的图斑 变化预测数字 孪生模型； 基于所述图斑 变化预测数字 孪生模型， 预测生产建 设活动扰动的持续扰动概 率。 2.根据权利要求1所述的一种基于数字孪生的生产建设活动扰动变化预测方法， 其特 征在于， 所述获取所述生产建设活动扰动图斑矢量在其外包矩形缓冲区内的卫星遥感影像 的步骤包括： 获取当前生产建设活动扰动图斑矢量在其外包矩形缓冲区内的第一卫星遥感影像和 第二卫星遥感影像； 其中， 所述第一卫星遥感影像包括历史最靠近当前时间节点的卫星遥 感影像数据， 所述第二卫星遥感影像包括历史次最靠近当前时间节点的卫星遥感影像数 据； 分别读取所述第 一卫星遥感影像和所述第 二卫星遥感影像中的数据， 生成对应的第 一 卫星遥感影 像数组和第二 卫星遥感影 像数组。 3.根据权利要求2所述的一种基于数字孪生的生产建设活动扰动变化预测方法， 其特 征在于， 所述根据所述卫星遥感影像， 建立基于数字孪生技术的图斑变化预测数字孪生模 型的步骤 包括： 分别计算所述第一卫星遥感影像数组的第一比值建筑指数数组和所述第二卫星遥感 影像数组的第二比值建筑指数 数组； 将所述第一比值建筑指数数组和所述第 二比值建筑指数数组分别进行二值化处理， 得 到对应的第一 二值化数组和第二 二值化数组； 对所述第一二值化数组和所述第 二二值化数组分别进行主成分分析， 得到对应的第 一 主成分值和第二主成分值； 根据所述第一主成分值和所述第二主成分值， 计算得到 趋势因子 。 4.根据权利要求3所述的一种基于数字孪生的生产建设活动扰动变化预测方法， 其特 征在于： 所述第一比值建筑指数 数组的计算公式为： 所述第二比值建筑指数 数组的计算公式为： 式中， I1为第一比值建筑指数数组； I2为第二比值建筑指数数组； B1_n为第一卫星遥感影 像数组中的近红外波段反射率值或DN值； B1_r为第一卫星遥感影像数组中的红波段反射率 值或DN值； B1_g为第一卫星遥感影像数组中的绿波段反射率值或DN值； B2_n为第二卫星遥感 影像数组中的近红外波段反射率值或DN值； B2_r为第二卫星遥感影像数组中的红波段反射 率值或DN值； B2_g为第二卫星遥感影像数 组中的绿波 段反射率值 或DN值； α 为土壤调节因子，权　利　要　求　书 1/3 页 2 CN 115222134 A 2取值范围为[0， 1]； β 为水体调节因子， 取值范围为[0， 1]。 5.根据权利要求3所述的一种基于数字孪生的生产建设活动扰动变化预测方法， 其特 征在于， 所述 二值化处理的计算公式为： 式中， BA1表示第一二值化数组； I1表示第一比值建筑指数数组； BA2表示第二二值化数 组； I2表示第二比值建筑指数 数组。 6.根据权利要求3所述的一种基于数字孪生的生产建设活动扰动变化预测方法， 其特 征在于： 所述第一主成分值的计算公式为： 其中， pc1表示第一主成分值； N1表示第一二值化数组的元素个数； BA1_i表示第一二值化 数组的第i个元 素； 所述第二主成分值的计算公式为： 其中， pc2表示第二主成分值； N2表示第二二值化数组的元素个数； BA2_j表示第二二值化 数组的第j个元 素。 7.根据权利要求3所述的一种基于数字孪生的生产建设活动扰动变化预测方法， 其特 征在于， 所述趋势因子的计算公式为： 式中， δ表示趋势因子； pc1表示第一主成分值； pc2表示第二主成分值。 8.根据权利要求1所述的一种基于数字孪生的生产建设活动扰动变化预测方法， 其特 征在于， 所述持续扰动概 率的计算公式为： 式中， P表示持续扰动概 率； δ表示趋势因子 。 9.根据权利要求1 ‑8任意一项所述的一种基于数字孪生的生产建设活动扰动变化预测 方法， 其特征在于， 所述基于所述图斑变化预测数字孪生模型， 预测生产建设活动扰动的持 续扰动概 率的步骤之后， 所述方法还 包括： 判断所述持续扰动概 率与预设的扰动阈值的关系； 若所述持续扰动概率大于预设的扰动阈值， 则提示该图斑为生产建设项目， 否则提示 该图斑为非生产建 设项目。权　利　要　求　书 2/3 页 3 CN 115222134 A 3