[发明专利]页岩CT图像重构方法及装置

摘要

本发明涉及图像处理技术领域，具体提供了一种页岩CT图像重构方法及装置，旨在解决如何在不损坏页岩样品的情况下，准确地获取其微观结构信息的技术问题。为此目的，本发明中的方法包括下述步骤：获取页岩样品的CT图像；依据所获取的CT图像，对预设的正则化模型进行迭代计算，得到所述页岩样品的重构图像。同时，本发明中的装置可以执行上述方法中的各步骤。本发明的技术方案，能够对页岩CT图像进行高分辨率重构，并消除重构伪影。

主权项

1.一种页岩CT图像重构方法，其特征在于，所述方法包括：获取页岩样品的CT图像d；依据所获取的CT图像d，采用预条件共轭梯度法对预设的正则化模型min(fα,β)进行迭代计算，得到所述页岩样品的重构图像m，具体为：对所述正则化模型min(fα,β)进行初始化；判断当前迭代目标函数fα,β对应的梯度值是否满足预设的迭代终止条件；若所述梯度值满足所述迭代终止条件，则停止迭代计算，并将当前正则化模型min(fα,β)中的m作为最终的重构图像；若所述梯度值不满足所述迭代终止条件，则执行下述操作：确定下一次迭代的迭代方向，具体为依据第k次迭代与第k‑1次迭代的目标函数的梯度值g(mk)与g(mk‑1)，并按照下式所示的方法计算所述第k‑1次迭代对应的正则化因子βk‑1：依据所计算的正则化因子βk‑1，并按照下式所示的方法计算第k次迭代对应的迭代方向hk＝‑g(mk)+βk‑1hk‑1；确定下一次迭代的迭代步长，并依据所确定的迭代步长以及迭代方向对所述正则化模型min(fα,β)进行数据更新；依据数据更新后的目标函数fα,β计算梯度值，并重新判断所计算的梯度值是否满足所述预设的迭代终止条件；其中，所述hk‑1为第k‑1次迭代对应的迭代方向；所述所述fα,β为需要最小化的目标函数，L表示离散形式的拉登变换，所述为Lm‑d的L2范数的平方，所述为m的L2范数的平方，所述||m||₁为m的L1范数，所述α与β均为预设的正则化因子。