CN111839568B - 一种新型大视场直线扫描ct***及图像重建方法 - Google Patents
一种新型大视场直线扫描ct***及图像重建方法 Download PDFInfo
- Publication number
- CN111839568B CN111839568B CN202010713021.3A CN202010713021A CN111839568B CN 111839568 B CN111839568 B CN 111839568B CN 202010713021 A CN202010713021 A CN 202010713021A CN 111839568 B CN111839568 B CN 111839568B
- Authority
- CN
- China
- Prior art keywords
- scanning
- detector
- ray source
- field
- view
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000001514 detection method Methods 0.000 claims abstract description 25
- 238000012545 processing Methods 0.000 claims abstract description 8
- 102000011990 Sirtuin Human genes 0.000 claims description 9
- 108050002485 Sirtuin Proteins 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 4
- 230000002238 attenuated effect Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 2
- 230000001105 regulatory effect Effects 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 abstract description 4
- 238000002591 computed tomography Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000007689 inspection Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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/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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4476—Constructional features of apparatus for radiation diagnosis related to motor-assisted motion of the source unit
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Optics & Photonics (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种新型大视场直线扫描CT***及图像重建方法，属于X射线成像检测技术领域。该***包括X射线源、探测器、检测对象、数据采集***、图像处理***、控制***和机械***；X射线源为锥束射线源，探测器为面阵平板探测器。扫描过程中，探测器固定，X射线源相对探测器做平行直线运动。本发明与传统CT***相比，增大了成像视场，具有结构简单，成本低，易于实现等优点。针对本发明大视场直线扫描CT***图像重建，本发明构建一种基于TV最小化的优化重建算法，可以直接对该***采集的投影数据进行重建，具有较好的抗噪声和去伪影能力。
Description
技术领域
本发明属于X射线成像检测技术领域，涉及一种新型大视场直线扫描CT***及图像重建方法。
背景技术
计算机断层成像技术(computed tomography，CT)自发明以来，广泛用于医疗、工业、安全检查等领域。在医疗领域，受CT机器结构限制，经常对一些肥胖病人进行CT扫描无能为力。在工业以及安全检查领域，对一些超出传统锥束CT视场的大物体的断层成像检测，如集装箱、发动机、大型涡轮叶片等，还存在亟待解决的问题。另一方面，在CT检测中，要求检测对象的大小必须在扫描视野范围内，如果检测对象过大，需要适当降低放大比，牺牲图像的分辨力来完成扫描，如何在保证高图像分辨力的同时，让扫描包含整个物体的支撑仍是一个未解决的问题。
最近，公开号为CN104809750A的专利提出了一种结构简单、低成本、可移动/便携的直线扫描CT***，***采用X射线源和探测器沿不同的方向平行移动的扫描方式，并且已经被证实在低成本CT扫描仪有巨大的潜力。为了实现对超出直线扫描CT视场的物体进行成像，公开号为CN106447740的专利提出了一种相对平行直线扫描CT感兴趣区域图像重建方法，该方法可对超出直线扫描CT视场的物体进行局部成像，但此方法不能对整个物体进行成像。为了满足对大物体的成像需求，本发明旨在提供一种新型直线扫描方法，可以增大扫描视场，同时继承直线扫描CT***结构简单、低成本、可移动/便携的优点。
发明内容
有鉴于此，本发明的目的之一在于提供一种新型大视场直线扫描CT***，该***采用探测器固定，射线源相对探测器做平行直线运动的扫描方式，从而扩大射线束扫描范围，以达到增大扫描视场的目的。
本发明的目的之二在于提供一种新型大视场直线扫描CT的图像重建方法，利用压缩感知在不完备数据CT图像重建中发挥的作用，构建针对新型大视场直线扫描CT***的TV最小化图像重建算法，为新型大视场直线扫描CT***提供相应的图像重建方法。
为达到上述目的之一，本发明提供如下技术方案：
一种新型大视场直线扫描CT***，包括X射线源(1)、检测对象(2)、探测器(3)、数据采集***(4)、图像处理***(5)、控制***(6)和机械***(7)；控制***(6)控制X射线源(1)发出X射线，穿透检测对象，到达探测器(3)；探测器(3)将探测到的衰减后的X射线转换成电信号，并由数据采集***(4)负责采集；采集后的信号传送给图像处理***(5)重建后显示；机械***(7)负责X射线源(1)的扫描运动、和检测物体(2)以及探测器(3)的位置调整运动。
在扫描过程中，探测器(3)固定，X射线源(1)相对于探测器做平行直线运动。
进一步，所述X射线源(1)为锥束射线源，探测器(3)为面阵平板探测器。
进一步，设探测器(3)长度为LD，其中DOD、SOD分别为探测器(3)和X射线源(1)到扫描中心点O的距离，则该***的扫描视场半径为
进一步，可得X射线源(1)的扫描行程为LS∈[LSmin,LSmax]，其中，
进一步，所述X射线源(1)的扫描方式采用等距扫描，其扫描间距ΔLS＝LS/(N-1)，其中N为一次直线扫描采样点总数。
进一步，该***采用T次圆周均布的直线扫描实现对检测对象的扫描，其中，
T＝[π/θ]
[π/θ]表示大于π/θ的最小整数。
进一步，以扫描物体为中心建立坐标系，平行于X射线源(1)运动轨迹的轴为x轴，垂直于X射线源(1)运动轨迹的轴为y轴，记点S1(-LSmin/2,-SOD)，对于任意扫描行程LS∈[LSmin,LSmax]，有点S2(-LS/2,-SOD)，记过点S1与视场边缘相切的直线为L1，过点S2与视场边缘相切的直线为L2，P(xP,yP)为直线L1与L2的交点，则θ为直线L1与L2形成的夹角，有
进一步，可得直线L1过点S1(-LSmin,-SOD)，以及点(-LD/2,DOD)，其方程为：
y＝k1(x+LD/2)+DOD,
其中，对于任意扫描行程LS∈[LSmin,LSmax]，可得过点S2(-LS/2,-SOD)与视场边缘相切的直线L2为：
y＝k2(x+LS/2)-SOD,
其中，则直线L1与L2的交点P(xP,yP)为：
进一步，该***对检测对象进行一次圆周扫描获得的投影数为I＝N×T。
为达到上述目的之二，本发明提供如下技术方案：
一种新型大视场直线扫描CT的图像重建方法，利用迭代方式进行图像重建，根据压缩感知构建基于全变分(total variation，TV)最小化的优化图像重建算法，优化图像质量。
进一步，构建的基于TV最小化的优化图像重建算法，表达式为：
其中，为***矩阵，/>为重建图像，/>为矫正后的投影数据，λ表示正则化参数，用于权衡数据保证项和正则项，TV(f)＝||▽f||1表示图像梯度L1范数。
本发明的有益效果在于：
(1)本发明提供的大视场直线扫描CT***相对于相对直线扫描CT相比，本发明可以扩大射线束扫描范围，从而增大扫描视场。
(2)本发明的另一个特点在于，扫描方式简单，成本低，运动简单易于实现。
(3)针对本发明的新型大视场直线扫描CT***，本发明引入压缩感知，利用图像在梯度变换域的稀疏性，构建基于TV最小化的重建方法，对投影数据进行直接图像重建，可以减少数据不一致性对重建图像的影响。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作优选的详细描述，其中：
图1为本发明所述的大视场直线扫描CT***模型；
图2为本发明***的一种结构示意图；
图3为本发明***的一种扫描运动示意图；
图4为本发明***的CT扫描二维几何模型；
图5为等距3次圆周均布直线扫描方式示意图；
图6为采用本发明***与传统CT***重建的效果对比图，其中(a)为实验采用的Shepp-Logan模型；(b)为采用本发明所述的图像重建方法针对新型大视场直线扫描CT重建的图像；(c)为相同扫描参数下，采用传统CT扫描，SIRT算法重建的图像；
其中，1为X射线源，2为检测对象，3为探测器，4为数据采集***，5为图像处理***，6为控制***，7为机械***，8为旋转转台，9为直线滑轨。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。需要说明的是，以下实施例中所提供的图示仅以示意方式说明本发明的基本构想，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。
请参阅图1～图6，图1为一种新型大视场直线扫描CT***的模型图，该***包括X射线源1、检测对象2、探测器3，以及数据采集***4、图像处理***5、控制***6和机械***7。控制***6控制X射线源1发出X射线，穿透检测对象2，到达探测器3；探测器3将探测到的衰减后的X射线转换成电信号，并由数据采集***4负责采集；采集后的信号交由图像处理***5重建后显示。在扫描过程中，探测器3固定，X射线源1相对于探测器3做平行直线运动。图2提供了一种参考结构示意，由X射线源1、旋转转台8、探测器3和直线滑轨9组成，检测对象可放置在旋转转台8上，直线滑轨9分别实现射线源1的扫描运动，探测器3与检测对象2、射线源1的距离调整。图3为对应的扫描运动示意图，在一次直线扫描运动过程中，检测对象2与探测器3固定，X射线源1做直线运动，其扫描轨迹平行于探测器3中心行。一次直线扫描完成后，检测对象2可旋转相应角度，实现第二次直线扫描。
参阅图4，对于确定的探测器3长度为LD，以及探测器3和X射线源1到扫描中心点O的距离，可以得到新型大视场直线扫描CT的扫描半径为：
其中，DOD，SOD分别为探测器(3)和X射线源(1)到扫描中心点O的距离。在相同扫描参数下，传统圆周CT扫描***以及相对平行直线扫描CT***的扫描视场为：
参阅图4可以直观看出，本发明的新型大视场直线扫描CT较传统CT视场明显增大。
对于确定的探测器(3)大小，视场半径，可确定X射线源(1)的扫描行程为LS∈[LSmin,LSmax]，其中,
对于一次直线扫描采样点数N，采用等距扫描方式，可得扫描间距为ΔLS＝LS/(N-1)。
根据Quinto对CT图像伪影分析的理论，在扫描视场圆周内，与射线相切的点能被精确重建。显然，一次直线运动扫描不能满足扫描视场圆周上所有的点都与射线相切，同时由于本发明扫描方式具有对称性，为此，本发明采用T次圆周均布的直线扫描实现对检测对象的扫描。其中：
T＝[π/θ] (5)
其中，以扫描物体为中心建立坐标系，平行于X射线源(1)运动轨迹的轴为x轴，垂直于X射线源(1)运动轨迹的轴为y轴，记点S1(-LSmin/2,-SOD)，对于任意扫描行程LS∈[LSmin,LSmax]，有点S2(-LS/2,-SOD)，记过点S1与视场边缘相切的直线为L1，过点S2与视场边缘相切的直线为L2，P(xP,yP)为直线L1与L2的交点，则θ表示直线L1与L2形成的夹角，有
由于直线L1过点S1(-LSmin,-SOD)，以及(-LD/2,DOD)，可得其方程为：
y＝k1(x+LD/2)+DOD, (7)
其中，对于任意扫描行程LS∈[LSmin,LSmax]，可得过点S2(-LS/2,-SOD)与视场边缘相切的直线L2为：
y＝k2(x+LS/2)-SOD, (8)
其中，进一步可得直线L1与L2的交点P(xP,yP)为：
由式(5)(6)(9)可知，当探测器大小，物体中心到射线源以及探测器距离确定，扫描次数T的大小与一次直线扫描射线源运动行程LS有关，对于LS∈[LSmin,LSmax]，LS越大，扫描次数T越小。
针对本发明的一种新型大视场直线扫描CT***，提供了相应的图像重建算法。超出传统CT视野物体的重建问题可以归纳为一种不完全数据重建问题。本发明引入压缩感知理论，利用图像在梯度变换域的稀疏性，构建基于TV最小化的优化图像重建算法，其表达式如下：
其中，为***矩阵，/>为重建图像，/>为矫正后的投影数据，λ表示正则化参数，用于权衡数据保真项和正则项，/>表示图像梯度L1范数。式(10)的求解可分为如下两步：
步骤1：凸集投影
(1)取初值f(k＝0)＝0，其中，k为迭代次数；
(2)利用SIRT算法实现数据保真；
其中，表示重建图像第j个元素第k次迭代的SIRT结果，Ai为***矩阵的第i行行向量，β为松弛因子；
(3)非负约束：
其中，表示重建图像第j个元素第k+1次非负约束结果；
步骤2：TV最小化
(1)fTV(m＝0)＝fPOCS(k+1)，其中m表示TV最小化的迭代次数；
(2)计算TV的梯度及梯度方向：
(3)沿TV梯度下降的方向迭代修正图像：
其中，α为调节因子；
(4)若m＞M，则
fSIRT(k)＝fTV(m+1) (15)
其中，M为TV最小化的最大迭代次数；
(5)若k＞K，则输出迭代结果。其中K为最大迭代次数。
验证实施例：
利用MATLAB工具，仿真本发明的大视场直线扫描CT及图像重建方法，其中直线扫描次数为3次，如图5所示，其仿真扫描参数如表1所示。
表1.仿真扫描参数
其重建图像大小为256×256，迭代次数SIRT迭代次数K＝500，TV最小化迭代次数M＝20。图6(b)为采用本发明所述的图像重建方法针对新型大视场直线扫描CT***重建的图像；图6(c)为相同扫描参数下，利用SIRT算法对传统CT扫描数据的重建结果。对比图6(a)采用的Shepp-Logan模型可以看出，本发明扫描CT***能有效增大扫描视场，并精确恢复图像结构信息。
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种新型大视场直线扫描CT***，其特征在于，该***包括X射线源(1)、探测器(3)、数据采集***(4)、图像处理***(5)、控制***(6)和机械***(7)；所述控制***(6)控制X射线源(1)发出X射线，穿透检测对象(2)，到达探测器(3)；所述探测器(3)将探测到的衰减后的X射线转换成电信号，并由数据采集***(4)负责采集；采集后的信号传送给图像处理***(5)重建后显示；机械***(7)负责X射线源(1)的扫描运动和检测物体(2)以及探测器(3)的位置调整运动；在扫描过程中，探测器(3)固定，X射线源(1)相对于探测器(3)做平行直线运动；
该***的扫描视场半径为：
其中，DOD、SOD分别为探测器(3)和X射线源(1)到扫描中心点O的距离，LD为探测器(3)长度；
X射线源(1)的扫描行程为LS∈[LSmin,LSmax]，其中，
所述X射线源(1)的扫描方式采用等距扫描，其扫描间距ΔLS＝LS/(N-1)，其中N为一次直线扫描采样点总数；
以检测对象(2)为中心建立坐标系，平行于X射线源(1)运动轨迹的轴为x轴，垂直于X射线源(1)运动轨迹的轴为y轴，记点S1(-LSmin/2,-SOD)，对于任意扫描行程LS∈[LSmin,LSmax]，有点S2(-LS/2,-SOD)，记过点S1与视场边缘相切的直线为L1，过点S2与视场边缘相切的直线为L2，P(xP,yP)为直线L1与L2的交点，则θ为直线L1与L2形成的夹角，有
该***采用T次圆周均布的直线扫描方式对检测进行扫描，则检测对象进行一次圆周扫描获得的投影数为I＝N×T，其中，T＝[π/θ]，[π/θ]为大于π/θ的最小整数，N为一次直线扫描采样点总数；
过点S1(-LSmin,-SOD)以及点(-LD/2,DOD)的直线L1为：
y＝k1(x+LD/2)+DOD
其中，对于任意扫描行程LS∈[LSmin,LSmax]，过点S2(-LS/2,-SOD)与视场边缘相切的直线L2为：
y＝k2(x+LS/2)-SOD
其中，则直线L1与L2的交点P(xP,yP)为：
2.根据权利要求1所述的一种新型大视场直线扫描CT***，其特征在于，所述X射线源(1)为锥束射线源，所述探测器(3)为面阵平板探测器。
3.根据权利要求1或2所述的新型大视场直线扫描CT***，其特征在于，针对该***的图像重建算法为：利用迭代方式进行图像重建，根据压缩感知构建基于TV最小化的优化图像重建算法，优化图像质量；所述TV是全变分；
构建的基于TV最小化的优化图像重建算法，表达式为：
其中，为***矩阵，/>为重建图像，/>为矫正后的投影数据，λ表示正则化参数，TV(f)＝||▽f||1表示图像梯度L1范数；
上式求解分为如下两步：
步骤1：凸集投影
(1)取初值f(k＝0)＝0，其中，k为迭代次数；
(2)利用SIRT算法实现数据保真；
其中，表示重建图像第j个元素第k次迭代的SIRT结果，Ai为***矩阵的第i行行向量，β为松弛因子；
(3)非负约束：
其中，表示重建图像第j个元素第k+1次非负约束结果；
步骤2：TV最小化
(1)fTV(m＝0)＝fPOCS(k+1)，其中m表示TV最小化的迭代次数；
(2)计算TV的梯度及梯度方向：
(3)沿TV梯度下降的方向迭代修正图像：
其中，α为调节因子；
(4)若m＞M，则
fSIRT(k)＝fTV(m+1)
其中，M为TV最小化的最大迭代次数；
(5)若k＞K，则输出迭代结果；其中K为最大迭代次数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010713021.3A CN111839568B (zh) | 2020-07-22 | 2020-07-22 | 一种新型大视场直线扫描ct***及图像重建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010713021.3A CN111839568B (zh) | 2020-07-22 | 2020-07-22 | 一种新型大视场直线扫描ct***及图像重建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111839568A CN111839568A (zh) | 2020-10-30 |
| CN111839568B true CN111839568B (zh) | 2023-12-12 |
Family
ID=72950255
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010713021.3A Active CN111839568B (zh) | 2020-07-22 | 2020-07-22 | 一种新型大视场直线扫描ct***及图像重建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111839568B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113533392B (zh) * | 2021-07-12 | 2022-08-26 | 重庆大学 | 一种组合扫描cl成像方法 |
| CN114758079B (zh) * | 2022-05-19 | 2024-02-20 | 国网河南省电力公司电力科学研究院 | 电气设备内部结构局部成像三维重建方法 |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008122971A1 (en) * | 2007-04-10 | 2008-10-16 | Arineta Ltd. | Cone-beam ct |
| CN101393145A (zh) * | 2008-10-22 | 2009-03-25 | 重庆大学 | 大尺寸物体的锥束双螺旋ct扫描成像方法 |
| CN201242531Y (zh) * | 2008-04-17 | 2009-05-20 | 清华大学 | 一种直线轨迹扫描成像*** |
| CN101516268A (zh) * | 2006-09-29 | 2009-08-26 | 皇家飞利浦电子股份有限公司 | 飞掠扫描 |
| WO2010030270A1 (en) * | 2008-09-10 | 2010-03-18 | Analogic Corporation | Ct scanning systems and methods using multi-pixel x-ray sources |
| CN102004111A (zh) * | 2010-09-28 | 2011-04-06 | 北京航空航天大学 | 一种倾斜多锥束直线轨迹ct成像方法 |
| CN103163165A (zh) * | 2013-02-28 | 2013-06-19 | 重庆大学 | 一种二代ct扫描成像方法 |
| CN104597062A (zh) * | 2015-02-02 | 2015-05-06 | 天津三英精密仪器有限公司 | 一种柱形束大视场x射线ct成像*** |
| CN104757988A (zh) * | 2015-05-04 | 2015-07-08 | 重庆大学 | 一种电子直线扫描微纳焦点ct扫描***及方法 |
| CN104809750A (zh) * | 2015-05-04 | 2015-07-29 | 重庆大学 | 一种直线扫描ct***及图像重建方法 |
| CN107764846A (zh) * | 2017-10-20 | 2018-03-06 | 重庆大学 | 一种正交直线扫描的cl成像***及分析方法 |
| CN107796834A (zh) * | 2017-10-20 | 2018-03-13 | 重庆大学 | 一种正交电子直线扫描cl成像***及方法 |
| CN109106390A (zh) * | 2018-06-28 | 2019-01-01 | 北京航空航天大学 | 一种旋转中心偏置的扇束短扫描ct重建方法及装置 |
| CN110179486A (zh) * | 2019-04-11 | 2019-08-30 | 清华大学 | 多能量ct成像***及其应用 |
| CN111077173A (zh) * | 2020-01-03 | 2020-04-28 | 重庆大学 | 一种倒置结构大视场Mirco-CT扫描成像*** |
| CN111265231A (zh) * | 2019-04-15 | 2020-06-12 | 清华大学 | 分布式光源ct图像重建方法与*** |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8571172B2 (en) * | 2008-04-14 | 2013-10-29 | Arineta Ltd. | CT cone beam scanner |
-
2020
- 2020-07-22 CN CN202010713021.3A patent/CN111839568B/zh active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101516268A (zh) * | 2006-09-29 | 2009-08-26 | 皇家飞利浦电子股份有限公司 | 飞掠扫描 |
| WO2008122971A1 (en) * | 2007-04-10 | 2008-10-16 | Arineta Ltd. | Cone-beam ct |
| CN201242531Y (zh) * | 2008-04-17 | 2009-05-20 | 清华大学 | 一种直线轨迹扫描成像*** |
| WO2010030270A1 (en) * | 2008-09-10 | 2010-03-18 | Analogic Corporation | Ct scanning systems and methods using multi-pixel x-ray sources |
| CN101393145A (zh) * | 2008-10-22 | 2009-03-25 | 重庆大学 | 大尺寸物体的锥束双螺旋ct扫描成像方法 |
| CN102004111A (zh) * | 2010-09-28 | 2011-04-06 | 北京航空航天大学 | 一种倾斜多锥束直线轨迹ct成像方法 |
| CN103163165A (zh) * | 2013-02-28 | 2013-06-19 | 重庆大学 | 一种二代ct扫描成像方法 |
| CN104597062A (zh) * | 2015-02-02 | 2015-05-06 | 天津三英精密仪器有限公司 | 一种柱形束大视场x射线ct成像*** |
| CN104757988A (zh) * | 2015-05-04 | 2015-07-08 | 重庆大学 | 一种电子直线扫描微纳焦点ct扫描***及方法 |
| CN104809750A (zh) * | 2015-05-04 | 2015-07-29 | 重庆大学 | 一种直线扫描ct***及图像重建方法 |
| CN107764846A (zh) * | 2017-10-20 | 2018-03-06 | 重庆大学 | 一种正交直线扫描的cl成像***及分析方法 |
| CN107796834A (zh) * | 2017-10-20 | 2018-03-13 | 重庆大学 | 一种正交电子直线扫描cl成像***及方法 |
| CN109106390A (zh) * | 2018-06-28 | 2019-01-01 | 北京航空航天大学 | 一种旋转中心偏置的扇束短扫描ct重建方法及装置 |
| CN110179486A (zh) * | 2019-04-11 | 2019-08-30 | 清华大学 | 多能量ct成像***及其应用 |
| CN111265231A (zh) * | 2019-04-15 | 2020-06-12 | 清华大学 | 分布式光源ct图像重建方法与*** |
| CN111077173A (zh) * | 2020-01-03 | 2020-04-28 | 重庆大学 | 一种倒置结构大视场Mirco-CT扫描成像*** |
Non-Patent Citations (3)
| Title |
|---|
| Fenglin Liu.Top-level design and pilot analysis of low-end CT scanners based on linear scanning for developing countries.《Journal of X-Ray Science and Technology》.2014,第22卷(第5期),673-686. * |
| In-line NDT with X-Ray CT combining sample rotation and translation;Thomas De Schryver等;《NDT&E international》;第84卷;89-98 * |
| Top-level design and pilot analysis of low-end CT scanners based on linear scanning for developing countries;Fenglin Liu;《Journal of X-Ray Science and Technology》;第22卷(第5期);673-686 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111839568A (zh) | 2020-10-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110807737B (zh) | 迭代图像重建框架 | |
| US8737711B2 (en) | X-ray CT image forming method and X-ray CT apparatus using the same | |
| US6907102B1 (en) | Iterative reconstruction methods for multi-slice computed tomography | |
| US8971599B2 (en) | Tomographic iterative reconstruction | |
| EP0492896A2 (en) | Reconstructing a 3-D image | |
| US20070116175A1 (en) | System and method for reconstructing image by using straight-line trajectory scan | |
| CN111265231B (zh) | 分布式光源ct图像重建方法与*** | |
| CN111839568B (zh) | 一种新型大视场直线扫描ct***及图像重建方法 | |
| CN107822652B (zh) | 用于重建光谱结果图像数据的方法 | |
| US20040258194A1 (en) | Fourier space tomographic image reconstruction method | |
| CN111739113B (zh) | 直线分布式光源和探测器的ct图像重建方法及装置 | |
| Ametova et al. | Software-based compensation of instrument misalignments for X-ray computed tomography dimensional metrology | |
| Sun et al. | Reduction of artifacts induced by misaligned geometry in cone-beam CT | |
| JP7206073B2 (ja) | 医用画像処理装置、医用画像生成装置、医用画像処理方法、およびプログラム | |
| CN112070856B (zh) | 基于非下采样轮廓波变换的有限角c型臂ct图像重建方法 | |
| CN112070704A (zh) | 一种基于紧小波框架的双正则化有限角ct图像重建方法 | |
| CN111583389A (zh) | 一种基于cad模型的不完全扫描ct图像重建方法 | |
| Deng et al. | An improved total variation minimization method using prior images and split-Bregman method in CT reconstruction | |
| Deng et al. | Limited-angle CT reconstruction with generalized shrinkage operators as regularizers | |
| Lu et al. | TV-based image reconstruction of multiple objects in a fixed source-detector geometry | |
| US9508164B2 (en) | Fast iterative image reconstruction method for 3D computed tomography | |
| CN110264536B (zh) | 一种在平行束超分重建中计算高低分辨率投影关系的方法 | |
| US9558569B2 (en) | Method and system for substantially reducing cone beam artifacts based upon image domain differentiation in circular computer tomography (CT) | |
| WO2023243503A1 (ja) | インテリアｃｔの画像再構成方法、画像再構成装置、及び、プログラム | |
| CN112529980B (zh) | 一种基于极大极小化的多目标有限角ct图像重建方法 |
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 |