CN112884862A - 基于质心投影轨迹拟合的锥束ct温度漂移校正方法及*** - Google Patents
基于质心投影轨迹拟合的锥束ct温度漂移校正方法及*** Download PDFInfo
- Publication number
- CN112884862A CN112884862A CN202110290331.3A CN202110290331A CN112884862A CN 112884862 A CN112884862 A CN 112884862A CN 202110290331 A CN202110290331 A CN 202110290331A CN 112884862 A CN112884862 A CN 112884862A
- Authority
- CN
- China
- Prior art keywords
- projection
- curve
- temperature drift
- image
- centroid
- 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
Links
- 238000012937 correction Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000013519 translation Methods 0.000 claims description 7
- 238000003702 image correction Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 6
- 238000002591 computed tomography Methods 0.000 description 42
- 238000004088 simulation Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 235000007688 Lycopersicon esculentum Nutrition 0.000 description 1
- 240000003768 Solanum lycopersicum Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010603 microCT Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- 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/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明属于CT图像重建技术领域,特别涉及一种基于质心投影轨迹拟合的锥束CT温度漂移校正方法及***,用于CT投影重建,首先依据CT原始投影图像数据,得到被投影物体质心投影轨迹曲线;然后通过局部区间长度比例及区间步长,利用该比例下的被投影物体局部区间质心投影估计未发生温度漂移时被投影物体整体投影轨迹曲线,通过做差获取温度漂移量;利用温度漂移量将原始投影图像亚像素级平移来实现CT图像偏移校正。本发明能够利用最小二乘拟合质心运动轨迹实现投影温度漂移校正,保证CT图像重建效果的同时,能够节约图像重建成本和资源,具有较好的应用前景。
Description
技术领域
本发明属于CT图像重建技术领域,特别涉及一种基于质心投影轨迹拟合的锥束CT温度漂移校正方法及***。
背景技术
X射线计算机断层扫描技术(CT)迅速发展,在工业领域得到了广泛的应用。对于微米CT和纳米CT,由于其高精度、高分辨率的特点,对射线源和机械结构的稳定性要求很高。微纳米CT在工作过程中扫描时间长,外部环境温度在扫描期间发生变化,由此引起射线源的偏移和物体热胀冷缩,使得投影在水平和垂直方向发生漂移,导致重建物体包含严重的运动伪影,表现为边缘轮廓位置的模糊。为了观察物体的精细结构,人们提出了消除伪影的校正算法。现有的温度漂移校正算法可以分为两类:基于短时参考的迭代补偿法和基于体模的补偿法,这些方法在一定程度上可以消去伪影,但是短时参考补偿法需要对物体进行额外的扫描,浪费了时间和计算资源。基于体模的补偿法需要额外制作校正体模,增加实验成本。
发明内容
为此,本发明提供一种基于质心投影轨迹拟合的锥束CT温度漂移校正方法及***,能够利用最小二乘拟合质心运动轨迹实现投影温度漂移校正,保证CT图像重建效果的同时,能够节约图像重建成本和资源。
按照本发明所提供的设计方案,提供一种基于质心投影轨迹拟合的锥束CT温度漂移校正方法,用于CT投影重建,包含如下内容:
依据CT原始投影图像数据,得到被投影物体质心投影轨迹曲线;
设置局部区间长度比例及区间步长,利用该比例下的被投影物体局部区间质心投影估计未发生温度漂移时被投影物体整体投影轨迹曲线;
依据估计得到的整体投影轨迹曲线和计算得到的投影轨迹曲线两者的曲线偏差获取温度漂移量;
通过温度漂移量将原始投影图像亚像素级平移来实现CT图像偏移校正。
作为本发明基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,利用CT原始投影图像,通过计算水平方向和垂直方向被投影物体原始质心来获取其投影轨迹曲线。
作为本发明基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,针对水平方向的温度漂移,利用二维投影图并通过固定纵坐标来遍历水平坐标获取一维数组,通过遍历投影角度得到水平方向的质心正弦图;针对垂直方向的温度漂移,利用二维投影图并依据设定的投影角度,遍历纵坐标对水平坐标的投影值获取一维曲线,计算一维曲线质心,通过遍历所有角度来获取垂直方向的质心正弦图。
作为本发明基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,通过设定投影间隔,使物体被扫描一周,并设定拟合曲线周期,利用先验知识的拟合曲线来估计整体投影轨迹曲线,其中,先验知识为水平方向和垂直方向投影总数拟合比例下的局部质心投影。
作为本发明基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,通过最小二乘法进行拟合来估计未发生温度漂移时被投影物体整体投影轨迹曲线,其中,拟合曲线方程表示为:r'(θ)=a+bsin(ωθ)+ccos(ωθ),θ为投影角度,a,b,c为待估计常数,w为已知常数。
作为本发明基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,温度漂移拟合比例设定为投影总数的后14%。
作为本发明基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,将拟合曲线和估计曲线在极值处对齐,获取对齐后的拟合曲线;通过对齐后的拟合曲线和估计曲线做差来得到温度漂移量。
作为本发明基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,将原始投影图像作为待平移图像,依据设定的图像插值倍数进行插值;针对插值后的图像,依据温度漂移量进行图像平移,通过图像采样来获取原始尺寸图像。
作为本发明基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,图像插值倍数大小设定为10。
进一步地,本发明还提供一种基于质心投影轨迹拟合的锥束CT温度漂移校正***,用于CT投影重建,包含:数据处理模块、曲线估计模块、温漂计算模块和图像校正模块,其中,
数据处理模块,用于依据CT原始投影图像数据,得到被投影物体质心投影轨迹曲线;
曲线估计模块,用于通过设置局部区间长度比例及区间步长,利用该比例下的被投影物体局部区间质心投影估计未发生温度漂移时被投影物体整体投影轨迹曲线;
温漂计算模块,用于依据估计得到的整体投影轨迹曲线和计算得到的投影轨迹曲线两者的曲线偏差获取温度漂移量;
图像校正模块,用于通过温度漂移量将原始投影图像亚像素级平移来实现CT图像偏移校正。
本发明的有益效果:
本发明通过质心轨迹的估计曲线和实际曲线之间的差异计算水平和垂直方向的漂移,可以实现图像在亚像素级的漂移校正;并进一步为了提高抗噪声性能,对实际曲线进行拟合,并和无漂移估计曲线对齐后计算各角度下漂移偏差,通过插值-采样的亚像素级像素的平移,将漂移校正后的投影实施重建,即可获得去除温度漂移的重建物体,提升物体重建质量。通过实际数据和仿真数据的测试,本案方案能够较好地校正温度漂移,使重建切片边缘更加尖锐,实现较好的效果,具有较好的应用前景。
附图说明:
图1为实施例中基于质心投影轨迹拟合的锥束CT温度漂移校正方法流程示意;
图2为实施例中锥束CT投影行质心正弦拟合温漂校正仿真效果图;
图3为束CT投影行质心正弦拟合温漂校正番茄种子伪影对比效果图。
具体实施方式:
为使本发明的目的、技术方案和优点更加清楚、明白,下面结合附图和技术方案对本发明作进一步详细的说明。
本发明实施例,提供一种基于质心投影轨迹拟合的锥束CT温度漂移校正方法,用于CT投影重建,包含如下内容:依据CT原始投影图像数据,得到被投影物体质心投影轨迹曲线;设置局部区间长度比例及区间步长,利用该比例下的被投影物体局部区间质心投影估计未发生温度漂移时被投影物体整体投影轨迹曲线;依据估计得到的整体投影轨迹曲线和计算得到的投影轨迹曲线两者的曲线偏差获取温度漂移量;通过温度漂移量将原始投影图像亚像素级平移来实现CT图像偏移校正。
参见图1所示,首先根据原始投影图像计算被投影物体质心的投影轨迹。然后利用局部质心位置估计未发生温度漂移时整体的质心轨迹,最终根据曲线间的偏差计算温度漂移量。经次方法计算的温度漂移量可以实现图像在亚像素级的漂移校正,将图像平移亚像素距离。将漂移校正后的投影实施重建,即可获得去除温度漂移的重建物体,提升图像重建效率。
作为本发明实施例中基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,利用CT原始投影图像,通过计算水平方向和垂直方向被投影物体原始质心来获取其投影轨迹曲线。进一步地,针对水平方向的温度漂移,利用二维投影图并通过固定纵坐标来遍历水平坐标获取一维数组,通过遍历投影角度得到水平方向的质心正弦图;针对垂直方向的温度漂移,利用二维投影图并依据设定的投影角度,遍历纵坐标对水平坐标的投影值获取一维曲线,计算一维曲线质心,通过遍历所有角度来获取垂直方向的质心正弦图。
基于CT图像投影图像,可直接对原始的投影数据进行处理。对于水平方向的温漂,首先获得固定某一y方向的正弦图(通常取中间切片的y位置),正弦图获取的方法是,投影Iθ(u,v)上固定纵坐标v0,遍历水平坐标u得到一维数组遍历投影角度θ,得到正弦图通过下式计算:
正弦图的质心记做ry(θ),通过下式计算:
在垂直方向,二维投影图Iθ(u,v)针对某一投影角度,遍历纵坐标v对水平坐标u投影值求和得到一维曲线tθ(v),计算一维曲线tθ(v)的质心,遍历所有角度θ,得到垂直方向的质心正弦曲线rx(θ)。
作为本发明实施例中基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,通过设定投影间隔,使物体被扫描一周,并设定拟合曲线周期,利用先验知识的拟合曲线来估计整体投影轨迹曲线,其中,先验知识为水平方向和垂直方向投影总数拟合比例下的局部质心投影。
使用最小二乘法估计未发生温度漂移的质心轨迹曲线r'(θ),由于正弦图的质心反映实际物体空间某一点的投影情况,可在进行实验时,使用全角度采集,也就是说通过设置投影间隔,使物体被扫描一周,在拟合曲线时可以加入先验知识,即正弦图质心的周期,根据以上分析,拟合曲线的周期应该为180,加入先验知识的拟合曲线,拟合误差更小。用于拟合的局部区间确定方法为:设置局部区间长度为投影数量的14%。以步长为1,将总投影划分为局部区间。将所有局部区间进行曲线拟合,其中包含二次谐波幅度最小的曲线对应的局部区间被认为是最优区间,使用此区间进行后续计算。可设拟合曲线方程为
r'(θ)=a+bsin(ωθ)+ccos(ωθ)
其中a,b,c,d为待估计常数,w为已知常数。
则A的计算式为
A=(UTU)-1UTR'
其中R'=[r'(θ1),r'(θ2)...r'(θn)],n为投影图像数。最终得到水平方向无温漂的估计曲线rx'(θ)和垂直方向估计曲线ry'(θ)。
作为本发明实施例中基于质心投影轨迹拟合的锥束CT温度漂移校正方法,进一步地,将拟合曲线和估计曲线在极值处对齐,获取对齐后的拟合曲线;通过对齐后的拟合曲线和估计曲线做差来得到温度漂移量。进一步地,将原始投影图像作为待平移图像,依据设定的图像插值倍数进行插值;针对插值后的图像,依据温度漂移量进行图像平移,通过图像采样来获取原始尺寸图像。
抖动误差可由原始测量曲线和拟合曲线的差值近似计算,由于微纳米CT的噪声较大,所以对原始质心轨迹rx(θ)和ry(θ)进行步骤3的拟合,得到原始质心轨迹的拟合曲线Rx(θ)和Ry(θ)。由于最小二乘拟合时周期存在误差,故需要分别将Rx(θ)和rx'(θ)、Ry(θ)和ry'(θ)在极值处对齐,Rx(θ±φ)和Ry(θ±φ)是分别与预测曲线对齐后的拟合曲线,分别记为Rx(θ)和Ry(θ)。将水平方向上的拟合曲线Rx(θ±φ)和估计曲线rx'(θ)做差,垂直方向上执行同样操作,得到漂移像素值,如下式所示:
Δ(θ)=R(θ)-r'(θ)
其中Δ(θ)为计算漂移值。
为实现图像的亚像素级平移,可将图像10倍插值,再计算漂移值扩大10倍平移,再对插值后的图像采样,得到原始图像尺寸。
进一步地,基于上述的方法,本发明实施例还提供一种基于质心投影轨迹拟合的锥束CT温度漂移校正***,用于CT投影重建,包含:数据处理模块、曲线估计模块、温漂计算模块和图像校正模块,其中,
数据处理模块,用于依据CT原始投影图像数据,得到被投影物体质心投影轨迹曲线;
曲线估计模块,用于通过设置局部区间长度比例及区间步长,利用该比例下的被投影物体局部区间质心投影估计未发生温度漂移时被投影物体整体投影轨迹曲线;
温漂计算模块,用于依据估计得到的整体投影轨迹曲线和计算得到的投影轨迹曲线两者的曲线偏差获取温度漂移量;
图像校正模块,用于通过温度漂移量将原始投影图像亚像素级平移来实现CT图像偏移校正。
为进一步验证本案方案有效性,参见图2和3所示,根据投影域局部质心位置通过最小二乘法拟合估计未产生漂移时质心轨迹曲线,通过计算预测曲线与实际曲线之间的差值,计算出投影在水平和垂直方向的温度漂移,通过图像插值和采样完成对图像的亚像素级平移。通过仿真数据和实际数据的测试完成所提出方法的验证,本案所提出的方案在保持重建切片细节的基础上,去除了切片边缘的漂移伪影,在视觉效果和定量指标均得到了满意的效果,能够快速准确的完成锥束CT温度漂移的校正。
除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对步骤、数字表达式和数值并不限制本发明的范围。
在这里示出和描述的所有示例中,任何具体值应被解释为仅仅是示例性的,而不是作为限制,因此,示例性实施例的其他示例可以具有不同的值。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种基于质心投影轨迹拟合的锥束CT温度漂移校正方法,用于CT投影重建,其特征在于,包含如下内容:
依据CT原始投影图像数据,得到被投影物体质心投影轨迹曲线;
设置局部区间长度比例及区间步长,利用该比例下的被投影物体局部区间质心投影估计未发生温度漂移时被投影物体整体投影轨迹曲线;
依据估计得到的整体投影轨迹曲线和计算得到的投影轨迹曲线两者的曲线偏差获取温度漂移量;
通过温度漂移量将原始投影图像亚像素级平移来实现CT图像偏移校正。
2.根据权利要求1所述的基于质心投影轨迹拟合的锥束CT温度漂移校正方法,其特征在于,利用CT原始投影图像,通过计算水平方向和垂直方向被投影物体原始质心来获取其投影轨迹曲线。
3.根据权利要求1或2所述的基于质心投影轨迹拟合的锥束CT温度漂移校正方法,其特征在于,针对水平方向的温度漂移,利用二维投影图并通过固定纵坐标来遍历水平坐标获取一维数组,通过遍历投影角度得到水平方向的质心正弦图;针对垂直方向的温度漂移,利用二维投影图并依据设定的投影角度,遍历纵坐标对水平坐标的投影值获取一维曲线,计算一维曲线质心,通过遍历所有角度来获取垂直方向的质心正弦图。
4.根据权利要求1所述的基于质心投影轨迹拟合的锥束CT温度漂移校正方法,其特征在于,通过设定投影间隔,使物体被扫描一周,并设定拟合曲线周期,利用先验知识的拟合曲线来估计整体投影轨迹曲线,其中,先验知识为水平方向和垂直方向投影总数拟合比例下的局部质心投影。
5.根据权利要求1或4所述的基于质心投影轨迹拟合的锥束CT温度漂移校正方法,其特征在于,通过最小二乘法进行拟合来估计未发生温度漂移时被投影物体整体投影轨迹曲线,其中,拟合曲线方程表示为:r'(θ)=a+b sin(ωθ)+c cos(ωθ),θ为投影角度,a,b,c为待估计常数,w为已知常数。
6.根据权利要求1或4所述的基于质心投影轨迹拟合的锥束CT温度漂移校正方法,其特征在于,温度漂移拟合比例设定为投影总数的后14%。
7.根据权利要求1所述的基于质心投影轨迹拟合的锥束CT温度漂移校正方法,其特征在于,将拟合曲线和估计曲线在极值处对齐,获取对齐后的拟合曲线;通过对齐后的拟合曲线和估计曲线做差来得到温度漂移量。
8.根据权利要求1所述的基于质心投影轨迹拟合的锥束CT温度漂移校正方法,其特征在于,将原始投影图像作为待平移图像,依据设定的图像插值倍数进行插值;针对插值后的图像,依据温度漂移量进行图像平移,通过图像采样来获取原始尺寸图像。
9.根据权利要求8所述的基于质心投影轨迹拟合的锥束CT温度漂移校正方法,其特征在于,图像插值倍数大小设定为10。
10.一种基于质心投影轨迹拟合的锥束CT温度漂移校正***,用于CT投影重建,其特征在于,包含:数据处理模块、曲线估计模块、温漂计算模块和图像校正模块,其中,
数据处理模块,用于依据CT原始投影图像数据,得到被投影物体质心投影轨迹曲线;
曲线估计模块,用于通过局部区间长度比例及区间步长,利用该比例下的被投影物体局部区间质心投影估计未发生温度漂移时被投影物体整体投影轨迹曲线;
温漂计算模块,用于依据估计得到的整体投影轨迹曲线和计算得到的投影轨迹曲线两者的曲线偏差获取温度漂移量;
图像校正模块,用于通过温度漂移量将原始投影图像亚像素级平移来实现CT图像偏移校正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110290331.3A CN112884862B (zh) | 2021-03-18 | 2021-03-18 | 基于质心投影轨迹拟合的锥束ct温度漂移校正方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110290331.3A CN112884862B (zh) | 2021-03-18 | 2021-03-18 | 基于质心投影轨迹拟合的锥束ct温度漂移校正方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112884862A true CN112884862A (zh) | 2021-06-01 |
CN112884862B CN112884862B (zh) | 2022-11-01 |
Family
ID=76041054
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110290331.3A Active CN112884862B (zh) | 2021-03-18 | 2021-03-18 | 基于质心投影轨迹拟合的锥束ct温度漂移校正方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112884862B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8023767B1 (en) * | 2008-03-10 | 2011-09-20 | University Of Rochester | Method and apparatus for 3D metal and high-density artifact correction for cone-beam and fan-beam CT imaging |
CN102426696A (zh) * | 2011-10-24 | 2012-04-25 | 西安电子科技大学 | 光学投影断层成像运动伪影校正方法 |
US20140270439A1 (en) * | 2013-03-14 | 2014-09-18 | Wisconsin Alumni Research Foundation | System and method for simulataneous image artifact reduction and tomographic reconstruction |
EP3009862A2 (en) * | 2014-10-16 | 2016-04-20 | Uniwersystet Slaski w Katowicach | Method of determining distortion and/or distortion correction for projection image obtained during computer tomography routine |
WO2018126335A1 (zh) * | 2017-01-03 | 2018-07-12 | 苏州海斯菲德信息科技有限公司 | 基于小球模体的锥束ct***几何参数评价及校正方法 |
CN109685877A (zh) * | 2018-12-27 | 2019-04-26 | 重庆大学 | 一种基于自适应投影图像特征区域匹配的微纳ct焦点漂移校正方法 |
CN111882624A (zh) * | 2020-06-19 | 2020-11-03 | 中国人民解放军战略支援部队信息工程大学 | 基于多采集序列的纳米ct图像运动伪影校正方法及装置 |
-
2021
- 2021-03-18 CN CN202110290331.3A patent/CN112884862B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8023767B1 (en) * | 2008-03-10 | 2011-09-20 | University Of Rochester | Method and apparatus for 3D metal and high-density artifact correction for cone-beam and fan-beam CT imaging |
CN102426696A (zh) * | 2011-10-24 | 2012-04-25 | 西安电子科技大学 | 光学投影断层成像运动伪影校正方法 |
US20140270439A1 (en) * | 2013-03-14 | 2014-09-18 | Wisconsin Alumni Research Foundation | System and method for simulataneous image artifact reduction and tomographic reconstruction |
EP3009862A2 (en) * | 2014-10-16 | 2016-04-20 | Uniwersystet Slaski w Katowicach | Method of determining distortion and/or distortion correction for projection image obtained during computer tomography routine |
WO2018126335A1 (zh) * | 2017-01-03 | 2018-07-12 | 苏州海斯菲德信息科技有限公司 | 基于小球模体的锥束ct***几何参数评价及校正方法 |
CN109685877A (zh) * | 2018-12-27 | 2019-04-26 | 重庆大学 | 一种基于自适应投影图像特征区域匹配的微纳ct焦点漂移校正方法 |
CN111882624A (zh) * | 2020-06-19 | 2020-11-03 | 中国人民解放军战略支援部队信息工程大学 | 基于多采集序列的纳米ct图像运动伪影校正方法及装置 |
Non-Patent Citations (4)
Title |
---|
SHINKOOK CHOI等: "A Hybrid Approach to Reduce Cone-Beam Artifacts for a Circular Orbit Cone-Beam CT System", 《IEEE ACCESS》 * |
张峰等: "锥束CT圆轨迹半覆盖扫描的几何校正", 《光学精密工程》 * |
杨民等: "三维CT重建几何参数的非线性最小二乘估计", 《北京航空航天大学学报》 * |
王珏等: "基于线框模型的锥束CT几何参数校正方法", 《仪器仪表学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112884862B (zh) | 2022-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101473538B1 (ko) | 엑스선 이미지들 내의 마커 위치들로부터의 환자 모션 벡터들의 추출 | |
US8731266B2 (en) | Method and system for correcting artifacts in image reconstruction | |
EP2523604B1 (en) | A computed tomography imaging process and system | |
CN107764846B (zh) | 一种正交直线扫描的cl成像***及分析方法 | |
US9826942B2 (en) | Correcting and reconstructing x-ray images using patient motion vectors extracted from marker positions in x-ray images | |
Sun et al. | An iterative projection‐based motion estimation and compensation scheme for head x‐ray CT | |
CN109685877B (zh) | 一种基于自适应投影图像特征区域匹配的微纳ct焦点漂移校正方法 | |
CN107978005B (zh) | 一种基于保边界扩散和平滑的有限角ct图像重建算法 | |
US9824468B2 (en) | Dictionary learning based image reconstruction | |
CN105225208A (zh) | 一种计算机断层成像金属伪影校正方法及装置 | |
CN110057847B (zh) | Tr层析扫描投影重排方法及装置 | |
Li et al. | Sparse CT reconstruction based on multi-direction anisotropic total variation (MDATV) | |
Niebler et al. | Projection‐based improvement of 3D reconstructions from motion‐impaired dental cone beam CT data | |
Ametova et al. | Software-based compensation of instrument misalignments for X-ray computed tomography dimensional metrology | |
CN112884862B (zh) | 基于质心投影轨迹拟合的锥束ct温度漂移校正方法及*** | |
Tan et al. | An interval subdividing based method for geometric calibration of cone-beam CT | |
CN110308165B (zh) | 使用Pi线优化的计算机断层扫描投影的几何对齐、样本移动校正以及强度归一化 | |
CN109360252B (zh) | 基于投影变换的锥束cl投影数据等效转换方法 | |
Jun et al. | Alignment theory of parallel-beam computed tomography image reconstruction for elastic-type objects using virtual focusing method | |
Powierza et al. | Fully experiment-based evaluation of few digital volume correlation techniques | |
Delgado-Friedrichs et al. | PI-line difference for alignment and motion-correction of cone-beam helical-trajectory micro-tomography data | |
Dennerlein et al. | Geometric jitter compensation in cone-beam CT through registration of directly and indirectly filtered projections | |
US20210012544A1 (en) | Method and apparatus for metal artifact reduction in industrial 3-dimensional cone beam computed tomography | |
CN107025634B (zh) | 负字典学习 | |
CN111899312A (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 |