专利名称：一种层析成像几何参数的校准方法
技术领域：
本发明涉及层析成像技术领域，具体涉及一种层析成像几何参数的校准方法，尤其适合于锥形束X射线层析成像系统的几何参数校准。
背景技术：
锥形束X射线层析成像(cone beam computed tomography, CBCT)是临床医学、生物医学等领域非常重要的诊断与研究工具，它包括乳腺CT、牙科CT、四肢CT、小动物CT等。 CBCT系统的几何参数对获取高分辨率、低伪影的重建图像起着决定性作用，几何参数的些许偏差将会带来成像质量的严重下降。因此，需要对CBCT系统的几何参数进行校准。
现有的CBCT系统几何参数的校准方法可以分为两类
一类是离线校准，即在对被测样品扫描之前，先对专用模体成像得到系统的几何参数，然后保证系统固定不动，再进行样品试验。这类方法中有些需要由几十个钢球构成的已知相对几何结构的精密模体(美国专利US_2005_0117708_A1 ；美国专利 US_2006_0245628_A1 ；中国发明专利:ZL 200510045796. 3 ；Physics In Medicine and Biology 54(2009) 1633-1660 ；Physics In Medicine and Biology 54(2009)7239-7261； Medical Physics 37(2010)3844-3854；Medical Physics 38(2011)2829-2840 ；Medical Physics 31(2004)3242-3266 ；Medical Physics 32 (2005) 968-983)，有些需要精密加工的图形模具(美国专利US_2007_0041508_A1 ；欧洲专利W0_2010_050970_A1)、圆杆 (中国发明专利=ZL 200910023137. 8)、线模体(中国发明专利=ZL 200910087131.7)、 已知相对垂直距离的两个滚珠模型(中国发明专利ZL 200910188615.0 ；中国发明专利ZL200910079277. 7 ；中国发明专利ZL 200610066252. X ；Physics In Medicine and Biology 45(2000)3489-3508 ；Medical Physics 33(2006) 1695-1706 ；IEEE Transactions on Information Technology in Biomedicine 15 (2011) 655-660，)。该类方法对模体的要求较高，这些模体的加工与制作精度关乎着几何校准的准确度，实际操作中很容易引入各种机械或者测量误差。
另一类是在线校准，即利用被测样品自身的投影图计算出CBCT系统的几何参数， 不需要制作任何模体。这是因为在实际系统中，由于机械结构受环境的影响会逐渐出现偏差，致使系统的几何参数在动态的变化，特别在高分辨系统中，微小的系统几何参数的变化都会带来性能的急剧下降；此外，还有部分系统需要频繁调整视场，放大倍率等，这使得离线校准的方法繁琐并且不可靠。中国发明专利ZL 200810136662. 6通过对投影图进行二值化处理后求图像重心的方法得到了旋转中心位置；中国发明专利ZL200810112193.4利用正弦图找出物体边缘的投影点在探测器阵列中的位置，从而求出了旋转中心位置；文章 Medical Physics 36 (2009) 48_58利用正弦图求解出了旋转轴位置和平面内的倾角；文章 Physics In Medicine and Biology53 Q008) 3841-3861利用正弦图求得与几何参数相关联的代价函数，然后采用单纯形算法求得旋转轴位置及探测器的两个倾角，但该方法不适合大锥角的系统，计算结果对感兴趣区域的选择较敏感，操作复杂，采用的单纯形算法易陷入局部极小从而得不到正确的结果；文章Medical Physics 38 (2011) 4934-4945是把几何参数代入到CT重建算法中，通过比较重建出的图像边缘的锐化程度得到了旋转轴位置、探测器的两个倾角及射线源到探测器的距离，该方法在求解过程中，是采用了穷举的方法，没有做优化，而且每一次调整几何参数都需要重建一次图像，过程复杂耗时。现有的在线校准方法的整体问题就是目标函数的建立非常困难、求解复杂，且精度较离线校准稍低。发明内容
有鉴于此，本发明的目的在于提供一种层析成像几何参数的校准方法，用于在不需要精密模体的情况下，实现高精度校准。
本发明提供了一种层析成像几何参数的校准方法，包括以下步骤
采集层析成像中的全角度投影图像，对所述采集到的全角度投影图像分别进行负对数运算，将每个角度的负对数运算后的图像进行叠加，得到叠加图像，利用所述叠加图像建立与几何参数相关联的目标函数，利用单纯形-模拟退火算法对所述目标函数进行全局最小值优化。
本发明利用了层析成像全角度投影图的叠加图像的对称性与几何参数的数学关系，采用了单纯形-模拟退火全局最优化算法精确求解出了图像重建所需要的关键几何参数，并利用求出的几何参数重建出高分辨率、低伪影的正确结果。本发明不需要任何精密的模体就可以实现多个几何参数的标定；由于利用了叠加图像，因此对噪声不敏感，抗噪声能力强；此外，本发明采用的单纯形-模拟退火全局最优化算法，可以避免陷入局部极小值， 保证了多个几何参数同时求解的精度与准确度。
图1为本发明实施例提供的层析成像几何参数校准系统的结构图2为本发明实施例提供的层析成像几何参数校准方法的流程图3为本发明实施例提供的建立目标函数的流程图4为本发明实施例中偏移探测器上像素点坐标与理想探测器平面上像素点坐标之间的对应关系示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面结合附图对本发明作进一步的详细描述。
本发明实施例提供的层析成像几何参数校准系统的结构及待求几何参数如图1 所示，射线源焦点ι到探测器2之间的距离为d，射线源焦点1在探测器2上的投影为(Utl， K)，其中Utl为旋转轴位置，％为中平面位置，探测器2绕y轴旋转的角度记为n (即探测器的倾斜角skew)，探测器2绕Z(V)轴旋转的角度记为θ (即探测器的滚转角roll)。射线源包括X射线，Y射线，可见光与近红外光。
本发明实施例提供的层析成像几何参数校准方法的流程如图2所示，该方法包括以下步骤
步骤201、定义右手笛卡尔坐标系x-y-z。其中ζ轴是旋转台的旋转轴，y轴垂直于旋转轴且穿过射线源焦点，χ轴垂直于y-ζ平面且穿过ι轴与ζ轴的交点，探测器平面上横向定义为U轴，纵向定义为V轴。
步骤202、采集层析成像中的投影图像，即物体或者机架每旋转一个角度，探测器就采集一幅投影图像，直至等角度旋转完360度。
步骤203、对层析成像中采集到的全角度投影图像分别进行负对数运算。全角度是指投影图像覆盖了 0 360度的旋转角度。运算方法由公式1给出
权利要求
1.一种层析成像几何参数的校准方法，其特征在于，包括以下步骤采集层析成像中的全角度投影图像，对所述采集到的全角度投影图像分别进行负对数运算，将每个角度的负对数运算后的图像进行叠加，得到叠加图像，利用所述叠加图像建立与几何参数相关联的目标函数，利用单纯形-模拟退火算法对所述目标函数进行全局最小值优化。
2.根据权利要求1所述的校准方法，其特征在于，该方法进一步包括以下步骤设定右手笛卡尔坐标系χ-y-z，ζ轴为旋转台的旋转轴，y轴垂直于旋转轴且穿过射线源焦点，χ轴垂直于y-ζ平面且穿过y轴与ζ轴的交点，探测器平面上横向设定为u轴，纵向设定为ν轴。
3.根据权利要求1所述的校准方法，其特征在于，所述层析成像所使用的射线源具体包括X射线、Y射线、可见光或近红外光。
4.根据权利要求1所述的校准方法，其特征在于，所述全角度投影图像为该投影图像覆盖了 0 360度的旋转角度。
5.根据权利要求1至4中任意一项所述的校准方法，其特征在于，所述建立目标函数的步骤具体包括建立空间中点的坐标与探测器平面上像素点之间的对应关系； 将所述叠加图像绕y轴旋转角度n； 再将所述叠加图像绕ζ轴旋转角度Θ ；根据获得的绕y轴和ζ轴旋后新的叠加图像，计算所述新的叠加图像关于旋转轴对称的最小二乘误差，得到最终的以几何参数为变量的多参数目标函数。
6.根据权利要求5所述的校准方法，其特征在于，所述目标函数为 f ( η，θ，u0，v0，d) = minimize (err)其中，η，θ，Utl，Vtl，d分别为探测器的倾斜角skew、探测器的滚转角roll、旋转轴位置、 中平面位置以及射线源焦点到探测器的距离，err为误差值。
7.根据权利要求6所述的校准方法，其特征在于，所述对目标函数进行全局最小值优化的方法具体包括设置初始条件为 IU0-UinI ^ 40pixels I V0-Vin I < 40pixels I n in| < 4degrees I d-dinI ^ IOOpixels I θ - θ in I ^ 4degrees其中，(uin，vin，nin，din，ein)表示输入的理想的初始几何参数。
8.根据权利要求7所述的校准方法，其特征在于，该方法适用于千伏特kv和兆伏特MV 的工业CT，单光子发射层析成像SPECT以及正电子发射层析成像PET。
全文摘要
本发明提供了一种层析成像几何参数的校准方法，涉及层析成像技术领域。该方法包括以下步骤采集层析成像中的全角度投影图像，对所述采集到的全角度投影图像分别进行负对数运算，将每个角度的负对数运算后的图像进行叠加，得到叠加图像，利用所述叠加图像建立与几何参数相关联的目标函数，利用单纯形-模拟退火算法对所述目标函数进行全局最小值优化。本发明不需要任何精密的模体就可以实现多个几何参数的标定；由于利用了叠加图像，因此对噪声不敏感，抗噪声能力强；此外，本发明采用的单纯形-模拟退火全局最优化算法，可以避免陷入局部极小值，保证了多个几何参数同时求解的精度与准确度。
文档编号A61B6/03GK102488528SQ20111040257
公开日2012年6月13日 申请日期2011年12月7日 优先权日2011年12月7日
发明者孟远征, 杨孝全, 骆清铭, 龚辉 申请人:华中科技大学