[发明专利]遥感数字试验场选址与评价方法

权利要求书

1.遥感数字试验场选址与评价方法，其特征在于，所述方法包括以下步骤：

S1、确定选址场地的地表类型，并初步确定需要进行选址评价的备选区域内的各个选址场地，并分别对各个选址场地进行编号；

S2、获取备选区域内的各个选址场地相应的中低分辨率卫星遥感影像地表反射率数据，在中低分辨率条件下计算备选区域内各个选址场地的空间均一性指标和时间稳定性指标，初步剔除不符合条件的选址场地，缩小数字试验场的备选区域的范围；

S3、利用中低分辨率卫星遥感影像对备选区域的范围进行初步筛选后，使用缩小后的备选区域内的各个选址场地对应的高分辨率遥感影像数据，再次计算备选区域内的各个选址场地的空间均一性指标与时间稳定性指标，进一步对备选区域内的各个选址场地进行筛选，得到数字试验场的候选范围；

S4、综合比较数字试验场的候选范围内的各目标场地的指标情况,筛选出数字试验场的选取结果；

所述S2中初步剔除不符合条件的选址场地，缩小数字试验场的备选区域的范围的方法包括以下步骤：

S2.1、在中低分辨率条件下计算备选区域内选址场地的空间均一性指标；

S2.2、在中低分辨率条件下计算备选区域内选址场地的时间稳定性指标；

S2.3、对比数据库中，获取所需数字试验场的空间均一性指标阈值区间B1及时间稳定性阈值区间B2；

S2.4、判断备选区域内每个选址场地对应的空间均一性指标是否均在B1内且对应的时间稳定性指标是否均在B2内，

当备选区域内选址场地对应的空间均一性指标存在不在B1内的情况或者对应的时间稳定性指标存在不在B2内的情况，则将该选址场地从备选区域内剔除，反之，则不对该选址场地进行处理；

S2.5、得到缩小后的数字试验场的备选区域的范围；

所述S2.1中在中低分辨率条件下计算备选区域内选址场地的空间均一性指标的方法包括以下步骤：

S2.1.1、获取当前时间对应的备区域内各个选址场地对应的遥感影像；

S2.1.2、选取一个选址场地，对该选址场地对应的遥感影像中的每个像素进行循环处理，以每个像素为窗口中心，根据第二窗口大小选择窗口内的所有像元，所述第二窗口为长为n1、宽为n2的矩形，则第二窗口内的像元个数为n1*n2，记为n；

S2.1.3、计算该窗口内所有像元相同波段像素值的均值，每个像元对应n3个波段的像素值，同一像元不同波段对应的像素值不同，将第k个像元的第i个波段对应的像素值记为DN_ik，将该窗口内所有像元第i个波段像素值的均值记为μi，所述

S2.1.4、计算该窗口内所有像元相同波段像素值的标准差，将该窗口内所有像元第i个波段像素值的标准差σi，

S2.1.5、根据空间均一性计算公式计算该选址场地对应的遥感影像中该像素不同波段分别对应的空间均一值，得到该选址场地对应的遥感影像中该像素对应的空间均一性数组，每个空间均一性数组包含该像素在n3个波段中分别对应的空间均一值；

S2.1.6、循环遍历整幅遥感影像，最终生成该选址场地空间均一性分布图，即该选址场地的空间均一性指标，所述选址场地空间均一性分布图中各像素对应一个空间均一性数组；

所述S2.2中在中低分辨率条件下计算备选区域内选址场地的时间稳定性指标的方法包括以下步骤：

S2.2.1、按照时间先后的顺序遍历备区域内每个选址场地对应所有遥感影像；

S2.2.2、选取一个选址场地，分别以该选址场地的每幅遥感影像对应的时间为中心，选取第一单位时间跨度内相应选址场地的所有遥感影像，将第一单位时间跨度内所有遥感影像的个数记为m，对应的时间为中心的遥感影像，其对应的时间与相应的第一单位时间跨度的中心对应的时间相同；

S2.2.3、计算第一单位时间跨度内第j幅遥感影像中像元数值的均值ρj；

S2.2.4、计算第一单位时间跨度内所有遥感影像中像元数值的均值

所述像元数值指像元对应位置的图像反射率数值；

S2.2.5、计算该第一单位时间跨度对应的遥感影像的时间稳定值δ，

所述

S2.2.6、循环遍历该选址场地的所有遥感影像分别对应的时间，得到该选址场地对应的时间稳定性数组，即该选址场地对应的时间稳定性指标；

所述S4中筛选出数字试验场的选取结果的方法包括以下步骤：

S4.1、获取目标场地对应的高分辨率遥感影像相应的空间均一性指标与时间稳定性指标；

S4.2、分别获取目标场地对应的空间均一性分布图中每个空间均一性数组中各个空间均一值的平均值，获取目标场地对应的空间均一性分布图中每个空间均一性数组对应平均值的中位数，将第r个目标场地的空间均一性指标对应的中位数记为Qr；

S4.3、获取目标场地对应的时间稳定性数组中各个时间稳定值的平均值，将第r个目标场地的时间稳定性指标对应的平均值记为Hr；

S4.4、对各个目标场地对应的空间均一性指标与时间稳定性指标进行归一化处理，将第r个目标场地对应的归一化处理结果记为Gr，所述Gr＝e1*Qr+e2*Hr，

其中，e1表示第一系数，e2表示第二系数；

S4.5、对比数据库，获取所需数字试验场的归一化处理结果预设值，记为R；

S4.6、计算r为不同值时，第r个目标场地与所需数字试验场的偏差值Pr，

所述其中r为正整数；

S4.7、获取r为不同值时，Pr为最小值时对应的r值，第r个目标场地即为数字试验场的选取结果。