CN103654818B - 一种实现多能量ct扫描重建的方法及装置 - Google Patents

一种实现多能量ct扫描重建的方法及装置 Download PDF

Info

Publication number
CN103654818B
CN103654818B CN201210324231.9A CN201210324231A CN103654818B CN 103654818 B CN103654818 B CN 103654818B CN 201210324231 A CN201210324231 A CN 201210324231A CN 103654818 B CN103654818 B CN 103654818B
Authority
CN
China
Prior art keywords
voltage
data
scan
energy
projection
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
CN201210324231.9A
Other languages
English (en)
Other versions
CN103654818A (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.)
Neusoft Medical Systems Co Ltd
Original Assignee
Neusoft Medical Systems 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 Neusoft Medical Systems Co Ltd filed Critical Neusoft Medical Systems Co Ltd
Priority to CN201210324231.9A priority Critical patent/CN103654818B/zh
Priority to US13/974,074 priority patent/US9232927B2/en
Publication of CN103654818A publication Critical patent/CN103654818A/zh
Application granted granted Critical
Publication of CN103654818B publication Critical patent/CN103654818B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • 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 or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating thereof
    • A61B6/582—Calibration
    • A61B6/583—Calibration using calibration phantoms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/408—Dual energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Optics & Photonics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Pulmonology (AREA)

Abstract

本发明公开了一种实现多能量CT数据采集的方法，在数据采集期间以周期性变化的方式对扫描对象进行扫描；还公开了一种实现多能量CT扫描重建的方法及装置，通过控制第一类电压的m个采样点的电压值对第一类扫描对象进行电压周期性变化的数据采集，获得n组多能量第一类扫描数据{yi}，又通过m个采样点的每个电压值分别对包含第一类扫描对象待分割的两种物质的双柱校正模进行数据采集，获得所述m个采样点的电压值分别对应的组合系数，进而得到所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数和又通过计算所述第一类扫描数据{yi}与组合投影数据之间差距的最小值，得到图像向量X第一和X第二，实现了由多能量扫描数据重建出组织分离的图像的目的。

Description

一种实现多能量CT扫描重建的方法及装置
技术领域
本发明涉及CT领域，特别涉及一种实现多能量CT扫描重建的方法及装置。
背景技术
现在，CT扫描被广泛的应用在在医学领域。利用不同能量的x射线在物质中衰减不同，可以实现对CT图像中的不同组织结构进行分离。
目前，不同能量数据的采集方法一般有两种，一种是通过高压交替地快速切换，使得不同方向的投影数据具有不同的电压，从而使得一次扫描中，不同投影方向交替获得不同能量的数据；另外一种是利用双球管技术，每个球管发出的射线具有不同的能量，从而也能利用一次扫描获得两种能量的数据。
对同一种物质，例如水或者骨，两种能量的x射线穿过其中以后的衰减不同，从而得到的图像的像素值也不同，利用这种差别可以将图像中的两种成分分离出来。例如：用两种能量扫描获得的数据进行CT图像重建，得到两幅图像I1和I2，可简单地通过如下2组方程式获得图像中各点水和骨的成分比例：
Water＝Coeff1*I1+Coeff2*I2
Bone＝Coeff3*I1+Coeff4*I2
若组合系数Coeff1，Coeff2，Coeff3，Coeff4已知，则很容易根据两种能量的图像I1和I2计算出水和骨的比例来，而这些系数的获得一般是通过事先扫描简单的校正模体，并通过一定的校正算法来得到。
但是，获取双能的数据对硬件要求较高，比如，上面的两种获取方式中，一种要求高压以很快的速度在不同投影方向间切换，另一种则要求安装两个球管。
发明内容
有鉴于此，本发明的主要目的在于提供一种实现多能量CT数据采集的方法，以解决获取双能的数据的数据采集方法对硬件要求较高的问题。
本发明提供了一种实现多能量CT数据采集的方法，包括：
利用CT扫描仪控制电压在数据采集期间以周期性变化的方式对扫描对象进行扫描，所述电压在数据采集期间内为完整周期的变化。。
优选地，所述控制电压在数据采集期间以周期性变化的方式对扫描对象进行扫描，具体为控制电压在最高电压与最低电压之间的m个不同采样点周期性地对扫描对象进行扫描，其中m≥3。
优选地，利用第一函数f(x)控制电压在数据采集期间以周期性变化的方式扫描，其中所述第一函数f(x)以2T为周期，并且在每个周期内，f(x)能取到最高电压和最低电压。
优选地，所述第一函数f(x)具体为：：
其中，x为整数，2T为电压变化的周期， MaxKV为最高电压，MinKV为最低电压。
本发明还提供一种实现多能量CT扫描重建的方法及装置以实现能够通过电压周期性变化的数据采集方式获得多能量扫描数据、由多能量扫描数据重建出组织分离的图像的目的，从而避免了通过获取双能的数据进行重建组织分离图像对硬件要求较高的问题。
本发明提供了一种实现多能量CT扫描重建的方法，该方法包括第一采集步骤、第二采集步骤、多能量校正步骤和多能量图像重建步骤：
所述第一采集步骤包括：
以由第一物质和第二物质组成的第一类扫描对象作为扫描对象，应用以上所述的一种实现多能量CT数据采集的方法，采用第一类电压进行数据采集，采集到n个投影角度下的n组第一类扫描数据{yi}，其中，i为整数，1≤i≤n，n≥2，所述第一类电压在数据采集期间周期性变化，所述第一类电压在最高电压与最低电压之间共有m个不同的采样点，所述采用第一类电压进行数据采集具体为采用所述m个采样点的电压值在数据采集期间周期性地进行数据采集；
所述第二采集步骤包括：
利用CT扫描仪，控制所述m个采样点中每个采样点的电压值，分别对双柱校正模进行所述n个投影角度的扫描，获得m个n组第二类扫描数据，其中所述双柱校正模中包含所述第一类扫描对象的两种待分割物质，一柱体包含第一物质，另一柱体包含第二物质；
所述多能量校正步骤包括：
利用所述m个n组第二类扫描数据重建并校正出m个第二类扫描图像，其中，所述每个n组第二类扫描数据重建并校正出对应的一个第二类扫描图像；
根据所述第一物质与第二物质在双柱校正模中的位置范围，得到所述第一物质的第三类图像向量和第二物质的第三类图像向量；
根据m组方程式计算m个采样点的电压值分别对应的组合系数与
所述多能量图像重建步骤包括：
通过计算所述第一类扫描数据{yi}与组合投影数据之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二；
其中，和为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数与获得；
为投影矩阵的第i行；
将满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二分别作为第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
本发明还提供一种实现多能量CT扫描重建的装置，包括：
第一采集单元：用于以由第一物质和第二物质组成的第一类扫描对象作为扫描对象，应用以上所述的一种实现多能量CT数据采集的方法，采用第一类电压进行数据采集，采集到n个投影角度下的n组第一类扫描数据{yi}，其中，i为整数，1≤i≤n，n≥2，所述第一类电压在数据采集期间周期性变化；所述第一类电压在最高电压与最低电压之间共有m个不同的采样点，所述采用第一类电压进行数据采集具体为采用所述m个采样点的电压值在数据采集期间周期性地进行数据采集，将所述n组第一类扫描数据{yi}发送至多能量图像重建单元；
第二采集单元：利用CT扫描仪，控制所述m个采样点中每个采样点的电压值，分别对双柱校正模进行所述n个投影角度的扫描，获得m个n组第二类扫描数据，其中所述双柱校正模中包含所述第一类扫描对象的两种待分割物质，一柱体包含第一物质，另一柱体包含第二物质，将所述m个n组第二类扫描数据发送至多能量校正单元；
多能量校正单元：用于利用所述m个n组第二类扫描数据重建并校正出m个第二类扫描图像，其中，所述每个n组第二类扫描数据重建并校正出对应的一个第二类扫描图像，根据所述第一物质与第二物质在双柱校正模中的位置范围，得到所述第一物质的第三类图像向量和第二物质的第三类图像向量，根据m组方程式计算出m个采样点的电压值分别对应的组合系数与将所述m个采样点的电压值分别对应的组合系数与发送至多能量图像重建单元；
多能量图像重建单元：用于通过计算所述第一类扫描数据{yi}与组合投影数据之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二，其中，和为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数与获得；为投影矩阵的第i行；将满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
可见本发明具有如下有益效果：
由于本发明提供的一种实现多能量CT数据采集的方法在数据采集期间以周期性变化的方式对扫描对象进行扫描，避免了获取双能的数据的数据采集方法对硬件要求较高的问题；而且本发明还提供一种实现多能量CT扫描重建的方法及装置，通过控制第一类电压的m个采样点的电压值对由第一物质和第二物质组成的第一类扫描对象进行电压周期性变化的数据采集，获得n组多能量第一类扫描数据{yi}，又通过控制m个采样点的每个电压值分别对包含第一类扫描对象待分割的两种物质的双柱校正模进行数据采集，获得所述m个采样点的电压值分别对应的组合系数与进而得到所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数和又通过计算所述第一类扫描数据{yi}与组合投影数据之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二，将满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置，可见，利用本发明的多能量CT扫描重建方法及装置能够通过电压周期性变化的方式获得多能量扫描数据、由多能量扫描数据重建出组织分离的图像，即第一物质的第一类图像向量和第二物质的第一类图像向量，投影之间的能量切换既可以在高电压和低电压之间快速进行，又可以通过逐渐变化的方式进行，使得电压从低压到高压的变化可以跨过多个投影，避免了获取双能的数据对硬件要求较高的问题。
附图说明
图1是本发明一种实现多能量CT扫描重建的方法步骤图例；
图2是本发明双柱校正模示意图例；
图3是本发明将X第一和X第二的进行循环迭代的步骤图例；
图4是本发明一种实现多能量CT扫描重建的装置的组成图例。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明实施例作进一步详细的说明。
本发明提供了一种实现多能量CT数据采集的方法，包括：
利用CT扫描仪控制电压在数据采集期间以周期性变化的方式对扫描对象进行扫描，所述电压在数据采集期间内为完整周期的变化。
在本发明的一个优选实施例中，所述控制电压在数据采集期间以周期性变化的方式对扫描对象进行扫描，具体为控制电压在最高电压与最低电压之间的m个不同采样点周期性地对扫描对象进行扫描，其中m≥3。
在本发明的又一优选实施例中，提出利用第一函数f(x)控制电压在数据采集期间以周期性变化的方式扫描，其中所述第一函数f(x)以2T为周期，并且在每个周期内，f(x)能取到最高电压和最低电压，具体地，所述第一函数f(x)在一个优选实施例中具体为：
其中，x为整数，2T为电压变化的周期， MaxKV为最高电压，MinKV为最低电压。
本发明提供了一种实现多能量CT扫描重建的方法，参见图1，该方法包括第一采集步骤、第二采集步骤、多能量校正步骤和多能量图像重建步骤：
其中，所述第一采集步骤包括：
S101、以由第一物质和第二物质组成的第一类扫描对象作为扫描对象，应用如以上所述的一种实现多能量CT数据采集的方法，采用第一类电压进行数据采集，采集到n个投影角度下的n组第一类扫描数据{yi}，其中，i为整数，1≤i≤n，n≥2，所述第一类电压在数据采集期间周期性变化，所述第一类电压在最高电压与最低电压之间共有m个不同的采样点，所述采用第一类电压进行数据采集具体为采用所述m个采样点的电压值在数据采集期间周期性地进行数据采集；
该步骤为对第一类扫描对象的数据采集过程，假设采集一圈的投影个数为n，采集到的n组第一类扫描数据{yi}为多能量数据，i表示第i个投影角度；yi表示利用CT扫描仪对第一类扫描对象应用如以上所述的一种实现多能量CT数据采集的方法以电压周期性变化的方式进行数据采集，采集到的第i个投影角度的第一类扫描数据；
所述第二采集步骤包括：
S102、利用CT扫描仪，控制所述m个采样点中每个采样点的电压值，分别对双柱校正模进行所述n个投影角度的扫描，获得m个n组第二类扫描数据，其中所述双柱校正模中包含所述第一类扫描对象的两种待分割物质，一柱体包含第一物质，另一柱体包含第二物质；
在该步骤中，利用CT扫描仪控制所述m个采样点中每个采样点的电压值，分别对双柱校正模进行所述n个投影角度的扫描，也就是分别用所述m个采样点中每个采样点的电压值各自对双柱校正模进行一圈的数据采集，即分别用所述m个采样点中每个采样点的电压值各自进行所述n个投影角度的扫描，以获得m个n组第二类扫描数据；
所述双柱校正模，具体可以采用如图2所示的双柱校正模体，其中的两个柱体的材料分别为将来要分割出的物质，例如：水和骨；
所述多能量校正步骤包括：
S103、利用所述m个n组第二类扫描数据重建并校正出m个第二类扫描图像，其中，所述每个n组第二类扫描数据重建并校正出对应的一个第二类扫描图像；
在该步骤中，所述每个n组第二类扫描数据重建并校正出对应的一个第二类扫描图像，具体可以采用本领域现有的重建、硬化校正方法，例如卷积反投影法、双能校正法；
S104、根据所述第一物质与第二物质在双柱校正模中的位置范围，得到所述第一物质的第三类图像向量和第二物质的第三类图像向量；
在该步骤中，可以根据第一物质和第二物质的几何位置和尺寸计算出图像中哪些像素是第一物质，哪些像素是第二物质，也可将步骤S103获得的任意一幅第二类扫描图像简单阈值分割得到第一物质的第三类图像向量和第二物质的第三类图像向量，其中，在第一物质的第三类图像中，除第一物质为1外，其他点都为0；在第二物质的第三类图像中，除第二物质为1外，其他点都为0；
S105、根据m组方程式计算出m个采样点的电压值分别对应的组合系数与
在该步骤中，具体可以采用最小二乘法计算出组合系数与具体地，为利用m个采样点中第H个采样点电压值进行数据采集时，第一物质对应的组合系数，为利用m个采样点中第H个采样点电压值时，第二物质对应的组合系数；
所述多能量图像重建步骤包括：
S106、通过计算所述第一类扫描数据{yi}与组合投影数据之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二；
其中，和为所述第一类电压在采集第i个投影角度的第一类扫描数据yi时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数与获得；
为投影矩阵的第i行；
需要说明的是：
(1)因为和为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，为第一物质对应的组合系数，为第二物质对应的组合系数，而且，所述第一类电压进行数据采集具体为采用所述m个采样点的电压值在数据采集期间周期性地进行数据采集，因此，和可以具体根据所述m个采样点的电压值分别对应的组合系数与获得，例如：
m个采样点的电压值分别对应的组合系数包括到 包括到
那么，可以得到在数据采集期间，所述第一类电压在采集第一个变化周期内所跨投影角度的第一类扫描数据时，所对应的组合系数为：
第一物质对应的组合系数为：到 到
第二物质对应的组合系数为：到 到
依次重复，可以得知所述第一类电压在采集后续变化周期内所跨投影角度的第一类扫描数据时对应的组合系数，进而得到和
(2)所述投影矩阵P可通过本领域现有的方法计算获得，例如，利用第一类扫描对象的医学图像各个像素的坐标、x射线光源的坐标、CT扫描仪各个检测器的坐标之间的几何关系计算出光源发出的x射线经过某个像素的长度获得所述投影矩阵P。
S107、将满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
通过以上具体实施方式可见，本发明提出的实现多能量CT扫描重建的方法，在投影之间的能量切换既可以在高电压和低电压之间快速进行，又可以通过逐渐变化的方式进行，使得电压从低压到高压的变化可以跨过多个投影，避免了获取双能的数据对硬件要求较高的问题。
优选地，在本发明的一个实施例中，投影之间的能量以逐渐变化的方式进行，在该实施例中，所述第一类电压在最高电压与最低电压之间m个不同采样点的个数m≥3。
下面对所述多能量图像重建步骤中所述计算第一类扫描数据{yi}与组合投影数据之间差距的最小值的步骤，进行详细说明，需要说明的是，以下所述多能量图像重建步骤既适用于投影间的能量高低压快速切换的情况，又适用于投影间的能量逐渐变化的情况：
在本发明的该优选实施例中，所述计算第一类扫描数据{yi}与组合投影数据之间差距的最小值具体通过计算目标函数的最小值实现，其中，是差距项，衡量{yi}与之间的差距，β(||X第一||TV+||X第二||TV)为正则项，β为折衷参数，合适的折衷参数β使得求出的目标图像向量X第一和X第二满足一定的光滑性，β越大，最终求得的图像越平滑，分辨率越低，β越小，最终求得的图像噪声越大，分辨率越高，在本发明具体实施时，应该根据需求，适当地调整β，使得图像能满足清晰度需求；||·||TV表示全变差。
在本发明的一个优选实施例中，所述计算目标函数的最小值具体通过对图像向量X第一和X第二进行循环迭代实现，参见图3，具体包括：
S301、首先设置迭代的初始条件：
以图像向量表示X第一；
以图像向量表示X第二；
定义r＝0，和为初始图像，具体可以是0图像；
利用投影矩阵P对进行正投影，获得初始的正投影数据利用投影矩阵P对进行正投影，获得初始的正投影数据
S302、从r＝0开始，进入迭代步骤：
若输出条件未满足，则重复以下迭代步骤，其中，本次迭代为第r+1次迭代，表示该第r+1次迭代得到的图像向量X第一，表示该第r+1次迭代前的图像向量X第一，表示该第r+1次迭代得到的图像向量X第二，表示该第r+1次迭代前的图像向量X第二，迭代步骤每执行完一次将r的值加1，所述迭代步骤包括图像向量的更新和图像向量的更新：
S302.1、所述的更新步骤包括：
重复以下步骤，直到所有的像素点都被更新：
S302.1A、选择图像向量中的任意一像素点且该像素点在该第r+1次迭代步骤中未被更新；
S302.1B、利用函数
将代入后，得到ψ关于的偏导等于0的方程：
S302.1C、假设：
其中Pij为投影矩阵P在第i行第j列上的投影系数；
根据等于投影数据的第i行，等于投影数据的第i行，计算获得与
S302.1D、根据利用二分法计算出完成像素点的更新；
在这里，计算出完成像素点的更新，即得到更新的图像向量
S302.1E、更新正投影数据
S302.2、所述的更新步骤包括：
重复以下步骤，直到所有的像素点都被更新：
S302.2A、选择图像向量中的任意一像素点且该像素点在该第r+1次迭代步骤中未被更新；
S302.2B、利用函数
将代入后，得到ψ关于的偏导等于0的方程：
S302.2C、假设：
其中Pij为投影矩阵P在第i行第j列上的投影系数；
根据等于投影数据的第i行，等于投影数据的第i行，计算获得与
S302.2D、根据利用二分法计算出完成像素点的更新；
在这里，计算出完成像素点的更新，即得到更新的图像向量
S302.2E、更新正投影数据
需要说明的是，所述迭代步骤循环的条件为输出条件未满足，所述输出条件具体可以设置为且其中ε为预设的小阈值，在理论上该预设的小阈值ε越小，图像质量越好，在本发明中，经过多次试验，该预设的小阈值ε至少应小于10-3，迭代优化最终得到的与图像质量较好。
本发明还包括一种实现多能量CT扫描重建的装置，参见图4，该装置包括：
第一采集单元401：用于以由第一物质和第二物质组成的第一类扫描对象作为扫描对象，应用如以上所述的一种实现多能量CT数据采集的方法，采用第一类电压进行数据采集，采集到n个投影角度下的n组第一类扫描数据{yi}，其中，i为整数，1≤i≤n，n≥2，所述第一类电压在数据采集期间周期性变化；所述第一类电压在最高电压与最低电压之间共有m个不同的采样点，所述采用第一类电压进行数据采集具体为采用所述m个采样点的电压值在数据采集期间周期性地进行数据采集，将所述n组第一类扫描数据{yi}发送至多能量图像重建单元404；
第二采集单元402：利用CT扫描仪，控制所述m个采样点中每个采样点的电压值，分别对双柱校正模进行所述n个投影角度的扫描，获得m个n组第二类扫描数据，其中所述双柱校正模中包含所述第一类扫描对象的两种待分割物质，一柱体包含第一物质，另一柱体包含第二物质，将所述m个n组第二类扫描数据发送至多能量校正单元403；
多能量校正单元403：用于利用所述m个n组第二类扫描数据重建并校正出m个第二类扫描图像，其中，所述每个n组第二类扫描数据重建并校正出对应的一个第二类扫描图像，根据所述第一物质与第二物质在双柱校正模中的位置范围，得到所述第一物质的第三类图像向量和第二物质的第三类图像向量，根据m组方程式计算出m个采样点的电压值分别对应的组合系数与将所述m个采样点的电压值分别对应的组合系数与发送至多能量图像重建单元404；
多能量图像重建单元404：用于通过计算所述第一类扫描数据{yi}与组合投影数据之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二，其中，和为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数与获得；为投影矩阵的第i行，为投影矩阵的第i行；将满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个......”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均包含在本发明的保护范围内。

Claims (13)

1.一种实现多能量CT扫描重建的方法，其特征在于，包括第一采集步骤、第二采集步骤、多能量校正步骤和多能量图像重建步骤：
所述第一采集步骤包括：
以由第一物质和第二物质组成的第一类扫描对象作为扫描对象，应用实现多能量CT数据采集的方法，采用第一类电压进行数据采集，采集到n个投影角度下的n组第一类扫描数据{yi}，其中，i为整数，1≤i≤n，n≥2，所述第一类电压在数据采集期间周期性变化，所述第一类电压在最高电压与最低电压之间共有m个不同的采样点，所述采用第一类电压进行数据采集具体为采用所述m个采样点的电压值在数据采集期间周期性地进行数据采集；其中，所述实现多能量CT数据采集的方法包括：利用CT扫描仪控制所述第一类电压在数据采集期间以周期性变化的方式对扫描对象进行扫描，所述第一类电压在数据采集期间内为完整周期的变化；
所述第二采集步骤包括：
利用CT扫描仪，控制所述m个采样点中每个采样点的电压值，分别对双柱校正模进行所述n个投影角度的扫描，获得m个n组第二类扫描数据，其中所述双柱校正模中包含所述第一类扫描对象的两种待分割物质，一柱体包含第一物质，另一柱体包含第二物质；
所述多能量校正步骤包括：
利用所述m个n组第二类扫描数据重建并校正出m个第二类扫描图像，其中，所述每个n组第二类扫描数据重建并校正出对应的一个第二类扫描图像；
根据所述第一物质与第二物质在双柱校正模中的位置范围，得到所述第一物质的第三类图像向量和第二物质的第三类图像向量；
根据m组方程式
计算出m个采样点的电压值分别对应的组合系数与
所述多能量图像重建步骤包括：
通过计算所述第一类扫描数据{yi}与组合投影数据之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二；
其中，和为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数与获得；
为投影矩阵的第i行；
将满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
2.根据权利要求1所述的方法，其特征在于，所述第一类电压在最高电压与最低电压之间m个不同采样点的个数m≥3。
3.根据权利要求1或2所述的方法，其特征在于，所述计算第一类扫描数据{yi}与组合投影数据之间差距的最小值，具体通过计算目标函数的最小值实现，其中，β为折衷参数。
4.根据权利要求3所述的方法，其特征在于，所述计算目标函数的最小值具体包括以下步骤：
以图像向量表示X第一；
以图像向量表示X第二；
定义r＝0，和为初始图像；
利用投影矩阵P对进行正投影，获得初始的正投影数据利用投影矩阵P对进行正投影，获得初始的正投影数据
若输出条件未满足，则重复以下迭代步骤，其中，本次迭代为第r+1次迭代，表示该第r+1次迭代得到的图像向量X第一，表示该第r+1次迭代前的图像向量X第一，表示该第r+1次迭代得到的图像向量X第二，表示该第r+1次迭代前的图像向量X第二，迭代步骤每执行完一次将r的值加1，所述迭代步骤包括图像向量的更新和图像向量的更新：
所述的更新步骤包括：
重复以下步骤，直到所有的像素点都被更新：
选择图像向量中的任意一像素点且该像素点在该第r+1次迭代步骤中未被更新；
利用函数将代入后，得到ψ关于的偏导等于0的方程：
定义：其中Pij为投影矩阵P在第i行第j列上的投影系数；
根据等于投影数据的第i行，等于投影数据的第i行，计算获得与
根据利用二分法计算出完成像素点的更新；
更新正投影数据
所述的更新步骤包括：
重复以下步骤，直到所有的像素点都被更新：
选择图像向量中的任意一像素点且该像素点在该第r+1次迭代步骤中未被更新；
利用函数将代入后，得到ψ关于的偏导等于0的方程：
定义：
根据等于投影数据的第i行，等于投影数据的第i行，计算获得与
根据利用二分法计算出完成像素点的更新
更新正投影数据
5.根据权利要求4所述的方法，其特征在于，所述输出条件具体为且其中ε为预设的小阈值。
6.根据权利要求5所述的方法，其特征在于，所述预设的小阈值ε至少小于10-3。
7.根据权利要求1所述的方法，其特征在于，所述根据m组方程式计算出m个采样点的电压值分别对应的组合系数与具体采用最小二分法进行计算。
8.根据权利要求1所述的方法，其特征在于，利用第一函数f(x)控制所述第一类电压在数据采集期间以周期性变化的方式扫描，其中所述第一函数f(x)以2T为周期，并且在每个周期内，f(x)能取到最高电压和最低电压。
9.根据权利要求8所述的方法，其特征在于，所述第一函数f(x)具体为：
其中，x为整数，2T为所述第一类电压变化的周期，MaxKV为最高电压，MinKV为最低电压。
10.一种实现多能量CT扫描重建的装置，其特征在于，包括：
第一采集单元：用于以由第一物质和第二物质组成的第一类扫描对象作为扫描对象，应用实现多能量CT数据采集的方法，采用第一类电压进行数据采集，采集到n个投影角度下的n组第一类扫描数据{yi}，其中，i为整数，1≤i≤n，n≥2，所述第一类电压在数据采集期间周期性变化；所述第一类电压在最高电压与最低电压之间共有m个不同的采样点，所述采用第一类电压进行数据采集具体为采用所述m个采样点的电压值在数据采集期间周期性地进行数据采集，将所述n组第一类扫描数据{yi}发送至多能量图像重建单元；其中，所述实现多能量CT数据采集的方法包括：利用CT扫描仪控制所述第一类电压在数据采集期间以周期性变化的方式对扫描对象进行扫描，所述第一类电压在数据采集期间内为完整周期的变化；
第二采集单元：利用CT扫描仪，控制所述m个采样点中每个采样点的电压值，分别对双柱校正模进行所述n个投影角度的扫描，获得m个n组第二类扫描数据，其中所述双柱校正模中包含所述第一类扫描对象的两种待分割物质，一柱体包含第一物质，另一柱体包含第二物质，将所述m个n组第二类扫描数据发送至多能量校正单元；
多能量校正单元：用于利用所述m个n组第二类扫描数据重建并校正出m个第二类扫描图像，其中，所述每个n组第二类扫描数据重建并校正出对应的一个第二类扫描图像，根据所述第一物质与第二物质在双柱校正模中的位置范围，得到所述第一物质的第三类图像向量和第二物质的第三类图像向量，根据m组方程式计算出m个采样点的电压值分别对应的组合系数与将所述m个采样点的电压值分别对应的组合系数与发送至多能量图像重建单元；
多能量图像重建单元：用于通过计算所述第一类扫描数据{yi}与组合投影数据之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二，其中，和为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数与获得；为投影矩阵的第i行；将满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
11.根据权利要求10所述的装置，其特征在于，所述第一类电压在最高电压与最低电压之间m个不同采样点的个数m≥3。
12.根据权利要求10所述的装置，其特征在于，利用第一函数f(x)控制所述第一类电压在数据采集期间以周期性变化的方式扫描，其中所述第一函数f(x)以2T为周期，并且在每个周期内，f(x)能取到最高电压和最低电压。
13.根据权利要求12所述的装置，其特征在于，所述第一函数f(x)具体为：
其中，x为整数，2T为所述第一类电压变化的周期，MaxKV为最高电压，MinKV为最低电压。
CN201210324231.9A 2012-09-04 2012-09-04 一种实现多能量ct扫描重建的方法及装置 Active CN103654818B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201210324231.9A CN103654818B (zh) 2012-09-04 2012-09-04 一种实现多能量ct扫描重建的方法及装置
US13/974,074 US9232927B2 (en) 2012-09-04 2013-08-23 Method of reconstruction from multi-energy CT scan and device thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210324231.9A CN103654818B (zh) 2012-09-04 2012-09-04 一种实现多能量ct扫描重建的方法及装置

Publications (2)

Publication Number Publication Date
CN103654818A CN103654818A (zh) 2014-03-26
CN103654818B true CN103654818B (zh) 2016-01-27

Family

ID=50187615

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210324231.9A Active CN103654818B (zh) 2012-09-04 2012-09-04 一种实现多能量ct扫描重建的方法及装置

Country Status (2)

Country Link
US (1) US9232927B2 (zh)
CN (1) CN103654818B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105758873B (zh) * 2015-03-04 2019-12-31 公安部第一研究所 一种ct探测装置及其数据处理方法
CN105212958B (zh) * 2015-10-21 2018-08-21 沈阳东软医疗系统有限公司 Ct图像的获得方法、装置及设备
CN106651753B (zh) * 2016-09-28 2020-03-17 东软医疗系统股份有限公司 提高ct图像显示效果的方法及装置
CN106780654B (zh) * 2017-01-24 2021-05-07 东软医疗系统股份有限公司 一种图像重建方法和装置
CN107884806B (zh) * 2017-10-19 2020-05-08 天津大学 面向双能ct成像的x射线能谱探测及重构解析方法
CN112697821B (zh) * 2020-12-02 2022-12-02 赛诺威盛科技(北京)股份有限公司 多能谱ct扫描方法、装置、电子设备和ct设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101292877A (zh) * 2007-04-23 2008-10-29 Ge医疗系统环球技术有限公司 X射线ct设备
CN101772324A (zh) * 2007-08-06 2010-07-07 株式会社日立医药 X射线ct装置
CN103202707A (zh) * 2012-01-16 2013-07-17 株式会社东芝 X射线ct装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8515147B2 (en) * 2006-12-15 2013-08-20 Koninklijke Philips N.V. Spectrally resolving X-ray imaging device
JP2009000293A (ja) * 2007-06-21 2009-01-08 Toshiba Corp デュアルエネルギシステム及びその画像収集方法
EP2243020B1 (en) * 2008-02-15 2017-09-06 Mayo Foundation For Medical Education And Research System and method for quantitative imaging of chemical composition to decompose more than two materials
US8818058B2 (en) * 2009-11-30 2014-08-26 Toshiba Medical Systems Corporation Method for determining a correction function for correcting computed tomographic numbers of a small target object in a CT image

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101292877A (zh) * 2007-04-23 2008-10-29 Ge医疗系统环球技术有限公司 X射线ct设备
CN101772324A (zh) * 2007-08-06 2010-07-07 株式会社日立医药 X射线ct装置
CN103202707A (zh) * 2012-01-16 2013-07-17 株式会社东芝 X射线ct装置

Also Published As

Publication number Publication date
CN103654818A (zh) 2014-03-26
US20140064441A1 (en) 2014-03-06
US9232927B2 (en) 2016-01-12

Similar Documents

Publication Publication Date Title
CN103654818B (zh) 一种实现多能量ct扫描重建的方法及装置
CN103186883B (zh) 一种ct图像重建中骨硬化伪影的校正方法
CN103413280B (zh) 一种低剂量x射线ct图像重建方法
CN103900931B (zh) 一种多能谱ct成像方法及成像系统
CN102254310B (zh) 用于减少ct图像数据中的图像伪影、尤其是金属伪影的方法
US20150228092A1 (en) Digital breast tomosynthesis reconstruction using adaptive voxel grid
CN104318536A (zh) Ct图像的校正方法及装置
CN104597061B (zh) 基于虚拟探测器的大视场锥束ct成像方法
Wu et al. Spatial-temporal total variation regularization (STTVR) for 4D-CT reconstruction
CN102946807B (zh) X射线ct装置及其控制方法
CN104751502B (zh) 一种用于扩大视野的ct图像重建方法
CN103622717B (zh) 由单源单次扫描x光ct图像生成双能x光ct图像方法
CN105225208A (zh) 一种计算机断层成像金属伪影校正方法及装置
WO2015200300A1 (en) Techniques for suppression of motion artifacts in medical imaging
CN103308537B (zh) 一种递变能量x射线成像图像融合方法
CN105488826A (zh) 一种基于fbp的能谱ct迭代成像方法和成像系统
Yang et al. Shading correction assisted iterative cone-beam CT reconstruction
CN103793890A (zh) 一种能谱ct图像的恢复处理方法
KR102348139B1 (ko) 이중 해상도의 관심 영역 내외 투영 데이터를 이용한 체내 단층 촬영 방법 및 시스템
Li et al. Clinical micro-CT empowered by interior tomography, robotic scanning, and deep learning
CN110428478A (zh) 交替光源扇束x射线ct采样方法及装置
CN103310471A (zh) Ct 图像生成装置及方法、ct 图像生成系统
Eldib et al. A motion artifact reduction method for dental CT based on subpixel-resolution image registration of projection data
CN110458908A (zh) 基于有限角度迭代重建超视野ct图像的方法
CN109523602A (zh) 迭代重建优化方法、装置和计算机可读介质

Legal Events

Date Code Title Description
PB01 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
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province

Patentee after: DongSoft Medical System Co., Ltd.

Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16

Patentee before: Dongruan Medical Systems Co., Ltd., Shenyang