CN103750854B - 高效迭代四维心脏锥束ct重建 - Google Patents
高效迭代四维心脏锥束ct重建 Download PDFInfo
- Publication number
- CN103750854B CN103750854B CN201310738121.1A CN201310738121A CN103750854B CN 103750854 B CN103750854 B CN 103750854B CN 201310738121 A CN201310738121 A CN 201310738121A CN 103750854 B CN103750854 B CN 103750854B
- Authority
- CN
- China
- Prior art keywords
- voxel
- phase volume
- phase
- volume
- projection
- 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.)
- Expired - Fee Related
Links
- 230000000747 cardiac effect Effects 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000012360 testing method Methods 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 12
- 230000005670 electromagnetic radiation Effects 0.000 claims description 11
- 238000002591 computed tomography Methods 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 4
- 238000004154 testing of material Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims 4
- 230000007423 decrease Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 16
- 230000000875 corresponding effect Effects 0.000 description 9
- 230000005855 radiation Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000004590 computer program Methods 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000008521 reorganization Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Medical Informatics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Mathematical Optimization (AREA)
- Pulmonology (AREA)
- Mathematical Analysis (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Algebra (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Mathematical Physics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Generation (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
如果使用锥束重建近似方法,那么当前的高端CT系统及未来的CT系统逐渐增大的锥角导致图像质量下降。根据本发明典型实施例,提供迭代四维心脏CT重建,其中从四维数据集中选择相位体积,每个在不同的相点具有同样的空间体积。然后把这些相位体积内相应的体素正向投影到同样的投影上。在计算不同投影后,更新这些体素。这可提供高效执行的迭代四维心脏锥束重建。
Description
本申请是申请号为200680033400.6的题为“高效迭代四维心脏锥束CT重建”的中国专利申请的分案申请。
本发明涉及X射线成像领域。尤其,本发明涉及用于检查感兴趣对象的检查装置、图像处理设备、检查感兴趣对象的方法、计算机可读介质及程序单元。
如果使用锥束重建近似方法,那么当前的高端CT系统及未来的CT系统的逐渐增大的锥角导致图像质量下降。对于心脏锥束CT的情况,精确重建是不可行的,并且因此,迭代重建方法对这种应用是潜在的候选。它们为大锥角提供出色的图像质量且它们能增加重建后的图像的信噪比(SNR)。
但是,这种迭代重建方法的应用由于它们的高计算消耗而受到限制。尤其是对于运动对象的四维重建而言,这个问题是严重的,因为其中必须针对大量相点重建整个体积。
因此,需要具有一种四维数据集的高效迭代图像重建。
根据本发明的典型实施例,提供用于检查感兴趣对象的检查装置,检查装置包括采集单元,用于采集感兴趣对象的四维数据集;计算单元,适于从四维数据集中选择第一相位体积和第二相位体积,基于第一组系数正向投影第一相位体积的第一体素,基于第二组系数正向投影第二相位体积的第二体素,基于第三组更新系数更新正向投影的第一体素,并且基于第四组更新系数更新正向投影的第二体素,其中第一体素和第二体素具有相同的空间坐标但相点不同。
因此,根据本发明的这个典型实施例,提供一种检查装置,其适于执行感兴趣对象的四维数据集的迭代重建,该重建可在硬件中高效地实现。因为第一体素和第二体素都具有同样的空间坐标而仅仅相点不同(因为它们是从不同的相位体积选择的),两种体素的正向投影以及两种体素的反向投影或更新可高效地执行。正向投影和更新体素具有同样的空间坐标(其中体素的正向投影可在第一步骤中执行,并且,其中两种体素的更新可在第二步骤中执行),这使得计算消耗下降。
根据本发明的另一个典型实施例,计算单元还适于计算加权函数,产生对应第一相位体积的第一权数和对应第二相位体积的第二权数。
然后,根据本发明的这个典型实施例,第一相位体积的第一权数和第二相位体积的第二权数可在一个计算步骤中得以计算(并且基本同时)。
根据本发明的另一个典型实施例,计算单元还适于计算用于将第一相位体积的第一体素映射到投影上的第一组系数,和用于将第二相位体积的第二体素映射到投影上的第二组系数。
因此,用于第一体素的系数和用于第二体素的系数可在同样的计算步骤中计算。因为两种体素具有同样的空间坐标,这个计算步骤是非常高效的。
根据本发明的另一个典型实施例,计算单元还适于基于第一投影和第二投影计算差分投影(differenceprojection)并基于该差分投影计算第一组更新系数和第二组更新系数。
因此,根据本发明的典型实施例,可对于相位相关投影的整个堆栈计算差分投影。
根据本发明的另一个典型实施例,用相应的心脏权数乘以该差分投影。
根据本发明的另一个典型实施例,检查装置还包括心电图单元,适于采集心电图数据,其中,加权函数为用于第一相位和第二相位的心脏加权函数。
因此,根据本发明的这个典型实施例,加权函数与感兴趣对象(其此时为患者)的心搏对应。这提供了运动伪影的高效减少。
根据本发明的另一个典型实施例,加权函数为来自心电图数据、扫描参数、视场大小、和视场位置的投影数量中的至少一个的函数。
因此,根据本发明的这个典型实施例,加权函数至少反映了心电图数据、数据采集过程的扫描参数、视场大小或位置。
根据本发明的另一个典型实施例,第一相位体积和第二相位体积具有相同的空间体积但相点不同。
例如,两种相位体积大小都为(nx,ny,nz),其中nx(ny,nz)描述了x(y,z)大小中的体素数量。但是在不同相点(时期)采集这两种相位体积。
根据本发明的另一个典型实施例,第一权数和第二权数大于零。
因此,根据本发明的这个典型实施例,仅仅从多相位数据集中选择那些相位体积,其对于该投影具有大于零的心脏权数并定义了相应的相位相关投影的集合。
根据本发明的另一个典型实施例,检查装置还包括电磁辐射源,适于向感兴趣对象发射电磁辐射;和准直器,布置在电磁辐射源与采集单元之间,其中,准直器适于使电磁辐射源发射的电磁辐射束准直,以形成扇束或锥束。
根据本发明的另一个典型实施例,检查装置作为心脏锥束计算机断层摄影装置。
而且,根据本发明的另一个典型实施例,四维数据集包括第一子集和第二子集,其中,第一和第二子集中每个的大小都与用于非门限迭代重建的子集大小乘以平均心动周期长度并除以平均门限窗口宽度对应。
这产生了子集大小,其提供与非门限重建相似的收敛速度。如果非门限重建的子集大小关于例如图像质量是优化的,其还适于门限重建。
根据本发明的另一个典型实施例,检查装置构造为包括行李检验装置、医疗应用装置、材料测试装置和材料科学分析装置的组中的一个。本发明的应用领域可以是材料科学分析,因为本发明限定的功能性允许快速、高效及高精度分析材料。
根据本发明的另一个典型实施例,提供检查感兴趣对象的方法,该方法包括采集感兴趣对象的四维数据集,从四维数据集中选择第一相位体积和第二相位体积,基于第一组系数正向投影第一相位体积的第一体素,基于第二组系数正向投影第二相位体积的第二体素,基于第三组更新系数更新正向投影的第一体素,基于第四组更新系数更新正向投影的第二体素,其中,第一体素和第二体素具有相同的空间坐标但相点不同。
根据本发明的另一个典型实施例,提供了检查感兴趣对象的图像处理设备,该图像处理设备包括存储器,用于存储感兴趣对象的四维数据集。而且,该图像处理设备包括计算单元,其适于执行上述方法步骤。
因此,提供图像处理设备,其适于基于四维数据集执行感兴趣对象的高效迭代重建。
根据本发明的另一个典型实施例,提供一种计算机可读介质,用于检查感兴趣对象的计算机程序存储在其中,当被处理器执行时,其适于执行上述方法步骤。
而且,本发明涉及检查感兴趣对象的程序单元,其存储在计算机可读介质上。该程序单元适于执行上述方法步骤。
程序单元优选地装载进数据处理器的工作存储器中。因此,数据处理器被装备以实现本发明方法的典型实施例。可用任意的编程语言编写计算机程序,例如C++并且存储在计算机可读介质上,例如CD-ROM。而且,计算机程序可从网络得到,例如WorldWideWeb,从中可将计算机程序下载到图像处理单元或处理器,或任何适合的计算机。
能够看出,作为本发明典型实施例的要点,从多相位数据集中选择相位体积,其中,每个相位体积具有相同的空间体积但是在不同的相点。在这些相位体积中的每个中,由于具有相同的空间坐标但是不同的相点,识别出一个体素。然后从不同的心脏相位体积将该体素集正向投影到定义为相位阵列的投影上。在计算了不同投影后,执行反向投影(再次对于整个体素集合,以同样的空间坐标并行)。
通过参照此后描述的实施例将清楚本发明的这些及其他方面。
将在下面参照附图描述本发明的典型实施例。
图1示出了代表根据本发明典型实施例的CT扫描器系统的简化视图;
图2示出了表示根据本发明典型实施例的方法的流程图;
图3示出了用于单个体素及所选的投影的投影几何结构的示意性表示;
图4示出了心脏加权函数的示意性表示;
图5示出了根据本发明的图像处理设备的典型实施例,其用于执行根据本发明的方法的典型实施例。
附图中的视图是示意性的。在不同的附图中,相似或相同的元件具有相同的参考标号。
图1示出了根据本发明典型实施例的CT扫描器系统的典型实施例。参照这个典型实施例,将针对医疗成像中的应用来描述本发明。但是,应当知道,本发明不限于这种应用,而是还可应用在行李检验,或其他工业应用领域,例如材料测试。
图1所示的计算机断层摄影装置100是锥束CT扫描器。图1所示的CT扫描器包括构架101,其绕旋转轴102旋转。构架101由马达103驱动。参考标号104代表辐射源,例如X射线源,其,根据本发明的一个方面,发射多波长辐射。
参考标号105代表孔径系统,其将从辐射源发射的辐射束整形为锥形辐射束106。锥束106被导向,从而其穿透布置在构架101的中央中,例如CT扫描器的检查区域中的感兴趣对象107,并入射到探测器108上。如从图1可知,探测器108布置在构架101上与辐射源104相对的一侧,从而探测器108的表面被锥束106覆盖。如图1所示,探测器108包括多个探测器元件123,每个元件能够以能量解析方式探测穿透感兴趣对象107的X射线或单独的光子。
在扫描感兴趣对象107期间,辐射源104、孔径系统105及探测器108沿构架101以箭头116指示的方向旋转。为了使带有辐射源104、孔径系统105及探测器108的构架101旋转,将马达103连接到马达控制单元117,其连接到计算或确定单元118。
在图1中,感兴趣对象107可以是放在传送带119上的患者或行李物品。在扫描感兴趣对象107期间,当构架101绕行李物品107旋转时,传送带119沿平行于构架101的旋转轴102的方向使感兴趣对象107移位。通过这种方式,可沿螺旋扫描路径扫描感兴趣对象107。在扫描期间还可停止传送带119由此测量单个切片。代替传送带119,例如,在感兴趣对象107为患者的医疗应用中,可采用可移动的平台。但是,应当知道,在所有描述的情况中,还可能通过以源-探测器装置的频率的两倍频率周期性地向后及向前移动平台,来执行其他扫描路径,例如鞍形轨迹。
探测器108连接到计算单元118。计算单元118接收探测结果,即来自探测器108的探测器元件123的读数,并且基于该读数确定扫描结果。而且,计算单元118与马达控制单元117通信,以便用马达103和120协调构架101与传送带119的运动。
根据本发明典型实施例,计算单元118适于迭代地重建感兴趣对象的图像。重建单元118产生的重建图像通过接口122输出到显示器(图1未示出)。
计算单元118由数据处理器实现,以处理来自探测器108的探测器元件的读数。
而且,由图1可知,重建单元118连接到扬声器121,例如,其用于在行李物品107中的探测到可疑材料时自动地输出警报。
用于检查感兴趣对象107的计算机断层摄影装置100包括探测器108,其具有以类似矩阵方式布置的多个探测器元件123,每个都适于探测X射线。而且,计算机断层摄影装置100包括确定单元或适于重建感兴趣对象107的图像的重建单元118。
计算机断层摄影装置100包括适于向感兴趣对象107发射X射线的X射线源104。设置在电磁辐射源104与探测元件123之间的准直器105适于准直从电磁辐射源104发射的电磁辐射束,以形成锥束。探测元件123构成多层探测器阵列108。
图2示出了根据本发明典型实施例的方法。由于有限的时间分辨率,在心脏锥束重建的情况下,所重建的用于产生四维数据集的心脏相位的数量大于心动周期的长度除以平均门限窗宽度。因此,不同心脏相位的心脏门限函数的非零加权值重叠。
心脏门限函数重叠的另一个原因是选择心脏相位,使得它们的距离与测量的时间分辨率对应并且用于整个心动周期的门限函数非零。
根据本发明的典型实施例,该方法从步骤1开始,将要重建的感兴趣的心脏体积f(x,y,z)定义为大小为(nx,ny,nz,nt)的多相位体积f(x,y,z,t)。此处,nx(ny,nz)表示在x(y,z)大小中的体素的数量,且nt表示相位的数量。
然后,在步骤2,将用于有序子集迭代重建的子集大小定义为用于非门限迭代重建的子集大小乘以平均心动周期长度与平均门限窗口宽度的比值。
然后,在步骤3,作为投影数量的函数,根据ECG值、扫描参数及视场大小和位置计算用于每个心脏相位的心脏加权函数。
之后,在步骤4,从多相位数据集选择那些体积,其对于这种投影具有大于零的心脏权数并定义了相应的相位相关投影的集合。
这样的一组体积201,202,203,204在图3中示出,其示出了用于单个体素及所选的投影的投影几何结构的示意性表示。每个体积201至204在不同的相点具有相同的空间体积。而且,如投影205,206,207,208所示,每个相位体积201至204与同样的投影205至208相关。
然后,在步骤5,计算出用于将在这些相位体积201,202,203,204中的位置(x,y,z)处的体素值映射到投影上的系数。
然后,在步骤6,应用这些计算出的系数将体素从全部不同的心脏相位体积201至204正向投影到定义为相位阵列的投影上。
重复步骤5和6,直到对于所考虑的全部相位,相位体积201到204的全部体素都被投影到特定投影上。
而且,对子集中的每个投影重复包括步骤4到7的循环。
然后,在步骤8,对这个子集的全部投影,计算相位相关投影的全部堆栈的不同投影并且乘以心脏权数。
对于子集中的每个投影执行后面的步骤9到11。
在步骤9,从这个投影计算更新在体积中的单个体素所需的系数。
在步骤10,使用这些系数更新在全部体积中的体素,该体素对于这个投影具有非零心脏权数。
然后,在步骤11,对每个体素重复步骤9和10,直至整个体积被更新。方法跳至步骤12。
在步骤13,对全部子集重复方法步骤4至12。
然后,方法跳至步骤14,在那里结束。
图4示出了心脏加权函数301的示意性表示。垂线302、303、304、305表示与图3中所述的四个相位体积对应的不同相点。如可从心脏加权函数301知道的,每个相位体积对应于不同的心脏权数。
图5示出了根据本发明的用于执行根据本发明方法典型实施例的图像处理设备的典型实施例。图5中的图像处理设备400包括中央处理单元(CPU)或图像处理器401,其连接到用于存储描述感兴趣对象(例如胸部或其它组织部分)的图像的存储器402。数据处理器401连接到例如CT设备的多个诊断设备的输入/输出网络。而且,数据处理器401连接到显示设备403,例如计算机监视器,用于显示数据处理器410中计算的或改编的信息或图像。操作者或用户通过键盘404和/或其它输出设备,其未在图5中示出,与数据处理器401交互。
而且,通过总线系统405,还可能把图像处理和控制处理器401连接到例如运动监视器,其监视感兴趣对象的运动。例如,运动传感器可以是呼气传感器或心电图单元。
本发明的典型实施例可以作为成像工作站的软件选项出售。
根据本发明典型实施例的检查装置、图像处理设备、方法、程序单元及计算机可读介质提供了改善的图像质量及更高速度的迭代心脏锥束CT重建。例如,重建50个心脏相位,并且为20%心动周期的平均门限窗口宽度再次导致10倍的重建速度。
应当注意到,术语“包括”不排除其它元件或步骤,并且“一个”或“一”不排除多个,并且单个处理器或系统可实现权利要求所述的多个装置或单元的功能。也可组合结合不同实施例描述的元件。
还应知道,权利要求中的参考标号不构成对权利要求范围的限制。
Claims (14)
1.一种用于检查感兴趣对象(107)的检查装置(100),所述检查装置(100)包括:
采集单元(108),用于采集所述感兴趣对象(107)的四维数据集;
计算单元(118),适于:
在步骤a)从所述四维数据集选择第一相位体积和第二相位体积;
在步骤b)计算用于将所述第一相位体积中的第一体素映射到投影上的系数;
然后,在步骤c)基于用于映射的所述系数将所述第一相位体积的第一体素和所述第二相位体积的第二体素正向投影;
重复步骤b)和c)直到所述第一相位体积和所述第二相位体积的所有体素都被正向投影;
然后,在步骤d)计算用于更新所述正向投影的第一体素的系数;
在步骤e)基于用于更新的所述系数更新所述正向投影的第一体素和所述正向投影的第二体素;
针对每个体素重复步骤d)和e)直到所述第一相位体积和所述第二相位体积被更新为止;
其中,所述第一相位体积和所述第二相位体积具有相同的空间体积但是处于不同的相点,并且是要被重建的体积;并且
其中,所述第一相位体积中的所述第一体素和所述第二相位体积中的所述第二体素具有相同的空间坐标但是不同的相点。
2.根据权利要求1所述的检查装置(100),
其中,所述计算单元(118)还适于:
计算加权函数,产生对应于所述第一相位体积的第一权数和对应于所述第二相位体积的第二权数。
3.根据权利要求1所述的检查装置(100),
其中,所述计算单元(118)还适于:
基于用于所述第一体素的投影和用于所述第二体素的投影计算差分投影;和
基于所述差分投影计算用于更新的所述系数。
4.根据权利要求3所述的检查装置(100),
其中,用对应的心脏权数乘以所述差分投影。
5.根据权利要求2所述的检查装置(100),还包括:
心电图单元,适于采集心电图数据;
其中,所述加权函数为用于所述第一相位和所述第二相位的心脏加权函数。
6.根据权利要求5所述的检查装置(100),
其中,所述加权函数为至少一个投影数量的函数,所述投影数量来自所述心电图数据、扫描参数、视场大小、和视场位置。
7.根据权利要求2所述的检查装置(100),
其中,所述第一权数和所述第二权数大于零。
8.根据权利要求1所述的检查装置(100),还包括:
电磁辐射源(104),适于向所述感兴趣对象(107)发射电磁辐射;和
准直器(105),布置在所述电磁辐射源(104)与所述采集单元(108)之间;
其中,所述准直器(105)适于使所述电磁辐射源(104)发射的电磁辐射束准直以形成扇束或锥束。
9.根据权利要求1所述的检查装置(100),
其中,所述检查装置(100)适于作为心脏锥束计算机断层摄影装置。
10.根据权利要求1所述的检查装置(100),
其中,所述四维数据集包括第一子集和第二子集;
其中,所述第一和第二子集中每个的大小都与用于非门限迭代重建的子集大小乘以平均心动周期长度并除以平均门限窗口宽度相对应。
11.根据权利要求1所述的检查装置(100),其中,所述检查装置(100)被配置作为行李检验装置、医疗应用装置、材料测试装置和材料科学分析装置中的一个。
12.一种用于检查感兴趣对象的图像处理设备,所述图像处理设备包括:
存储器,用于存储所述感兴趣对象(107)的四维数据集;
计算单元(108)适于:
在步骤a)从所述四维数据集选择第一相位体积和第二相位体积;
在步骤b)计算用于将所述第一相位体积中的第一体素映射到投影上的系数;
然后,在步骤c)基于用于映射的所述系数将所述第一相位体积的第一体素和所述第二相位体积的第二体素正向投影;
重复步骤b)和c)直到所述第一相位体积和所述第二相位体积的所有体素都被正向投影;
然后,在步骤d)计算用于更新所述正向投影的第一体素的系数;
在步骤e)基于用于更新的所述系数更新所述正向投影的第一体素和所述正向投影的第二体素;
针对每个体素重复步骤d)和e)直到所述第一相位体积和所述第二相位体积被更新为止;
其中,所述第一相位体积和所述第二相位体积具有相同的空间体积但是处于不同的相点,并且是要被重建的体积;并且
其中,所述第一相位体积中的所述第一体素和所述第二相位体积中的所述第二体素具有相同的空间坐标但是不同的相点。
13.一种用于检查感兴趣对象(107)的方法,该方法如下步骤:
采集所述感兴趣对象(107)的四维数据集;
在步骤a)从所述四维数据集选择第一相位体积和第二相位体积;
在步骤b)计算用于将所述第一相位体积中的第一体素映射到投影上的系数;
然后,在步骤c)基于用于映射的所述系数将所述第一相位体积的第一体素和所述第二相位体积的第二体素正向投影;
重复步骤b)和c)直到所述第一相位体积和所述第二相位体积的所有体素都被正向投影;
然后,在步骤d)计算用于更新所述正向投影的第一体素的系数;
在步骤e)基于用于更新的所述系数更新所述正向投影的第一体素和所述正向投影的第二体素;
针对每个体素重复步骤d)和e)直到所述第一相位体积和所述第二相位体积被更新为止;
其中,所述第一相位体积和所述第二相位体积具有相同的空间体积但是处于不同的相点,并且是要被重建的体积;并且
其中,所述第一相位体积中的所述第一体素和所述第二相位体积中的所述第二体素具有相同的空间坐标但是不同的相点。
14.一种用于检查感兴趣对象(107)的装置,该装置包括:
用于采集所述感兴趣对象(107)的四维数据集的单元;
用于在步骤a)从所述四维数据集选择第一相位体积和第二相位体积的单元;
用于在步骤b)计算用于将所述第一相位体积中的第一体素映射到投影上的系数的单元;
用于然后,在步骤c)基于用于映射的所述系数将所述第一相位体积的第一体素和所述第二相位体积的第二体素正向投影的单元;
用于重复步骤b)和c)直到所述第一相位体积和所述第二相位体积的所有体素都被正向投影的单元;
用于然后,在步骤d)计算用于更新所述正向投影的第一体素的系数的单元;
用于在步骤e)基于用于更新的所述系数更新所述正向投影的第一体素和所述正向投影的第二体素的单元;
用于针对每个体素重复步骤d)和e)直到所述第一相位体积和所述第二相位体积被更新为止的单元;
其中,所述第一相位体积和所述第二相位体积具有相同的空间体积但是处于不同的相点,并且是要被重建的体积;并且
其中,所述第一相位体积中的所述第一体素和所述第二相位体积中的所述第二体素具有相同的空间坐标但是不同的相点。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP05108408.5 | 2005-09-13 | ||
EP05108408 | 2005-09-13 | ||
CN200680033400.6A CN101262821A (zh) | 2005-09-13 | 2006-09-01 | 高效迭代四维心脏锥束ct重建 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200680033400.6A Division CN101262821A (zh) | 2005-09-13 | 2006-09-01 | 高效迭代四维心脏锥束ct重建 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103750854A CN103750854A (zh) | 2014-04-30 |
CN103750854B true CN103750854B (zh) | 2016-03-23 |
Family
ID=37865323
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310738121.1A Expired - Fee Related CN103750854B (zh) | 2005-09-13 | 2006-09-01 | 高效迭代四维心脏锥束ct重建 |
CN200680033400.6A Pending CN101262821A (zh) | 2005-09-13 | 2006-09-01 | 高效迭代四维心脏锥束ct重建 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200680033400.6A Pending CN101262821A (zh) | 2005-09-13 | 2006-09-01 | 高效迭代四维心脏锥束ct重建 |
Country Status (6)
Country | Link |
---|---|
US (1) | US7676018B2 (zh) |
EP (1) | EP1933708B1 (zh) |
CN (2) | CN103750854B (zh) |
AT (1) | ATE493072T1 (zh) |
DE (1) | DE602006019292D1 (zh) |
WO (1) | WO2007031899A2 (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007072279A1 (en) * | 2005-12-19 | 2007-06-28 | Philips Intellectual Property & Standards Gmbh | Iterative image reconstruction of a moving object from projection data |
ITVA20070065A1 (it) * | 2007-07-25 | 2009-01-26 | St Microelectronics Srl | Metodo e dispositivo di pilotaggio per motore brushless con profilo di tensione predisposto per una commutazione progressiva e automatica da un pilotaggio di tipo sinusoidale trifase ad un pliotaggio trifase ad onda quadra |
US8787521B2 (en) | 2009-12-23 | 2014-07-22 | General Electric Company | System and method of iterative image reconstruction for computed tomography |
CN103403702A (zh) | 2010-12-24 | 2013-11-20 | 澳大利亚国立大学 | 动态多维图像数据的重构 |
CN103236058B (zh) * | 2013-04-25 | 2016-04-13 | 内蒙古科技大学 | 获取四维心脏图像感兴趣体积的方法 |
US9626779B1 (en) * | 2015-10-20 | 2017-04-18 | International Business Machines Corporation | Efficient back-projection operation using precomputed table |
US11200709B2 (en) * | 2016-12-27 | 2021-12-14 | Canon Medical Systems Corporation | Radiation image diagnostic apparatus and medical image processing apparatus |
CN108182720A (zh) * | 2017-12-28 | 2018-06-19 | 上海联影医疗科技有限公司 | 心脏图像的重建方法及医疗影像系统 |
US11138770B2 (en) | 2017-11-06 | 2021-10-05 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for medical imaging |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7058440B2 (en) | 2001-06-28 | 2006-06-06 | Koninklijke Philips Electronics N.V. | Dynamic computed tomography imaging using positional state modeling |
US7020511B2 (en) * | 2002-02-22 | 2006-03-28 | Ge Medical Systems Global Technology Company, Llc | Method for three dimensional cine EBA/CTA imaging |
US6904118B2 (en) * | 2002-07-23 | 2005-06-07 | General Electric Company | Method and apparatus for generating a density map using dual-energy CT |
US7042975B2 (en) * | 2002-10-25 | 2006-05-09 | Koninklijke Philips Electronics N.V. | Four-dimensional helical tomographic scanner |
US7545903B2 (en) | 2003-07-16 | 2009-06-09 | Koninklijke Philips Electronics N.V. | Reconstruction of an image of a moving object from volumetric data |
DE10354214A1 (de) | 2003-11-20 | 2005-06-02 | Siemens Ag | Verfahren zur Erzeugung von tomographischen Schnittbildern eines sich periodisch bewegenden Objektes mit mehreren Fokus-Detektor-Kombinationen |
US20060198491A1 (en) * | 2005-03-04 | 2006-09-07 | Kabushiki Kaisha Toshiba | Volumetric computed tomography system for imaging |
-
2006
- 2006-09-01 EP EP06795860A patent/EP1933708B1/en not_active Not-in-force
- 2006-09-01 CN CN201310738121.1A patent/CN103750854B/zh not_active Expired - Fee Related
- 2006-09-01 US US12/065,615 patent/US7676018B2/en not_active Expired - Fee Related
- 2006-09-01 AT AT06795860T patent/ATE493072T1/de not_active IP Right Cessation
- 2006-09-01 WO PCT/IB2006/053055 patent/WO2007031899A2/en active Application Filing
- 2006-09-01 CN CN200680033400.6A patent/CN101262821A/zh active Pending
- 2006-09-01 DE DE602006019292T patent/DE602006019292D1/de active Active
Non-Patent Citations (2)
Title |
---|
Transmission maximum-likelihood reconstruction with ordered subsets for cone beam CT;S H Manglos et al;《PHYSICS IN MEDICINE AND BIOLOGY》;19951231;第40卷;第1225-1241页 * |
Truncation artifact suppression in cone-beam radionulide transmission CT using maximum likelihood techniques: evaluation with human subjects;Manglos S H et al;《PHYSICS IN MEDICINE AND BIOLOGY》;19920301;第47卷(第3期);第550页 * |
Also Published As
Publication number | Publication date |
---|---|
DE602006019292D1 (de) | 2011-02-10 |
US20080267342A1 (en) | 2008-10-30 |
ATE493072T1 (de) | 2011-01-15 |
EP1933708B1 (en) | 2010-12-29 |
EP1933708A2 (en) | 2008-06-25 |
CN103750854A (zh) | 2014-04-30 |
WO2007031899A2 (en) | 2007-03-22 |
CN101262821A (zh) | 2008-09-10 |
WO2007031899A3 (en) | 2008-04-17 |
US7676018B2 (en) | 2010-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103750854B (zh) | 高效迭代四维心脏锥束ct重建 | |
CN101313334B (zh) | 高对比度对象的运动补偿ct重建 | |
CN1926578B (zh) | 运动补偿 | |
CN101262819B (zh) | Ct散射的直接测量及校正 | |
US6665370B2 (en) | Computed tomography method and apparatus for acquiring images dependent on a time curve of a periodic motion of the subject | |
US7221728B2 (en) | Method and apparatus for correcting motion in image reconstruction | |
US20080267455A1 (en) | Method for Movement Compensation of Image Data | |
US20070248253A1 (en) | Artifact Reduction | |
CN101505660A (zh) | 多管x射线探测 | |
EP1814459A1 (en) | Computer tomography apparatus and method for examining an object of interest | |
JP2007512034A (ja) | 発散ビームスキャナのための画像再構成方法 | |
KR20190038574A (ko) | 시공간적으로 오버랩하는 x선들로부터 3d 영상을 재구성하는 방법과 장치 | |
CN100573580C (zh) | 多焦点采集 | |
US20080226016A1 (en) | Reconstruction Algorithm for Object Point Outside the Scan-Field-of-View | |
US8135197B2 (en) | Automatic adaptive soft tissue thresholding for two-pass CT cone-beam artifact reduction | |
JP2010504160A (ja) | ヘリカル・コーン・ビームデータの射線整合方式による再構成 | |
CN109416833A (zh) | 根据谱ct数据确定钙含量 | |
US8854355B2 (en) | System and method of visualizing features in an image | |
CN101478920A (zh) | 对狭窄的局部运动补偿重建 | |
EP1694211B1 (en) | Computer tomography method for periodically moving objects | |
JP2009537814A (ja) | CSCTのq範囲を拡大する装置及び方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160323 Termination date: 20200901 |
|
CF01 | Termination of patent right due to non-payment of annual fee |