阅读说明：本技术 一种复合能谱ct成像方法 (Composite energy spectrum CT imaging method ) 是由 曾凯 徐丹 于 2019-10-08 设计创作，主要内容包括：本发明提出一种复合能谱CT成像方法,该方法基于传统CT系统,扫描过程中通过高压发生器按照指定频率切换低电压和高电压,获得复合扫描数据；对低电压扫描数据和高电压扫描数据分别进行图像重建,得到低电压图像和高电压图像；对重建后的低电压图像和高电压图像进行图像域的物质分解得到基物质图像。本发明基于传统CT系统实现,电压切换的速度可以根据现有硬件的条件来调整,不需要对硬件进行任何改动,因此实现成本极低。由于切换速度的降低,电压能保持和单能扫描一样的稳定性,因此不需要进行特殊的校准。相对于GE的高速切换方式,本发明主要有成本低,易实现,能量分离大,剂量低的优点。(The invention provides a composite energy spectrum CT imaging method, which is based on a traditional CT system, wherein in the scanning process, a high-voltage generator switches low voltage and high voltage according to specified frequency to obtain composite scanning data; respectively carrying out image reconstruction on the low-voltage scanning data and the high-voltage scanning data to obtain a low-voltage image and a high-voltage image; and performing image domain material decomposition on the reconstructed low-voltage image and the reconstructed high-voltage image to obtain a base material image. The invention is realized based on the traditional CT system, the voltage switching speed can be adjusted according to the conditions of the existing hardware, and the hardware does not need to be changed, so the realization cost is extremely low. Due to the reduced switching speed, the voltage can maintain the same stability as a single energy scan, and therefore no special calibration is required. Compared with a GE high-speed switching mode, the method mainly has the advantages of low cost, easiness in implementation, large energy separation and low dosage.)

一种复合能谱CT成像方法

技术领域

本发明涉及能谱成像技术领域，尤其是一种复合能谱CT成像方法。

背景技术

能谱成像技术对医学影像诊断有着重要的意义，他能够分离物质不同能量的信息，显著抑制射线硬化伪影，给临床诊断带来更多的依据。但是要达到能谱成像必须要采用比传统CT更加先进的硬件系统，比如西门子的双射线源双探测器CT系统，通用电气的高速切换CT扫描系统，飞利浦的双层探测器CT系统。这些系统都是依赖他们自己专有的硬件技术来扫描获得不同能量的物质信息。因此传统的CT系统难以在实际临床应用中进行能谱成像。

能谱CT自从西门子在2006年推出双源双探测器的解决方案开始，已经在科研和临床上得到了广泛应用。传统CT成像提供的是被扫描物体在某个扫描kV下对X射线的有效吸收系数图像，与物体的大小、所使用的射线过滤器等等都相关。即使在同一扫描条件下，在不同CT扫描系统上的CT值也有差异。因此在传统CT上，定量分析相对比较困难。能谱CT通过在不同能量下对物体进行扫描，能感知物体对不同能量X射线的衰减情况，因此能分辨出被扫描物体中物质的组成，而且得到的图像受扫描条件的影响小，能更准确地提供定量的分析。

能谱CT是通过测量被扫描物体(病人)对不同能量的x射线的吸收情况来达到物质分解的目的。为了获得更好的图像质量，一般来说，这在两个或多个能量下的测量要求有好的同时性，以及大的能量分离。目前来说市面上已有产品中有表1所示的几种实现方式：一是通用电气公司采用的能谱快速切换，二是西门子公司的双源双探测器，第三是飞利浦公司的双层探测器。另外还有光子计数CT正在研发中，市场还没有能用于临床的产品。

表1现有能谱CT技术

目前的技术水平来说，光子计数CT的同时性和能量分离是最好的，但是其造价也是相当昂贵。而且光子计数探测器在高剂量下的稳定性问题一直没有得到很好的解决，因此在商业化的道路上走得比较艰难。快速能谱切换虽然有很好的同时性以及低廉的成本，但是其切换速度也是有极限的。当切换速度越来越快时，由于球管电压存在一定的上升和下降时间，切换的能谱之间的差异将越来越小，导致重建图像的效果变差。双源双探测器的成本也居高不下，因此一直只在西门子的最高端CT产品中才有能谱扫描技术。双层探测器虽然有很好的同时性，但是它的能量分辨能力是所有技术中最差的，而且探测器本身的造价也较高。

发明内容

发明目的：为了降低CT成像对于能谱CT系统的硬件要求，本发明提出一种复合能谱CT成像方法，该方法基于传统CT系统而提出的，主要是通过扫描控制、适合的图像重建算法和能谱分解技术来得到能谱图像，在硬件上不需要进行任何更改。

技术方案：为实现上述目的，本发明提出的技术方案为：

一种复合能谱CT成像方法，包括步骤：

(1)在扫描过程中通过高压发生器按照指定频率切换低电压和高电压，获得复合扫描数据，具体包括在一个扫描周期内进行a、b、c所述的任意一种切换方式：

a.射线源每旋转一圈切换一次高低电压；

b.射线源每旋转N圈切换一次高低电压，0.5≤N<1，且保证一个扫描周期内至少分别获得半圈高电压扫描数据和低电压扫描数据，前一次高压/低压扫描和后一次低压/高压扫描之间关闭或不关闭射线；

c.在每一圈扫描过程中切换多次，相邻两圈扫描过程中电压切换方式互补，使获得的每一圈扫描数据中高电压扫描数据区域和低电压扫描数据区域交替分布，而相邻两圈扫描数据中，前一圈的高电压扫描数据区域与后一圈的低电压扫描数据区域在位置上一一对应；

(2)对低电压扫描数据和高电压扫描数据分别进行图像重建，得到低电压图像和高电压图像；

(3)对重建后的低电压图像和高电压图像进行图像域的物质分解得到基物质图像。

进一步的，所述扫描获得扫描包括轴扫描和螺旋扫描。

进一步的，当采用螺旋扫描时，对于b所述的切换方式，扫描的螺距p满足：

前一次高压/低压扫描和后一次低压/高压扫描之间的射线关闭时间为

进一步的，所述图像重建的方法包括：滤波反投影法、迭代重建法和压缩感知法。

进一步的，当采用螺旋扫描方式进行扫描时，所述步骤(2)中采用去伪影的图像重建方法对低电压扫描数据和高电压扫描数据分别进行图像重建，具体步骤包括：

(51)构建训练数据：获得清晰的无伪影CT图像作为目标图像；通过对与目标图像相应的体模进行螺旋CT扫描得到螺旋扫描锥形束投影数据，对螺旋扫描锥形束投影数据通过现有的图像重建方法进行图像重建，得到初始图像；将初始图像经过数值模拟过程产生对应的模拟锥形束投影数据；对模拟锥形束投影数据采用与初始图像相同的图像重建方法重建得到二次图像；

(52)搭建卷积神经网络，在卷积神经网络的输入层增加一个输入通道，即构建好的卷积神经网络具有两个输入通道；将初始图像和二次图像分别送入两个输入通道，通过卷积神经网络提取初始图像和二次图像的灰度信息、结构信息，并根据提取出的灰度信息、结构信息进行自主学习后估计真实图像；

(53)构建关于目标图像和估计真实图像的损失函数，利用梯度下降法训练所述卷积神经网络；所述损失函数为：

其中，Img表示卷积神经网络的输出图像，Img_ture表示与Img对应的目标图像，Img_k表示图像Img中第k个像素点的表示像素值，Img_ture,k表示图像Img_ture中第k个像素点的表示像素值；

(54)对于步骤(1)中得到的低电压螺旋扫描锥形束投影数据和高电压螺旋扫描锥形束投影数据，分别采用步骤(51)所述的方法得到初始图像和二次图像；将得到的初始图像和二次图像送入训练好的卷积神经网络，得到相应的清晰无伪影的高/电压重建图像。

具体的，所述卷积神经网络包括：CNN、ResNet和Unet。

有益效果：与现有技术相比，本发明具有以下优势：

本发明在实现方式上虽然与GE一样是通过高压切换来实现双能扫描，但是本发明不需要像GE的系统一样在每次采样都进行切换。切换的速度可以根据现有硬件的条件来调整，不需要对硬件进行任何改动，因此实现成本极低。由于切换速度的降低，电压能保持和单能扫描一样的稳定性，因此不需要进行特殊的校准。曝光电流也容易修改，而且有能力实现自动电流调制。总体来说，相对于GE的高速切换方式，本发明主要有成本低，易实现，能量分离大，剂量低的优点。

附图说明

图1为本发明的流程图；

图2为采用轴扫描方式时，分别在低电压和高电压扫描完整的两圈数据示意图；

图3为采用轴扫描方式时，在同一圈内切换多次高低电压，得到的相邻两圈扫描数据的示意图；

图4为采用轴扫描方式时，在高、低电压分别扫描两个半圈数据的示意图；

图5为螺旋模式下每隔一圈切换低电压和高电压的扫描数据示意图；

图6为螺旋模式下每隔0.75圈切换低电压和高电压的扫描数据示意图；

图7为螺旋模式下在特定的方向上降低电流或完全关闭射线的扫描数据示意图；

图8为卷积神经网络的结构图。

具体实施方式

下面结合附图对本发明作更进一步的说明。

本发明所提出的能谱CT成像方法是基于传统CT系统而提出的，通过常规CT系统实现更加灵活的电压切换的扫描方式来采集数据，并利用图像重建和能谱分解技术来得到能谱图像。该方法在硬件上不需要进行任何更改。

对于能谱成像来说，最关键的就是要获得不同能谱的物质信息。本发明提出了一种利用现有CT硬件，在只改变扫描控制方式的基础上，在低电压和高电压之间进行切换，得到一组复合能量下的扫描数据，然后利用高级重建技术，实现双能扫描的方法。CT系统的一次扫描中，高压控制器在低电压和高电压中以现有硬件支持的速度反复切换，达到双能扫描的目的，扫描方式可以是轴扫，也可以是螺旋扫描。扫描过程中，高压发生器在低电压和高电压之间来回切换。切换的频率可以调整，保证在每个能量下，能有足够的数据分别独立重建低电压和高电压的图像。

扫描通常有轴扫描模式和螺旋扫描模式。

在轴扫模式下，病床不动，为了让曝光后在低电压和高电压都能有足够的数据重建图像，一般来说每个电压下面都需要有至少半圈完整的数据。这种模式下，可以设计多种切换方式。

图2至4给出了在扫描周期为2圈的情况下的3种轴扫模式下的扫描数据示意图。其中，图2为采用轴扫描方式时，分别在低电压和高电压扫描一圈完整的数据。图3为采用轴扫描方式时，扫描两圈完整的数据，其中每圈内切换数次，但是第一圈和第二圈的切换方式相反，即第一圈内的高电压区域与第二圈的低电压区域互相重合，第一圈的低电压区域与第二圈的高电压区域互相重合。图4为采用轴扫描方式时，在高、低电压分别扫描两个半圈数据的示意图。在这些切换模式下，低电压和高电压都分别有足够的数据重建出图像。

在螺旋扫描模式下，我们可以通过控制螺距和切换间隔，来保证低电压和高电压下能分别重建出完整的图像。在螺旋扫描模式下，螺距是p(每圈的移动距离/扫描宽度)。每个扫描切换的周期可以按照1/p圈来规定，比如p＝0.5，每两圈就是一个扫描周期。对于每一个扫描周期，图5至7给出了3种螺旋扫描模式下的扫描数据示意图。其中，图5是螺旋模式下射线源每旋转一圈切换一次电压；图6是螺旋模式下射线源每旋转0.75圈切换一次电压；图7是螺旋模式下射线源以高电压扫描0.75圈后切换低电压再扫描0.75圈，但是前一次扫描和后一次扫描之间引入时间间隔(约为0.25圈)，在这个时间间隔内降低电流或者完全关闭射线。

如图7这种情况，射线源每旋转N圈切换一次高低电压，0.5≤N<1，只需要保证一个扫描周期内至少分别获得半圈高电压扫描数据和低电压扫描数据即可，前一次高压/低压扫描和后一次低压/高压扫描之间关闭或不关闭射线。这个N也可以取其他满足条件的值，比如如p＝0.5,N＝0.75圈，这样情况下如果在前后两次扫描之间关闭射线的话，关闭时间

此处就是vof＝0.5圈；或者可以取p＝0.75，N＝0.75圈，这样vof＝0圈。

上述的轴扫描和螺旋双能扫描方式，低电压和高电压的曝光数据都能独立重建出整个扫描范围内的图像。在获得低电压和高电压图像后，可以直接在图像域进行物质分解。

下面通过具体实施方式来进一步说明本发明的技术方案，实施例提出的复合能谱CT成像方法流程如图1所示，包括步骤：

Step1:基于图2至7中举例说明的电压切换形式，利用现有CT系统进行扫描。电压的切换间隔可以采用系统硬件能够达到的速度。比如扫描转速是0.5秒一圈，采用0.75圈切换一次，低电压是80kVp，高电压是140kVp。

Step2:通过扫描采集数据后，对低电压和高电压数据分别进行校正和预处理，利用CT图像重建技术对两种电压的图像分别进行重建。

对于轴扫描数据，其重建算法包括但不限于：滤波反投影，迭代重建，压缩感知等。

对于螺旋扫描数据，在分别重建低电压和高电压的扫描数据的时候，由于切换的频率比较低，再加上重建误差和锥角(正比于探测器的排数)的平方成正比，所以当探测器排数增大到128甚至256排的时候就会带来很大的误差，具体体现在重建后的图像中存在较为严重的伪影。为了消除图像伪影，这里优选提供了一种改进的螺旋重建方法来减少锥形束伪影对图像质量带来的影响，具体包括以下步骤：

1)构建训练数据：获得清晰的无伪影CT图像作为目标图像；通过对与目标图像相应的体模进行螺旋CT扫描得到螺旋扫描锥形束投影数据，对螺旋扫描锥形束投影数据通过现有的图像重建方法进行图像重建，得到初始图像；将初始图像经过数值模拟过程产生对应的模拟锥形束投影数据；对模拟锥形束投影数据采用与初始图像相同的图像重建方法重建得到二次图像；

2)搭建如图8所示的卷积神经网络，在卷积神经网络的输入层增加一个输入通道，即构建好的卷积神经网络具有两个输入通道；将初始图像和二次图像分别送入两个输入通道，通过卷积神经网络提取初始图像和二次图像的灰度信息、结构信息，并根据提取出的灰度信息、结构信息进行自主学习后估计真实图像；

3)构建关于目标图像和估计真实图像的损失函数，利用梯度下降法训练所述卷积神经网络；所述损失函数为：

其中，Img表示卷积神经网络的输出图像，Img_ture表示与Img对应的目标图像，Img_k表示图像Img中第k个像素点的表示像素值，Img_ture,k表示图像Img_ture中第k个像素点的表示像素值；

4)对于步骤(1)中得到的低电压螺旋扫描锥形束投影数据和高电压螺旋扫描锥形束投影数据，分别采用步骤1)所述的方法得到初始图像和二次图像；将得到的初始图像和二次图像送入训练好的卷积神经网络，得到相应的清晰无伪影的高/电压重建图像。

上述改进的螺旋重建方法中，所述的现有的图像重建方法包括但不限于滤波反投影法、迭代重建法和压缩感知法，所述卷积神经网络包括：CNN、ResNet和Unet。

Step3：得到的两个电压的重建图像可以直接在图像域进行物质分解。物质分解时，图像域的多项式的分解函数F可以通过下面的方法来求解得到：

其中，F_wat、F_iod分别表示待求解的图像阈物质分解的函数，通常用多项式表示；j表示输入的高电压图像/低电压图像的每个像素的索引值；Img_low、Img_high分别表示图像重建得到的低电压图像和高电压图像，Img_wat、Img_iod分别表示真实的水基图像和碘基图像，通常可以通过体模测量或数值仿真提前设定得到。

目前能通过高压切换进行双能扫描的商业CT系统是GE的高速切换能谱CT。该产品在硬件上需要改变高压发生器的设计，同时降低对稳压的要求来达到高速切换的目的。这种实现方式在硬件上有成本投入，系统的能谱扫描需要进行特殊的校准，而且扫描协议中不能随意修改曝光电流，导致无法进行电流调制来降低辐射剂量。

本发明在实现方式上虽然与GE一样是通过高压切换来实现双能扫描，但是本发明不需要像GE的系统一样在每次采样都进行切换。切换的速度可以根据现有硬件的条件来调整，不需要对硬件进行任何改动，因此实现成本极低。由于切换速度的降低，电压能保持和单能扫描一样的稳定性，因此不需要进行特殊的校准。曝光电流也容易修改，而且有能力实现自动电流调制。总体来说，相对于GE的高速切换方式，本发明主要有成本低，易实现，能量分离大，剂量低的优点。

以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。