CN102763138B - 辐射治疗中的运动校正 - Google Patents
辐射治疗中的运动校正 Download PDFInfo
- Publication number
- CN102763138B CN102763138B CN201080051809.7A CN201080051809A CN102763138B CN 102763138 B CN102763138 B CN 102763138B CN 201080051809 A CN201080051809 A CN 201080051809A CN 102763138 B CN102763138 B CN 102763138B
- Authority
- CN
- China
- Prior art keywords
- image data
- motion
- dissection
- motion model
- interest
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/254—Analysis of motion involving subtraction of images
-
- 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]
-
- 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/10108—Single photon emission computed tomography [SPECT]
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Abstract
一种诊断成像系统包括断层摄影扫描器10,所述断层摄影扫描器生成解剖和功能图像数据集。适配单元50基于在多个运动阶段内采集的运动平均的体积图像表示使运动模型适配所述感兴趣对象的几何性质。在所述多个运动阶段利用运动模型从所述解剖投影图像数据模拟虚拟图像数据。比较单元54确定实际和虚拟解剖图像数据之间的差异。如果该差异满足停止标准,使用运动模型校正所采集的功能图像数据,并从其重建经校正的功能图像。如果不是,迭代地更新运动模型,直到所述差异满足停止标准。
Description
技术领域
本申请涉及诊断成像领域。本发明尤其适于与组合式X射线计算机断层摄影(CT)扫描器和发射断层摄影扫描器结合使用,所述发射断层摄影诸如是正电子发射断层摄影(PET)和单光子发射计算机断层摄影(SPECT)。
背景技术
在诊断核成像中,随着放射性核素通过患者的血流来研究放射性核素分布,以对循环系统成像,或者对积聚所注入的放射性药剂的特定器官成像。在单光子发射计算机断层摄影(SPECT)中,使用一个或多个辐射探测器,通称为伽马(gamma)照相机,通过由放射性衰变事件导致的辐射发射来探测放射性药剂。通常,每个伽马照相机都包括辐射探测器阵列和设置于辐射探测器阵列前方的蜂房型准直器。蜂房型准直器界定线性或小角度圆锥视线,使得所探测的辐射包括投影数据。如果在一定范围的视角上,例如在180°或360°角范围上移动伽马照相机,那么可以利用滤波反向投影、期望值最大化或者另一种成像技术将所得的投影数据重建成患者体内放射性药剂分布的图像。有利地,可以设计放射性药剂以使其集中在选定组织中,从而提供那些选定组织的优先成像。
在正电子发射断层摄影(PET)中,放射性药剂的放射性衰变事件产生正电子。每个正电子与电子交互作用以产生正电子-电子湮灭事件,该湮灭事件发射两个反向的伽马射线。使用符合(coincidence)探测电路,围绕成像患者的辐射探测器的环形阵列探测与正电子-电子湮灭相对应的符合反向伽马射线事件。连接两个符合探测的响应线(LOR)包含正电子-电子湮灭事件的位置。这样的响应线类似于投影数据,并且能够被重建以产生二维或三维图像。在飞行时间PET(TOF-PET)中,使用两个符合γ射线事件的探测之间的小时间差异来沿着LOR(响应线)定位湮灭事件。
SPECT和PET两种成像技术都有一个问题,即,放射性核素与探测器之间由患者的解剖结构造成的光子吸收和散射使所得的图像失真。为了获得更为准确的核图像,利用透射计算机断层摄影技术进行直接透射辐射测量。透射数据用于构造整个身体内的密度差的衰减图并用于校正所发射光子的吸收。在过去,放射性同位素线或点源与探测器相对放置,使得探测器能够收集透射数据。在有和无患者时,两个值的比率用于校正不均匀密度,这可能导致图像噪声、图像伪影、图像失真,并可能掩盖重要特征。
另一种技术使用X射线CT扫描数据生成更准确的衰减图。由于与软组织相比,X射线和伽马射线两者都受到硬组织,诸如骨骼甚至合成材料植入物的更强的衰减,所以可以使用CT数据来估计放射性药剂所发射的伽马射线的衰减图。通常,在所发射伽马射线的适当能量下,使用取决于能量的缩放因子将CT像素值(单位Hounsfield(HU))转换成线性衰减系数(LAC)。
在过去,核扫描器和CT扫描器以固定的关系彼此相邻地永久安装,并共享公共的患者支撑物。将患者从CT扫描器的检查区域平移到核扫描器的检查区域。然而,由于患者的潜在移动或者CT扫描器和核扫描器之间的重新定位,这种技术在核图像和CT图像之间的对准方面带来不确定性。
为了消除对准问题,当前的系统将CT和核成像系统安装到公共的扫描架。然而,该设计意味着扫描架的速度限于每次旋转数十秒。如果患者在CT采集期间屏息,可以在CT数据中消除或减少运动。这出现了一个问题,即核成像采集时间比屏息以生成充分多数据更长。因此,患者自由呼吸。在屏息CT扫描期间患者的几何性质与自由呼吸的核扫描的几何性质不匹配。这会导致重建伪影,因为在衰减图与采集核数据的若干分钟内所采集的发射数据之间,尤其是在运动更大的区域中,例如隔膜、心壁等中,存在不匹配。
发明内容
本申请提供了一种对核成像中运动对象的衰减和散射校正的新的改进的方法和设备,其克服了上述问题和其他问题。
根据一个方面,提供了一种用于生成运动模型的方法。在感兴趣对象的多个运动阶段期间采集解剖投影图像数据集。将所采集的解剖投影图像数据集重建成运动平均的解剖图像表示。基于所述运动平均的体积图像表示使运动模型的几何性质适配所述感兴趣对象的几何性质。在所述多个运动阶段利用运动模型从所述运动平均的解剖图像表示来模拟解剖投影图像数据。基于所采集的解剖投影图像数据集和所模拟的解剖图像数据之间的差异来更新所述运动模型。
根据另一方面,一种被配置成执行用于生成运动模型的方法的处理器。
根据另一方面,一种诊断成像系统,包括断层摄影扫描器,所述断层摄影扫描器连续生成解剖和功能图像数据集。所述诊断成像系统包括被编程控制以执行生成运动模型的方法的一个或多个处理器。
根据另一方面,一种诊断成像系统,包断层摄影扫描器,所述断层摄影扫描器生成感兴趣对象的解剖和功能图像数据集。解剖重建单元将所述解剖投影图像数据集重建成运动平均的解剖图像表示。适配单元基于所述运动平均的体积图像表示使运动模型适配所述感兴趣对象的几何性质。模拟单元在所述多个运动阶段利用运动模型从所述运动平均的解剖图像表示来模拟解剖投影图像数据。比较单元确定所采集的解剖投影图像数据集与所模拟的解剖图像数据之间的差异。运动模型更新单元基于由所述比较单元确定的差异来更新所述运动模型。
一个优点在于,可以在多个运动阶段内采集感兴趣对象的图像数据。
另一优点在于,针对采集运动中的感兴趣对象的图像数据,改善了信噪比(SNR)。
另一优点在于,可以在断层摄影扫描器的扫描架旋转期间采集感兴趣对象的图像数据。
另一优点在于,在采集投影数据期间减少了受检者的辐射暴露。
另一优点在于,可以针对感兴趣对象的各个运动阶段采集用于校正发射数据的校正数据。
附图说明
在阅读和理解下文的详细描述之后,本领域技术人员将认识到本发明的其他优点。
本发明可以采用各种部件和部件布置,以及各种步骤和步骤安排的形式。附图仅用于图示说明优选实施例,而不应被解释为限制本发明。
图1是具有运动建模单元的组合式SPECT/CT单扫描架系统的示意图;以及
图2是用于生成运动模型的方法的流程图。
具体实施方式
参考图1,诊断成像系统10同时和/或独立地执行X射线计算机断层摄影(XCT)和核成像,诸如PET或SPECT。成像系统10包括界定患者接收膛14的静止外壳12。由外壳12支撑的可旋转扫描架16布置在膛的周围,以界定公共检查区域18。沿纵向和/或垂直地调节支撑待成像和/或检查的患者或受检者22的患者支撑物20,以实现患者在检查区域中的期望定位。
为了提供XCT成像能力,安装在可旋转扫描架16上的X射线组件24包括诸如X射线管的X射线源26和准直器或光阀(shutter)组件28。准直器将来自X射线源26的辐射准直成锥形或楔形射束、一个或多个基本平行的扇形射束等。光阀使射束选通和关闭。X射线探测器30,诸如固态、平板探测器,安装在可旋转扫描架16上,并与辐射组件24相对。在扫描架旋转时,X射线组件24和探测器30一致绕检查区域18旋转以采集跨越半转、完整360°旋转、多转或更小弧的XCT投影数据。每个XCT投影都指示沿X射线组件24和X射线探测器30之间的线性路径的X射线衰减。所采集的XCT投影数据被存储在数据缓存器32中并由XCT重建处理器34处理成XCT图像表示,然后存储在XCT图像存储器36中。结合到一起,X射线源、准直器/光阀组件、探测器和重建处理器界定了用于生成解剖图像的模块(means)。
为了提供功能性核成像能力,将至少两个核探测器头40a、40b,诸如单光子发射断层摄影(SPECT)探测器,可活动地安装到旋转扫描架16上。安装X射线组件24和核探测器头40a、40b允许在无需移动患者22的情况下由两种模态对检查区域18成像。在一个实施例中,探测器头由安装到旋转扫描架16的机器人组件(未示出)可移动地支撑。机器人组件能够绕着患者22定位探测器头以采集跨越不同角范围,例如90°偏移、彼此相对的180°等的视图。每个SPECT探测器头都包括准直器,使得已知每个探测到的辐射事件都沿着可识别的线性或小角度圆锥形视线发生,从而所采集的辐射包括投影数据。所采集的SPECT投影数据被存储在数据缓存器42中并由SPECT重建处理器44将其处理成SPECT图像表示,然后存储在SPECT图像存储器46中。结合在一起,SPECT探测器头和SPECT重建处理器界定了用于生成功能图像的模块。
在另一实施例中,功能性成像模块包括正电子发射断层摄影(PET)探测器。绕着患者接收膛14布置PET探测器的一个或多个环,以从其中接收伽马辐射。所探测到的符合辐射事件对界定PET投影数据,PET投影数据存储在数据缓存器中并由PET重建处理器处理成PET图像表示,然后存储在PET图像存储器中。结合在一起,PET探测器环和PET重建处理器界定了用于生成功能图像的模块。
通常,在功能性核成像中,从受检者的透射数据生成衰减图。衰减图用于校正所采集的功能性投影数据的衰减,即本来会包括在功能图像中的光子,由于密度更大的组织吸收更多的发射光子,造成图像变化。在多扫描架系统中,在屏息采集期间从解剖成像系统采集透射数据。然后将受检者重新定位到功能性成像系统中,功能性成像系统通常与解剖成像系统相邻并共享同一患者支撑物。
即使在两个成像系统彼此非常接近时,也可能发生重新定位误差,这降低了衰减图的精确度。此外,功能性成像时间充分长,持续若干个呼吸周期。另一方面,可以在充分短的时间内生成解剖图像,可以在单次屏息期间生成解剖图像。然而,因为在呼吸阶段的整个范围内生成功能图像数据;而解剖图像数据是在单次呼吸阶段中生成的,所以解剖和功能图像表示并非在所有呼吸阶段中都匹配。这会导致图像伪影。为了克服这些问题,从解剖图像数据生成感兴趣对象的运动模型。利用运动模型生成针对感兴趣对象的每个运动阶段的衰减图。
继续参考图1,由控制器50操作诊断成像扫描器以执行成像序列。在将受检者定位在检查区域18中之后,在对象经历多个阶段的呼吸或其他运动的同时,例如在对象经历呼吸周期的同时,利用解剖图像生成模块,成像序列在多个投影角采集感兴趣对象的投影成像数据集。在数据缓存器32中存储所采集的解剖图像投影数据集。解剖重建处理器34从所采集的解剖投影图像数据集重建至少一个运动平均的解剖体积表示。在解剖图像存储器36中存储重建的运动平均的解剖体积表示。由于解剖图像投影数据是在多个运动阶段期间采集的,所以所得的运动平均的体数据表示是感兴趣对象的模糊图像。例如,如果感兴趣对象是位于一个肺中的肿瘤,它将由于呼吸而经历周期运动。与在单次屏息中旋转扫描架以收集完整的数据集的屏息成像序列不同,本装置允许受检者在采集期间自由呼吸以适应单次旋转长于典型屏息的扫描架16。
从运动平均的体数据表示,感兴趣对象的模糊表面或边界表示感兴趣对象的运动阶段。因此,界定用于适配的模块的适配单元50基于运动平均的体数据表示自动或半自动地使运动模型适配感兴趣对象的几何性质。适配单元包括一般运动模型的库,例如,基于非均匀有理B样条(NURBS)的核计算机轴向断层摄影(NCAT)和X射线计算机轴向断层摄影(XCAT)计算体模,基于感兴趣对象的几何性质从其确定最佳匹配。利用已知的分割和/或拟合方法,诸如针对三维(3D)区域的多边形网格或点云(CoP)拟合方案,将所确定的最佳匹配运动模型拟合到感兴趣对象的几何性质,适配单元利用来自运动平均的解剖图像表示、解剖成像扫描的持续时间和/或与解剖图像投影数据相关联的时间戳的其模糊边界确定感兴趣对象的运动阶段。
界定了用于模拟的模块的模拟单元52,基于运动模型生成虚拟解剖投影图像数据。在本领域中已知用于生成3D患者图像或模型的二维(2D)解剖投影数据的模拟方法,例如基于蒙特卡洛(MC)的方法,包括康普顿和/或瑞利散射建模等。
界定了用于比较的模块的比较单元54,通过基于解剖图像的虚拟二维(2D)投影与在已知呼吸阶段中在对应角度实际采集的2D解剖投影图像数据之间的差异生成在每个投影角的变形场,来比较虚拟和实际采集的解剖投影图像数据。通过分析虚拟和所采集的解剖图像或投影之间的差异,比较单元导出每个投影角的2D变形场。比较可以以基于界标的变形计算或2D弹性对准计算为基础,在基于界标的变形计算中,针对每个投影角计算每个界标的运动的两个分量,2D弹性对准计算计算每个投影角的2D变形矢量场。
界定了用于几何校正的模块的几何校正单元56,组合在所有投影角的2D变形场以形成相容(consistent)的3D变形场。几何校正单元执行的组合通过导出最可能的3D变形场,可以基于最大似然(ML)运动模型,最可能的3D变形场最好地解释了所观测到的2D变形,或者可以基于纯几何方法,这种方法在不同视角中求解个体界标投影线的3D交点。几何校正单元确定每个运动阶段对运动模型的几何校正,以便使所采集的解剖投影图像数据与模拟的投影图像数据之间的差异最小化。适配单元50施加几何校正,使得运动模型与感兴趣对象的几何性质一致。
结合在一起,适配单元50、模拟单元52、比较单元54和几何校正单元56界定了用于生成运动模型的模块。迭代地重复生成运动模型,直到达到预选的品质因数或停止标准。
一旦生成了达标的运动模型,扫描器控制器继续成像序列以在对象经历多个运动阶段的同时利用功能图像生成模块采集感兴趣对象的功能性成像数据集。或者,可以与解剖图像投影数据同时生成功能性成像数据并加以存储,直到生成3D运动模型。通常,为待成像的受检者注射一种或多种放射性药剂或放射性同位素示踪剂。这样的示踪剂的范例是Tc-99m、Ga-67、In-111和I-123。感兴趣对象之内存在示踪剂从感兴趣对象生成发射辐射事件,由核探测器头40a、40b探测到这种事件。在数据缓存器42中存储所采集的功能图像数据集。界定了用于感测运动的模块的运动感测装置60,在采集所述功能图像数据集期间生成运动信号。运动信号指示正在采集功能图像数据的同时感兴趣对象的当前运动阶段。运动感测装置的范例包括呼吸带、光学跟踪系统、心电图(ECG)、脉搏计等。所生成的运动信号用于将所采集的功能图像数据分装(bin)到相等患者几何性质,即相同运动阶段的集合中。
利用来自运动模型生成模块的运动模型和所生成的运动信号,界定了用于校正的模块的校正单元62,针对感兴趣对象的每个运动阶段校正该功能图像数据集。校正类型的范例包括衰减校正、散射校正、分体积校正等。为了校正衰减,校正单元基于所生成的运动模型为感兴趣对象的每个运动阶段生成衰减图。利用衰减图校正功能图像数据的每种分装,该衰减图对应于与该分装相关联的运动阶段。因此,校正单元基于所生成的运动模型为感兴趣对象的每个运动阶段生成散射校正函数。利用散射校正函数校正功能图像数据的每种分装,该散射校正函数对应于与该分装相关联的运动阶段。校正单元基于所生成的运动模型为感兴趣对象的每个运动阶段生成标准摄取值(SUV)校正因数。利用SUV校正因数校正功能图像数据的种分装,该SUV校正因数对应于与该分装相关联的运动阶段。应当认识到,还想到了用于衰减、散射和分体积校正的其他方法。
在更为具体的范例中,运动模型是四维(4D)模型,即针对每次呼吸或其他运动阶段的一堆3D衰减图。在探测器头收集数据时,利用在探测器头上的位置、探测器头的角位置和运动阶段对每个辐射事件进行编码。在重建期间,针对对应的运动阶段,通过运动阶段对数据进行分装,并利用衰减图加以校正。
功能重建处理器44从经校正的功能图像数据集重建至少一个功能图像表示。重建的功能图像表示存储在功能图像存储器46中。工作站或图形用户接口70包括显示装置和用户输入装置,医生能够使用它们选择扫描序列和规程、显示图像数据等。
任选的图像组合器72将解剖图像表示和功能图像表示组合成一个或多个组合图像表示,以同时加以显示。例如,可以以不同颜色中叠加图像,可以在解剖图像表示上叠加功能图像表示的轮廓或特征,可以在功能图像表示上叠加解剖图像表示的分割解剖结构的轮廓或特征,可以以共同的比例尺并排显示功能和解剖图像表示,等等。组合图像存储在组合图像存储器74中。
参考图1,扫描器控制器50包括利用计算机程序编程程控的处理器(所述计算机程序被存储在计算机可读介质上),以执行根据图示的流程图的方法,该方法可以包括,但不限于,控制功能和解剖成像模块,即光子发射断层摄影扫描器和X射线断层摄影扫描器。适当的计算机可读介质包括光学、磁性或固态存储器,例如CD、DVD、硬盘、软盘、RAM、闪存等。
根据图2,用于生成运动模型的方法包括采集解剖图像数据。将所采集的解剖图像数据重建成解剖图像表示。使运动模型适配感兴趣对象,在解剖图像表示中突出显示。通过利用多个运动阶段的运动模型模拟所采集的解剖图像数据来生成虚拟解剖图像数据。实际采集的解剖图像数据是要虚拟解剖图像数据。如果实际和虚拟解剖图像数据之间的差异低于阈值或者满足停止标准,使用运动模型来校正正确的功能图像数据,并从其重建功能图像表示。如果实际和虚拟解剖图像数据之间的差异不低于阈值或者不满足停止标准,基于差异来更新运动模型,并迭代地重复模拟直到生成适当的运动模型。
已经参考优选实施例描述了本发明。他人在阅读和理解以上详细描述之后可能想到修改和变更。应当将本发明解释为包括所有这样的修改和变更,只要它们在所附权利要求或其等价要件的范围之内。
Claims (21)
1.一种用于生成运动模型的方法,包括:
在感兴趣对象的多个运动阶段期间采集解剖投影图像数据集;
将所述解剖投影图像数据集重建成运动平均的解剖体积图像表示;
基于所述运动平均的体积图像表示使运动模型的几何性质适配所述感兴趣对象的几何性质;
在所述多个运动阶段利用所述运动模型从所述运动平均的解剖图像表示来模拟所述解剖投影图像数据;以及
基于所采集的解剖投影图像数据集与所模拟的解剖图像数据之间的差异来更新所述运动模型。
2.根据权利要求1所述的方法,还包括:
迭代地重复模拟所述解剖投影图像数据然后更新所述运动模型的步骤,直到达到停止标准。
3.根据权利要求1或2所述的方法,其中,所述解剖投影图像数据集是在多个投影角中的每个处采集的。
4.根据权利要求3所述的方法,其中,更新所述运动模型的步骤还包括:
基于在所述投影角中的每个处的所述解剖投影图像数据集与所模拟的解剖投影图像数据集之间的差异来生成在对应投影角处的变形场;
组合在每个投影角处的所述变形场以形成三维(3D)变形场;以及
基于所述三维(3D)变形场更新所述运动模型的几何性质。
5.根据权利要求4所述的方法,还包括:
在所述感兴趣对象的所述多个运动阶段期间采集功能图像数据集;
针对每个运动阶段基于所述运动模型校正所述功能图像数据集;以及
将经校正的功能图像数据集重建成所述感兴趣对象的至少一个经校正的功能图像表示。
6.根据权利要求5所述的方法,还包括:
在采集所述功能图像数据集期间从运动感测装置采集运动信号,所述运动信号表征所述感兴趣对象的每个运动阶段。
7.根据权利要求6所述的方法,其中,校正所述功能图像数据集的步骤还包括:
根据所采集的运动信号,基于针对每个运动阶段的所述三维(3D)变形场来生成衰减图;以及
根据针对每个运动阶段的所述衰减图校正所述功能图像数据集的衰减和散射。
8.一种用于生成运动模型的装置,包括:
用于在感兴趣对象的多个运动阶段期间采集解剖投影图像数据集的模块;
用于将所述解剖投影图像数据集重建成运动平均的解剖体积图像表示的模块;
用于基于所述运动平均的体积图像表示使运动模型的几何性质适配所述感兴趣对象的几何性质的模块;
用于在所述多个运动阶段利用所述运动模型从所述运动平均的解剖图像表示来模拟所述解剖投影图像数据的模块;以及
用于基于所采集的解剖投影图像数据集与所模拟的解剖图像数据之间的差异来更新所述运动模型的模块。
9.根据权利要求8所述的装置,还包括:
用于迭代地重复模拟所述解剖投影图像数据然后更新所述运动模型的步骤,直到达到停止标准的模块。
10.根据权利要求8或9所述的装置,其中,所述解剖投影图像数据集是在多个投影角中的每个处采集的。
11.根据权利要求10所述的装置,其中,用于更新所述运动模型的模块还包括:
用于基于在所述投影角中的每个处的所述解剖投影图像数据集与所模拟的解剖投影图像数据集之间的差异来生成在对应投影角处的变形场的模块;
用于组合在每个投影角处的所述变形场以形成三维(3D)变形场的模块;以及
用于基于所述三维(3D)变形场更新所述运动模型的几何性质的模块。
12.根据权利要求11所述的装置,还包括:
用于在所述感兴趣对象的所述多个运动阶段期间采集功能图像数据集的模块;
用于针对每个运动阶段基于所述运动模型校正所述功能图像数据集的模块;以及
用于将经校正的功能图像数据集重建成所述感兴趣对象的至少一个经校正的功能图像表示的模块。
13.根据权利要求12所述的装置,还包括:
用于在采集所述功能图像数据集期间从运动感测装置采集运动信号的模块,所述运动信号表征所述感兴趣对象的每个运动阶段。
14.根据权利要求13所述的装置,其中,用于校正所述功能图像数据集的模块还包括:
用于根据所采集的运动信号,基于针对每个运动阶段的所述三维(3D)变形场来生成衰减图的模块;以及
用于根据针对每个运动阶段的所述衰减图校正所述功能图像数据集的衰减和散射的模块。
15.一种诊断成像系统,包括:
断层摄影扫描器(10),其连续生成解剖和功能图像数据集;以及
一个或多个处理器,其被编程控制以执行根据权利要求1-7中的任一项所述的方法的步骤。
16.一种诊断图像扫描器,包括:
断层摄影扫描器(10),其在感兴趣对象的多个运动阶段期间采集解剖投影图像数据集;
解剖重建单元(34),其将所述解剖投影图像数据集重建成运动平均的解剖图像表示;
适配单元(50),其基于所述运动平均的体积图像表示使运动模型适配所述感兴趣对象的几何性质;
模拟单元(52),其在所述多个运动阶段利用所述运动模型从所述运动平均的解剖图像表示来模拟解剖投影图像数据;以及
比较单元(54),其确定所采集的解剖投影图像数据集与所模拟的解剖图像数据之间的差异;以及
运动模型更新单元(56),其基于所述比较单元(54)所确定的所述差异来更新所述运动模型。
17.根据权利要求16所述的诊断图像扫描器,其中:
所述模拟单元(52)迭代地重复利用经更新的运动模型对所述解剖投影图像数据的所述模拟,直到达到停止标准。
18.根据权利要求16或17所述的诊断图像扫描器,其中,所述断层摄影扫描器(10)在每个投影角处都采集一次所述解剖投影图像数据集。
19.根据权利要求18所述的诊断图像扫描器,其中:
所述比较单元(54)基于在所述投影角的每个处的所述解剖投影图像数据集与所模拟的解剖投影图像数据之间的差异来生成在对应投影角处的变形场;并且
所述运动模型更新单元(56)组合在每个投影角处的所述变形场以形成三维(3D)变形场并基于所述三维(3D)变形场更新所述运动模型的几何性质。
20.根据权利要求19所述的诊断图像扫描器,其中,所述断层摄影扫描器(10)在所述感兴趣对象的所述多个运动阶段期间采集功能图像数据集,所述诊断图像扫描器还包括:
校正单元(62),其针对每个运动阶段基于所述运动模型校正所述功能图像数据集;以及
功能重建单元(44),其将经校正的所述功能图像数据集重建成所述感兴趣对象的至少一个经校正的功能图像表示。
21.根据权利要求20所述的诊断图像扫描器,还包括:
运动感测装置(60),其在采集所述功能图像数据集期间采集运动信号,所述运动信号表征所述感兴趣对象的每个运动阶段;并且
其中,所述校正单元(62)根据所采集的运动信号基于针对每个运动阶段的所述三维(3D)变形场来生成衰减图;并且
所述校正单元(62)根据针对每个运动阶段的所述衰减图来校正所述功能图像数据集的衰减和散射。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US26217209P | 2009-11-18 | 2009-11-18 | |
US61/262,172 | 2009-11-18 | ||
PCT/IB2010/054665 WO2011061644A1 (en) | 2009-11-18 | 2010-10-14 | Motion correction in radiation therapy |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102763138A CN102763138A (zh) | 2012-10-31 |
CN102763138B true CN102763138B (zh) | 2016-02-17 |
Family
ID=43501165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201080051809.7A Expired - Fee Related CN102763138B (zh) | 2009-11-18 | 2010-10-14 | 辐射治疗中的运动校正 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20120278055A1 (zh) |
EP (1) | EP2502204A1 (zh) |
CN (1) | CN102763138B (zh) |
RU (1) | RU2012124998A (zh) |
WO (1) | WO2011061644A1 (zh) |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5400546B2 (ja) * | 2009-09-28 | 2014-01-29 | 株式会社日立メディコ | X線ct装置 |
RU2582887C2 (ru) * | 2010-11-23 | 2016-04-27 | Конинклейке Филипс Электроникс Н.В. | Калибровка рет с изменяющимися интервалами совпадений |
JP5947813B2 (ja) * | 2011-01-05 | 2016-07-06 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | ゲート信号を持つリストモードpetデータにおいて運動を検出し補正する方法及び装置 |
US20150150525A1 (en) * | 2012-06-11 | 2015-06-04 | Nassir Navab | Dynamic nuclear emission and x-ray imaging device and respective imaging method |
DE102012213551A1 (de) * | 2012-08-01 | 2014-02-06 | Siemens Aktiengesellschaft | Verfahren zur bewegungsgemittelten Schwächungskorrektur und Magnetresonanz-Anlage |
KR101461099B1 (ko) * | 2012-11-09 | 2014-11-13 | 삼성전자주식회사 | 자기공명영상장치 및 기능적 영상획득방법 |
EP2760028B1 (en) * | 2013-01-23 | 2018-12-12 | Samsung Electronics Co., Ltd | Radiation generator |
US9443346B2 (en) * | 2013-07-23 | 2016-09-13 | Mako Surgical Corp. | Method and system for X-ray image generation |
CN106030657B (zh) * | 2014-02-19 | 2019-06-28 | 皇家飞利浦有限公司 | 医学4d成像中的运动自适应可视化 |
JP6722652B2 (ja) | 2014-07-28 | 2020-07-15 | インテュイティブ サージカル オペレーションズ, インコーポレイテッド | 手術中のセグメンテーションについてのシステム及び方法 |
US9763631B2 (en) | 2014-09-17 | 2017-09-19 | General Electric Company | Systems and methods for imaging plural axial locations |
DE102015206362B3 (de) * | 2015-04-09 | 2016-07-21 | Siemens Healthcare Gmbh | Multizyklische dynamische CT-Bildgebung |
US9965875B2 (en) * | 2016-06-21 | 2018-05-08 | Carestream Health, Inc. | Virtual projection image method |
JP6799292B2 (ja) * | 2017-07-06 | 2020-12-16 | 株式会社島津製作所 | 放射線撮影装置および放射線画像検出方法 |
CN108389232B (zh) * | 2017-12-04 | 2021-10-19 | 长春理工大学 | 基于理想视点的非规则表面投影图像几何校正方法 |
US11568581B2 (en) | 2018-01-27 | 2023-01-31 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for correcting mismatch induced by respiratory motion in positron emission tomography image reconstruction |
US10504250B2 (en) | 2018-01-27 | 2019-12-10 | Uih America, Inc. | Systems and methods for correcting mismatch induced by respiratory motion in positron emission tomography image reconstruction |
US10492738B2 (en) | 2018-02-01 | 2019-12-03 | Siemens Medical Solutions Usa, Inc. | Motion detection for nuclear medicine imaging |
EP3547262A1 (en) * | 2018-03-28 | 2019-10-02 | Koninklijke Philips N.V. | Tomographic x-ray image reconstruction |
EP3628225B1 (de) * | 2018-09-26 | 2021-03-31 | Siemens Healthcare GmbH | Verfahren zur aufnahme von bilddaten und medizinisches bildgebungssystem |
JP7330833B2 (ja) * | 2019-09-20 | 2023-08-22 | 株式会社日立製作所 | 放射線撮像装置および放射線治療装置 |
CN110842918B (zh) * | 2019-10-24 | 2020-12-08 | 华中科技大学 | 一种基于点云伺服的机器人移动加工自主寻位方法 |
US11410354B2 (en) * | 2020-02-25 | 2022-08-09 | Uih America, Inc. | System and method for motion signal recalibration |
CN111476897B (zh) * | 2020-03-24 | 2023-04-18 | 清华大学 | 基于同步扫描条纹相机的非视域动态成像方法及装置 |
US11222447B2 (en) * | 2020-05-06 | 2022-01-11 | Siemens Medical Solutions Usa, Inc. | Inter-frame motion correction in whole-body direct parametric image reconstruction |
EP3961567A1 (en) * | 2020-08-27 | 2022-03-02 | Koninklijke Philips N.V. | Apparatus, method and computer program for registering pet images |
WO2022170607A1 (zh) * | 2021-02-10 | 2022-08-18 | 北京大学 | 一种定位图像的转化系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101542530A (zh) * | 2006-11-22 | 2009-09-23 | 皇家飞利浦电子股份有限公司 | 基于有限数据集合的图像生成 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007015199A2 (en) * | 2005-08-04 | 2007-02-08 | Koninklijke Philips Electronics, N.V. | Motion compensation in functional imaging |
US20080095414A1 (en) | 2006-09-12 | 2008-04-24 | Vladimir Desh | Correction of functional nuclear imaging data for motion artifacts using anatomical data |
-
2010
- 2010-10-14 CN CN201080051809.7A patent/CN102763138B/zh not_active Expired - Fee Related
- 2010-10-14 US US13/503,933 patent/US20120278055A1/en not_active Abandoned
- 2010-10-14 RU RU2012124998/08A patent/RU2012124998A/ru unknown
- 2010-10-14 EP EP10777106A patent/EP2502204A1/en not_active Withdrawn
- 2010-10-14 WO PCT/IB2010/054665 patent/WO2011061644A1/en active Application Filing
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101542530A (zh) * | 2006-11-22 | 2009-09-23 | 皇家飞利浦电子股份有限公司 | 基于有限数据集合的图像生成 |
Non-Patent Citations (2)
Title |
---|
PET/CT成像原理概述;李广义等;《医学影像学杂志》;20040831;第14卷(第8期);第681-684页 * |
Toward time resolved 4D cardiac CT imaging with patient dose reduction: estimating the global heart motion;Katsuyuki Taguchi et al.;《Proceedings of SPIE》;20060302;第6142卷;第61420J-1页至第61420J-9页 * |
Also Published As
Publication number | Publication date |
---|---|
RU2012124998A (ru) | 2013-12-27 |
EP2502204A1 (en) | 2012-09-26 |
CN102763138A (zh) | 2012-10-31 |
US20120278055A1 (en) | 2012-11-01 |
WO2011061644A1 (en) | 2011-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102763138B (zh) | 辐射治疗中的运动校正 | |
US7813783B2 (en) | Methods and systems for attenuation correction in medical imaging | |
EP2210238B1 (en) | Apparatus and method for generation of attenuation map | |
JP6192542B2 (ja) | Spect/ctシステムのための反復コーンビームct再構成のための打ち切り補正 | |
US8761478B2 (en) | System and method for tomographic data acquisition and image reconstruction | |
CN103315760B (zh) | 用于核医疗成像中基于发射数据进行衰减补偿的系统和方法 | |
CN104252714B (zh) | 时变数据的重建 | |
US7729467B2 (en) | Methods and systems for attentuation correction in medical imaging | |
US7991450B2 (en) | Methods and systems for volume fusion in diagnostic imaging | |
JP5254810B2 (ja) | リストモードデータに基づく局所動き補償 | |
CN105989621B (zh) | 用于在图像重构中执行联合估计技术的方法和系统 | |
US20140056500A1 (en) | Generating Attenuation Correction Maps for Combined Modality Imaging Studies and Improving Generated Attenuation Correction Maps Using MLAA and DCC Algorithms | |
CN103533892A (zh) | 运动补偿成像 | |
JP5975748B2 (ja) | 核画像におけるアーティファクト除去 | |
CN108242068A (zh) | 用于正电子发射断层成像中的散射校正的方法和系统 | |
Feng et al. | Use of three-dimensional Gaussian interpolation in the projector/backprojector pair of iterative reconstruction for compensation of known rigid-body motion in SPECT | |
EP3804625A1 (en) | Internal dose tomography | |
CN107886554A (zh) | 流数据的重构 | |
James | Visual Dissection and Motion Field Mosaics for Breathing Motion Correction in Nuclear Medicine Imaging | |
Schuh et al. | PET-/SPECT-MRI Attenuation Correction using Image Registration | |
Qiao | Approaches for respiratory motion compensation in positron emission tomography imaging |
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: 20160217 Termination date: 20161014 |
|
CF01 | Termination of patent right due to non-payment of annual fee |