CN106251313B - 医学成像方法及系统 - Google Patents

医学成像方法及系统 Download PDF

Info

Publication number
CN106251313B
CN106251313B CN201610668696.4A CN201610668696A CN106251313B CN 106251313 B CN106251313 B CN 106251313B CN 201610668696 A CN201610668696 A CN 201610668696A CN 106251313 B CN106251313 B CN 106251313B
Authority
CN
China
Prior art keywords
correction
correction coefficient
scanning
normalized
acquiring
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
CN201610668696.4A
Other languages
English (en)
Other versions
CN106251313A (zh
Inventor
董晓霞
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.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201610668696.4A priority Critical patent/CN106251313B/zh
Publication of CN106251313A publication Critical patent/CN106251313A/zh
Priority to PCT/CN2017/096642 priority patent/WO2018033003A1/en
Priority to US15/569,403 priority patent/US10429525B2/en
Priority to US16/586,701 priority patent/US11105939B2/en
Application granted granted Critical
Publication of CN106251313B publication Critical patent/CN106251313B/zh
Priority to US17/446,461 priority patent/US11874412B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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)
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • 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/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/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Medical Informatics (AREA)
  • Data Mining & Analysis (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Optics & Photonics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Nuclear Medicine (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明提供了一种医学成像方法及系统，所述方法包括：获取扫描数据；获取校正系数，所述校正系数包括归一化校正系数和总计数定量一致性校正系数；使用所述校正系数对所述扫描数据进行校正，获得校正后的扫描数据；使用所述校正后的扫描数据进行图像重建，获得重建图像。本发明方法可以实现用户在对扫描数据进行归一化校正处理后，重建图像的SUV值保持一致。

Description

医学成像方法及系统
技术领域
本申请涉及医学成像技术领域，特别是涉及一种医学成像方法及系统。
背景技术
正电子发射断层成像(positron emission tomography,PET)系统的探测器是由许多模块组成的，众多晶体块的光子产额及光电倍增管的增益难免存在差异，这种性能的非均一性会造成各探测器模块探测效率的不一致。另外，探测器的灵敏度与探头的几何形状、几何尺寸有关。对造成探测器灵敏度差异的因素进行校正的过程，称为探测器的归一化。
校正的方法是应用归一化系数对每条响应线进行归一化运算。归一化系数是通过对已知强度均匀的符合源长时间采集，获得探测视野内所有符合计数的平均值和每一对探测器的归一化系数。归一化系数可分解为两种影像因素(探测器效率和几何影响因素)的乘积。探测器效率会随时间而发生变化，须定期校正；几何影响因素不会随时间变化，经一次校正采集和计算后即可确定。在PET出厂时及以后定期用旋转线源或均匀的柱源测试系统探测效率的不均匀性，并且将测试结果转换成校正系数图供数据处理使用。
目前，归一化校正方法主要有自归一化(self-normalization)校正方法和组分归一化(component-based normalization)校正方法，其中，组分归一化校正方法考虑不同因子对探测器灵敏度差异的影响。采用组分归一化校正方法对探测器进行归一化校正时，容易导致用于重建的响应线总计数不同，表现在图像上即图像的SUV(Standardized UptakeValue，标准摄取值)定量不一致。
发明内容
本申请要解决的问题是采用组分归一化校正方法对探测器进行归一化校正时，容易导致用于重建的响应线总计数不同。
为实现上述目的，本申请提供了一种医学成像方法，包括：获取扫描数据；获取校正系数，所述校正系数包括归一化校正系数和总计数定量一致性校正系数；使用所述校正系数对所述扫描数据进行校正，获得校正后的扫描数据；使用所述校正后的扫描数据进行图像重建，获得重建图像。
在本申请的一些实施例中，上述归一化校正系数包括通过组分归一化校正方法计算得到的归一化校正系数。
在本申请的一些实施例中，上述获取校正系数包括：获取扫描时的计数率；获取所述扫描时的计数率对应的校正系数。
在本申请的一些实施例中，上述获取扫描时的计数率对应的校正系数包括：获取计数率与校正系数的对应关系；根据所述扫描时的计数率和所述计数率与校正系数的对应关系，确定所述扫描时的计数率对应的校正系数。
在本申请的一些实施例中，上述获取计数率与校正系数的对应关系包括：获取第一查找表，所述第一查找表中记录计数率与归一化校正系数的对应关系。
在本申请的一些实施例中，上述获取计数率与校正系数的对应关系包括：获取第二查找表，所述第二查找表中记录计数率与总计数定量一致性校正系数的对应关系。
为实现上述目的，本发明还提供了一种医学成像系统，包括扫描设备和处理器，所述扫描设备与处理器通讯，用于响应处理器发出的扫描指令，进行扫描数据采集，并将所采集的扫描数据发送给处理器；所述处理器包括接收模块、存储模块及数据处理模块，所述接收模块用于获取扫描数据和校正系数，所述存储模块用于存储扫描数据、多种校正系数，所述多种校正系数包括归一化校正系数和总计数定量一致性校正系数，所述数据处理模块用于使用所述校正系数对所述扫描数据进行校正，获得校正后的扫描数据，使用校正后的扫描数据进行图像重建，获得重建图像。
在本申请的一些实施例中，上述存储模块进一步用于存储第一查找表，所述第一查找表中记录计数率与归一化校正系数的对应关系。
在本申请的一些实施例中，上述存储模块进一步用于存储第二查找表，所述第二查找表中记录计数率与总计数定量一致性校正系数的对应关系。
在本申请的一些实施例中，上述归一化校正系数包括通过组分归一化校正方法计算得到的归一化校正系数。
本申请的一部分附加特性可以在下面的描述中进行说明。通过对以下描述和相应附图的检查或者对实施例的生产或操作的了解，本申请的一部分附加特性对于本领域技术人员是明显的。本披露的特性可以通过对以下描述的具体实施例的各种方面的方法、手段和组合的实践或使用得以实现和达到。
附图说明
图1是根据本申请的一些实施例所示的医学成像系统的示意图；
图2是根据本申请的一些实施例所示的处理器的示意图；
图3是根据本申请的一些实施例所示的医学成像方法的流程图；
图4是根据本申请的一些实施例所示的获取扫描时的计数率对应的归一化校正系数的流程图；
图5是根据本申请的一些实施例所示的对扫描数据进行总计数定量一致性校正的流程图；
图6是根据本申请的一些实施例所示的获取扫描时的计数率对应的总计数定量一致性校正系数的流程图；
图7是根据本申请的一些实施例所示的获取计数率与总计数定量一致性校正系数对应关系的流程图；
图8是根据本申请的一些实施例所示的第一归一化校正图像和计数率的对应关系曲线与第二归一化校正图像与计数率的对应关系曲线的对比示意图；
图9是根据本申请的一些实施例所示的拟合得到的计数率与总计数定量一致性校正系数的关系曲线的示意图。
具体实施方式
为了更清楚地说明本申请的实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单的介绍。显而易见地，下面描述中的附图仅仅是本申请的一些示例或实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图将本申请应用于其他类似情景。除非从语言环境中显而易见或另做说明，图中相同标号代表相同结构或操作。
如本申请和权利要求书中所示，除非上下文明确提示例外情形，“一”、“一个”、“一种”和/或“该”等词并非特指单数，也可包括复数。一般说来，术语“包括”与“包含”仅提示包括已明确标识的步骤和元素，而这些步骤和元素不构成一个排它性的罗列，方法或者设备也可能包含其他的步骤或元素。
虽然本申请对根据本申请的实施例的系统中的某些模块做出了各种引用，然而，任何数量的不同模块可以被使用并运行在客户端和/或服务器上。所述模块仅是说明性的，并且所述系统和方法的不同方面可以使用不同模块。
本申请中使用了流程图用来说明根据本申请的实施例的系统所执行的操作。应当理解的是，前面或下面操作不一定按照顺序来精确地执行。相反，可以按照倒序或同时处理各种步骤。同时，也可以将其他操作添加到这些过程中，或从这些过程移除某一步或数步操作。
图1是根据本申请的一些实施例所示的医学成像系统的示意图。医学成像系统100可以包含一个处理器120、一个网络130和一个扫描设备110。处理器120对收集到的信息(例如数据等)进行处理。处理器120可以是一个实体的电子设备，也可以是一个服务器。所述电子设备可以包括计算机、平板、手机、智能终端设备等。处理器120可以是集中式的，例如数据中心；也可以是分布式的，例如一个分布式系统。处理器120可以是本地的，也可以是远程的。在一些实施例中，所述信息可以是通过扫描或其他方式获得的一个或多个对象的扫描数据、成像参数等。
在一些实施例中，处理器120可以包括中央处理器(Central Processing Unit，CPU)、专门应用集成电路(Application Specific Integrated Circuit，ASIC)、专用指令处理器(Appl ication Specific Instruction Set Processor，ASIP)、物理处理器(Physics Processing Unit，PPU)、数字信号处理器(Digital Processing Processor，DSP)、现场可编程逻辑门阵列(Field-Programmable Gate Array，FPGA)、可编程逻辑器件(Programmable Logic Device，PLD)、处理器、微处理器、控制器、微控制器等中的一种或几种的组合。
网络130可以是单个网络，也可以是多个不同网络的组合。例如，网络130可能是一个局域网(Local Area Network，LAN)、广域网(Wide Area Network，WAN)、公用网络、私人网络、专有网络、公共交换电话网(Public Switched Telephone Network，PSTN)、互联网、无线网络、虚拟网络、或者上述网络的任何组合。网络130也可以包括多个网络接入点。有线网络可以包括利用金属电缆、混合电缆、一个或多个接口等一种或多种组合的方式。无线网络可以包括利用蓝牙、区域局域网(LAN)、广域局域网(WAN)、无线个域网(WPAN)、近源场通信(Near Field Communication，NFC)等一种或多种组合的方式。网络130可以适用于本申请所描述的范围内，但并不局限于所述描述。
扫描设备110可以包括对一个或多个对象进行扫描的一个或多个设备，进一步地，所述用于扫描的设备可以被用在但不仅限于医学领域的应用，例如医学检测等。在一些实施例中，医学检测可以包括正电子发射断层成像(PET)、单光子发射计算机断层显像(SPECT)、正电子发射断层成像-计算机断层成像(PET-CT)、正电子发射断层成像-磁共振成像(PET-MR)或者上述一种或多种医学检测的组合。在一些实施例中，所述对象可以是器官、机体、物体、机能障碍、肿瘤等一种或多种的组合。在一些实施例中，所述对象可以是头部、胸腔、器官、骨骼、血管等一种或多种的组合。在一些实施例中，扫描设备110可以由一个或多个检测模块拼接而成。进一步地，所述一个或多个检测模块的探测器可以连续地放置在所述对象的周围。
在一些实施例中，扫描设备110和处理器120可以是一体的。在一些实施例中，扫描设备110可以通过网络130发送信息到处理器120。在一些实施例中，扫描设备110也可以直接发送信息到处理器120。在一些实施例中，处理器120也可以包含处理本身存储的信息。
图2是根据本申请的一些实施例所示的处理器的示意图。处理器120可以包含一个或多个显示模块210、一个或多个接收模块230、一个或多个存储模块220、一个或多个数据处理模块240。
显示模块210具有显示界面，可以显示用户输入的信息、数据处理过程信息及数据处理结果。在一些实施例中，显示模块210可以显示扫描对象的基本信息、设定的扫描参数、重建参数、校正模式、重建图像等信息。在一些实施例中，显示模块210可以显示归一化校正选项、死时间校正选项、衰减校正选项、散射校正选项等中的一种或多种的组合。在一些实施例中，显示模块210可以显示用户选定的归一化校正选项对应的归一化校正系数、选定的死时间校正选项对应的死时间校正系数。在一些实施例中，显示模块210可以包括一个显示设备，如显示屏等。在一些实施例中，用户可以根据显示模块240所显示的图像进行设置和/或操作，如设置为二维图像的显示、设置为三维图像的显示、显示扫描数据对应的图像、显示控制界面、显示输入界面、显示不同区域的图像、显示图像重建的过程、显示图像重建的结果，接受到用户的输入后对显示图像进行放大处理、缩小处理、设置多个图像同时显示等一种或几种设置和/或操作的组合。在一些实施例中，显示模块210可以包括一个或多个输入设备，如键盘、触屏、触板、鼠标、远程控制、语音识别设备等一个或多个。在一些实施例中，用户可以通过所述一个或多个输入设备输入一些原始参数和/或设置对应图像显示和/或处理的初始化条件。在一些实施例中，用户通过所述一个或多个输入设备输入选定的对扫描数据进行的校正项，例如包括归一化校正、死时间校正、衰减校正等中的一种或多种的组合。
存储模块220可以存储数据。所述存储的数据可以来自扫描设备110，网络130，和/或处理器120中的其他模块/单元(显示模块210或其他相关模块(未示出))。在一些实施例中，存储模块220可以存储扫描数据、归一化校正系数、死时间校正系数、衰减校正系数等各种校正系数。在一些实施例中，存储模块220可以存储查找表，例如，记录计数率与死时间校正系数的对应关系的查找表、记录计数率与总计数定量一致性校正系数的对应关系的查找表、记录计数率与归一化校正系数的对应关系的查找表等。存储模块220可以是利用电能方式存储信息的设备，例如各种存储器，如随机存取存储器(Random Access Memory(RAM))、只读存储器(Read Only Memory(ROM))等。其中随机存储器可以包括十进计数管、选数管、延迟线存储器、威廉姆斯管、动态随机存储器(DRAM)、静态随机存储器(SRAM)、晶闸管随机存储器(T-RAM)、零电容随机存储器(Z-RAM)等中的一种或几种的组合。只读存储器可以包括磁泡存储器、磁钮线存储器、薄膜存储器、磁镀线存储器、磁芯内存、磁鼓存储器、光盘驱动器、硬盘、磁带、早期非易失存储器(NVRAM)、相变化内存、磁阻式随机存储式内存、铁电随机存储内存、非易失SRAM、闪存、电子抹除式可复写只读存储器、可擦除可编程只读存储器、可编程只读存储器、屏蔽式堆读内存、浮动连接门随机存取存储器、纳米随机存储器、赛道内存、可变电阻式内存、可编程金属化单元等中的一种或几种的组合。存储模块220可以是利用磁能方式存储信息的设备，例如硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘、闪存等。存储模块220可以是利用光学方式存储信息的设备，例如CD或DVD等。存储模块220可以是利用磁光方式存储信息的设备，例如磁光盘等。存储模块220的存取方式可以是随机存储、串行访问存储、只读存储等中的一种或几种的组合。存储模块220可以是非永久记忆存储器，或永久记忆存储器。
存储模块220可以与一个或多个显示模块210、接收模块230、数据处理模块240或其他相关模块(未示出)关联。在一些实施例中，存储模块220可以通过网络130选择性地关联一个或多个虚拟存储资源，例如云盘存储(cloud storage)、虚拟私人网络(a virtualprivate network)和/或其他虚拟存储资源。存储的数据可以是各种形式的数据，例如数值、信号、图像、既定目标的相关信息、命令、算法、程序等一种或多种组合。
接收模块230可以以一种或多种方式收集所需要的信息。所述收集信息的方式可以包括扫描一个对象(例如通过成像设备110获取一个对象的信息)，通过收集预先存储的信息(例如通过收集存储模块220中的信息或通过网络130获得的远程信息)等。信息的种类可以包括体素数据、计数、矩阵、图像、向量、向量库等。
数据处理模块240可以处理接收模块230所收集的信息。在一些实施例中，数据处理模块可以对扫描数据进行处理，例如，在图像重建前进行一系列的数据处理前的准备工作，包括计算衰减系数、弓形几何校正位置、散射校正滤波函数、衰减校正系数和重建平滑滤波函数等；对扫描数据进行校正计算，包括探测器效率归一化、散射符合、衰减、死时间、弓形几何失真和辐射衰变等；使用校正后的数据进行图像重建运算。在一些实施例中，数据处理模块可以对重建图像进行处理，例如计算参数，包括刻度系数(scale factor)、像素尺寸(pixel size)和标定系数(calibration factor)等。
对于本领域的专业人员来说，在了解本申请的医学成像方法及系统的原理后，可能在不背离这一原理的情况下，对各个模块进行任意组合，或者构成子系统与其他模块连接，对实施上述方法和系统的应用领域形式和细节上的各种修正和改变，但是这些修正和改变仍在以上描述的范围之内。例如，上述模块可以是体现在一个系统中的不同模块，也可以是一个模块实现上述的两个或两个以上模块的功能。比如，在本申请的一些实施例中，存储模块220可以被包含在任何一个或多个所述模块中。
图3是根据本申请的一些实施例所示的医学成像方法的流程图。在一些实施例中，过程300可以由处理器120实现。如步骤301所示，获取扫描数据。处理器120可以获取一个对象的扫描数据。在一些实施例中，处理器120可以通过成像设备110获取所述扫描数据。进一步地，所述成像设备110可以包括PET成像设备。在一些实施例中，所述扫描数据可以从存储模块220中获得。在一些实施例中，所述扫描数据也可以通过网络130从远程存储模块(如云盘)中获得。在获取扫描数据后，处理器120可以对扫描数据执行如步骤303所示的校正操作及步骤304所示的重建图像操作。
如步骤302所示，获取校正系数。在一些实施例中，该校正系数可以包括归一化校正系数、死时间校正系数、衰减校正系数、总计数定量一致性校正系数等中的一种或多种的组合。数据处理模块240可以由存储模块220、接收模块230、或系统100的其他装置等中获取归一化校正系数。在一些实施例中，存储模块220可以存储一种或多种归一化校正系数。在一些实施例中，可以从存储模块220获取归一化校正系数和总计数定量一致性校正系数。在其他实施例中，可以从显示模块或系统100的其他存储器中获取归一化校正系数，例如可以从扫描设备110中的存储器中获取。在一些实施例中，存储模块220可以存储一种或多种归一化校正系数。在一些实施例中，归一化校正系数可以是通过自归一化(self-normalization)校正方法计算得到的归一化校正系数。在一些实施例中，归一化校正系数可以是通过组分归一化(component-based normalization)校正方法计算得到的归一化校正系数。在一些实施例中，归一化校正系数可以预先通过实验获取，例如，可以通过对已知强度均匀的符合源长时间采集，获取探测视野内所有符合计数的平均值和每一对探测器的符合的计数值，二者之比即为各对探测器的归一化校正系数。
在一些实施例中，可以使用现有的组分归一化校正方法计算得到归一化校正系数，例如Phys.Med.Biol.44(1999)571–594公开了一种使用组分归一化校正方法计算得到的归一化校正系数可以表述为公式1的形式：
NCuivj＝εuiεvjcuimodDcvjmodDduvrkfuvguvr (1)
其中，D表示沿探测器一边的晶体数量，r表示响应线(LOR)的径向位置，k表示探测器模块中响应线的相对位置，εui、εvj表示位于探测器环u、v中的晶体i、j的晶体效率，cuimodD、cvjmodD表示探测器模块内晶体相对位置不同引起的效率变化因子，duvrk表示晶体干扰因子，fuv表示轴向几何因子，guvr表示径向几何因子。
在一些实施例中，组分归一化校正系数包含的多种归一化因子中的部分归一化因子可以与计数率有关，例如，计数率不同，其对应的归一化因子不同。图4是根据本申请的一些实施例所示的获取扫描时的计数率对应的归一化校正系数的流程图。过程400可以由处理器120实现。如步骤401所示，获取扫描时的计数率。在一些实施例中，可以在扫描时，统计计数率。如步骤402所示，获取计数率与归一化校正系数的对应关系。在一些实施例中，数据处理模块240可以从存储模块220获取计数率与归一化校正系数的对应关系，例如，存储模块220可以存储第一查找表，该第一查找表中记录计数率与归一化校正系数的对应关系。在一些实施例中，计数率与归一化校正系数的对应关系可以预先通过实验确定，例如，通过对已知强度均匀的符合源长时间采集，获得探测视野内所有符合计数的平均值和每一对探测器的符合的计数值，二者之比即为各对探测器的归一化因子。最后如步骤403所示，根据扫描时的计数率和计数率与归一化校正系数的对应关系，可以确定扫描时的计数率对应的归一化校正系数。
如步骤303所示，使用所述校正系数对所述扫描数据进行校正，获得校正后的扫描数据。在一些实施例中，该校正操作可以由数据处理模块240实现。在一些实施例中，数据处理模块240从存储模块220中获取选定的归一化校正系数和/或总计数定量一致性校正系数，对扫描数据进行归一化校正处理和/或总计数定量一致性校正处理。
如步骤304所示，使用所述校正后的扫描数据进行图像重建，获得重建图像。重建图像可以通过显示模块210呈现给用户。
本申请中，用户在使用获取的归一化校正系数对扫描数据进行归一化校正时，可能导致用于重建的响应线总计数和/或计数率不同，表现在重建图像上即图像的SUV定量不一致。为解决该问题，在一些实施例中，可以使用总计数定量一致性校正系数对扫描数据进行校正。仅仅作为示例，获取组分归一化校正系数，使用组分归一化校正系数对扫描数据进行校正，为确保用于重建的响应线总计数和/或计数率保持一致，可以使用总计数定量一致性校正系数对扫描数据作进一步校正，得到SUV定量一致的重建图像。
图5是根据本申请的一些实施例所示的对扫描数据进行总计数定量一致性校正的流程图。过程500可以由处理器120实现。如步骤501所示，获取扫描数据。所述扫描数据可以来自接收模块230、存储模块220、或系统100中的其他存储器等中的一种或多种的组合。在一些实施例中，所述扫描数据可以是数据处理模块240中经归一化校正后的数据。
如步骤502所示，获取总计数定量一致性校正系数。在一些实施例中，总计数定量一致性校正系数可以来自存储模块220、接收模块230、或系统100的其他存储器等中的一种或多种的组合。在一些实施例中，存储模块220可以存储一种或多种总计数定量一致性校正系数。在一些实施例中，获取的总计数定量一致性校正系数可以与选定的归一化校正模式有关，例如，选定的归一化校正模式不同，其对应的获取的总计数定量一致性校正系数可以不同。在一些实施例中，选定的总计数定量一致性校正系数可以与选定的归一化校正模式无关，例如，选定的归一化校正模式不同，其对应的总计数定量一致性校正系数可以相同。在一些在本申请的一些实施例中实施例中，总计数定量一致性校正系数可以与计数率有关，例如，获取的总计数定量一致性校正系数可以是扫描时的计数率对应的总计数定量一致性校正系数。
图6是根据本申请的一些实施例所示的获取扫描时的计数率对应的总计数定量一致性校正系数的流程图。过程600可以由处理器120实现。如步骤601所示，获取扫描时的计数率。在一些实施例中，可以在扫描时，统计计数率。如步骤602所示，获取计数率与总计数定量一致性校正系数的对应关系。在一些实施例中，数据处理模块240可以从存储模块220获取计数率与总计数定量一致性校正系数的对应关系，例如，存储模块220可以存储第二查找表，该第二查找表中记录计数率与总计数定量一致性校正系数的对应关系。在一些实施例中，计数率与总计数定量一致性校正系数的对应关系可以预先通过实验确定。最后如步骤603所示，根据所述扫描时的计数率和所述计数率与总计数定量一致性校正系数的对应关系，确定所述扫描时的计数率对应的总计数定量一致性校正系数。
图7是根据本申请的一些实施例所示的获取计数率与总计数定量一致性校正系数对应关系的流程图。过程700可以由处理器120实现。如步骤701所示，可以对均匀模体进行扫描，获取属于多个时间间隔的扫描数据。所述模体可以是造影剂均匀分布的水模、或其他均匀模体。获取属于多个时间间隔的扫描数据，该属于多个时间间隔的扫描数据对应所述均匀模体从高活度到低活度的一系列衰变数据。
请参考图8，如步骤702所示，使用通过自归一化校正方法计算得到的归一化校正系数校正所述多个时间间隔的扫描数据，使用校正后的所述多个时间间隔的扫描数据进行图像重建，获得所述多个时间间隔中的每一个时间间隔对应的第一归一化校正图像。如步骤703所示，使用通过组分归一化校正方法计算得到的归一化校正系数校正所述多个时间间隔的扫描数据，使用校正后的所述多个时间间隔的扫描数据进行图像重建，获得所述多个时间间隔中的每一个时间间隔对应的第二归一化校正图像。
如步骤704至706所示，计算所述多个时间间隔中的每一个时间间隔对应的第一归一化校正图像和第二归一化校正图像的比值；根据所述多个时间间隔中的每一个时间间隔对应的第一归一化校正图像和第二归一化校正图像的比值确定所述多个时间间隔中的每一个时间间隔对应的总计数定量一致性校正系数；获取所述多个时间间隔中的每一个时间间隔对应的扫描时的计数率；根据所述多个时间间隔中的每一个时间间隔对应的总计数定量一致性校正系数及扫描时的计数率，建立计数率与总计数定量一致性校正系数的对应关系。在一些实施例中，计数率与总计数定量一致性校正系数的对应关系可以是计数率与总计数定量一致性校正系数的关系曲线。进一步地，可以对计数率与总计数定量一致性校正系数的关系曲线进行拟合操作，例如，二次拟合或线性拟合等。图9所示为执行二次拟合操作后得到的总计数定量一致性校正系数的关系曲线。
上文已对基本概念做了描述，显然，对于本领域技术人员来说，上述发明披露仅仅作为示例，而并不构成对本申请的限定。虽然此处并没有明确说明，本领域技术人员可能会对本申请进行各种修改、改进和修正。该类修改、改进和修正在本申请中被建议，所以该类修改、改进、修正仍属于本申请示范实施例的精神和范围。
同时，本申请使用了特定词语来描述本申请的实施例。如“一个实施例”、“一实施例”、和/或“一些实施例”意指与本申请至少一个实施例相关的某一特征、结构或特点。因此，应强调并注意的是，本说明书中在不同位置两次或多次提及的“一实施例”或“一个实施例”或“一替代性实施例”并不一定是指同一实施例。此外，本申请的一个或多个实施例中的某些特征、结构或特点可以进行适当的组合。
此外，本领域技术人员可以理解，本申请的各方面可以通过若干具有可专利性的种类或情况进行说明和描述，包括任何新的和有用的工序、机器、产品或物质的组合，或对他们的任何新的和有用的改进。相应地，本申请的各个方面可以完全由硬件执行、可以完全由软件(包括固件、常驻软件、微码等)执行、也可以由硬件和软件组合执行。以上硬件或软件均可被称为“数据块”、“模块”、“引擎”、“单元”、“组件”或“系统”。此外，本申请的各方面可能表现为位于一个或多个计算机可读介质中的计算机产品，该产品包括计算机可读程序编码。
例如，在本申请的一些实施例中，提供了一种医学成像系统，包括扫描设备和处理器，所述扫描设备与处理器通讯，用于响应处理器发出的扫描指令，进行扫描数据采集，并将所采集的扫描数据发送给处理器；所述处理器包括接收模块、存储模块及数据处理模块，所述接收模块用于获取扫描数据和校正系数，所述存储模块用于存储扫描数据、多种校正系数，所述多种校正系数包括归一化校正系数和总计数定量一致性校正系数，所述数据处理模块用于使用所述校正系数对所述扫描数据进行校正，获得校正后的扫描数据，使用校正后的扫描数据进行图像重建，获得重建图像。在本申请的一些实施例中，上述存储模块进一步用于存储第一查找表，所述第一查找表中记录计数率与归一化校正系数的对应关系。在本申请的一些实施例中，上述存储模块进一步用于存储第二查找表，所述第二查找表中记录计数率与总计数定量一致性校正系数的对应关系。
计算机可读信号介质可能包含一个内含有计算机程序编码的传播数据信号，例如在基带上或作为载波的一部分。该传播信号可能有多种表现形式，包括电磁形式、光形式等等、或合适的组合形式。计算机可读信号介质可以是除计算机可读存储介质之外的任何计算机可读介质，该介质可以通过连接至一个指令执行系统、装置或设备以实现通讯、传播或传输供使用的程序。位于计算机可读信号介质上的程序编码可以通过任何合适的介质进行传播，包括无线电、电缆、光纤电缆、RF、或类似介质、或任何上述介质的组合。
本申请各部分操作所需的计算机程序编码可以用任意一种或多种程序语言编写，包括面向对象编程语言如Java、Scala、Smalltalk、Eiffel、JADE、Emerald、C++、C#、VB.NET、Python等，常规程序化编程语言如C语言、Visual Basic、Fortran 2003、Perl、COBOL 2002、PHP、ABAP，动态编程语言如Python、Ruby和Groovy，或其他编程语言等。该程序编码可以完全在用户计算机上运行、或作为独立的软件包在用户计算机上运行、或部分在用户计算机上运行部分在远程计算机运行、或完全在远程计算机或服务器上运行。在后种情况下，远程计算机可以通过任何网络形式与用户计算机连接，比如局域网(LAN)或广域网(WAN)，或连接至外部计算机(例如通过因特网)，或在云计算环境中，或作为服务使用如软件即服务(SaaS)。
此外，除非权利要求中明确说明，本申请所述处理元素和序列的顺序、数字字母的使用、或其他名称的使用，并非用于限定本申请流程和方法的顺序。尽管上述披露中通过各种示例讨论了一些目前认为有用的发明实施例，但应当理解的是，该类细节仅起到说明的目的，附加的权利要求并不仅限于披露的实施例，相反，权利要求旨在覆盖所有符合本申请实施例实质和范围的修正和等价组合。例如，虽然以上所描述的系统组件可以通过硬件设备实现，但是也可以只通过软件的解决方案得以实现，如在现有的服务器或移动设备上安装所描述的系统。
同理，应当注意的是，为了简化本申请披露的表述，从而帮助对一个或多个发明实施例的理解，前文对本申请实施例的描述中，有时会将多种特征归并至一个实施例、附图或对其的描述中。但是，这种披露方法并不意味着本申请对象所需要的特征比权利要求中提及的特征多。实际上，实施例的特征要少于上述披露的单个实施例的全部特征。
一些实施例中使用了描述成分、属性数量的数字，应当理解的是，此类用于实施例描述的数字，在一些示例中使用了修饰词“大约”、“近似”或“大体上”来修饰。除非另外说明，“大约”、“近似”或“大体上”表明所述数字允许有±20％的变化。相应地，在一些实施例中，说明书和权利要求中使用的数值参数均为近似值，该近似值根据个别实施例所需特点可以发生改变。在一些实施例中，数值参数应考虑规定的有效数位并采用一般位数保留的方法。尽管本申请一些实施例中用于确认其范围广度的数值域和参数为近似值，在具体实施例中，此类数值的设定在可行范围内尽可能精确。
针对本申请引用的每个专利、专利申请、专利申请公开物和其他材料，如文章、书籍、说明书、出版物、文档、物件等，特此将其全部内容并入本申请作为参考。与本申请内容不一致或产生冲突的申请历史文件除外，对本申请权利要求最广范围有限制的文件(当前或之后附加于本申请中的)也除外。需要说明的是，如果本申请附属材料中的描述、定义、和/或术语的使用与本申请所述内容有不一致或冲突的地方，以本申请的描述、定义和/或术语的使用为准。
最后，应当理解的是，本申请中所述实施例仅用以说明本申请实施例的原则。其他的变形也可能属于本申请的范围。因此，作为示例而非限制，本申请实施例的替代配置可视为与本申请的教导一致。相应地，本申请的实施例不仅限于本申请明确介绍和描述的实施例。

Claims (8)

1.一种医学成像方法，其特征在于，包括：
获取扫描数据；
获取校正系数，所述校正系数包括归一化校正系数和总计数定量一致性校正系数，所述归一化校正系数包括通过组分归一化校正方法计算得到的归一化校正系数；
使用所述校正系数对所述扫描数据进行校正，获得校正后的扫描数据；
使用所述校正后的扫描数据进行图像重建，获得重建图像。
2.根据权利要求1所述的方法，其特征在于，所述获取校正系数包括：
获取扫描时的计数率；
获取所述扫描时的计数率对应的校正系数。
3.根据权利要求2所述的方法，其特征在于，所述获取扫描时的计数率对应的校正系数包括：
获取计数率与校正系数的对应关系；
根据所述扫描时的计数率和所述计数率与校正系数的对应关系，确定所述扫描时的计数率对应的校正系数。
4.根据权利要求3所述的方法，其特征在于，所述获取计数率与校正系数的对应关系包括：获取第一查找表，所述第一查找表中记录计数率与归一化校正系数的对应关系。
5.根据权利要求3所述的方法，其特征在于，所述获取计数率与校正系数的对应关系包括：获取第二查找表，所述第二查找表中记录计数率与总计数定量一致性校正系数的对应关系。
6.一种医学成像系统，其特征在于，包括扫描设备和处理器，
所述扫描设备与处理器通讯，用于响应处理器发出的扫描指令，进行扫描数据采集，并将所采集的扫描数据发送给处理器；
所述处理器包括接收模块、存储模块及数据处理模块，所述接收模块用于获取扫描数据和校正系数，所述存储模块用于存储扫描数据、多种校正系数，所述多种校正系数包括归一化校正系数和总计数定量一致性校正系数，所述归一化校正系数包括通过组分归一化校正方法计算得到的归一化校正系数，所述数据处理模块用于使用所述校正系数对所述扫描数据进行校正，获得校正后的扫描数据，使用校正后的扫描数据进行图像重建，获得重建图像。
7.根据权利要求6所述的系统，其特征在于，所述存储模块进一步用于存储第一查找表，所述第一查找表中记录计数率与归一化校正系数的对应关系。
8.根据权利要求6所述的系统，其特征在于，所述存储模块进一步用于存储第二查找表，所述第二查找表中记录计数率与总计数定量一致性校正系数的对应关系。
CN201610668696.4A 2016-08-15 2016-08-15 医学成像方法及系统 Active CN106251313B (zh)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CN201610668696.4A CN106251313B (zh) 2016-08-15 2016-08-15 医学成像方法及系统
PCT/CN2017/096642 WO2018033003A1 (en) 2016-08-15 2017-08-09 System and method for pet data correction
US15/569,403 US10429525B2 (en) 2016-08-15 2017-08-09 System and method for PET data correction
US16/586,701 US11105939B2 (en) 2016-08-15 2019-09-27 System and method for PET data correction
US17/446,461 US11874412B2 (en) 2016-08-15 2021-08-30 System and method for PET data correction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610668696.4A CN106251313B (zh) 2016-08-15 2016-08-15 医学成像方法及系统

Publications (2)

Publication Number Publication Date
CN106251313A CN106251313A (zh) 2016-12-21
CN106251313B true CN106251313B (zh) 2020-06-26

Family

ID=57592771

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610668696.4A Active CN106251313B (zh) 2016-08-15 2016-08-15 医学成像方法及系统

Country Status (3)

Country Link
US (3) US10429525B2 (zh)
CN (1) CN106251313B (zh)
WO (1) WO2018033003A1 (zh)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106251313B (zh) 2016-08-15 2020-06-26 上海联影医疗科技有限公司 医学成像方法及系统
JP2019007825A (ja) * 2017-06-23 2019-01-17 株式会社島津製作所 核医学診断装置および核医学診断装置の日常保守および点検方法
CN117017334A (zh) * 2017-10-31 2023-11-10 上海联影医疗科技股份有限公司 Pet校正系统和方法
CN115153608A (zh) 2017-11-08 2022-10-11 上海联影医疗科技股份有限公司 校正投影图像的系统和方法
WO2019090541A1 (en) 2017-11-08 2019-05-16 Shenzhen United Imaging Healthcare Co., Ltd. Systems and methods for correcting projection images in computed tomography image reconstruction
WO2019090533A1 (en) 2017-11-08 2019-05-16 Shenzhen United Imaging Healthcare Co., Ltd. System and method for diagnositic and treatment
CN108287361B (zh) * 2018-01-03 2019-08-27 东软医疗系统股份有限公司 一种单事件死时间的检测方法及装置
CN108186043B (zh) * 2018-02-09 2021-03-30 东软医疗系统股份有限公司 校正pet设备灵敏度的方法和装置
US11867627B2 (en) * 2018-10-12 2024-01-09 Washington University Compact guided diffuse optical tomography system for imaging a lesion region
CN110084800B (zh) * 2019-04-28 2021-03-16 上海海事大学 一种用于四肢软组织肉瘤病人的肺转移预测方法
CN109965897B (zh) * 2019-05-10 2022-03-22 上海联影医疗科技股份有限公司 Pet扫描仪校正方法、装置、计算机设备和可读存储介质
CN110687585B (zh) * 2019-09-23 2021-05-11 上海联影医疗科技股份有限公司 获取晶体效率的方法、装置、计算机设备和存储介质
CN110706175B (zh) * 2019-09-27 2022-07-29 上海联影医疗科技股份有限公司 Pet校正系数的生成方法、系统、可读存储介质和设备
US11763498B2 (en) * 2019-11-12 2023-09-19 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for image reconstruction
CN111882499B (zh) * 2020-07-15 2024-04-16 上海联影医疗科技股份有限公司 Pet图像的降噪方法、装置以及计算机设备
CN111815734B (zh) * 2020-08-21 2024-04-16 上海联影医疗科技股份有限公司 一种pet扫描数据修复方法、装置和计算机设备
CN112817035B (zh) * 2021-01-28 2022-11-25 上海联影医疗科技股份有限公司 数据补偿方法、装置、计算机设备和存储介质
CN113030814B (zh) * 2021-03-09 2021-12-21 电子科技大学 一种磁共振k空间数据校正的方法
WO2023241722A1 (en) * 2022-06-17 2023-12-21 Shanghai United Imaging Healthcare Co., Ltd. Methods and systems for image reconstruction

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1627100A (zh) * 2003-12-08 2005-06-15 西门子公司 运行具有改进的线性性的计数辐射检测器的方法
CN104183012A (zh) * 2013-10-31 2014-12-03 上海联影医疗科技有限公司 一种pet三维图像重建方法和装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
PT103200B (pt) * 2004-09-30 2006-08-24 Taguspark-Soc. Prom.Desenv.Parq.Ci.Tec.Area Lisboa Sistema de tomografia por emissão de positrões (pet)
CN100508891C (zh) * 2007-09-04 2009-07-08 陈武凡 在pet成像中最大后验优化图像重建方法
US8193505B2 (en) * 2008-09-29 2012-06-05 Siemens Medical Solutions Usa, Inc. System and method for scatter normalization of PET images
CN102813527A (zh) * 2011-06-10 2012-12-12 北京大基康明医疗设备有限公司 正电子发射计算机断层扫描系统
CN102324089B (zh) * 2011-07-13 2013-04-03 南方医科大学 基于广义熵与mr先验的pet图像最大后验重建方法
US9645261B2 (en) * 2013-01-11 2017-05-09 Siemens Medical Solutions Usa, Inc. Normalization coefficients in PET continuous bed motion acquisition
CN103279964B (zh) * 2013-04-23 2015-10-28 浙江大学 一种基于prca的pet图像动态重建方法及系统
JP6462262B2 (ja) * 2013-08-30 2019-01-30 キヤノンメディカルシステムズ株式会社 Ｘ線コンピュータ断層撮影装置及びフォトンカウンティングｃｔ装置
CN105389811B (zh) 2015-10-30 2017-04-05 吉林大学 一种基于多级阈值分割的多模态医学图像处理方法
CN106251313B (zh) * 2016-08-15 2020-06-26 上海联影医疗科技有限公司 医学成像方法及系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1627100A (zh) * 2003-12-08 2005-06-15 西门子公司 运行具有改进的线性性的计数辐射检测器的方法
CN104183012A (zh) * 2013-10-31 2014-12-03 上海联影医疗科技有限公司 一种pet三维图像重建方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Does partial volume corrected maximum SUV based on count recovery coefficient in 3D-PET/CT correlate with clinical aggressiveness of non-Hodgkin’s lymphoma?;T Tsujikawa, H Otsuka et al.;《Annals of Nuclear Medicine》;20080131;23-30 *
正电子发射断层扫描仪PET中的数据校正常用方法;刘力，王燕芳;《中国体视学与图像分析》;20070630;147-151 *

Also Published As

Publication number Publication date
CN106251313A (zh) 2016-12-21
US20200041669A1 (en) 2020-02-06
US11874412B2 (en) 2024-01-16
US20210389484A1 (en) 2021-12-16
US20190049605A1 (en) 2019-02-14
US11105939B2 (en) 2021-08-31
US10429525B2 (en) 2019-10-01
WO2018033003A1 (en) 2018-02-22

Similar Documents

Publication Publication Date Title
CN106251313B (zh) 医学成像方法及系统
US10692212B2 (en) System and method for image reconstruction
US11331066B2 (en) Imaging method and system for determining a scan area
US11574716B2 (en) System and method for generating imaging report
US20170301066A1 (en) System and method for image correction
CN107330952B (zh) 电子计算机断层扫描图像数据处理系统
CN107330951B (zh) 图像重建系统及方法
US11393139B2 (en) System and method for MPR streak reduction
CN107392876B (zh) 获取归一化校正因子的方法、装置及医学成像方法
CN107909618B (zh) 图像重建系统和方法
CN106485680A (zh) 图像校正方法及装置
US11080896B2 (en) Systems and methods for computed tomography image reconstruction
US11361480B2 (en) System and method for 3D image reconstruction from axial step-and-shoot CT
US10679817B2 (en) Method and system for adjusting focal point position
US11900602B2 (en) System and method for medical imaging
US11576627B2 (en) Systems and methods for data collection in a medical device
CN111080737B (zh) 图像重建方法、装置及pet扫描系统
WO2023087260A1 (en) Systems and methods for data processing
CN117635511A (zh) 医学图像处理方法、装置、计算机设备和存储介质
CN115176280A (zh) 用于基于深度学习的图像衰减校正的方法和装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258

Patentee after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258

Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

CP01 Change in the name or title of a patent holder