CN106780651B - 一种心脏螺旋ct图像重建方法及装置 - Google Patents

一种心脏螺旋ct图像重建方法及装置 Download PDF

Info

Publication number
CN106780651B
CN106780651B CN201611251221.1A CN201611251221A CN106780651B CN 106780651 B CN106780651 B CN 106780651B CN 201611251221 A CN201611251221 A CN 201611251221A CN 106780651 B CN106780651 B CN 106780651B
Authority
CN
China
Prior art keywords
projection data
data range
range
motion estimation
reconstruction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201611251221.1A
Other languages
English (en)
Other versions
CN106780651A (zh
Inventor
王鑫
全国涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201611251221.1A priority Critical patent/CN106780651B/zh
Publication of CN106780651A publication Critical patent/CN106780651A/zh
Application granted granted Critical
Publication of CN106780651B publication Critical patent/CN106780651B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography

Abstract

本发明实施例提供了一种心脏螺旋CT图像重建方法及装置,该方法包括:根据目标相位与门控曲线确定用于重建心脏的第一投影数据范围;根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围;根据所述第二投影数据范围与所述第三投影数据范围确定运动估计向量MVF;根据所述第一投影数据范围与所述运动估计向量MVF,得到心脏重建图像。本发明实施例通过采用上述技术方案,使得重建心脏时所需处理的数据量减少,提高图像的重建速度,并且运动估计更加准确,有效减少重建图像的运动伪影。

Description

一种心脏螺旋CT图像重建方法及装置
技术领域
本发明实施例涉及图像处理技术领域,尤其涉及一种心脏螺旋CT图像重建方法及装置。
背景技术
图像重建是指将物体外部测量的数据,通过数字处理获得物体形状信息的技术。现有的图像重建算法,在对目标物体进行成像时,若目标物体处于运动状态,直接成像会产生较大的运动伪影,伪影是指原本被扫描物体并不存在而在图像上却出现的各种形态的影像。
以进行冠脉造影检查为例,通过CT设备,以螺旋扫描的方式对人体进行扫描,得到CT扫描数据;进而采用图像重建算法对CT扫描数据进行图像处理,重建得到心脏各截面的影像,供医生使用。但是,由于心脏是持续跳动的,所以扫描过程中并不是对静态的心脏器官进行扫描,而是对运动中的心脏进行扫描,这使得各角度获得的CT扫描数据并不是针对静态物体扫描得到的。由此直接成像会存在较大的运动伪影。
为了抑制运动伪影,通常通过门控曲线选择心脏运动相对平稳时刻的数据进行重建;具体是,在CT扫描过程中会同步检测得到心脏的心电信号,反映心脏跳动的情况。采用心脏门控技术选择心脏运动相对平稳的时刻,进而将平稳时刻扫描得到的CT扫描数据进行图像重建,得到心脏图像。
但是,门控曲线选择的数据仍然存在一定的运动成分,在某些情况下,甚至无法找到运动相对平稳时刻的数据进行重建,无法有效减少重建图像中的运动伪影。
发明内容
有鉴于此,本发明实施例提供一种心脏螺旋CT图像重建方法及装置,以解决现有图像重建算法在对运动的目标进行成像时,无法有效减少运动伪影的技术缺陷。
第一方面,本发明实施例提供了心脏螺旋CT图像重建方法,包括:
根据目标相位与门控曲线确定用于重建心脏的第一投影数据范围;
根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围,其中所述第二投影数据范围在所述目标相位对应的目标投影数据位置之前采集,所述第三投影数据范围在所述目标相位对应的目标投影数据位置之后采集,并且所述第二投影数据范围与第三投影数据范围对应的重建图像在Z轴方向具有公共重建范围;
根据所述第二投影数据范围与所述第三投影数据范围确定运动估计向量MVF;
根据所述第一投影数据范围与所述运动估计向量MVF,得到心脏重建图像。
第二方面,本发明实施例还提供了一种心脏螺旋CT图像重建装置,包括:
重建数据确定模块,用于根据目标相位与门控曲线确定用于重建心脏的第一投影数据范围;
运动估计数据确定模块,用于根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围,其中所述第二投影数据范围在所述目标相位对应的目标投影数据位置之前采集,所述第三投影数据范围在所述目标相位对应的目标投影数据位置之后采集,并且所述第二投影数据范围与第三投影数据范围对应的重建图像在Z轴方向具有公共重建范围;
运动估计向量确定模块,用于根据所述第二投影数据范围与所述第三投影数据范围确定运动估计向量MVF;
重建模块,用于根据所述第一投影数据范围与所述运动估计向量MVF,得到心脏重建图像。
本发明实施例提供的技术方案,通过门控曲线确定出用于重建心脏的投影数据范围,之后选取运动估计的数据范围并确定运动估计向量MVF,在根据门控曲线确定的投影数据进行心脏重建时,利用运动估计向量进行运动估计的校正,得到最终的心脏重建图像。本发明实施例通过采用上述技术方案,使得重建心脏时所需处理的数据量减少,提高图像的重建速度,并且运动估计更加准确,有效减少重建图像的运动伪影。
附图说明
图1a为本发明实施例所适用的螺旋CT扫描设备的结构示意图;
图1b为本发明实施例一提供的一种心脏螺旋CT图像重建方法的流程示意图;
图1c为本发明实施例一提供的心动周期与门控曲线的示意图;
图2为本发明实施例二提供的一种心脏螺旋CT图像重建方法的流程示意图;
图3为本发明实施例三提供的一种心脏螺旋CT图像重建方法的流程示意图;
图4为本发明实施例四提供的一种心脏螺旋CT图像重建装置的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
在更加详细地讨论示例性实施例之前应当提到的是,一些示例性实施例被描述成作为流程图描绘的处理或方法。虽然流程图将各步骤描述成顺序的处理,但是其中的许多步骤可以被并行地、并发地或者同时实施。此外,各步骤的顺序可以被重新安排。当其操作完成时所述处理可以被终止,但是还可以具有未包括在附图中的附加步骤。所述处理可以对应于方法、函数、规程、子例程、子程序等等。
为清楚介绍本发明实施例的技术方案,首先介绍本发明实施例所适用的螺旋CT扫描设备及其扫描原理。如图1a所示为本发明实施例所适用的螺旋CT扫描设备的结构示意图,该设备包括扫描床和扫描架构成,扫描架内设置有扫描仪,扫描仪呈环形,围绕在扫描床外侧,按照设定扫描需要,扫描床以均匀速度或者可变速度,沿着扫描床的长度方向(人的身高方向,称为Z轴)水平移动。在水平移动的同时,扫描仪内部会绕Z轴旋转,形成图1a所示的螺旋扫描轨迹。
实施例一
图1b为本发明实施例一提供的一种心脏螺旋CT图像重建方法的流程示意图,该方法适用于心脏螺旋CT图像重建的情况,该方法可以由心脏螺旋CT图像重建装置执行,其中该装置可由软件和/或硬件实现,一般可集成CT成像设备中。如图1b所示,该方法包括:
步骤101、根据目标相位与门控曲线确定用于重建心脏的第一投影数据范围。
示例性的,本实施例提供的心脏螺旋CT图像重建具体为回顾式心脏扫描的重建,即螺旋CT扫描图像重建,得到的投影数据为心脏的投影数据,相位是心脏的心电图(electrocardiogram,ECG)信号中的术语,表示心脏处于当前心动周期的位置,对应心脏的运动状态。心脏的ECG信号可以被R波分成一个一个的周期,即为心动周期,若相位处于当前心动周期的40%左右,通常心脏处于收缩末期,若相位处于当前心动周期的75%左右,通常心脏处于舒张中期。目标相位是指所要重建的心脏所在心动周期的相位,可以通过扫描协议参数确定,也可以是用户通过界面设置得到的。
示例性的,在进行螺旋CT扫描时,通过球管发射的X射线来扫描人体从而采集到投影数据,球管发射一次射线得到的所有射束信息可以用一个view来表示,其中,射束信息可以包括:采集数据的时间、采集到的投影数据、扫描角度、Z轴位置等信息,当确定目标相位后,则可以根据目标相位与射束信息在时间上的对应关系,确定与目标相位对应的view,将目标相位记录为P,则其对应的射束信息为可以表示为viewP,将viewP加入门控算法中进行计算,从而确定用于重建心脏的第一投影数据的范围,可以记录为D1:[View1 i,View2 i],第一投影数据的范围是与目标相位有设定对应关系。
具体的,图1c为本发明实施例一提供的心动周期与门控曲线的示意图,如图1c所示,ECG为采集心电信号得到的心电图,ECG图中的每一个峰值位置为R波位置,即图中以“x”标识的位置,相邻两个R波之间为一个心动周期,虚线与ECG相交得到的任一交点可以作为上述目标相位,如将相位P作为目标相位,图中Gweight表示将目标相位加入门控算法得到的门控曲线,Gweight曲线中根据目标相位确定的门控范围,即Gweight曲线中的凸起位置范围W,对应的投影数据区域即为上述第一投影数据区域,虚线与Gweight门控曲线交点的Z轴位置对应上述目标投影数据位置。
步骤102、根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围。
其中,第二投影数据范围在目标相位对应的目标投影数据位置之前采集,第三投影数据范围在目标相位对应的目标投影数据位置之后采集,并且第二投影数据范围与第三投影数据范围对应的重建图像在Z轴方向具有公共重建范围。
示例性的,运动估计是指基于心脏在不同时刻的图像来估计心脏的运动趋势,根据目标相位P对应的投影数据的射束信息viewP,确定在viewP之前采集到的第二投影数据范围和在viewP之后采集到的第三投影数据范围,并且分别根据第二投影数据范围与第三投影数据范围重建的图像在Z轴上应具有公共的重建范围,从而能够根据第二投影数据范围与第三投影数据范围重建上述公共重建范围内的图像,得到相同体素在不同时刻的运动情况,即心脏的相同位置在不同时刻的运动情况,以用于后续运动估计,可以理解的是,不同时刻对应不同的扫描角度。
步骤103、根据第二投影数据范围与第三投影数据范围确定运动估计向量MVF。
示例性的,根据第二投影数据范围与第三投影数据范围对应的重建图像组,确定出运动估计的向量MVF,依据MVF来估计心脏的运动趋势。
步骤104、根据第一投影数据范围与运动估计向量MVF,得到心脏重建图像。
示例性的,根据第一投影数据范围内的投影数据进行心脏重建时,依据得到的运动估计向量MVF对重建的图像进行运动估计校正,从而得到校正后的心脏重建图像。
在计算机断层扫描成像中,通常是根据360度的扫描范围内得到的所有投影数据进行运动估计与校正,从而得到的人体断层图像。本实施例提供的技术方案,通过门控曲线确定出用于重建心脏的投影数据范围,之后选取运动估计的数据范围并确定运动估计向量MVF,在根据门控曲线确定的投影数据进行心脏重建时,利用运动估计向量进行运动估计的校正,得到最终的心脏重建图像。通过将门控曲线与运动估计校正进行结合,使得重建心脏时所需处理的数据量减少,提高图像的重建速度,并且运动估计更加准确,有效减少重建图像的运动伪影。
可选的,根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围包括:以上述目标投影数据位置为中心,按照第一预设步长分别确定位于目标投影数据位置两侧的第二投影数据范围中心和第三投影数据范围中心;根据第二投影数据范围中心和第三投影数据范围中心,以及作为区域半径值的第二预设步长,确定第二投影数据范围与第三投影数据范围。
示例性的,第一预设步长可以为与目标投影数据位置对应的扫描角度相差90度对应的投影数据的长度,可以表示为V90,则第二投影数据范围中心为viewP-V90,则第三投影数据范围中心为viewP+V90;第二预设步长可以为预设扫描角度差值对应的投影数据的距离,该预设扫描角度差值可以为60度~180度,优选的为90度,其对应的投影数据的长度可以表示为V0,则第二投影数据范围可以表示为D2:[viewp-V90-V0/2,viewp-V90+V0/2],第三投影数据范围可以表示为D3:[viewp+V90-V0/2,viewp+V90+V0/2],其中V0的取值可以通过样本统计来确定,V0越小,则灵敏度越高,运动估计越准确,抗干扰能力越弱,V0越大,则灵敏度越低,运动估计越不准确,但抗干扰能力越强,本实施例对V0取值不作具体限定,可以根据实际使用需求,平衡灵敏度与抗干扰能力来确定合适的V0,选取合适的预设扫描角度差值。
实施例二
图2为本发明实施例二提供的一种心脏螺旋CT图像重建方法的流程示意图,本实施例在上述实施例一的基础上,对确定运动估计向量MVF进行优化。如图2所示,该方法包括:
步骤201、根据目标相位与门控曲线确定用于重建心脏的第一投影数据范围。
步骤202、根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围。
其中第二投影数据范围在目标相位对应的目标投影数据位置之前采集,第三投影数据范围在目标相位对应的目标投影数据位置之后采集,并且第二投影数据范围与第三投影数据范围对应的重建图像在Z轴方向具有公共重建范围。
步骤203、分别根据第二投影数据范围与第三投影数据范围确定图像组I1与图像组I2。
示例性的,根据第二投影数据范围与第三投影数据范围进行图像重建,得到相应的图像组I1与图像组I2。
可选的,计算第二投影数据范围与第三投影数据范围对应的重建图像在Z轴方向的公共重建范围;分别重建第二投影数据范围与第三投影数据范围在公共重建范围内的图像组I1与图像组I2。
示例性的,根据螺旋反投影重建的几何映射关系,可以计算出上述D2与D3数据范围对应的重建图像的范围,D2与D3对应的重建图像在Z轴方向的重建范围分别记作Z1与Z2,Z1与Z2重合的重建范围记作Zi,则Zi即为D2与D3对应的重建图像在Z轴方向的公共范围。重建第二投影数据范围D2在Zi范围内的图像,得到图像组I1,重建第三投影数据范围D3在Zi范围内的图像,得到图像组I2。
步骤204、根据图像组I1与图像组I2确定运动估计向量MVF。
可选的,对图像组I1与图像组I2进行三维非刚性配准,确定运动估计向量MVF。
示例性的,得到的运动估计向量MVF为空间向量,该运动估计向量表示了D2与D3数据范围内的重建图像所包含的运动伪影信息。
步骤205、根据第一投影数据范围与运动估计向量MVF,得到心脏重建图像。
本实施例提供的技术方案,根据第二投影数据范围与第三投影数据范围的重建图像在Z轴方向上的公共重建范围,重建得到图像组I1与图像组I2,进而确定运动估计向量MVF,通过对用于运动估计的数据范围进行了筛选,使得运动估计更加准确,在根据门控曲线确定的投影数据进行心脏重建时,利用运动估计向量对重建图像进行校正,得到最终的心脏重建图像,通过将门控曲线与运动估计校正进行结合,使得重建心脏时所需处理的数据量减少,提高图像的重建速度,并且运动估计更加准确,有效减少重建图像的运动伪影。
实施例三
图3为本发明实施例三提供的一种心脏螺旋CT图像重建方法的流程示意图,本实施例在上述各实施例的基础上,对“根据第一投影数据范围与运动估计向量MVF,得到心脏重建图像”进行优化,如图3所示,该方法包括:
步骤301、根据目标相位与门控曲线确定用于重建心脏的第一投影数据范围。
步骤302、根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围。
其中第二投影数据范围在目标相位对应的目标投影数据位置之前采集,第三投影数据范围在目标相位对应的目标投影数据位置之后采集,并且第二投影数据范围与第三投影数据范围对应的重建图像在Z轴方向具有公共重建范围;
步骤303、根据第二投影数据范围与第三投影数据范围确定运动估计向量MVF。
步骤304、根据运动估计向量MVF,对第一投影数据范围内所包含投影数据对应的扫描角度进行插值计算,得到各扫描角度对应的运动向量。
示例性的,根据某一投影数据位置能够确定该投影数据位置对应的射束信息,即确定投影数据位置对应的view,根据射束信息则能够确定出该投影数据对应的扫描角度,假设根据上述D2与D3确定的在公共重建范围的重建图像组I1与图像组I2的对应的扫描角度相差180度,即扫描角度相差180度时对应的运动估计向量为MVF,则第一投影数据范围内的,扫描角度为45度时对应的view的运动估计向量记作MVF45,则MVF45计算公式如下:
Figure BDA0001196973240000111
在计算其他扫描角度对应的view的运动估计向量时,依据上述公式进行线性插值类推即可。
步骤305、根据各扫描角度对应的运动向量,对依据第一投影数据范围重建的图像进行运动补偿,得到心脏重建图像。
示例性的,根据第一投影数据范围重建心脏图像时,即根据第一投影数据范围D1:[View1 i,View2 i]内包含的view对应的投影数据来重建心脏图像,根据上述插值计算方法,确定第一投影数据范围内包含的各个viewi对应的运动估计向量MVFi,在进行图像重建时,将MVFi加入到投影过程中,对每个viewi对应的射线投影到图像上的位置利用MVFi进行修正,利用修正后的投影数据进行反投影重建,从而消除重建图像中的运动伪影,得到最终的心脏重建图像,其中,每个viewi对应的射线投影到图像上的位置可以以坐标形式(x,y,z)进行表示,i可以取不同的数值,如i可以为1、2、3……,来表示不同的view。
本实施例提供的技术方案,根据第二投影数据范围与第三投影数据范围确定运动估计向量MVF之后,对运动估计向量MVF进行插值计算,得到第一投影数据范围内所包含投影数据对应的各扫描角度对应的运动向量,进而对第一投影数据范围重建的图像进行运动补偿,得到最终的心脏重建图像,使得对每个扫描角度相应的重建数据均进行运动补偿,进一步确保有效减少重建图像的运动伪影。
实施例四
图4为本发明实施例四提供的一种心脏螺旋CT图像重建装置的结构示意图,该装置可由软件和/或硬件实现,一般可集成在CT成像设备中,可通过执行心脏螺旋CT图像重建方法来进行心脏螺旋CT图像重建。如图4所示,该装置包括:
重建数据确定模块401,用于根据目标相位与门控曲线确定用于重建心脏的第一投影数据范围;
运动估计数据确定模块402,用于根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围,其中第二投影数据范围在目标相位对应的目标投影数据位置之前采集,第三投影数据范围在目标相位对应的目标投影数据位置之后采集,并且第二投影数据范围与第三投影数据范围对应的重建图像在Z轴方向具有公共重建范围;
运动估计向量确定模块403,用于根据第二投影数据范围与第三投影数据范围确定运动估计向量MVF;
重建模块404,用于根据第一投影数据范围与运动估计向量MVF,得到心脏重建图像。
本实施例提供的技术方案,通过门控曲线确定出用于重建心脏的投影数据范围,之后选取运动估计的数据范围并确定运动估计向量MVF,在根据门控曲线确定的投影数据进行心脏重建时,利用运动估计向量进行运动估计的校正,得到最终的心脏重建图像。通过将门控曲线与运动估计校正进行结合,使得重建心脏时所需处理的数据量减少,提高图像的重建速度,并且运动估计更加准确,有效减少重建图像的运动伪影。
在上述实施例的基础上,运动估计数据确定模块402包括:
数据中心确定单元,用于以目标投影数据位置为中心,按照第一预设步长分别确定位于目标投影数据位置两侧的第二投影数据范围中心和第三投影数据范围中心;
运动估计数据确定单元,用于根据第二投影数据范围中心和第三投影数据范围中心,以及作为区域半径值的第二预设步长,确定第二投影数据范围与第三投影数据范围。
在上述实施例的基础上,运动估计向量确定模块403包括:
图像组确定单元,用于分别根据第二投影数据范围与第三投影数据范围确定图像组I1与图像组I2;
运动估计向量确定单元,用于根据图像组I1与图像组I2确定运动估计向量MVF。
在上述实施例的基础上,图像组确定单元包括:
公共范围确定子单元,用于计算第二投影数据范围与第三投影数据范围对应的重建图像在Z轴方向的公共重建范围;
图像组确定子单元,用于分别重建第二投影数据范围与第三投影数据范围,在公共重建范围内的图像组I1与图像组I2。
在上述实施例的基础上,重建模块404包括:
插值单元,用于根据运动估计向量MVF,对第一投影数据范围内所包含投影数据对应的扫描角度进行插值计算,得到各扫描角度对应的运动向量;
重建单元,用于根据各扫描角度对应的运动向量,对依据第一投影数据范围重建的图像进行运动补偿,得到心脏重建图像。
上述实施例中提供的心脏螺旋CT图像重建装置可执行本发明任意实施例所提供的心脏螺旋CT图像重建方法,具备执行该方法相应的功能模块和有益效果。未在上述实施例中详尽描述的技术细节,可参见本发明任意实施例所提供的心脏螺旋CT图像重建方法。
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。

Claims (8)

1.一种心脏螺旋CT图像重建方法,其特征在于,包括:
根据目标相位与门控曲线确定用于重建心脏的第一投影数据范围;
根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围,其中所述第二投影数据范围在所述目标相位对应的目标投影数据位置之前采集,所述第三投影数据范围在所述目标相位对应的目标投影数据位置之后采集,并且所述第二投影数据范围与所述第三投影数据范围对应的重建图像在Z轴方向上具有公共重建范围;
根据所述第二投影数据范围与所述第三投影数据范围确定运动估计向量MVF;
根据所述第一投影数据范围与所述运动估计向量MVF,得到心脏重建图像;
其中,根据所述第一投影数据范围与所述运动估计向量MVF,得到心脏重建图像包括:
根据所述运动估计向量MVF,对所述第一投影数据范围内所包含投影数据对应的扫描角度进行插值计算,得到各扫描角度对应的运动向量;
根据所述各扫描角度对应的运动向量,对依据所述第一投影数据范围重建的图像进行运动补偿,得到心脏重建图像。
2.根据权利要求1所述的方法,其特征在于,根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围包括:
以所述目标投影数据位置为中心,按照第一预设步长分别确定位于所述目标投影数据位置两侧的第二投影数据范围中心和第三投影数据范围中心;
根据所述第二投影数据范围中心和第三投影数据范围中心,以及作为区域半径值的第二预设步长,确定所述第二投影数据范围与所述第三投影数据范围。
3.根据权利要求1所述的方法,其特征在于,根据所述第二投影数据范围与所述第三投影数据范围确定运动估计向量MVF包括:
分别根据所述第二投影数据范围与所述第三投影数据范围确定图像组I1与图像组I2;
根据所述图像组I1与所述图像组I2确定运动估计向量MVF。
4.根据权利要求3所述的方法,其特征在于,分别根据所述第二投影数据范围与所述第三投影数据范围确定图像组I1与图像组I2包括:
计算所述第二投影数据范围与所述第三投影数据范围对应的重建图像在Z轴方向的公共重建范围;
分别重建所述第二投影数据范围与所述第三投影数据范围在所述公共重建范围内的所述图像组I1与所述图像组I2。
5.一种心脏螺旋CT图像重建装置,其特征在于,包括:
重建数据确定模块,用于根据目标相位与门控曲线确定用于重建心脏的第一投影数据范围;
运动估计数据确定模块,用于根据目标相位确定用于运动估计的第二投影数据范围与第三投影数据范围,其中所述第二投影数据范围在所述目标相位对应的目标投影数据位置之前采集,所述第三投影数据范围在所述目标相位对应的目标投影数据位置之后采集,并且所述第二投影数据范围与所述第三投影数据范围对应的重建图像在Z轴方向具有公共重建范围;
运动估计向量确定模块,用于根据所述第二投影数据范围与所述第三投影数据范围确定运动估计向量MVF;
重建模块,用于根据所述第一投影数据范围与所述运动估计向量MVF,得到心脏重建图像;
其中,所述重建模块包括:
插值单元,用于根据所述运动估计向量MVF,对所述第一投影数据范围内所包含投影数据对应的扫描角度进行插值计算,得到各扫描角度对应的运动向量;
重建单元,用于根据所述各扫描角度对应的运动向量,对依据所述第一投影数据范围重建的图像进行运动补偿,得到心脏重建图像。
6.根据权利要求5所述的装置,其特征在于,所述运动估计数据确定模块包括:
数据中心确定单元,用于以所述目标投影数据位置为中心,按照第一预设步长分别确定位于所述目标投影数据位置两侧的第二投影数据范围中心和第三投影数据范围中心;
运动估计数据确定单元,用于根据所述第二投影数据范围中心和第三投影数据范围中心,以及作为区域半径值的第二预设步长,确定所述第二投影数据范围与所述第三投影数据范围。
7.根据权利要求5所述的装置,其特征在于,所述运动估计向量确定模块包括:
图像组确定单元,用于分别根据所述第二投影数据范围与所述第三投影数据范围确定图像组I1与图像组I2;
运动估计向量确定单元,用于根据所述图像组I1与所述图像组I2确定运动估计向量MVF。
8.根据权利要求7所述的装置,其特征在于,所述图像组确定单元包括:
公共范围确定子单元,用于计算所述第二投影数据范围与所述第三投影数据范围对应的重建图像在Z轴方向的公共重建范围;
图像组确定子单元,用于分别重建所述第二投影数据范围与所述第三投影数据范围在所述公共重建范围内的所述图像组I1与所述图像组I2。
CN201611251221.1A 2016-12-29 2016-12-29 一种心脏螺旋ct图像重建方法及装置 Active CN106780651B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611251221.1A CN106780651B (zh) 2016-12-29 2016-12-29 一种心脏螺旋ct图像重建方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611251221.1A CN106780651B (zh) 2016-12-29 2016-12-29 一种心脏螺旋ct图像重建方法及装置

Publications (2)

Publication Number Publication Date
CN106780651A CN106780651A (zh) 2017-05-31
CN106780651B true CN106780651B (zh) 2021-03-09

Family

ID=58928117

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611251221.1A Active CN106780651B (zh) 2016-12-29 2016-12-29 一种心脏螺旋ct图像重建方法及装置

Country Status (1)

Country Link
CN (1) CN106780651B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107348969B (zh) * 2017-07-11 2020-07-14 上海联影医疗科技有限公司 一种pet数据处理方法、系统及pet成像设备
CN109345606B (zh) * 2018-09-27 2023-01-03 上海联影医疗科技股份有限公司 心脏图像重建方法、装置、计算机设备和可读存储介质
US10950016B2 (en) 2018-06-11 2021-03-16 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for reconstructing cardiac images
CN111150419B (zh) * 2020-02-14 2022-12-09 赛诺威盛科技(北京)股份有限公司 螺旋ct扫描重建图像的方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102144927A (zh) * 2010-02-10 2011-08-10 清华大学 基于运动补偿的ct设备和方法
CN104331914A (zh) * 2014-09-18 2015-02-04 东莞南方医大松山湖科技园有限公司 重建心脏相位图像的方法和装置
CN106232009A (zh) * 2014-02-21 2016-12-14 三星电子株式会社 断层扫描设备和由断层扫描设备重构断层扫描图像的方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2002950210A0 (en) * 2002-07-11 2002-09-12 Mediaware Solutions Pty Ltd Mosaic construction from a video sequence
US7630528B2 (en) * 2004-03-02 2009-12-08 Koninklijke Philips Electronics N.V. Motion compensation
US9129426B2 (en) * 2010-08-31 2015-09-08 General Electric Company Motion compensation in image processing

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102144927A (zh) * 2010-02-10 2011-08-10 清华大学 基于运动补偿的ct设备和方法
CN106232009A (zh) * 2014-02-21 2016-12-14 三星电子株式会社 断层扫描设备和由断层扫描设备重构断层扫描图像的方法
CN104331914A (zh) * 2014-09-18 2015-02-04 东莞南方医大松山湖科技园有限公司 重建心脏相位图像的方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CT成像技术的发展;石明国 等;《中国医学装备》;20070430;第4卷(第4期);第56-60页 *
Enhanced temporal resolution at cardiac CT with a novel CT image reconstruction algorithm: Initial patient experience;Apfaltrer Paul 等;《European Journal of Radiology》;20130228;第82卷(第2期);第270-274页 *

Also Published As

Publication number Publication date
CN106780651A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
US7426256B2 (en) Motion-corrected three-dimensional volume imaging method
CN110751702B (zh) 图像重建方法、系统、装置及存储介质
US8184883B2 (en) Motion compensated CT reconstruction of high contrast objects
CN110009709B (zh) 医学图像成像方法和系统
US9576391B2 (en) Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus
CN106780651B (zh) 一种心脏螺旋ct图像重建方法及装置
CN106920265B (zh) 计算机断层扫描图像重建方法及装置
EP3107457B1 (en) Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus
US20080267455A1 (en) Method for Movement Compensation of Image Data
JP4828920B2 (ja) 3次元画像処理装置
JP2012233896A (ja) 医用画像を運動補正する方法及び装置
CN108876730B (zh) 校正运动伪影的方法、装置及设备和存储介质
JP2005296652A (ja) 運動する対象の構造データの取得装置
US20140126685A1 (en) Angiographic examination method
JP5209210B2 (ja) X線ct装置
US8855391B2 (en) Operating method for an imaging system for the time-resolved mapping of an iteratively moving examination object
US20230130015A1 (en) Methods and systems for computed tomography
JP4777164B2 (ja) 心拍位相決定装置、プログラム及びx線診断装置
JP2007260398A (ja) 心臓のコンピュータ断層撮影表示における階段状アーチファクトの低減方法およびメモリ媒体
CN110728730A (zh) 图像重建方法、装置、ct设备及ct系统
US20120275656A1 (en) Method for supplying a 3d x-ray image data record for a moving object with highly absorbent material
CN115067980A (zh) 扫描时相的确定方法、装置、计算机设备和存储介质
CN112330800A (zh) 图像重建方法、装置、控制台设备及ct系统
US10537297B2 (en) Radiation tomography apparatus and program for cardiac-gated imaging
Sindel et al. Respiratory motion compensation for C-arm CT liver imaging

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

GR01 Patent grant
GR01 Patent grant