CN103329168A - 针对spect/ct系统的迭代锥形射束ct重建的截断补偿 - Google Patents

针对spect/ct系统的迭代锥形射束ct重建的截断补偿 Download PDF

Info

Publication number
CN103329168A
CN103329168A CN2012800064587A CN201280006458A CN103329168A CN 103329168 A CN103329168 A CN 103329168A CN 2012800064587 A CN2012800064587 A CN 2012800064587A CN 201280006458 A CN201280006458 A CN 201280006458A CN 103329168 A CN103329168 A CN 103329168A
Authority
CN
China
Prior art keywords
fov
data
experimenter
projection
transverse axis
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
Application number
CN2012800064587A
Other languages
English (en)
Other versions
CN103329168B (zh
Inventor
E·S·汉西斯
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN103329168A publication Critical patent/CN103329168A/zh
Application granted granted Critical
Publication of CN103329168B publication Critical patent/CN103329168B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computerised tomographs
    • A61B6/037Emission tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computerised tomographs
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/488Diagnostic techniques involving pre-scan acquisition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5294Devices using data or image processing specially adapted for radiation diagnosis involving using additional data, e.g. patient information, image labeling, acquisition parameters
    • 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computerised tomographs
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/432Truncation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection

Abstract

一种多模态成像系统(10)包括锥形射束计算机断层摄影扫描器(24、30)和核成像扫描器,所述锥形射束计算机断层摄影扫描器采集检查区域(18)中的受试者(22)的CT投影数据,所述核成像扫描器同时/随后采集所述检查区域中的所述受试者的核投影数据。CT重建处理器(34)被编程为执行以下步骤:在所述CT投影数据中,定义视场(FoV),所述视场具有横截轴方向的体素网格;确定所述受试者的最大横截轴跨度;通过将所述FoV的所述体素网格扩展到所述FoV外部的至少一个扩展区域,生成扩展FoV,所述扩展区域至少囊括所确定的最大横截轴跨度和所述横截轴方向中的所有衰减;以及将所述CT投影数据迭代地重建为所述扩展FoV的衰减图。至少一个核重建处理器(44)被编程为基于迭代重建的衰减图校正所采集的核投影数据。

Description

针对SPECT/CT系统的迭代锥形射束CT重建的截断补偿
技术领域
本申请涉及医学成像领域。其具体应用于计算机断层摄影(CT)成像的重建方案。然而,其还一般地应用于其他成像模态,例如但不限于,单光子发射计算机断层摄影(SPECT)和正电子发射断层摄影(PET)。
背景技术
在诊断核成像中,在放射性核素通过患者的血流时,针对对循环系统进行成像或者针对对积聚所注射放射性药物的特定器官进行成像,而对放射性进行研究。有利地,所述放射性药物可以被设计为集中于所选组织中,以提供对那些所选组织的优选成像。
在单光子发射计算机断层摄影(SPECT)中,使用一个或多个通常被称作伽马相机的辐射探测器,以通过由放射性衰变事件造成的辐射发射,来探测放射性药物。通常,每个伽马相机包括辐射探测器阵列和设置在所述辐射探测器阵列前方的准直器。所述准直器限定线形或小角度锥形视线,使得探测到的辐射包括投影数据。如果在角度视场范围,例如180°或360°角度范围上移动所述伽马相机,则得到的投影数据可以被重建为患者中的放射性药物分布的图像。
在正电子发射断层摄影(PET)中,放射性药物的放射性衰变事件产生正电子。每个正电子与电子相互作用,以产生正电子-电子湮灭事件,其发射两个反向的伽马射线。使用符合探测电路,围绕成像患者的辐射探测器的环形阵列探测与所述正电子-电子湮灭相对应的符合反向伽马射线事件。连接两个符合探测的响应线(LOR)包括正电子-电子湮灭事件的位置。这种响应线可以被重建,以产生所述放射性药物分布的图像。
在飞行时间PET(TOF-PET)中,所述两个符合γ射线事件的探测时间之间小的时间差被用于沿所述LOR(响应线)定位所述湮灭事件。
在计算机断层摄影(CT)成像中,辐射源辐照成像受试者;并且设置在所述成像受试者相对侧的辐射探测器阵列探测所传输的辐射。归因于由所述成像受试者中的组织对辐射的不同衰减,所探测的辐射可以被重建为描绘所述成像受试者中的辐射吸收结构的图像。
针对核发射断层摄影和X射线透射型断层摄影的重建算法包括滤波后向投影(FBP)或卷积后向投影法,以及重建发射和透射投影数据所必须的相关联的滤波。FBP是采用到达探测器的信号(signal-to-detector)的方案的分析方法。FBP容易实现、计算快速,并且一般为线性的。然而,FBP有一些限制,值得一提的的是未提供对源自于探测器的低光子计数的噪声的建模,这抑制了信噪比(SNR)。降低使用FBP重建的图像噪声的一种方式是使用更高的辐射剂量,但这是关乎患者安全的问题。
另一类重建技术是使用复杂的迭代算法的迭代重建方法,所述迭代算法细化并建立图像表示,并且能够提供针对衰减、噪声、散射、探测器效率、空载时间、解剖模型、患者运动,以及放射性示踪剂动力学,等等的校正。
然而,迭代重建方法容易有截断误差。原理上,迭代重建方法通过将对重建图像的模拟采集与所测量的数据进行比较而进行。例如,在CT的情况中,所述方法计算所重建的衰减图像的线积分,并将它们与所测量的衰减进行比较。根据模拟数据与测量数据之间的差,以迭代方式细化所述重建图像,直到所述模拟数据匹配所述测量数据。在未囊括完整成像目标的截断重建FoV(重建体积)的情况中,所述测量数据可以包含最初位于所述重建FoV外部的衰减。然而,所述迭代方法,在给定固定的FoV在匹配模拟的和测量的衰减中,将该衰减放在所述FoV内部。这可能导致强的伪影,例如拖尾和阴影。截断问题对于成像目标的横截轴(trans-axial)截断尤其严重。增大重建FoV,以囊括完整的成像目标,去除截断伪像。类似的论证应用于发射断层摄影的情况,在发射断层摄影中重建示踪剂放射性,而非衰减。
有关SPECT和PET成像技术两者的一个问题在于,由放射性核素与探测器之间的患者的解剖结构的光子吸收和散射扭曲了得到的图像。一种技术使用X射线CT扫描数据来生成衰减图。由于X射线和伽马射线两者均被硬组织(例如,骨或甚至是合成植入物)更强地衰减(与软组织相比),所述CT数据可以被用于估计针对由放射性药物发射的伽马射线的衰减图。如果所重建的衰减图遭遇截断误差的问题,所述误差将被带到经校正的SPECT/PET投影数据中,造成不准确的重建放射性值。
本申请提供新的改进的系统和方法,其克服了上文提到的问题以及其他问题。
发明内容
根据一个方面,提出了一种用于诊断成像的方法。所述方法包括使用第一成像模态采集处于检查区域中的受试者的投影数据;定义视场(FoV),所述视场具有横截轴方向的体素网格;确定所述受试者的最大横截轴跨度;通过将所述FoV的所述体素网格扩展到所述FoV外部的至少一个扩展区域,生成扩展FoV,所述扩展区域至少囊括所确定的最大横截轴跨度和所述横截轴方向中的衰减;以及将所采集的投影数据迭代地重建为所述扩展FoV的图像表示。
根据另一方面,提出了一种成像扫描器。所述成像扫描器包括第一成像扫描器,所述第一成像扫描器使用第一成像模态采集处于检查区域中的受试者的投影数据。至少一个处理器被编程为执行以下步骤:定义视场(FoV),所述视场具有在横截轴方向中的体素网格;确定所述受试者的最大横截轴跨度;通过将所述FoV的所述体素网格扩展到所述FoV外部的至少一个扩展区域,生成扩展FoV,所述扩展区域至少囊括所确定的最大横截轴跨度和所述横截轴方向中的衰减;以及将所采集的投影数据迭代地重建为所述扩展FoV的图像表示。
根据另一方面,提出了一种多模态成像系统。所述系统包括CT扫描器和核成像扫描器,所述CT扫描器采集检查区域中受试者的CT投影数据,所述核成像扫描器采集所述检查区域中的受试者的核投影数据。至少一个CT重建处理器被编程为执行以下步骤:在所述CT投影数据中定义视场(FoV),所述视场具有在横截轴方向中的体素网格;确定所述受试者的最大横截轴跨度;通过将所述FoV的所述体素网格扩展到所述FoV外部的至少一个扩展区域,生成扩展FoV,所述扩展区域至少囊括所确定的最大横截轴跨度和所述横截轴方向中的衰减;将所述CT投影数据迭代地重建为所述扩展FoV的衰减图;以及基于所述迭代地重建的衰减图校正所述和核投影数据。
一个优点在于视场截断误差得以减小。
另一优点在于可以更准确地针对衰减和散射效应校正发射断层摄影数据。
本领域普通技术人员在阅读和理解以下详细描述时将认识到本发明更进一步的优点。
附图说明
本发明可以采取各种部件和部件的布置,以及各种步骤和步骤的安排的形式。附图仅出于图示优选的实施例的目的,并且不应被解释为限制本发明。
图1是组合SPECT/CT单机架系统的图解视图;
图2A是在未在横截轴方向中囊括受试者的视场中的迭代重建的图像表示,其造成强的伪影;
图2B是迭代重建的图像表示,对其以扩展视场进行所述重建以减小截断误差;
图2C是用户选择的视场和扩展视场的图像表示;
图3是用于针对迭代重建方法的截断补偿的方法的流程图;
图4A至图4D是用于确定扩展视场的可选方法的流程图;以及
图5是用于扩展视场的体素网格为扩展视场的方法的流程图。
具体实施方式
参考图1,诊断成像系统10同时地和/或独立地执行计算机断层摄影(CT)和核成像,例如PET或SPECT。所述成像系统10包括限定患者接收膛14的固定壳体12。由所述壳体12支撑的可旋转机架16被布置为绕所述膛旋转,以限定共用检查区域18。纵向和/或垂直地调节支撑要被成像和/或检查的患者或受试者22的患者支撑体20,以实现对所述患者在所述检查区域中的期望定位。
为了提供CT成像能力,安装在所述可旋转机架16上的X射线组件24包括诸如X射线管的X射线源26,以及过滤器、准直器,和/或快门组件28。所述准直器将来自所述X射线源26的辐射准直为锥形射束或楔形射束,一个或多个基本上平行的扇形射束,等等。所述快门门控所述射束打开或关闭。X射线探测器30,例如固态平板探测器,安装在所述可旋转机架16上与所述X射线组件24相对。在图示的实施例中,所述探测器面板相对于辐射的投影中心偏移,或者在横截轴平面中从辐射的中心横向平移。偏移探测器的几何位置是合乎期望的,因为它们允许增大的FoV,或者允许较小的探测器尺寸。较大的探测器倾向与更为复杂、造价昂贵,可能限制总体系统设计,并且可能限制探测器定位或患者接触等等。也预期跨越完整视场的具有宽的、非偏移探测器的X射线系统。
随着所述机架旋转,所述X射线组件24和所述X射线探测器30一起绕所述检查区域18旋转,以跨完整360°回转、多重回转,或较小的弧度,采集CT投影数据。每个CT投影指示沿所述X射线组件24与所述X射线探测器30的探测元件之间的线性路径的X射线衰减。所采集的CT投影数据被存储在CT数据缓存32中,并且由CT重建处理器34处理为CT图像表示,并然后被存储在CT图像存储器单元36中。所述X射线源、所述准直器/快门组件、所述探测器,以及所述重建处理器一起,定义了用于生成解剖学、CT、X射线或第一图像的系统或装置。
为了提供功能性核成像能力,将一个或多个核探测头40a、40b,例如单光子发射断层摄影(SPECT)探测器,以能够移动的方式安装到所述旋转机架16。在同一机架上安装所述X射线组件24和所述核探测头40a、40b允许所述检查区域18通过两种模态被成像,而无需移动所述患者22。在一个实施例中,所述探测头由安装到所述旋转机架16的机器人组件(未示出)以能够移动的方式支撑。所述机器人组件使得能够关于所述患者22以可选的偏移量,例如90°偏移、180°彼此相对等等,定位所述探测头。每个SPECT探测头均包括准直器,使得已知每个所探测的辐射事件已沿可识别的线性或小角度锥形视线发生,以使所采集的辐射包括投影数据。所采集的SPECT投影数据被存储在数据缓存42中,并被SPECT重建处理器44处理成SPECT图像表示,并且被存储在SPECT图像存储器单元46中。所述SPECT探测头和所述SPECT重建处理器一起,定义了用于生成核、功能或第二图像的系统装置。
在另一实施例中,所述功能成像系统或装置包括正电子发射断层摄影(PET)探测器。PET探测器的一个或多个环关于所述患者接收膛14布置,以接收来自其的伽马辐射。探测到的符合辐射事件的对定义了LOR,所述LOR以列表模式存储在数据缓存中,并被PET重建处理器处理为PET图像表示,并且存储在PET图像存储器单元中。所述(一个或多个)PET探测器环和所述PET重建处理器一起,定义了用于生成所述功能图像的系统或装置。
通常采用滤波卷积后向投影(FBP)来重建发射和透射投影数据。FBP是采用到达探测器的信号方案的分析方法,即,其假设来自X射线管上的焦斑的点源辐射作为单一细线通过体素的中心,去往每个探测器元的中心。FBP容易实现并且计算快速。然而,FBP具有一些限制,值得一提的是,没有提供对来源于低光子计数的噪声的建模,这抑制了信噪比(SNR)。为了克服使用FBP重建的图像噪声,可以使用较高的辐射剂量,但这是关乎患者安全的问题。
可选的技术是迭代重建方法,其使用复杂的迭代算法,例如最概然期望最大化法(MLEM)或者有序子集期望最大化(OSEM),其细化并建立图像表示,并且能够提供针对衰减、噪声、散射、探测器效率、空载时间、解剖模型、患者运动,以及放射性示踪剂动力学,等等的校正。在迭代算法中,做出对FoV中的放射性或衰减的初始估计。换言之,例如,从FBP或假设均匀的信号分布,生成初始图像表示。计算投影数据,或从所述初始图像前向投影所述投影数据,并将所述投影数据与相应的采集投影数据进行比较。将所计算的投影数据与实际投影数据之间的差后向投影到所述初始图像表示上,以更新所述图像表示。迭代所述算法,直到所述差最小化,或者最佳方案可用。
尽管迭代重建方法具有吸引人的性质,例如在解决噪声数据时它们能够生成具有与FBP相比较低噪声的重建,但它们易于产生截断误差。
由用户或临床医师选择具有选定的体素网格的FoV,并且确定包括所述FoV外部的扩展区域的扩展FoV(扩展FoV)。为了补偿截断误差,所述CT重建处理器34在重建期间自动地扩展所述视场(FoV),使得其在所述横截轴方向中包括完整的受试者,以计入所述CT投影数据中存在的所有衰减,或者所述PET/SPECT投影数据中存在的所有放射性。
在一个用于确定所述扩展FoV的实施例中,将所述FoV前向投影在所采集的X射线投影数据上,以定义所述FoV“覆盖区(footprint)”。所述横截轴方向中所述FoV覆盖区外部实质的衰减或放射性的量的存在指示截断。然后选择所述扩展FoV足够大,使得其囊括所述横截轴方向中所有投影中存在的所有衰减或放射性。如果所述受试者在所述横截轴方向延伸超过所述探测器,则该方法可能不适用,即所述受试者延伸超过所述扫描器FoV,在该情况下,可以使用以下方法中的一种用于确定所述扩展FoV。
在另一实施例中,从所述受试者的体重和身高结合关于典型患者体型的统计数据,估计所述受试者的横截轴跨度。测量所述受试者的身高和体重,并然后将它们与患者汇集或数据库进行比较,所述患者汇集或数据库根据关于身高和体重的统计数据,把最大横截轴跨度编入索引。然后确定所述扩展FoV,以囊括该估计的患者尺寸。类似地,可以从患者汇集估计典型的最大横截轴跨度。然后将该最大跨度用于确定所述扩展FoV。
在用于确定所述扩展FoV的跨度的另一实施例中,可以从其他可用图像数据估计所述受试者的跨度,例如从之前采集的CT、MRI、SPECT、PET或类似的图像表示。类似地,可以从同时或相继采集的被成像的患者区域没有移出所述检查区域18的投影数据或重建的图像数据,估计所述受试者的跨度。例如,在用所述成像系统10进行SPECT/CT或PET/CT采集时,所述SPECT/PET数据或从其重建的图像表示可以被用于估计所述受试者的轮廓,并且由此估计所述受试者的跨度。或者,可以采集额外的CT投影数据,用于确定所述受试者的跨度。
在用于确定所述扩展FoV的另一实施例中,在非常大的FoV进行第一重建(可能是以低的分辨率),已知所述非常大的FoV囊括完整的受试者。例如,所述FoV可以被选为至少如所述扫描器膛那样大。然后可以在该第一重建中探测真实患者轮廓,并且可以相应地选择针对所扩展的重建的所述扩展FoV。
一旦所述扩展FoV得以确定,则根据所述重建的使用目的和计算硬件限制,以多种方式扩展所述FoV以囊括所述扩展FoV。在一个实施例中,所述FoV的所述体素网格被扩展到所述FoV外部的至少一个扩展区域,直到所述扩展FoV的尺寸。换言之,增加所述FoV的所述体素网格中体素的数目,同时每个体素的尺寸保持不变。在另一实施例中,如果由于处理存储器或处理时间上的限制而使得无法达到所述FoV体素网格的扩展,则增加所述FoV的每个体素的尺寸或维度,同时所述网格中体素的数目保持恒定,由此扩大所述体素网格的空间跨度。
在确定了所述扩展FoV的所述体素网格的位置之后,基于所述扩展FoV,例如使用迭代重建方法,至少生成所述FoV的图像表示,并将所述图像表示存储到CT图像存储器36中,用于在图形用户界面50上显示。所述图形用户界面50还包括用户输入设备,临床医师或用户通过用户输入设备与扫描控制器52交互,以选择扫描序列和协议,等等。
在一个实施例中,所述扩展FoV被显示有所述用户选择的FoV标签。额外地,标签可以指示完全采样的FoV以及可能的欠采样区域。或者,在所述FoV的所述体素网格被扩展的实施例中,可以剪切重建的图像表示,以仅显示用户定义的FoV。在所述FoV的所述体素网格被扩大或扩展并扩大的实施例中,或者所述体素的一些或全部的尺寸被放大的实施例中,重建的图像表示被再次采样,回到所述用户选择的FoV的体素尺寸。
在另一实施例中,重建的扩展FoV被用于生成衰减图,用于校正用所述核成像扫描器采集的投影数据,例如SPECT或PET发射数据。对于精确的衰减校正,重要的是所述衰减图表示整个目标,以抵消采集期间的所有衰减。尽管对于用户定义的FoV的不是这种情况,通过使用所描述的方法以及从所述扩展FoV重建生成衰减图,可以确保完全覆盖。例如,所述SPECT重建处理器44接收所述扩展FoV的重建的图像表示,并且从其生成衰减图。所述衰减图被用于校正由所述受试者的解剖结构对光子的吸收和散射。
尽管针对在计算机断层摄影图像重建中的使用描述了利用扩展FoV以补偿截断误差的迭代重建方法,但是所述方法还可应用于其他成像模态中,也可以针对所述其他成像模态执行迭代重建方法。例如,所述扩展FoV迭代重建方法可以被应用于单层CT、PET、SPECT,等等。至于所述成像系统10,所述SPECT重建处理器44可以通过利用单层CT结合迭代重建算法,补偿所采集的SPECT发射数据中的截断误差,以生成SPECT图像表示。
参考图3,提供一种用于迭代重建方法中的截断补偿的方法。从所述检查区域18采集所述受试者的投影数据S100。使用图形用户界面52,临床医师定义或选择视场S102。确定受试者的最大横截轴跨度S104,然后将其用于确定扩展的视场S106,所述扩展视场至少囊括所确定的最大横截轴跨度以及横截轴方向中的所有衰减或放射性。使用所采集的投影数据迭代地重建所述扩展的视场S108。
可以以多种方式显示所重建的图像表示S110。
在CT投影数据为CT透射数据的实施例中,可以将重建的所述扩展FoV的CT图像表示用于生成衰减图S110,所述衰减图被用于校正核投影数据(例如SPECT/PET数据)S112,然后将所述核投影数据,例如迭代地,重建为衰减校正的核图像表示S114。
可以以多种方式确定所述受试者的横截轴跨度。在图4A中所图示的实施例中,所述扩展FoV的位置、所述FoV被前向投影在所采集的X射线投影数据上,以定义所述FoV“覆盖区”S200。所述横截轴方向中所述FoV覆盖区外部实质衰减的量的存在指示截断。所述扩展FoV然后被选择为足够大S202,以使其囊括所述横截轴方向中所有投影中存在的所有衰减。如果所述受试者在所述横截轴方向扩展超过所述探测器,即,所述受试者扩展超过所述扫描器FoV,则该方法可能不适用,在这种情况下,可以使用以下方法中的一种用于确定所述扩展FoV。
参考图4B,在另一实施例中,从所述受试者的体重和身高,结合有关典型患者体型的统计数据,估计他/她的横截轴跨度S304。测量所述受试者的身高和体重S300,并且将其交叉引用到患者汇集或数据库,并且根据有关身高和体重的统计数据估计所述最大横截轴跨度302。然后确定所述扩展FoV304,以囊括该估计的患者尺寸。类似地,可以从患者汇集估计典型的最大横截轴跨度。然后将该最大跨度用于确定所述扩展FoV。
参考图4C,在另一实施例中,从其他可获得的图像数据,例如从之前采集的CT、SPECT、PET、等等图像表示,估计所述受试者的跨度。通过所述重建处理器34、44获得之前的采集S400,以及生成所述受试者的轮廓S402。然后从所述受试者轮廓确定所述最大跨度S404。类似地,可以从同时或相继采集的投影数据,估计所述受试者的跨度,在所述投影数据中,被成像的患者区域保留在所述检查区域18内。
参考图4D,使用低分辨率和/或诸如FBP的快速算法,在所选的非常大的FoV S500中执行第一重建S502,已知所述FoV S500囊括完整的受试者。例如,所述FoV被选为至少如所述扫描器膛那样大。在该第一重建中确定真实的患者轮廓S504,并且根据所确定的轮廓的最大横截轴跨度,选择针对所扩展的重建的所述扩展FoV。
参考图5,根据所述重建的使用目的和计算硬件限制,可以以多种方式将所述FoV扩展以囊括所述扩展FoV。一在个实施例中,所述FoV的体素网格被扩展到所述FoV外部的至少一个扩展区域,直到扩展FoV的尺寸。换言之,增加所述FoV的体素网格中体素的数目,而每个体素的尺寸保持相同。在另一实施例中,如果所述FoV体素网格的扩展可能由于处理存储器或处理时间而不可达到,则增加所述FoV的每个体素的尺寸或维度,而网格中体素的数目保持恒定,由此扩大所述体素网格的空间跨度。
可以以各种方式实现本文中公开的各种计算和控制部件,例如通过计算机或其他设备,包括被编程的数字处理器,或者包括执行所公开的处理的固件,或者通过被配置为执行所公开的处理或其部分的混合或模拟电路,等等。在一些实施例中,所述计算部件可以通过具有适当固件或编程的计算机来实现。本文中所公开的技术可以通过处理器或其他硬件得以实现,和/或可以被实现为存储指令的存储介质,所述指令在被这种处理器或其他硬件运行时,执行所公开的方法。这样的存储介质可以通过一种或多种类型的存储介质来实现,例如以下的一种或多种:磁盘、光盘、闪存或其他静电存储器;随机存取存储器(RAM);只读存储器(ROM);DVD;等等。
己参考优选的实施例描述了本发明。他人在阅读和理解前面的详细描述时可以想到多种修改和变型。本发明旨在被解释为包含所有这样的修改和变型,只要它们落入权利要求书或其等价方案的范围之内。

Claims (20)

1.一种成像方法,包括:
利用第一成像模态,采集处于检查区域中的受试者的投影数据;
定义视场(FoV),所述视场具有横截轴方向的体素网格;
确定所述受试者的最大横截轴跨度;
通过将所述FoV的所述体素网格扩展到所述FoV外部的至少一个扩展区域,来生成扩展FoV,所述扩展区域至少囊括所确定的最大横截轴跨度和所述横截轴方向上的衰减;以及
将采集的投影数据重建为所述扩展FoV的图像表示。
2.根据权利要求1所述的方法,其中,所述重建步骤包括将所述采集的投影数据迭代地重建为所述图像表示。
3.根据权利要求1和2中任一项所述的方法,还包括:
基于重建的所述扩展FoV的图像表示,生成衰减图;
生成利用第二成像模态采集的所述受试者的第二投影数据集;以及
重建所述第二投影图像数据,包括使用所述衰减图。
4.根据权利要求1和2中任一项所述的方法,其中,与所述第二成像模态的所述投影数据同时地或在时间上相近地,采集来自所述第一成像模态的所述投影数据。
5.根据权利要求1至4中任一项所述的方法,其中,根据所述受试者的体重、所述受试者的身高和患者汇集数据库中的至少一个,估计所述受试者的跨度。
6.根据权利要求1至4中任一项所述的方法,其中,从之前采集的所述受试者的图像,确定所述受试者的跨度。
7.根据权利要求1至4中任一项所述的方法,其中,确定所述受试者的跨度包括:
在已知大于所述受试者的横截轴跨度的FoV中定义粗略体素网格;
将大的FoV重建为粗略图像表示,所述粗略图像表示的分辨率比所述重建的图像表示的分辨率更低;并且
基于所述粗略图像表示,确定所述受试者的横截轴跨度。
8.根据权利要求1至4中任一项所述的方法,其中,确定所述扩展FoV包括:
将所述FoV前向投影到所述采集的投影数据上,以定义FoV覆盖区;并且
扩展所述FoV覆盖区,使得其包括所述投影数据中存在的所有衰减或放射性。
9.根据权利要求1至8中任一项所述的方法,其中,扩展所述FoV的步骤包括以下中的至少一个:
增加所述FoV中体素的数目,以囊括所述至少一个扩展区域;或者
增大所述FoV中每个体素的尺寸,以囊括所述至少一个扩展区域。
10.根据权利要求9所述的方法,其中,所述FoV的所述体素和所述至少一个扩展区域的所述体素具有不同的尺寸。
11.根据权利要求9所述的方法,其中,所述重建的图像表示包括所述FoV和所述扩展区域,所述方法还包括:
对所述重建的图像表示重新采样,以将所述体素尺寸减小到所述FoV体素网格的尺寸;以及
显示重新采样的图像表示,使得定义的FoV可见。
12.一种计算机介质,其承载有软件以控制一个或多个处理器执行根据权利要求1至11中任一项所述的方法。
13.一种处理器,其被编程为执行根据权利要求1至11中任一项所述的方法。
14.一种成像扫描器(10),包括:
第一成像扫描器(24、30、40a、40b),其使用第一成像模态采集处于检查区域(18)中的受试者的投影数据;
至少一个处理器(34、44),其被编程为执行根据权利要求1至11中任一项所述的方法。
15.根据权利要求14所述的成像扫描器,其中,所述第一成像扫描器包括锥形射束计算机断层摄影扫描器(24、30),所述锥形射束计算机断层摄影扫描器(24、30)包括偏移探测器(30)。
16.根据权利要求15所述的成像扫描器(10),还包括:
核成像扫描器(24、30、40a、40b),其采集处于所述检查区域(18)中的所述受试者的核数据;并且
其中,所述处理器还被编程为:
基于迭代地重建的所述扩展FoV的图像表示,生成衰减图;并且使用所述衰减图重建所述核数据。
17.根据权利要求14至16中任一项所述的成像扫描器(10),其中,所述处理器被编程为从以下中的至少一项,确定所述受试者的最大横截轴:
从所述受试者的体重、所述受试者的身高或者来自之前采集的图像的患者汇集数据库;或者
通过定义已知大于所述受试者的横截轴跨度的粗略FoV,将所述粗略FoV重建为粗略图像表示,所述粗略图像表示具有比所述迭代地重建的图像表示更低的分辨率,并且基于所述粗略图像表示确定所述受试者的横截轴跨度;或者
将所述FoV前向投影到采集的投影数据上,以定义覆盖区,并且扩展所述覆盖区,以包括所述采集的投影数据中存在的所有衰减或放射性。
18.根据权利要求14至17中任一项所述的成像扫描器(10),其中,所述处理器被编程为通过以下中的至少一项,扩展所述FoV:
增加所述FoV中体素的数目,以囊括所述至少一个扩展区域;以及
增大所述FoV中每个体素的尺寸,以囊括所述至少一个扩展区域。
19.一种成像扫描器,包括:
用于定义视场(FoV)的器件,所述视场具有横截轴方向的体素网格;
用于确定所述受试者的最大横截轴跨度的器件;
用于通过将所述FoV的所述体素网格扩展到所述FoV外部的至少一个扩展区域来生成扩展FoV的器件,所述扩展区域至少囊括所确定的最大横截轴跨度以及所述横截轴方向中的所有衰减;以及
用于将采集的投影数据迭代地重建为所述扩展FoV的图像表示的器件。
20.一种多模态成像系统(10),包括:
扩展视场计算机断层摄影扫描器(24、30),其包括采集检查区域(18)中受试者的CT投影数据的偏移探测器(30);
核成像扫描器(40a、40b),其采集所述检查区域中所述受试者的核投影数据;
至少一个CT重建处理器(34),其被编程为执行以下步骤:
在所述CT投影数据中定义视场(FoV),所述视场具有横截轴方向的体素网格;
确定所述受试者的最大横截轴跨度;
通过将所述FoV的所述体素网格扩展到所述FoV外部的至少一个扩展区域,来生成扩展FoV,所述扩展区域至少囊括所确定的最大横截轴跨度以及所述横截轴方向中的所有衰减;
将所述CT投影数据,例如迭代地,重建为所述扩展FoV的衰减图;以及
至少一个核重建处理器(44),其被编程为基于迭代地重建的衰减图,校正采集的所述核投影数据。
CN201280006458.7A 2011-01-27 2012-01-18 针对spect/ct系统的迭代锥形射束ct重建的截断补偿 Active CN103329168B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161436714P 2011-01-27 2011-01-27
US61/436,714 2011-01-27
PCT/IB2012/050239 WO2012101548A2 (en) 2011-01-27 2012-01-18 Truncation compensation for iterative cone-beam ct reconstruction for spect/ct systems

Publications (2)

Publication Number Publication Date
CN103329168A true CN103329168A (zh) 2013-09-25
CN103329168B CN103329168B (zh) 2016-10-26

Family

ID=45607309

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280006458.7A Active CN103329168B (zh) 2011-01-27 2012-01-18 针对spect/ct系统的迭代锥形射束ct重建的截断补偿

Country Status (6)

Country Link
US (1) US9462988B2 (zh)
EP (1) EP2668639B1 (zh)
JP (1) JP6192542B2 (zh)
CN (1) CN103329168B (zh)
RU (1) RU2606561C2 (zh)
WO (1) WO2012101548A2 (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105578963A (zh) * 2013-09-27 2016-05-11 皇家飞利浦有限公司 用于组织剂量估计的图像数据z轴覆盖范围延伸
CN106296765A (zh) * 2016-08-02 2017-01-04 上海联影医疗科技有限公司 图像重建方法及系统
CN108447553A (zh) * 2018-03-07 2018-08-24 上海联影医疗科技有限公司 一种医学成像方法、装置、系统及计算机可读存储介质
CN109961491A (zh) * 2019-04-12 2019-07-02 上海联影医疗科技有限公司 多模态图像截断补偿方法、装置、计算机设备和介质
CN110215228A (zh) * 2019-06-11 2019-09-10 上海联影医疗科技有限公司 Pet重建衰减校正方法、系统、可读存储介质和设备
CN110693513A (zh) * 2018-10-27 2020-01-17 上海联影医疗科技有限公司 多模态医学系统的控制方法、系统及存储介质

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8977026B2 (en) * 2012-05-30 2015-03-10 General Electric Company Methods and systems for locating a region of interest in an object
US10846860B2 (en) * 2013-03-05 2020-11-24 Nview Medical Inc. Systems and methods for x-ray tomosynthesis image reconstruction
WO2014192831A1 (ja) * 2013-05-28 2014-12-04 株式会社東芝 医用画像診断装置及び制御方法
US9600913B2 (en) 2013-07-31 2017-03-21 Koninklijke Philips N.V. Iterative CT image reconstruction of a ROI with objects outside the scan FoV
US9980693B2 (en) 2014-01-23 2018-05-29 Ermi, Inc. Imaging proxy
US10188358B2 (en) * 2014-05-15 2019-01-29 General Electric Company System and method for subject shape estimation
US9763631B2 (en) 2014-09-17 2017-09-19 General Electric Company Systems and methods for imaging plural axial locations
KR102379067B1 (ko) 2014-12-01 2022-03-25 삼성전자주식회사 의료 영상 장치 및 의료 영상 처리 방법
US9936926B2 (en) * 2015-02-02 2018-04-10 Palodex Group Oy System and method of small field of view X-ray imaging
US11308662B2 (en) 2016-08-02 2022-04-19 Shanghai United Imaging Healthcare Co., Ltd. System and method for image reconstruction
US10347014B2 (en) 2016-08-02 2019-07-09 Shanghai United Imaging Healthcare Co., Ltd. System and method for image reconstruction
US10395353B2 (en) * 2016-08-31 2019-08-27 Siemens Medical Solutions Usa, Inc. Model-based scatter in multi-modality multi-energy SPECT reconstruction
JP2018105734A (ja) * 2016-12-27 2018-07-05 キヤノンメディカルシステムズ株式会社 核医学診断装置
CN111373448B (zh) 2017-09-22 2023-12-08 尼维医疗公司 使用机器学习正则化器的图像重建
US20210012546A1 (en) * 2018-03-26 2021-01-14 Koninklijke Philips N.V. Automatic fault detection in hybrid imaging
EP3886703B1 (en) * 2018-11-30 2023-12-20 Accuray, Inc. Method and apparatus for improving scatter estimation and correction in imaging
US11357467B2 (en) 2018-11-30 2022-06-14 Accuray, Inc. Multi-pass computed tomography scans for improved workflow and performance
EP3998106A1 (en) 2018-12-17 2022-05-18 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for determining a region of interest of a subject
CN113081010A (zh) * 2021-04-12 2021-07-09 北京肿瘤医院(北京大学肿瘤医院) 基于待测物体体表引导的ct成像系统、方法及其应用
US11647975B2 (en) 2021-06-04 2023-05-16 Accuray, Inc. Radiotherapy apparatus and methods for treatment and imaging using hybrid MeV-keV, multi-energy data acquisition for enhanced imaging
US11605186B2 (en) 2021-06-30 2023-03-14 Accuray, Inc. Anchored kernel scatter estimate
US11794039B2 (en) 2021-07-13 2023-10-24 Accuray, Inc. Multimodal radiation apparatus and methods
US11854123B2 (en) 2021-07-23 2023-12-26 Accuray, Inc. Sparse background measurement and correction for improving imaging

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1504960A (zh) * 2002-10-04 2004-06-16 GEҽҩϵͳ����Ƽ���˾ 用于截断补偿的方法和装置
US6907102B1 (en) * 2002-12-16 2005-06-14 Ken Sauer Iterative reconstruction methods for multi-slice computed tomography
US20050267348A1 (en) * 2004-06-01 2005-12-01 Wollenweber Scott D Methods and apparatus for automatic protocol selection
US20070116344A1 (en) * 2005-11-23 2007-05-24 General Electric Company Method and apparatus for field-of-view expansion of volumetric CT imaging

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2071725C1 (ru) * 1993-06-22 1997-01-20 Научно-производственная коммерческая фирма "Ренси Лтд." Вычислительный томограф (варианты)
EP0747728B1 (en) * 1995-05-11 2003-01-02 Koninklijke Philips Electronics N.V. Improved gamma camera imaging system
JP4170449B2 (ja) 1998-07-07 2008-10-22 株式会社東芝 トランスミッションctのトランケーション補正装置、核医学診断装置及びトランケーション補正方法
JP2002209880A (ja) 2000-12-26 2002-07-30 Ge Medical Systems Global Technology Co Llc X線ctシステム、並びに、操作コンソール及びその制御方法及び記憶媒体
JP2004073432A (ja) 2002-08-15 2004-03-11 Ge Medical Systems Global Technology Co Llc X線ct装置および画像処理方法、ならびにプログラム
CN100573588C (zh) 2004-04-21 2009-12-23 皇家飞利浦电子股份有限公司 使用截短的投影和在先采集的3d ct图像的锥形束ct设备
EP1828808B1 (en) 2004-12-17 2010-04-21 Koninklijke Philips Electronics N.V. Truncation compensation algorithm for iterative reconstruction
JP2006312027A (ja) 2005-04-05 2006-11-16 Toshiba Corp 放射線診断装置
US7502440B2 (en) 2005-04-05 2009-03-10 Kabushiki Toshiba Radiodiagnostic apparatus
DE102006014630B4 (de) 2006-03-29 2014-04-24 Siemens Aktiengesellschaft Verfahren zur Korrektur von Trunkierungsartefakten
JP5337031B2 (ja) 2006-08-17 2013-11-06 コーニンクレッカ フィリップス エヌ ヴェ Ct画像取得
EP2067055B1 (en) * 2006-09-21 2016-03-23 Koninklijke Philips N.V. Cardiac spect system with trajectory optimization
US7737406B2 (en) 2006-09-27 2010-06-15 Siemens Medical Solutions Usa, Inc. Compensating for truncated CT images for use as attenuation maps in emission tomography
DE102006062277A1 (de) 2006-12-22 2008-06-26 Siemens Ag Verfahren zur Korrektur von Trunkierungsartefakten in einem Rekonstruktionsverfahren für computertomographische Aufnahmen
US8184889B2 (en) * 2007-09-28 2012-05-22 Siemens Medical Solutions Usa, Inc. Reconstruction support regions for improving the performance of iterative SPECT reconstruction techniques
US8600139B2 (en) * 2009-06-08 2013-12-03 Koninklijke Philips N.V. Time-of-flight positron emission tomography reconstruction using image content generated event-by-event based on time-of-flight information

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1504960A (zh) * 2002-10-04 2004-06-16 GEҽҩϵͳ����Ƽ���˾ 用于截断补偿的方法和装置
US6907102B1 (en) * 2002-12-16 2005-06-14 Ken Sauer Iterative reconstruction methods for multi-slice computed tomography
US20050267348A1 (en) * 2004-06-01 2005-12-01 Wollenweber Scott D Methods and apparatus for automatic protocol selection
US20070116344A1 (en) * 2005-11-23 2007-05-24 General Electric Company Method and apparatus for field-of-view expansion of volumetric CT imaging

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
E HANSIS等: "《High-quality 3-D coronary artery imaging on an interventional C-arm x-ray system》", 《MEDICAL PHYSICS》 *
HANSIS, E等: "《Nuclear Science Symposium Conference Record (NSS/MIC), 2010 IEEE》", 30 October 2010 *
MAASS, C等: "《Nuclear Science Symposium Conference Record (NSS/MIC), 2010 IEEE》", 30 October 2010, article "《Simple ROI cone-beam computed tomography》" *
WIEGERT J等: "《Projection Extension for Region Of Interest Imaging in Cone-Beam CT》", 《ACADEMIC RADIOLOGY》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105578963A (zh) * 2013-09-27 2016-05-11 皇家飞利浦有限公司 用于组织剂量估计的图像数据z轴覆盖范围延伸
CN105578963B (zh) * 2013-09-27 2019-07-12 皇家飞利浦有限公司 用于组织剂量估计的图像数据z轴覆盖范围延伸
CN106296765A (zh) * 2016-08-02 2017-01-04 上海联影医疗科技有限公司 图像重建方法及系统
CN106296765B (zh) * 2016-08-02 2020-06-02 上海联影医疗科技有限公司 图像重建方法及系统
CN108447553A (zh) * 2018-03-07 2018-08-24 上海联影医疗科技有限公司 一种医学成像方法、装置、系统及计算机可读存储介质
CN108447553B (zh) * 2018-03-07 2021-07-30 上海联影医疗科技股份有限公司 一种医学成像方法、装置、系统及计算机可读存储介质
CN110693513A (zh) * 2018-10-27 2020-01-17 上海联影医疗科技有限公司 多模态医学系统的控制方法、系统及存储介质
CN109961491A (zh) * 2019-04-12 2019-07-02 上海联影医疗科技有限公司 多模态图像截断补偿方法、装置、计算机设备和介质
CN110215228A (zh) * 2019-06-11 2019-09-10 上海联影医疗科技有限公司 Pet重建衰减校正方法、系统、可读存储介质和设备
CN110215228B (zh) * 2019-06-11 2023-09-05 上海联影医疗科技股份有限公司 Pet重建衰减校正方法、系统、可读存储介质和设备

Also Published As

Publication number Publication date
WO2012101548A3 (en) 2013-01-03
WO2012101548A2 (en) 2012-08-02
US9462988B2 (en) 2016-10-11
EP2668639B1 (en) 2016-11-02
JP6192542B2 (ja) 2017-09-06
RU2013139541A (ru) 2015-03-10
CN103329168B (zh) 2016-10-26
EP2668639A2 (en) 2013-12-04
RU2606561C2 (ru) 2017-01-10
JP2014507988A (ja) 2014-04-03
US20130294570A1 (en) 2013-11-07

Similar Documents

Publication Publication Date Title
CN103329168A (zh) 针对spect/ct系统的迭代锥形射束ct重建的截断补偿
CN101849247B (zh) 用于生成衰减图的设备和方法
CN107111867B (zh) 多模态成像系统及方法
CN102763138B (zh) 辐射治疗中的运动校正
US8008625B2 (en) Method and apparatus for high-sensitivity single-photon emission computed tomography
US9747701B2 (en) Systems and methods for emission tomography quantitation
CN101080651B (zh) 成像系统及成像方法
CN101842806B (zh) 脏同位素pet重建
Wu et al. Absolute quantitative total-body small-animal SPECT with focusing pinholes
CN102982510A (zh) Spect图像的运动校正
CN103315760A (zh) 用于核医疗成像中基于发射数据进行衰减补偿的系统和方法
CN106659452B (zh) 在定量单光子发射计算机断层扫描中利用多个光电峰的重构
CN101646957B (zh) 各向同性分辨率的图像重建
HU231327B1 (hu) Többemissziós energiák egyfotonos emissziós komputertomográfiában
Rahman et al. Fisher information analysis of list-mode SPECT emission data for joint estimation of activity and attenuation distribution
Ghani et al. An empirical method for geometric calibration of a photon counting detector-based cone beam CT system
Rahman Development and Objective Task-Based Evaluation of Computational Methods for Clinical SPECT
Eclov On the Feasibility of Positron Emission Tomography Imaging Systems Utilizing Mixed Spatial Resolution Detectors
Francis 3D Volumetric Reconstruction for Light-Field SPECT
JP2023141790A (ja) 核医学診断装置および吸収係数画像推定方法
Du et al. Principles of brain single-photon emission computed tomography imaging
Pointoin Model-based randoms correction for 3D positron emission tomography
NUYTS et al. HIGH QUALITY IMAGE FORMATION OF SMALL ANIMALS USING EMISSION TOMOGRAPHY
Kolthammer Quantitative Positron Emission Tomography for Estimation of Absolute Myocardial Blood Flow

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant