[发明专利]一种CT图像环形伪影的迭代后处理去除方法

权利要求书

1.一种CT图像环形伪影的迭代后处理去除方法，其特征在于包括下述步骤：

(1)利用插值方法将带环形伪影的CT图像S转换为极坐标CT图像S_polar；

(2)利用图像结构提取算法从S_polar中提取结构图像S_structure，并从S_polar中减去S_structure生成纹理图像S_texture；

(3)从S_texture中提取条纹伪影图像S_stripe；

(4)利用插值方法将S_stripe转化为环形伪影图像S_ring,k，其中k为当前迭代的次数；

(5)从S中减去S_ring,k，得到校正的CT图像S_correct；

(6)判定迭代停止条件||S_ring,k||₁/||S_ring,1||₁≤Q₁是否成立，成立则环形伪影校正完成，不成立则用S_correct替换S并跳转至步骤(1)进行迭代计算，其中Q₁为迭代停止阈值，||S_ring,k||₁为S_ring,k的L₁范数，||S_ring,1||₁为第1次迭代得到的环形伪影图像S_ring,1的L₁范数。

2.根据权利要求1所述的一种CT图像环形伪影的迭代后处理去除方法，其特征在于：在所述步骤(3)中，从S_texture中提取S_stripe时需逐行进行相同操作，以从第i行纹理中提取第i行条纹伪影为例，具体操作步骤包括：

①计算的F(p,q)取最小值时所对应的p和q，计算公式为：

其中M为该行像素总数，C₁和C₂为该行的两类像素，C₁为该行第p个至第q个像素，p和q应满足M·T≤(q-p)≤M·(1-T)，T为阈值，C₂为该行除C₁外的其它所有像素，分别为C₁、C₂中的像素灰度均值，R_i为C₁中第i个像素的灰度值，R_j为C₂中第j个像素的灰度值；

②取C₁的条纹伪影强度为C₂的条纹伪影强度为由此可得该行条纹伪影其中n为当前迭代的次数；

③从中减去得到该行校正纹理

④判定迭代停止条件是否成立，成立则执行下一步，不成立则用替换并跳转至步骤①进行迭代计算，其中Q₂为迭代停止阈值，为的L₁范数，为第1次迭代得到的当前行条纹伪影的L₁范数；

⑤原始减去即得到该行条纹伪影