[发明专利]CT图像中校正金属伪影的方法和系统

摘要

本发明涉及一种CT图像中校正金属伪影的方法和系统，上述CT图像中校正金属伪影的方法包括对未校正CT图像中的金属区域进行分割，得到只含有金属信息的图像；对未校正CT图像和只含有金属信息的图像进行限值处理，得到未校正限值图像和金属信息限值图像；对未校正限值图像和金属信息限值图像进行前向投影，得到p_original和p_metal；利用p_metal定位出p_original中的金属区域，利用p_original中金属区域两端外的信息插值估计金属区域内的信息，根据插值估计后的投影数据重建CT图像，对重建CT图像进行滤波处理，构造先验图像；对先验图像进行前向投影，得到p_prior，根据p_prior对p_original中的金属区域进行插值校正，得到p_correct2；对p_correct2进行重建，得到I_correct2，将I_correct2与只含有金属信息的图像融合，校正CT图像。

主权项

一种CT图像中校正金属伪影的方法，其特征在于，包括如下步骤：对未校正CT图像中的金属区域进行分割，得到只含有金属信息的图像；分别对未校正CT图像和只含有金属信息的图像进行限值处理，得到未校正限值图像和金属信息限值图像；分别对未校正限值图像和金属信息限值图像进行前向投影，得到p_original和p_metal；其中，p_original为未校正限值图像对应的投影数据，p_metal为金属信息限值图像对应的投影数据；利用p_metal定位出p_original中的金属区域，并利用p_original中金属区域两端外的信息插值估计金属区域内的信息，根据插值估计后的投影数据重建CT图像，并对重建CT图像进行滤波处理，构造先验图像；利用p_metal定位出p_original中的金属区域，并利用p_original中金属区域两端外的信息插值估计金属区域内的信息，根据插值估计后的投影数据重建CT图像，并对重建CT图像进行滤波处理，构造先验图像的过程包括：利用非零判断运算计算出投影数据p_metal中的金属区域；根据p_metal中的金属区域定位出p_original中的金属区域；其中，p_original为M行N列的二维矩阵，p_original中第i行为p_original的第i个投影，p_original的第i个投影下的金属区域起点为i_start，终点为i_end，i为大于1且小于M的整数；对p_original的第i行，利用第i_start‑1个元素和第i_end+1个元素，对第i_start个元素到第i_end个元素之间的数据进行线性插值，得到线性插值后的投影数据p_line；计算出p_metal中每行中金属信息的起点和终点之间的距离，获取距离最小的第min行，获取第min行中金属信息的起点min_start，终点min_end；其中，p_metal为M行N列的二维矩阵；计算出投影数据p_line中第min行的数据之和sum_min；计算p_line中除第min行外的第i行数据的和sum_i，计算出第i行的sum_i与第min行的sum_min的差diff_sum_i；对于除第min行外的其它第i行，以i_start为起点以i_end为终点，构造叠加函数；其中，所述叠加函数为：fi(j)=0,ifj<i_start0,ifj>i_enddiff_sum_i(i_end-i_start)2πsin(π(j-i_start)i_end-i_start),otherwise,]]>其中，j为大于等于1，且小于等于N的整数；分别将函数fi(j)的N个值一一对应叠加至投影数据p_line中除第min行外的第i行数据上，形成修正后的投影数据p_line；对修正后的投影数据p_line每一行的i_start和i_end处进行高斯光滑滤波，得到用于重建的投影数据p_correct1；对投影数据p_correct1重建，得到图像I_correct1，对所述图像I_correct1进行具有保护图像结构边缘功能的滤波处理，得到先验图像；对先验图像进行前向投影，得到投影数据p_prior，根据p_prior对p_original中的金属区域进行插值校正，得到投影数据p_correct2；对p_correct2进行重建，得到重建图像I_correct2，将I_correct2与只含有金属信息的图像融合，得到校正后的CT图像。