一种基于运动估计的肺4d-ct图像的超分辨率冠矢状面图像重建方法
【专利摘要】一种基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，依次包括：（1）读取由多个不同相位的肺部3D图像组成的肺部4D-CT图像数据；（2）根据肺部4D-CT图像数据，对每个相位提取同一肺部部位对应的冠矢状面图像；（3）采用基于完全搜索块匹配算法估计不同“帧”肺部冠矢状面图像之间的运动矢量场；（4）以步骤（3）得到的运动矢量场为基础，采用迭代反投影法重建高分辨率肺4D-CT冠矢状面图像。本发明获得的肺4D-CT图像的超分辨率冠矢状面重建图像的分辨率有明显的提高，局部放大图中，肺实质中的血管和周边组织的亮度和清晰度有明显增强,可以克服由采集时间和放射剂量造成的图像分辨率低的限制，能够有效引导肺癌的精确放射治疗。
【专利说明】—种基于运动估计的肺4D-CT图像的超分辨率冠矢状面图 像重建方法
【技术领域】
[0001]本发明涉及医学图像处理【技术领域】，具体涉及一种基于运动估计的肺4D-CT图像 的超分辨率冠矢状面图像重建方法。
【背景技术】
[0002]肺4D-CT图像能够提供全面的高精度放射治疗呼吸运动表征。在肺4D-CT图像数 据中，由于有多个相位的图像，通常是10-20个，通过各相位图像有助于获取肺部呼吸运动 信息，是放射治疗目标精确定位的关键，因此肺4D-CT技术在肺肿瘤精确放射治疗中发挥 着越来越重要的作用。
[0003]肺4D-CT数据的获取，通常是根据床位和肺体积，将多个自由呼吸的3D-CT数据段 排序而得。然而，由于CT固有的高剂量照射，沿纵向(通常命名为Z轴方向)的密集采样往 往是不实际的，从而导致所获得的肺4D-CT数据的层间分辨率远低于层内分辨率，造成数 据显著的各项异性。
[0004]因此，在对各相位3D数据进行冠矢状面观察时，为了获得正确比例的图像，需要 根据3D数据的层间分辨率和层内分辨率的比例，沿Z轴方向进行插值放大。常用的插值方 法为最近邻或双线性插值法，但是，这些方法都会导致图像模糊，尤其是当层间分辨率与 层内分辨率比例差别比较大时图像模糊更严重。
[0005]因此，针对现有技术不足，提供一种基于运动估计的肺4D-CT图像的超分辨率冠 矢状面图像重建方法以克服现有技术不足甚为必要。
【发明内容】
[0006]本发明的目的在于针对现有技术不足，提供一种基于运动估计的肺4D-CT图像的 超分辨率冠矢状面图像重建方法，该方法能够提高肺4D-CT图像的冠矢状面图像的分辨率。
[0007]本发明的上述目的通过如下技术方案实现。
[0008]一种基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，依次包括 如下步骤，
(1)读取由多个不同相位的肺部3D图像组成的肺部4D-CT图像数据；
(2)根据肺部4D-CT图像数据，对每个相位提取同一肺部部位对应的冠矢状面图像；
(3)估计不同“帧”肺部冠矢状面图像之间的运动矢量场；
(4)以步骤(3)得到的运动矢量场为基础，重建高分辨率肺4D-CT冠矢状面图像。
[0009]上述步骤(3)是采用基于完全搜索块匹配算法估计不同“帧”肺部冠矢状面图像之 间的运动矢量场。
[0010]上述步骤(3 )具体包括:
(3.1)在当前帧中选取一子块，根据最小绝对误差匹配准则，在参考帧的给定的搜索区域内找出与当前帧中的当前块最相似的块作为匹配块；
根据匹配块与当前块的相对位置计算运动位移作为当前块的运动矢量，所述运动矢量亦为全局最优的相对运动矢量；
(3.2)所述最小绝对误差匹配准则如下:
【权利要求】
1.一种基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，其特征在于:依次包括如下步骤，(1)读取由多个不同相位的肺部3D图像组成的肺部4D-CT图像数据；(2)根据肺部4D-CT图像数据，对每个相位提取同一肺部部位对应的冠矢状面图像；(3)估计不同“帧”肺部冠矢状面图像之间的运动矢量场；(4)以步骤(3)得到的运动矢量场为基础，重建高分辨率肺4D-CT冠矢状面图像。
2.根据权利要求1所述的基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，其特征在于:所述步骤(3)是采用基于完全搜索块匹配算法估计不同“帧”肺部冠矢状面图像之间的运动矢量场。
3.根据权利要求1所述的基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，其特征在于:上述步骤(3)具体包括:(3.1)在当前帧中选取一子块，根据最小绝对误差匹配准则，在参考帧的给定的搜索区域内找出与当前帧中的当前块最相似的块作为匹配块；根据匹配块与当前块的相对位置计算运动位移作为当前块的运动矢量，所述运动矢量亦为全局最优的相对运动矢量；(3.2)所述最小绝对误差匹配准则如下:
4.根据权利要求1所述的基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，其特征在于:所述步骤(4)具体是采用迭代反投影法重建高分辨率肺4D-CT冠矢状面图像。
5.根据权利要求4所述的基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，其特征在于:所述步骤(4)具体是包括:(4.1)将需要重建的原始低分辨率图像與插值放大为初始高分辨率图像p>s),n为迭代次数；(4.2)根据退化模型将初始高分辨率图像/?'模拟成像过程得到低分辨图像的集合c4wh，k = i K &表示原始序列低分辨率图像的数量； 所述退化模型具体为: I1=IlIfti, i=yjPji； 其中:1?表示IB幅低分辨率图像中的第fc幅，/表示需要重建的初始高分辨率图像； I?表示下采样矩阵； 办是系统加性噪声； Mfc表示几何变换,为由运动估计求得的运动矢量； Sk表示模糊矩阵，是由光学系统本身、成像系统与原始场景的相对运动，以及低分辨率传感器的点扩散函数而造成的； 具体的，在第次迭代过程中，Irin)的成像过程由退化模型模拟得到: TlA 其中:咖表示第次迭代过程中假设的高分辨率图像di0表示次迭代后由退化模 fTlTl型得到的低分辨率图像；rfc表示从/到Λ的二维几何变换，即为步骤(3)获得的运动矢量；A是高斯模糊算子；A S是下采样算子； (4.3)判断误差? 是否达到最小值，如果达到最小值，停止迭代，以当前估计的高分辨率图像/(R)作为最终所求的超分辨率图像； 如果误差未达到最小值，则进入步骤(4.4)； (4.4)根据误差对当前高分辨率图像进行更新，更新过程具体如下式:
i户叫=/w + 泛 ιΤ.— sf) τ s)?p)其中 表示上采样算子；P表示背投影算子，由It.和决定； (4.5)将更新后的高分辨率图像作为初始高分辨率图像，进入步骤(4.2)。
6.根据权利要求5所述的基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，其特征在于:所述步骤(4.3)中判断误
7.根据权利要求5所述的基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，其特征在于:所述步骤(4.3)中判断误差gk — j是否达到最小值具体是通过判断是否达到最大迭代次数M来判定的,当迭代次数I达到M时，判定误差达到最小值，否则判定误差未达到最小值。
8.根据权利要求5所述的基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，其特征在于:最大迭代次数
9.根据权利要求5所述的基于运动估计的肺4D-CT图像的超分辨率冠矢状面图像重建方法，其特征在于:最大迭代次数M设置为3。
【文档编号】G06T7/20GK103440676SQ201310350915
【公开日】2013年12月11日 申请日期:2013年8月13日 优先权日:2013年8月13日
【发明者】张煜, 肖珊 申请人:南方医科大学