[发明专利]一种医学锥束CT图像中运动伪影修正方法

权利要求书

1.一种医学锥束CT图像中运动伪影修正方法，其特征在于，包括如下步骤：

S1：初始图像重建，利用病人或者模体锥束CT投影数据重建初始图像f₀；

S2：区域分割提取，利用图像分割技术和图像特征提取技术对所述初始图像f₀进行处理，获取气泡区域f_gas、骨性组织区域f_bone、阴影伪影与条状伪影存在的主要区域f_r；

S3：阴影伪影消除，利用线性插值技术对所述气泡区域f_gas、所述阴影伪影与条状伪影存在的主要区域f_r进行处理，获取阴影伪影修正后的第一锥束CT图像；

所述步骤S3具体包括如下步骤：

S31：将局部Radon变换作用于所述主要区域f_r得到第一虚拟投影，即将图像域转换到虚拟投影域；

S32：将Radon变换作用于所述气泡区域f_gas，获取含有阴影伪影污染的虚拟投影数据；

S33：在所述第一虚拟投影域中，利用线性插值技术替换含有阴影伪影污染的虚拟投影数据；

S34：重建替换后的虚拟投影数据，对重建后的虚拟投影数据进行低通滤波后得到f_mask1，将f_mask1加到所述初始图像f₀中，获取阴影伪影修正后的第一锥束CT图像，即实现锥束CT图像的阴影伪影修正；

S4：条状伪影消除，利用频率分割技术对阴影伪影修正后的第一锥束CT图像、所述气泡区域f_gas、所述骨性组织区域f_bone进行处理，获取条状伪影修正后的第二锥束CT图像，实现锥束CT图像的运动伪影消除；

所述步骤S4具体包括如下步骤：

S41：对步骤S3中获取的第一锥束CT图像进行频率分割，获取第一锥束CT图像的中频成份；

S42：将Radon变换作用于所述中频成份、所述气泡区域f_gas 和所述骨性组织区域f_bone ，获取第二虚拟投影；

S43：选取中频成份投影中的气泡成份和骨性组织成份，选取到的成份为条状伪影成份；

S44：重建第二虚拟投影中选取的条状伪影成份，得到fmask2，将步骤S3中获取的第一锥束CT图像减去fmask2，获取条状伪影修正后的第二锥束CT图像，实现锥束CT图像的运动伪影消除。

2.根据权利要求1所述的医学锥束CT图像中运动伪影修正方法，其特征在于，所述步骤S1中：采用滤波反投影技术重建所述初始图像f₀。

3.根据权利要求1所述的医学锥束CT图像中运动伪影修正方法，其特征在于，所述步骤S1中：所述初始图像f₀中的运动伪影包括阴影伪影和条状伪影。

4.根据权利要求1所述的医学锥束CT图像中运动伪影修正方法，其特征在于，所述步骤S2中：所述图像分割技术采用阈值分割技术，所述图像特征提取技术采用图像腐蚀和图像膨胀技术。

5.根据权利要求1所述的医学锥束CT图像中运动伪影修正方法，其特征在于，所述步骤S34中：采用高斯滤波器对所述投影数据进行低通滤波处理。

6.根据权利要求1所述的医学锥束CT图像中运动伪影修正方法，其特征在于，所述步骤S41中：所述第一锥束CT图像的中频成份包括条状伪影、图像边界以及图像的纹理特征。