CN101128153A - X射线ct设备 - Google Patents

Abstract

一种X射线CT设备，包括：X射线源，向受检者照射X射线；X射线探测器，以使受检者位于X射线探测器和X射线源之间的形式与X射线源相对放置，并且用于探测透射过受检者的X射线，作为投影数据；旋转装置，使X射线源和X射线探测器旋转；控制装置，通过从转动装置对X射线管和X射线探测器的转动所获得的多个角度方向收集投影数据，对收集到的投影数据进行重建计算来产生受检者的层析图，并控制X射线源和旋转装置；以及显示装置，显示产生的层析图。所述X射线CT设备进一步包括：投影数据分析装置，从投影数据重建受检者的分析用成像部位的层析图，并通过对重建层析图像的重新投影产生控制曲线；管电流控制装置，根据产生的控制曲线控制馈入X射线管的电流值。

Description

X射线CT设备

技术领域

本发明涉及对X射线CT设备的改进，使控制馈入X射线源(在下文中称为X射线管)的X射线源电流(在下文中称为X射线管电流)成为可能，此电流依据受检者身体轴线方向的位置确定该受检者的成像部位中的X射线照射密度。

本申请要求基于日本专利法下的日本专利申请No.2005-051497的巴黎公约优先权，且通过引用而享有日本专利申请No.2005-051497的权益。

背景技术

X射线CT设备得到的图像质量由诸如空间分辨率和图像所含的噪声之类的因素决定。在上述因素中，图像噪声主要取决于X射线透射强度。X射线透射强度由X射线的照射强度和跨过受检者成像部位横截面的X射线透射长度决定。X射线的照射强度由馈入X射线管的电流决定。若把人体的构成假设为椭圆，则依据透射强度在前面和后面方向(在下文中可称为前后方向)与侧面到侧面方向(在下文中可称为右到左方向)的差来确定受检者成像部位的横截面的吸收系数。透射X射线强度的减弱增加噪声与透射X射线强度信号的比值。

例如，专利文献1中披露了改进透射强度减弱的技术。在专利文件1中，透射X射线强度减弱的改进是通过以下步骤来实现的。首先，在对用于诊断的X射线层析图像的测量(亦称作真实扫描)之前所获得的扫描图图像数据进行分析，并利用分析的扫描图图像数据创建受检者的类似三维透射长度模型。随后，基于先前创建的类似三维透射长度模型和在真实扫描的扫描条件，创建依赖于受检者成像部位的X射线管电流的变化模式。

专利文件1：JP-A-2002-263097

发明内容

本发明要解决的问题：

然而，现有技术仅限于根据创建的扫描图来设置X射线管电流控制值，并不考虑扫描图中X射线吸收量对于受检者的不同组织而不同这一事实，例如骨骼对X射线的衰减大，而软组织小。例如，在如肩部这一多骨骼的部位，透射X射线强度不足。在透射X射线强度不足的情况下，噪声量相对于透射X射线强度信号增加。即，即使已经采用现有技术，由于透射X射线强度不足导致噪声量增加，因此在受检者具有多骨骼的部位层析图像质量下降这个问题仍然没有解决。

本发明的目的是提供一种X射线CT设备，使在受检者透射X射线强度变化大的部位能够得到高质量的层析图像。

解决问题的措施：

根据本发明的X射线CT设备包括：X射线源，向受检者照射X射线；X射线探测器，以使受检者位于X射线源和X射线探测器之间的方式与X射线源相对放置，并探测通过受检者的透射X射线作为投影数据；旋转装置，旋转X射线源和X射线探测器；控制装置，收集通过旋转装置旋转X射线源和X射线探测器所获得的多个角度方向的投影数据，对这些收集的投影数据进行重建计算来产生受检者的层析图像，并控制X射线源和旋转装置；显示装置，显示产生的层析图像，其特征在于：上述X射线CT设备还包括投影数据分析装置，从投影数据重建受检者成像部位的分析用层析图像，并通过将重建的层析图像重新投影产生控制曲线(control profile)；以及管电流控制装置，根据产生的控制曲线控制馈入X射线管的电流值。

本发明的优点

根据本发明，甚至在受检者透射X射线强度变化大的部位也可以得到高质量的层析图像。

附图说明

图1是本发明所应用的X射线CT设备的示意组成图。

图2是解释投影数据分析设备的第一实施例的方框图，该投影数据分析设备在本发明中应用于管电流控制单元。

图3是解释本发明中应用的管电流控制流程的流程图。

图4(A)是分析用重建层析图像，它是图2示出的根据本发明的投影数据分析设备基于具有相对较少骨骼的受检者成像部位的投影数据而重建的。

图4(B)是分析用重建层析图像在阈值处理之后的图像，以及解释通过在X方向重新投影该图像得到的X方向的骨骼衰减曲线B的简图。

图5(A)是分析用重建层析图像，它是图2示出的根据本发明投影数据分析设备基于具有相对较多骨骼的受检者成像部位的投影数据而重建的。

图5(B)是分析用重建层析图像在阈值处理之后的图像，以及解释通过在X方向重新投影该图像得到的X方向的骨骼衰减曲线B的简图。

图6是解释根据本发明的X射线CT设备开始实际扫描时，如图2所示的根据本发明的投影数据分析设备的各个处理过程的时序图。

图7是模型简图，示出在X方向和y方向的重新投影结果，它是沿着受检者身体轴线方向针对三种组织：软组织、骨骼和肺部，根据由投影数据分析设备得到的分析用受检者重建图像以及受检者相应部位的扫描图所得到的。

具体实施方式

下文将参考附图解释本发明X射线CT设备的实施例。根据本实施例的X射线CT设备1在根据X射线管位置(θ，z)控制X射线管电流时，测量受检者的投影数据。进而，在用于解释本发明实施例的所有附图中，具有相同功能的附图被标以相同的参考号，它的重复解释将被略去。

图1示出了根据本实施例的X射线CT设备的组成。此X射线CT设备1由扫描器10、连接于扫描器10的主计算机20、连接于主机算计10的操作单元24和显示单元30组成。

首先解释扫描器10组成单元。

X射线管11向受检者照射X射线。X射线探测器12与X射线管11相对放置，并且探测透射受检者的X射线。数据测量设备(DAS：数据采集系统)13对由X射线探测器12探测到的透射X射线执行预定的处理，并计算投影数据。X射线高压设备14是馈给X射线管11的电源，设计成能够被操作单元24设置电压、电流和供电时间(相当于X射线照射时间)，这将在下文解释。电源提供的电压称为管电压，电源提供的电流称为管电流。数据传输设备15在旋转系统和静止系统之间传递数据，并且包括滑动环和电刷，或者旋转变压器。扫描控制设备16控制旋转板(扫描器)的旋转量，X射线管11和X射线探测器12附于该旋转板上。投影数据分析设备17根据据由数据测量设备13计算的投影数据数重建受检者的层析图像。管电流控制单元18控制X射线管11的管电流。这些组成单元被分别置于可绕受检者旋转的旋转体上。

现在，主机算机20对X射线CT设备1执行总体控制，并且包括一下组成单元。控制设备21由例如CPU和存储器组成，用于对X射线CT设备1进行总体控制。数据接收设备22，从扫描器10中的透射设备15接收测量数据。图像处理设备23根据通过数据接收设备22接收到的测量数据进行图像重建计算，产生上述层析图像。操作单元24是例如轨迹球、鼠标和键盘，用于针对管电流控制模式的开/关(ON/OFF)来提供输入指令。

接着，显示设备30具有显示产生的层析图像的功能。尽管在图示中被忽略，X射线CT设备1具有病人台(patient table)，在扫描时放置受检者。在根据本实施例的X射线CT设备1中，管电流控制模式开/关可根据操作单元24的输入指令来进行。在管电流控制模式关时，实际扫描在恒定X射线照射强度下进行。X射线探测器12测量的透射X射线强度数据被DAS13转换成数字数据后，从旋转系统中的数据传输设备15传递到静止系统。在静止系统，投影数据被数据接收设备22中的接收单元获取，发送到图像处理设备23，并且在被重建为层析图像后，图像被显示在显示设备30上来进行图像判读。

将在下文参考图2来描述根据第一实施例的投影数据分析设备17的组成。

投影数据分析设备17为要流入X射线管11的管电流确定控制值。所述数据分析设备17包括分析用图像重建单元17a、重建图像分析单元17b、管电流控制值计算单元17c。当投影数据输入时，分析用图像重建单元17a开始对投影数据进行重建处理，并针对每一个预定的观察时间间隔来计算层析图像。重建图像分析单元17b分析重建的图像，并计算例如横截面的骨骼和软组织各自的投影最大值，所述投影最大值是通过将这些和包括图像的轴向和垂轴宽度的已转换透射长度进行加和而获得的。管电流控制值计算单元17c根据获得的已转换透射长度和来计算最优管电流控制值，并且将包括新管电流值的成像条件输入到X射线高压设备14。

进而，投影数据分析设备17可配备有部位权重确定单元17d。

更具体的，如图2的虚线区域所示，可以加入部位权重确定单元17d。关于器官部位的信息通过对重新投影值的绝对量、比例或者变化进行估计来获得。部位权重确定单元17d利用所获得的器官部位信息针对每一个器官部位修改电流控制值计算单元17c获得的电流控制值，并向X射线高电压设备14输出新的电流控制值。例如，尽管针对每个器官部位来预先确定权重系数，但是通过将计算的电流值乘以权重来得到新的电流控制值。例如，当受检者为女性时，对于骨盆腔(腹下部)的权重被设置为一个小值，照射量被控制到比从透射长度获得的照射量低。

在X射线CT设备1中的管电流控制单元18，管电流通过以下步骤来得以控制。首先，在实际扫描中，在任意时间从测量的投影数据(测量数据)来重建分析用层析图像。在任意时间重建的分析用层析图像被分析，而且基于该分析来确定新管电流值。新确定的管电流值被直接馈回X射线管11。在这些处理流程序列中，在对扫描图成像和确定成像范围与X射线条件初始值后，过程进入到实际扫描，以开始每个视角(view)的投影数据测量(视角数据测量：view data measurement)。在实际扫描中，在任意时间更新包括上述管电流值的X射线状态。

下文将参照图3详细解释管电流控制处理流程。

首先，在步骤S1，用X射线CT设备1进行扫描图成像。在步骤S2，基于步骤S1拍摄的扫描图来设置成像范围。在步骤S3，设置包括管电流的X射线成像条件。在此设置的管电流值是成像条件的初始值。在步骤S4，根据步骤S2和S3设置的成像范围和成像条件进行视角数据测量。在步骤S5，通过在投影数据分析设备17的输入级配备的数据压缩设备，进行投影数据的压缩处理。

在步骤S6，将在步骤S5压缩的投影数据输入到投影数据分析设备17的分Q析用图像重建单元17a。分析用图像重建单元17a为分析用的图像进行重建处理。在步骤S7，基于步骤S6创建的分析用图像，重建图像分析单元17b分析重建的图像，并根据一个预定参考材料(例如水)，计算出分析用图像的转换透射长度。在步骤S8，管电流控制值计算单元17c根据获得的已转换透射长度并且利用例如透射长度-控制值表来计算最优管电流值。在步骤S9，根据与初始状态的偏差来判断是否基于管电流控制值计算单元17c计算的最优管电流值来修改X射线高压设备14的管电流值。当管电流不被修改，处理回到步骤S4，根据当时的管电流值继续进行视角数据测量。当管电流被修改时，处理进入步骤S10，包括管电流的X射线条件被重置，新重置的管电流值被输入到X射线高压设备14。接着处理回到步骤S4，在新的电流控制值下，根据X射线条件进行观察数据测量。

在步骤S11，根据在步骤S4测量的投影数据，图像处理设备23进行图像重建处理。步骤S11和S5～10可以并行进行。通过此并行处理，可以提高测量和数据处理序列的速度。在步骤S4获得的投影数据通过旋转系统中的数据传输设备15输出给主计算机20。在步骤S11，主计算机20中的数据传输设备22接收投影数据，并将其输出给图像处理设备23。图像处理设备23根据投影数据进行图像重建处理，并输出层析图像给显示设备30。在步骤S12，判断测量是否完成，并且当结果是“是”时，测量结束。当结果是“不”时，处理回到步骤S4。

现在，重建图像分析单元17b将在下文进行更详细解释。重建图像分析单元17b主要用来对分析用重建图像进行重新投影处理，并且能够在重新投影方向计算转换透射长度信息。进而，根据本实施例的重建图像分析单元17b能够在重新投影时进行阈值处理。例如，能够计算在图4和图5中示出的重新投影数据B。

图4(A)示出了层析图40的例子，包含骨骼部位(B)41、42、43、44、46和49，软组织(T)45、46和47。在图4(A)中，用预定CT值对层析图40进行阈值处理，并且例如骨骼和软组织的数个区域被抽出。接着，从X方向抽出层析图40区域的X射线透射曲线被计算出来作为4A。在计算出的曲线中，重新投影值的最大值Bmax被确定。依照同样办法，在除X方向之外的重新投影值的最大值Bmax也被确定，并且受检者周围所有的值都加在一起从而获得控制曲线，且被反映给X射线管电流。

图5是一个例子，其中得到了在与图4不同的另一部位的已转换透射长度。类似于图4，示出了一个层析图50例子，包含骨骼部位(B)51，52，53，54，56和59，和软组织(T)55，56和57。在此例中，X射线透射曲线被计算作为5A。当将图4中的4A和图5中的5A作比较时，可以发现重新投影值的最大值大体根据受检者的测量位置变化。

管电流控制值计算单元17c将获得的控制曲线转换为已转换透射长度，并从已转换透射长度计算最优管电流控制值。当计算管电流控制值时，受检者被近似为一个充满了均匀材料(例如，水)的椭圆，转换通过利用具有椭圆近似结构的材料(例如，水)的等效透射长度作为参考来进行。例如，由于∑F(j)是整个CT值在重新投影方向(例如X方向)的积分值，水的等效透射长度Dw[m]由以下公式(1)确定。其中，空气的CT值假设为0，水为1000。

Dw＝∑F(j)*p/N*1000         ...(1)

在此，p是像素大小[mm]，N是重新投影束数量，它覆盖受检者(对应于当从重新投影方向看去时受检者的宽度)。管电流的电流控制值Q预先作为Dw的函数(或者以表格的形式)给出，并且根据公式(2)确定。由于使用分析受检者重建图像而不是投影数据，在受检者横截面的估计精度增加。

Q＝f(Dw)                   ...(2)

当利用阈值进行重新投影处理进行后，软组织和骨髂的透射长度可以根据公式(3)和公式(4)计算。其中Nt和Nb是分别包括软组织和骨骼的重新投影束的数量，骨骼的CT值假定为2000，软组织的CT值假定1040。

Db＝∑F(j)*p/N*2000        ...(3)

Dt＝∑F(j)*p/N*1040        ...(4)

在根据公式(5)分别对组织和骨骼的透射长度进行加权时，在确定了水的透射长度之后，管电流控制值根据公式(2)确定。

尽管权重系数Wt和Wb可以通常假设为1.0，然而，由于当管电压低时，骨骼的影响增加。通过增大骨骼权重Wb能够将管电流设置的相对高。

进而，在婴儿情况下，由于骨骼的影响小，优选的将管电流到尽可能小，以抑制X射线曝光。在此情况下，骨骼权重Wb可被设置为低。

进而，Wt和Wb可用于给每个成像器官部位的计算权重，

Dw＝Wt*Dt+Wb*Dt         ...(5)

进而，在不使用权重系数的方法中，事先确定软组织的透射长度与骨骼的透射长度之间的关系，以及最优管电流值。在此情况下，对于最优管电流的电流控制值Q表示为公式(6)。

Q＝f(Dt，Db)            ...(6)

在此情况下，尽管可以分析医学数据，可以通过统计确定一个二维表。

在任何情况下，由于能够计算骨骼的透射长度被计入的受检者的横截面，，进一步提高了受检者横截面的测量精度。

重新投影处理能够在多个方向上基于一幅分析用重建图像而进行，并且当所述处理在多于一个方向进行时，在X方向(时钟3点方向)和Y方向(时钟0点方向)，能够估计在时钟6点方向和9点方向的超前相位的透射长度。当以正弦波或者任意函数控制管电流时，如果给出其最大值和最小值就是足够的，可以进行X方向和Y方向的交替计算来减少计算量。

现在描述X射线CT设备1的操作。当根据本实施例的X射线CT设备1中X射线探测器12包括例如1024通道的探测器单元时，每转动一周获得1024次视角数据。进一步，利用来自DAS 13的投影数据，在多片(multi-slice)探测器中间两行的均值被建立，并通过在投影数据分析设备17的输入级提供的数据压缩设备输入到投影数据分析设备17。进一步，在视角方向和通道方向两个方向增加8个采样数据，128通道的视角数据在转动一周输入128次。若假定初始数据的采样间隔是0.5mm，能够重建具有分辨率4mm的分析用图像。进一步，若假设分析用图像矩阵为128像素，最大视场(感兴趣区域)在4mm分辨率将是512mm。若与读片用重建图像相比，由于重建矩阵是1/16，并且视角数目是1/8，其计算量将是约1/128。在此情况下，通过模糊校正滤波器进行带有滤波的校正处理，所述模糊校正滤波器是通过进行256点傅立叶变换而预先存储的。进一步，在应用了模糊校正滤波之后进行反投影滤波。当完成重建所需的视角数据反向投影处理后，重建分析用层析图像。另外，尽管应当与水的数据一起保持空气数据和位移数据，由于例如空气数据和位移数据的量对于一个视角足够，所以对应的数据所需的小存储容量是足够的。

在多片CT中，当需要排列的行的数目增加时，需要使用例如Feldkamp方法的三维反投影计算，然而，根据本实施例，即使如果只有特定的行以其目前利用的现有二维重建方法重建，结果也显示出足够的实用性。

进一步，如果除了只对接近中心行处图像进行重建之外，还对处在边缘行的图像一起重建，并且在多个点获得身体轴线方向的已转换透射长度，则用于评估的横截面的已转换透射长度的估计精度能够提高。

在实际扫描之前，输入X射线条件的初始值X0和控制命令的时间间隔Δt。X0可以是从扫描图的分析结果得到的最优值或者操作者经验得到的值。

进一步，尽管控制命令的时间间隔Δt可任意确定，然而，如果管电流值经常更新且X射线曝光没有优化，命令的更新是无意义的。因此，时间间隔是根据例如控制系统的反应时间来确定的，并且，当然，时间间隔可以被定义为如视角的数目和角度间隔。例如，当假设视角数目在转动一周是128次，管电流值每32个视角输出一次，指令可以每转一周更新4次。

图6示出了扫描中各个过程的时序图。当在t0扫描开始，X射线开始曝光，与之同步数据收集操作开始。此时管电流初始值被假设为X0。全部测量数据被送入数据传输设备15并且传输到静止系统。在主算计机20获取的数据在图像处理设备23中被重建为读片用层析图。

在另一方面，投影数据分析设备17进行分析用图像重建、重建图像分析和最优管电流控制值计算过程，来确定管电流控制值。分析用图像重建、重建图像分析和最优管电流控制值计算过程以流水线方式分别进行。通过最优管电流控制值计算得到的新的管电流控制值送入X射线高压设备14，并且反馈在t6首次起效。反馈起效的时刻可以被用作一个实际成像起始点。

此处获得的管电流值是基于由t1-t2时刻之间获得的投影数据重建的分析用图像所设置的，并且包括对应的阶段(角度)信息。在此处，X高压设备14在对应阶段t6改变实际控制值。

在如现有技术一样通过直接处理投影数据确定管电流控制值时，当例如前置放大器增益、log转换增益和焦点探测器距离变化时，管电流的参数必须复查。然而，在本实施例中，由于管电流值是基于表示绝对值的CT值而设置的，就可以得到一个稳定的控制。

图7示出了分别对于软组织、骨骼和肺部三种组织的在X和Y方向重新投影值。水平轴表示身体轴线方向，对应于在背景示出的扫描图的位置。这里，当表示软组织的重新投影结果时，只有肺部存在位置在X和Y两个方向示出了高重新投影值。

在另一方面，在接近肩膀处骨骼的重新投影值增加。以此方式，可以理解结果分别示出了成像部位的特点。从而，根据本实施例，各个部位可以通过分别利用重新投影值进行判断。

在本实施例中，尽管曝光剂量对于不同部位可以自动控制，通过由操作员利用例如扫描图将一个特别需要的部位设置为低照射剂量，可以利用这个信息进行照射剂量的控制。然而，若为了根据层析图的测量部位设置管电流值，从用户友好度来看，这些是进一步需要的。

进而，在本实施例中，描述了透射长度转换成水对X射线的吸收系数的例子。然而，作为X射线吸收系数的参数，材料不局限于水，任何材料例如骨骼和软组织都可以选择。因为如果受检者能够转换成为预定材料的透射长度的形式，那就是足够的。

Claims (10)

1.一种X射线CT设备，包括：

X射线源(11)，向受检者照射X射线，

X射线探测器(12)，以将受检者放置在X射线探测器(12)和X射线源(11)两者之间的方式与X射线源(11)相对设置，并且用于探测通过受检者的透射X射线，作为投影数据，

旋转装置(19)，旋转X射线源(11)和X射线探测器(12)，

控制装置(20)，从旋转装置(19)对X射线源(11)和X射线探测器(12)的旋转得到的多个角度方向收集投影数据，对收集的投影数据进行重建计算来产生受检者的层析图，并且控制X射线源(11)和旋转装置(19)，以及

显示装置(30)，显示产生的层析图，

其特征在于，所述X射线CT设备进一步包括：

投影数据分析装置(17)，从投影数据重建受检者的分析用成像部位的层析图，并且通过对重建层析图进行重新投影来产生控制曲线，以及

管电流控制装置(18)，根据产生的控制曲线控制要馈入X射线管的电流值。

2.根据权利要求1的X射线CT设备，其特征在于，当从预定方向照射X射线，所述投影数据分析装置(17)通过对由具有不同X射线吸收量并存在于透射方向的部位组成的X射线吸收量曲线进行重新投影，来产生控制曲线。

3.根据权利要求1的X射线CT设备，其特征在于，投影数据分析装置(17)确定围绕受检者整个周边重新投影的X射线吸收量曲线，并且通过添加整个周边的X射线吸收量曲线来产生控制曲线。

4.根据权利要求1到3之一的X射线CT设备，其特征在于，投影数据分析装置(17)包括：

分析用图像重建装置(17a)，从投影数据的一部分来重建分析用层析图，

重建图像分析装置(17b)，针对重建的分析用层析图，计算具有预定X射线吸收的参考材料的X射线透射长度，作为已转换透射长度，以及

管电流控制值计算装置(17c)，根据计算的已转换透射长度计算要馈入X射线管的电流值。

5.根据权利要求4的X射线CT设备，其特征在于，投影数据分析装置(17)进一步包括部位权重确定装置(17d)，利用与每个器官部位相对应的器官部位信息，修改由管电流控制值计算装置(17c)计算的管电流控制值。

6.根据权利要求4的X射线CT设备，其特征在于，投影数据分析装置(17)获取受检者的扫描图，根据扫描图为X射线成像设置包括成像范围和管电流值的初始条件，根据确定的成像范围和成像条件测量投影数据，利用测量的投影数据进行分析用图像的重建处理产生进行了重建处理的分析用图像的预定方向上的控制曲线，根据产生的控制曲线确定已转换透射长度，并且由确定的已转换透射长度计算管电流值。

7.根据权利要求6的X射线CT设备，其特征在于，投影数据分析装置(17)利用参考材料从已转换透射长度计算管电流值。

8.根据权利要求6的X射线CT设备，其特征在于，投影数据分析装置(17)向骨骼透射长度和软组织透射长度提供预定的权重，并计算已转换透射长度。

9.根据权利要求8的X射线CT设备，其特征在于，在基于分析用层析图指定受检者器官部位之后，根据指定的器官部位来执行预定权重处理。

10.根据权利要求6到9之一的X射线CT设备，其特征在于，要分析的已转换透射长度是通过包括水等效透射长度、骨骼透射长度和软组值透射长度的透射长度转换的。

Images