CT图像伪影去除方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN112785520A(43)申请公布日2021.05.11(21)申请号202110059221.6(22)申请日2021.01.14(71)申请人明峰医疗系统股份有限公司地址310016浙江省杭州市经济技术开发区益丰路129号5幢1-2层(72)发明人徐亦飞陈伟潘华素(74)专利代理机构绍兴市越兴专利事务所(普通合伙)33220代理人蒋卫东(51)Int.Cl.G06T5/00(2006.01)G06T5/50(2006.01)权利要求书1页说明书3页附图6页(54)发明名称CT图像伪影去除方法(57)摘要本发明提供了一种CT图像伪影去除方法，涉及医学影像技术领域，通过对输入图像进行滤波得到带有边缘保护区域的图像纹理，获得图像纹理在频率阈的相位分布和幅值分布，计算幅值分布中相同频率成分的均值并做归一化，得到归一化的幅值分布，并根据设定的阈值和低频掩膜得到峰值的分布；对峰值分布进行高斯滤波得到滤波矩阵，得到无混叠的纹理频域图像；结合图像纹理在频率阈的相位分布，对无混叠的纹理频域图像处理，得到无混叠的纹理空间域图像；然后将图像纹理与无混叠的纹理空间域图像相减，得到混叠空间域图像；最后将输入图像减去混叠伪影图像，得到消除混叠伪影之后的图像。本发明方法在去除混叠伪影的同时，可保留更多的图像细节。CN112785520ACN112785520A权利要求书1/1页1.一种CT图像伪影去除方法，其特征在于包括如下步骤：（1）输入图像，对输入图像进行滤波得到带有边缘保护区域的图像纹理；（2）获得图像纹理在频率阈的相位分布和幅值分布；（3）计算幅值分布中相同频率成分的均值并做归一化，得到归一化的幅值分布，并根据设定的阈值和低频掩膜得到峰值的分布；（4）对峰值分布进行高斯滤波得到滤波矩阵；将滤波矩阵与纹理幅值分布做矩阵点乘，得到无混叠的纹理频域图像；（5）结合步骤（2）图像纹理在频率阈的相位分布，对无混叠的纹理频域图像处理，得到无混叠的纹理空间域图像；然后将步骤（1）得到的图像纹理与无混叠的纹理空间域图像相减，得到混叠空间域图像；（6）将步骤（1）输入图像减去混叠伪影图像，得到消除混叠伪影之后的图像。2.如权利要求1所述的一种CT图像伪影去除方法，其特征在于：所述步骤（1）中，首先对输入图像分别做中值滤波和拉普拉斯滤波，将输入图像与中值滤波的结果相减得到图像纹理分布；输入图像通过拉普拉斯滤波可以得到图像的边缘分布；根据拉普拉斯滤波的结果，将图像纹理分布中较强的边缘设为保护区域，得到带有保护区域的图像纹理分布。3.如权利要求1所述的一种CT图像伪影去除方法，其特征在于：所述步骤（2）中，对带有边缘保护区域的图像纹理分布做二维傅里叶变换，得到图像纹理分布在频率域的相位分布和幅值分布。4.如权利要求1所述的一种CT图像伪影去除方法，其特征在于：所述步骤（4）中，滤波核为51x51像素的高斯形凹陷。5.如权利要求3所述的一种CT图像伪影去除方法，其特征在于：或者，所述步骤（5）为：将步骤（1）得到的图像纹理减去无混叠的纹理频域图像得到混叠频域图像；结合步骤（2）图像纹理在频率阈的相位分布，对混叠频域图像处理得到混叠空间域图像。6.如权利要求5所述的一种CT图像伪影去除方法，其特征在于：所述步骤（5）中的图像处理具体为对图像做二维傅里叶逆变换。2CN112785520A说明书1/3页CT图像伪影去除方法技术领域[0001]本发明涉及一种CT伪影校正方法，尤其涉及一种CT图像伪影去除方法，属于医学影像技术领域。背景技术[0002]混叠伪影是CT图像重建中非常常见的一种伪影。该伪影的产生原因一般为帧采样不足或者探测器采样不足。其在图像上的表现为图像边缘或者高密度物质的边缘出现放射状条纹。混叠伪影的存在轻则降低图像质量，重则影响临床诊断。[0003]为了解决帧采样不足导致的CT混叠伪影，常见的方法是提高投影方向的采样率。为了解决探测器采样不足导致的混叠伪影，在系统硬件上常见的方法有以下几种：一种是使用1/4探测器偏移。这种方法通过使用对向角度采集的数据来提高分辨率，减少混叠伪影。但是使用该技术时只能保证探测器视野中心的位置有较好的效果。另一种是飞焦点技术，在采集投影时焦点在两个固定位置上来回摆动。该技术也可以有效的提高分辨率，较少混叠。但是该技术要求使用高成本的球管。[0004]除了在硬件上使用1/4探测器偏移、飞焦点等技术外，还可以通过各种图像处理技术来消除混叠伪影。常见的图像处理方法包括以下几种：第一类是在投影域对数据进行处理。在投影域的处理一般是通过使用更平滑的卷积核来实现的。使用平滑的卷积核可以明显的压制投影数据中的高频成分，减少混叠伪影。该方法的缺点是会导致图像细节丢失，变