一种岩石CT图像目标分割方法技术

本发明专利技术公开了一种岩石CT图像分割方法，主要分为5个步骤：(1)对岩石CT图像进行自适应阈值分割，所谓自适应阈值分割是根据CT图像中相邻图像相关性，在分割过程中自动调节分割阈值；(2)对分割结果进行视觉判断，判断是否满意，是否需要进一步修复；(3)如果需要修复，则对分割结果自动判别分割优劣，进行自动区域生长修复；(4)对自动修复结果再此进行人眼判断，判断是否有个别感兴趣的目标需要进一步手工修复；(5)如果需要进一步手工修复，则进行手工选取特定目标的自动修复。通过本发明专利技术分割岩石CT图像的效率很高，整个过程中几乎没有人为参与，对大规模的CT图像分割效果很好。

【技术实现步骤摘要】
一种岩石CT图像目标分割方法
本专利技术设计一种CT图像分割方法，尤其涉及一种岩石CT图像目标分割方法，属于CT图像分割
技术背景在石油地质分析过程中，CT成像技术(ComputedTomography，电子计算机X射线断层扫描技术)得到越来越广泛的应用，为了分析演示孔隙和颗粒结构，得到连通性、渗透率和形状因子等微观信息，通常采用对岩石进行CT扫描得到二维序列图，再根据二维序列图进行相应目标(孔隙和颗粒)分割提取，然后重建岩石的三维立体模型的方法。二维图像分割结果的好坏，直接关系到重建之后三维模型的真实程度。图像分割是目标识别和特征提取的基础，是近年来计算机视觉研究的热点之一。图像分割技术发展到现在，产生了很多从经典到前沿的分割方法。主要分为以下几种：阈值分割，基于区域的分割，基于边缘的分割，基于特定数学理论的分割(包括数学形态学，模糊理论，神经网络等)。其中经典的分割技术当属阈值分割。根据特定的阈值来判断目标像素，从而得到感兴趣的目标。而根据阈值的选取的个数不同，阈值分割又分为单阈值，双阈值，多阈值等分割方法。双阈值分割方法简单有效，速度较快。在岩石图像中，感兴趣的目标往往是孔隙和颗粒，其边缘比较模糊，需要通过预览手动调节合适的阈值来进行目标分割。但对于作为三维重建的岩石CT图像，图像的数量巨大，手动对每一张图像设定合适阈值进行分割的工作量非常大。针对这一问题，目前的一种解决方案是选择其中一幅图像，设定合适的阈值，其他所有图像都根据这个阈值进行分割，或者将CT图像分几个小组，每个小组选择一个合适的阈值进行分割。另一种解决方案是针对每一幅图像使用合适的分割方法，每一幅图像独立分割。第一种解决方案中分组分割是在分割效果和效率中间取一个平衡点。这种方法一方面平衡点的寻找需要多次的试探，工作量本身就很大，另一方面如果CT图像亮度变化比较剧烈，这种解决方案可能需要对每一幅图进行单独分割，工作量不会得到有效的降低。第二种解决方案理论上是最好的方法，整个过程几乎不需要进行人为干预。但是分割技术发展到现在，还没有一种分割方法对亮度不均匀的所有CT图像都能根据人的视觉判断进行理想的分割。况且如果分割方法过度复杂，对于大规模的CT图像进行处理的效率将会很低。而且在图像分割过程中，经常会有对分割结果进行局部手工修复的需求。并且期望手工修复的时候可以尽量降低工作量。因为当规模较大的时候，对每一幅图像进行手工修复的工作量是巨大的。然而对所有图像采用同一种固定修复方法，又很难达到预期的修复效果。针对这种需求，目前的研究情况还没有提出一种可以对岩石CT图像进行自适应，高效的分割和修复方法。本专利技术针对这种问题和需求提出并实现了一种可以根据人眼视觉判断，进行的自适应，并且高效的分割方法，并且利用序列图的相关性提出并实现了一种自适应的修复方法。可以在不加大工作量的前提下，进而对分割效果进一步提高。
技术实现思路
本专利技术的目的就是为了满足上述需求而提出并实现一种具有实用性的，高效的岩石CT图像目标分割方法。本专利技术提供的岩石CT图像目标分割方法，是基于以下思想提出的：岩石CT图像具有规模大，灰度对比度低、纹理单一、噪声较多、且结构不规则等特征。由于CT扫描技术的原因，通常会还会造成CT图像前后亮度不均匀，亮度差别较大。岩石内部结构比较复杂，直接进行分割难度较大。但是CT图像有一个很重要的特征，就是相邻CT图的特征信息基本保持一致，包括亮度，结构，统计特性等。而且同一地质环境下的CT图像，孔隙和颗粒在岩石中的分布有很大相关性，特征变化缓慢。合理利用这些特性可以解决很多目前比较难以解决的问题。而且大规模图像处理过程中，图像处理的效果和效率，以及人所参与的工作量，是判断整个图像处理方法优劣的重要指标。均衡这些指标也正是图像处理方法研究的关键所在。图像分割作为图像处理中的代表应用技术不可能避开这些指标。在效果方面，优秀的图像分割技术要保证能将所有感兴趣的目标分割出来，并且非感兴趣目标要被抑制。由于CT图像虽然具有很大相关性，但是图像之间亮度有一定差别，采用固定阈值分割的时候，无法保证每一幅CT图像都达到预期的分割效果。而根据相关性信息，自动调节分割阈值，可以达到预期的分割效果。但即便阈值自动调节，分割效果得到很大改善，还会难免有少量图像达不到人眼视觉判断的分割效果，存在一定的欠分割现象。所谓欠分割就是没有分割彻底，较大的目标只被分割出少量区域。以当前的分割结果，自动进行区域生长，将目标区域进一步修复完善，便可以达到较理想的分割效果。在效率方面，工作量(包括人和计算机资源)要在可接受范围内尽可能降低。自动阈值调节和自适应区域生长就是在改善分割效果的同时兼顾效率的方法。即便在最终对个别目标进行人工修复的过程中，针对相关目标自适应进行区域生长修复也是将效率作为重要指标。只需对一幅图像中的目标进行手工修复，其他CT图中的目标将自动进行判别修复。综合要求分割方法简单有效，分割方法不能太复杂，需尽量使整个过程中没有人参与或者人参与成份很少。本专利技术提供岩石CT图像目标分割方法，包括以下步骤：步骤1：对岩石CT图像进行自动调节阈值分割，把感兴趣的目标提取出来；步骤2：对步骤1目标提取的结果进行人眼视觉判断，分割效果是否达到要求，如果未达到要求，需要进一步对分割结果修复，进行步骤3，否则分割结束；步骤3：对步骤1分割结果进行自动区域生长修复，使提取的目标更加完善；步骤4：对步骤3目标修复的结果进行人眼视觉判断，目标修复效果是否达到要求，如果没有达到要求，需要进一步对分割修复的结果进行手工修复，进行步骤5，否则分割结束；步骤5：对步骤3分割修复结果进行手工选取特定目标自动修复，以修复的结果作为最终提取的目标。在本专利技术上述技术方案中，所述步骤1对岩石CT图像进行自动调节阈值分割，优先采用由下述方法进行分割：(1)选定CT图像中的第一幅图像，通过双阈值二值化方法，手动调阈值，根据人眼视觉对感兴趣目标区域的判断，选定合适的阈值(thL，thU),其中灰度值在[thl,thU]范围内的像素点属于分割区域，使当前图像的分割效果较为理想。(2)计算图像的直方图H0(x)，根据H0(x)记录用于分割的下阈值点(thL，H0(thL))和与其相邻的左右各5个点组成下阈值匹配模版Tl。同样记录用于分割的上阈值点(thU，H0(thU)和与其相邻的左右各5个点组成上阈值匹配模版Tu，Tl和Tu的中心点分别为(thL，H0(thL))和(thU，H0(thU))。(3)计算当前图像的直方图H(x),使用Tl和Tu在H(x)上进行匹配。计算Tl在H(x)上对应位置差值绝对值τ，找到使τ最小的中心点thli，将thli用于当前图像分割的下阈值点。使用同样的方法找到Tu最匹配的中心点thui，将thui用于当前图像分割的上阈值点。(4)对当前图像使用(thli，thui)阈值进行分割，并用thli更新thL，用thui更新thU，用H(x)更新H0(x)。(5)重复步骤(2)，(3)，(4)，直到所有图像分割完毕。在本专利技术上述技术方案中，所述步骤3对分割后的岩石CT图像进行自动区域生长修复，优先采用由下述方法进行修复：(1)根据人眼视觉判断，选择一幅分割比较理想的图片(通常直本文档来自技高网...

【技术保护点】
一种岩石CT图像目标分割方法，其特征在于：包括以下步骤：步骤1：对岩石CT图像进行自动调节阈值分割，把感兴趣的目标提取出来；步骤2：对步骤1目标提取的结果通过人眼视觉进行判断，判断分割效果是否达到要求，如果未达到要求，需要进一步对分割结果修复，进行步骤3，否则分割结束；步骤3：对步骤1分割结果进行自动区域生长修复，使提取的目标更加完善；步骤4：对步骤3目标修复的结果通过人眼视觉进行判断，判断目标修复效果是否达到要求，如果没有达到要求，需要进一步对分割修复的结果进行手工修复，进行步骤5，否则分割结束；步骤5：对步骤3分割修复结果进行手工选取特定目标自动修复，以修复的结果作为最终提取的目标。

【技术特征摘要】
1.一种岩石CT图像目标分割方法，包括以下步骤：步骤1：对岩石CT图像进行自动调节阈值分割，把感兴趣的目标提取出来；步骤2：对步骤1目标提取的结果通过人眼视觉进行判断，判断分割效果是否达到要求，如果未达到要求，需要进一步对分割结果修复，进行步骤3，否则分割结束；步骤3：对步骤1分割结果进行自动区域生长修复，使提取的目标更加完善；步骤4：对步骤3目标修复的结果通过人眼视觉进行判断，判断目标修复效果是否达到要求，如果没有达到要求，需要进一步对分割修复的结果进行手工修复，进行步骤5，否则分割结束；步骤5：对步骤3分割修复结果进行手工选取特定目标自动修复，以修复的结果作为最终提取的目标；步骤1所述的对岩石CT图像进行自动调节阈值分割，采取下述方法进行：(1)选定CT图像中的一幅图像，作为分割开始的基准图像，通过双阈值二值化分割方法，根据人眼视觉对感兴趣目标区域的判断，手动选定阈值区间[thL，thU],其中灰度值在[thL，thU]范围内的像素点属于分割区域；(2)计算基准图像的直方图H0(x)，根据H0(x)记录用于分割的下阈值点(thL，H0(thL))和与其相邻的左右各5个点组成下阈值匹配模版Tl，同样记录用于分割的上阈值点(thU，H0(thU)和与其相邻的左右各5个点组成上阈值匹配模版Tu，Tl和Tu的中心点分别为(thL，H0(thL))和(thU，H0(thU))；(3)计算当前图像的直方图H(x),使用Tl和Tu在H(x)上进行匹配，计算Tl在H(x)上对应位置差值绝对值和τ，找到使τ最小的中心点thli，将thli用于当前图像分割的下阈值点，使用同样的方法找到Tu最匹配的中心点thui，将thui用于当前图像分割的上阈值点；(4)对当前图像使用(thli，thui)阈值进行分割，并用thli更新thL，用thui更新thU，用H(x)更新H0(x)；(5)重复步骤(2)，(3)，(4)，直到所有图像分割完毕。2.根据权利要求1所述的岩石CT图像目标分割方法，其特征在于，步骤3对自动调节阈值分割结果进行自动区域生长修复，采取下述方法进行：(1)根据人眼视觉判断，选择一幅分割较理想的图像作为基准图像，计算其分割区域的覆盖率δ0，以δ0作为分割优劣的标准；(2)计算CT图像中每一幅图已分割的覆盖率δ，以及δ和δ0的比值ω，根据ω将当前图像的分割程度归类到A1,A2,A3,A4；(3)对CT图像进行扫描，通过区域生长方法对已分割区域进行识别，得到已分割目标区域的像素点集Ω，并计算Ω的平均灰度值Ψ；(4)计算判别像素点P的8邻域灰度均值和Ψ的差值Δ以及8邻域灰度值与判别像素点P灰度的平均差值ε，将Δ和ε作为区域生长的判别准则,并根据A1,A2,A3,A4不同类别自动设定特定的Δ和ε判别值，对目标区域进行区域生长；(5)在Ω的基础上根据步骤(4)中的判别规则对目标区域进行改进的区域生长方法修复，使目标区域分割更加完整；(6)重复步骤(4)和(5)直到图像中所有目标都修复完毕。3.根据权利要求1所述的岩石CT图像目标分割方法，其特征在于，步骤5对分割结果进行手工选取特定目标自动修复，采取下述方法进行：(1)对CT图像中分割不理想的目标进行手工修复，得到修复区域的点集Φ；(2)将Φ所有像素点作为种子点，计算Φ区域的加权平均灰度值φ和计算判别像素点P的8邻域灰度均值与φ的差值Δ，并计算P的8邻域灰度值λi与判别像素点P的灰度值λ4的平均差值ε，这两个值在一定范围内的像素点将归并到当前区域中，最终得到区域A；(3)计算区域A的面积s0和像素加权平均值g0，并运用拉普拉斯变换在图形层上提取A的边缘D；(4)将区域A投影到下一幅相邻图像，通过对应位置图像灰度值的比较，进行像素点取舍，得到投影点集B；以点集B代替Φ作为种子点，重复步骤(2)得到生长完成的区域C；(5)将当前区域生长完成的区域C与区域A做比较，将C中存在而在A中不存在的像素点集标记为区域F；(6)计算边缘D的总线长d0，以及区域F与边缘D的交线线长d...

【专利技术属性】
技术研发人员：滕奇志，徐永进，何小海，袁浩，吴小强，王正勇，
申请(专利权)人：四川大学，
类型：发明
国别省市：四川;51