一种岩芯柱状图生成方法、装置、终端设备及存储介质与流程

1.本发明涉及地质勘察和图像处理技术领域，尤其涉及一种岩芯柱状图生成方法、装置、终端设备及存储介质。


背景技术：

2.在地质勘察领域，需基于现场钻孔岩芯图片生成岩芯柱状图，以使勘察人员基于岩芯柱状图进一步绘制底下全貌；现有的岩芯柱状图生成方法主要由设计人员根据现场钻孔岩芯图片中管道排列情况，人为识别管道间分割线，根据识别的分割线利用制图软件和图像合成软件对管道进行切割，根据切割后的各管道图片拼接成岩芯柱状图；采用上述的方法存在生成效率较低的问题。


技术实现要素：

3.本发明实施例提供一种岩芯柱状图生成方法、装置、终端设备及存储介质，能提高岩芯柱状图生成的效率。
4.本发明一实施例提供一种岩芯柱状图生成方法，包括：
5.获取待处理钻孔岩芯图；
6.对待处理钻孔岩芯图做直线检测；
7.在确定待处理钻孔岩芯图中的直线数量不等于第一阈值时，对待处理钻孔岩芯图的直线进行聚类，生成修复后的钻孔岩芯图，将修复后的钻孔岩芯图中所存在的直线作为分割线，根据修复后的钻孔岩芯图中的分割线对修复后的钻孔岩芯图进行切割，生成若干钻孔岩芯管道图；其中，所述第一阈值根据待处理钻孔岩芯图中管道的数量进行确定；
8.在确定待处理钻孔岩芯图中的直线数量等于第一阈值时，根据待处理钻孔岩芯图的直线对待处理钻孔岩芯图进行分割，生成若干钻孔岩芯管道图；
9.将各钻孔岩芯管道图进行拼接生成岩芯柱状图。
10.进一步地，在对待处理钻孔岩芯图做直线检测之前，还包括：
11.对待处理钻孔岩芯图做灰度化处理，生成第一钻孔岩芯图；
12.对第一钻孔岩芯图做滤波处理，生成第二钻孔岩芯图；
13.对第二钻孔岩芯图做直方图均衡化处理，生成第三钻孔岩芯图；
14.对第三钻孔岩芯图做图像开运算和图像闭运算处理，生成第四钻孔岩芯图；
15.所述对待处理钻孔岩芯图做直线检测，包括：
16.对所述第四钻孔岩芯图做直线检测。
17.进一步地，所述对待处理钻孔岩芯图做直线检测，包括：
18.获取待处理钻孔岩芯图中的所有像素点；
19.对于任意一个像素点，计算待处理钻孔岩芯图中经过像素点的所有直线的正弦函数；
20.根据所有像素点的正弦函数生成正弦函数图像；
21.将正弦函数图像中，正弦函数相交的像素点作为同一直线上的像素点。
22.进一步地，所述对待处理钻孔岩芯图的直线进行聚类，包括：
23.根据k-means算法对待处理钻孔岩芯图的直线进行聚类。
24.进一步地，对第一钻孔岩芯图做滤波处理，生成第二钻孔岩芯图，包括：
25.对第一钻孔岩芯图做高斯滤波处理，生成第二钻孔岩芯图。
26.在上述方法项实施例的基础上，本发明对应提供了装置项实施例；
27.本发明一实施例对应提供了一种岩芯柱状图生成装置，包括：数据获取模块、直线检测模块和直方图生成模块；
28.所述数据获取模块，用于获取待处理钻孔岩芯图；
29.所述直线检测模块，用于对待处理钻孔岩芯图做直线检测；
30.所述直方图生成模块，用于在确定待处理钻孔岩芯图中的直线数量不等于第一阈值时，对待处理钻孔岩芯图的直线进行聚类，生成修复后的钻孔岩芯图，将修复后的钻孔岩芯图中所存在的直线作为分割线，根据修复后的钻孔岩芯图中的分割线对修复后的钻孔岩芯图进行切割，生成若干钻孔岩芯管道图；其中，所述第一阈值根据待处理钻孔岩芯图中管道的数量进行确定；在确定待处理钻孔岩芯图中的直线数量等于第一阈值时，根据待处理钻孔岩芯图的直线对待处理钻孔岩芯图进行分割，生成若干钻孔岩芯管道图；将各钻孔岩芯管道图进行拼接生成岩芯柱状图。
31.本发明另一实施例提供了一种终端设备，包括处理器、存储器以及存储在所述存储器中且被配置为由所述处理器执行的计算机程序，所述处理器执行所述计算机程序时实现上述发明实施例所述的一种岩芯柱状图生成方法。
32.本发明另一实施例提供了一种存储介质，所述存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述存储介质所在设备执行上述发明实施例所述的一种岩芯柱状图生成方法。
33.通过实施本发明具有如下有益效果：本发明提供了一种岩芯柱状图生成方法、装置、终端设备及存储介质，该方法通过获取待处理钻孔岩芯图，对所述钻孔岩芯图中用于分割各管道的直线进行检测，在检测到直线数量不等于由待处理钻孔岩芯图中管道数量确定的第一阈值时，对待处理钻孔岩芯图中的直线进行修复，并以修复后钻孔岩芯图中所存在的直线作为分割线，根据修复后的钻孔岩芯图中的分割线对修复后的钻孔岩芯图进行切割，根据切割图像拼接生成岩芯柱状图；在检测到直线数量等于第一阈值时，根据直线对钻孔岩芯图做切割，并根据切割图像拼接生成岩芯柱状图；通过实施本发明的方法，可以对钻孔岩芯图中的直线进行检测和修复，继而根据直线对钻孔岩芯图中各管道进行分割后拼接生成岩芯柱状图，避免人为处理图像，提高了岩芯柱状图生成的效率。
附图说明
34.图1是本发明一实施例提供的一种岩芯柱状图生成方法的流程示意图。
35.图2是本发明一实施例提供的灰度化处理后的图像。
36.图3是本发明一实施例提供的高斯滤波后的图像。
37.图4是本发明一实施例提供的直方图均衡化处理后的图像。
38.图5是本发明一实施例提供的图像开运算处理后的图像。
39.图6是本发明一实施例提供的sobel滤波器处理后的图像。
40.图7是本发明一实施例提供的图像闭运算处理后的图像。
41.图8是本发明一实施例提供的直线检测时的图像。
42.图9是本发明一实施例提供的一种岩芯柱状图生成装置的结构示意图。
具体实施方式
43.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
44.如图1所示，本发明一实施例提供的一种岩芯柱状图生成方法，包括：
45.步骤s1：获取待处理钻孔岩芯图；
46.步骤s2：对待处理钻孔岩芯图做直线检测；
47.步骤s3：在确定待处理钻孔岩芯图中的直线数量不等于第一阈值时，对待处理钻孔岩芯图的直线进行聚类，生成修复后的钻孔岩芯图，将修复后的钻孔岩芯图中所存在的直线作为分割线，根据修复后的钻孔岩芯图中的分割线对修复后的钻孔岩芯图进行切割，生成若干钻孔岩芯管道图；其中，所述第一阈值根据待处理钻孔岩芯图中管道的数量进行确定；
48.步骤s4：在确定待处理钻孔岩芯图中的直线数量等于第一阈值时，根据待处理钻孔岩芯图的直线对待处理钻孔岩芯图进行分割，生成若干钻孔岩芯管道图；
49.步骤s5：将各钻孔岩芯管道图进行拼接生成岩芯柱状图。
50.需要说明的是，钻孔岩芯图中，通常以管道纵向排布，每一管道包含一定深度的岩芯；例如：在表示地下0-5米的钻孔岩芯图中，包含0-1米岩芯的管道、1-2米岩芯的管道、2-3米岩芯的管道、3-4米岩芯的管道和4-5米岩芯的管道；每一管道以此纵向排列；分割线包括各管道之间的直线以及钻孔岩芯图图像中第一个管道上边缘与最后一个管道下边缘。
51.对于步骤s1、待处理钻孔岩芯图的获取可以通过用户拍摄现场钻孔岩芯图后得到。
52.在一个优先的实施例中，在对待处理钻孔岩芯图做直线检测之前，还包括：对待处理钻孔岩芯图做灰度化处理，生成第一钻孔岩芯图；对第一钻孔岩芯图做滤波处理，生成第二钻孔岩芯图；对第二钻孔岩芯图做直方图均衡化处理，生成第三钻孔岩芯图；对第三钻孔岩芯图做图像开运算和图像闭运算处理，生成第四钻孔岩芯图；所述对待处理钻孔岩芯图做直线检测，包括：对所述第四钻孔岩芯图做直线检测。
53.具体的，将待处理钻孔岩芯图做灰度化处理，将待处理钻孔岩芯图从彩色图像转换为灰度图像。与彩色图像相比，灰度图像能够在保存图片纹理信息的同时，以更小的图像单位信息进行保存，提高后续处理效率。本实施例使用y亮度灰度图像处理方法对待处理钻孔岩芯图做灰度化处理；该灰度化处理方法相对其他灰度化方法，如gamma校正灰度化或最大值灰度化，具有计算量小，处理速度快的优点，能节省对图片做灰度化处理的时间。
54.可通过以下公式对待处理钻孔岩芯图做灰度化处理：
55.y1＝0.2r+0.59g+0.11b
56.其中，y1为灰度化处理后的图像(即上述第一钻孔岩芯图)，一般的彩色图片rgb图片可用一个三维矩阵表示，r为待处理钻孔岩芯图的三维矩阵的红色维度，g为待处理钻孔岩芯图的三维矩阵的绿色维度，b为待处理钻孔岩芯图的三维矩阵的蓝色维度；经过上述y亮度灰度图像处理方法对待处理钻孔岩芯图做灰度化处理后，可得到如图2所示的灰度化处理后的图像(即上述第一钻孔岩芯图)。
57.图像滤波能在尽量保留图像细节特征的条件下对目标图像的噪声像素进行抑制，除去拍摄时产生的自然噪声像素点；为提高后续直线检测的准确性，在得到灰度化处理后的图像后，需对该图像做进一步的滤波处理，以防止图像噪声影响后续直线检测。
58.在一个优选的实施例中，对第一钻孔岩芯图做滤波处理，生成第二钻孔岩芯图，包括：对第一钻孔岩芯图做高斯滤波处理，生成第二钻孔岩芯图。
59.具体的，由于在数字图像采集的过程中最容易出现的是高斯噪声，它是由于不良照明或高温引起的传感器噪声。一般的图像滤波方法有高斯滤波、均值滤波、中值滤波、维纳滤波，但在实际试验中发现高斯滤波其实在处理高斯噪声方面卓有成效，而中值滤波比较适合处理椒盐形式的噪声。同时逆滤波比较适合处理无噪声的运动模糊图像，如果图像既有噪声又有运动模糊，那么维纳滤波是一个很好的选择。本发明的待处理钻孔岩芯图的拍摄场景存在不够明亮，分布不均衡的情况，这种拍摄场景下拍摄的钻孔岩芯图会存在高斯噪声，因此，在本发明中采用针对处理高斯噪声的高斯滤波器对灰度化处理后的钻孔岩芯图做高斯滤波处理；本发明采用卷积核为(2n+1)
×
(2n+1)的高斯滤波器对灰度化处理后的钻孔岩芯图做高斯滤波处理的公式如下：
[0060][0061]
其中，y2(x,y)为对坐标为(x,y)的像素点做高斯滤波处理后的像素点坐标，卷积核大小为2n+1，n可以为任意大小，σ为高斯分布的标准差；σ越大，权重分布越均匀，滤波效果越好，同时图像越模糊；σ越小，权重分布越偏向于窗口中心点，滤波效果越差，同时图像越能保留其原有清晰度。
[0062]
通过上述高斯滤波处理后，可得到如图3所示的高斯滤波后的图像(即上述第二钻孔岩芯图)。
[0063]
在得到高斯滤波后的图像后，需将滤波后的图像做直方图均衡化处理，以灰度块计算的方式将滤波后的图像变换为灰度直方图，并进一步将其构造成相对均匀分布的新图像，从而提高图像在各像素块上的对比度，继而凸显出直线的位置，使后续直线检测结果更准确；可通过以下公式对滤波后的图像做直方图均衡化处理：
[0064][0065]
其中，y3代表当前灰度级经过累积分布函数映射后的图像，n为图像中像素的总和，nj为当前灰度级的像素个数，l是图像中的灰度级总数；
[0066]
通过上述直方图均衡化处理后，可得到如图4所示的直方图图像(即上述第三钻孔岩芯图)。
[0067]
为使后续直线检测的结果更准确，还需对直方图图像(即上述第三钻孔岩芯图)做图像开运算和图像闭运算处理；图像开运算处理对图像先腐蚀后膨胀，能除去图像中孤立的小点和毛刺，同时保持图像中岩芯形状轮廓基本不变，通过图像开运算处理能使得处理后得到的图像在后续进行直线检测时将图像中的直线连在一起；可通过以下公式对直方图图像做图像开运算处理：
[0068]
y4＝dilate(enrode(y3,element))
[0069]
其中，y4为开运算处理后得到的图像，element代表卷积核，enrode为腐蚀操作，dilate为图像的膨胀操作；
[0070]
通过开运算处理后得到的图像如图5所示。
[0071]
由于图像在各个局部的纹理强弱、稀疏度不同，所以在全局上很难找到合适的阈值去分割出边缘，进而准确找到直线的边缘；因此，首先要在各个局部找到边缘线，而后将边缘线连接在一起。具体操作：采用滑动窗口的方法，在每个小窗口内，采用canny算子，取得边缘，而这些边缘会组成多个连通域(连通域是指图像中具有相同像素值并且位置相邻的像素组成的区域)。canny边缘检测的第一步需要用sobel滤波器的步骤来计算每个像素的梯度和方向。sobel滤波器需要两个卷积核，g
x
和gy：
[0072][0073][0074][0075][0076]
其中，每个像素点的梯度大小为g，梯度方向为θ；
[0077]
遍历图像中的所有像素点，判断当前像素点是否是周围像素点中具有相同方向梯度的最大值。如果是梯度最大的像素点，就保留，否则就抑制这一步的目的是将模糊的边界变得清晰，从而实现保留了每个像素点上梯度强度的极大值，继而删掉其他的值。对处理后的图片的梯度大小进行筛选，如果它的梯度大小大于最高阈值，则它会作为边界的像素保留；若它的值低于最低阈值则会把它排除；若它处于最高阈值和最低阈值之间，则它只会在连接到一个高于最高阈值时会被保留。一般在最高阈值和最低阈值的比例为3:1到2:1之间。经过上述处理后可得到如图6所示的图像。
[0078]
此外，由于钻孔岩芯图中的泥土或石头会在边缘检测后产生小的连通域，连通域一般是指图像中具有相同像素值且位置相邻的前景像素点组成的图像区域，这些小的连通域会干扰对钻孔岩芯图中各管道间边界的判断，因此，需要通过对连通域的大小进行筛选，从而除去图像中的小的闭环区域。对上述开运算处理后的图像做图像闭运算操作，图像闭运算操作即为先膨胀后腐蚀，以单通道的黑白图像为例，先通过膨胀操作，将白色区域膨胀扩大，再通过腐蚀操作，将黑色区域腐蚀扩大。一个小的黑色像素点，先通过膨胀操作，被周
围的白色覆盖掉，再进行腐蚀操作，去除小的黑色像素点。图像闭运算能将图像恢复到待处理钻孔岩芯图原图像边缘的粗细，进而提取出整幅图像内小窗口内的边缘；经过上述闭运算处理后得到的图像如图7所示。
[0079]
对于步骤s2、以经过上述一系列处理后的图像作为直线检测时的待处理钻孔岩芯图。
[0080]
在一个优选的实施例中，所述对待处理钻孔岩芯图做直线检测，包括：获取待处理钻孔岩芯图中的所有像素点；对于任意一个像素点，计算待处理钻孔岩芯图中经过像素点的所有直线的正弦函数；根据所有像素点的正弦函数生成正弦函数图像；将正弦函数图像中，正弦函数相交的像素点作为同一直线上的像素点。
[0081]
具体的，在进行直线检测时，采用霍夫直线检测算法对图7闭运算处理后的图像中的直线数量进行检测；霍夫直线检测的原理是建立在任何一个在平面坐标系上的点都可以转化为极坐标内的一个如下正弦函数，例如：经过某一个坐标为(x0,y0)的像素点的所有直线的正弦函数可以表示为：
[0082]
r(θ)＝x0·
sinθ+y0·
cosθ
[0083]
其中，r是极坐标的极径，θ为极角，r(θ)函数是一个以θ为自变量的正弦函数；
[0084]
通过上述正弦函数的求导可将图像中所有的点在同一个坐标中绘制出不同的正弦函数图像。此外，直线检测的依据为：如果有多个像素点的正弦函数相交于一点，则就代表这些像素点它们通过同一条直线。因此，将图7闭运算处理后的图像中的像素点都用于绘制出各自的正弦函数图像，即可将正弦函数图像中，正弦函数相交的像素点作为同一直线上的像素点。通过上述直线检测时的图像如图8所示。
[0085]
需要补充的是，对于一张待处理钻孔岩芯图，经过霍夫直线检测后将得到长短不一，斜率不同的各类直线。为了满足切割需要，需根据长度、斜率等筛选出斜率为0的直线(即近似为是岩芯装载体的直线，岩心图像中直线的特征都是水平的)，具体做法是依次逐步取出控制斜率为0的固定范围直线，然后将这些直线在整幅图像检测完毕后，组成局部直线集合；此时得到待处理钻孔岩芯图中的多条水平直线。
[0086]
对于步骤s3、对待处理钻孔岩芯图中的水平直线数量进行检测，在直线数量不等于第一阈值时，则表明待处理钻孔岩芯图中的水平直线数量与所需分割线数量不相等，则需对待处理钻孔岩芯图的直线进行聚类，以得到聚类后的分割线；需要说明的是，上述第一阈值等于待处理钻孔岩芯图中纵向排列的管道数量加一。
[0087]
在一个优选的实施例中，所述对待处理钻孔岩芯图的直线进行聚类，包括：根据k-means算法对待处理钻孔岩芯图的直线进行聚类。
[0088]
在直线检测后，在待处理钻孔岩芯图中的各管道的边缘检测到的直线并非都是只有一条的，这是由于边缘直线受到泥土遮挡后，完整的直线会被分割成断断续续的多条直线，因此，对直线进行聚类就是将这些本属于同一边缘直线的分散线段聚类到一起，把它们归为同一条直线，以实现直线的修复。
[0089]
在本发明中采用k-means算法对待处理钻孔岩芯图的直线进行聚类，k-means算法过程简单而且主要应用于对于点的分类，本过程中只需要对于线段的纵坐标进行分类，所以使用k-means算法适用于本步骤；对于直线的聚类还可使用dbscan、optics等，但由于上述两种聚类方法相对复杂，在考虑处理速度的前提下，本发明采用k-means算法能提高图像
的处理速度。此外，k-means算法的自动聚类能忽略没有检测到的线条，能预测和弥补霍夫直线检测中检测不到的线条。通过k-means算法对直线进行聚类的方式为：先找到每类直线大概的位置，而后在直线的所有位置范围内找到当前直线集合中的直线元素的中值，最终确定直线的位置，作为分割图像的依据。具体的，k-means算法首先随机选取k个簇中心，然后定义代价函数：
[0090][0091]
其中，xi代表第i个样本,ci是xi所属于的簇，代表簇对应的中心点,m是样本总数；
[0092]
重复下面两个步骤直至代价函数j收敛：
[0093]
对于每一个样本xi，将其分配到距离最近的簇：
[0094][0095]
对于每一个类簇k，重新计算该类簇的中心：
[0096][0097]
需要说明的是，k-means算法在迭代时，若当前j没有达到最小值，那么首先固定簇中心μk，调整每个样例xi所属的类别ci，来让j函数减少，然后固定ci，调整簇中心μk使j减小。这两个过程交替循环，j单调递减，当j递减到最小值时，μk和ci也同时收敛。
[0098]
通过上述聚类后，可得到最终分割线。其中，该钻孔岩芯图中的分割线数量等于第一阈值。根据钻孔岩芯图中的分割线对钻孔岩芯图进行切割，并生成若干钻孔岩芯管道图。
[0099]
对于步骤s4、在直线检测后得到直线数量等于第一阈值时，表明此时图像中不存在分割线被泥土遮挡或缺失的情况，此时，之间基于检测到直线对待处理钻孔岩芯图进行分割，并生成若干钻孔岩芯管道图。
[0100]
对于步骤s5、对于步骤s3或步骤s4切割后得到的图像，将其按照管道中岩芯深度的次序依次进行拼接得到岩芯柱状图；需要补充的是，由于切割后的钻孔岩芯管道图的图像尺寸可能不一样，因此在进行图像依次拼接前，还需对图像进行缩放以使图像间的宽度和高度合适，继而生成更为美观的岩芯柱状图。此外，为使所生成的岩芯柱状图更便于查看，在生成岩芯柱状图时，对岩芯柱状图中每一管道所代表的岩芯深度进行标识。
[0101]
进一步地，可将上述方法实施在用于辅助岩土勘测信息软件上，用于满足用户对不同项目不同类型图片的保存。用户只需将拍摄的图像和文件上传到对应的媒体服务器，通过软件后台执行上述岩芯柱状图生成方法；同时，结合软件内部处理流程，可在生成岩芯柱状图时，同步生成每一管道对应的描述信息。
[0102]
通过实施本发明，有以下有益效果：
[0103]
1、能根据钻孔岩芯图自动生成岩芯柱状图，减少了认为处理时产生的误差，提高了岩芯柱状图生成的准确率和效率；
[0104]
2、在将钻孔岩芯图做直线检测前，采用契合本发明实施场景的图像处理方法对图
像做灰度化、滤波等处理，提高了直线检测的准确性，进一步提高了岩芯柱状图生成的准确率和效率。
[0105]
在上述方法项实施例的基础上，本发明对应提供了装置项实施例。
[0106]
如图9所示，本发明一实施例提供了一种岩芯柱状图生成装置，包括：数据获取模块、直线检测模块和直方图生成模块；
[0107]
所述数据获取模块，用于获取待处理钻孔岩芯图；
[0108]
所述直线检测模块，用于对待处理钻孔岩芯图做直线检测；
[0109]
所述直方图生成模块，用于在确定待处理钻孔岩芯图中的直线数量不等于第一阈值时，对待处理钻孔岩芯图的直线进行聚类，生成修复后的钻孔岩芯图，将修复后的钻孔岩芯图中所存在的直线作为分割线，根据修复后的钻孔岩芯图中的分割线对修复后的钻孔岩芯图进行切割，生成若干钻孔岩芯管道图；其中，所述第一阈值根据待处理钻孔岩芯图中管道的数量进行确定；在确定待处理钻孔岩芯图中的直线数量等于第一阈值时，根据待处理钻孔岩芯图的直线对待处理钻孔岩芯图进行分割，生成若干钻孔岩芯管道图；将各钻孔岩芯管道图进行拼接生成岩芯柱状图。
[0110]
需说明的是，以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。另外，本发明提供的装置实施例附图中，模块之间的连接关系表示它们之间具有通信连接，具体可以实现为一条或多条通信总线或信号线。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
[0111]
所述领域的技术人员可以清楚地了解到，为了方便和简洁，上述描述的装置的具体工作过程，可参考前述方法实施例中对应的过程，在此不再赘述。
[0112]
在上述方法项实施例的基础上，本发明对应提供了终端设备项实施例。
[0113]
本发明一实施例提供了一种终端设备，包括处理器、存储器以及存储在所述存储器中且被配置为由所述处理器执行的计算机程序，所述处理器执行所述计算机程序时实现本发明中任意一项所述的一种岩芯柱状图生成方法。
[0114]
终端设备可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述终端设备可包括，但不仅限于，处理器、存储器。
[0115]
所称处理器可以是中央处理单元(central processing unit，cpu)，还可以是其他通用处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、现成可编程门阵列(field-programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等，所述处理器是所述终端设备的控制中心，利用各种接口和线路连接整个终端设备的各个部分。
[0116]
所述存储器可用于存储所述计算机程序，所述处理器通过运行或执行存储在所述存储器内的计算机程序，以及调用存储在存储器内的数据，实现所述终端设备的各种功能。所述存储器可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少
一个功能所需的应用程序等；存储数据区可存储根据手机的使用所创建的数据等。此外，存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡(smart media card,smc)，安全数字(secure digital,sd)卡，闪存卡(flash card)、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
[0117]
在上述方法项实施例的基础上，本发明对应提供了存储介质项实施例。
[0118]
本发明一实施例提供了一种存储介质，所述存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述存储介质所在设备执行本发明中任意一项所述的一种岩芯柱状图生成方法。
[0119]
所述存储介质为计算机可读存储介质，所述计算机程序存储在所述计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、u盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括电载波信号和电信信号。
[0120]
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。

技术特征：
1.一种岩芯柱状图生成方法，其特征在于，包括：获取待处理钻孔岩芯图；对待处理钻孔岩芯图做直线检测；在确定待处理钻孔岩芯图中的直线数量不等于第一阈值时，对待处理钻孔岩芯图的直线进行聚类，生成修复后的钻孔岩芯图，将修复后的钻孔岩芯图中所存在的直线作为分割线，根据修复后的钻孔岩芯图中的分割线对修复后的钻孔岩芯图进行切割，生成若干钻孔岩芯管道图；其中，所述第一阈值根据待处理钻孔岩芯图中管道的数量进行确定；在确定待处理钻孔岩芯图中的直线数量等于第一阈值时，根据待处理钻孔岩芯图的直线对待处理钻孔岩芯图进行分割，生成若干钻孔岩芯管道图；将各钻孔岩芯管道图进行拼接生成岩芯柱状图。2.如权利要求1所述的一种岩芯柱状图生成方法，其特征在于，在对待处理钻孔岩芯图做直线检测之前，还包括：对待处理钻孔岩芯图做灰度化处理，生成第一钻孔岩芯图；对第一钻孔岩芯图做滤波处理，生成第二钻孔岩芯图；对第二钻孔岩芯图做直方图均衡化处理，生成第三钻孔岩芯图；对第三钻孔岩芯图做图像开运算和图像闭运算处理，生成第四钻孔岩芯图；所述对待处理钻孔岩芯图做直线检测，包括：对所述第四钻孔岩芯图做直线检测。3.如权利要求1任意一项所述的一种岩芯柱状图生成方法，其特征在于，所述对待处理钻孔岩芯图做直线检测，包括：获取待处理钻孔岩芯图中的所有像素点；对于任意一个像素点，计算待处理钻孔岩芯图中经过像素点的所有直线的正弦函数；根据所有像素点的正弦函数生成正弦函数图像；将正弦函数图像中，正弦函数相交的像素点作为同一直线上的像素点。4.如权利要求3所述的一种岩芯柱状图生成方法，其特征在于，所述对待处理钻孔岩芯图的直线进行聚类，包括：根据k-means算法对待处理钻孔岩芯图的直线进行聚类。5.如权利要求2所述的一种岩芯柱状图生成方法，其特征在于，对第一钻孔岩芯图做滤波处理，生成第二钻孔岩芯图，包括：对第一钻孔岩芯图做高斯滤波处理，生成第二钻孔岩芯图。6.一种岩芯柱状图生成装置，其特征在于，包括：数据获取模块、直线检测模块和直方图生成模块；所述数据获取模块，用于获取待处理钻孔岩芯图；所述直线检测模块，用于对待处理钻孔岩芯图做直线检测；所述直方图生成模块，用于在确定待处理钻孔岩芯图中的直线数量不等于第一阈值时，对待处理钻孔岩芯图的直线进行聚类，生成修复后的钻孔岩芯图，将修复后的钻孔岩芯图中所存在的直线作为分割线，根据修复后的钻孔岩芯图中的分割线对修复后的钻孔岩芯图进行切割，生成若干钻孔岩芯管道图；其中，所述第一阈值根据待处理钻孔岩芯图中管道的数量进行确定；在确定待处理钻孔岩芯图中的直线数量等于第一阈值时，根据待处理钻
孔岩芯图的直线对待处理钻孔岩芯图进行分割，生成若干钻孔岩芯管道图；将各钻孔岩芯管道图进行拼接生成岩芯柱状图。7.一种终端设备，其特征在于，包括处理器、存储器以及存储在所述存储器中且被配置为由所述处理器执行的计算机程序，所诉处理器执行所述计算机程序时实现如权利要求1至5中任意一项所述的一种岩芯柱状图生成方法。8.一种存储介质，其特征在于，所述存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述存储介质所在设备执行如权利要求1至5中任意一项所述的一种岩芯柱状图生成方法。

技术总结
本发明公开了一种岩芯柱状图生成方法、装置、终端设备及存储介质，该方法包括：获取待处理钻孔岩芯图；对待处理钻孔岩芯图做直线检测；在确定待处理钻孔岩芯图中的直线数量不等于第一阈值时，则对待处理钻孔岩芯图的直线进行聚类，生成修复后的钻孔岩芯图，将修复后的钻孔岩芯图中所存在的直线作为分割线，根据修复后的钻孔岩芯图中的分割线对修复后的钻孔岩芯图进行切割，生成若干钻孔岩芯管道图；其中，所述第一阈值根据待处理钻孔岩芯图中管道的数量进行确定；在确定待处理钻孔岩芯图中的直线数量等于第一阈值时，根据待处理钻孔岩芯图的直线对待处理钻孔岩芯图进行分割，生成若干钻孔岩芯管道图；将各钻孔岩芯管道图进行拼接生成岩芯柱状图。接生成岩芯柱状图。接生成岩芯柱状图。


技术研发人员：徐利成 冀盼彧 关汉辉 梁志鹏 于泽江 李儒挺 姚治远 关俊威 陆敏凤 曹进锋 林俊龙 刘轩 任沁枭 翁金塔
受保护的技术使用者：广州市交通设计研究院有限公司
技术研发日：2023.05.26
技术公布日：2023/8/24