专利名称：基于序列相似性插值的部分k空间图像重构方法
技术领域：
本发明属于图像处理技术领域，涉及医学图像处理方法，可用于MRI序列图像重构。
背景技术：
部分K空间图像重构是为了加快磁共振图像成像速度而提出的一种减少数据采集量来重构高分辨率图像的问题。为了解决这一问题，有很多经典的方法被提出:第一种是补零法，即未采集的K空间数据用零填补，然后做傅里叶反变换得到图像空间的成像方法，这种成像方法可以提高成像速度，缺陷是图像中有伪影。第二种是相位校正法，这类方法假设磁共振图像空间的相位呈缓慢变化状态。它用部分低频K空间数据重构的图像来估计相位，并用于相位校正，从而达到利用对称性来补充未采集的K空间数据的目的。Stark，H.提出的POCS方法就是这类方法中最常见的方法。但是由于磁共振图像相位缓慢变化的条件在整个图像空间中常常难以满足，导致相位估计误差大，引起较大的重构误差，以至到目前为止无法在医学临床中应用。第三种是信号估计法，这种方法运用信号估计理论，利用部分K空间数据的插值、外推及多目标优化等方法来获取未采集的K空间数据。M.Funderer于1989年提出了图像最大似然的方法，E.M.Haacke在同年提出的约束方法，都是是这类方法的典型方法，这类方法成像效果远劣于相位校正方法。第四种是奇异谱模型法，这类方法的思想是基于任何信号都可以用奇异函数的加权和表示的事实，建立新的图像表达模型，由部分K空间数据提取模型参数，然后按模型和参数重构出完整的K空 间，这种方法多数情况下优于相位校正法和信号估计法。第五种是压缩感知的方法，该方法利用小波，有限差分，字典学习等对重构图像进行稀疏表示，这类方法的效果多数情况下都优于其他方法。上述部分K空间重构方法通常都需要采样率达到30%以上，才能获得较好的重构效果。但是采样率越高，所需采集时间越长，被成像者要长时间滞留在成像仪器内，由于被成像者的运动将会使图像产生运动模糊。然而，进一步降低采样率会造成重构图像产生伪影和图像分辨率低的问题。
发明内容
本发明的目的在于针对上述已有技术的不足，提出一种基于序列相似性插值的部分K空间图像重构方法，以降低数据采样率，提高重构图像的质量。为实现上述目的，本发明包括如下步骤:(I)采集15幅K空间数据Qi, i = 1，2，…，15，对Qi做10倍下采样，得到对应的部分K空间数据Pi = Q1.*M，M为下采样矩阵作傅里叶反变换，得到低分辨率图像Li，并将其作为样本图像；(2)输入待重构的部分K空间数据Pt，对该部分K空间数据Pt采用补零法处理，得到初始重构图像Lt ；
11/,-/..1l(3)分别求初始重构图像Lt和15幅样本图像Li的误差:
权利要求
1.一种基于序列相似性插值的部分K空间图像重构方法，包括如下步骤: (1)采集15幅K空间数据Qi,i = 1,2, -,IS,对Qi做10倍下采样，得到对应的部分K空间数据Pi = Q1.*M，M为下采样矩阵；对Pi作傅里叶反变换，得到低分辨率图像Li,并将其作为样本图像； (2)输入待重构的部分K空间数据Pt，对该部分K空间数据Pt采用补零法处理，得到初始重构图像Lt ； (3)分别求初始重构图像Lt和15幅样本图像Li的误差
2.根据权利要求1所述的基于序列相似性插值的部分K空间图像重构方法，其中所述步骤(I)中的下采样矩阵M是一个随机矩阵，矩阵的大小与K空间数据Qi的大小相同，其中的元素只包含O和1,0和I的个数比例为9:1,且这些元素随机分布。
3.根据权利要求1所述的基于序列相似性插值的部分K空间图像重构方法，其中步骤(2)所述的对部分K空间数据Pt采用补零法处理，是对未采集的数据补零后，再进行傅里叶反变换，得到重构图像。
全文摘要
本发明公开了一种基于序列相似性插值的部分K空间图像重构方法，主要解决现有方法在10倍下采样的情况下重构图像质量下降比较严重的问题。其主要步骤为采集15幅K空间数据Qi，i＝1,2,…，15;对Qi做10倍下采样，再做傅里叶反变换，得到样本图像Li；输入待重构的部分K空间序列数据Pt，对其用补零法处理得到初始重构图像Lt；求初始重构图像Lt与样本图像Li的误差；用误差较小的五幅样本图像Lj对应的K空间数据Qj分别对输入的部分K空间数据Pt进行补充，取平均值，再做傅里叶反变换，得到重构图像。本发明在10倍下采样的条件下能够有效提高重构图像质量，可用于多个部位的MRI序列图像重构。
文档编号G06T5/50GK103236049SQ201310161378
公开日2013年8月7日 申请日期2013年5月5日 优先权日2013年5月5日
发明者缑水平, 刘芳, 唐晓, 焦李成, 盛珂, 吴建设, 王爽, 马文萍, 马晶晶 申请人:西安电子科技大学