CN102805635A - 一种优化扫描流程方法和装置 - Google Patents
一种优化扫描流程方法和装置 Download PDFInfo
- Publication number
- CN102805635A CN102805635A CN2011101452106A CN201110145210A CN102805635A CN 102805635 A CN102805635 A CN 102805635A CN 2011101452106 A CN2011101452106 A CN 2011101452106A CN 201110145210 A CN201110145210 A CN 201110145210A CN 102805635 A CN102805635 A CN 102805635A
- Authority
- CN
- China
- Prior art keywords
- angle
- projection
- probe access
- data
- projection angle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及医疗设备领域,特别是一种优化扫描流程的方法,所述方法包括:根据原始投影数据组来判断待检对象发生运动的位置;对所述位置再次进行扫描。本发明还提供一种优化扫描流程的装置。采用本发明能自动判断出待检对象在本次扫描中发生运动的具***置,并仅对所述具***置处再次进行扫描,从而更加精确,并减少了待检对象接收的X射线辐射剂量。
Description
技术领域
本发明涉及医疗设备领域,特别是一种优化扫描流程的方法和装置。
背景技术
随着医疗水平的进步,X射线计算机断层成像(Computed Tomography,CT)设备得到了广泛的应用。当对待检对象进行CT扫描时,旋转机架绕Z轴旋转360度,并同时带动包括X射线源和探测器的数据采集***同步围绕待检对象旋转360°,以从[0,360°]的投影角度来获取X射线投影数据,其中Z轴为CT检查床进出旋转机架的水平方向。在每个投影角度处,探测器的M个通道同时接收X射线投影数据,并将其传送至CT的图像重建***用于生成重建图像。然而由于人体不可避免会发生运动,所以在对人体进行CT扫描时,重建图像会因人体的运动而存在运动伪影,从而影响重建图像的质量。
在目前的CT扫描中,CT图像重建***在每次重建图像后可以显示该次重建的所有图像,据此观察者(如医生或CT技术人员)会查看所有的这些重建图像上,用肉眼观察每幅重建图像中是否有运动伪影。由于运动伪影导致重建图像模糊,对观察者观察病灶有很大影响,所以通常在观察到运动伪影后就要对待检对象重新进行扫描,以得到没有运动伪影的重建图像。为了减少待检对象接收的扫描剂量,观察者会根据自己的经验来判断待检对象发生运动的位置,并将扫描线设置在该位置处,自此处开始再次对待检对象进行扫描。不过这种方式存在人为误差。
发明内容
本发明的目的在于提供一种优化扫描流程的方法和装置,其可自动判断出待检对象发生运动的具***置,从而更加精确;并仅对所述具***置处再次进行扫描,以减少待检对象接收的辐射剂量。
有鉴于此,本发明提出一种优化扫描流程的方法,所述方法包括:根据原始投影数据组来判断待检对象发生运动的位置;对所述位置再次进行扫描。以通过自动方式获取待检对象发生运动的位置,减少待检对象接收的辐射剂量。
根据本发明的一个实施例,所述对所述位置再次进行扫描包括:若为螺旋扫描,则先根 据所述位置来计算需再次扫描的范围,再对所述范围再次进行扫描;若为顺序扫描,则直接对所述位置再次进行扫描。以通过自动方式获取对待检对象进行螺旋扫描的再次扫描范围,从而更加精确。
根据本发明的一个实施例,所述根据原始投影数据组来判断待检对象发生运动的位置包括:获取所述原始投影数据组的补偿数据组;计算所述原始投影数据组与所述补偿数据组的距离矩阵;取所述距离矩阵的某一列,找出该列中最小值所对应的行值,若该行值与该列值的差值的绝对值大于一预设阈值,则根据该行值得到所述待检对象发生运动的具***置。可见本发明能直接从原始投影数据组中就判断出待检对象在本次扫描中发生运动的具***置,所以准确率和可靠性更高。
根据本发明的一个实施例,所述距离矩阵为欧几里得距离矩阵。
根据本发明的再一个实施例,所述获取补偿数据组包括:读取所述原始投影数据组中与一投影数据对应的第一投影角度和第一探测器通道;根据所述第一投影角度和所述第一探测器通道来计算一第二投影角度和第二探测器通道,其中所述第二投影角度下的第二探测器通道与所述第一投影角度下的第一探测器通道位于同一X射线传播路径上;将所述第二投影角度下第二探测器通道所采集到的投影数据作为所述补偿数据组中的数据。
根据本发明的又一个实施例,所述计算第二投影角度和第二探测器通道包括:根据所述第一探测器通道按如下公式来计算所述X射线传播路径偏离X射线束中心的角度:θ=β-n1*Δβ,其中,n1为所述第一探测器通道,Δβ为所述探测器通道的角度步距,θ为所述X射线传播路径与X射线束中心的夹角,β为所述X射线束的一半张角;根据所述夹角按如下公式来计算所述第二探测器通道:n2=(β+θ)/Δβ,其中,n2为所述第二探测器通道,且所述第二探测器通道与所述第一探测器通道关于所述探测器的中心呈镜像对称;根据所述夹角按如下公式来计算所述第二投影角度:α2=α1±pi-2θ,其中,α2为所述第二投影角度,α1为所述第一投影角度。
根据本发明的又一个实施例,所述第一投影角度的范围为[0,180°+2β]。这样数据采集***只需围绕待检对象进行[0,180°+2β]的投影即可得到所述原始投影数据,相比360°的投影方式减少了待检对象接收的辐射剂量,同时减少了原始投影数据组的数据量,减少了数据存储空间。
根据本发明的又一个实施例,所述预设阈值为20。
本发明还提供一种计算机程序,使计算机执行所述优化扫描流程的方法。
本发明还提供一种计算机可读的记录介质,存储了所述的计算机程序。
本发明还提供一种优化扫描流程的装置,所述装置包括:一运动位置判断组件,用于根据原始投影数据组来判断待检对象发生运动的位置;一再扫描执行组件;用于对所述位置再次进行扫描。以通过自动方式获取待检对象发生运动的位置,减少待检对象接收的辐射剂量。
根据本发明的一个实施例,所述再扫描执行组件包括:一再扫描范围计算模块,用于在螺旋扫描时根据所述位置来计算需再次扫描的范围;一再扫描执行模块,用于对所述范围再次进行扫描。以通过自动方式获取对待检对象进行螺旋扫描的再次扫描范围,从而更加精确。
根据本发明的一个实施例,所述运动位置判断组件包括:一补偿数据获取模块,用于获取所述原始投影数据组的补偿数据组;一距离矩阵计算模块,用于计算所述原始投影数据组与所述补偿数据组的距离矩阵;一运动位置判断模块,取所述距离矩阵的某一列,找出该列中最小值所对应的行值,若该行值与该列值的差值的绝对值大于一预设阈值,则根据该行值得到所述待检对象发生运动的具***置。可见本发明能直接从原始投影数据组中就判断出待检对象在本次扫描中发生运动的具***置,所以准确率和可靠性更高。
根据本发明的再一个实施例,所述补偿数据获取模块包括:一读取单元,用于读取所述原始投影数据组中与一投影数据对应的第一投影角度和第一探测器通道;一计算单元,用于根据所述第一投影角度和所述第一探测器通道来计算一第二投影角度和一第二探测器通道,其中所述第二投影角度下的第二探测器通道与所述第一投影角度下的第一探测器通道位于同一X射线传播路径上;一补偿数据获取单元,其将所述第二投影角度下第二探测器通道所采集到的投影数据作为所述补偿数据组中的数据。
根据本发明的又一个实施例,所述计算单元包括:一角度计算子单元,用于根据所述第一探测器通道按如下公式来计算所述X射线传播路径偏离X射线束中心的角度:θ=β-n1*Δβ,其中,n1为所述第一探测器通道,Δβ为所述探测器通道的角度步距,θ为所述X射线传播路径与X射线束中心的夹角,β为所述X射线束的一半张角;一第二通道计算子单元,用于根据所述夹角按如下公式来计算所述第二探测器通道:n2=(β+θ)/Δβ,其中,n2为所述第二探测器通道,且所述第二探测器通道与所述第一探测器通道关于所述探测器的中心呈镜像对称;一第二角度计算子单元,用于根据所述夹角按如下公式来计算所述第二投影角度:α2=α1±pi-2θ,其中,α2为所述第二投影角度,α1为所述第一投影角度。
从上述方案中可以看出,由于本发明能自动判断出待检对象在本次扫描中发生运动的 具***置,从而更加精确。而且本发明仅对所述具***置处再次进行扫描,减少了待检对象接收的X射线辐射剂量。同时本发明的原始投影数据为数据采集***围绕待检对象进行[0,180°+2β]的投影,相比360°的投影方式减少了待检对象接收的辐射剂量,也减少了原始投影数据组的数据量,减少了数据存储空间。
附图说明
下面将通过参照附图详细描述本发明的实施例,使本领域技术人员更清楚本发明的上述及其它特征和优点,附图中:
图1a为本发明优化顺序扫描流程的方法示意图。
图1b为本发明优化螺旋扫描流程的方法示意图。
图2为本发明获取补偿数据的原理示意图。
图3a为待检对象未发生运动时本发明距离矩阵中某列数据的示意图。
图3b为待检对象发生运动时本发明距离矩阵中某列数据的示意图,图中同时示意了待检对象发生运动具***置。
图4为本发明优化顺序扫描流程的装置示意图。
图5为本发明运动位置判断组件的示意图。
图6为本发明补偿数据获取模块的示意图。
图7为本发明再扫描执行组件的示意图。
附图标记
步骤101~110优化扫描流程的步骤
1X射线源
3探测器
10本发明优化扫描流程的装置
11运动位置判断组件
12再扫描执行组件
111补偿数据获取模块
112距离矩阵计算模块
113运动位置判断模块
221读取单元
222计算单元
223补偿数据获取单元
121再扫描范围计算模块
122再扫描执行模块
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下举实施例对本发明进一步详细说明。
以下以CT设备为例,具体说明本发明优化扫描流程的方法和装置。
对目前的双排CT而言,通常数据采集***围绕待检对象同步旋转一周的时间少于1s,如为0.6s,而人体发生运动频率通常为一次1秒,所以在扫描过程中,并非每次扫描得到的原始投影数据组中都含有待检对象的运动数据。此外,每次扫描过程中,数据采集***是根据一定的角度间隔来采集X射线投影数据的,即读片(reading),如每次扫描需要读片2000次,由于人体在扫描过程中并非一直在动,所以并非每次读片的投影数据都含有运动数据,除非这个人患有帕金森综合症,这种症状的病人的运动频率非常快,需要扫描时间很短的高端CT来对其进行扫描。
对于顺序扫描而言,现有的扫描流程如下:
首先,根据重建图像来人为判断这些重建图像中是否有运动伪影,例如在一幅重建图像pic1中发现运动伪影。
接着,CT***根据pic1获得与pic1对应的检查床位置pos1。由于顺序扫描是在指定的复数个检查床位置处依次进行扫描,所以能根据pic1及检查床的进床速度获得与pic1对应的检查床位置pos1。
然后,CT***根据输入的指令再增加一次扫描,并在在上次扫描结束位置附近设置一条扫描指示线,然而这并不是应该执行再扫描的准确位置。
再然后,CT***根据输入指令移动扫描指示线至pos1处。这样就得到另一幅重建图像pic2。
最后,并用pic2来替换pic1用于进行图像分析。
但是该流程中需要人为发现图像中的运动伪影,比较费时,也存在人为误差。
而螺旋扫描的扫描流程则与顺序扫描的扫描流程有所不同。在该扫描流程中CT***根据pic1获得与pic1对应的检查床位置pos1和X射线源投影角度δ1,然后在pos1和δ1的位置之前和之后都需要进行再扫描,其中再扫描的起点线相同,只是扫描范围box1和box2不同。需要将这两个扫描范围box1和box2的终点线都移动到pos1和δ1的位置,然后进行再扫描,并将再扫描得到的图像用来替换有运动伪影的原图像pic1。
然而实际应用中,观察者并非会按上述方式操作。例如对于有10幅图像的顺序扫描,如果其中3幅图像有运动伪影,那么为了保证重建图像的质量,一般会重新扫描这10幅图像。而实际上没有运动伪影的7幅图像是无需再次扫描的,这样待检对象就接收了额外的X射线辐射。
有鉴于此,本发明提出一种优化扫描流程的方法。用于优化再扫描的流程,以提高判断图像中运动伪影的准确率,并减少待检对象接收的X射线剂量。
如图1a所示,对顺序扫描而言,所述方法包括:
步骤101,读取原始投影数据组R中与一投影数据v1对应的第一投影角度α1和第一探测器通道n1。
本发明是基于如图2所示的数据补偿的原理来进行的。图2所示数据采集***包括一种X射线源1(F1或F2)和探测器3,探测器3的通道数为M,其起始探测器通道为C-1,结束探测器通道为C-M。其中,X射线源1发射-X射线束,X射线束张角为2β,探测器3与X射线源1一起围绕一待检对象顺时针同步旋转,即X射线源从F1处顺时针旋转到F2处,以从复数个视点来采集X射线投影数据。理论上,在同一X射线传播路经F1F2上,由于X射线从F1到F2方向经过待检对象的衰减与从F2到F1方向经过待检对象的衰减相同,所以在投影角度为α1(从X轴正半轴旋转到F1点的旋转角度)时,探测器3的一个通道n1上接收到的X射线投影数据v1应该等于投影角度为α2(从X轴正半轴旋转到F2点的旋转角度)时,探测器3的另一个通道n2接收的X射线投影数据v2,即v1=v2。本发明据此来获取补充数据组R_C的数据。
步骤102,根据第一探测器通道n1按如下公式(1)来计算X射线传播路径F1F2偏离X射线束中心的角度θ:
θ=β-n1*Δβ (1)
其中,Δβ为探测器3中相邻两个通道间的角度步距,β为所述X射线束的一半张角。
步骤103,根据夹角θ按如下公式(2)来计算第二探测器通道n2:
n2=(β+θ)/Δβ (2)
这里,可以将公式(1)变形为如下公式(3):
n1=(β-θ)/Δβ (3)
比较公式(2)和(3)可以得出:第一探测器通道n1与第二探测器通道n2关于探测器3的中心是镜像对称的。
步骤104,根据夹角θ按如下公式(4)来计算第二投影角度α2:
α2=α1±pi-φ (4)
由于φ=2θ,所以
α2=α1±pi-2θ (5)
图2的实施例中α2=α1-pi-2θ。
步骤105,将在第二投影角度α2下第二探测器通道n2所采集到的投影数据v2作为原始投影数据组R的补偿数据组R_C中的数据。
采用上述步骤101~105,无需包括数据采集***围绕待检对象旋转一周(360°),而只需进行[0,180°+2β]的投影即可获得补偿数据组R_C中的全部数据。于是,本发明中第一投影角度范围为[0,180°+2β],这样根据公式(5)可以得到第二投影角度范围也为[0,180°+2β],这样就减少了待检对象接收的辐射剂量,同时减少了原始投影数据的数据量,减少了数据存储空间。
步骤106,计算原始投影数据组R与补偿数据组R_C的距离矩阵S。
根据本发明的一个实施例,所述距离矩阵为欧几里得距离矩阵。
在CT扫描中,数据采集***是根据一定的角度间隔来采集X射线投影数据的,即读片。如在一次扫描过程中本发明将[0,180°+2β]分为N个角度间隔,如N=2000,则数据采集***围绕待检对象每旋转((180°+2β)/2000)°,探测器3就读取一次X射线投影数据,总共需要读片2000次。
假设距离矩阵S为一N行N列的矩阵,其中第i行第j列的数据可以是原始投影数据组R中第i次读片的数据与补偿数据组R_C中第j次读片的数据的欧几里得(Euclidean)距离。 也就是说,距离矩阵S中第i行的数据是原始投影数据组R中第i次读片的数据分别与补偿数据组R_C中第1~N次读片的数据的欧几里得距离,距离矩阵S中第j列的数据是补偿数据组R_C中第j次读片的数据分别与原始投影数据组R中第1~N次读片的数据的欧几里得距离。
当然距离矩阵S第i行第j列的数据也可以是补偿数据组R_C中第i次读片的数据与原始投影数据组R中第j次读片的数据的欧几里得(Euclidean)距离。
步骤107,取距离矩阵S的某一列j,找出该列中最小值所对应的行值i,若该行值i与该列值j的差值的绝对值大于一预设阈值,则根据该行值i得到待检对象发生运动的具***置。
假设j=1000,则距离矩阵S中第1000列数据为补偿数据组R_C中第1000次读片的数据分别与原始投影数据组R中第1~N次读片的数据的欧几里得距离。理论上,距离矩阵S中第j列的最小数据为Sj,j,即i=j。然而实际上,由于各种原因如机械误差(例如探测器通道的大小、X射线源焦点到扫描中心的距离、X射线源焦点到探测器的距离、X射线束张角等的误差)、条状伪影检测的准确率、飞焦点、旋转机架的速度等,使得S中第j列的最小数据并不一定为Sj,j,而是在Sj,j附近的一个预设阈值范围内。根据本发明的一个实施例,所述预设阈值为20。如图3a所示为待检对象未发生运动时S中第1000列数据的示意图,其中横坐标为S的行值i,纵坐标为Si,j的相对大小,图中在i=1000附近有一个圆圈,该圆圈示意了一预设阈值范围,即当S中第1000列的最小值在该圆圈内,则认为在该次扫描过程中待检对象没有发生运动。反之,若第1000列的最小值不在该圆圈内,如图3b所示,则认为该次扫描过程中待检对象发生了运动。
如图3b所示S中第1000列最小值为i=1250左右,其相比i=1000附近的圆圈有很大的偏移,可知在本次扫描过程中约第1250次读片时待检对象发生了运动,且纵坐标有一个明显的跳跃(sharp break)。可见本发明能自动判断出待检对象在本次扫描中发生运动的具***置,相比现有顺序扫描再次扫描的方法而言,避免了人工误差,更加精确。
步骤108,对所述位置再次进行扫描。
根据步骤107得到的行值i得到在本次扫描过程中约第i次读片时待检对象发生了运动,据此得到检查床的位置posi。将检查床移动到posi处,并使X射线源处于第i次读片的投影角度δi=(180°+2β)/N*i,对该处待检对象再次进行扫描。相比现有顺序扫描再次扫描的方法而言,减少待检对象接收的X射线剂量。
图1b所示为本发明优化螺旋扫描流程的方法:其包括其中步骤101~107与图1a相同,不同的是,在根据步骤107得到待检对象在本次扫描中发生运动的具***置后,螺旋扫描的步骤为:
步骤109,根据所述位置来计算需再次扫描的范围。
根据步骤107得到的行值i来获得第i次读片时检查床位置posi和X射线源投影角度δi,将其作为两个再次扫描范围box3和box4的起点。
步骤110,对所述范围再次进行扫描。启动再次扫描,使box3和box4的终点线分别都移动到posi和δ1的位置,这样就实现了再扫描。
本发明还提供一种计算机程序,使计算机执行所述优化扫描流程的方法。
本发明还提供一种计算机可读的记录介质,存储用于使一机器执行如本文所述的优化扫描流程的方法指令。具体地,可以提供配有存储介质的***或者装置,在该存储介质上存储着实现上述实施例中任一实施例的功能的软件程序代码,且使该***或者装置的计算机(或CPU或MPU)读出并执行存储在存储介质中的程序代码。这种情况下,从存储介质读取的程序代码本身可实现上述实施例中任何一项实施例的功能,因此程序代码和存储程序代码的存储介质构成了本发明的一部分。
用于提供程序代码的存储介质实施例包括软盘、硬盘、磁光盘、光盘(如CD-ROM、CD-R、CD-RW、DVD-ROM、DVD-RAM、DVD-RW、DVD+RW)、磁带、非易失性存储卡和ROM。可选择地,可以由通信网络从服务器计算机上下载程序代码。
此外,应该清楚的是,不仅可以通过执行计算机所读出的程序代码,而且可以通过基于程序代码的指令使计算机上操作的操作***等来完成部分或者全部的实际操作,从而实现上述实施例中任意一项实施例的功能。
此外,可以理解的是,将由存储介质读出的程序代码写到***计算机内的扩展板中所设置的存储器中或者写到与计算机相连接的扩展单元中设置的存储器中,随后基于程序代码的指令使安装在扩展板或者扩展单元上的CPU等来执行部分和全部实际操作,从而实现上述实施例中任一实施例的功能。
本发明还提供一种优化扫描流程的装置,如图4所示,所述装置10包括:一运动位置判断组件11和一再扫描执行组件12。其中运动位置判断组件11,用于根据原始投影数据组来判断待检对象发生运动的位置,如图5所示,其包括:
一补偿数据获取模块111,用于获取原始投影数据组R的补偿数据组R_C。
具体地,如图6所示补偿数据获模块11件包括:
一读取单元221,用于读取原始投影数据组R中与一投影数据v1对应的第一投影角度α1和第一探测器通道n1。
一计算模块222,用于根据第一投影角度α1和第一探测器通道n1来计算一第二投影角度α2和第二探测器通道n2,其中在第二投影角度α2下的第二探测器通道n2与在第一投影角度α1下的第一探测器通道n1位于同一X射线传播路径F1F2上。具体地包括:
一角度计算子单元,用于根据第一探测器通道n1按公式(1)来计算X射线传播路径F1F2偏离X射线束中心的角度θ;
一第二通道计算子单元,用于根据夹角θ按公式(2)来计算第二探测器通道n2;
一第二角度计算子单元,用于根据夹角θ按公式(5)来计算第二投影角度α2。
一补偿数据获取模块223,其将计算单元222得到的第二投影角度α2下第二探测器通道n2所采集到的投影数据v2作为补偿数据组R_C中的数据。
一距离矩阵计算模块112,用于计算原始投影数据组R与补偿数据组R_C的距离矩阵S。
一运动位置判断模块113,取所述距离矩阵中的第j列,找出该列中最小值所对应的行值i,若该行值i与该列值j的差值的绝对值大于一预设阈值,则根据该行值得到所述待检对象发生运动的具***置。根据本发明的一个实施例,预设阈值为20。具体判断方式已在本发明的方法中做了详细说明,这里不再赘述。
如图7所示,再扫描执行组件12包括:
一再扫描范围计算模块121,用于在螺旋扫描时根据所述位置来计算需再次扫描的范围;
一再扫描执行模块122,用于对所述范围再次进行扫描。
本发明涉及医疗设备领域,特别是一种根据原始投影数据组来判断待检对象是否发生运动的方法,所述方法包括:获取所述原始投影数据组的补偿数据组;计算所述原始投影数据组与所述补偿数据组的距离矩阵;取所述距离矩阵的某一列,找出该列中最小值所对应的行值,若该行值与该列值的差值的绝对值大于一预设阈值,则所述原始投影数据组中存在运动伪影。本发明还提供一种根据原始投影数据组来判断待检对象是否发生运动的装置。采用本发明能直接从[0,180°+2β]的原始投影数据组中就判断出待检对象在本次扫描中是否发生运动,并进一步得到发生运动的具***置,同时减少了待检对象接收的X射线辐射剂量。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (15)
1.一种优化扫描流程的方法,所述方法包括:
根据原始投影数据组来判断待检对象发生运动的位置;
对所述位置再次进行扫描。
2.根据权利要求1所述的方法,其特征在于,所述对所述位置再次进行扫描包括:
若为螺旋扫描,则先根据所述位置来计算需再次扫描的范围,
再对所述范围再次进行扫描;
若为顺序扫描,则直接对所述位置再次进行扫描。
3.根据权利要求1所述的方法,其特征在于,所述根据原始投影数据组来判断待检对象发生运动的位置包括:
获取所述原始投影数据组的补偿数据组;
计算所述原始投影数据组与所述补偿数据组的距离矩阵;
取所述距离矩阵的某一列,找出该列中最小值所对应的行值,若该行值与该列值的差值的绝对值大于一预设阈值,则根据该行值得到所述待检对象发生运动的具***置。
4.根据权利要求3所述的方法,其特征在于,所述距离矩阵为欧几里得距离矩阵。
5.根据权利要求3所述的方法,其特征在于,所述获取补偿数据组包括:
读取所述原始投影数据组中与一投影数据对应的第一投影角度和第一探测器通道;
根据所述第一投影角度和所述第一探测器通道来计算一第二投影角度和第二探测器通道,其中所述第二投影角度下的第二探测器通道与所述第一投影角度下的第一探测器通道位于同一X射线传播路径上;
将所述第二投影角度下第二探测器通道所采集到的投影数据作为所述补偿数据组中的数据。
6.根据权利要求5所述的方法,其特征在于,所述计算第二投影角度和第二探测器通道包括:
根据所述第一探测器通道按如下公式来计算所述X射线传播路径偏离X射线束中心的角度:
θ=β-n1*Δβ,
其中,n1为所述第一探测器通道,Δβ为所述探测器通道的角度步距,θ为所述X射线传播路径与X射线束中心的夹角,β为所述X射线束的一半张角;
根据所述夹角按如下公式来计算所述第二探测器通道:
n2=(β+θ)/Δβ,
其中,n2为所述第二探测器通道,且所述第二探测器通道与所述第一探测器通道关于所述探测器的中心呈镜像对称;
根据所述夹角按如下公式来计算所述第二投影角度:
α2=α1±pi-2θ,
其中,α2为所述第二投影角度,α1为所述第一投影角度。
7.根据权利要求6所述的方法,其特征在于,所述第一投影角度的范围为[0,180°+2β]。
8.根据权利要求3所述的方法,其特征在于,所述预设阈值为20。
9.一种计算机程序,使计算机执行如权利要求1-8中任一权利要求所述的方法。
10.一种计算机可读的记录介质,存储了如权利要求9所述的计算机程序。
11.一种优化扫描流程的装置,所述装置(10)包括:
一运动位置判断组件(11),用于根据原始投影数据组来判断待检对象发生运动的位置;
一再扫描执行组件(12);用于对所述位置再次进行扫描。
12.根据权利要求11所述的装置,其特征在于,所述再扫描执行组件(12)包括:
一再扫描范围计算模块(121),用于在螺旋扫描时根据所述位置来计算需再次扫描的范围;
一再扫描执行模块(122),用于对所述范围再次进行扫描。
13.根据权利要求11所述的装置,其特征在于,所述运动位置判断组件(11)包括:
一补偿数据获取模块(111),用于获取所述原始投影数据组的补偿数据组;
一距离矩阵计算模块(112),用于计算所述原始投影数据组与所述补偿数据组的距离矩阵;
一运动位置判断模块(113),取所述距离矩阵的某一列,找出该列中最小值所对应的行值,若该行值与该列值的差值的绝对值大于一预设阈值,则根据该行值得到所述待检对象发生运动的具***置。
14.根据权利要求3所述的装置,其特征在于,所述补偿数据获取模块(111)包括:
一读取单元(221),用于读取所述原始投影数据组中与一投影数据对应的第一投影角度和第一探测器通道;
一计算单元(222),用于根据所述第一投影角度和所述第一探测器通道来计算一第二投影角度和一第二探测器通道,其中所述第二投影角度下的第二探测器通道与所述第一投影角度下的第一探测器通道位于同一X射线传播路径上;
一补偿数据获取单元(223),其将所述第二投影角度下第二探测器通道所采集到的投影数据作为所述补偿数据组中的数据。
15.根据权利要求14所述的装置,其特征在于,所述计算单元(222)包括:
一角度计算子单元,用于根据所述第一探测器通道按如下公式来计算所述X射线传播路径偏离X射线束中心的角度:
θ=β-n1*Δβ,
其中,n1为所述第一探测器通道,Δβ为所述探测器通道的角度步距,θ为所述X射线传播路径与X射线束中心的夹角,β为所述X射线束的一半张角;
一第二通道计算子单元,用于根据所述夹角按如下公式来计算所述第二探测器通道:
n2=(β+θ)/Δβ,
其中,n2为所述第二探测器通道,且所述第二探测器通道与所述第一探测器通道关于所述探测器的中心呈镜像对称;
一第二角度计算子单元,用于根据所述夹角按如下公式来计算所述第二投影角度:
α2=α1±pi-2θ,
其中,α2为所述第二投影角度,α1为所述第一投影角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110145210.6A CN102805635B (zh) | 2011-05-31 | 2011-05-31 | 一种优化扫描流程方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110145210.6A CN102805635B (zh) | 2011-05-31 | 2011-05-31 | 一种优化扫描流程方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102805635A true CN102805635A (zh) | 2012-12-05 |
CN102805635B CN102805635B (zh) | 2015-04-08 |
Family
ID=47229631
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110145210.6A Active CN102805635B (zh) | 2011-05-31 | 2011-05-31 | 一种优化扫描流程方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102805635B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107481226A (zh) * | 2017-07-27 | 2017-12-15 | 沈阳东软医疗***有限公司 | 去除异常扫描数据的方法、装置及pet*** |
CN115381471A (zh) * | 2022-10-26 | 2022-11-25 | 南方医科大学南方医院 | 一种基于运动检测的ct扫描辅助***及方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040136490A1 (en) * | 2002-07-23 | 2004-07-15 | Edic Peter Michael | Method and apparatus for correcting motion in image reconstruction |
CN1565001A (zh) * | 2002-07-23 | 2005-01-12 | 通用电气公司 | 从投影数据得出运动信息的方法和设备 |
US20050175140A1 (en) * | 2004-02-05 | 2005-08-11 | Canon Kabushiki Kaisha | X-ray image processing apparatus and x-ray image processing method |
JP2006000376A (ja) * | 2004-06-17 | 2006-01-05 | Konica Minolta Medical & Graphic Inc | 医用画像撮影システム、情報処理装置 |
CN1915170A (zh) * | 2005-07-18 | 2007-02-21 | 西门子公司 | 对生物的运动身体部位产生图像的方法和放射线诊断装置 |
CN1918605A (zh) * | 2004-02-13 | 2007-02-21 | 皇家飞利浦电子股份有限公司 | 运动假影补偿 |
CN102805634A (zh) * | 2011-05-31 | 2012-12-05 | 上海西门子医疗器械有限公司 | 一种判断待检对象发生运动的方法和装置 |
CN102894989A (zh) * | 2011-07-28 | 2013-01-30 | 上海西门子医疗器械有限公司 | 一种减少x射线扫描剂量的方法和装置 |
-
2011
- 2011-05-31 CN CN201110145210.6A patent/CN102805635B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040136490A1 (en) * | 2002-07-23 | 2004-07-15 | Edic Peter Michael | Method and apparatus for correcting motion in image reconstruction |
CN1565001A (zh) * | 2002-07-23 | 2005-01-12 | 通用电气公司 | 从投影数据得出运动信息的方法和设备 |
US20050175140A1 (en) * | 2004-02-05 | 2005-08-11 | Canon Kabushiki Kaisha | X-ray image processing apparatus and x-ray image processing method |
CN1918605A (zh) * | 2004-02-13 | 2007-02-21 | 皇家飞利浦电子股份有限公司 | 运动假影补偿 |
JP2006000376A (ja) * | 2004-06-17 | 2006-01-05 | Konica Minolta Medical & Graphic Inc | 医用画像撮影システム、情報処理装置 |
CN1915170A (zh) * | 2005-07-18 | 2007-02-21 | 西门子公司 | 对生物的运动身体部位产生图像的方法和放射线诊断装置 |
CN102805634A (zh) * | 2011-05-31 | 2012-12-05 | 上海西门子医疗器械有限公司 | 一种判断待检对象发生运动的方法和装置 |
CN102894989A (zh) * | 2011-07-28 | 2013-01-30 | 上海西门子医疗器械有限公司 | 一种减少x射线扫描剂量的方法和装置 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107481226A (zh) * | 2017-07-27 | 2017-12-15 | 沈阳东软医疗***有限公司 | 去除异常扫描数据的方法、装置及pet*** |
CN115381471A (zh) * | 2022-10-26 | 2022-11-25 | 南方医科大学南方医院 | 一种基于运动检测的ct扫描辅助***及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102805635B (zh) | 2015-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100545745C (zh) | 放射线图像摄影装置及其控制方法 | |
CN102884550B (zh) | 用于动态灌注ct的图像数据配准 | |
CN102754128B (zh) | X射线图像中的标志识别和处理 | |
CN101953691B (zh) | X射线诊断装置以及x射线照射区域控制方法 | |
CN102711617B (zh) | 扫描规划视场调整器、确定器和/或质量评估器 | |
US6793390B2 (en) | Method for automatic arrangement determination of partial radiation images for reconstructing a stitched full image | |
CN102753097B (zh) | 在存在移动时用于在3d空间中的x射线标志定位的方法 | |
JP3999176B2 (ja) | X線ct装置、情報処理方法ならびに記憶媒体、プログラム | |
CN105377141B (zh) | X射线ct装置 | |
CN102598059A (zh) | 利用偏移探测器几何结构生成计算断层摄影图像的方法和设备 | |
CN103892862B (zh) | 用于侦查图像的自动扫描定位装置 | |
CN102177430A (zh) | 使用经移位的几何结构改进ct图像采集的方法和设备 | |
US20080019474A1 (en) | X-ray ct scanner and data processing method of x-ray ct scanner | |
CN102376084B (zh) | 使用各向异性噪声模型对ct图像的迭代图像滤波 | |
US8279999B2 (en) | X-ray CT apparatus and a method of controlling the X-ray CT apparatus | |
CN106232009A (zh) | 断层扫描设备和由断层扫描设备重构断层扫描图像的方法 | |
US8948337B2 (en) | Computed tomography image reconstruction | |
US8953739B2 (en) | Image processing apparatus and X-ray CT system | |
CN102805634B (zh) | 一种判断待检对象发生运动的方法和装置 | |
CN102805635A (zh) | 一种优化扫描流程方法和装置 | |
JP4982091B2 (ja) | 放射線断層撮影装置 | |
CN103501703B (zh) | X射线ct装置 | |
US10217248B2 (en) | Method for removing streak from detector cell with performance difference | |
CN102894989B (zh) | 一种减少x射线扫描剂量的方法和装置 | |
CN104039231A (zh) | X射线ct装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |