[发明专利]一种基于证据权和分形的全月着陆区选址方法

权利要求书

1.一种基于证据权和分形的全月着陆区选址方法，其特征在于，所述选址方法包括：

S1、对全月范围进行网格划分，计算每个单元网格中所含有的撞击坑个数，并通过计算网格均值获得每个证据因子所对应的网络均值层；

所述网格均值为单元网格内的像素点的均值；

所述证据因子包括月壳厚度、粗糙度图、坡度、数字高程模型、重力梯度、氧化铁含量分布以及月壤光学成熟度；

S2、通过二态变量表示每个网络均值层和目标图层，计算先验几率O_先验，以及每个证据因子在每个单元网格中的权重和后验概率P_后验；

所述目标图层包括Apollo 11、12、13、14、15、16、17号，Luna 16、17、20、21、23、24号，Surveyor 1、2、3、5、6、7号，Chang′E3、Ranger6、7、8、9号，CRAIL及LADEE共50个已着陆点；

S3、基于分形理论，统计出全部单元网格中含有多于S个撞击坑的单元网格的总数A，通过log(S)-log(A)双对数图拟合来确定第一阈值，根据所述第一阈值确定每个单元网格内撞击坑的分布复杂度；

其中，S＝0,1,...,k，以及A＝A₀,A₁,...,A_k；

S4、基于第一单元网格和第二单元网格的索引值，降低第一单元网格的后验概率P_后验，升高第二单元网格的后验概率P_后验，以获得修正后验概率P_修正后验；

所述第一单元网格为撞击坑个数S大于拐点值T的单元网格；

所述第二单元网格为撞击坑个数为0的单元网格；

所述索引值为基于全月范围进行网格划分时，所设定的用于区分不同网格的数值；

S5、通过对所述修正后验概率P_修正后验进行颜色渲染，获得后验概率图，将所述后验概率图与专家筛选图层、撞击坑数据库进行结合，进行全月着陆区的选址。

2.根据权利要求1所述的选址方法，其特征在于，所述S1包括：

S11、将每个证据图层的像素点读成二维矩阵；

所述证据图层为每个证据因子所对应的图层；

S12、将每个证据图层的像素矩阵行列数的公约数作为所划分网格的尺寸，对全月范围内的区域进行网格划分；

所述像素矩阵为单元网格内的像素点所形成的矩阵；

S13、若证据图层的像素点值不能被网格尺寸所整除，则对所述证据图层的边缘像素点进行裁剪；

S14、基于函数cal_avi，遍历各证据图层中每个单元网格内的像素点，计算每个单元网格像素点均值，以获得每个证据因子所对应的网格均值层。

3.根据权利要求1所述的选址方法，其特征在于，所述S2包括：

S21、将每个证据因子所对应的网格均值层表示为0-1变量，以及基于已着陆点坐标，确定每个已着陆点所在的单元网格，用二态变量表示所确定的单元网格；

S22、根据总单元网格的个数和含已着陆点单元网格的个数确定先验概率P_先验，通过证据权重法计算先验几率O_先验，以及每个证据因子在每个单元网格内的权重和后验概率P_后验。

4.根据权利要求3所述的选址方法，其特征在于，所述S21包括：

S211、遍历每个网格均值层中的网格均值，根据第二阈值进行判断，若网格均值大于/小于第二阈值，则将网格均值所对应的单元网格设置为“0/1”或“1/0”；

所述第二阈值为基于专家知识判断所确定的阈值；

S212、将已着陆点坐标范围从-90°至90°调整为0°至180°，将-180°至180°调整为0°至360°；

S213、基于已着陆点坐标，遍历所划分网格，判断单元网格里是否含有已着陆点；通过将存在已着陆点的单元网格设置为“1”，将不含有已着陆点的单元网格设置为“0”，获得已着陆点网格层。