CN108986035B - 一种含有两根岩芯的数字图像校正方法 - Google Patents
一种含有两根岩芯的数字图像校正方法 Download PDFInfo
- Publication number
- CN108986035B CN108986035B CN201710397114.8A CN201710397114A CN108986035B CN 108986035 B CN108986035 B CN 108986035B CN 201710397114 A CN201710397114 A CN 201710397114A CN 108986035 B CN108986035 B CN 108986035B
- Authority
- CN
- China
- Prior art keywords
- image
- correction
- core
- point
- axis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000003702 image correction Methods 0.000 title claims abstract description 19
- 239000011435 rock Substances 0.000 title claims abstract description 19
- 238000012937 correction Methods 0.000 claims abstract description 36
- 238000009499 grossing Methods 0.000 claims abstract description 8
- 238000005520 cutting process Methods 0.000 claims abstract description 6
- 150000001875 compounds Chemical class 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 3
- 239000011707 mineral Substances 0.000 abstract description 3
- 238000005286 illumination Methods 0.000 abstract description 2
- 238000003860 storage Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明属于矿产勘查技术领域,具体涉及一种含有两根岩芯的数字图像校正方法。本发明包括以下步骤:(1)从存储介质中读取采集的岩芯图像数据;(2)对图像进行裁剪,对有效区域建立直角坐标系;(3)对图像进行几何校正和亮度校正;(4)对图像进行平滑处理。本发明解决了现有方法每次只能对含有一根岩芯的数字图像进行校正、且没有考虑光源照射对于图像亮度的影响的技术问题。本发明可以对含有两根岩芯的数字图像进行校正,提高了工程效率,且考虑了照明光源对图像亮度的影响,使得校正后的图像亮度更加均匀,本发明的方法具有普遍适用性,可以简单类推出含有三根岩芯或更多根岩芯的数字图像校正方法。
Description
技术领域
本发明属于矿产勘查技术领域,具体涉及一种含有两根岩芯的数字图像校正方法。
背景技术
岩芯是研究地下岩石层结构、了解地下油气层和矿产情况最重要的实物地质资料。岩芯图像采集系统运用光学成像技术,把岩芯表面的光学影像转化为数字图像信息存储到计算机中。与保存岩芯实物相比,岩芯数字图像具有存取方便、不会随着时间推移而发生变化、能保证岩芯的原有特性等优点。
目前用于岩芯数字图像扫描的面阵采集头在使用中还存在一些问题。面阵采集头采集图像的理想位置是垂直于拍摄平面,但由于岩芯是圆柱状不能满足该要求。并且面阵采集头采用的是广角镜头,由于光学系统设计和制造的缺陷及安装误差,不可避免地会使采集得到的图像存在失真。失真图像不能反映岩芯真实特征信息,会影响地质专家对岩芯的分析及后期对岩芯图像的处理。因此,在运用此系统时,必须对失真图像进行校正,使其能够比较真实地反映岩芯的特征信息,满足后期对岩芯图像分析处理的要求。
现有的岩芯图像校正方法,每次只能对含有一根岩芯的数字图像进行校正,因此在获取岩芯数字图像过程中,需要对岩芯进行单根扫描,耗时较多;此外,现有的岩芯图像校正方法,没有考虑光源照射对于图像亮度的影响,校正后的图像亮度不够均匀。因此,亟需研究一种新的岩芯图像校正方法,能够对含有多根岩芯的数字图像进行校正,且充分考虑照明光源对图像亮度的影响,实现提高工程效率、使校正后图像亮度更加均匀的目的。
发明内容
本发明需要解决的技术问题为:现有的岩芯图像校正方法,每次只能对含有一根岩芯的数字图像进行校正,且没有考虑光源照射对于图像亮度的影响,存在工程效率低、校正后的图像亮度不均匀的缺陷。
本发明的技术方案如下所述:
一种含有两根岩芯的数字图像校正方法,包括以下步骤:步骤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):
对于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)为:
由此,可以求出图像各个新点的像素值,并且得到平滑的图像。
同理,对位于图像右侧的一根岩芯采用相同的方法进行几何校正、亮度校正和平滑处理。
以上的具体描述,对发明的目的、技术方案和有益效果进行了详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种含有两根岩芯的数字图像校正方法,其特征在于:包括以下步骤:步骤S1:读取图像;步骤S2:裁剪图像,建立坐标系;步骤S3:几何校正和亮度校正;步骤S4:平滑处理;
步骤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中,进行亮度校正的步骤如下所述:
通过下式进行亮度校正:
L'=L·μ
式中,
L为原图像某点的亮度;
μ为该点亮度校正因子;
L'为该点校正后的亮度。
2.根据权利要求1所述的一种含有两根岩芯的数字图像校正方法,其特征在于:步骤S2中,根据亮度差异裁剪有效区域,然后对含有两根岩芯的数字图像数据建立直角坐标系。
3.根据权利要求2所述的一种含有两根岩芯的数字图像校正方法,其特征在于:步骤S2中,所述直角坐标系建立方法如下:以数字图像上与两根岩芯上表面均相切的直线作为X轴方向,Y轴方向与X轴方向相垂直、且Y轴与两根岩芯横截面圆心的距离相等。
7.根据权利要求6所述的一种含有两根岩芯的数字图像校正方法,其特征在于:步骤S4中,根据原图像各点在校正后图像的位置以及它们的像素值,采用位置权重法计算校正后各新点的像素值,得到平滑的图像。
9.根据权利要求7所述的一种含有两根岩芯的数字图像校正方法,其特征在于:对位于图像右侧的一根岩芯进行几何校正、亮度校正和平滑处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710397114.8A CN108986035B (zh) | 2017-05-31 | 2017-05-31 | 一种含有两根岩芯的数字图像校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710397114.8A CN108986035B (zh) | 2017-05-31 | 2017-05-31 | 一种含有两根岩芯的数字图像校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108986035A CN108986035A (zh) | 2018-12-11 |
CN108986035B true CN108986035B (zh) | 2022-01-07 |
Family
ID=64501127
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710397114.8A Active CN108986035B (zh) | 2017-05-31 | 2017-05-31 | 一种含有两根岩芯的数字图像校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108986035B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104574420A (zh) * | 2015-01-29 | 2015-04-29 | 中国石油大学(华东) | 一种纳米级泥页岩数字岩心构建方法 |
CN106530335A (zh) * | 2015-08-06 | 2017-03-22 | 通用电气公司 | 在没有校准靶标的情况下的计算机断层扫描自校准 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8731248B2 (en) * | 2012-02-16 | 2014-05-20 | Arcsoft (Hangzhou) Multimedia Technology Co., Ltd. | Method of performing eye circle correction an image and related computing device |
-
2017
- 2017-05-31 CN CN201710397114.8A patent/CN108986035B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104574420A (zh) * | 2015-01-29 | 2015-04-29 | 中国石油大学(华东) | 一种纳米级泥页岩数字岩心构建方法 |
CN106530335A (zh) * | 2015-08-06 | 2017-03-22 | 通用电气公司 | 在没有校准靶标的情况下的计算机断层扫描自校准 |
Also Published As
Publication number | Publication date |
---|---|
CN108986035A (zh) | 2018-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110807809B (zh) | 基于点线特征和深度滤波器的轻量级单目视觉定位方法 | |
Bianco et al. | A comparative analysis between active and passive techniques for underwater 3D reconstruction of close-range objects | |
CN113112496B (zh) | 一种基于自适应阈值的亚像素轴类零件尺寸测量方法 | |
CN105783786A (zh) | 一种基于结构光视觉的零件倒角测量方法与装置 | |
CN107845073B (zh) | 一种基于深度图的局部自适应三维点云去噪方法 | |
EP3373239B1 (en) | Method and system for image correction using a quasiperiodic grid | |
CN109035170B (zh) | 基于单网格图分段映射的自适应广角图像校正方法及装置 | |
CN103500722B (zh) | 一种电子束缺陷检测方法 | |
CN111985507B (zh) | 一种岩体三维点云节理迹线提取方法 | |
CN111915517A (zh) | 一种适用于室内光照不利环境下rgb-d相机全局定位方法 | |
CN108445009A (zh) | 一种太阳能电池板裂纹检测方法 | |
CN111693028A (zh) | 一种基于投影影像获取数字水深模型的方法 | |
CN108986035B (zh) | 一种含有两根岩芯的数字图像校正方法 | |
CN107993193B (zh) | 基于光照均衡化和改进surf算法的隧道衬砌图像拼接方法 | |
CN108921936A (zh) | 一种基于光场模型的水下激光条纹匹配和立体重建方法 | |
CN105809657A (zh) | 一种角点检测方法和装置 | |
CN102339464A (zh) | 线搜索式角点检测方法 | |
CN115790539B (zh) | 一种合作靶标水下摄影测量方法 | |
CN114998397B (zh) | 一种多视角卫星影像立体像对优化选择方法 | |
Wu et al. | Defect inspection for underwater structures based on line-structured light and binocular vision | |
CN114396895B (zh) | 一种隧洞衬砌混凝土管片表面粗糙度的测量方法 | |
Ziqiang et al. | Research of the algorithm calculating the length of bridge crack based on stereo vision | |
CN111047635B (zh) | 基于深度图像的平面触摸方法、装置与触摸系统 | |
CN110211053B (zh) | 用于三维测量的快速精确相位匹配方法 | |
CN112907656A (zh) | 机器人的位置检测方法、检测装置、处理器和电子设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |