一种浅层岩性地质图生成方法、装置、储存介质及设备与流程

本发明涉及地质勘测技术领域，尤其涉及一种浅层岩性地质图生成方法、装置、储存介质及设备。

背景技术：

我国覆盖区的面积约占全国陆地面积的1/3以上，其中相当一部分地区是覆盖层厚度不超过100m的浅覆盖区(第四系覆盖层厚度小于100m，覆盖层面积占图幅面积50％或50％以上的地区)。在浅覆盖区地质填图工作中，常规填图方法由于基岩露头少而受到限制，隐伏的岩体、地层、构造、矿产(化)以及地表地质现象在深部发生的变化情况等信息无法直接观察研究，利用转石或局部零散露头填制的地质图可信度低、信息量少、整体质量不高。

目前，通过单一勘查方法进行地质填图的工作效率地，准确度及精度较差，尤其是通过单一方法勘查覆盖区，得到的区域地质图件的信息承载量较少，因此，如何提高图件的质量是一件亟待解决的问题。

技术实现要素：

(一)发明目的

本发明的目的是提供一种浅层岩性地质图生成方法、装置、储存介质及设备以解决现有技术中地质图件绘制效率低，图件准确度及精度较差的问题。

(二)技术方案

为解决上述问题，本发明的第一方面提供了一种浅层岩性地质图生成方法，包括：

基于模糊c均值聚类算法对目标区域的多种类型遥感数据进行岩性识别，得到所述多种类型遥感数据所对应的识别结果；

分析所述识别结果的数据之间的相关性，筛选出所述识别结果中最能体现目标区域中岩性的岩性地质图，并反向计算得到与之对应的最能体现目标区域中岩性的遥感数据；

将所述最能体现目标区域中岩性的遥感数据重新进行多类型遥感数据聚类岩性识别，得到目标区域岩性地质图。

进一步地，所述基于模糊c均值聚类算法对目标区域的多种类型遥感数据进行岩性识别，得到所述多种类型遥感数据所对应的识别结果具体为：

假设将遥感数据每个采样点视为一个采集向量，所述向量维度等于遥感数据通道数，每幅遥感数据有n个向量xi(i＝1,2,…,n)，将分为c个模糊组，并求每组的聚类中心，使得非相似性指标的价值函数达到最小，与引入模糊划分相适应，隶属矩阵u允许有取值在0，1间的元素，同时加上归一化规定，一个数据集的隶属度的和总等于1：

那么，基于模糊c均值聚类算法的目标函数为：

这里uij介于0，1间；ci为模糊组i的聚类中心，dij＝||ci-xj||为第i个聚类中心与第j个数据点间的欧几里德距离；且m∈[1,∞)是一个加权指数；

构造如下新的目标函数，可求得使(3.2)式达到最小值的必要条件：

这里λj，j＝1到n，是1)式的n个约束式的拉格朗日乘子，对所有输入参量求导，使式(2)达到最小的必要条件为：

和

由上述两个必要条件，模糊c均值聚类算法的计算步骤如下：

步骤1：用值在0，1间的随机数初始化隶属矩阵u，使其满足式(1)中的约束条件；

步骤2：用式(4)计算c个聚类中心ci，i＝1,…,c；

步骤3：根据式(2)计算价值函数。如果它小于某个确定的阀值，或它相对上次价值函数值的改变量小于某个阀值，则算法停止。

进一步地，所述分析所述识别结果的数据之间的相关性，筛选出所述识别结果中最能体现目标区域中岩性的岩性地质图，并反向计算得到与之对应的最能体现目标区域中岩性的遥感数据包括：

对所述识别结果进行数据间的相关性分析，通过计算得到的多幅岩性地质图的相关性，设得到岩性地质图综述为m，筛选出相关性较差的前n幅图件，所述相关性较差为协方差较大：

其中，x,y表示两幅岩性地质图，xi，yi表示岩性地质图中的第i个点，表示不同地质图的均值，

筛选出预设数量相关性低的岩性地质图组成最能体现目标区域中岩性的岩性地质图；

基于所述最能体现目标区域中岩性的岩性地质图反向计算得到与之对应的最能体现目标区域中岩性的遥感数据。

进一步地，还包括：

将多种类型遥感数据进行预处理后进行拼接，得到目标区域的多种类型遥感数据。

进一步地，所述预处理包括：岩性碎片的均色处理、岩性碎片的边界位置插值和羽化处理。

进一步地，还包括：

对所述目标区域岩性地质图进行空间滤波处理。

根据本发明的另一个方面，提供一种浅层岩性地质图生成装置，包括：

模糊c均值聚类模块，用于基于模糊c均值聚类算法对目标区域的多种类型遥感数据进行岩性识别，得到所述多种类型遥感数据所对应的识别结果；

数据筛选模块，用于分析所述识别结果的数据之间的相关性，筛选出所述识别结果中最能体现目标区域中岩性的岩性地质图，并反向计算得到与之对应的最能体现目标区域中岩性的遥感数据；

聚类模块，用于将所述最能体现目标区域中岩性的遥感数据重新进行多类型遥感数据聚类岩性识别，得到目标区域岩性地质图。

进一步地，所述基于模糊c均值聚类算法对目标区域的多种类型遥感数据进行岩性识别，得到所述多种类型遥感数据所对应的识别结果具体为：

假设将遥感数据每个采样点视为一个采集向量，所述向量维度等于遥感数据通道数，每幅遥感数据有n个向量xi(i＝1,2,…,n)，将分为c个模糊组，并求每组的聚类中心，使得非相似性指标的价值函数达到最小，与引入模糊划分相适应，隶属矩阵u允许有取值在0，1间的元素，同时加上归一化规定，一个数据集的隶属度的和总等于1：

那么，基于模糊c均值聚类算法的目标函数为：

这里uij介于0，1间；ci为模糊组i的聚类中心，dij＝||ci-xj||为第i个聚类中心与第j个数据点间的欧几里德距离；且m∈[1,∞)是一个加权指数；

构造如下新的目标函数，可求得使(3.2)式达到最小值的必要条件：

这里λj，j＝1到n，是1)式的n个约束式的拉格朗日乘子，对所有输入参量求导，使式(2)达到最小的必要条件为：

和

由上述两个必要条件，模糊c均值聚类算法的计算步骤如下：

步骤1：用值在0，1间的随机数初始化隶属矩阵u，使其满足式(1)中的约束条件；

步骤2：用式(4)计算c个聚类中心ci，i＝1,…,c；

步骤3：根据式(2)计算价值函数。如果它小于某个确定的阀值，或它相对上次价值函数值的改变量小于某个阀值，则算法停止。

进一步地，所述分析所述识别结果的数据之间的相关性，筛选出所述识别结果中最能体现目标区域中岩性的岩性地质图，并反向计算得到与之对应的最能体现目标区域中岩性的遥感数据包括：

相关性分析单元，用于对所述识别结果进行数据间的相关性分析；

筛选单元，用于筛选出预设数量相关性低的岩性地质图组成最能体现目标区域中岩性的岩性地质图；

反向计算模块，用于基于所述最能体现目标区域中岩性的岩性地质图反向计算得到与之对应的最能体现目标区域中岩性的遥感数据。

进一步地，还包括：

拼接模块，用于将多种类型遥感数据进行预处理后进行拼接，得到目标区域的多种类型遥感数据。

进一步地，所述预处理包括：岩性碎片的均色处理、岩性碎片的边界位置插值和羽化处理。

进一步地，还包括：

空间滤波模块，用于对所述目标区域岩性地质图进行空间滤波处理。

根据本发明的又一方面，提供一种计算机存储介质，所述存储介质上存储有计算机程序，所述程序被处理器执行时实现上述技术方案中任意一项所述方法的步骤。

根据本发明的又一方面，提供一种设备，包括存储器、显示器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述程序时实现上述技术方案中任意一项所述方法的步骤。

(三)有益效果

本发明的上述技术方案具有如下有益的技术效果：

通过多种勘查结果的综合应用，筛选各种方法的优质结果进行拼接，从而形成高质量的地质图，提高地质填图的工作效率，丰富地质图，尤其是覆盖区区域地质图件的信息承载量，提高图件的质量和基础地质研究的水平。

附图说明

图1是根据本发明第一实施方式的浅层岩性地质图生成方法流程图；

图2是根据本发明一可选实施方式的高分二号遥感图像机器识别的对比图；

图3是根据本发明一可选实施方式的spot数据及聚类结果图；

图4是根据本发明一可选实施方式的对spot遥感数据进行机器识别处理后的数据及聚类结果图；

图5是根据本发明一可选实施方式的aster传感器数据及聚类结果图；

图6是根据本发明一可选实施方式的对aster遥感数据进行机器识别处理后的数据及聚类结果图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明了，下面结合具体实施方式并参照附图，对本发明进一步详细说明。应该理解，这些描述只是示例性的，而并非要限制本发明的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本发明的概念。

显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

此外，下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。

如图1所示，在本发明实施例的第一方面，提供了一种浅层岩性地质图生成方法，包括：

s1：基于模糊c均值聚类算法对目标区域的多种类型遥感数据进行岩性识别，得到所述多种类型遥感数据所对应的识别结果；

s2：分析所述识别结果的数据之间的相关性，筛选出所述识别结果中最能体现目标区域中岩性的岩性地质图，根据所述岩性地质图筛选出对应的遥感数据，形成新的最能体现目标区域中岩性的遥感数据；此步骤筛选出的遥感数据，可以形成一个新的聚类分析样本，具体通过协方差相关的方法，筛选出相关度较低的遥感数据，形成新的聚类分析样本集。

s3：将所述最能体现目标区域中岩性的遥感数据重新进行多类型遥感数据聚类岩性识别，得到目标区域岩性地质图。

上述实施例方法通过多种类型遥感数据的综合应用，筛选各种方法的优质结果进行拼接，从而形成高质量的地质图，提高地质填图的工作效率，丰富地质图，尤其是覆盖区区域地质图件的信息承载量，提高图件的质量和基础地质研究的水平。

可选的，所述基于模糊c均值聚类算法对目标区域的多种类型遥感数据进行岩性识别，得到所述多种类型遥感数据所对应的识别结果具体为：

假设将遥感数据每个采样点视为一个采集向量，所述向量维度等于遥感数据通道数，每幅遥感数据有n个向量xi(i＝1,2,…,n)，将分为c个模糊组，并求每组的聚类中心，使得非相似性指标的价值函数达到最小，与引入模糊划分相适应，隶属矩阵u允许有取值在0，1间的元素，同时加上归一化规定，一个数据集的隶属度的和总等于1：

那么，基于模糊c均值聚类算法的目标函数为：

这里uij介于0，1间；ci为模糊组i的聚类中心，dij＝||ci-xj||为第i个聚类中心与第j个数据点间的欧几里德距离；且m∈[1,∞)是一个加权指数；

构造如下新的目标函数，可求得使(3.2)式达到最小值的必要条件：

这里λj，j＝1到n，是1)式的n个约束式的拉格朗日乘子，对所有输入参量求导，使式(2)达到最小的必要条件为：

和

由上述两个必要条件，模糊c均值聚类算法的计算步骤如下：

步骤1：用值在0，1间的随机数初始化隶属矩阵u，使其满足式(1)中的约束条件；

步骤2：用式(4)计算c个聚类中心ci，i＝1,…,c；

步骤3：根据式(2)计算价值函数。如果它小于某个确定的阀值，或它相对上次价值函数值的改变量小于某个阀值，则算法停止。

模糊聚类算法是一种基于函数最优方法的聚类算法，使用微积分计算技术求最优代价函数。在基于概率算法的聚类方法中将使用概率密度函数，为此要假定合适的模型，模糊聚类算法的向量可以同时属于多个聚类，从而摆脱上述问题。在模糊聚类算法中，定义了向量与聚类之间的近邻函数，并且聚类中向量的隶属度由隶属函数集合提供。对模糊方法而言，在不同聚类中的向量隶属函数值是相互关联的。

fcm把n个向量xi(i＝1,2,…,n)分为c个模糊组，并求每组的聚类中心，使得非相似性指标的价值函数达到最小。fcm与hcm的主要区别在于fcm用模糊划分，使得每个给定数据点用值在0，1间的隶属度来确定其属于各个组的程度。与引入模糊划分相适应，隶属矩阵u允许有取值在0，1间的元素。不过，加上归一化规定，一个数据集的隶属度的和总等于1：

那么，fcm的价值函数(或目标函数)就是：

这里uij介于0，1间；ci为模糊组i的聚类中心，dij＝||ci-xj||为第i个聚类中心与第j个数据点间的欧几里德距离；且m∈[1,∞)是一个加权指数。

构造如下新的目标函数，可求得使(5.2)式达到最小值的必要条件：

这里λj，j＝1到n，是(5.1)式的n个约束式的拉格朗日乘子。对所有输入参量求导，使式(5.2)达到最小的必要条件为：

和

由上述两个必要条件，模糊c均值聚类算法是一个简单的迭代过程。在批处理方式运行时，fcm用下列步骤确定聚类中心ci和隶属矩阵u[1]：

步骤1：用值在0，1间的随机数初始化隶属矩阵u，使其满足式(5.1)中的约束条件

步骤2：用式(5.4)计算c个聚类中心ci，i＝1,…,c。

步骤3：根据式(5.2)计算价值函数。如果它小于某个确定的阀值，或它相对上次价值函数值的改变量小于某个阀值，则算法停止。

步骤4：用(5.5)计算新的u矩阵。返回步骤2。

上述算法也可以先初始化聚类中心，然后再执行迭代过程。由于不能确保fcm收敛于一个最优解。算法的性能依赖于初始聚类中心。因此，我们要么用另外的快速算法确定初始聚类中心，要么每次用不同的初始聚类中心启动该算法，多次运行fcm。

设被分类的对象的集合为：x＝{x1，x2，…，xn}，其中每一个对象xk有n个特性指标，设为xk＝(x1k，x2k，…，xnk)t,如果要把x分成c类，则它的每一个分类结果都对应一个c×n阶的boolean矩阵u＝[uik]c×n，对应的模糊c划分空间为：

在此空间上，模糊c均值算法如下：

(dik)²＝(xk-pi^(l))^ta(xk-pi),1≤i≤c,1≤k≤n

1≤k≤n

(dik)²＞0；i＝1，2，…，c

||u^(l)-u^(l-1)||<ε

进一步地，所述分析所述识别结果的数据之间的相关性，筛选出所述识别结果中最能体现目标区域中岩性的岩性地质图，并反向计算得到与之对应的最能体现目标区域中岩性的遥感数据包括：

对所述识别结果进行数据间的相关性分析，通过计算得到的多幅岩性地质图的相关性，设得到岩性地质图综述为m，筛选出相关性较差的前n幅图件，所述相关性较差为协方差较大：

其中，x,y表示两幅岩性地质图，xi，yi表示岩性地质图中的第i个点，表示不同地质图的均值，

筛选出预设数量相关性低的岩性地质图组成最能体现目标区域中岩性的岩性地质图；

基于所述最能体现目标区域中岩性的岩性地质图反向计算得到与之对应的最能体现目标区域中岩性的遥感数据。

可选的，分析岩性地质图集，筛选出最能体现目标区域中岩性的岩性地质图，得到目标区域的岩性地质图优选集包括：对岩性地质图集中每套地质图进行解析，综合多套地质图分析目标区域所包含的岩性；筛选出最能体现目标区域所包含的岩性的岩性地质图，得到目标区域的岩性地质图优选集。

可选的，将岩性地质图优选集中岩性地质图进行拼接，得到目标区域岩性地质图包括：将岩性地质图优选集中每幅图所最能体现岩性的部分进行裁剪，得到目标区域的岩性碎片集；将岩性碎片集中的岩性碎片进行预处理，使岩性碎片色谱统一，得到统一色谱岩性碎片集；将统一色谱岩性碎片集中的岩性碎片按照目标区域的岩性分布进行拼接、融合，得到目标区域岩性地质图。

可选的，预处理包括：岩性碎片的均色处理、岩性碎片的边界位置插值和羽化处理。

可选的，还包括：对目标区域岩性地质图进行空间滤波处理，减低地表覆盖物对浅层岩性识别的影响。

在本发明一具体实施例中，基于多遥感数据的岩性识别分析,通过上述模糊聚类算法应用，提取遥感数据中反映岩性的信息并对识别结果进行分类。

1)高分二号遥感数据聚类分析

如图2所示，对高分二号遥感数据进行机器识别处理，并进行聚类分析得到如下分类结果：

(1)class5：城镇等建筑物的典型特征。遥感图像为亮白色(如图5所示)，说明高分二号卫星全波段的反射均较好，即使通过低通滤波处理(提高深层地层分辨能力)，可见光波段数据仍然很难消除此类地表干扰的影响。

(2)class2和class3：植被覆盖的典型特征。遥感图像为绿色(如图5所示)，分类结果以浅绿色和深绿色表示。从分类结果可以看到，道路的干扰信息已经被有效的消除。

(3)class4：聚类分析结果中白色点虚线圈画的深蓝色区域在遥感图像和遥感处理图像上表现为植被覆盖，但是区别于class2和class3的分类结果，实际近地表岩性显示为黄土或黏土为主。对比现有地质填图结果，此区域多为湖相黏土和风成黄土，聚类结果预测效果较好。

(4)class1：在遥感图像上为植被覆盖，区别于其他植被的分类结果，聚类为红色，对比现有地质填图结果，此区域多为花岗闪长岩或砂砾岩为主，聚类结果预测效果较好。

高分二号数据能够很好的识别运城地区北部花岗岩分布规律，对北部风成黄土和中部砂岩分布条带均有较好识别，同时，对区内黏土的分布位置也较为准确，但单一遥感数据存在一定的多解性，需要通过多套数据的综合分析对近地表岩性进行更为精细的分析。

2)spot卫星遥感数据聚类分析

spot卫星数据受到云层干扰，如图3所示，对于不同波段的聚类分析结果具有不同程度的影响，通过对比，近红外遥感数据的聚类分析结果能够更好的减弱云层遮蔽的影响。

如图5-6所示，对spot遥感数据进行机器识别处理，并进行聚类分析得到如下分类结果，如图4所示

(1)近红外波段的聚类结果能够较好的反映工区北部(1)花岗岩和(2)砂砾石冲积水道及冲积扇的分布范围，对比现有地质填图，白垩系花岗岩及沉积水道岩性与聚类结果吻合，且具有较高分辨能力，尤其对冲积水道分布位置及形态表征优于高分二号遥感数据分类结果。

(2)近红外波段的聚类结果能够较好的反映出北部冲积扇体的平面分布，与现有的地质填图所揭示的湖相沉积泥岩结果吻合较好。工区东北部水道向南沉积，形成冲积扇相沉积，与西北部的冲积扇相被分为一类(深蓝色)，说明近红外波段的聚类结果能够较好的反映河流相及冲积扇相沉积岩。

(3)近红外波段的聚类结果能够较好的反映出工区的中西部和南部风成黄土(红色和黄色混合)，与现有的地质填图所揭示的结果较为一致，但是此聚类特征存在一定的多解性，有待补充数据进行进一步分析。

(4)近红外波段的聚类结果能够较好的反映工区中部粉砂岩条带分布特征，与现有地质填图结果一致。聚类结果中条带中西部为黄色(岩性为粉砂岩)，东部为蓝色，钻井显示为取土井，岩性在条带上的变化与聚类结果一致，说明聚类结果较为准确。

spot数据近红外波段聚类结果受地表植被的影响较小，相比于高分二号数据，能够更好的反映出工区北部河流相和冲积扇相岩性特征；对工区中部砂岩条带及其中的岩性变化识别较为准确；对工区的中西部和南部的风成黄土有较好的识别，聚类识别结果与现有的地质填图结果较为一致。

3)aster传感器遥感数据聚类分析

aster传感器遥感数据波段范围较宽，包括：(1)可见光/近红外；(2)短波红外；(3)热红外，通过对比不同波段的遥感数据，如图5所示，热红外波段与全波段数据聚类分析结果基本一致，全波段能反映出更多细节和岩性信息。

如图6所示，对spotaster遥感数据进行机器识别处理，并进行聚类分析得到如下分类结果：

(1)全波段遥感数据包括近红外短波长信号，聚类结果能够较好的反映出工区北部花岗岩、河流相和冲积扇相岩性分布，短波长信号可以反映出更为细节的岩性信息，此预测结果与spot遥感数据和现有岩性填图结果一致。

(2)全波段遥感数据包含的热红外波段波场数据波场较长，能够更好的反映更为深层的岩性信息。聚类结果能够较好的反映出工区中部偏北的砂体(红色)，此预测结果与现有的岩性填图结果吻合较好。同时，聚类结果能够较好反映工区中部条带状砂体分布(红色)，在条带东部取土井位置，聚类结果表现为绿色，说明岩性横向变化，此结果与现有的岩性填图结果一致。

(3)全波段遥感数据能够较好的反映出工区南部湖相沉积黏土。此举类结果与现有岩性填图结果分布范围较为一致，说明长波长数据能够较好的压制地表植被的影响。

aster传感器数据能够很好的识别运城地区北部花岗岩分布规律，对北部河流相及冲击扇相沉积体有较好的聚类效果。同时，聚类结果能够反映出中部偏北砂岩及中部条带分布砂体的特征。从全波段数据聚类结果分析，aster传感器热红外波段对浅地表水体的分布较为敏感，会造成工区南北湖相沉积泥岩的聚类特征不同。

遥感数据进行岩性识别，需要对数据进行预处理，项目中提出的均色处理及无缝拼接处理方法能够较好的完成数据的一致性处理。应用遥感数据进行近地表岩性识别前，需要针对不同类型遥感数据进行多尺度空间滤波处理，使得数据能够较好的反映近地表岩性特征，减弱地表干扰。项目提出的遥感数据模糊聚类方法能够较好的识别出近地表岩性信息，不同类型遥感数据所反映的岩性信息不同，所以通过多套遥感数据的综合应用，浅地表岩性分布情况可以得到更为准确的预测。

在本发明实施例的另一个方面，提供一种浅层岩性地质图生成装置，包括：

模糊c均值聚类模块，用于基于模糊c均值聚类算法对目标区域的多种类型遥感数据进行岩性识别，得到所述多种类型遥感数据所对应的识别结果；

数据筛选模块，用于分析所述识别结果的数据之间的相关性，筛选出所述识别结果中最能体现目标区域中岩性的岩性地质图，并反向计算得到与之对应的最能体现目标区域中岩性的遥感数据；

聚类模块，用于将所述最能体现目标区域中岩性的遥感数据重新进行多类型遥感数据聚类岩性识别，得到目标区域岩性地质图。

可选的，所述基于模糊c均值聚类算法对目标区域的多种类型遥感数据进行岩性识别，得到所述多种类型遥感数据所对应的识别结果具体为：

假设将遥感数据每个采样点视为一个采集向量，所述向量维度等于遥感数据通道数，每幅遥感数据有n个向量xi(i＝1,2,…,n)，将分为c个模糊组，并求每组的聚类中心，使得非相似性指标的价值函数达到最小，与引入模糊划分相适应，隶属矩阵u允许有取值在0，1间的元素，同时加上归一化规定，一个数据集的隶属度的和总等于1：

那么，基于模糊c均值聚类算法的目标函数为：

这里uij介于0，1间；ci为模糊组i的聚类中心，dij＝||ci-xj||为第i个聚类中心与第j个数据点间的欧几里德距离；且m∈[1,∞)是一个加权指数；

构造如下新的目标函数，可求得使(3.2)式达到最小值的必要条件：

这里λj，j＝1到n，是1)式的n个约束式的拉格朗日乘子，对所有输入参量求导，使式(2)达到最小的必要条件为：

和

由上述两个必要条件，模糊c均值聚类算法的计算步骤如下：

步骤1：用值在0，1间的随机数初始化隶属矩阵u，使其满足式(1)中的约束条件；

步骤2：用式(4)计算c个聚类中心ci，i＝1,…,c；

步骤3：根据式(2)计算价值函数。如果它小于某个确定的阀值，或它相对上次价值函数值的改变量小于某个阀值，则算法停止。

可选的，所述分析所述识别结果的数据之间的相关性，筛选出所述识别结果中最能体现目标区域中岩性的岩性地质图，并反向计算得到与之对应的最能体现目标区域中岩性的遥感数据包括：

相关性分析单元，用于对所述识别结果进行数据间的相关性分析；

筛选单元，用于筛选出预设数量相关性低的岩性地质图组成最能体现目标区域中岩性的岩性地质图；

反向计算模块，用于基于所述最能体现目标区域中岩性的岩性地质图反向计算得到与之对应的最能体现目标区域中岩性的遥感数据。

可选的，还包括：

拼接模块，用于将多种类型遥感数据进行预处理后进行拼接，得到目标区域的多种类型遥感数据。

可选的，所述预处理包括：岩性碎片的均色处理、岩性碎片的边界位置插值和羽化处理。

可选的，还包括：

空间滤波模块，用于对所述目标区域岩性地质图进行空间滤波处理。

在本发明实施例的又一方面，提供一种计算机存储介质，存储介质上存储有计算机程序，程序被处理器执行时实现上述实施例中任意一项方法的步骤。

在本发明实施例的又一方面，提供一种设备，包括存储器、显示器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行程序时实现上述实施例中任意一项方法的步骤。

本发明旨在保护一种浅层岩性地质图生成方法，包括：基于模糊c均值聚类算法对目标区域的多种类型遥感数据进行岩性识别，得到目标区域岩性地质图集；分析岩性地质图集，筛选出最能体现目标区域中岩性的岩性地质图，得到目标区域的岩性地质图优选集；将岩性地质图优选集中岩性地质图进行拼接，得到目标区域岩性地质图。该方法通过多种类型遥感数据的综合应用，筛选各种方法的优质结果进行拼接，从而形成高质量的地质图，提高地质填图的工作效率，丰富地质图，尤其是覆盖区区域地质图件的信息承载量，提高图件的质量和基础地质研究的水平。

应当理解的是，本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理，而不构成对本发明的限制。因此，在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。此外，本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。