本发明专利技术提供了针对非理想圆形轨迹且运动重复性差的CBCT校准方法，包括：在模体表面上镶嵌钢珠，记录各钢珠的几何位置，将模体放置在CBCT中心，对模体进行照射，采集图像，得到二维平片；通过图像处理算法将平片进行处理，自动将金属点的坐标标记出来；根据模体中已知的钢珠标记点的坐标值和第二步中自动标记出来的钢珠标记点序号，得到各金属标记点在模体坐标系中的坐标，并找到对应的探测器平板坐标系中的坐标值；计算矩阵；对每个Gantry角度θ计算修正矩阵；将各角度的采集的投影图像用修正矩阵进行修正，然后再进行三维重建。本发明专利技术提供的方法可以对非理想圆形轨迹且运动重复性差的CBCT进行校准，实现了CT图像精确重建。实现了CT图像精确重建。实现了CT图像精确重建。
【技术实现步骤摘要】
针对非理想圆形轨迹且运动重复性差的CBCT校准方法
[0001]本专利技术涉及CBCT校准
，具体为针对非理想圆形轨迹且运动重复性差的CBCT校准方法。
技术介绍
[0002]计算机断层成像技术是通过对物体进行不同角度的射线投影测量获取物体横截面信息的成像技术。按照射线束采集方式，可分为扇形束CT和锥形束CT两大类。其中扇形束螺旋CT目前已经在医疗领域广泛应用，锥形束CT(Cone
‑
beamCT，CBCT)由于剂量低，空间分辨率高和扫描速度快等优点已成为行业热点。临床应用中，CBCT图像重建通常采用FDK算法来对图像进行重建。FDK算法有两个前提条件：第一、射线源焦点与探测器中心的连线必须通过旋转中心线，且垂直于探测器平面；第二、旋转轴应与探测器的中心列方向平行。然而在实际搭建CBCT系统时，由于机械误差的存在，CBCT系统不能严格满足上述2个前提条件，因此存在一定的几何偏移。如果不对偏移进行校正，重建的图像会产生严重的伪影，影响图像质量。因此，几何参数的精确求解(即几何校正)是实现CT图像精确重建的前提。
[0003]目前业内使用的几何校正方法需要设计校正模体，通过多角度或者一个角度下获取的校正体模的投影数据，利用解析方法获得CT系统的几何参数，然后根据该参数进行后续图像重建。各种CBCT校正方法需要系统运动为标准的圆周运动，且系统的运动具有高度重复性(即运动完全稳定)。在此前提下，校正算法对系统进行刚性位移校正。然而在一些实际应用场景中，受限于机械精度，系统运动不规则，且重复性较差，即系统运动具有随机派偏差，如在运动颠簸状态下的影像扫描采集。针对这类特殊应用场景，目前业内尚未有成熟的解决方案。
技术实现思路
[0004]本专利技术所解决的技术问题在于提供针对非理想圆形轨迹且运动重复性差的CBCT校准方法，以解决技术背景中所提到的问题。为实现上述目的，本专利技术提供如下技术方案：针对非理想圆形轨迹且运动重复性差的CBCT校准方法，具体包括以下内容：
[0005](1)在圆柱形模体表面上镶嵌钢珠，记录各钢珠的几何位置，将模体放置在CBCT中心，对圆柱形模体进行照射，采集图像，得到二维平片；
[0006](2)通过图像处理算法将平片进行处理，自动将金属点的坐标标记出来，而测量结果为探测器平板坐标系中的坐标值
[0007](3)根据模体中已知的钢珠标记点的坐标值和第二步中自动标记出来的钢珠标记点序号，得到第二步中各金属标记点在模体坐标系中的坐标p
i
，并找到p
i
对应的探测器平板坐标系中的坐标值q
i
；
[0008](4)计算矩阵M
t
：设坐标在平板上的投影在平板坐标系中为则
[0009][0010]其中为模体坐标系到CBCT坐标系的转换矩阵，M
p
为CBCT坐标系中的点到平板坐标系的投影矩阵，即
[0011][0012]设为从平片图像中提取的在平板上的投影坐标，则解
[0013][0014]可得对于模体坐标系中任意点p
i
，M
t
·
p
i
为该点在CBCT坐标系中的坐标；
[0015](5)对每个Gantry角度θ计算M(θ)；
[0016](6)将各角度的采集的投影图像用修正矩阵M(θ)进行修正，然后再进行三维重建。
[0017]优选地，所述(1)中圆柱形模体表面嵌有16个钢珠，且呈螺旋状分布。
[0018]优选地，所述(1)中采集图像的方法为：在0
°‑
359
°
每个角度分别对圆柱形模体进行照射，采集图像，得到360张二维平片。
[0019]优选地，所述(2)中图像处理算法为霍夫变换。
[0020]优选地，所述(5)的计算方法为：则在Gantry角度θ下，可得方程：
[0021]q
i
＝M
p
·
M(θ)
·
M
t
·
p
i
[0022]设为从平片图像中提取的在平板上的投影坐标，则解
[0023][0024]得M(θ)，即为在Gantry为θ角度下的修正矩阵。
[0025]与现有技术相比，本专利技术提供的方法可以对非理想圆形轨迹且运动重复性差的CBCT进行校准，实现了CT图像精确重建。
附图说明
[0026]图1为实施例中圆柱形模体表面上镶嵌钢珠示意图；
[0027]图2为实施例中图像采集示意图。
[0028]为了更清楚地说明本专利技术的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本专利技术的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
具体实施方式
[0029]为了使本专利技术的技术手段、创作特征、工作流程、使用方法达成目的与功效易于明白了解，下面将结合本专利技术实施例，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，及
进行非创造性的扩展而得出的其它结论，都属于本专利技术保护的范围。
[0030]实施例
[0031]针对非理想圆形轨迹且运动重复性差的CBCT校准方法，具体包括以下内容：
[0032](1)在圆柱形模体表面嵌有16个钢珠，且呈螺旋状分布，如图1所示，记录各钢珠的几何位置，将模体放置在CBCT中心，在0
°‑
359
°
每个角度分别对圆柱形模体进行照射，采集图像，如图2所示，得到360张二维平片；
[0033](2)通过图像处理算法(霍夫变换)将平片进行处理，自动将金属点的坐标标记出来，而测量结果为探测器平板坐标系中的坐标值
[0034](3)根据模体中已知的钢珠标记点的坐标值和第二步中自动标记出来的钢珠标记点序号，得到第二步中各金属标记点在模体坐标系中的坐标p
i
，并找到p
i
对应的探测器平板坐标系中的坐标值q
i
；
[0035](4)计算矩阵M
t
：设坐标在平板上的投影在平板坐标系中为则
[0036][0037]其中为模体坐标系到CBCT坐标系的转换矩阵，M
p
为CBCT坐标系中的点到平板坐标系的投影矩阵，即
[0038][0039]设为从平片图像中提取的在平板上的投影坐标，则解
[0040][0041]可得对于模体坐标系中任意点p
i
，M
t
·
p
i
为该点在CBCT坐标系中的坐标；
[0042](5)本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.针对非理想圆形轨迹且运动重复性差的CBCT校准方法，其特征在于，具体包括以下内容：(1)在圆柱形模体表面上镶嵌钢珠，记录各钢珠的几何位置，将模体放置在CBCT中心，对圆柱形模体进行照射，采集图像，得到二维平片；(2)通过图像处理算法将平片进行处理，自动将金属点的坐标标记出来，而测量结果为探测器平板坐标系中的坐标值(3)根据模体中已知的钢珠标记点的坐标值和第二步中自动标记出来的钢珠标记点序号，得到第二步中各金属标记点在模体坐标系中的坐标p
i
，并找到p
i
对应的探测器平板坐标系中的坐标值q
i
；(4)计算矩阵M
t
：设坐标在平板上的投影在平板坐标系中为则其中为模体坐标系到CBCT坐标系的转换矩阵，M
p
为CBCT坐标系中的点到平板坐标系的投影矩阵，即设为从平片图像中提取的在平板上的投影坐标，则解可得对于模体坐标系中任意点p
i
，M
t
·
p
i
为该点在CBCT坐标系中的坐标；(5)对每个Gantry角度θ计算M(θ)；(6)将各角...
【专利技术属性】
技术研发人员：陈蔚海，周云，马丽，
申请(专利权)人：太丛信息科技上海有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。