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

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

Info

Publication number
CN103654818A
CN103654818A CN201210324231.9A CN201210324231A CN103654818A CN 103654818 A CN103654818 A CN 103654818A CN 201210324231 A CN201210324231 A CN 201210324231A CN 103654818 A CN103654818 A CN 103654818A
Authority
CN
China
Prior art keywords
data
voltage
scan
projection
image vector
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
CN201210324231.9A
Other languages
English (en)
Other versions
CN103654818B (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

Images

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个投影角度的第一类扫描数据时对应的组合系数
Figure DDA00002098257000011
和
Figure DDA00002098257000012
又通过计算所述第一类扫描数据{yi}与组合投影数据
Figure DDA00002098257000013
之间差距的最小值，得到图像向量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)具体为：：
Figure BDA00002098256800021
其中，x为整数，2T为电压变化的周期，
Figure BDA00002098256800023
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个采样点的电压值分别对应的组合系数
Figure BDA00002098256800032
与
Figure BDA00002098256800033
所述多能量图像重建步骤包括：
通过计算所述第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800034
之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800035
之间差距为最小值的图像向量X第一和X第二；
其中，
Figure BDA00002098256800036
和
Figure BDA00002098256800037
为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数
Figure BDA00002098256800038
与
Figure BDA00002098256800039
获得；
Figure BDA000020982568000310
为投影矩阵的第i行；
将满足所述第一类扫描数据{yi}与组合投影数据
Figure BDA000020982568000311
之间差距为最小值的图像向量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组方程式
Figure BDA00002098256800041
计算出m个采样点的电压值分别对应的组合系数
Figure BDA00002098256800042
与
Figure BDA00002098256800043
将所述m个采样点的电压值分别对应的组合系数
Figure BDA00002098256800044
与
Figure BDA00002098256800045
发送至多能量图像重建单元；
多能量图像重建单元：用于通过计算所述第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800046
之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800047
之间差距为最小值的图像向量X第一和X第二，其中，和
Figure BDA00002098256800049
为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数与
Figure BDA00002098256800052
获得；
Figure BDA00002098256800053
为投影矩阵的第i行；将满足所述第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800054
之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
可见本发明具有如下有益效果：
由于本发明提供的一种实现多能量CT数据采集的方法在数据采集期间以周期性变化的方式对扫描对象进行扫描，避免了获取双能的数据的数据采集方法对硬件要求较高的问题；而且本发明还提供一种实现多能量CT扫描重建的方法及装置，通过控制第一类电压的m个采样点的电压值对由第一物质和第二物质组成的第一类扫描对象进行电压周期性变化的数据采集，获得n组多能量第一类扫描数据{yi}，又通过控制m个采样点的每个电压值分别对包含第一类扫描对象待分割的两种物质的双柱校正模进行数据采集，获得所述m个采样点的电压值分别对应的组合系数与
Figure BDA00002098256800056
进而得到所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数
Figure BDA00002098256800057
和
Figure BDA00002098256800058
又通过计算所述第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800059
之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据
Figure BDA000020982568000510
之间差距为最小值的图像向量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为电压变化的周期，
Figure BDA00002098256800073
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组方程式
Figure BDA00002098256800081
计算出m个采样点的电压值分别对应的组合系数
Figure BDA00002098256800082
与
Figure BDA00002098256800083
在该步骤中，具体可以采用最小二乘法计算出组合系数与
Figure BDA00002098256800085
具体地，
Figure BDA00002098256800086
为利用m个采样点中第H个采样点电压值进行数据采集时，第一物质对应的组合系数，
Figure BDA00002098256800087
为利用m个采样点中第H个采样点电压值时，第二物质对应的组合系数；
所述多能量图像重建步骤包括：
S106、通过计算所述第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800091
之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800092
之间差距为最小值的图像向量X第一和X第二；
其中，和
Figure BDA00002098256800094
为所述第一类电压在采集第i个投影角度的第一类扫描数据yi时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数
Figure BDA00002098256800095
与
Figure BDA00002098256800096
获得；
为投影矩阵的第i行；
需要说明的是：
(1)因为
Figure BDA00002098256800098
和为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，
Figure BDA000020982568000910
为第一物质对应的组合系数，
Figure BDA000020982568000911
为第二物质对应的组合系数，而且，所述第一类电压进行数据采集具体为采用所述m个采样点的电压值在数据采集期间周期性地进行数据采集，因此，
Figure BDA000020982568000912
和
Figure BDA000020982568000913
可以具体根据所述m个采样点的电压值分别对应的组合系数
Figure BDA000020982568000914
与
Figure BDA000020982568000915
获得，例如：
m个采样点的电压值分别对应的组合系数
Figure BDA000020982568000916
包括到
Figure BDA000020982568000918
包括
Figure BDA000020982568000920
到
Figure BDA000020982568000921
那么，可以得到在数据采集期间，所述第一类电压在采集第一个变化周期内所跨投影角度的第一类扫描数据时，所对应的组合系数为：
第一物质对应的组合系数为：
Figure BDA00002098256800101
到
Figure BDA00002098256800102
Figure BDA00002098256800103
到
Figure BDA00002098256800104
第二物质对应的组合系数为：
Figure BDA00002098256800105
到
Figure BDA00002098256800106
Figure BDA00002098256800107
到
Figure BDA00002098256800108
依次重复，可以得知所述第一类电压在采集后续变化周期内所跨投影角度的第一类扫描数据时对应的组合系数，进而得到
Figure BDA00002098256800109
和
Figure BDA000020982568001010
(2)所述投影矩阵P可通过本领域现有的方法计算获得，例如，利用第一类扫描对象的医学图像各个像素的坐标、x射线光源的坐标、CT扫描仪各个检测器的坐标之间的几何关系计算出光源发出的x射线经过某个像素的长度获得所述投影矩阵P。
S107、将满足所述第一类扫描数据{yi}与组合投影数据
Figure BDA000020982568001011
之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
通过以上具体实施方式可见，本发明提出的实现多能量CT扫描重建的方法，在投影之间的能量切换既可以在高电压和低电压之间快速进行，又可以通过逐渐变化的方式进行，使得电压从低压到高压的变化可以跨过多个投影，避免了获取双能的数据对硬件要求较高的问题。
优选地，在本发明的一个实施例中，投影之间的能量以逐渐变化的方式进行，在该实施例中，所述第一类电压在最高电压与最低电压之间m个不同采样点的个数m≥3。
下面对所述多能量图像重建步骤中所述计算第一类扫描数据{yi}与组合投影数据
Figure BDA000020982568001012
之间差距的最小值的步骤，进行详细说明，需要说明的是，以下所述多能量图像重建步骤既适用于投影间的能量高低压快速切换的情况，又适用于投影间的能量逐渐变化的情况：
在本发明的该优选实施例中，所述计算第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800111
之间差距的最小值具体通过计算目标函数
Figure BDA00002098256800112
的最小值实现，其中，
Figure BDA00002098256800113
是差距项，衡量{yi}与
Figure BDA00002098256800114
之间的差距，β(||X第一||TV+||X第二||TV)为正则项，β为折衷参数，合适的折衷参数β使得求出的目标图像向量X第一和X第二满足一定的光滑性，β越大，最终求得的图像越平滑，分辨率越低，β越小，最终求得的图像噪声越大，分辨率越高，在本发明具体实施时，应该根据需求，适当地调整β，使得图像能满足清晰度需求；||·||TV表示全变差。
在本发明的一个优选实施例中，所述计算目标函数
Figure BDA00002098256800115
的最小值具体通过对图像向量X第一和X第二进行循环迭代实现，参见图3，具体包括：
S301、首先设置迭代的初始条件：
以图像向量
Figure BDA00002098256800116
表示X第一；
以图像向量
Figure BDA00002098256800117
表示X第二；
定义r＝0，和
Figure BDA00002098256800119
为初始图像，具体可以是0图像；
利用投影矩阵P对
Figure BDA000020982568001110
进行正投影，获得初始的正投影数据
Figure BDA000020982568001111
利用投影矩阵P对
Figure BDA000020982568001112
进行正投影，获得初始的正投影数据
Figure BDA000020982568001113
S302、从r＝0开始，进入迭代步骤：
若输出条件未满足，则重复以下迭代步骤，其中，本次迭代为第r+1次迭代，
Figure BDA000020982568001114
表示该第r+1次迭代得到的图像向量X第一，
Figure BDA000020982568001115
表示该第r+1次迭代前的图像向量X第一，
Figure BDA00002098256800121
表示该第r+1次迭代得到的图像向量X第二，
Figure BDA00002098256800122
表示该第r+1次迭代前的图像向量X第二，迭代步骤每执行完一次将r的值加1，所述迭代步骤包括图像向量
Figure BDA00002098256800123
的更新和图像向量
Figure BDA00002098256800124
的更新：
S302.1、所述的更新步骤包括：
重复以下步骤，直到
Figure BDA00002098256800126
所有的像素点都被更新：
S302.1A、选择图像向量
Figure BDA00002098256800127
中的任意一像素点
Figure BDA00002098256800128
且该像素点在该第r+1次迭代步骤中未被更新；
S302.1B、利用函数
将
Figure BDA000020982568001210
代入后，得到ψ关于
Figure BDA000020982568001211
的偏导等于0的方程：
Figure BDA000020982568001212
S302.1C、假设：
Figure BDA000020982568001213
Figure BDA000020982568001214
其中Pij为投影矩阵P在第i行第j列上的投影系数；
根据
Figure BDA000020982568001215
等于投影数据
Figure BDA000020982568001216
的第i行，
Figure BDA000020982568001217
等于投影数据
Figure BDA000020982568001218
的第i行，计算获得
Figure BDA000020982568001219
与
Figure BDA000020982568001220
S302.1D、根据
Figure BDA000020982568001221
利用二分法计算出
Figure BDA000020982568001222
完成像素点
Figure BDA000020982568001223
的更新；
在这里，计算出
Figure BDA000020982568001224
完成像素点
Figure BDA000020982568001225
的更新，即得到更新的图像向量
Figure BDA000020982568001226
S302.1E、更新正投影数据
Figure BDA00002098256800131
S302.2、所述
Figure BDA00002098256800132
的更新步骤包括：
重复以下步骤，直到
Figure BDA00002098256800133
所有的像素点都被更新：
S302.2A、选择图像向量
Figure BDA00002098256800134
中的任意一像素点
Figure BDA00002098256800135
且该像素点在该第r+1次迭代步骤中未被更新；
S302.2B、利用函数
Figure BDA00002098256800136
将代入后，得到ψ关于
Figure BDA00002098256800138
的偏导等于0的方程：
S302.2C、假设：
Figure BDA000020982568001310
Figure BDA000020982568001311
其中Pij为投影矩阵P在第i行第j列上的投影系数；
根据
Figure BDA000020982568001312
等于投影数据
Figure BDA000020982568001313
的第i行，
Figure BDA000020982568001314
等于投影数据
Figure BDA000020982568001315
的第i行，计算获得与
S302.2D、根据
Figure BDA000020982568001318
利用二分法计算出
Figure BDA000020982568001319
完成像素点
Figure BDA000020982568001320
的更新；
在这里，计算出
Figure BDA000020982568001321
完成像素点
Figure BDA000020982568001322
的更新，即得到更新的图像向量
Figure BDA000020982568001323
S302.2E、更新正投影数据
需要说明的是，所述迭代步骤循环的条件为输出条件未满足，所述输出条件具体可以设置为
Figure BDA00002098256800141
且
Figure BDA00002098256800142
其中ε为预设的小阈值，在理论上该预设的小阈值ε越小，图像质量越好，在本发明中，经过多次试验，该预设的小阈值ε至少应小于10-3，迭代优化最终得到的与
Figure BDA00002098256800144
图像质量较好。
本发明还包括一种实现多能量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组方程式
Figure BDA00002098256800145
计算出m个采样点的电压值分别对应的组合系数与
Figure BDA00002098256800147
将所述m个采样点的电压值分别对应的组合系数与
Figure BDA00002098256800149
发送至多能量图像重建单元404；
多能量图像重建单元404：用于通过计算所述第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800151
之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二，其中，
Figure BDA00002098256800153
和
Figure BDA00002098256800154
为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数
Figure BDA00002098256800155
与
Figure BDA00002098256800156
获得；为投影矩阵的第i行，
Figure BDA00002098256800158
为投影矩阵的第i行；将满足所述第一类扫描数据{yi}与组合投影数据
Figure BDA00002098256800159
之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个......”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均包含在本发明的保护范围内。

Claims (12)

1.一种实现多能量CT数据采集的方法，其特征在于，包括：
利用CT扫描仪控制电压在数据采集期间以周期性变化的方式对扫描对象进行扫描，所述电压在数据采集期间内为完整周期的变化。
2.根据权利要求1所述的方法，其特征在于，所述控制电压在数据采集期间以周期性变化的方式对扫描对象进行扫描，具体为控制电压在最高电压与最低电压之间的m个不同采样点周期性地对扫描对象进行扫描，其中m≥3。
3.根据权利要求1所述的方法，其特征在于，利用第一函数f(x)控制电压在数据采集期间以周期性变化的方式扫描，其中所述第一函数f(x)以2T为周期，并且在每个周期内，f(x)能取到最高电压和最低电压。
4.根据权利要求2所述的方法，其特征在于，所述第一函数f(x)具体为：
Figure FDA00002098256700011
其中，x为整数，2T为电压变化的周期，
Figure FDA00002098256700012
Figure FDA00002098256700013
MaxKV为最高电压，MinKV为最低电压。
5.一种实现多能量CT扫描重建的方法，其特征在于，包括第一采集步骤、第二采集步骤、多能量校正步骤和多能量图像重建步骤：
所述第一采集步骤包括：
以由第一物质和第二物质组成的第一类扫描对象作为扫描对象，应用如权利要求1-4所述的一种实现多能量CT数据采集的方法，采用第一类电压进行数据采集，采集到n个投影角度下的n组第一类扫描数据{yi}，其中，i为整数，1≤i≤n，n≥2，所述第一类电压在数据采集期间周期性变化，所述第一类电压在最高电压与最低电压之间共有m个不同的采样点，所述采用第一类电压进行数据采集具体为采用所述m个采样点的电压值在数据采集期间周期性地进行数据采集；
所述第二采集步骤包括：
利用CT扫描仪，控制所述m个采样点中每个采样点的电压值，分别对双柱校正模进行所述n个投影角度的扫描，获得m个n组第二类扫描数据，其中所述双柱校正模中包含所述第一类扫描对象的两种待分割物质，一柱体包含第一物质，另一柱体包含第二物质；
所述多能量校正步骤包括：
利用所述m个n组第二类扫描数据重建并校正出m个第二类扫描图像，其中，所述每个n组第二类扫描数据重建并校正出对应的一个第二类扫描图像；
根据所述第一物质与第二物质在双柱校正模中的位置范围，得到所述第一物质的第三类图像向量和第二物质的第三类图像向量；
根据m组方程式
Figure FDA00002098256700021
，计算出m个采样点的电压值分别对应的组合系数与
所述多能量图像重建步骤包括：
通过计算所述第一类扫描数据{yi}与组合投影数据之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据
Figure FDA00002098256700025
之间差距为最小值的图像向量X第一和X第二；
其中，
Figure FDA00002098256700026
和
Figure FDA00002098256700027
为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数
Figure FDA00002098256700028
与
Figure FDA00002098256700029
获得；
为投影矩阵的第i行；
将满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
6.根据权利要求5所述的方法，其特征在于，所述第一类电压在最高电压与最低电压之间m个不同采样点的个数m≥3。
7.根据权利要求5或6所述的方法，其特征在于，所述计算第一类扫描数据{yi}与组合投影数据
Figure FDA000020982567000212
之间差距的最小值，具体通过计算目标函数
Figure FDA00002098256700031
的最小值实现，其中，β为折衷参数。
8.根据权利要求7所述的方法，其特征在于，所述计算目标函数的最小值具体包括以下步骤：
以图像向量
Figure FDA00002098256700033
表示X第一；
以图像向量表示X第二；
定义r＝0，
Figure FDA00002098256700035
和
Figure FDA00002098256700036
为初始图像；
利用投影矩阵P对
Figure FDA00002098256700037
进行正投影，获得初始的正投影数据
Figure FDA00002098256700038
利用投影矩阵P对
Figure FDA00002098256700039
进行正投影，获得初始的正投影数据
Figure FDA000020982567000310
若输出条件未满足，则重复以下迭代步骤，其中，本次迭代为第r+1次迭代，表示该第r+1次迭代得到的图像向量X第一，
Figure FDA000020982567000312
表示该第r+1次迭代前的图像向量X第一，
Figure FDA000020982567000313
表示该第r+1次迭代得到的图像向量X第二，
Figure FDA000020982567000314
表示该第r+1次迭代前的图像向量X第二，迭代步骤每执行完一次将r的值加1，所述迭代步骤包括图像向量
Figure FDA000020982567000315
的更新和图像向量
Figure FDA000020982567000316
的更新：
所述
Figure FDA000020982567000317
的更新步骤包括：
重复以下步骤，直到
Figure FDA00002098256700041
所有的像素点都被更新：
选择图像向量
Figure FDA00002098256700042
中的任意一像素点
Figure FDA00002098256700043
且该像素点在该第r+1次迭代步骤中未被更新；
利用函数
Figure FDA00002098256700044
将
Figure FDA00002098256700045
代入后，得到ψ关于
Figure FDA00002098256700046
的偏导等于0的方程：
Figure FDA00002098256700047
定义：
Figure FDA00002098256700049
其中Pij为投影矩阵P在第i行第j列上的投影系数；
根据等于投影数据
Figure FDA000020982567000411
的第i行，等于投影数据
Figure FDA000020982567000413
的第i行，计算获得与
Figure FDA000020982567000415
根据
Figure FDA000020982567000416
利用二分法计算出
Figure FDA000020982567000417
完成像素点
Figure FDA000020982567000418
的更新；
更新正投影数据
Figure FDA000020982567000419
所述
Figure FDA000020982567000420
的更新步骤包括：
重复以下步骤，直到所有的像素点都被更新：
选择图像向量
Figure FDA000020982567000422
中的任意一像素点
Figure FDA000020982567000423
且该像素点在该第r+1次迭代步骤中未被更新；
利用函数
Figure FDA00002098256700051
将
Figure FDA00002098256700052
代入后，得到ψ关于
Figure FDA00002098256700053
的偏导等于0的方程：
Figure FDA00002098256700054
定义：
Figure FDA00002098256700055
根据
Figure FDA00002098256700056
等于投影数据
Figure FDA00002098256700057
的第i行，
Figure FDA00002098256700058
等于投影数据
Figure FDA00002098256700059
的第i行，计算获得与
Figure FDA000020982567000511
根据
Figure FDA000020982567000512
利用二分法计算出
Figure FDA000020982567000513
完成像素点
Figure FDA000020982567000514
的更新
更新正投影数据
Figure FDA000020982567000515
9.根据权利要求8所述的方法，其特征在于，所述输出条件具体为
Figure FDA000020982567000516
其中ε为预设的小阈值。
10.根据权利要求9所述的方法，其特征在于，所述预设的小阈值ε至少小于10-3。
11.根据权利要求5所述的方法，其特征在于，所述根据m组方程式
Figure FDA000020982567000517
计算出m个采样点的电压值分别对应的组合系数
Figure FDA000020982567000518
与具体采用最小二分法进行计算。
12.一种实现多能量CT扫描重建的装置，其特征在于，包括：
第一采集单元：用于以由第一物质和第二物质组成的第一类扫描对象作为扫描对象，应用如权利要求1-4所述的一种实现多能量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组方程式
Figure FDA00002098256700061
计算出m个采样点的电压值分别对应的组合系数
Figure FDA00002098256700062
与
Figure FDA00002098256700063
将所述m个采样点的电压值分别对应的组合系数
Figure FDA00002098256700064
与发送至多能量图像重建单元；
多能量图像重建单元：用于通过计算所述第一类扫描数据{yi}与组合投影数据
Figure FDA00002098256700066
之间差距的最小值，得到满足所述第一类扫描数据{yi}与组合投影数据之间差距为最小值的图像向量X第一和X第二，其中，和
Figure FDA00002098256700069
为所述第一类电压在采集第i个投影角度的第一类扫描数据时对应的组合系数，具体根据所述m个采样点的电压值分别对应的组合系数
Figure FDA000020982567000610
与
Figure FDA000020982567000611
获得；
Figure FDA000020982567000612
为投影矩阵的第i行；将满足所述第一类扫描数据{yi}与组合投影数据
Figure FDA000020982567000613
之间差距为最小值的图像向量X第一和X第二分别作为待求的第一物质的第一类图像向量的转置和第二物质的第一类图像向量的转置。
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 true CN103654818A (zh) 2014-03-26
CN103654818B 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)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105212958A (zh) * 2015-10-21 2016-01-06 沈阳东软医疗系统有限公司 Ct图像的获得方法、装置及设备
CN105758873A (zh) * 2015-03-04 2016-07-13 公安部第研究所 一种ct探测装置及其数据处理方法
CN106651753A (zh) * 2016-09-28 2017-05-10 沈阳东软医疗系统有限公司 提高ct图像显示效果的方法及装置
CN106780654A (zh) * 2017-01-24 2017-05-31 沈阳东软医疗系统有限公司 一种图像重建方法和装置
CN112697821A (zh) * 2020-12-02 2021-04-23 赛诺威盛科技（北京）有限公司 多能谱ct扫描方法、装置、电子设备和ct设备

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107884806B (zh) * 2017-10-19 2020-05-08 天津大学 面向双能ct成像的x射线能谱探测及重构解析方法

Citations (4)

* 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设备
JP2009000293A (ja) * 2007-06-21 2009-01-08 Toshiba Corp デュアルエネルギシステム及びその画像収集方法
CN101772324A (zh) * 2007-08-06 2010-07-07 株式会社日立医药 X射线ct装置
CN103202707A (zh) * 2012-01-16 2013-07-17 株式会社东芝 X射线ct装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008072175A1 (en) * 2006-12-15 2008-06-19 Koninklijke Philips Electronics N.V. Spectrally resolving x-ray imaging device
WO2009102996A2 (en) * 2008-02-15 2009-08-20 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 (4)

* 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设备
JP2009000293A (ja) * 2007-06-21 2009-01-08 Toshiba Corp デュアルエネルギシステム及びその画像収集方法
CN101772324A (zh) * 2007-08-06 2010-07-07 株式会社日立医药 X射线ct装置
CN103202707A (zh) * 2012-01-16 2013-07-17 株式会社东芝 X射线ct装置

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105758873A (zh) * 2015-03-04 2016-07-13 公安部第研究所 一种ct探测装置及其数据处理方法
CN105758873B (zh) * 2015-03-04 2019-12-31 公安部第一研究所 一种ct探测装置及其数据处理方法
CN105212958A (zh) * 2015-10-21 2016-01-06 沈阳东软医疗系统有限公司 Ct图像的获得方法、装置及设备
CN106651753A (zh) * 2016-09-28 2017-05-10 沈阳东软医疗系统有限公司 提高ct图像显示效果的方法及装置
CN106651753B (zh) * 2016-09-28 2020-03-17 东软医疗系统股份有限公司 提高ct图像显示效果的方法及装置
CN106780654A (zh) * 2017-01-24 2017-05-31 沈阳东软医疗系统有限公司 一种图像重建方法和装置
CN112697821A (zh) * 2020-12-02 2021-04-23 赛诺威盛科技（北京）有限公司 多能谱ct扫描方法、装置、电子设备和ct设备
CN112697821B (zh) * 2020-12-02 2022-12-02 赛诺威盛科技(北京)股份有限公司 多能谱ct扫描方法、装置、电子设备和ct设备

Also Published As

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

Similar Documents

Publication Publication Date Title
CN103654818B (zh) 一种实现多能量ct扫描重建的方法及装置
CN104318536B (zh) Ct图像的校正方法及装置
CN103150712B (zh) 一种基于投影序列数据相似性的图像去噪方法
CN103186883B (zh) 一种ct图像重建中骨硬化伪影的校正方法
CN105225208B (zh) 一种计算机断层成像金属伪影校正方法及装置
CN103413280B (zh) 一种低剂量x射线ct图像重建方法
CN104156917A (zh) 基于双能谱的x射线ct图像增强方法
CN111265231B (zh) 分布式光源ct图像重建方法与系统
CN110742635B (zh) 一种复合能谱ct成像方法
CN109801343A (zh) 基于重建前后图像的环形伪影校正方法、ct控制系统
CN103308537B (zh) 一种递变能量x射线成像图像融合方法
Wang et al. Design, analysis and simulation for development of the first clinical micro-CT scanner1
CN104751502A (zh) 一种用于扩大视野的ct图像重建方法
US6351548B1 (en) Fast hierarchical reprojection algorithm for tomography
CN105488826A (zh) 一种基于fbp的能谱ct迭代成像方法和成像系统
CN105796121A (zh) 一种ct和x射线激发荧光双模同步断层成像方法
CN103793890A (zh) 一种能谱ct图像的恢复处理方法
CN110428478A (zh) 交替光源扇束x射线ct采样方法及装置
JP5570716B2 (ja) Ｘ線ｃｔ装置およびプログラム
CN110458908A (zh) 基于有限角度迭代重建超视野ct图像的方法
CN109816747A (zh) 一种头部ct图像的金属伪影校正方法
CN104361615B (zh) 一种采用圆轨道扇形束x射线ct扫描机快速重建断层图像的方法
CN105222730B (zh) 一种基于图像复原的工业ct几何尺寸测量方法
CN109448071A (zh) 一种能谱图像重建方法及系统
CN104545964A (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

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

CP03 Change of name, title or address