(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111641581.3 (22)申请日 2021.12.2 9 (71)申请人 中科星通 （廊坊） 信息技 术有限公司 地址 065001 河北省廊坊市经济与技 术开 发区云鹏 道20号 (72)发明人 臧文乾　刘其悦　陶炳成　黄祥志　 赵亚萌　王更科　盛磊　张丽丽　 (74)专利代理 机构 北京华际知识产权代理有限 公司 11676 专利代理师 王超 (51)Int.Cl. G06F 30/20(2020.01) G06Q 10/06(2012.01) (54)发明名称 遥感数字试验场选 址与评价方法 (57)摘要 本发明公开了遥感数字试验场选址与评价 方法， 获取备选区域内的各个选址场地相应的中 低分辨率卫星遥感影像地表反射率数据， 在中低 分辨率条件下计算备选区域内各个选址场地的 空间均一性指标和时间稳定性指标， 初步剔除不 符合条件的选址场地， 缩小数字试验场的备选区 域的范围。 本发 明为了更好的进行卫星影像间的 交叉验证， 在考虑场地均匀性、 稳定性和地面特 征的基础上， 在全国范围内选择多个具有大面积 均匀目标作为试验场， 以此为空间范围， 作为选 择交叉验证影像的依据； 在需要进行场地提取的 目标区域， 在经纬网格范围内计算像元尺度上的 空间一致性与时间稳定性， 通过大范围的计算与 对比， 遴选得到最优的均匀数字试验场。 权利要求书3页 说明书11页 附图4页 CN 114462199 A 2022.05.10 CN 114462199 A 1.遥感数字试验场选 址与评价方法， 其特 征在于， 所述方法包括以下步骤： S1、 确定选址场地的地表类型， 并初步确定需要进行选址评价的备选区域内的各个选 址场地， 并分别对各个选 址场地进行编号； S2、 获取备选区域内的各个选址场地相应的中低分辨率卫星遥感影像地表反射率数 据， 在中低分辨率条件下计算备选区域内各个选址场地的空间均一性指标和时间稳定性指 标， 初步剔除不符合条件的选 址场地， 缩小数字试验场的备选区域的范围； S3、 利用中低分辨率卫星遥感影像对备选区域的范围进行初步筛选后， 使用缩小后的 备选区域内的各个选址场地对应的高分辨率遥感影像数据， 再次计算备选区域内的各个选 址场地的空间均一性指标与时间稳定性指标， 进一步对备选区域内的各个选址场地进 行筛 选， 得到数字试验场的候选范围； S4、 综合比较数字试验场的候选范围内的各目标场地的指标情况,筛选出数字试验场 的选取结果。 2.根据权利要求1所述的遥感数字试验场选址与评价方法， 其特征在于： 数字试验场的 候选范围内的目标场地的获取 方法包括以下步骤： S‑A4.1、 获取所需数字试验场的最小面积及相应的形状， 记为第一窗口， 获取数字试验 场的候选范围； A4.2、 分别在数字试验场的候选范围内各个选址场地中变换第一窗口的位置， 将第一 窗口的位置对应的场地记为目标场地， 不同的目标场地在数字试验场的候选范围中对应选 址场地中对应第一窗口 的位置不同， 将目标场地用{x1， y1， z1}表示， 当数字试验场的候选范围内选址场 地对应的面积小于第 一窗口对应的面积时， 则直接 将该选址场地从数字试验场的候选范围内剔除， 所述x1表示目标场地对应的选 址场地编号， 所述y1表示目标场地中心点 位置对应的经度， 所述z1表示目标场地中心点 位置对应的纬度； 数字试验场的候选范围内的各目标场地的指标情况指的是数字试验场的候选范围内 的各目标场地分别对应的空间均一 性指标及时间稳定性指标。 3.根据权利要求1所述的遥感数字试验场选址与评价方法， 其特征在于： 所述S1中初步 确定需要 进行选址评价的备选区域内的选 址场地的方法包括以下步骤： S1.1、 利用地 面实测照片， 分别得到各个选 址场地的地表类型； S1.2、 对比数据库， 获取 所需数字试验场的地表类型； S1.3、 将S1.1中获取的各个选址场地的地表类型分别与S1.2中所需的数字试验场的地 表类型进行比较， 将满足S1.2中所需的数字试验场的地表类型的选址场地记为需要进行选址评价的备 选区域内的选 址场地。 4.根据权利要求1所述的遥感数字试验场选址与评价方法， 其特征在于： 所述S2中初步 剔除不符合条件的选 址场地， 缩小数字试验场的备选区域的范围的方法包括以下步骤： S2.1、 在中低分辨 率条件下计算备选区域内选 址场地的空间均一 性指标； S2.2、 在中低分辨 率条件下计算备选区域内选 址场地的时间稳定性指标； S2.3、 对比数据库中， 获取所需数字试验场的空间均一性指标阈值区间B1及时间稳定权　利　要　求　书 1/3 页 2 CN 114462199 A 2性阈值区间B2； S2.4、 判断备选区域内每个选址场地对应的空间均 一性指标是否均在B1内且对应的时 间稳定性指标 是否均在 B2内， 当备选区域内选址场地对应的空间均一性指标存在不在B1内的情况或者对应的时间 稳定性指标存在不在B2内的情况， 则将该选址场地 从备选区域内剔除， 反之， 则不对该选址 场地进行处 理； S2.5、 得到缩小后的数字试验场的备选区域的范围。 5.根据权利要求4所述的遥感数字试验场选址与评价方法， 其特征在于： 所述S2.1中在 中低分辨 率条件下计算备选区域内选 址场地的空间均一 性指标的方法包括以下步骤： S2.1.1、 获取当前时间对应的备区域内各个选 址场地对应的遥感影 像； S2.1.2、 选取一个选址场地， 对该选址场地对应的遥感影像中的每个像素进行循环处 理， 以每个像素为窗口中心， 根据第二窗口大小选择窗口内的所有像元， 所述第二窗口为长 为n1、 宽为n2的矩形， 则第二窗口内的像元个数为 n1*n2， 记为 n； S2.1.3、 计算该窗口内所有像元相同波段像素值的均值， 每个像元对应n3个波段的像 素值， 同一像元不同波段对应的像素值不同， 将第k个像元 的第i个波段对应的像素值记为 DNik， 将该窗口内所有像元第i个波段像素值的均值记为 μi， 所述 S2.1.4、 计算该窗口内所有像元相同波段像素值的标准差， 将该 窗口内所有像元第i个 波段像素值的标准差σ i， S2.1.5、 根据空间均一性计算公式 计算该选址 场地对应的遥感影像中该像素不同波段分别对应的空间均一值， 得到该选址场地对应的遥 感影像中该像素对应的空间均一性数 组， 每个空间均一性数 组包含该元素在n3个波段中分 别对应的空间均一 值； S2.1.6、 循环遍历整幅遥感影像， 最终生成该选址场地空间均 一性分布图， 即该选址场 地的空间均一性指标， 所述选址场地空间均一性分布图中各像素对应一个空间均一性数 组。 6.根据权利要求5所述的遥感数字试验场选址与评价方法， 其特征在于： 所述S2.2中在 中低分辨 率条件下计算备选区域内选 址场地的时间稳定性指标的方法包括以下步骤： S2.2.1、 按照时间先后的顺序遍历备区域内每 个选址场地对应所有遥感影 像； S2.2.2、 选取一个选址场地， 分别以该选址场地的每幅遥感影像对应的时间为 中心， 选 取第一单位时间跨度内相应选址场地的所有遥感影像， 将第一单位时间跨度内所有遥感影 像的个数记为m， 对应的时间为中心的遥感影像， 其对应的时间与相应的第一单位时间跨度 的中心对应的时间相同； S2.2.3、 计算第一单位时间跨度内第j幅遥感影 像中像元 数值的均值ρj； S2.2.4、 计算第一单位时间跨度内所有遥感影 像中像元 数值的均值 所述像元 数值指像元对应位置的图像反射 率数值； S2.2.5、 计算该第一单位时间跨度对应的遥感影 像的时间稳定值 δ，权　利　要　求　书 2/3 页 3 CN 114462199 A 3