CN114581556B - 一种区域地质调查中的数字填图方法 - Google Patents

Abstract

本发明公开了一种区域地质调查中的数字填图方法，包括以下步骤：步骤S1、获取目标区域原有的地质调查报告作为基础地质报告，以及获取多个调查组在目标区域获得的多个现有地质调查报告作为多个更新地质报告，并在所述基础地质报告和多个更新地质报告中提取出填图数据；步骤S2、将多个更新地质报告的填图数据分别与基础地质报告的填图数据进行差异性量化，并基于所述差异性的量化结果确定目标区域的更新概率，再基于所述更新概率确定多个更新地质报告的填图数据的更新权重。本发明对加权后的基础地质报告的填图数据和多个更新地质报告的填图数据进行求和得到高准确性的目标区域的更新填图数据，实现在更新目标区域填图数据时提高泛化力。

Description

一种区域地质调查中的数字填图方法

技术领域

本发明涉及地质填图技术领域，具体涉及一种区域地质调查中的数字填图方法。

背景技术

地质剖面图是地质内容最基本的表达形式之一，是将各种勘查工程(钻探、坑探等)及不同岩层、岩体、构造形态等地质和矿产现象沿不同方位二维投影，并用不同的岩性花纹、符号及各种注记绘制的二维图形。地质剖面图是沿某一选定方向，垂直于水平面的二维平面图(或投影)。地质平剖面图对分析地层结构具有重要意义，是整个地质工作成果的基础和先导图件，可以整体地反映所调查区域的地质状况，发挥着指导地质勘查、成矿预测、矿山设计与生产决策等的重要作用。传统的地质信息的表达方式主要有两种：一种是采用平面图和剖面图来表达,其实质就是将三维地质环境中的地质现象投影到某一平面(XY平面、XZ平面或YZ平面)上进行表达；另一种是采用透视和轴侧投影原理,对三维地质环境中的地质现象进行透视制图,或是将它们投影到两个以上的平面上进行组合表达,以增强三维视觉效果，提高人们对目标体的三维理解。而二维平面图或剖面图是构建三维立体地质信息表达的基础，所以针对地质剖面图的填图研究对后续研究三维建模大有裨益。

近年来，自动化成图的研究愈发得到研究者们的青睐，地质剖面图的自动数字化也在不断发展中。当前的地质剖面图数字化填图多以钻孔数据为主，数据收集较为困难，且生成的图像不能编辑，不能进行一些必要的空间分析。目前，地质图件构图多以数据驱动构建，数据驱动的建模方法是指根据建模数据直接建模的方法。按照数据源的不同，数据驱动的建模方法可分为基于钻孔、剖面、平面地质图和多源数据融合的建模方法。相较于数据驱动的建模方法，知识驱动的建模方法在建模过程中增加了地质知识的应用，现有技术的数据驱动方式，在更新现有填图数据时缺乏更新泛化性，更新精度低。

发明内容

本发明的目的在于提供一种区域地质调查中的数字填图方法，以解决现有技术中在更新现有填图数据时缺乏更新泛化性，更新精度低的技术问题。

为解决上述技术问题，本发明具体提供下述技术方案：

一种区域地质调查中的数字填图方法，包括以下步骤：

步骤S1、获取目标区域原有的地质调查报告作为基础地质报告，以及获取多个调查组在目标区域获得的多个现有地质调查报告作为多个更新地质报告，并在所述基础地质报告和多个更新地质报告中提取出填图数据；

步骤S2、将多个更新地质报告的填图数据分别与基础地质报告的填图数据进行差异性量化，并基于所述差异性的量化结果确定目标区域的更新概率，再基于所述更新概率确定多个更新地质报告的填图数据的更新权重；

步骤S3、设定所述基础地质报告的填图数据的基础权重，将所述基础权重和更新权重分别加权至所述基础地质报告的填图数据和多个更新地质报告的填图数据上，并对加权后的基础地质报告的填图数据和多个更新地质报告的填图数据进行求和得到高准确性的目标区域的更新填图数据；

步骤S4、根据更新填图数据提取目标区域的地质剖面线，通过反距离权重法确定地质剖面线上内插点的高程值，并用最小二乘法多项式对所述高程值进行拟合出所述地质剖面线，根据填图数据和地质剖面线确定出各岩层界面线，根据地质剖面线和各岩层界面线依次对象化各岩层得到多边形岩层，以及构建出多边形岩层的拓扑关系，再在所述多边形岩层添加更新填图数据，以完成地质剖面图的填图表达。

作为本发明的一种优选方案，所述在所述基础地质报告和多个更新地质报告中提取出填图数据，包括：

依次截取基础地质报告、更新地质报告中所有描述目标区域地质特征的语义文本，并对所述语义文本利用正则表达式过滤噪声文本，所述正则表达式由填图数据的撰写格式决定，所述填图数据包括目标区域的经纬度数据，以及各岩层的属性数据、倾角数据和厚度数据。

作为本发明的一种优选方案，所述将多个更新地质报告的填图数据分别与基础地质报告的填图数据进行差异性量化，包括：

对所述多个更新地质报告的填图数据和基础地质报告的填图数据进行标准化处理以消除量纲误差，并依次计算多个更新地质报告的填图数据与基础地质报告的填图数据的欧式距离；

基于所述欧式距离构建表征所述更新地质报告的填图数据与基础地质报告差异性的量化公式，所述差异性的量化公式为：

式中，J_i表征为第i个更新地质报告的填图数据与基础地质报告的填图数据的差异性，x_i,j表征为第i个更新地质报告的填图数据中第j个数据分量，x_o,j表征为基础地质报告的填图数据中第j个数据分量，n_i表征为第i个更新地质报告的填图数据中包含的数据分量总数，i，j均表征为计量常数。

作为本发明的一种优选方案，所述基于所述差异性的量化结果确定目标区域的更新概率，包括：

设定差异性阈值，并将每个更新地质报告的填图数据与基础地质报告的填图数据的差异性与差异性阈值比较，其中，

若所述差异性高于或等于差异性阈值，将对应的更新地质报告标记为有效更新报告；

若所述差异性低于差异性阈值，将对应的更新地质报告标记为非有效更新报告；

统计有效更新报告总数，并将有效更新报告总数与更新地质报告总数的比值作为所述目标区域的更新概率，所述更新概率的计算公式为：

式中，P表征为更新概率，k表征为有效更新报告总数，m表征为更新地质报告总数。

作为本发明的一种优选方案，所述基于所述更新概率确定多个更新地质报告的填图数据的更新权重，包括：

设定更新概率阈值，其中，

若更新概率高于或等于更新概率阈值，则将多个更新地质报告的填图数据的更新权重设置为：

式中，p_i表征为第i个更新地质报告的填图数据的更新权重；

若更新概率低于更新概率阈值，则将多个更新地质报告的填图数据的更新权重设置为0。

作为本发明的一种优选方案，所述设定所述基础地质报告的填图数据的基础权重，包括：

若多个更新地质报告的填图数据的更新权重均为0，则将所述基础地质报告的填图数据的基础权重设置为1；

若多个更新地质报告的填图数据的更新权重不均为0，则将所述基础地质报告的填图数据的基础权重设置为：

式中，p_o表征为所述基础权重；

优选的，基于所述基础权重对所述更新权重进行归一化处理，所述归一化处理公式为：

式中，p_i′表征为归一化处理后的第i个更新地质报告的填图数据的更新权重。

作为本发明的一种优选方案，所述将所述基础权重和更新权重分别加权至所述基础地质报告的填图数据和多个更新地质报告的填图数据上，包括：

将所述基础权重加权至所述基础地质报告的填图数据上得到加权后的基础地质报告的填图数据为：

x_o，j′＝p₀*x_o，j；

将所述更新权重加权至所述更新地质报告的填图数据上得到加权后的更新地质报告的填图数据为：

x_i，j′＝p_i′*x_i，j；

式中，x_o,j′表征为加权后的基础地质报告的填图数据中第j个数据分量，x_i,j′表征为加权后的第i个更新地质报告的填图数据中第j个数据分量。

作为本发明的一种优选方案，所述对加权后的基础地质报告的填图数据和多个更新地质报告的填图数据进行求和得到高准确性的目标区域的更新填图数据，包括：

所述更新填图数据的计算公式为：

式中，y_o,j表征为所述更新填图数据中第j个数据分量。

作为本发明的一种优选方案，所述根据更新填图数据提取目标区域的地质剖面线，通过反距离权重法确定地质剖面线上内插点的高程值，包括：

利用目标区域的经纬度数据获取目标区域的DEM数据，直线连接相应经纬度点作为目标区域的地质剖面图的起点A和终点B，并将起点A和终点B间的直线等距离划分得到多个端点D_r(r∈[1,M])，M表征为端点总数；

通过反距离权重插值法获得各个端点D_r(r∈[1,M])的高程值。

作为本发明的一种优选方案，所述用最小二乘法多项式对所述高程值进行拟合出所述地质剖面线，包括：

以相邻俩端点间的距离为x_r坐标，所述端点的高程值作为y_r坐标，通过最小二乘法多项式拟合坐标为{(x_r,y_r)|r∈[1,M]}的多个散点得到最优拟合函数；

将最优拟合函数绘制出的曲线作为所述地质剖面线。

本发明与现有技术相比较具有如下有益效果：

本发明获取目标区域原有的地质调查报告作为基础地质报告，以及获取多个调查组在目标区域获得的多个现有地质调查报告作为多个更新地质报告，并在所述基础地质报告和多个更新地质报告中提取出填图数据，确定多个更新地质报告的填图数据的更新权重，对加权后的基础地质报告的填图数据和多个更新地质报告的填图数据进行求和得到高准确性的目标区域的更新填图数据，实现在更新目标区域填图数据时提高泛化力，增强填图精度。

附图说明

为了更清楚地说明本发明的实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是示例性的，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图引伸获得其它的实施附图。

图1为本发明实施例提供的区域地质调查中的数字填图方法流程图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

如图1所示，本发明提供了一种区域地质调查中的数字填图方法，包括以下步骤：

步骤S1、获取目标区域原有的地质调查报告作为基础地质报告，以及获取多个调查组在目标区域获得的多个现有地质调查报告作为多个更新地质报告，并在基础地质报告和多个更新地质报告中提取出填图数据；

在更新目标区域的地质填图时，需要建立在目标区域承认度最高的地质调查报告，比如地质局官方使用的目标区域的地质调查报告，即本实施例中的目标区域原有的地质调查报告，随着调查技术发展或地质变动，需要对目标区域的地址填图进行更新时可派遣多组调查组去重新获取目标的地质调查报告，即本实施例中的多个现有地质调查报告，实现在原有的精度较高的地质调查报告的基础上，融合进多个现有地质调查报告实现地质填图的更新，即本实施例在基础地质报告的融合多个更新地质报告，具体方法如下：

在基础地质报告和多个更新地质报告中提取出填图数据，包括：

依次截取基础地质报告、更新地质报告中所有描述目标区域地质特征的语义文本，并对语义文本利用正则表达式过滤噪声文本，正则表达式由填图数据的撰写格式决定，填图数据包括目标区域的经纬度数据，以及各岩层的属性数据、倾角数据和厚度数据。

在基础地质报告和多个更新地质报告中提取出填图数据，包括：

依次截取基础地质报告、更新地质报告中所有描述目标区域地质特征的语义文本，并对语义文本利用正则表达式过滤噪声文本，正则表达式由填图数据的撰写格式决定，填图数据包括目标区域的经纬度数据，以及各岩层的属性数据、倾角数据和厚度数据。

步骤S2、将多个更新地质报告的填图数据分别与基础地质报告的填图数据进行差异性量化，并基于差异性的量化结果确定目标区域的更新概率，再基于更新概率确定多个更新地质报告的填图数据的更新权重；

将多个更新地质报告的填图数据分别与基础地质报告的填图数据进行差异性量化，包括：

对多个更新地质报告的填图数据和基础地质报告的填图数据进行标准化处理以消除量纲误差，并依次计算多个更新地质报告的填图数据与基础地质报告的填图数据的欧式距离；

基于欧式距离构建表征更新地质报告的填图数据与基础地质报告差异性的量化公式，差异性的量化公式为：

式中，J_i表征为第i个更新地质报告的填图数据与基础地质报告的填图数据的差异性，x_i,j表征为第i个更新地质报告的填图数据中第j个数据分量，x_o,j表征为基础地质报告的填图数据中第j个数据分量，n_i表征为第i个更新地质报告的填图数据中包含的数据分量总数，i，j均表征为计量常数。

进行差异性量化来判定目标区域的基础地质报告是否存在较大的变动，更新地质报告的填图数据与基础地质报告的填图数据的差异性越大，则说明更新地质报告的变动越大，而更新地质报告的填图数据与基础地质报告的填图数据的差异性越小，则说明更新地质报告的变动越小。

基于差异性的量化结果确定目标区域的更新概率，包括：

设定差异性阈值，并将每个更新地质报告的填图数据与基础地质报告的填图数据的差异性与差异性阈值比较，其中，

若差异性高于或等于差异性阈值，将对应的更新地质报告标记为有效更新报告；

若差异性低于差异性阈值，将对应的更新地质报告标记为非有效更新报告；

统计有效更新报告总数，并将有效更新报告总数与更新地质报告总数的比值作为目标区域的更新概率，更新概率的计算公式为：

式中，P表征为更新概率，k表征为有效更新报告总数，m表征为更新地质报告总数。

差异性阈值由使用者进行自行设定，当差异性高于或等于差异性阈值的更新地质报告数量越多，即有效更新报告越多，则说明多数调查组的调查结果是目标区域的地质变动或原有地质调查报告精度不足，因此需要对目标区域的地址填图进行更新的越有必要，原有地质填图不可使用，本实施例利用更新概率进行更新必要性的量化，即更新概率越高；

当差异性高于或等于差异性阈值的更新地质报告数量越小，即有效更新报告越少，则说明多数调查组的调查结果是目标区域的未地质变动或原有地质调查报告精度足够，因此需要对目标区域的地址填图进行更新的越无必要，原有地质填图结果可以继续使用，本实施例利用更新概率进行更新必要性的量化，即更新概率越低。

基于更新概率确定多个更新地质报告的填图数据的更新权重，包括：

设定更新概率阈值，其中，

若更新概率高于或等于更新概率阈值，则将多个更新地质报告的填图数据的更新权重设置为：

式中，p_i表征为第i个更新地质报告的填图数据的更新权重；

若更新概率低于更新概率阈值，则将多个更新地质报告的填图数据的更新权重设置为0。

步骤S3、设定基础地质报告的填图数据的基础权重，将基础权重和更新权重分别加权至基础地质报告的填图数据和多个更新地质报告的填图数据上，并对加权后的基础地质报告的填图数据和多个更新地质报告的填图数据进行求和得到高准确性的目标区域的更新填图数据；

设定基础地质报告的填图数据的基础权重，包括：

若多个更新地质报告的填图数据的更新权重均为0，则将基础地质报告的填图数据的基础权重设置为1；

若多个更新地质报告的填图数据的更新权重不均为0，则将基础地质报告的填图数据的基础权重设置为：

式中，p_o表征为基础权重；

基础权重的值域为(0.5,1)，因此将基础权重用于对基础地质报告的填图数据进行加权可保证更新填图数据中基础地质报告的填图数据的占比范围在(0.5,1)，有效的保障了更新填图数据中对基础地质报告的填图数据对目标区域中基础地质报告的填图数据进行记忆力保留，而且在有效更新报告总数大时，则说明地质填图需要更新的部分越多，即增强了更新地质报告的填图数据在进行地质填图的参考性，削弱了基础地质报告的填图数据在进行地质填图的参考性，因此基础权重降低，并将降低的部分转至更新权重，以覆盖所有更新地质报告的填图数据，使得更新填图数据中包含所有更新地质报告的填图数据的填图记忆力，提高更新填图数据的泛化力，在有效报告总数小时，则说明地质填图需要更新的部分越少，即削弱了更新地质报告的填图数据在进行地质填图的参考性，增强了基础地质报告的填图数据在进行地质填图的参考性，因此更新权重增加，并将降低的部分转至基础权重，提高更新填图数据的记忆力，实现记忆力和泛化力的动态兼顾，得到精度高的更新填图数据进行地质填图。

优选的，基于基础权重对更新权重进行归一化处理，归一化处理公式为：

式中，p_i′表征为归一化处理后的第i个更新地质报告的填图数据的更新权重。

将基础权重和更新权重分别加权至基础地质报告的填图数据和多个更新地质报告的填图数据上，包括：

将基础权重加权至基础地质报告的填图数据上得到加权后的基础地质报告的填图数据为：

x_o，j′＝P_o*x_o，j；

将更新权重加权至更新地质报告的填图数据上得到加权后的更新地质报告的填图数据为：

x_i，j′＝p_i′*x_i，j；

式中，x_o,j′表征为加权后的基础地质报告的填图数据中第j个数据分量，x_i,j′表征为加权后的第i个更新地质报告的填图数据中第j个数据分量。

对加权后的基础地质报告的填图数据和多个更新地质报告的填图数据进行求和得到高准确性的目标区域的更新填图数据，包括：

更新填图数据的计算公式为：

式中，y_o,j表征为更新填图数据中第j个数据分量。

步骤S4、根据更新填图数据提取目标区域的地质剖面线，通过反距离权重法确定地质剖面线上内插点的高程值，并用最小二乘法多项式对高程值进行拟合出地质剖面线，根据填图数据和地质剖面线确定出各岩层界面线，根据地质剖面线和各岩层界面线依次对象化各岩层得到多边形岩层，以及构建出多边形岩层的拓扑关系，再在多边形岩层添加更新填图数据，以完成地质剖面图的填图表达。

根据更新填图数据提取目标区域的地质剖面线，通过反距离权重法确定地质剖面线上内插点的高程值，包括：

利用目标区域的经纬度数据获取目标区域的DEM数据，直线连接相应经纬度点作为目标区域的地质剖面图的起点A和终点B，并将起点A和终点B间的直线等距离划分得到多个端点D_r(r∈[1,M])，M表征为端点总数；

通过反距离权重插值法获得各个端点D_r(r∈[1,M])的高程值。

用最小二乘法多项式对高程值进行拟合出地质剖面线，包括：

以相邻俩端点间的距离为x_r坐标，端点的高程值作为y_r坐标，通过最小二乘法多项式拟合坐标为{(x_r,y_r)|r∈[1,M]}的多个散点得到最优拟合函数；

将最优拟合函数绘制出的曲线作为地质剖面线。

为各岩层界面线自定义设置一个统一长度L，其次，假设存在这样一个点，同时满足以下两个条件：(1)该点到相邻岩层界面线(左侧)的距离，使其满足读取到的岩层真厚度；(2)该点需满足拟合后的地质剖面线多项式方程。如果同时满足条件，则记录该点，否则，更改制图比例尺后重新执行该步骤。然后，依次读取岩层厚度H_l和岩层倾角数据

l表示岩层实体类别编号，并结合设置的岩层界面线长度L绘制该点处的界面线，得到新的各岩层界面线。其中，点到直线距离公式如下：

其中a,b,c分别为前一条界面线(左侧)的直线方程系数，x₀和y₀为满足上述两个条件时记录的点的横坐标和纵坐标。

通过以上步骤后，已构建完整的地质剖面线和各岩层界面线。为实现岩层的对象化编辑和分析，将各岩层界面线终点连接起来，从而形成封闭的多边形连接线，最后，通过矢量线转面的形式实现各岩层形成独立的多边形，并构建拓扑关系以及检查是否存在拓扑错误，若有拓扑错误，则修改拓扑错误后为各个对象化岩层添加倾角信息，以及各岩层的描述信息和厚度信息，从而完成几何形态和属性信息结合的地质剖面图近似表达。

本发明获取目标区域原有的地质调查报告作为基础地质报告，以及获取多个调查组在目标区域获得的多个现有地质调查报告作为多个更新地质报告，并在基础地质报告和多个更新地质报告中提取出填图数据，确定多个更新地质报告的填图数据的更新权重，对加权后的基础地质报告的填图数据和多个更新地质报告的填图数据进行求和得到高准确性的目标区域的更新填图数据，实现在更新目标区域填图数据时提高泛化力，增强填图精度。

以上实施例仅为本申请的示例性实施例，不用于限制本申请，本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内，对本申请做出各种修改或等同替换，这种修改或等同替换也应视为落在本申请的保护范围内。

Claims (10)

1.一种区域地质调查中的数字填图方法，其特征在于，包括以下步骤：

步骤S1、获取目标区域原有的地质调查报告作为基础地质报告，以及获取多个调查组在目标区域获得的多个现有地质调查报告作为多个更新地质报告，并在所述基础地质报告和多个更新地质报告中提取出填图数据；

步骤S2、将多个更新地质报告的填图数据分别与基础地质报告的填图数据进行差异性量化，并基于所述差异性的量化结果确定目标区域的更新概率，再基于所述更新概率确定多个更新地质报告的填图数据的更新权重；

步骤S3、设定所述基础地质报告的填图数据的基础权重，将所述基础权重和更新权重分别加权至所述基础地质报告的填图数据和多个更新地质报告的填图数据上，并对加权后的基础地质报告的填图数据和多个更新地质报告的填图数据进行求和得到高准确性的目标区域的更新填图数据；

步骤S4、根据更新填图数据提取目标区域的地质剖面线，通过反距离权重法确定地质剖面线上内插点的高程值，并用最小二乘法多项式对所述高程值进行拟合出所述地质剖面线，根据填图数据和地质剖面线确定出各岩层界面线，根据地质剖面线和各岩层界面线依次对象化各岩层得到多边形岩层，以及构建出多边形岩层的拓扑关系，再在所述多边形岩层添加更新填图数据，以完成地质剖面图的填图表达。

2.根据权利要求1所述的一种区域地质调查中的数字填图方法，其特征在于：所述在所述基础地质报告和多个更新地质报告中提取出填图数据，包括：

依次截取基础地质报告、更新地质报告中所有描述目标区域地质特征的语义文本，并对所述语义文本利用正则表达式过滤噪声文本，所述正则表达式由填图数据的撰写格式决定，所述填图数据包括目标区域的经纬度数据，以及各岩层的属性数据、倾角数据和厚度数据。

3.根据权利要求2所述的一种区域地质调查中的数字填图方法，其特征在于：所述将多个更新地质报告的填图数据分别与基础地质报告的填图数据进行差异性量化，包括：

对所述多个更新地质报告的填图数据和基础地质报告的填图数据进行标准化处理以消除量纲误差，并依次计算多个更新地质报告的填图数据与基础地质报告的填图数据的欧式距离；

基于所述欧式距离构建表征所述更新地质报告的填图数据与基础地质报告差异性的量化公式，所述差异性的量化公式为：

式中，J_i表征为第i个更新地质报告的填图数据与基础地质报告的填图数据的差异性，x_i,j表征为第i个更新地质报告的填图数据中第j个数据分量，x_o,j表征为基础地质报告的填图数据中第j个数据分量，n_i表征为第i个更新地质报告的填图数据中包含的数据分量总数，i，j均表征为计量常数。

4.根据权利要求3所述的一种区域地质调查中的数字填图方法，其特征在于：所述基于所述差异性的量化结果确定目标区域的更新概率，包括：

设定差异性阈值，并将每个更新地质报告的填图数据与基础地质报告的填图数据的差异性与差异性阈值比较，其中，

若所述差异性高于或等于差异性阈值，将对应的更新地质报告标记为有效更新报告；

若所述差异性低于差异性阈值，将对应的更新地质报告标记为非有效更新报告；

统计有效更新报告总数，并将有效更新报告总数与更新地质报告总数的比值作为所述目标区域的更新概率，所述更新概率的计算公式为：

式中，P表征为更新概率，k表征为有效更新报告总数，m表征为更新地质报告总数。

5.根据权利要求4所述的一种区域地质调查中的数字填图方法，其特征在于：所述基于所述更新概率确定多个更新地质报告的填图数据的更新权重，包括：

设定更新概率阈值，其中，

若更新概率高于或等于更新概率阈值，则将多个更新地质报告的填图数据的更新权重设置为：

式中，p_i表征为第i个更新地质报告的填图数据的更新权重；

若更新概率低于更新概率阈值，则将多个更新地质报告的填图数据的更新权重设置为0。

6.根据权利要求5所述的一种区域地质调查中的数字填图方法，其特征在于：所述设定所述基础地质报告的填图数据的基础权重，包括：

若多个更新地质报告的填图数据的更新权重均为0，则将所述基础地质报告的填图数据的基础权重设置为1；

若多个更新地质报告的填图数据的更新权重不均为0，则将所述基础地质报告的填图数据的基础权重设置为：

式中，p_o表征为所述基础权重；

基于所述基础权重对所述更新权重进行归一化处理，所述归一化处理公式为：

式中，p_i’表征为归一化处理后的第i个更新地质报告的填图数据的更新权重。

7.根据权利要求6所述的一种区域地质调查中的数字填图方法，其特征在于，所述将所述基础权重和更新权重分别加权至所述基础地质报告的填图数据和多个更新地质报告的填图数据上，包括：

将所述基础权重加权至所述基础地质报告的填图数据上得到加权后的基础地质报告的填图数据为：

x_o，j′＝p_o*x_o，j；

将所述更新权重加权至所述更新地质报告的填图数据上得到加权后的更新地质报告的填图数据为：

x_i，j′＝p_i′*x_i，j；

式中，x_o,j’表征为加权后的基础地质报告的填图数据中第j个数据分量，x_i,j’表征为加权后的第i个更新地质报告的填图数据中第j个数据分量。

8.根据权利要求7所述的一种区域地质调查中的数字填图方法，其特征在于，所述对加权后的基础地质报告的填图数据和多个更新地质报告的填图数据进行求和得到高准确性的目标区域的更新填图数据，包括：

所述更新填图数据的计算公式为：

式中，y_o,j表征为所述更新填图数据中第j个数据分量。

9.根据权利要求8所述的一种区域地质调查中的数字填图方法，其特征在于，所述根据更新填图数据提取目标区域的地质剖面线，通过反距离权重法确定地质剖面线上内插点的高程值，包括：

利用目标区域的经纬度数据获取目标区域的DEM数据，直线连接相应经纬度点作为目标区域的地质剖面图的起点A和终点B，并将起点A和终点B间的直线等距离划分得到多个端点D_r(r∈[1,M])，M表征为端点总数；

通过反距离权重插值法获得各个端点D_r(r∈[1,M])的高程值。

10.根据权利要求9所述的一种区域地质调查中的数字填图方法，其特征在于，所述用最小二乘法多项式对所述高程值进行拟合出所述地质剖面线，包括：

以相邻俩端点间的距离为x_r坐标，所述端点的高程值作为y_r坐标，通过最小二乘法多项式拟合坐标为{(x_r,y_r)|r∈[1,M]}的多个散点得到最优拟合函数；

将最优拟合函数绘制出的曲线作为所述地质剖面线。

Images