一种确定c型臂断层成像系统中几何标定参数的方法
【专利摘要】本发明公开了一种确定C型臂断层成像系统中几何标定参数的方法，具体而言是用螺旋线模体直接确定C-arm?CT系统中9个几何标定参数，属于医疗成像领域。该方法利用螺旋线的几何特性，将其巧妙地转换为一系列的平行四边形。然后获得这些平行四边形对角线的交点，用来确定模体坐标系坐标轴的投影，最后将其集成到相应的标定算法中。本发明基于广泛应用的螺旋线模体，不需要额外加工标定模体，并且是一种解析的方法，充分利用螺旋线模体上标志点的信息，具有较高的鲁棒性及准确性。
【专利说明】一种确定C型臂断层成像系统中几何标定参数的方法
【技术领域】
[0001]本发明涉及一种C型臂断层成像（c-arm CT)系统中的离线标定方法，更具体地说，它涉及一种基于螺旋线模体直接确定C-arm CT系统中几何标定参数的方法，属于医疗成像领域。
【背景技术】
[0002]集成了平板探测器的C型臂断层成像（C-arm CT)系统，通过获得一系列二维X射线投影数据去重建物体的断层图像。C型臂成像系统主要由X射线源，平板探测器，C型臂构架组成。相对于传统CT而言，c-arm CT具有移动方便，价格便宜等优点，尤其是其开放性设计，能广泛地应用在介入治疗中的术前规划，术中检测及结果评估中。
[0003]但是在实际应用中，因为机械制造，系统集成及重力下垂等因素的影响，C-arm CT系统的实际轨道会偏离预设轨道。这些不规则的偏差会降低图像分辨率，或导致在重建图像中出现严重伪影，影响医生判断。
[0004]目前针对此问题，螺旋线模体由于具有良好的3维空间信息，被广泛地应用在实际的几何标定中。 此模体主体为中空的圆柱形的树脂，上面嵌有按螺旋线均匀排列的一系列小钢珠，这些小钢珠即为标志点。研究者采用的方法大多基于投影矩阵，即利用模体上标志点的3维空间信息及投影图像的2维信息的对应关系来标定系统。但此方法没有对成像系统进行具体地描述，如X射线源位置，探测器的倾斜及旋转角度，源到探测器的距离等，并且在系统偏离较大的情况下，以此方法为基础得到重建图像质量可能会不理想。
[0005]另外一些研究者能获得系统的具体标定参数，但是此时螺旋线模体不再适用，需要设计专门的标定模体，这样加大了额外工作量。
【发明内容】
[0006]本发明的目的是为了解决上述问题，提出一种基于螺旋线模体直接确定C型臂断层成像系统中几何标定参数的方法。该方法利用常用的螺旋线模体，用解析方法获取了成像系统的全部9个几何标定参数，具有较高的鲁棒性及准确性。
[0007]本发明利用螺旋线的几何特性，将螺旋线上的标志点按不同的分类方法分类，对每一组而言，每四个标志点定义一个平行四边形。在不考虑投影点重叠的情况下，这四个标志点的投影点确定一个四边形。这个投影得到的四边形对角线的交点即为平行四边形的对角线的交点的投影。具体的技术方案如下：
[0008]第一步，将螺旋线模体放在C型臂置物台上，用C型臂扫描，获取投影数据。
[0009]第二步，确定螺旋线上标志点在平板探测器上的投影位置。
[0010]第三步，利用螺旋线的几何特性，确定模体坐标系z轴，X轴及两条平行于I轴的投影。
[0011]第四步，确定模体坐标系原点在平板探测器上的投影。
[0012]第五步，获取中间参数：a、b、c ;[0013]第六步，通过中间参数获得系统的9个几何标定参数：探测器的旋转角度（3个参数：俯仰角，偏转角，倾斜角)，x射线源到探测器距离（I个参数)，X射线源位置（3个参数：x方向，y方向，Z方向），模体坐标系原点在探测器上的投影（2个参数：水平方向和竖直方向)。
[0014]本发明的优点在于：
[0015](I)该方法能获得系统的全部9个几何标定参数，从而能具体地描述此成像系统，对机械制造有指导意义。并且获得独立的几何标定参数，为后期将其集成到相应的重建算法中提供了可能性；
[0016](2)该方法基于广泛应用的螺旋线模体，不需要额外加工标定模体；
[0017](3)该方法是一种解析方法，不基于数学约化。并且充分利用螺旋线模体上标志点的信息，保证了方法的鲁棒性及准确性。
【专利附图】
【附图说明】
[0018]图I是螺旋线模体示意图。
[0019]图2是确定几何标定参数的流程图。
[0020]图3是确定坐标系z轴投影的示意图。
[0021]图中：【具体实施方式】
[0022]下面将结合附图和实施例对本发明作进一步的详细说明。
[0023]本发明是一种确定C型臂断层成像系统中几何标定参数的方法，如图2所示，包括以下几个步骤：
[0024]第一步，将螺旋线模体放在C型臂断层成像系统中的置物台上，用C型臂扫描，获取螺旋线模体在平板探测器上的投影图像。
[0025]螺旋线模体如图I所示，本发明中模体镶嵌小钢珠（即标志点）的树脂部份跟通用的不大一致，本发明采用螺旋管道代替了常用的中空圆柱形，这样就降低了制造模体所需要的树脂，从而降低了制造此模体的成本(采用快速成型的加工方法)。小钢珠的排列跟通用的螺旋线状模体是一致的，这样不会影响本方法的通用性。
[0026]螺旋线模体的放置要保证螺旋线的中心轴与C型臂系统旋转轴的方向大概一致。用C型臂连续扫描一周，角度范围设定根据机型的不同，适用范围的不同调整。满足进行断层成像要求的最小角度范围为180° +最大扇角的一半。
[0027]最后，得到X射线穿透螺旋线模体在平板探测器上的投影图像。
[0028]第二步，确定螺旋线上标志点在平板探测器上的投影位置。
[0029]为了确定标志点投影在探测器上的位置，首先用canny边界检测的方法获得标志点的边界，接着采用基于奇异值分解的椭圆拟合算法确定投影图像上标志点中心。本发明采用椭圆的中心作为投影中心，根据相关文献，这两者间的误差可以忽略。
[0030]第三步，利用螺旋线的几何特性，确定模体坐标系z轴、X轴及两条平行于I轴的轴的投影。
[0031]螺旋线上标志点的排列满足半径为R，螺距为2 h的螺旋线方程，第I个标志点位置的定义如下：
[0032]r1 = (Rsin Ψ1, Rcos Ψ1, h Ψ1)
[0033]Ψ1 = 1XΔΨ(1 = -N1, —N1+1, ? ? ?，N2-1, N2)
[0034]其中：A Ψ是螺旋线上标志点的角度间隔，由方程Δ Ψ = Ji /N给出。N的取值为整数，理论上来说，N的取值范围可以从1到无穷大。但在实际应用中，当螺旋线长度保持固定，随着N取值的增大，标志点的数量也随着增多，此时标志点的投影可能会在投影图像上出现重叠，不利于后期处理。当N取值太小时，标志点数目过少，会影响算法精度。所以，在实际应用中，N的取值范围为3≤N≤18是比较合理的。
[0035]参数N1和N2的取值为整数，它们的选取跟螺旋线中轴线的总长度和参考点的选取有关。假设均匀排列的标志点构成的螺旋线中轴线的长度为L，那么N1+N2 = L/(h Δ Ψ）-1。一般选取位于螺旋线中部或接近中部的标志点做为参考点，此参数点的角度为0,即Ψt0 =
O。在给出螺旋线中轴线的总长度及参考点的选取位置后，可以确定出N1和N2的取值。
[0036]( 1)确定坐标系z轴投影的方法，如图3所示，从螺旋线上选取4个标志点，形成平行四边形Pkz如下：
[0037]
【权利要求】
1.一种确定C型臂断层成像系统中几何标定参数的方法，包括以下几个步骤： 第一步，将螺旋线模体放在C型臂断层成像系统中的置物台上，用C型臂扫描，获取螺旋线模体在平板探测器上的投影图像； 第二步，确定螺旋线上标志点在平板探测器上的投影位置； 第三步，利用螺旋线的几何特性，确定模体坐标系Z轴、X轴及两条平行于I轴的轴的投影； 螺旋线上标志点的排列满足半径为R，螺距为2 Jih的螺旋线方程，第I个标志点位置为：
2.根据权利要求1所述的一种确定C型臂断层成像系统中几何标定参数的方法，所述的螺旋线模体采用螺旋管道，小钢珠位于螺旋管道中。
3.根据权利要求1所述的一种确定C型臂断层成像系统中几何标定参数的方法，第二步中，首先用canny边界检测的方法获得标志点的边界，然后采用基于奇异值分解的椭圆拟合算法确定投影图像上标志点中心。
【文档编号】A61B6/03GK103549971SQ201310551531
【公开日】2014年2月5日 申请日期:2013年11月7日 优先权日:2013年11月7日
【发明者】李德玉, 徐敏, 樊瑜波, 蒲放, 张弛 申请人:北京航空航天大学