[发明专利]钻孔数据抽稀方法和装置、设备及储存介质

说明书

本申请公开了一种钻孔数据抽稀方法，其方法包括如下步骤：获取钻孔数据集；其中，钻孔数据集包括多个钻孔数据及每个钻孔数据所对应的显示图层；显示图层的个数为m；对钻孔数据进行经纬度坐标换算，得到钻孔数据的像素索引；基于每一显示图层下钻孔数据的像素索引，生成对应显示图层下的抽稀格网；其中，抽稀格网的行数为p行，列数为q列；遍历每一显示图层下的抽稀格网中的每个子格，若子格内存在多个钻孔数据时，由多个钻孔数据中选取一个作为保留数据，得到抽稀处理后的钻孔数据。本公开的方法提高了处理速度，抽稀后的钻孔数据分布均匀，并且提高了抽稀后的钻孔数据的真实性。

技术领域

本公开涉及地质工程领域，尤其涉及一种钻孔数据抽稀方法和装置、设备及储存介质。

背景技术

近年来，我国各城市都开展了大量的城建工程、地下工程、地下水资源和矿产资源的勘查工作，积累了海量的地下地质结构与资源数据。随着数据规模越来越大，建模精度要求也越来越高，其中某些城市或地区单是工程地质钻孔数据就可达到百万个，面对如此庞大的数据，系统现有的海量数据的管理与操作能力效率始终不尽如人意，难以满足实际应用中的实时绘制和交互处理的需求。不能保留关键地形特征，真实性低。

发明内容

有鉴于此，本公开提出了一种钻孔数据抽稀方法，包括：

获取钻孔数据集；其中，所述钻孔数据集包括多个钻孔数据及每个所述钻孔数据所对应的显示图层；所述显示图层的个数为m；

对所述钻孔数据进行经纬度坐标换算，得到所述钻孔数据的像素索引；

基于每一显示图层下所述钻孔数据的像素索引，生成对应显示图层下的抽稀格网；其中，所述抽稀格网的行数为p行，列数为q列；

遍历每一显示图层下的所述抽稀格网中的每个子格，若所述子格内存在多个钻孔数据时，由所述多个钻孔数据中选取一个作为保留数据，得到抽稀处理后的钻孔数据。

在一种可能实现的方式中，所述m的取值范围为：0＜m≤22。

在一种可能实现的方式中，对所述钻孔数据进行经纬度坐标换算，得到所述钻孔数据的像素索引，包括：

读取x像素索引转换公式：

x＝{(lon+180)×256×(2ⁿ-1)/360}

利用x像素索引转换公式将每一所述显示图层下的钻孔数据经度坐标换算为像素索引；

读取y像素索引转换公式：

利用y像素索引转换公式将每一所述显示图层下的钻孔数据纬度坐标换算为像素索引；

其中，lon表征经度坐标，lat表征纬度坐标，n表征缩放的图层，0≤n≤m。

在一种可能实现的方式中，基于每一显示图层下所述钻孔数据的像素索引，生成对应显示图层下的抽稀格网，包括：

由每一所述显示图层下所述钻孔数据的像素索引中查找最小横向索引值x1、最大横向索引值x2、最小纵向索引值y1和最大纵向索引值y2；

以(x1，y1)作为左上角坐标，(x2，y2)作为右下角坐标生成每一所述显示图层所对应的包络矩形；

将每一所述显示图层所对应的所述包络矩形划分为p行，q列后，得到每一所述显示图层所对应的所述抽稀格网。

在一种可能实现的方式中，将每一所述显示图层所对应的所述包络矩形划分为p行，q列时，根据以下计算公式获取行数p和列数q：

其中，n为缩放的图层。