一种含有两根岩芯的数字图像校正方法与流程

本发明属于矿产勘查技术领域，具体涉及一种含有两根岩芯的数字图像校正方法。

背景技术

岩芯是研究地下岩石层结构、了解地下油气层和矿产情况最重要的实物地质资料。岩芯图像采集系统运用光学成像技术，把岩芯表面的光学影像转化为数字图像信息存储到计算机中。与保存岩芯实物相比，岩芯数字图像具有存取方便、不会随着时间推移而发生变化、能保证岩芯的原有特性等优点。

目前用于岩芯数字图像扫描的面阵采集头在使用中还存在一些问题。面阵采集头采集图像的理想位置是垂直于拍摄平面，但由于岩芯是圆柱状不能满足该要求。并且面阵采集头采用的是广角镜头，由于光学系统设计和制造的缺陷及安装误差，不可避免地会使采集得到的图像存在失真。失真图像不能反映岩芯真实特征信息，会影响地质专家对岩芯的分析及后期对岩芯图像的处理。因此，在运用此系统时，必须对失真图像进行校正，使其能够比较真实地反映岩芯的特征信息，满足后期对岩芯图像分析处理的要求。

现有的岩芯图像校正方法，每次只能对含有一根岩芯的数字图像进行校正，因此在获取岩芯数字图像过程中，需要对岩芯进行单根扫描，耗时较多；此外，现有的岩芯图像校正方法，没有考虑光源照射对于图像亮度的影响，校正后的图像亮度不够均匀。因此，亟需研究一种新的岩芯图像校正方法，能够对含有多根岩芯的数字图像进行校正，且充分考虑照明光源对图像亮度的影响，实现提高工程效率、使校正后图像亮度更加均匀的目的。

技术实现要素：

本发明需要解决的技术问题为：现有的岩芯图像校正方法，每次只能对含有一根岩芯的数字图像进行校正，且没有考虑光源照射对于图像亮度的影响，存在工程效率低、校正后的图像亮度不均匀的缺陷。

本发明的技术方案如下所述：

一种含有两根岩芯的数字图像校正方法，包括以下步骤：步骤s1：读取图像；步骤s2：裁剪图像，建立坐标系；步骤s3：几何校正和亮度校正；步骤s4：平滑处理。

作为优选方案：步骤s2中，根据亮度差异裁剪有效区域，然后对含有两根岩芯的数字图像数据建立直角坐标系。步骤s2中，所述直角坐标系建立方法如下：以数字图像上与两根岩芯上表面均相切的直线作为x轴方向，y轴方向与x轴方向相垂直、且y轴与两根岩芯横截面圆心的距离相等。

作为优选方案：步骤s3中，进行几何校正的步骤如下所述：

光源点p为y轴上某点，表示探测器和光源的位置；o为xy坐标系原点；o1’为位于图像左侧的一根岩芯在xy平面横截面的圆心，c为o1’在x轴的投影；o1为直线po1’与该岩芯横截面外表面的交点，b为o1在x轴的投影；a为直线po1’与x轴的交点，即a为o1所成的像，根据几何投影关系，此点没有失真；x1为该岩芯横截面外表面一点；x1’为直线px1与x轴的交点，即x1’为x1所成的像；x1”为x1’经校正后的在x轴的新位置，线段ax1”的长度与线段o1x1对应的弧长l1相等；β11为∠x1’pa的角度；β12为∠apo的角度；α1为∠x1o1’o1的角度；a1'为线段oa的长度；a1为线段ox1’的长度；h为线段op的长度；r1为该岩芯半径；d1为y轴距离该岩芯横截面外表面的最短距离；l1为线段o1x1对应的该岩芯横截面外表面上的弧长；

对于x1’经校正后的在x轴的新位置x1”，而线段ax1”的长度与线段o1x1对应的该岩芯横截面外表面上的弧长相等，即|ax1”|＝l1，采用下式计算l1：

式中，

(x1,y1)为o1点坐标；

(x,y)为x1点坐标；

采用以上步骤对采集的岩芯横截面外表面上各点进行计算，确定其校正后的新位置，完成几何校正。

步骤s3中，通过下式计算o1点坐标(x1,y1)：

步骤s3中，通过下式计算x1点坐标(x,y)：

步骤s3中，进行亮度校正的步骤如下所述：

通过下式进行亮度校正：

l'＝l·μ

式中，

l为原图像某点的亮度；

μ为该点亮度校正因子；

l'为该点校正后的亮度。

步骤s3中，

作为优选方案：步骤s4中，根据原图像各点在校正后图像的位置以及它们的像素值，采用位置权重法计算校正后各新点的像素值，得到平滑的图像。步骤s4中，所述位置权重法的具体计算方法如下：

假设原始图像同一行相邻两个点i和i+1，它们的像素值分别为img(i)和img(i+1)，校正后它们在该行的位置分别变为p和q，则该行它们之间的新点n，q≥n≥p，像素值img(n)为：

作为优选方案：对位于图像右侧的一根岩芯采用相同的方法进行几何校正、亮度校正和平滑处理。

本发明的有益效果为：

(1)本发明的一种含有两根岩芯的数字图像校正方法，可以对含有两根岩芯的数字图像进行校正，在获取岩芯数字图像过程中可以同时扫描两根岩芯进行数字图像数据采集，提高了工程效率，节省了时间，在工程上具有很高的实用价值；

(2)本发明的一种含有两根岩芯的数字图像校正方法，可以推广到同时采集三根或者更多根岩芯，进一步提高工程效率、节省扫描时间；

(3)本发明的一种含有两根岩芯的数字图像校正方法，考虑了照明光源对图像亮度的影响，使得校正后的图像亮度更加均匀；

(4)本发明的一种含有两根岩芯的数字图像校正方法，具有普遍适用性，可以类推出含有三根岩芯或更多根岩芯的数字图像校正方法。

附图说明

图1为本发明的一种含有两根岩芯的数字图像校正方法流程图；

图2为本发明的一种含有两根岩芯的数字图像校正方法坐标示意图。

具体实施方式

下面结合附图和实施例对本发明的一种含有两根岩芯的数字图像校正方法进行详细说明。

实施例1

如图1所示，本实施例的一种含有两根岩芯的数字图像校正方法，包括如下步骤：步骤s1：读取图像；步骤s2：裁剪图像，建立坐标系；步骤s3：几何校正和亮度校正；步骤s4：平滑处理。

步骤s1：读取图像

从存储介质中读取含有两根岩芯的数字图像数据。

步骤s2：裁剪图像，建立坐标系

根据亮度差异裁剪有效区域，然后对含有两根岩芯的数字图像数据建立直角坐标系。所述直角坐标系的建立方法如下所述：如图2所示，以数字图像上与两根岩芯上表面均相切的直线作为x轴方向，y轴方向与x轴方向相垂直、且y轴与两根岩芯横截面圆心的距离相等(即两根岩芯相对于y轴左右对称)，xy平面与岩芯横截面平行。

步骤s3：几何校正和亮度校正

选取合适的步长进行采点，进行几何校正和亮度校正。

如图2中所示，光源点p为y轴上某点，表示探测器和光源的位置；o为xy坐标系原点；o1’为位于图像左侧的一根岩芯在xy平面横截面的圆心，c为o1’在x轴的投影；o1为直线po1’与该岩芯横截面外表面的交点，b为o1在x轴的投影；a为直线po1’与x轴的交点，即a为o1所成的像，根据几何投影关系，此点没有失真；x1为该岩芯横截面外表面一点；x1’为直线px1与x轴的交点，即x1’为x1所成的像；x1”为x1’经校正后的在x轴的新位置，线段ax1”的长度与线段o1x1对应的弧长l1相等。β11为∠x1’pa的角度；β12为∠apo的角度；α1为∠x1o1’o1的角度；a1'为线段oa的长度；a1为线段ox1’的长度；h为线段op的长度；r1为该岩芯半径；d1为y轴距离该岩芯横截面外表面的最短距离；l1为线段o1x1对应的该岩芯横截面外表面上的弧长。

由于△oap∽△o1’ac∽△o1ab，通过下式计算o1点坐标(x1,y1)：

由于直线px1的方程为x1所在的圆的方程为(x+r1+d1)²+(y+r1)²＝r1²，通过下式计算x1点坐标(x,y)：

对于x1’经校正后的在x轴的新位置x1”，而线段ax1”的长度与线段o1x1对应的该岩芯横截面外表面上的弧长相等，即|ax1”|＝l1，采用下式计算l1：

采用以上步骤对采集该岩芯横截面外表面上各点进行计算，确定其校正后的新位置，完成几何校正。

通过下式进行亮度校正：

l'＝l·μ

其中，

式中，

l为原图像某点的亮度；

μ为该点亮度校正因子；

l'为该点校正后的亮度。

步骤s4：平滑处理

使用相同步长确定校正后图像的新点，根据原图像各点在校正后图像的位置以及它们的像素值，采用位置权重法计算各新点的像素值，得到平滑的图像。

所述位置权重法的具体计算方法如下：

假设原始图像同一行相邻两个点i和i+1，它们的像素值分别为img(i)和img(i+1)，校正后它们在该行的位置分别变为p和q，则该行它们之间的新点n，q≥n≥p，像素值img(n)为：

由此，可以求出图像各个新点的像素值，并且得到平滑的图像。

同理，对位于图像右侧的一根岩芯采用相同的方法进行几何校正、亮度校正和平滑处理。

以上的具体描述，对发明的目的、技术方案和有益效果进行了详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。