CN111183371B - 处理数字正电子发射断层摄影的探测器像素性能变化 - Google Patents

处理数字正电子发射断层摄影的探测器像素性能变化 Download PDF

Info

Publication number
CN111183371B
CN111183371B CN201880061097.3A CN201880061097A CN111183371B CN 111183371 B CN111183371 B CN 111183371B CN 201880061097 A CN201880061097 A CN 201880061097A CN 111183371 B CN111183371 B CN 111183371B
Authority
CN
China
Prior art keywords
detector
pixels
pixelated
dead
detector pixels
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201880061097.3A
Other languages
English (en)
Other versions
CN111183371A (zh
Inventor
白传勇
A·安德烈耶夫
S·林
张滨
M·A·米勒
宋犀云
叶京汉
S·德维韦迪
胡志强
Y-L·谢
I·布罗德斯基
T·C·布尔格林
Y-M·朱
D·B·麦克奈特
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 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 NV filed Critical Koninklijke Philips NV
Publication of CN111183371A publication Critical patent/CN111183371A/zh
Application granted granted Critical
Publication of CN111183371B publication Critical patent/CN111183371B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T7/00—Details of radiation-measuring instruments
    • G01T7/005—Details of radiation-measuring instruments calibration techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G01T1/2914—Measurement of spatial distribution of radiation
    • G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
    • 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/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
    • A61B6/582—Calibration
    • A61B6/585—Calibration of detector units
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Optics & Photonics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

一种非瞬态计算机可读介质存储指令，所述指令能由包括至少一个电子处理器(20)的工作站(18)读取并运行以执行质量控制(QC)方法(100)。所述方法包括：接收由像素化探测器(14)采集的当前QC数据集和由所述像素化探测器采集的一个或多个先前QC数据集；根据所述当前QC数据集和所述一个或多个先前QC数据集来确定所述像素化探测器的探测器像素(16)随着时间的稳定性水平；当针对所述探测器像素确定的所述稳定性水平在稳定性阈值范围之外时，将所述像素化探测器的探测器像素标记为死的；以及在与所述工作站可操作地连接的显示设备(24)上显示对被标记为死的所述探测器像素的识别(28)。

Description

处理数字正电子发射断层摄影的探测器像素性能变化
技术领域
以下大致涉及医学成像领域、医学图像解读领域、图像重建领域、医学成像设备维护领域、以及相关的领域。
背景技术
数字正电子发射断层摄影(PET)探测器包括大量像素的组件。在一种像素化探测器设计中，每个探测器像素是被切成期望尺寸的小闪烁晶体，并且具有相关联的闪烁光探测单元和电子设备以在PET扫描中探测511keV伽马射线。晶体制备、探测器组装过程等被保持为尽可能系统性的，使得大多数像素的性能是相对类似的且可预测的，由此形成平均像素的主要子集。
探测器像素的小部分可以由于许多原因(诸如晶体不均匀性、制造加工变化、光子探测单元(例如，光电二极管)性能波动、组装工艺不一致、电子设备变化等)而显著不同于平均像素。这样的像素中的一些具有比平均像素低得多的灵敏度。这样的像素中的一些具有比平均像素高得多的灵敏度；并且还有一些像素可能具有不稳定的灵敏度。即，灵敏度不时急剧地或连续地变化。
对探测器像素性能的评价通常通过探测器像素校准来执行，所述探测器像素校准涉及从标准均匀体模采集事件数据，所述标准均匀体模优选地被放置在扫描器等中心处以便到探测器环的探测器的距离相等。由于所涉及的时间和努力，像素校准通常很少(例如仅在主要维护等之后)被执行。校准通常还包括确定用于探测器像素的规范化因子以考虑个体像素的灵敏度的差异。
在一些现有PET系统中，具有极其低的灵敏度(例如，为平均像素的灵敏度的20％或更小的下限)的像素和具有不自然地高的灵敏度(例如，高于平均像素40％的上限)的像素被视为是死像素。当像素被识别为死的时，它们从数据处理中排除，好像它们不存在。例如，与死像素相关联的事件从系统性能评估(诸如美国电气制造商协会(NEMA)灵敏度、计数率性能评估)中排除，以及从图像重建过程(例如列表模式迭代重建)中排除。在重建中，用于处理死像素的方法有时用来最小化死像素在图像质量(例如，分辨率、伪影等)和定量准确性(例如，病变强度、SUV等)中的负面影响。
用于处理非典型探测器像素的这些方法具有一些缺点。例如，具有随着时间在上限与下限之间波动的灵敏度的不稳定像素可能难以识别。如果像素具有在例如平均像素的20％与140％之间变化的灵敏度，那么对PET系统的校准可能不能探测到这样的像素，因为校准很少被执行并且不以高时间分辨率跟踪性能。规范化可以隐含地补偿这种变化，但是发生在两个规范化之间的任何像素性能变化不能被探测和补偿。
如果一组聚集像素(例如，瓦片、模块等)具有非常低的灵敏度(例如，平均的50％)或其灵敏度与其他像素大不相同，那么它们可能在某些情况下引入显著的图像伪影和定量误差。在灵敏度变化的根本原因处于瓦片或模块水平(例如影响瓦片或模块的所有探测器的关于瓦片或模块电子设备的问题)的情况下，这种情况会出现。规范化过程可以隐含地补偿这种瓦片或模块水平变化，但是规范化很少被执行。如果瓦片或模块的灵敏度变化发生在两个规范化之间，那么患者数据可能受损，并且得到的图像可能具有伪影和定量误差。
将像素标记为死的现有方法也可能是涵盖过广的，致使排除有用的成像数据。如果被标记为死的像素在扫描的持续时间内是稳定的，那么与这些死像素相关联的计数仍然可以是有用的。实际上，这样的计数在低剂量研究、短扫描或动态扫描中可以是非常有价值的。然而，现有方法从被使用中排除来自死像素的计数。
当死像素的数量增加和/或他们变得聚集时，用户可能担心系统使用起来是否仍然是安全的。服务呼叫的指导方针对那些熟练的用户来说可能是不足的，可能导致不必要或过早的维护呼叫。相反，注意力不集中的用户可能不能识别死像素的数量正在变得高到足以不利地影响临床成像质量。
以下公开了克服这些问题的新的且改善的系统。
发明内容
在一个公开的方面中，一种非瞬态计算机可读介质存储指令，所述指令能由包括至少一个电子处理器的工作站读取并运行以执行质量控制(QC)方法。所述方法包括：接收由像素化探测器采集的当前QC数据集和由所述像素化探测器采集的一个或多个先前QC数据集；根据所述当前QC数据集和所述一个或多个先前QC数据集来确定所述像素化探测器的探测器像素随着时间的稳定性水平；当针对所述探测器像素确定的所述稳定性水平在稳定性阈值范围之外时，将所述像素化探测器的探测器像素标记为死的；以及在与所述工作站可操作地连接的显示设备上显示对被标记为死的所述探测器像素的识别。
在另一公开的方面中，一种成像系统包括图像采集设备和显示设备，所述图像采集设备配置为采集患者的成像数据，所述图像采集设备包括具有多个探测器像素的像素化探测器。至少一个电子处理器被编程为：接收由所述像素化探测器采集的当前质量控制(QC)数据集和由所述像素化探测器采集的一个或多个先前QC数据集；根据所述当前QC数据集和所述一个或多个先前QC数据集来确定所述探测器像素随着时间的稳定性水平；当针对所述探测器像素确定的所述稳定性水平在稳定性阈值范围之外时，将所述探测器像素中的一个或多个标记为死的；并且控制所述显示设备以显示对被标记为死的所述探测器像素的识别。
在另一公开的方面中，一种成像系统包括：正电子发射断层摄影(PET)设备，其配置为采集患者的成像数据。所述图像采集设备包括具有多个PET探测器像素的PET探测器环。所述成像系统还包括显示设备。至少一个电子处理器被编程为：接收由像素化探测器采集的当前质量控制(QC)数据集和由所述像素化探测器采集的一个或多个先前QC数据集；根据所述当前QC数据集来确定所述探测器像素的灵敏度水平；根据所述当前QC数据集和所述一个或多个先前QC数据集来确定所述探测器像素随着时间的稳定性水平；当所述灵敏度水平在最大灵敏度阈值以上时，将所述像素化探测器的探测器像素标记为死的；当所述灵敏度水平在最小灵敏度阈值以下时，将所述像素化探测器的探测器像素标记为冷的；并且控制所述显示设备以显示对被标记为死的或冷的所述探测器像素的识别。
一个优点在于不稳定探测器像素的更有效且及时的识别。
另一优点在于识别已经重新稳定的之前不稳定的探测器像素。
另一优点在于重建包括使用具有低但稳定的灵敏度水平的探测器像素采集的数据的成像数据。
另一优点在于为临床医生提供对被标记为死的探测器像素的识别。
另一优点在于更有效地通知临床医生死探测器像素对临床图像质量的影响。
另一优点在于减少对于PET探测器环维护问题的维护呼叫。
给定的实施例可以提供前述优点中的零个、一个、两个、更多个或全部，并且/或者可以提供在阅读并理解本公开后对于本领域普通技术人员来说将会变得显而易见的其他优点。
附图说明
本公开可以采取各种部件和部件的布置，以及各种步骤和步骤的安排的形式。附图仅出于说明优选实施例的目的，并且不得被解释为对本发明的限制。
图1图解性地示出了根据一个方面的图像重建系统。
图2示出了图1的系统的示范性流程图操作；以及
图3在图1的系统的显示器上说明性地示出了对探测器像素的识别。
具体实施方式
在现有的数字PET成像设备维护中，探测器像素规范化校准很少被执行，例如在主要维护之后或按照在校准之间具有长(例如多个月)间隔的日程计划。规范化涉及从标准均匀体模采集数据。在规范化校准期间，具有在最小阈值以下的灵敏度的任何像素(“冷”像素)被标记为死的，并且类似地，具有在最大阈值以上的灵敏度的任何像素(“热”像素)也被标记为死像素。根本原因是，冷像素最可能丢失许多计数，而热像素产生许多虚假计数。使用被标记为死像素的像素采集的计数被忽略，并且在迭代图像重建中使用的灵敏度矩阵被调整以考虑丢失的像素。
此外，日常质量控制(QC)程序通常使用标准点源体模(例如22Na点源)来执行。QC检查各种探测器像素参数(诸如能量分辨率、探测器不均匀性)，并且还探测死像素。从QC程序生成的信息不用来调整任何扫描器设置，而是仅仅用来验证扫描器正在可接受的包络内操作。如果例如QC探测到太多的死像素，那么这可能导致维护呼叫。
本文中公开的实施例有利地利用现有QC程序来评价探测器像素的每天稳定性。在一些预见到的实施方式中，该额外的QC信息是信息性的，即死像素基于灵敏度阈值或基于不稳定性来探测，并且如果存在太多的死像素，那么这可以触发维护呼叫或新的规范化校准。
在本文中公开的其他实施例中，从QC数据获得的信息可以用来将不稳定的探测器像素标记为死的，和/或用来将当前被标记为死的像素重新标记为活像素。在一个说明性过程中，标记包括：(1)即使其灵敏度在下限/上限阈值边界内，也将任何不稳定的像素标记为死的；(2)将其灵敏度超过阈值的任何像素标记为死的；以及(3)将其灵敏度低于阈值的任何稳定像素标记为活的。
在另一公开的方面中，新识别的死像素的影响可以被模拟以通知用户这些死像素的实际临床影响。这种模拟能够以简单的方式被完成：现有体模或临床图像与其重建利用由于所移除的新识别的死像素的任何计数被重复(并且在图像重建中使用的灵敏度矩阵被适当地调整)的原始数据集进行比较。除了模拟死像素的现有群体之外，预见到模拟针对额外死像素的预测来完成，以通知用户从临床视角死像素计数何时变得有问题。如果这针对实际临床图像来完成，那么影响是为临床医生提供死像素对在医院处实际执行的类型的临床成像的影响的真实预测。
在一些情况下，采集的数据能够被保存在性能图中。例如，在日常的QC数据收集中，像素性能图能够使用符合计数来获得。对于现有PET系统，为了避免几何响应变化和与真实死像素相反的正常晶体的遮蔽，能够备选地使用每个像素中的单个(即，单个511keV伽马)来获得像素性能图。由于像素性能图每日获得，所以对像素性能的分析也每日执行。这种图在数据处理中的应用能够提供每天采集的数据的可信赖性的保证。因此，即使日内变化不能被探测和处理，风险也被最小化。这种每天过程使每日像素灵敏度规范化校准程序不必要，由此对于每日操作显著地节省时间和费用(例如，体模、源等)。
通过现有PET系统实施的典型图像重建过程使用探测器规范化结果来隐含地处理像素性能变化。然而，规范化不是每日执行，因为由于对使均匀体模充满活动并且执行扫描和分析的需要，它是耗时的过程。在本文中描述的方法和系统中，在相继的规范化之间发生的像素性能变化在每日基础上被探测，并且将得到的性能图应用于图像重建中能够改善数据保真度，并且最小化数据受损的风险。在这些过程中，间歇的(即不稳定的)像素从数据处理中排除。证明随着时间稳定的个体像素灵敏度中的任何调整被包括于迭代重建算法或建模/包括重建中的系统响应的其他算法中的投影/反向投影模型中。
在一些实施例中，提供了用于在系统性能、图像质量和量化上评估像素性能变化的QC工具。QC工具提供了将所测量的像素性能图加载在工作站系统上，并且使用该图来运行关于系统性能、图像质量和量化的典型分析。QC工具还提供了用于手动地修改像素性能图以针对数字PET模拟不同的像素性能变化影响情景。QC工具能够使用模拟的数据或在系统上采集的数据用于系统性能、图像质量和量化分析，包括计数性能、分辨率、对比、噪声、均匀性、SUV等。QC工具为用户提供定性且定量的结果以对照要求评估系统性能，帮助用户理解当一些大像素性能变化被探测到时系统是否仍然适合于使用。
虽然参考数字PET进行描述，但是以下方法可应用于采用具有探测器像素的探测器的任何类型的成像，例如数字PET、透射型计算机断层摄影(CT)、数字单光子发射计算机断层摄影(SPECT)或数字放射摄影(DR)(即，具有平坦像素化的数字x-射线探测器的平板x-射线)。
参考图1，示出了说明性医学成像系统10。如图1中示出的，系统10包括图像采集设备12。在一个范例中，图像采集设备12能够包括发射成像设备(例如，正电子发射断层摄影(PET)设备、用于在单光子发射计算机断层摄影(SPECT)中使用的伽马相机等)；然而，将意识到，任何其他合适的成像模态(例如，透射型计算机断层摄影(CT)、X-射线等、以及混合系统(诸如PET/CT))可以被使用。图像采集设备12包括具有多个探测器像素16(在图1中被示为插图A)的像素化探测器14，所述多个探测器像素被布置为从被设置在检查区域17中的患者收集成像数据。取决于图像采集设备12的模态，像素化探测器14能够是PET设备的探测器环(例如，整个PET探测器环或其一部分，诸如探测器瓦片，探测器模块等等)；CT设备的探测器阵列；被配置为执行SPECT的伽马相机的探测器；以及数字放射摄影设备(诸如X-射线机器)的数字探测器阵列。
系统10还包括具有典型部件的计算机或工作站或其他电子数据处理设备18，诸如至少一个电子处理器20、至少一个用户输入设备(例如，鼠标、键盘、跟踪球等)22和显示设备24。在一些实施例中，显示设备24能够是与计算机18分开的部件。工作站18还能够包括一个或多个数据库26(被存储在非瞬态存储介质(诸如RAM或ROM、磁盘等等)中)，并且/或者工作站能够与一个或多个数据库27(例如，电子医学记录(EMR)数据库、图片归档与通信系统(PACS)数据库等)电子通信。如本文中描述的，数据库27是PACS数据库。
至少一个电子处理器20与存储指令的非瞬态存储介质(未示出)可操作地连接，所述指令能由至少一个电子处理器20读取并运行以执行所公开的操作，包括执行质量控制(QC)方法或过程100。非瞬态存储介质可以例如包括：硬盘驱动器、RAID或其他磁性存储介质；固态驱动器、闪烁驱动器、电可擦除只读存储器(EEROM)或其他电子存储器；光盘或其他光学存储设备；其各种组合；等等。在一些范例中，QC方法或过程100可以通过云处理来执行。与探测器规范化过程相比，QC方法或过程100在相对频繁的基础上被执行。例如，在一些实施例中，QC方法或过程100在每日基础上(例如在成像设备的启动期间在早晨)被执行。在执行QC过程100之后，成像设备然后用来执行对患者的临床成像，每个患者进而被装载到检查区域17中以进行成像并且然后被卸载以准许接纳下一临床患者。
参考图2，QC方法100的说明性实施例被图解性地示为流程图。为了开始该过程，标准点源体模被装载到检查区域17中，优选地在成像区域17的等中心处。为了PET校准，典型的标准点源体模是NEMA-标准22Na点源。在102处，至少一个电子处理器20被编程为接收由像素化探测器14从点源采集的当前质量控制(QC)数据集。QC方法100还具有到在早前的时候(例如在过去几天)由像素化探测器采集的一个或多个先前QC数据集的入口。先前QC数据集可以例如被存储在数据库26中。所采集的点源的QC数据(例如，符合数据或单个数据)指示探测器像素16的个体性能，并且能够在每日基础上被收集。由于点源具有例如以Bq单位测量的标准放射性，由所有探测器采集的计数应当在考虑从点源到探测器像素的距离之后是近似相同的，假设它们具有完全相同的灵敏度。例如，在多环PET探测器中，在图像采集设备12的等中心处的点源导致到不同像素化探测器14处的探测器像素16的不同距离，导致系统性灵敏度差异。同时，由于典型PET探测器环中关于视场的中心和源位置存在许多几何对称性，预期等效定位的探测器像素具有类似的灵敏度。因此，当前QC数据集中的计数的差异指示不同的像素灵敏度。关于探测器像素性能的其他信息也能够被确定。探测器像素16的性能数据能够包括灵敏度、能量分辨率、飞行时间(TOF)漂移(即探测器时间延迟)等中的一个或多个。代替提到的22Na点源，在一些实施例中，从像素化探测器14采集包括线辐射源的参考源的当前QC数据和先前QC数据集。线源能够为多环PET探测器提供更均匀的辐照。一旦被接收，所采集的QC数据(当前和先前两者)能够被保存在数据库26中。
在104处，至少一个电子处理器20被编程为根据当前QC数据集和一个或多个先前QC数据集来确定像素化探测器14的探测器像素16的灵敏度和稳定性水平。为了这样做，探测器像素16的灵敏度水平针对每个QC数据集(例如当前QC数据集和每个先前QC数据集)基于由探测器像素采集的实际计数与(基于点源或线源的反射性和/或所有探测器的平均计数来预期的)预期计数的比来确定。处理器20确定所确定的探测器像素16的灵敏度水平是否在灵敏度阈值范围之外。“热”像素可以被定义为具有高于上限阈值的灵敏度；而“冷”像素可以被定义为具有在下限阈值以下的灵敏度。具有在下限阈值与上限阈值之间的灵敏度的像素被视为是正常像素。
探测器像素稳定性被确定为当前和先前QC数据集中的探测器像素的灵敏度的随着时间的变化。由于QC数据集通常在每日基础上被采集，以“每天”时间间隔确定稳定性(即，根据时间的灵敏度)是简单的。不稳定的探测器像素是其随着时间的灵敏度变化是不可接受地大的探测器像素。例如，识别具有根据时间的大灵敏度变化(例如，一天灵敏度是平均值的10％，而另一天平均值的60％等等)的探测器像素16。应当注意，在如根据当前QC数据集确定的其灵敏度在下限与上限阈值之间的意义上，探测器像素可以被分类为“正常的”，但是如果如根据过去和当前QC数据集确定的其灵敏度每天显著变化(即使该每天变化保持在下限和上限阈值内)，则被分类为“不稳定的”。
在106处，至少一个电子处理器20被编程为，当针对探测器像素确定的灵敏度水平在灵敏度阈值范围之外(例如，具有小于平均灵敏度的20％的下限和具有大于平均灵敏度的140％的上限)时将像素化探测器14的探测器像素16中的至少一个标记为“死的”，并且其余探测器像素被标记为“活的”。例如，当所确定的灵敏度水平在灵敏度阈值范围的最大灵敏度阈值以上时，“热”探测器像素16被标记为死的。当所确定的灵敏度水平在灵敏度阈值范围的最小灵敏度阈值以下时，“冷”探测器像素16被标记为死的。在又一范例中，当之前被标记为死的一个或多个探测器像素的所探测到的稳定性水平在稳定性阈值范围内时，之前被标记为死的一个或多个探测器像素16能够被重新标记为活的。
此外，在操作106中，至少一个电子处理器20被编程为当针对探测器像素确定的稳定性水平在稳定性阈值范围之外时将像素化探测器14的探测器像素16中的至少一个标记物为“死的”。例如，作为非限制性说明性范例，稳定性阈值范围可以是在过去五天内不超过20％的灵敏度的变化。使用该范例，如果探测器像素在过去五天(包括当天)内已经测量到以下灵敏度：70％；65％；60％；62％；72％，那么该像素被视为满足20％稳定性阈值范围。相比之下，在过去五天(包括当天)内已经测量到以下灵敏度：52％，70％；75％；79％；71％的探测器像素将会被视为不稳定的，因为其灵敏性已经在52％与79％之间的范围内变动(27％的范围，大于20％稳定性阈值范围)。该像素将会由于其不稳定性而被标记为死的，即使所测量的灵敏度都落在灵敏度阈值范围(20％-140％)内。备选地，可以存在能够在确定它是死的之前被采取的一些第二遍统计步骤；例如，QC数据中的单个异常值可以被排除，并且/或者第二遍评价可以包括增加所使用的QC数据点的数量、考虑统计方差、其组合等等。
在之前的范例中，如果它不满足灵敏度阈值范围或稳定性阈值范围，则像素被标记为死的。在另一实施例中，稳定性评价用来保留通常将会被标记为死的一些稳定像素。在一个这种范例中，如果它们满足稳定性阈值范围，则冷的探测器像素16被标记为活的，并且被包括在进一步的数据处理操作中。相比之下，在该范例中，即使它们是稳定的，热像素也被标记为死的。该方法的根本原因是，稳定的冷像素仍然提供对辐射探测事件的有用计数；然而，稳定的热像素是不可靠的，因为高灵敏度可能归因于与辐射探测事件不相关的高暗计数率。
冷的但稳定的探测器像素的保留具有实质优点，尤其是在低计数成像的情况下。例如，对于具有非常小的焦距的点源扫描或患者扫描，如果来自这样的冷像素的计数被排除，那么像素填充方法可以被采用以估计与这样的像素相关联的计数。然而，如果填充使用相邻像素的平均或外插，那么像素填充方法可能将误差/偏差引入到数据。在像素是聚集的情况下尤为如此，对于其空间分辨率损失也将会是显著的。使用来自具有对应性能图的冷的(但稳定的)探测器像素16的计数保持数据(包括空间分辨率)的准确性，并且克服与冷像素相关联的稍微更高的噪声水平。
在104、106处对当前和先前QC数据集操作的探测器像素评价的输出可以以各种方式被使用，如接下来描述的。
在108处，在一个应用中，至少一个电子处理器20被编程为控制显示设备22以显示对被标记为死的探测器像素16的识别28。识别能够以任何合适的格式被显示，任何合适的格式诸如列表、探测器图等。在一些范例中，识别28能够显示探测器14的图，其中被标记为死的像素使用填充方框来标识，而活像素使用未填充方框来标识。额外地或备选地，可以提供死像素的列表。在一些实施例中，图可以不同地标识由于在灵敏度阈值范围之外而被标记为死的那些像素与由于在稳定性阈值范围之外而被标记为死的那些像素。备选地，(针对在灵敏度范围之外和针对不稳定的像素)可以示出两个不同的探测器图。虽然这种类型的呈现可以是有用的，但是它具有潜在的缺点，因为临床医生可能具有理解绘制的死像素实际上如何影响临床图像的困难。由于关于死像素的实际临床影响的这种不确定性，临床医生可能在(尚)未需要时预约服务呼叫。相反，如果临床医生低估死像素的实际临床影响，那么当对探测器14进行服务将会是更优选的时，成像系统可能继续用于临床成像。
在110处，在另一应用中，呈现死像素对实际临床图像的影响。为此目的，至少一个电子处理器20被编程为在具有或没有没有死像素的情况下重建采集的或模拟的成像数据，以便展示死像素在临床图像空间中的影响。在一个范例中，处理器20被编程为例如使用模拟来模拟将由像素化探测器14针对成像对象采集的成像数据。两个图像重建被执行。第一个使用所有模拟的成像数据(这里被称为“参考图像”)被执行。参考图像能够根据先前QC数据集中的数据(包括利用之前标记的死像素采集的数据)来生成，图像通过省去之前标记的死像素来重建，并且图像质量和量化已经被视为满意的。第二个使用所有模拟的成像数据来执行，除了被标记为死的探测器像素16被处置为不提供数据(这里被称为“当前QC图像”)。被模拟为已经通过这些死像素采集的数据从进一步数据处理中排除，并且重建过程基于这些排除的像素来调整。得到的当前QC图像和参考图像在显示器中并排或以一些其他容易感知到的可比较格式(例如响应于用户按压按钮而在两幅图像之间切换)被呈现，使得临床医生能够看见死像素对临床图像的影响)。
对于更高级的用户，前述模拟可以替代地针对点体模或线体模来执行，其中模拟的数据再次被重建以形成当前QC和参考图像。点体模或线体模的简化的几何形状可以以PET成像数据采集和图像重建过程的更好理解为用户提供更多信息性的比较图像。
应当注意，成像数据的模拟能够被执行一次，并且被存储在工作站18的数据库26中。同样地，使用所有成像数据的参考图像重建可以被完成一次，并且被存储在工作站18的数据库26中。此后，执行操作110仅仅需要在当前死像素集被省去的情况下使用存储的模拟数据来执行当前QC图像重建。
在108、110处的(一个或多个)呈现可以被动地被呈现给用户。在另一方法中，QC过程100可以额外地对这些结果执行主动诊断，以提供维护建议。因此，在112处，至少一个电子处理器20被编程为基于被标记为死的探测器像素16来确定是指示对像素化探测器14的维护还是校准。例如，如果大量探测器像素16被标记为死的，那么像素化探测器14可能需要被重新校准。如果在校准过程之后大量探测器像素16仍然被标记为死的，那么维护预约能够被请求以更换像素化探测器14。在这种情况下，显示设备20显示根据对是指示对像素化探测器14的维护还是校准的确定来指示执行对像素化探测器14的维护的建议或执行对像素化探测器14的校准的建议的消息。该操作允许用户系统性地评估死探测器像素16的影响，包括NEMA研究、患者研究中的系统计数性能、分辨率、图像质量和量化等。
在另一范例中，处理器20被编程为使用在110处产生的重建图像来执行建议操作112。基于当前QC图像与参考图像的比较，处理器20确定是对像素化探测器14的维护还是校准被指示为是需要的。在这种情况下，在这种情况下，显示设备20显示根据对是指示对像素化探测器14的维护还是校准的确定来指示执行对像素化探测器14的维护的建议或执行对像素化探测器14的校准的建议的消息。
图3示出了显示对探测器像素16的识别28的指示的显示设备20的范例。识别28示出了探测器像素16的不同分类之间的灵敏度与时间关系。第一曲线1示出了探测器像素16的灵敏度在最大灵敏度阈值以上。在一个实施例中，该像素被标记为“死的”，因为在高灵敏度可能归因于与辐射探测事件不相关的高暗计数率的根本原因下，不管稳定性如何，所有热像素都被标记为“死的”。第二曲线2示出了其性能在原始校准问题之后稳定至可接受灵敏度阈值范围内的探测器像素16。第三曲线3示出了在灵敏度阈值范围内但是其灵敏度在不同时间显著变化的探测器像素16。如果该像素因此被确定为是不稳定的，那么它被标记为“死的”，即使其灵敏度在可接受范围内。第四曲线4示出了在最小灵敏度阈值极限以下(即“冷”)但是其性能稳定的探测器像素16。在稳定的冷像素仍然提供辐射探测事件的有用计数的根本原因下，该稳定的探测器像素16也可以被使用在重建操作中。第五曲线5示出了具有显著降低的灵敏度的探测器像素16。该探测器像素16被识别为处于风险中，并且可能继续在另一校准过程之后失效。因为该像素由于其持续降低的灵敏度而表现出初始不稳定性，并且最终降至下限阈值以下(转变为“冷”像素)，所以它最终被重新标记为“死的”。在这种情况下，在日程计划的维护之前触发警报以进行预防性检查/维护，由此改善系统可靠性。服务可以包括更换模块，或在死像素正以群集出现的情况下，探测器瓦片或模块可以被重新布置以打破这些群集并且改善总体性能。
已经参考优选实施例描述了本公开。他人在阅读和理解以上具体实施方式的情况下可能想到修改和更改。本文旨在将本发明解释为包括所有这种修改和更改，只要它们落入随附权利要求书及其等价方案的范围之内。

Claims (18)

1.一种非瞬态计算机可读介质，存储指令，所述指令能由包括至少一个电子处理器(20)的工作站(18)读取并运行以执行质量控制(QC)方法(100)，所述方法包括：
接收由像素化探测器(14)采集的当前QC数据集和由所述像素化探测器采集的一个或多个先前QC数据集；
根据所述当前QC数据集和所述一个或多个先前QC数据集来确定所述像素化探测器的探测器像素(16)随着时间的稳定性水平，其中，所述稳定性水平被确定为在所述当前QC数据集和所述一个或多个先前QC数据集中的所述像素化探测器的所述探测器像素的灵敏度水平随着时间的变化；
当针对所述探测器像素确定的所述稳定性水平在稳定性阈值范围之外时，将所述像素化探测器的探测器像素标记为死的；以及
在与所述工作站可操作地连接的显示设备(24)上显示对被标记为死的所述探测器像素的识别(28)
所述方法还包括：
在由被标记为死的所述探测器像素(16)采集的数据被省去的情况下重建由所述像素化探测器(14)采集的成像数据集以生成当前QC图像；以及
在所述显示设备(24)上显示所述当前QC图像与利用所述一个或多个先前QC数据集生成的参考图像的比较。
2.根据权利要求1所述的非瞬态计算机可读介质，其中，所述方法(100)还包括：
根据所述当前QC数据集来确定所述像素化探测器(14)的所述探测器像素(16)的灵敏度水平；以及
当所述灵敏度水平在最大灵敏度阈值以上时，将所述像素化探测器的探测器像素标记为死的。
3.根据权利要求2所述的非瞬态计算机可读介质，其中，所述方法(100)还包括：
当所述灵敏度水平在最小灵敏度阈值以下时，将所述像素化探测器(14)的探测器像素(16)标记为死的。
4.根据权利要求2所述的非瞬态计算机可读介质，其中，所述方法(100)还包括：
当所述灵敏度水平在最小灵敏度阈值以下并且具有随着时间的稳定的灵敏度水平时，将所述像素化探测器(14)的探测器像素(16)标记为冷的；
其中，所述显示还包括对被标记为冷的所述探测器像素的识别(28)。
5.根据权利要求1-4中的任一项所述的非瞬态计算机可读介质，其中，所述方法(100)还包括：
基于被标记为死的所述探测器像素(16)来确定是指示对所述像素化探测器(14)的维护还是校准；
其中，在所述显示设备(24)上的所述显示还包括根据对是指示对所述像素化探测器的维护还是校准的所述确定来指示执行对所述像素化探测器的维护的建议或执行对所述像素化探测器的校准的建议的消息，所述消息还指示将被标记为死的探测器像素的群集重新布置在所述像素化探测器中的建议。
6.根据权利要求1-4中的任一项所述的非瞬态计算机可读介质，其中，当针对所述探测器像素确定的所述稳定性水平在稳定性阈值范围之外时将所述像素化探测器(14)的探测器像素(16)标记为死的操作还包括：
当对之前被标记为死的一个或多个探测器像素的探测到的稳定性水平在所述稳定性阈值范围内时，将所述之前被标记为死的一个或多个探测器像素重新标记为活的。
7.根据权利要求1-4中的任一项所述的非瞬态计算机可读介质，其中，所述方法(100)还包括：
在被标记为死的所述探测器像素(16)被处置为不提供数据的情况下模拟由所述像素化探测器(14)针对成像对象采集的成像数据；以及
重建所模拟的成像数据以生成模拟的图像。
8.根据权利要求1所述的非瞬态计算机可读介质，其中，所述方法(100)还包括：
比较所述当前QC图像与所述参考图像；
基于所述当前QC图像与所述参考图像的所述比较来确定是指示对所述像素化探测器(14)的维护还是校准；
其中，在所述显示设备(24)上的所述显示还包括根据对是指示对所述像素化探测器的维护还是校准的所述确定来指示执行对所述像素化探测器的维护的建议或执行对所述像素化探测器的校准的建议的消息。
9.一种成像系统(10)，包括：
图像采集设备(12)，其配置为采集患者的成像数据，所述图像采集设备包括具有多个探测器像素(16)的像素化探测器(14)；
显示设备(24)；以及
至少一个电子处理器(20)，其被编程为：
接收由所述像素化探测器采集的当前质量控制(QC)数据集和由所述像素化探测器采集的一个或多个先前QC数据集；
根据所述当前QC数据集和所述一个或多个先前QC数据集来确定所述探测器像素随着时间的稳定性水平，其中，所述稳定性水平被确定为在所述当前QC数据集和所述一个或多个先前QC数据集中的所述探测器像素的灵敏度水平随着时间的变化；
当针对所述探测器像素确定的所述稳定性水平在稳定性阈值范围之外时，将所述探测器像素中的一个或多个标记为死的；并且
控制所述显示设备以显示对被标记为死的所述探测器像素的识别(28)，
其中，所述至少一个电子处理器(20)还被编程为：
在由被标记为死的所述探测器像素(16)采集的数据被省去的情况下重建由所述像素化探测器(14)采集的成像数据集以生成当前QC图像；以及
控制所述显示设备(24)显示所述当前QC图像与利用所述一个或多个先前QC数据集生成的参考图像的比较。
10.根据权利要求9所述的成像系统(10)，其中，所述至少一个电子处理器(20)还被编程为：
根据所述当前QC数据集来确定所述像素化探测器(14)的所述探测器像素(16)的灵敏度水平；
当所述灵敏度水平在最大灵敏度阈值以上时，将所述像素化探测器的探测器像素标记为死的；
当所述灵敏度水平在最小灵敏度阈值以下时，将所述像素化探测器的探测器像素标记为冷的；
控制所述显示设备(24)以显示对作为死的或冷的所述像素的所述识别(28)。
11.根据权利要求9和10中的任一项所述的成像系统(10)，其中，所述至少一个电子处理器(20)还被编程为：
基于被标记为死的所述探测器像素(16)来确定是指示对所述像素化探测器(14)的维护还是校准；
其中，在所述显示设备(24)上的所述显示还包括根据对是指示对所述像素化探测器的维护还是校准的所述确定来指示执行对所述像素化探测器的维护的建议或执行对所述像素化探测器的校准的建议的消息。
12.根据权利要求9-10中的任一项所述的成像系统(10)，其中，所述至少一个电子处理器(20)还被编程为：
当之前被标记为死的一个或多个探测器像素的探测到的稳定性水平在所述稳定性阈值范围内时，将所述之前被标记为死的一个或多个探测器像素重新标记为活的。
13.根据权利要求9-10中的任一项所述的成像系统(10)，其中，所述至少一个电子处理器(20)还被编程为：
在被标记为死的所述探测器像素(16)被处置为不提供数据的情况下模拟由所述像素化探测器(14)针对成像对象采集的成像数据；并且
重建所模拟的成像数据以生成模拟的图像。
14.根据权利要求9所述的成像系统(10)，其中，所述至少一个电子处理器(20)还被编程为：
比较所述当前QC图像与所述参考图像；
基于对所述当前QC图像与所述参考图像的所述比较来确定是指示对所述像素化探测器(14)的维护还是校准；
其中，在所述显示设备(24)上的所述显示还包括根据对是指示对所述像素化探测器的维护还是校准的所述确定来指示执行对所述像素化探测器的维护的建议或执行对所述像素化探测器的校准的建议的消息。
15.根据权利要求9-10中的任一项所述的成像系统(10)，其中，所述像素化探测器(14)是以下中的一个：
正电子发射断层摄影(PET)设备的探测器环；透射型计算机断层摄影(CT)设备的探测器阵列；被配置为执行单光子发射计算机断层摄影(SPECT)的伽马相机的探测器；以及数字放射摄影设备的数字探测器阵列。
16.根据权利要求9-10中的任一项所述的成像系统(10)，其中，包括点辐射源或线辐射源的参考源的所述当前QC数据集和所述一个或多个先前QC数据集被采集。
17.一种成像系统(10)，包括：
正电子发射断层摄影(PET)设备(12)，其配置为采集患者的成像数据，所述正电子发射断层摄影(PET)设备包括具有多个PET探测器像素(16)的PET探测器环(14)；
显示设备(24)；以及
至少一个电子处理器(20)，其被编程为：
接收由像素化探测器采集的当前质量控制(QC)数据集和由所述像素化探测器采集的一个或多个先前QC数据集；
根据所述当前QC数据集来确定所述探测器像素的灵敏度水平；
根据所述当前QC数据集和所述一个或多个先前QC数据集来确定所述探测器像素随着时间的稳定性水平，其中，所述稳定性水平被确定为在所述当前QC数据集和所述一个或多个先前QC数据集中的所述探测器像素的灵敏度水平随着时间的变化；
当所述灵敏度水平在最大灵敏度阈值以上时，将所述像素化探测器的探测器像素标记为死的；
当所述灵敏度水平在最小灵敏度阈值以下时，将所述像素化探测器的探测器像素标记为冷的；并且
控制所述显示设备以显示对被标记为死的或冷的所述探测器像素的识别(28)，
其中，所述至少一个电子处理器(20)还被编程为：
在由被标记为死的所述探测器像素(16)采集的数据被省去的情况下重建由所述像素化探测器(14)采集的成像数据集以生成当前QC图像；以及
控制所述显示设备(24)显示所述当前QC图像与利用所述一个或多个先前QC数据集生成的参考图像的比较。
18.根据权利要求17所述的成像系统(10)，其中，所述至少一个电子处理器(20)还被编程为：
当之前被标记为死的一个或多个探测器像素的探测到的稳定性水平在最大稳定性阈值和最小稳定性阈值内时，将所述之前被标记为死的一个或多个探测器像素(16)重新标记为活的。
CN201880061097.3A 2017-09-22 2018-09-12 处理数字正电子发射断层摄影的探测器像素性能变化 Active CN111183371B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201762561706P 2017-09-22 2017-09-22
US62/561,706 2017-09-22
PCT/EP2018/074546 WO2019057572A1 (en) 2017-09-22 2018-09-12 MANAGING DETECTOR PIXEL PERFORMANCE VARIATION IN DIGITAL POSITRON EMISSION TOMOGRAPHY

Publications (2)

Publication Number Publication Date
CN111183371A CN111183371A (zh) 2020-05-19
CN111183371B true CN111183371B (zh) 2024-02-13

Family

ID=63579346

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201880061097.3A Active CN111183371B (zh) 2017-09-22 2018-09-12 处理数字正电子发射断层摄影的探测器像素性能变化

Country Status (5)

Country Link
US (1) US11175418B2 (zh)
EP (1) EP3685189A1 (zh)
JP (1) JP7041252B6 (zh)
CN (1) CN111183371B (zh)
WO (1) WO2019057572A1 (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110687585B (zh) * 2019-09-23 2021-05-11 上海联影医疗科技股份有限公司 获取晶体效率的方法、装置、计算机设备和存储介质
CN111624646B (zh) * 2020-06-11 2021-07-30 北京卫星环境工程研究所 多模块辐射探测器的重复事件甄别和处理方法
CN111728631A (zh) * 2020-07-30 2020-10-02 上海联影医疗科技有限公司 Pet系统探测数据修正方法、装置和计算机设备
KR102608895B1 (ko) * 2020-09-10 2023-11-30 신춘관 상하 수직류 순환에어방식의 부상여재의 컴팩트 소형화 여과장치 및 이를 이용한 최적화여재재생시스템
US20240151838A1 (en) * 2021-03-31 2024-05-09 Exo Imaging, Inc. Processing circuitry, system and method to test pixels in an ultrasonic imaging device

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003001243A2 (en) * 2001-06-23 2003-01-03 Koninklijke Philips Electronics Nv Fault-tolerant detector array for gamma ray imaging
EP1487193A1 (en) * 2003-06-13 2004-12-15 Canon Kabushiki Kaisha Method and apparatus for correcting defect pixels in radiation imaging, computer program and computer-readable recording medium
CN102449504A (zh) * 2009-06-01 2012-05-09 皇家飞利浦电子股份有限公司 具有改进的量化能力的pet探测器系统
CN102934526A (zh) * 2010-04-13 2013-02-13 卡尔斯特里姆保健公司 使用数字射线摄影探测器的曝光控制
CN104124256A (zh) * 2014-07-01 2014-10-29 上海奕瑞光电子科技有限公司 一种像素aec平板探测器
CN105700004A (zh) * 2014-12-10 2016-06-22 西门子股份公司 Ct图像中伪影的校正中的噪声抑制
CN105793732A (zh) * 2013-12-04 2016-07-20 皇家飞利浦有限公司 成像探测器自我诊断电路
CN107110983A (zh) * 2014-12-15 2017-08-29 皇家飞利浦有限公司 基于像素的死时间校正

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6457861B1 (en) * 2000-11-15 2002-10-01 Ge Medical Systems Global Technology Company, Llc Method and apparatus for correcting electronic offset and gain variations in a solid state X-ray detector
DE10122876C2 (de) * 2000-11-24 2003-08-21 Siemens Ag Verfahren zum Betrieb eines Bildsystems einer bildgebenden medizinischen Untersuchungseinrichtung und medizinische Untersuchungseinrichtung
JP4995193B2 (ja) 2006-04-28 2012-08-08 株式会社日立メディコ Ｘ線画像診断装置
US8158951B2 (en) * 2009-02-03 2012-04-17 General Electric Company Apparatus and methods for evaluating operation of pixelated detectors
US9201149B2 (en) * 2011-05-16 2015-12-01 Cmt Medical Technologies Ltd. X-ray radiation detector with automatic exposure control
JP5919717B2 (ja) 2011-10-07 2016-05-18 コニカミノルタ株式会社 動態医用画像生成システム
JP6116899B2 (ja) * 2012-01-16 2017-04-19 東芝メディカルシステムズ株式会社 医用画像診断装置及び制御プログラム
BR112014032249A2 (pt) 2012-06-27 2017-06-27 Koninklijke Philips Nv sistema para correção de energia de dados de um evento de tomografia por emissão de pósitron (pet), e método para correção de energia de dados de um evento de tomografia por emissão de pósitron (pet)
JP6338689B2 (ja) * 2014-03-28 2018-06-06 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 陽電子放出断層撮影においけるデッドピクセル識別
JP6383186B2 (ja) * 2014-06-12 2018-08-29 キヤノン株式会社 画像処理装置および画像処理方法、画像処理システム

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003001243A2 (en) * 2001-06-23 2003-01-03 Koninklijke Philips Electronics Nv Fault-tolerant detector array for gamma ray imaging
EP1487193A1 (en) * 2003-06-13 2004-12-15 Canon Kabushiki Kaisha Method and apparatus for correcting defect pixels in radiation imaging, computer program and computer-readable recording medium
CN102449504A (zh) * 2009-06-01 2012-05-09 皇家飞利浦电子股份有限公司 具有改进的量化能力的pet探测器系统
CN102934526A (zh) * 2010-04-13 2013-02-13 卡尔斯特里姆保健公司 使用数字射线摄影探测器的曝光控制
CN105793732A (zh) * 2013-12-04 2016-07-20 皇家飞利浦有限公司 成像探测器自我诊断电路
CN104124256A (zh) * 2014-07-01 2014-10-29 上海奕瑞光电子科技有限公司 一种像素aec平板探测器
CN105700004A (zh) * 2014-12-10 2016-06-22 西门子股份公司 Ct图像中伪影的校正中的噪声抑制
CN107110983A (zh) * 2014-12-15 2017-08-29 皇家飞利浦有限公司 基于像素的死时间校正

Also Published As

Publication number Publication date
CN111183371A (zh) 2020-05-19
JP7041252B6 (ja) 2022-05-31
US11175418B2 (en) 2021-11-16
WO2019057572A1 (en) 2019-03-28
JP2020534539A (ja) 2020-11-26
EP3685189A1 (en) 2020-07-29
US20200301032A1 (en) 2020-09-24
JP7041252B2 (ja) 2022-03-23

Similar Documents

Publication Publication Date Title
CN111183371B (zh) 处理数字正电子发射断层摄影的探测器像素性能变化
CN109564692B (zh) 使用局部修改的飞行时间（tof）内核进行tof pet图像重建
US6934590B2 (en) Quality control system for medical diagnostic apparatus
EP3555669B1 (en) Dead pixel correction for digital pet reconstruction
CN101849834A (zh) 放射线成像设备及其暗电流校正方法
EP3033637B1 (en) A hybrid method based on simulation and experimental data to normalize pet data
US10210635B2 (en) Reconstruction quality assessment with local non-uniformity in nuclear imaging
US20180120459A1 (en) Solving outside-field of view scatter correction problem in positron emission tomography via digital experimentation
US9480450B2 (en) Scatter component estimating method
US10786222B2 (en) Methods, systems, and apparatus for automatically assessing quality of imaging systems
US7915578B2 (en) Method and apparatus for correcting scattering in SPECT imaging
Hamer et al. Amorphous silicon, flat-panel, x-ray detector versus storage phosphor-based computed radiography: contrast-detail phantom study at different tube voltages and detector entrance doses
Pongnapang Practical guidelines for radiographers to improve computed radiography image quality
Brambilla et al. Impact of target‐to‐background ratio, target size, emission scan duration, and activity on physical figures of merit for a 3D LSO‐based whole body PET/CT scanner
Willis et al. Automated quality control assessment of clinical chest images
US10733770B2 (en) System and method for performing fault-tolerant reconstruction of an image
JPH1048341A (ja) 画像生成方法と装置
Kweon et al. Evaluation of the radiation dose to a phantom for various X-ray exposure factors performed using the dose area product in digital radiography
US11759163B2 (en) Radiation imaging system, control apparatus, control method of radiation imaging system, and storage medium
CN113759412B (zh) 获取束流形状和能量探测单元响应特征的方法、装置
US20220022838A1 (en) Evaluation apparatus, evaluation method, and evaluation system
Theodorakis et al. Prospective PET image quality gain calculation method by optimizing detector parameters
CN115153606A (zh) 探测器的死时间校正方法、系统、装置及存储介质
Karam Radiation-based quantitative bioimaging at the national institute of standards and technology
Kelley Quality control in computed tomography

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
GR01 Patent grant
GR01 Patent grant