发明内容
本发明提供一种基于CT定位片确定病人尺寸、位置,可以利用CT定位片获得准确的病人尺寸、位置信息,消除因病床投影造成的误差。
本发明还提供一种基于CT定位片确定X射线经过病人衰减的方法,可以准确获得X射线经过病人的衰减,有效排除病床投影对X射线衰减精度造成的误差。
为实现上述目的,本发明的基于CT定位片确定病人尺寸、位置的方法通过如下技术方案实现:一种基于CT定位片确定病人尺寸、位置的方法,包括如下步骤:a)对病人进行CT定位扫描获取CT定位片;b)去除CT定位片中病床与空气的投影值获得病人的投影值;c)根据病人的投影值获得病人尺寸和位置。
进一步地,所述CT定位片按照其上像素点的像素值大小分为仅含空气的投影区域、含有空气和病床的投影区域以及含有空气、病床和病人的投影区域。
进一步地,所述仅含空气的投影区域的像素点的像素值不小于零且不大于第一阈值T1,所述含有空气和病床的投影区域的像素点的像素值大于第一阈值T1,且不大于第二阈值T2,所述含有空气、病床和病人的投影区域的像素点的像素值大于第二阈值T2。
进一步地,所述获得病人投影值的过程包括如下步骤:将CT定位片中含空气和病床的投影区域的各像素点的投影值进行平滑滤波后拟合获得拟合后的投影值;将CT定位片中含空气、病床和病人的投影区域的各像素点的投影值与拟合后的投影值进行差值计算获得病人的投影值。
进一步地,所述获得病人尺寸、位置的过程包括如下步骤:根据病人的横截面获得病人横截面所在曲线;根据病人的投影值获得病人的投影值对应的间隔最大的两个检测器通道的位置;根据球管焦点的位置、两个检测器通道的位置获得通过该两个检测器通道的X射线的路径;根据X射线路径与病人横截面所在曲线获得病人的尺寸与位置。
进一步地,所述获得病人投影值对应的间隔最大的两个检测器通道的位置的过程包括:确定病人的投影值中大于第三阈值T3的检测器通道范围;从检测器通道范围中选出检测器通道编号最小的一个检测器通道m与检测器通道编号最大的一个检测器通道n;根据检测器通道m、检测器通道n获得该两个检测器通道的位置。
进一步地,所述第三阈值T3等于0.1。
进一步地,所述通过检测器通道m的X射线路径为其中,(xed1,yed1)为检测器通道m的位置,(xf,yf)为球管焦点的位置,所述通过检测器通道n的X射线路径为其中,(xed2,yed2)为检测器通道n的位置,(xf,yf)为球管焦点的位置。
进一步地,所述病人横截面所在曲线为其中,a为病人平躺在病床上时身体的宽度尺寸;x0为病人平躺在病床上时病人横截面中心至旋转中心在X轴方向上的偏移距离;b为病人平躺在病床上时病人的高度且b=0.5*P/μh,P为所有检测器通道上最大的病人投影值,μh为人体平均的X射线线性衰减系数;y0为病人平躺在病床上时病人横截面中心至旋转中心在Y轴方向上的偏移距离,且y0=yb-b,yb为当前床高时病床上表面的坐标。
进一步地,所述获得病人尺寸、位置的过程包括:根据通过检测器通道m的X射线路径与病人的横截面曲线仅具有一个交点以及通过检测器通道n的X射线路径与病人的横截面曲线仅具有一个交点获得病人尺寸、位置。
本发明基于CT定位片确定X射线经过病人衰减的方法通过如下技术方案实现:一种基于CT定位片确定X射线经过病人衰减的方法,包括:确定X射线经过病人的衰减面积、X射线经过各个检测器通道的路径、病床上表面所在曲线以及病床下表面所在曲线;根据X射线经过病人的衰减面积以及根据上述方法获得的病人的尺寸获得X射线经过病人的衰减系数;根据经过各个检测器通道的X射线路径、病人的横截面曲线获得各X射线穿过病人的路径长度;根据经过各检测器通道的X射线路径、病床上表面所在曲线、病床下表面所在曲线获得各X射线穿过病床的路径长度;根据X射线经过病人的衰减系数、各X射线穿过病人的路径长度、各X射线穿过病床的路径长度获得各X射线经过病人后到达检测器通道的衰减。
进一步地,所述X射线经过病人的衰减系数μ通过如下公式获得:其中,Sa为X射线经过病人的衰减面积,a为病人平躺在病床上时身体的宽度尺寸,b为病人平躺在病床上时病人的高度且b=0.5*P/μh,P为所有检测器通道上最大的病人投影值,μh为人体平均的X射线线性衰减系数。
进一步地,所述X射线经过各检测器通道的路径是根据球管焦点位置与各检测器通道位置获得的。
进一步地,所述X射线经过各检测器通道的路径为其中,(xedi,yedi)为第i个检测器通道的位置,i取m至n之间的整数,(xf,yf)为球管焦点的位置。
进一步地,所述获得某一X射线穿过病人的路径长度包括:根据该X射线经过的检测器通道的路径与病人的横截面曲线获得该X射线经过检测器通道的路径与病人横截面曲线的两个交点位置;根据两个交点的位置获得两个交点之间的距离,所述该两个交点之间的距离即为该X射线穿过病人的路径长度。
进一步地,所述病床上表面所在曲线为x=Rbcosθ,y=(yb-Rb)+Rbsinθ,其中,Rb为病床的半径,yb为病床的床高,αd为病床的扇角。
进一步地,所述病床下表面所在曲线为x=(Rb+Lb)cosθ,y=(yb+Lb-Rb)+Rbsinθ,其中,Rb为病床的半径,yb为病床的床高,Lb为病床的厚度,αd为病床的扇角。
进一步地,所述获得某一X射线穿过病床的路径长度的过程包括:根据该X射线经过的检测器通道的路径与病床的上表面所在曲线获得该X射线与病床的上表面的交点位置;根据该X射线经过的检测器通道的路径与病床的下表面所在曲线获得该X射线与病床的下表面的交点位置;根据两交点位置获得两交点之间的距离,所述该X射线穿过病床的路径长度即为两交点之间的距离。
进一步地,所述某一X射线经过病人后到达检测器通道的衰减Ai通过如下公式获得:Ai=μDi+μbdi,其中,μ为X射线经过病人的衰减系数,Di为该X射线经过病人的路径长度,μb为病床对X射线的平均衰减系数,di为该X射线经过病床的路径长度。
本发明对比现有技术有如下的有益效果:本发明基于CT定位片确定病人的尺寸、位置的方法通过从CT定位片去除病床与空气的投影值获得的病人的投影值、进而获得X射线经过病人的投影值对应的间隔最大的两个检测器通道的路径,从而根据该两条X射线路径与病人横截面所在曲线分别仅有一个交点获得病人的尺寸与位置,可以消除现有技术因CT定位片中因含有病床投影而造成的误差;本发明基于CT定位片确定X射线经过病人衰减的方法根据上述方法获得的病人的尺寸确定X射线经过病人的衰减,可以提高X射线衰减的准确度,有效排除病床投影对X射线衰减造成的误差。
具体实施方式
下面结合附图和实施例对本发明作进一步的描述。
图1为本发明基于CT定位片确定病人尺寸、位置的方法的流程示意图。
请参见图1,本发明提供的确定病人尺寸、位置的方法包括如下步骤:
步骤S1:对病人进行CT定位扫描获取CT定位片;
在实际的CT定位片扫描中,病人躺在CT机病床上,因此采集的定位片中往往包含了病床和空气的投影值,如图2所示。
步骤S2:去除CT定位片中病床与空气的投影值获得病人的投影值;
为了得到准确的病人尺寸和衰减信息,需要首先将CT定位片中空气和病床的投影去除,得到真正的病人投影。当病人的尺寸较小时,在检测器采集到的投影数据中,有可能存在三种情况:1)检测器没有采集到病人和病床的投影,即X射线传播路径只经过了空气,测量值为空气的投影;2)检测器的投影测量值为病床的投影;3)检测器的投影测量值为病人和病床的投影之和。因此,投影数据并不能直接反映病人对X射线的衰减信息,需要从投影测量值中去除空气和病床的影响。
请继续参见图2,所述CT定位片中空气、病床与病人的投影关系通过如下方式确定:所述CT定位片按照其上像素点的像素值大小分为仅含空气的投影区域、含有空气和病床的投影区域以及含有空气、病床和病人的投影区域;像素点的像素值不小于零且不大于第一阈值T1的部分为仅含空气的投影值区域,如曲线1所示;像素点的像素值大于第一阈值T1且不大于第二阈值T2的部分为含有空气和病床的投影区域,如曲线2所示;像素点的像素值大于第二阈值T2的部分为含有空气、病床和病人的投影区域,如曲线3所示。
所述CT定位片中空气和病床的投影值去除过程如下:由于CT病床的投影值变化比较平缓,先将CT定位片中含空气和病床的投影区域的各像素点的投影值进行平滑滤波后拟合获得拟合后的投影值,请参图6与图7所示;然后将CT定位片中含空气、病床和病人的投影区域的各像素点的投影值与拟合后的投影值进行差值计算获得病人的投影值,请参图3中曲线4所示。
由于CT机能够在X射线球管角度为0度,90度,180度和270度时采集定位片,上述病床的投影值确定方法可适用于球管角度为0度和180度采集的定位片,即球管位于病床的上方或下方。本发明还提供另一种确定病床投影值的方法,该方法适用于球管角度在以上任意一个角度的情况,图6为球管角度为180度时不同床高的病床投影值,图7为球管角度为90度时不同床高的病床投影值;请参见图6与图7所示,在本实施例中,需要在每个球管角度下,预先采集多个典型床高时病床的投影,进行平滑滤波并拟合后得到病床的投影,其他床高时的病床投影值可以通过插值得到,用于从包含病床投影的定位片中去除病床的投影影响,即将定位片的投影减去插值得到的病床投影,得到的病人投影可用于病人尺寸、衰减和位置信息的计算。
步骤S3:根据病人的投影值获得病人尺寸和位置。请参图8所示,所述获得病人尺寸、位置的过程包括如下步骤:
S30、根据病人的横截面获得病人横截面所在曲线;
S31、根据病人的投影值获得病人的投影值对应的间隔最大的两个检测器通道的位置;
S32、根据球管焦点的位置、两个检测器通道的位置获得通过该两个检测器通道的X射线的路径;
S33、根据X射线路径与病人横截面所在曲线获得病人的尺寸与位置。
对前面得到的病人投影值,通过有效投影值即第三阈值T3,检测器获得病人投影值对应的间隔最大的两个检测器通道的位置,从而判断出有效检测器通道范围m~n,具体过程如下:确定病人的投影值中大于第三阈值T3的检测器通道范围,所述第三阈值T3等于0.1;从检测器通道范围中选出检测器通道编号最小的一个检测器通道m与检测器通道编号最大的一个检测器通道n;根据检测器通道m、检测器通道n获得该两个检测器通道的位置。如图4中的ch(m)~ch(n)所示,其中m、n为正整数,通道m和n的坐标分别为(xed1,yed1)和(xed2,yed2);球管焦点6的坐标为(xf,yf),则通过通道m和n的X射线路径的方程分别为:
其中A1和B1分别为经过通道m的X射线路径方程的斜率和截距;A2和B2分别为经过通道n的X射线路径方程的斜率和截距。
根据通过检测器通道m的X射线路径与病人的横截面曲线仅具有一个交点以及通过检测器通道n的X射线路径与病人的横截面曲线仅具有一个交点可获得病人尺寸、位置。请参见图4,病人的横截面7所在曲线为长短半轴分别为a和b的椭圆,其中a轴与X轴平行,b轴与Y轴平行,a即为病人平躺在病床上时身体的宽度尺寸,b即为病人平躺在病床上时病人的高度,该病人对X射线的线性衰减系数为μ。则该病人横截面所在曲线为:
其中x0为病人平躺在病床上时病人横截面中心至旋转中心在X轴方向上的偏移距离;y0为病人平躺在病床上时病人横截面中心至旋转中心在Y轴方向上的偏移距离,,即x0和y0分别为该病人横截面所在曲线中心在X轴方向和Y轴方向上的坐标值。短轴b的长度可以根据病人投影值计算得到,公式为:
b=0.5*P/μh (4)
其中P为所有检测器通道的最大病人投影值,μh为人体的平均X射线线性衰减系数。另外在CT扫描中,床高是已知的,因此可以求得病人横截面所在曲线的中心坐标y0为:
y0=yb-b (5)
其中yb为当前床高时病床上表面的y坐标。
将方程(1)和(2)分别代入方程(3),可以求得中心坐标x0和病人横截面所在曲线长轴a。
根据上述方法获病人的尺寸a和以及病人横截面中心至旋转中心在X轴方向上的偏移距离x0。
请参见图5所示,本发明还提供一种基于CT定位片确定X射线经过病人衰减的方法,包括如下步骤:
S40、确定X射线经过病人的衰减面积、X射线经过各个检测器通道的路径、病床上表面所在曲线以及病床下表面所在曲线;
S41、根据X射线经过病人的衰减面积以及根据前述方法获得的病人的尺寸获得X射线经过病人的衰减系数;
S42、根据经过各个检测器通道的X射线路径、病人的横截面曲线获得各X射线穿过病人的路径长度;
S43、根据经过各检测器通道的X射线路径、病床上表面所在曲线、病床下表面所在曲线获得各X射线穿过病床的路径长度;
S44、根据X射线经过病人的衰减系数、各X射线穿过病人的路径长度、各X射线穿过病床的路径长度获得各X射线经过病人后到达检测器通道的衰减。
步骤S40中,所述X射线经过病人的衰减面积Sa通过如下公式获得:
其中,m、n为病人投影值对应的间隔最大的两个检测器通道,Pi为第i个检测器通道上的最大的病人投影值,Pi+1为第i+1个检测器通道上最大的病人投影值,Δd为相邻两个检测器通道之间的距离,SID为球管焦点6到旋转中心的距离,SDD为球管焦点6到检测器9的距离。
步骤S41中,所述X射线经过病人的衰减系数μ通过如下公式获得:其中,Sa为X射线经过病人的衰减面积,a为病人平躺在病床上时身体的宽度尺寸,b为病人平躺在病床上时病人的高度且b=0.5*P/μh,P为所有检测器通道上最大的病人投影值,μh为人体平均的X射线线性衰减系数。
所述病床上表面所在曲线为:
x=Rbcosθ,
y=(yb-Rb)+Rbsinθ (7)
其中,Rb为病床的半径,yb为病床8的床高,αd为病床的扇角。
所述病床下表面所在曲线为:
x=(Rb+Lb)cosθ
y=(yb+Lb-Rb)+Rbsomθ (8)
其中,Rb为病床的半径,yb为病床的床高,Lb为病床的厚度,αd为病床的扇角。
步骤S42中,所述X射线经过各检测器通道的路径是根据球管焦点位置与各检测器通道位置获得的。所述X射线经过各检测器通道的路径为: 其中,(xedi,yedi)为第i个检测器通道的位置,i取m至n之间的整数,(xf,yf)为球管焦点6的位置。
请参图9所示,步骤S42中,所述获得某一X射线穿过病人的路径长度包括如下步骤:
S50、根据该X射线经过的检测器通道的路径与病人的横截面曲线获得该X射线经过检测器通道的路径与病人横截面曲线的两个交点位置;
S51、根据两个交点的位置获得两个交点之间的距离,所述该两个交点之间的距离即为该X射线穿过病人的路径长度。
请参图10所示,步骤S43中,所述获得某一X射线穿过病床的路径长度的过程包括:
S60、根据该X射线经过的检测器通道的路径与病床的上表面所在曲线获得该X射线与病床的上表面的交点位置;
S61、根据该X射线经过的检测器通道的路径与病床的下表面所在曲线获得该X射线与病床的下表面的交点位置;
S62、根据两交点位置获得两交点之间的距离,所述该X射线穿过病床的路径长度即为两交点之间的距离。
在每个旋转角度下,球管焦点和每个检测器通道的空间坐标都是已知的,因此可以确定到达每个检测器通道的X射线路径:
其中,(xedi,yedi)为第i个检测器通道的位置,i取m至n之间的整数,(xf,yf)为球管焦点的位置。
将方程(9)与病人横截面所在曲线方程(3)联立,可求得X射线与病人横截面所在曲线的两个交点,进而求得X射线经过人体的路径长度Di。将方程(9)分别与方程(7)和(8)联立,可求得X射线与病床上表面、下表面的交点,进而获得X射线经过病床的路径长度di。所述某一X射线经过病人后到达检测器通道的衰减Ai通过如下公式获得:Ai=μDi+μbdi,其中,μ为X射线经过病人的衰减系数,Di为该X射线经过病人的路径长度,μb为病床对X射线的平均衰减系数,di为该X射线经过病床的路径长度。由此,可以确定任何旋转角度下,到达所有检测器通道的X射线的衰减。而在CT机的剂量调制算法中,一个基本的思想就是对X射线衰减强的身体部位和投影方向使用大的剂量,而对X射线衰减弱的身体部位和投影方向使用小的剂量。因此,本发明所提出的计算方法可以作为CT机剂量调制方法的依据。
虽然本发明已以较佳实施例揭示如上,然其并非用以限定本发明,任何本领域技术人员,在不脱离本发明的精神和范围内,当可作些许的修改和完善,因此本发明的保护范围当以权利要求书所界定的为准。