[发明专利]基于CT图像的岩心目标提取方法

权利要求书

1.一种基于CT图像的岩心目标提取方法，其特征在于：包括以下步骤：

S1.读取岩心的CT图像，并将所得CT图像记为I₀：f(x,y)=g，其中，(x,y)为图像像素坐标，g为图像像素坐标(x,y)对应的灰度值，g∈[0,1]；

S2.岩心CT图像的阈值变换：具体包括以下步骤：

S21.设灰度值变换的阈值为0；

S22.遍历整张岩心CT图像的像素，如果像素的灰度值为0，则保持其灰度值不变，如果像素的灰度值大于0，则将其灰度值设置为1，相应的表达式为f1(x,y)=0g=01g>0,]]>所得二值图像的结果记为I₁；

S3.对图像I₁进行轮廓提取；

S4.对S3所得的每个轮廓进行填充，即将每个轮廓内部的像素的灰度值都设为1，得到图像为I₂：f₂(x,y)；

S5.对S4所得每个连通区域进行像素数目统计，记录每个连通区域包含的像素数目，并对统计结果进行比对，获得具有最大像素数目的连通区域；

S6.将除具有最大像素数目的连通区域外的连通区域的像素的灰度值都设置为零，得到图像I₃：f₃(x,y)；

S7.对图像I₃进行轮廓提取，获得目标轮廓并标记为I₄：f₄(x,y)；

S8.将图像I₃的每个像素的灰度值与原始图像I₀中相应位置的像素的灰度值相乘，即得基于CT图像的岩心目标图像，记为I₅：f₅(x,y)＝f₀(x,y)×f₃(x,y)；

S9.将目标图像I₅与轮廓图像I₄进行加运算，即得基于CT图像的岩心目标及岩心目标轮廓。

2.根据权利要求1所述的基于CT图像的岩心目标提取方法，其特征在于：S3所述的轮廓提取具体包括以下步骤：

S31.遍历I₁的每个像素，对于灰度值为1的像素，如果该像素的八个邻域点的灰度值都为1，则将该像素的灰度值设置为零，即可提取I₁中所有轮廓；

S32.按照从左到右、从下到上的顺依次对S31所得轮廓进行八邻域点追踪，获得单值化的轮廓。

3.根据权利要求2所述的基于CT图像的岩心目标提取方法，其特征在于：所述八邻域点追踪具体如下：找到第一个灰度值为1的像素，假设为最左下方的边界点，将该点记为第一边界点，该像素的右、右上、上、左上四个邻域中至少有一个是边界点，记为第二边界点，从第二边界点开始，按右、右上、上、左上、左、左下、下和右下的顺序找相邻点中的边界点，记为第三边界点，如果第三边界点就是第一边界点，则表明闭合轮廓追踪完成；否则，从第三边界点继续，直到找到第一边界点为止。