CBCT图像环形伪影去除方法技术

技术编号:15705295 阅读:145 留言:0更新日期:2017-06-26 12:44
本发明专利技术公开了一种CBCT图像环形伪影去除方法,步骤包括:对采集的原始数据进行预处理,然后重建出原始图像;对所述原始图像进行分割得到金属图像;对所述金属图像进行强度平滑性的低通滤波并归一化得到权值矩阵;对所述原始图像进行分频处理得到原始图像的高频部分;对所述原始图像进行归一化校正处理得到校正图像;对所述校正图像进行分频处理得到校正图像的高频部分和低频部分;将所述原始图像的高频部分、权值矩阵、以及校正图像的高频部分和低频部分进行加权融合得到最终校正图像。本发明专利技术不仅能够有效解决插值法会带来新伪影的问题,还能解决插值法会导致金属附近区域信息和骨骼信息丢失的问题。

Removal of annular artifacts in CBCT images

The invention discloses a method for removing CBCT image, ring artifact comprises the following steps: the original data pretreatment, and then reconstruct the original image; the original image is segmented by the metal image; the metal image intensity of the smoothness of the low-pass filter and the normalized weight matrix; frequency treatment of high frequency part of original image of the original image; the normalized correction correction image of the original image; the corrected image frequency to get high frequency and low frequency of the corrected image; the original image of the high-frequency part, the weight matrix and correction of image high frequency and low frequency part of the final corrected image fusion. The invention not only can effectively solve the problem that the interpolation method brings new artifacts, but also solves the problem that the interpolation method leads to the loss of the area information and the skeleton information near the metal.

【技术实现步骤摘要】
CBCT图像环形伪影去除方法
本专利技术涉及医学仪器成像领域,具体涉及一种CBCT图像环形伪影去除方法。
技术介绍
近年来,锥形束计算机断层成像系统(Cone-beamComputedTomography)在临床诊断、3D植入物等医疗领域得到了广泛的应用。作为断层成像的一种,锥形束CT具有成像精准、对射线利用率高等众多优点。但是,由于系统工艺、设计以及重建方法等多重因素的限制,锥形束CT图像中往往存在有一系列以重建中心为圆心且灰度值区别于周围像素的同心圆环,即环形伪影。环形伪影的出现严重的影响了图像的质量与真实性,干扰了进一步的临床诊断与处理。由此,如何在不影响图像分辨率的情况下对环形伪影进行去除成为了目前业界研究的重点。现有技术中,去除环形伪影的方法包括:(1)基于硬件的矫正方法,在该方法中,比较常见的为平场校正方法中的“两点校正法”。在这个过程中,将进行两次图像的采样。首先得到一张无样本的系统暗本底CBCT图像,在无偏差的情况下,这张平场图像应该是均匀的,由此这张图像实际代表了每个探测元之间的偏移,即不一致性。然后,再在均匀的光照条件下对被测物体进行一次成像。最后,将均匀光场图像与暗场图像做差,再用相对标定的方法对图像进行校正。(2)基于投影正弦图的校正方法,在早期的伪影去除领域,业界提出了基于传统傅里叶变换的滤波方式。其基本思想为在空间域上对被测物体的投影正弦图做离散傅里叶变换,由于在频域中,图像的原始信息多表现为低频信号,而伪影和部分细节则表现为高频信号,由此只要通过低频滤波的方式就可以去除图像中的高频信息。其后再将处理后图像变换回空间域,就达到了去除投影正弦图中条纹型伪影的目的。(3)基于重建图像的校正方法,为了节省存储空间,提高伪影去除的效率,使图像的矫正结果有着更加直观的体现。近年来,业界逐渐出现了多种针对于重建图像的后处理方法。与基于投影正弦图的前处理方法相似,为了改变环形伪影的几何特性,减少伪影检测与去除过程中的难度,需首先对重建后图像进行极坐标变换。根据此类变换的数学性质,笛卡尔坐标系中的环状伪影将在极坐标系中表现为条纹型伪影,经过后续的一系列处理以及逆坐标变换,就可以得到环形伪影校正后的重建CT图像。现有的环形伪影去除方法存在如下缺点:(1)基于硬件的矫正方法,由于方法本身的局限性,在现代的CT成像系统中,这个方法很难对环形伪影进行有效的去除,而另一种方法则比较特殊,其通过在CT成像过程中改变单个探测器的排列顺序以达到补偿探测器响应不一致的问题,虽然也起到了一定的效果,但对硬件加工工艺有着一定的要求,由此造价比较昂贵;(2)基于投影正弦图的校正方法,这种方法的缺点在于进行平滑低通滤波的同时也去除了图像的部分细信息,从而造成重建后图像的模糊与观察质量的明显下降。
技术实现思路
本专利技术的专利技术目的是提供一种CBCT图像环形伪影去除方法,不仅能够有效解决插值法会带来新伪影的问题,还能解决插值法会导致金属附近区域信息和骨骼信息丢失的问题。为达到上述专利技术目的,本专利技术采用的技术方案是:一种CBCT图像环形伪影去除方法,包括如下步骤:步骤一、通过全变分自适应滤波对采集的原始数据进行预处理,然后重建出含金属伪影的原始图像;步骤二、对所述原始图像进行分割得到金属图像;步骤三、对所述金属图像进行强度平滑性的低通滤波并归一化得到权值矩阵;步骤四、对所述原始图像进行分频处理得到原始图像的高频部分;步骤五、对所述原始图像进行归一化校正处理得到校正图像;步骤六、对所述校正图像进行分频处理得到校正图像的高频部分和低频部分;步骤七、将所述原始图像的高频部分、权值矩阵、以及校正图像的高频部分和低频部分进行加权融合得到最终校正图像。进一步地,所述步骤五中的对原始图像进行归一化校正处理包括如下步骤:(1)对所述原始图像进行均值聚类得到先验图像模型再投影;(2)对所述原始图像进行投影,通过利用原始图像投影除以先验图像投影得到归一化投影;(3)对所述金属图像再投影,由金属图像投影确定插值区域;(4)对归一化的投影进行线性插值后通过乘以先验图像投影做去归一化处理得到校正投影;(5)对校正投影进行重建得到校正图像。优选地,步骤(1)中所述的均值聚类为k均值聚类。由于上述技术方案运用,本专利技术与现有技术相比具有下列优点:本专利技术通过结合归一化方法解决了插值法会带来新伪影的问题,在重建图像之前通过自适应滤波滤除原始数据的噪声成分以减少重建后图像的条纹伪影,并加权了重建后图像的高频成分解决了插值法会导致金属附近区域信息和骨骼信息丢失的问题。附图说明图1是本专利技术实施例一的算法流程图。图2是本专利技术实施例二中去除环形伪影前的图像示意。图3是本专利技术实施例二中去除环形伪影后的图像示意。具体实施方式下面结合附图及实施例对本专利技术作进一步描述:实施例一:参见图1所示,一种CBCT图像环形伪影去除方法,包括如下步骤:步骤一、通过全变分自适应滤波对采集的原始数据进行预处理,然后重建出含金属伪影的原始图像;步骤二、对所述原始图像进行分割得到金属图像;步骤三、对所述金属图像进行强度平滑性的低通滤波并归一化得到权值矩阵;步骤四、对所述原始图像进行分频处理得到原始图像的高频部分;步骤五、对所述原始图像进行归一化校正处理得到校正图像;步骤六、对所述校正图像进行分频处理得到校正图像的高频部分和低频部分;步骤七、将所述原始图像的高频部分、权值矩阵、以及校正图像的高频部分和低频部分进行加权融合得到最终校正图像。本实施例中,所述步骤五中的对原始图像进行归一化校正处理包括如下步骤:(1)对所述原始图像进行均值聚类得到先验图像模型再投影;(2)对所述原始图像进行投影,通过利用原始图像投影除以先验图像投影得到归一化投影;(3)对所述金属图像再投影,由金属图像投影确定插值区域;(4)对归一化的投影进行线性插值后通过乘以先验图像投影做去归一化处理得到校正投影;(5)对校正投影进行重建得到校正图像。其中,步骤(1)中所述的均值聚类为k均值聚类。k均值聚类算法具体为:先随机选取k个对象作为初始的聚类中心。然后计算每个对象与各个种子聚类中心之间的距离,把每个对象分配给距离它最近的聚类中心。聚类中心以及分配给他们的对象就代表一个聚类。一旦全部对象都被分配了,每个聚类的聚类中心会根据聚类中现有的对象被重新计算。这个过程将不断重复直到满足某个终止条件。终止条件可以是以下任何一个:1、没有(或最小数目)对象被重新分配给不同的聚类;2、没有(或最小数目)聚类中心再发生变化;3、误差平方和局部最小。本专利技术基于分频去除CT金属伪影,通过结合归一化方法解决了插值法会带来新伪影的问题,在重建图像之前通过自适应滤波滤除原始数据的噪声成分以减少重建后图像的条纹伪影,并加权了重建后图像的高频成分解决了插值法会导致金属附近区域信息和骨骼信息丢失的问题。实施例二:通过实测数据验证本专利技术的环形伪影去除算法的有效性和可靠性,具体为,扫描参数为120kv,层厚为0.6mm,重建FOV为500mm,滤波反投影重建的图像大小为512×512,环形伪影去除前后对比参见图2和3所示,可以看出,图3中恢复了清晰的边缘信息和结构信息,本专利技术能够在完全去除金属伪影的同时保留金属区域更多的细节信息。因此,本专利技术对金属伪影的去除取得了本文档来自技高网
...
CBCT图像环形伪影去除方法

【技术保护点】
一种CBCT图像环形伪影去除方法,其特征在于,包括如下步骤:步骤一、通过全变分自适应滤波对采集的原始数据进行预处理,然后重建出含金属伪影的原始图像;步骤二、对所述原始图像进行分割得到金属图像;步骤三、对所述金属图像进行强度平滑性的低通滤波并归一化得到权值矩阵;步骤四、对所述原始图像进行分频处理得到原始图像的高频部分;步骤五、对所述原始图像进行归一化校正处理得到校正图像;步骤六、对所述校正图像进行分频处理得到校正图像的高频部分和低频部分;步骤七、将所述原始图像的高频部分、权值矩阵、以及校正图像的高频部分和低频部分进行加权融合得到最终校正图像。

【技术特征摘要】
1.一种CBCT图像环形伪影去除方法,其特征在于,包括如下步骤:步骤一、通过全变分自适应滤波对采集的原始数据进行预处理,然后重建出含金属伪影的原始图像;步骤二、对所述原始图像进行分割得到金属图像;步骤三、对所述金属图像进行强度平滑性的低通滤波并归一化得到权值矩阵;步骤四、对所述原始图像进行分频处理得到原始图像的高频部分;步骤五、对所述原始图像进行归一化校正处理得到校正图像;步骤六、对所述校正图像进行分频处理得到校正图像的高频部分和低频部分;步骤七、将所述原始图像的高频部分、权值矩阵、以及校正图像的高频部分和低频部分进行加权融合得到最终校正图...

【专利技术属性】
技术研发人员:王燕妮
申请(专利权)人:江苏美伦影像系统有限公司
类型:发明
国别省市:江苏,32

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1