[发明专利]一种实用的板状构件层析扫描装置射束倾角标定方法

摘要

本发明公开了一种实用的板状构件层析扫描装置射束倾角标定方法。本方法基于层析扫描装置扫描台上的一点旋转360度时，其在探测器成像面的投影轨迹为一椭圆的原理，获取一目标体360度旋转范围内的投影序列，利用图像处理与最小二乘拟合法得到目标体椭圆投影轨迹的短半轴与长半轴之比，从而得到射束倾角的值。本发明标定方法具有结果唯一、抗噪性强、精度高的优点，从而有效保证了板状结构层析的扫描成像装置对输出层析图像的精度，减少重建伪影的产生。

主权项

一种实用的板状构件层析扫描装置射束倾角标定方法，所述板状构件层析扫描装置包括有射线源(1)、扫描平台(6)、平板探测器(8)、弧形导轨(9)；其特征在于该标定方法包括有下列实施步骤：步骤一：选取扫描平台(6)与射线源(1)之间的间距D1‑D2大于球形目标体与扫描平台的旋转中心点的距离的4倍；D1表示第一距离，D2表示第二距离；步骤二：将一球形目标体固定在扫描平台(6)上，且所述球形目标体与扫描平台(6)的旋转中心点O6的距离为10～100mm；步骤三：启动射线源(1)和扫描平台(6)，并使扫描平台(6)绕旋转轴(4)在360度范围内旋转，在每间隔1度～5度的旋转角度下，平板探测器(8)采集得到球状目标体的DR投影图像；步骤四：对每幅DR投影进行对数变换得到变换投影pj(x,y)：pj(x,y)=lnmax(pj0(x,y))pj0(x,y)]]>j＝1,2,3…Tx表示成像坐标系XOY下X轴上的坐标变量；y表示成像坐标系XOY下Y轴上的坐标变量；表示第j幅原始DR投影；表示第j幅原始投影的最大值；pj(x,y)表示第j幅对数变换后的变换投影；T表示第三步采集到的球形目标体的总投影幅数；步骤五：在图像重建与图像处理单元中对第四步对数变换后的每幅DR图像进行图像分割，分割后的图像背景灰度值为0，球状目标体投影的灰度值保持原来的值；步骤六：在图像重建与图像处理单元中求取第五步中图像分割后的每幅DR图像中球状目标体投影的质心坐标值(xj,yj)：xj=Σx=1mΣy=1nx×fj(x,y)ΣxΣyfj(x,y)yj=Σx=1mΣy=1ny×fj(x,y)ΣxΣyfj(x,y)]]>m表示DR投影的长度；n表示DR投影的高度；fj(x,y)表示DR图像的二维函数，其中，x表示成像坐标系XOY下X轴上的坐标变量，y表示成像坐标系XOY下Y轴上的坐标变量；xj表示第j幅DR图像中球状目标体投影质心在X轴上的坐标；yj表示第j幅DR图像中球状目标体投影质心在Y轴上的坐标；步骤七：在图像重建与图像处理单元中对第六步中得到的所有质心坐标进行最小二乘拟合，拟合曲线为一椭圆，椭圆方程可表示为：Ax2+Bxy+Cy2+Dx+Ey+F＝0x表示成像坐标系XOY下X轴上的坐标变量；y表示成像坐标系XOY下Y轴上的坐标变量；A表示椭圆方程中变量x的二次项系数；B表示椭圆方程中变量xy的二次项系数；C表示椭圆方程中变量y的二次项系数；D表示椭圆方程中变量x的一次项系数；E表示椭圆方程中变量y的一次项系数；F表示椭圆方程的常数项；步骤八：建立误差函数f(A,B,C,D,E,F)：f(A,B,C,D,E,F)=Σj=1n(Axj2+Bxjyj+Cyj2+Dxj+Eyj+F)2]]>n表示球状目标体DR投影质心的数目；f(A,B,C,D,E,F)代表误差函数；xj表示第j幅DR图像中球状目标体投影质心在X轴上的坐标；yj表示第j幅DR图像中球状目标体投影质心在Y轴上的坐标；对误差函数f(A,B,C,D,E,F)进行求偏导，并令可得到一线性方程组，求解该方程组，即可得到椭圆方程的系数A、B、C、D、E、F的值，利用这些值，得出椭圆的长半轴与短半轴的计算表达式：Ra=-2FA+C-B2+(A-CF)2Rb=-2FA+C+B2+(A-CF)2]]>步骤九：求取椭圆短半轴与长半轴之比该比值等于射束倾角α的余弦值即为板状构件层析扫描装置的射束倾角α的近似值为