[发明专利]岩心三维图像修复方法

说明书

技术领域

本发明涉及一种图像修复方法，尤其涉及一种岩心三维图像修复方法，属于三维图像修复技术领域。

背景技术

在石油地质分析过程中，通常需要获得三维岩心结构来定量地研究渗流的微观机理，分析储集层渗流性质和运移规律。随着CT（computed tomography，电子计算机X射线断层扫描技术）分辨率的提高，CT机已较多应用到岩石三维成像中，但CT的分辨率与样本尺寸是相矛盾的。为了得到准确的分析结果，需要获取高分辨率的孔隙结构图像，但扫描的样本尺寸受到很大限制，使岩样的代表性有所欠缺。因此很难得到可供分析的岩心高分辨率三维孔隙结构图像。针对这一问题，目前的主要解决方案是以二维薄片为基础进行三维重建，得到具有与薄片特性相近的三维图像。也有基于三维图像的三维重建，其原理与基于二维图像的三维重建相似。我们统一将此种基于局部特性重建出三维图像的方法称为三维数学建模。然而三维数学建模是基于局部特性，重建结果缺乏一定的实际代表性，同时存在随机性，且对于多相或非均质的岩心三维图像的重建存在较大的困难，重建效果与实际三维情况存在较大差距。

将大尺寸岩心样本切割为多个小样本，用CT机分别对小样本进行扫描，将获取的小样本图像拼接成较大的三维图像，能有效地解决这一问题。然而大尺寸样本在切割的过程中存在三维图像信息的缺失，且缺损信息较大，在小样本图像拼接成大样本图像的过程中需要恢复缺损的信息，实现大范围信息缺失的三维图像修复。因此本发明的目的就是提供一种可修复大范围信息缺失的岩心三维图像的方法，保持原有岩心三维图像的形态结构。

目前对于三维图像修复的研究还较少，大部分还是针对二维图像的修复。图像修复就是对图像的信息缺损区域进行信息填充的过程,其目的就是为了对有信息缺损的图像进行恢复,并且要使观察者无法察觉到图像曾经缺损或已被修复。目前的三维图像修复，主要是针对于医学MRI、CT图像小范围信息缺失的三维图像修复以及视频的修复。Satoru Morita（Three Dimensional Image Inpainting[C].In proceeding of:Advances in Natural Computation,Second International Conference,Xi'an,China,2006,4222:752-761）将基于偏微分方程PDE的方法拓展为三维修复算法，并将该算法应用于医学图像和视频的修复，但是该方法只是针对小范围信息缺失的三维图像修复，且主要是对于结构的修复较好，而对于纹理部分却存在模糊现象。由于视频中帧内的相关度往往大于帧间的相关度，帧内信息较为重要，所以对于视频的修复主要是逐帧进行，较少考虑帧之间的关系。而对于岩心三维图像，其在各个方向上的相关性都一样，需同等考虑各个方向的信息。

发明内容

本发明的目的就在于为了解决上述问题而提供一种复原度极高的缺失信息的岩心三维图像修复方法。

本发明通过以下技术方案来实现上述目的：

1、本发明所述岩心三维图像修复方法包括以下步骤：

（1）提取缺失信息的岩心三维图像即原始三维图像中待修复区域Ω的初始边界δΩ⁰，并设定模板大小；

（2）对原始三维图像的已知区域Φ进行区域划分；

（3）计算原始三维图像的边界上所有体素q的边界模板的优先级值，所述体素q指三维数据场中的采样点；

（4）搜索最大优先级值的边界模板Ψ_p，并保存边界体素P的坐标；

（5）判断边界体素P的区域类型，并设置边界模板Ψ_p的匹配块的全局搜索区域；

（6）根据边界体素P的匹配中心因子MCF和最大距离MD在相应的范围内搜索最优匹配块Ψ_q，并进行相应体素的灰度值填充；

（7）更新或设置以下信息：置信度、边界面δΩ、新产生的边界体素P′的匹配中心因子MCF和最大距离MD，以及新产生的边界模板Ψ_p′的优先级值；

（8）重复步骤（4）～（7）直至填充完待修复区域Ω，即

上述方法的基本原理如下：