CN110473269A - 一种图像重建方法、***、设备及存储介质 - Google Patents
一种图像重建方法、***、设备及存储介质 Download PDFInfo
- Publication number
- CN110473269A CN110473269A CN201910730387.9A CN201910730387A CN110473269A CN 110473269 A CN110473269 A CN 110473269A CN 201910730387 A CN201910730387 A CN 201910730387A CN 110473269 A CN110473269 A CN 110473269A
- Authority
- CN
- China
- Prior art keywords
- data
- scanning data
- raw scanning
- orthographic projection
- correction
- 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.)
- Granted
Links
Classifications
-
- 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/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明实施例公开了一种图像重建方法、***、设备及存储介质,通过获取目标对象的原始扫描数据,根据原始扫描数据确定正投影数据,计算原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正,对校正后的原始扫描数据进行重建,由于校正运动场对原始扫描数据进行了校正,因而得到伪影更少的目标重建图像,如果不满足要求,将目标重建图像作为初始图像,重新对初始图像进行正投影以及计算校正运动场以及对原始扫描数据进行校正,直至获得较为理想的目标重建图像,解决了现有技术中在图像域上对初始重建图像校正效果较差的问题,实现了在投影域上对原始扫描数据进行校正以减少重建图像的运动伪影的效果。
Description
技术领域
本发明实施例涉及医疗器械技术,尤其涉及一种图像重建方法、***、设备及存储介质。
背景技术
一般情况下,采用医学影像设备(如CT(Computed Tomography,电子计算机断层扫描)、PET(Positron Emission Computed Tomography,正电子发射型计算机断层显像))对被检测者的某一扫描区域扫描过程中,被检测者可能存在自主或不自主的运动(如被检测者自主的呼吸运动、不自主的心脏搏动和胃肠蠕动等),这些自主或不自主的运动就会在重建图像上形成运动伪影,降低图像质量甚至影响诊断。
为此,现有技术常采用运动校正方法消除形成的运动伪影,如选取初始重建图像的不同范围的投影数据,或者选取数据范围相同但权重不同的投影数据,将这些投影数据进行重建得到多个投影数据的重建图像,根据多个投影数据的重建图像计算运动场,并根据运动场对多个投影数据的重建图像进行校正。
可见,上述方法是在图像域上进行运动校正,在具体实现时,如果初始重建图像的运动伪影较明显,则根据选择的多个投影数据的重建图像计算得到的运动场很不准确,导致重建图像的校正效果较差。
发明内容
本发明实施例提供了一种图像重建方法、***、设备及存储介质,以实现在投影域进行伪影校正,大大减少重建图像的运动伪影。
第一方面,本发明实施例提供了一种图像重建方法,包括:
获取目标对象的原始扫描数据,根据所述原始扫描数据确定与所述原始扫描数据对应的正投影数据;
计算所述原始扫描数据与所述正投影数据之间的校正运动场,根据所述校正运动场对所述原始扫描数据进行校正;
对校正后的原始扫描数据进行重建,得到目标重建图像。
第二方面,本发明实施例还提供了一种图像重建***,该***包括:
获取模块,用于获取目标对象的原始扫描数据,根据所述原始扫描数据确定与所述原始扫描数据对应的正投影数据;
计算模块,用于计算所述原始扫描数据与所述正投影数据之间的校正运动场,根据所述校正运动场对所述原始扫描数据进行校正;
重建模块,用于对校正后的原始扫描数据进行重建,得到目标重建图像。
第三方面,本发明实施例还提供了一种图像重建设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其中,所述处理器执行所述计算机程序时实现如第一方面任一项所述的图像重建方法。
第四方面,本发明实施例还提供了一种包含计算机可执行指令的存储介质,其中,所述计算机可执行指令在由计算机处理器执行时实现如第一方面任一项所述的图像重建方法。
本发明实施例的技术方案,通过获取目标对象的原始扫描数据,根据原始扫描数据确定与原始扫描数据对应的正投影数据,并计算原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正,然后对校正后的原始扫描数据进行重建,得到目标重建图像,由于校正运动场对原始扫描数据进行了校正,因而校正后的原始扫描数据更加准确,并得到伪影更少的目标重建图像,解决了现有技术中在图像域上对初始重建图像校正效果较差的问题,实现了在投影域上对原始扫描数据进行校正以减少重建图像的运动伪影的效果。
附图说明
图1为本发明实施例一提供的一种图像重建方法的流程示意图;
图2为本发明实施例二提供的一种图像重建方法的流程示意图;
图3为本发明实施例三提供的一种图像重建方法的流程示意图;
图4为本发明实施例三提供的一种图像重建方法的结构示意图;
图5为本发明实施例四提供的一种图像重建***的结构示意图;
图6为本发明实施例五提供的一种图像重建设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
图1为本发明实施例一提供的一种图像重建方法的流程示意图,本实施例可适用于在投影域对获取的原始扫描图像进行伪影校正的情况,该方法可以由图像重建***来执行,其中该***可由软件和/或硬件实现,并一般集成在终端中。具体参见图1所示,该方法可以包括如下步骤:
S110,获取目标对象的原始扫描数据,根据原始扫描数据确定与原始扫描数据对应的正投影数据。
其中,目标对象可以理解为被扫描对象,示例性地,被扫描对象可以人体、动物体或者一些仿真实验体等。目标对象也可以是指被扫描个体的局部待扫描区域。例如,可以是被扫描的人体的头部区域或胸部区域等。
其中,原始扫描数据可以为CT扫描数据、PET扫描数据以及核磁共振扫描数据等,原始扫描数据可以是运动伪影较强的扫描数据段,也可以是运动伪影较弱的扫描数据段,本实施例不做具体限定。
可选地,可以通过确定与原始扫描数据对应的初始图像,对初始图像进行正投影,得到与原始扫描数据对应的正投影数据。
具体地,对原始扫描数据进行重建,得到与原始扫描数据对应的初始图像;或者,对目标对象的历史扫描数据进行重建,得到与原始扫描数据对应的初始图像。即,初始图像可以是目标对象当前的图像,也可以是目标对象先验的图像。
示例性地,如果初始图像是目标对象当前的图像,对当前的原始扫描数据进行重建,便得到该初始图像。类似地,如果初始图像是目标对象先验的图像,对先验的历史扫描数据进行重建,便得到该初始图像。可以理解的是,可以对原始扫描数据进行一次重建或者若干次重建得到目标对象当前的图像,类似地,也可以对先验的历史扫描数据通过一次重建或者若干次重建得到目标对象先验的图像。
S120,计算原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正。
其中,原始扫描数据可以包括数据值、采集位置以及角度等。相应的,与原始扫描数据对应的正投影数据也可以包括数据值、采集位置以及角度等。因而,原始扫描数据与正投影数据之间的校正运动场可以包括原始扫描数据的与正投影数据之间的相对运动方向以及相对距离等。
可以理解的是,原始扫描数据中可能包括较多的伪影数据,而相对于原始扫描数据,与原始扫描数据对应的正投影数据往往可能包括较少的伪影数据,因此,可以根据原始扫描数据与正投影数据之间的校正运动场,对原始扫描数据的校正。
S130,对校正后的原始扫描数据进行重建,得到目标重建图像。
由上可知,通过原始扫描数据与正投影数据就可得到校正后的原始扫描数据,其中,校正后的原始扫描数据为运动伪影较少的投影数据,相应的,得到的目标重建图像也为运动伪影较少的重建图像。
可选地,采用滤波反投影方法对校正后的原始扫描数据进行重建,得到目标重建图像;也可以采用迭代重建方法对校正后的原始扫描数据进行重建,得到目标重建图像。
本发明实施例提供了一种图像重建方法,通过获取目标对象的原始扫描数据,根据原始扫描数据确定与原始扫描数据对应的正投影数据,并计算原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正,然后对校正后的原始扫描数据进行重建,得到目标重建图像,由于校正运动场对原始扫描数据进行了校正,因而校正后的原始扫描数据更加准确,并得到伪影更少的目标重建图像,解决了现有技术中在图像域上对初始重建图像校正效果较差的问题,实现了在投影域上对原始扫描数据进行校正以减少重建图像的运动伪影的效果。
实施例二
图2为本发明实施例二提供的一种图像重建方法的流程示意图。本实施例的技术方案在上述实施例的基础上进行了细化,可选地,所述根据所述原始扫描数据确定与所述原始扫描数据对应的正投影数据,包括:确定与所述原始扫描数据对应的初始图像,对所述初始图像进行正投影,得到与所述原始扫描数据对应的正投影数据。
在上述各技术方案的基础上,可选地,所述计算所述原始扫描数据与所述正投影数据之间的校正运动场,还可包括:对所述原始扫描数据与所述正投影数据进行配准,得到所述原始扫描数据与所述正投影数据之间的校正运动场。
具体参见图2所示,本实施例的方法可以包括如下步骤:
S210,获取目标对象的原始扫描数据,根据原始扫描数据确定与原始扫描数据对应的正投影数据。
可选地,获取目标对象的原始扫描数据之后,对原始扫描数据进行重建得到初始图像,然后对初始图像进行正投影,得到与原始扫描数据对应的正投影数据。
可选地,正投影数据可以通过以下方式获取:遍历初始图像的每个体素;针对每个体素的至少一个角度做正投影,得到与原始扫描数据对应的正投影数据。
可选地,如果初始图像对应的是CT图像,可以将该CT图像离散化成规则的网格,每个网格单元可以作为一个体素。
示例性地,如果正投影的角度为ω个,体素个数为N*N*N个,如果在每个角度下做正投影可以得到M*M个数据,则对每个体素的至少一个角度做正投影可以得到[M*M-M*M*ω]个投影数据,由此,对N*N*N个体素的至少一个角度做正投影就可得到大量的正投影数据,由此可以得到与原始扫描数据对应的正投影数据。
S220,对原始扫描数据与正投影数据进行配准,得到原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正。
其中,对原始扫描数据与正投影数据进行配准的配准方法可以是刚性配准方法也可以是非刚性配准方法。
其中,计算校正运动场的原始扫描数据可以是当前的扫描数据,也可以是对当前的扫描数据进行插值处理或者重排处理之后的扫描数据。
其中,校正运动场可以是连续区域也可以是不连续区域。
可选地,基于将原始扫描数据的所有体素将原始扫描数据与正投影数据进行配准,得到原始扫描数据与正投影数据的之间的校正运动场;或者,基于原始扫描数据的部分体素将原始扫描数据与正投影数据进行配准,得到原始扫描数据与正投影数据之间的初始运动场,对初始运动场进行插值,得到原始扫描数据与正投影数据之间的校正运动场。
可以理解的是,可以根据计算的校正运动场对原始扫描数据进行修改,以使修改后的原始扫描数据更加准确,再根据校正运动场对修改后的原始扫描数据进行校正,或者,可以给原始扫描数据加上一个修正项,以使修改后的原始扫描数据更加准确,再根据校正运动场对加上修正项的原始扫描数据进行校正。
S230,对校正后的原始扫描数据进行重建,得到目标重建图像。
本发明实施例提供了一种图像重建方法,通过获取目标对象的原始扫描数据,确定与原始扫描数据对应的初始图像,对初始图像进行正投影,得到与原始扫描数据对应的正投影数据,并对原始扫描数据与正投影数据进行配准,得到原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正,解决了现有技术中在图像域上对初始重建图像校正效果较差的问题,实现通过对原始扫描数据与正投影数据进行配准以对原始扫描数据进行校正,并根据校正后的原始扫描数据进行重建以得到伪影更少的目标重建图像的效果。
实施例三
图3为本发明实施例三提供的一种图像重建方法的流程示意图。本实施例的技术方案在上述实施例的基础上进行了细化,可选地,所述根据所述原始扫描数据确定与所述原始扫描数据对应的正投影数据,包括:确定与所述原始扫描数据对应的初始图像,对所述初始图像进行正投影,得到与所述原始扫描数据对应的正投影数据。
在上述技术方案的基础上,可选地,所述对校正后的原始扫描数据进行重建,得到目标重建图像,包括:对校正后的原始扫描数据进行重建,得到当前重建图像;判断所述当前重建图像是否满足预设收敛条件;如果否,将所述当前重建图像作为初始图像,重复执行对所述初始图像进行正投影,得到与所述原始扫描数据对应的正投影数据的步骤,直至所述当前重建图像满足预设收敛条件;如果是,则将所述当前重建图像作为目标重建图像。
可以理解的是,如果得到的目标对象的运动较强,通过上述实施例得到的目标重建图像可能不能满足要求,还可以将目标重建图像作为当前重建图像,重新执行上述实施例的步骤。具体参见图3所示,本实施例的方法可以包括如下步骤:
S310,获取目标对象的原始扫描数据。
S320,确定与原始扫描数据对应的初始图像。
S330,对初始图像进行正投影,得到与原始扫描数据对应的正投影数据。
S340,计算原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正。
S350,对校正后的原始扫描数据进行重建,得到当前重建图像。
S360,判断当前重建图像是否满足预设收敛条件,如果是,执行S370,如果否,执行S380。
其中,预设收敛条件可以是预设的收敛阈值,如果当前重建图像更新幅度小于预设的收敛阈值,则说明当前重建图像满足预设收敛条件;如果当前重建图像的更新幅度与大于预设的收敛阈值,则说明当前重建图像不满足预设收敛条件。
S370,将当前重建图像作为目标重建图像。
S380,将当前重建图像作为初始图像,返回继续执行S330。
可选的,如果当前重建图像更新幅度大于收敛阈值,将当前重建图像作为初始图像,并返回继续执行S330以进行迭代,重新计算原始扫描数据与正投影数据之间的校正运动场,并对原始扫描数据进行校正并且进行图像重建的步骤,直至当前重建图像满足预设收敛条件。
示例性地解释上述方法,如图4所示,获取到原始扫描数据之后,可以对原始扫描数据进行重建,得到初始图像,然后采用上述实施例提供的正投影方法对得到的初始图像进行正投影,得到正投影数据,接着,可以采用上述实施例提供的配准方法对原始扫描数据与正投影数据进行配准,得到校正运动场,并根据校正运动场对原始扫描数据进行校正,并获得运动校正后的扫描数据,然后对运动校正后的扫描数据进行重建得到目标重建图像。可以理解的是,如果上述目标重建图像不满足预设收敛条件,可以将重建得到的目标重建图像作为初始图像,重复执行对初始图像进行正投影,得到与原始扫描数据对应的正投影数据的步骤,直至当前重建图像满足预设收敛条件。
本发明实施例公开了一种图像重建方法,通过获取目标对象的原始扫描数据,根据原始扫描数据确定与原始扫描数据对应的正投影数据,计算原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正,对校正后的原始扫描数据进行重建,得到目标重建图像,如果不满足要求,将当前的目标重建图像作为初始图像,重复执行对初始图像进行正投影,得到与原始扫描数据对应的正投影数据的步骤,直至当前重建图像满足预设收敛条件,通过以上步骤,可以得到较为理想的目标重建图像,解决了现有技术中在图像域上对初始重建图像校正效果较差的问题,实现了在投影域上对原始扫描数据进行校正以减少重建图像的运动伪影的效果。
实施例四
图5为本发明实施例四提供的一种图像重建***的结构示意图。参见图5所示,该***包括:获取模块51、计算模块52以及重建模块53。
其中,获取模块51,用于获取目标对象的原始扫描数据,根据原始扫描数据确定与原始扫描数据对应的正投影数据;计算模块52,用于计算原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正;重建模块53,用于对校正后的原始扫描数据进行重建,得到目标重建图像。
在上述各技术方案的基础上,获取模块51,还用于确定与原始扫描数据对应的初始图像,对初始图像进行正投影,得到与原始扫描数据对应的正投影数据。
在上述各技术方案的基础上,获取模块51,还用于:
对原始扫描数据进行重建,得到与原始扫描数据对应的初始图像;或者,
对目标对象的历史扫描数据进行重建,得到与原始扫描数据对应的初始图像。
在上述各技术方案的基础上,获取模块51,还用于:
遍历初始图像的每个体素;
针对每个体素的至少一个角度做正投影,得到与原始扫描数据对应的正投影数据。
在上述各技术方案的基础上,计算模块52还用于:
对原始扫描数据与正投影数据进行配准,得到原始扫描数据与正投影数据之间的校正运动场。
在上述各技术方案的基础上,计算模块52还用于:
基于将原始扫描数据的所有体素将原始扫描数据与正投影数据进行配准,得到原始扫描数据与正投影数据的之间的校正运动场;
或者,
基于原始扫描数据的部分体素将原始扫描数据与正投影数据进行配准,得到原始扫描数据与正投影数据之间的初始运动场,
对初始运动场进行插值,得到原始扫描数据与正投影数据之间的校正运动场。
在上述各技术方案的基础上,重建模块53还用于:
对校正后的原始扫描数据进行重建,得到当前重建图像;
判断当前重建图像是否满足预设收敛条件;如果否,将当前重建图像作为初始图像,重复执行对初始图像进行正投影,得到与原始扫描数据对应的正投影数据的步骤,直至当前重建图像满足预设收敛条件;
如果是,则将当前重建图像作为目标重建图像。
本发明实施例提供了一种图像重建***,通过获取目标对象的原始扫描数据,根据原始扫描数据确定与原始扫描数据对应的正投影数据,并计算原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正,然后对校正后的原始扫描数据进行重建,得到目标重建图像,由于校正运动场对原始扫描数据进行了校正,因而校正后的原始扫描数据更加准确,并得到伪影更少的目标重建图像,解决了现有技术中在图像域上对初始重建图像校正效果较差的问题,实现了在投影域上对原始扫描数据进行校正以减少重建图像的运动伪影的效果。
实施例五
图6为本发明实施例五提供的一种图像重建设备的结构示意图。图6示出了适于用来实现本发明实施方式的示例性图像重建设备12的框图。图6显示的图像重建设备12仅仅是一个示例,不应对本发明实施例的功能和使用范围带来任何限制。
如图6所示,图像重建设备12以通用计算设备的形式表现。图像重建设备12的组件可以包括但不限于:一个或者多个处理器或者处理单元16,***存储器28,连接不同***组件(包括***存储器28和处理单元16)的总线18。
总线18表示几类总线结构中的一种或多种,包括存储器总线或者存储器控制器,***总线,图形加速端口,处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说,这些体系结构包括但不限于工业标准体系结构(ISA)总线,微通道体系结构(MAC)总线,增强型ISA总线、视频电子标准协会(VESA)局域总线以及***组件互连(PCI)总线。
图像重建设备12典型地包括多种计算机***可读介质。这些介质可以是任何能够被图像重建设备12访问的可用介质,包括易失性和非易失性介质,可移动的和不可移动的介质。
***存储器28可以包括易失性存储器形式的计算机***可读介质,例如随机存取存储器(RAM)30和/或高速缓存存储器32。图像重建设备12可以进一步包括其它可移动/不可移动的、易失性/非易失性计算机***存储介质。仅作为举例,存储***34可以用于读写不可移动的、非易失性磁介质(图6未显示,通常称为“硬盘驱动器”)。尽管图6中未示出,可以提供用于对可移动非易失性磁盘(例如“软盘”)读写的磁盘驱动器,以及对可移动非易失性光盘(例如CD-ROM,DVD-ROM或者其它光介质)读写的光盘驱动器。在这些情况下,每个驱动器可以通过一个或者多个数据介质接口与总线18相连。存储器28可以包括至少一个程序产品,该程序产品具有一组(例如图像重建***的图像获取模块61和存储模块62)程序模块,这些程序模块被配置以执行本发明各实施例的功能。
具有一组(例如图像重建***的获取模块51、计算模块52和重建模块53)程序模块46的程序/实用工具44,可以存储在例如存储器28中,这样的程序模块46包括但不限于操作***、一个或者多个应用程序、其它程序模块以及程序数据,这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块46通常执行本发明所描述的实施例中的功能和/或方法。
图像重建设备12也可以与一个或多个外部设备14(例如键盘、指向设备、显示器24等)通信,还可与一个或者多个使得用户能与该图像重建设备12交互的设备通信,和/或与使得图像重建设备12能与一个或多个其它计算设备进行通信的任何设备(例如网卡,调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口22进行。并且,图像重建设备12还可以通过网络适配器20与一个或者多个网络(例如局域网(LAN),广域网(WAN)和/或公共网络,例如因特网)通信。如图所示,网络适配器20通过总线18与图像重建设备12的其它模块通信。应当明白,尽管图中未示出,可以结合图像重建设备12使用其它硬件和/或软件模块,包括但不限于:微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、RAID***、磁带驱动器以及数据备份存储***等。
处理单元16通过运行存储在***存储器28中的程序,从而执行各种功能应用以及数据处理,例如实现本发明实施例所提供的一种图像重建方法,该方法包括:
获取目标对象的原始扫描数据,根据原始扫描数据确定与原始扫描数据对应的正投影数据;
计算原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正;
对校正后的原始扫描数据进行重建,得到目标重建图像。
处理单元16通过运行存储在***存储器28中的程序,从而执行各种功能应用以及数据处理,例如实现本发明实施例所提供的一种图像重建方法。
当然,本领域技术人员可以理解,处理器还可以实现本发明任意实施例所提供的一种图像重建方法的技术方案。
实施例六
本发明实施例六还提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如本发明实施例所提供的一种图像重建方法,该方法包括:
获取目标对象的原始扫描数据,根据原始扫描数据确定与原始扫描数据对应的正投影数据;
计算原始扫描数据与正投影数据之间的校正运动场,根据校正运动场对原始扫描数据进行校正;
对校正后的原始扫描数据进行重建,得到目标重建图像。
当然,本发明实施例所提供的一种计算机可读存储介质,其上存储的计算机程序不限于如上的方法操作,还可以执行本发明任意实施例所提供的一种图像重建方法中的相关操作。
本发明实施例的计算机存储介质,可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的***、***或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行***、***或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在原始扫描数据、正投影数据、校正运动场等,其中承载了计算机可读的程序代码。这种传播的原始扫描数据、正投影数据、校正运动场等形式。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行***、***或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括——但不限于无线、电线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码,程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++,还包括常规的过程式程序设计语言—诸如”C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
值得注意的是,上述图像重建***的实施例中,所包括的各个模块只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。
Claims (10)
1.一种图像重建方法,其特征在于,包括:
获取目标对象的原始扫描数据,根据所述原始扫描数据确定与所述原始扫描数据对应的正投影数据;
计算所述原始扫描数据与所述正投影数据之间的校正运动场,根据所述校正运动场对所述原始扫描数据进行校正;
对校正后的原始扫描数据进行重建,得到目标重建图像。
2.根据权利要求1所述的方法,其特征在于,根据所述原始扫描数据确定与所述原始扫描数据对应的正投影数据,包括:
确定与所述原始扫描数据对应的初始图像,对所述初始图像进行正投影,得到与所述原始扫描数据对应的正投影数据。
3.根据权利要求2所述的方法,其特征在于,确定与所述原始扫描数据对应的初始图像,包括:
对所述原始扫描数据进行重建,得到与所述原始扫描数据对应的初始图像;或者,
对所述目标对象的历史扫描数据进行重建,得到与所述原始扫描数据对应的初始图像。
4.根据权利要求2所述的方法,其特征在于,对所述初始图像进行正投影,得到与所述原始扫描数据对应的正投影数据,包括:
遍历所述初始图像的每个体素;
针对所述每个体素的至少一个角度做正投影,得到与所述原始扫描数据对应的正投影数据。
5.根据权利要求1所述的方法,其特征在于,计算所述原始扫描数据与所述正投影数据之间的校正运动场,包括:
对所述原始扫描数据与所述正投影数据进行配准,得到所述原始扫描数据与所述正投影数据之间的校正运动场。
6.根据权利要求5所述的方法,其特征在于,对所述原始扫描数据与所述正投影数据进行配准,得到所述原始扫描数据与所述正投影数据的校正运动场,包括:
基于将所述原始扫描数据的所有体素将所述原始扫描数据与所述正投影数据进行配准,得到所述原始扫描数据与所述正投影数据的之间的校正运动场;
或者,
基于所述原始扫描数据的部分体素将所述原始扫描数据与所述正投影数据进行配准,得到所述原始扫描数据与所述正投影数据之间的初始运动场,对所述初始运动场进行插值,得到所述原始扫描数据与所述正投影数据之间的校正运动场。
7.根据权利要求2所述的方法,其特征在于,对校正后的原始数据进行重建,得到目标重建图像,包括:
对校正后的原始扫描数据进行重建,得到当前重建图像;
判断所述当前重建图像是否满足预设收敛条件;如果否,将所述当前重建图像作为初始图像,重复执行对所述初始图像进行正投影,得到与所述原始扫描数据对应的正投影数据的步骤,直至所述当前重建图像满足预设收敛条件;
如果是,则将所述当前重建图像作为目标重建图像。
8.一种图像重建***,其特征在于,包括:
获取模块,用于获取目标对象的原始扫描数据,根据所述原始扫描数据确定与所述原始扫描数据对应的正投影数据;
计算模块,用于计算所述原始扫描数据与所述正投影数据之间的校正运动场,根据所述校正运动场对所述原始扫描数据进行校正;
重建模块,用于对校正后的原始扫描数据进行重建,得到目标重建图像。
9.一种图像重建设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1-7中任一项所述的图像重建方法。
10.一种包含计算机可执行指令的存储介质,其特征在于,所述计算机可执行指令在由计算机处理器执行时实现如权利要求1-7中任一项所述的图像重建方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910730387.9A CN110473269B (zh) | 2019-08-08 | 2019-08-08 | 一种图像重建方法、***、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910730387.9A CN110473269B (zh) | 2019-08-08 | 2019-08-08 | 一种图像重建方法、***、设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110473269A true CN110473269A (zh) | 2019-11-19 |
CN110473269B CN110473269B (zh) | 2023-05-26 |
Family
ID=68510460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910730387.9A Active CN110473269B (zh) | 2019-08-08 | 2019-08-08 | 一种图像重建方法、***、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110473269B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111612867A (zh) * | 2020-06-01 | 2020-09-01 | 上海联影医疗科技有限公司 | 运动伪影校正方法、装置、计算机设备和可读存储介质 |
CN111798534A (zh) * | 2020-07-17 | 2020-10-20 | 沈阳先进医疗设备技术孵化中心有限公司 | 图像重建方法、装置、控制台设备及ct*** |
CN113313649A (zh) * | 2021-06-01 | 2021-08-27 | 上海联影医疗科技股份有限公司 | 图像重建方法及装置 |
CN113628293A (zh) * | 2020-05-06 | 2021-11-09 | 美国西门子医疗***股份有限公司 | 全身直接参数化图像重建中的帧间运动校正 |
WO2022022535A1 (en) * | 2020-07-28 | 2022-02-03 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for image processing |
Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6570951B1 (en) * | 2002-05-14 | 2003-05-27 | Ge Medical Systems Global Technology Company, Llc | Image space compensation scheme for reducing artifacts |
US6647134B1 (en) * | 2000-03-30 | 2003-11-11 | Mayo Foundation For Medical Education And Research | Autocorrection of MR projection images |
US20040136490A1 (en) * | 2002-07-23 | 2004-07-15 | Edic Peter Michael | Method and apparatus for correcting motion in image reconstruction |
CN101126722A (zh) * | 2007-09-30 | 2008-02-20 | 西北工业大学 | 基于配准模型仿真的锥束ct射束硬化校正方法 |
US20110123081A1 (en) * | 2009-11-25 | 2011-05-26 | David Sebok | Correcting and reconstructing x-ray images using patient motion vectors extracted from marker positions in x-ray images |
US20130051644A1 (en) * | 2011-08-29 | 2013-02-28 | General Electric Company | Method and apparatus for performing motion artifact reduction |
CN102982510A (zh) * | 2011-06-21 | 2013-03-20 | 通用电气公司 | Spect图像的运动校正 |
CN103455989A (zh) * | 2013-09-24 | 2013-12-18 | 南京大学 | 一种结合超声图像提高有限角度ct成像质量的方法 |
CN103679642A (zh) * | 2012-09-26 | 2014-03-26 | 上海联影医疗科技有限公司 | 一种ct图像金属伪影校正方法、装置及ct设备 |
US20140212011A1 (en) * | 2013-01-28 | 2014-07-31 | Koninklijke Philips Electronics N.V. | Spect motion-correction |
CN104424659A (zh) * | 2013-08-30 | 2015-03-18 | 西门子公司 | 图像数据的基于图像的运动补偿 |
CN104574292A (zh) * | 2014-11-26 | 2015-04-29 | 沈阳东软医疗***有限公司 | 一种ct图像的校正方法和装置 |
CN104867168A (zh) * | 2015-04-28 | 2015-08-26 | 南京邮电大学 | 基于p范数的压缩感知计算机断层扫描图像重建方法 |
US20160104279A1 (en) * | 2014-10-10 | 2016-04-14 | Cedars-Sinai Medical Center | Method and system for "push-button" comprehensive cardiac mr examination using continuous self-gated 3d radial imaging |
US20160110893A1 (en) * | 2014-10-21 | 2016-04-21 | Shenyang Neusoft Medical Systems Co., Ltd. | Correcting a ct scan image |
CN106683146A (zh) * | 2017-01-11 | 2017-05-17 | 上海联影医疗科技有限公司 | 一种图像重建方法和图像重建算法的参数确定方法 |
CN106683144A (zh) * | 2016-12-30 | 2017-05-17 | 上海联影医疗科技有限公司 | 一种图像迭代重建方法及装置 |
CN106960429A (zh) * | 2017-02-16 | 2017-07-18 | 中国科学院苏州生物医学工程技术研究所 | 一种ct图像金属伪影校正方法及装置 |
CN107845121A (zh) * | 2017-11-03 | 2018-03-27 | 中国工程物理研究院应用电子学研究所 | 一种探测器偏置扫描中加权伪影的校正方法 |
US20180137657A1 (en) * | 2016-11-16 | 2018-05-17 | University Of Central Florida Research Foundation, Inc. | System and method for motion estimation and compensation in helical computed tomography |
CN108876730A (zh) * | 2018-05-24 | 2018-11-23 | 沈阳东软医疗***有限公司 | 校正运动伪影的方法、装置及设备和存储介质 |
US20190180482A1 (en) * | 2016-07-04 | 2019-06-13 | Hitachi, Ltd. | Image reconstruction device, x-ray ct device, and image reconstruction method |
-
2019
- 2019-08-08 CN CN201910730387.9A patent/CN110473269B/zh active Active
Patent Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6647134B1 (en) * | 2000-03-30 | 2003-11-11 | Mayo Foundation For Medical Education And Research | Autocorrection of MR projection images |
US6570951B1 (en) * | 2002-05-14 | 2003-05-27 | Ge Medical Systems Global Technology Company, Llc | Image space compensation scheme for reducing artifacts |
US20040136490A1 (en) * | 2002-07-23 | 2004-07-15 | Edic Peter Michael | Method and apparatus for correcting motion in image reconstruction |
CN101126722A (zh) * | 2007-09-30 | 2008-02-20 | 西北工业大学 | 基于配准模型仿真的锥束ct射束硬化校正方法 |
US20110123081A1 (en) * | 2009-11-25 | 2011-05-26 | David Sebok | Correcting and reconstructing x-ray images using patient motion vectors extracted from marker positions in x-ray images |
CN102982510A (zh) * | 2011-06-21 | 2013-03-20 | 通用电气公司 | Spect图像的运动校正 |
US20130051644A1 (en) * | 2011-08-29 | 2013-02-28 | General Electric Company | Method and apparatus for performing motion artifact reduction |
CN103679642A (zh) * | 2012-09-26 | 2014-03-26 | 上海联影医疗科技有限公司 | 一种ct图像金属伪影校正方法、装置及ct设备 |
US20140212011A1 (en) * | 2013-01-28 | 2014-07-31 | Koninklijke Philips Electronics N.V. | Spect motion-correction |
CN104424659A (zh) * | 2013-08-30 | 2015-03-18 | 西门子公司 | 图像数据的基于图像的运动补偿 |
CN103455989A (zh) * | 2013-09-24 | 2013-12-18 | 南京大学 | 一种结合超声图像提高有限角度ct成像质量的方法 |
US20160104279A1 (en) * | 2014-10-10 | 2016-04-14 | Cedars-Sinai Medical Center | Method and system for "push-button" comprehensive cardiac mr examination using continuous self-gated 3d radial imaging |
US20160110893A1 (en) * | 2014-10-21 | 2016-04-21 | Shenyang Neusoft Medical Systems Co., Ltd. | Correcting a ct scan image |
CN104574292A (zh) * | 2014-11-26 | 2015-04-29 | 沈阳东软医疗***有限公司 | 一种ct图像的校正方法和装置 |
CN104867168A (zh) * | 2015-04-28 | 2015-08-26 | 南京邮电大学 | 基于p范数的压缩感知计算机断层扫描图像重建方法 |
US20190180482A1 (en) * | 2016-07-04 | 2019-06-13 | Hitachi, Ltd. | Image reconstruction device, x-ray ct device, and image reconstruction method |
US20180137657A1 (en) * | 2016-11-16 | 2018-05-17 | University Of Central Florida Research Foundation, Inc. | System and method for motion estimation and compensation in helical computed tomography |
CN106683144A (zh) * | 2016-12-30 | 2017-05-17 | 上海联影医疗科技有限公司 | 一种图像迭代重建方法及装置 |
CN106683146A (zh) * | 2017-01-11 | 2017-05-17 | 上海联影医疗科技有限公司 | 一种图像重建方法和图像重建算法的参数确定方法 |
CN106960429A (zh) * | 2017-02-16 | 2017-07-18 | 中国科学院苏州生物医学工程技术研究所 | 一种ct图像金属伪影校正方法及装置 |
CN107845121A (zh) * | 2017-11-03 | 2018-03-27 | 中国工程物理研究院应用电子学研究所 | 一种探测器偏置扫描中加权伪影的校正方法 |
CN108876730A (zh) * | 2018-05-24 | 2018-11-23 | 沈阳东软医疗***有限公司 | 校正运动伪影的方法、装置及设备和存储介质 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113628293A (zh) * | 2020-05-06 | 2021-11-09 | 美国西门子医疗***股份有限公司 | 全身直接参数化图像重建中的帧间运动校正 |
CN111612867A (zh) * | 2020-06-01 | 2020-09-01 | 上海联影医疗科技有限公司 | 运动伪影校正方法、装置、计算机设备和可读存储介质 |
CN111612867B (zh) * | 2020-06-01 | 2023-08-15 | 上海联影医疗科技股份有限公司 | 运动伪影校正方法、装置、计算机设备和可读存储介质 |
CN111798534A (zh) * | 2020-07-17 | 2020-10-20 | 沈阳先进医疗设备技术孵化中心有限公司 | 图像重建方法、装置、控制台设备及ct*** |
CN111798534B (zh) * | 2020-07-17 | 2024-03-08 | 东软医疗***股份有限公司 | 图像重建方法、装置、控制台设备及ct*** |
WO2022022535A1 (en) * | 2020-07-28 | 2022-02-03 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for image processing |
CN113313649A (zh) * | 2021-06-01 | 2021-08-27 | 上海联影医疗科技股份有限公司 | 图像重建方法及装置 |
CN113313649B (zh) * | 2021-06-01 | 2022-09-16 | 上海联影医疗科技股份有限公司 | 图像重建方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110473269B (zh) | 2023-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110473269A (zh) | 一种图像重建方法、***、设备及存储介质 | |
CN1809841B (zh) | 运动补偿的重建方法、设备与*** | |
US7948503B2 (en) | Techniques for 3-D elastic spatial registration of multiple modes of measuring a body | |
CN109961491A (zh) | 多模态图像截断补偿方法、装置、计算机设备和介质 | |
US10213179B2 (en) | Tomography apparatus and method of reconstructing tomography image | |
US20110262015A1 (en) | Image processing apparatus, image processing method, and storage medium | |
JP2020506742A (ja) | 断層撮影再構成に使用するためのデータのディープラーニングに基づく推定 | |
US20110044559A1 (en) | Image artifact reduction | |
Kim et al. | Cardiac motion correction based on partial angle reconstructed images in x‐ray CT | |
CN102783966A (zh) | 用于运动校正的医学图像的方法和器件 | |
AU2019271915A1 (en) | Method and system for motion correction in CT imaging | |
CN108876730B (zh) | 校正运动伪影的方法、装置及设备和存储介质 | |
CN104252714A (zh) | 时变数据的重建 | |
US10013778B2 (en) | Tomography apparatus and method of reconstructing tomography image by using the tomography apparatus | |
CN107784684B (zh) | 一种锥束ct三维重建方法及*** | |
CN112150571A (zh) | 一种图像运动伪影消除方法、装置、设备及存储介质 | |
Müller et al. | Evaluation of interpolation methods for surface‐based motion compensated tomographic reconstruction for cardiac angiographic C‐arm data | |
CN108074219A (zh) | 一种图像校正方法、装置及医疗设备 | |
Zhang et al. | Comparison of different attenuation correction methods for dual gating myocardial perfusion SPECT/CT | |
CN111476854B (zh) | 图像重建方法、装置、终端及存储介质 | |
Sun et al. | Quantification of local reconstruction accuracy for helical CT with motion correction | |
CN110853113B (zh) | 基于bpf的tof-pet图像重建算法及重建*** | |
Feng et al. | Automated least-squares calibration of the coregistration parameters for a micro PET-CT system | |
Schäfer et al. | Limited angle C-arm tomography and segmentation for guidance of atrial fibrillation ablation procedures | |
US11712215B2 (en) | Devices, systems, and methods for motion-corrected medical 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 |
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. |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |