发明名称 扇束2D-CT扫描系统转台旋转中心偏移量的测量方法
摘要 本发明公开了一种扇束2D-CT扫描系统转台旋转中心偏移量的测量方法,是将扇束投影重排为平行束投影,根据投影角相差π的平行束投影相似的原理,将重排后的平行束正弦图中[0°,180°]和[180°,360°]对应的投影信息取出,获得第一正弦图和第二正弦图;将第二正弦图水平翻转后形成的第三正弦图与第一正弦图进行互相关运算,从而求出转台旋转中心偏移量τ0。本发明的方法不需要制作专门的校正试件,直接利用被扫描断层的正弦图进行测量。有效保证了扇束2D-CT系统转台旋转中心偏移量的测量精度,提高了重建图像的质量。
申请公布号 CN101825433B 申请公布日期 2011.09.28
申请号 CN201010139235.0 申请日期 2010.03.31
申请人 北京航空航天大学 发明人 杨民;刘永瞻;梁丽红
分类号 G01B11/02(2006.01)I 主分类号 G01B11/02(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 李有浩
主权项 1.一种扇束2D-CT扫描系统转台旋转中心偏移量的测量方法,其特征在于包括有下列测量步骤:步骤一:在2D-CT扫描系统中,线阵列探测器(6)和射线源焦点(1)不动,让转台(4)绕旋转中心点O作圆周转动;所述转台(4)绕旋转中心点O的转动,使得线阵列探测器(6)能够采集到不同投影角下的被扫描断层(5)的投影数据信息,将所有投影角下的投影数据合成为一幅扇束射线扫描正弦图p<sup>fan</sup>(θ,s),θ表示转台的旋转角度,也代表投影角,S代表探测器单元坐标;步骤二:将扇束射线扫描正弦图p<sup>fan</sup>(θ,s)重排为平行束扫描正弦图p(θ,s),p(θ,s)为二维矩阵,该矩阵的每一行代表每一投影角下的投影数据,θ=β时的投影数据记为p<sub>β</sub>(s),p<sub>β</sub>(s)为p(θ,s)的某一行;步骤三:将正弦图p(θ,s)中θ∈[0°180°]范围内的对应投影数据取出合成第一正弦图<img file="FSB00000514424800011.GIF" wi="145" he="63" />将正弦图p(θ,s)中θ∈[180°360°]范围内的对应投影数据取出合成第二正弦图<img file="FSB00000514424800012.GIF" wi="289" he="64" />和<img file="FSB00000514424800013.GIF" wi="118" he="63" />均为二维矩阵,<img file="FSB00000514424800014.GIF" wi="118" he="63" />为p(θ,s)的前半部分,<img file="FSB00000514424800015.GIF" wi="119" he="64" />为p(θ,s)的后半部分;步骤四:将<img file="FSB00000514424800016.GIF" wi="118" he="64" />进行水平翻转,得到第三正弦图<img file="FSB00000514424800017.GIF" wi="300" he="64" />与<img file="FSB00000514424800018.GIF" wi="117" he="64" />的关系为<img file="FSB00000514424800019.GIF" wi="402" he="63" />N为线阵列探测器(6)中探测单元的总数目;步骤五:将第一正弦图<img file="FSB000005144248000110.GIF" wi="117" he="64" />与第三正弦图<img file="FSB000005144248000111.GIF" wi="118" he="64" />的对应行进行互相关运算,那么对应行的互相关函数最大值对应的坐标就是<img file="FSB000005144248000112.GIF" wi="118" he="64" />与<img file="FSB000005144248000113.GIF" wi="117" he="63" />对应行之间的位移差数列记为d(i),i=1~M,M表示<img file="FSB000005144248000114.GIF" wi="118" he="74" />或<img file="FSB000005144248000115.GIF" wi="118" he="64" />的高度;步骤六:求位移差数列d(i)的均值τ′,该均值τ′的1/2等于中心射线(2)偏移转台(4)的旋转中心点O的偏移量τ<sub>0</sub>,即<img file="FSB000005144248000116.GIF" wi="349" he="121" />
地址 100191 北京市海淀区学院路37号