一种双能CT图像重建方法、系统及计算机可读存储介质与流程

一种双能ct图像重建方法、系统及计算机可读存储介质
技术领域
1.本发明涉及医学影像领域，尤其涉及一种双能ct图像重建方法、系统及计算机可读存储介质。


背景技术：

2.第三代ct系统的主要组成部分包括tube(球管)，collimator(限束器)，detector(探测器)。x射线球管发出x光，经过限束器限制后形成一个锥形的光束。锥形光束照射在探测器上经过探测器转化为电信号并经过数据采集和转换单元转换为数字信息存储在图像处理系统中。图像处理系统经过一系列的校正算法和图像重建算法生成图像显示在显示器上，从成本和技术程度来说，主流探测器都采用了由很多个探测器模块排列在一个弧形或者多边形的面上来上形成探测器。每个探测器模块中规则排列一个探测器单元形成的矩阵。
3.在现有技术当中，能谱ct的主要方式包括1)双层探测器：双层探测器通常是在一块电路基板的两侧贴两种探测器。离x射线焦点近的一层能吸收低能x射线，而高能x射线不能吸收。在电路基板的另一层探测器能够吸收剩余的高能x射线。两层探测器可以进行合成和分解生成能谱ct图像。2)单一球管高低压切换：采用能够瞬间快速切换球管kv的高压发生器，在每个奇数采样使用一个kv值而在偶数采样使用另外一个kv值。通过对两个kv所对应的数据进行合成分解生成能谱ct图像。3)双球管采用不同的kv：两个相距90度的球管分别采用不同的kv数值同时曝光，球管对应的两个探测器同时接收数据。两个探测器的数据进行重建生成能谱ct成像。4)光子计数探测器：能够在不同的x射线能量间隔内对x射线光子进行计数生成不同能量区间间隔的采样数据，进一步重建生成能谱ct图像。5)改变x射线滤过方式来改变能谱：通过改变x射线滤过来机械式的改变能谱，不同采样时刻得到不同能谱对应的采样数据。6)两个多半圈的ct扫描方式：两个半圈扫描分别采用两个不同的kv进行扫描。
4.在上述的多种方式的能谱ct来说，在实际图像重建中一般把高低能扫描所得到的投影数据转换为康普顿散射和光电效应组合，或者转换为两种基物质如水和碘的组合或水和钙的组合，在现有计算式中，经过一系列转换和校正后得出计算式为：
[0005][0006][0007]
其中的多项式系数c
nm
和d
nm
可以经过校正方法得到。一种得到这些系数的方法为通过输入x射线源的光谱、探测器的能谱响应曲线、射线路径上的各种滤过的能谱响应到物理模型中，在扫描路径上增加各种不同长度组合的水和碘物质进行模拟，每个水和碘物质对的吸收(p
water，j
，p
iodine，j
)，j＝1,2...都可以得到，然后通过上式拟合计算得到系数cnm和dnm。但通过此分解方法存在一个问题，当扫描物体是纯水物质时，由于d
nm
拟合误差的存在，不能保证计算出的碘吸收p
iodine
为零。进行了射束硬化处理后纯水物质的p1和p2相等，c
nm
的
拟合误差必然会导致即计算出的水不等于实际水吸收。此外，在被扫描物体的密度较小接近空气密度时，由于拟合误差的存在，会导致p
water
和p
iodine
是符号相反且绝对值大小接近的数值。
[0008]
综上所述，本发明实际解决的技术问题是：现有技术在双能ct基物质分解时采用非约束的多项式拟合，由于测量和拟合偏差不能保证纯水物质时碘的分量完全为零，因而也不能保证扫描纯水时水物质本身的准确度。


技术实现要素：

[0009]
为了克服上述技术缺陷，本发明的目的在于提供一种双能ct图像重建方法、系统及计算机可读存储介质，此方法在纯水模体的情况下p1和p2相等，因而能保证p
iodine
为零，以及p
water
＝p1＝p2。
[0010]
本发明公开了一种双能ct图像重建方法，包括步骤：
[0011]
s100：基于ct成像系统分别获取管电压值为kv1和kv2时扫描空气和扫描患者时探测器接收到的光信号量iair
kv1
、iair
kv2
和i
kv1
、i
kv2
；
[0012]
s200：根据iair
kv1
、iair
kv2
和i
kv1
、i
kv2
计算出患者在管电压值kv1和kv2下的吸收值p
kv1
和p
kv2
；
[0013][0014][0015]
s300：将p
kv1
和p
kv2
校正为单色光p1和p2；
[0016]
将p
kv1
和p
kv2
校正为单色光p1和p2的计算式为：
[0017][0018][0019]
其中an和bn为经过校正得到的多项式系数。
[0020]
s400：根据p1和p2经过多项式分解后分别得到两种基物质水和碘的吸收值p
water
和p
iodine
；
[0021]
s500：根据p
water
和p
iodine
重建基物质图，并根据基物质图计算出所需的能谱图像。
[0022]
优选地，根据p1和p2经过多项式分别得到两种基物质水和碘的吸收值p
water
和p
iodine
的计算式为：
[0023][0024][0025]
优选地，计算式中的幂指数a、b以及拟合n、m的阶数以及系数fnm和gnm由拟合偏差最小原则定出。
[0026]
优选地，根据p
water
和p
iodine
重建基物质图时，通过fbp和/或fdk的近似重建算法或解析重建算法。
[0027]
优选地，根据p
water
和p
iodine
重建基物质图时，还可以通过迭代重建算法重建出基物
质图。
[0028]
有鉴于此，本发明的目的之二在于提供一种双能谱图像重建系统，包括处理器，以及与处理器连接的存储器，
[0029]
存储器上存储有双能谱图像重建程序，双能谱图像重建被处理器执行时实现如上述双能ct图像重建方法。
[0030]
有鉴于此，本发明的目的之三在于提供一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序被处理器执行时实现如上述的双能ct图像重建方法。
[0031]
采用了上述技术方案后，与现有技术相比，本发明的有益效果在于在能谱分解时增加了运约束减少了拟合偏差；在纯水模体的情况下p1和p2相等，因而保证了p
iodine
为零，以及p
water
＝p1＝p2，且能够保证p
water
≥0，在接近空气密度的物质被扫描时，不存在水和碘的吸收正负相消的情况。
附图说明
[0032]
图1为本发明一种双能ct图像重建方法的流程示意图；
[0033]
图2为本发明一种双能ct图像重建方法的ct系统形式布局示意图。
具体实施方式
[0034]
以下结合附图与具体实施例进一步阐述本发明的优点。
[0035]
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。
[0036]
在本公开使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本公开。在本公开和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
[0037]
应当理解，尽管在本公开可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本公开范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在
……
时”或“当
……
时”或“响应于确定”。
[0038]
在本发明的描述中，需要理解的是，术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
[0039]
在本发明的描述中，除非另有规定和限定，需要说明的是，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是机械连接或电连接，也可以是两个元件内部的连通，可
以是直接相连，也可以通过中间媒介间接相连，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。
[0040]
在后续的描述中，使用用于表示元件的诸如“模块”、“部件”或“单元”的后缀仅为了有利于本发明的说明，其本身并没有特定的意义。因此，“模块”与“部件”可以混合地使用。
[0041]
参阅图1所示，本实施例提供一种双能ct图像重建方法，包括步骤：
[0042]
s100：基于ct成像系统分别获取管电压值为kv1和kv2时扫描空气和扫描患者时探测器接收到的光信号量iair
kv1
、iair
kv2
和i
kv1
、i
kv2
；
[0043]
s200：根据iair
kv1
、iair
kv2
和i
kv1
、i
kv2
计算出患者在管电压值kv1和kv2下的吸收值p
kv1
和p
kv2
；
[0044]
s300：将p
kv1
和p
kv2
校正为单色光p1和p2；
[0045]
s400：根据p1和p2经过多项式分解后分别得到两种基物质水和碘的吸收值p
water
和p
iodine
；
[0046]
s500：根据p
water
和p
iodine
重建基物质图，并根据基物质图计算出所需的能谱图像。
[0047]
需要说明的是，在上述步骤s100中，iair
kv1
、iair
kv2
分别表示以管电压值kv1和kv2扫描空气时所获取的光信号量，具体如图2所示，当球管发出x射线，经过限束器形成一个锥形光束后，穿透空气照射在探测器上，探测器将接到的一个信号量，当球管的管电压值为kv1时，x射线穿过空气后探测器接收到的光信号量为iair
kv1
，当球管的管电压值为kv2时，x射线穿过空气后探测器接收到的光信号量为iair
kv2
。i
kv1
、i
kv2
分别表示以管电压值kv1和kv2扫描患者时所获取的光信号量，如图2所示，当球管发出x射线，经过限束器形成一个锥形光束后，穿透患者照射在探测器上，探测器将接到一个信号量，当球管的管电压值为kv1时，x射线穿过患者后探测器接收到的光信号量为i
kv1
，当球管的管电压值为kv2时，x射线穿过患者后探测器接收到的光信号量为i
kv2
。
[0048]
需要说明的是，在上述步骤s200中，根据iair
kv1
、iair
kv2
和i
kv1
、i
kv2
计算出患者在管电压值kv1和kv2下的吸收值p
kv1
和p
kv2
，计算式为：
[0049][0050][0051]
需要说明的是，上述步骤s300中，由于球管发出的x射线并不是单色光，而是有若干种波长组成的x射线的组合光，因此在管电压值kv1和管电压值kv2下扫描均存在射束硬化现象，因此首先对计算出的吸收p
kv1
和p
kv2
进行射束硬化校正过的吸收值p1和p2，计算式为：
[0052][0053][0054]
需要说明的是，计算式(5)和(6)中an的bn为经过校准得到的多项式系数。
[0055]
需要说明的是，上述步骤s400中根据p1和p2经过多项式分解后分别得到两种基物质水和碘的吸收值p
water
和p
iodine
的计算式为：
[0056]
[0057][0058]
需要说明的是，上述计算式(7)和(8)中的a、b可以根据拟合的准确长度来取不同的数值，f
nm
和g
nm
可以通过理论模拟或者采用实验测量进行拟合，例如，一种特殊情况下，可以取a＝1，b＝1，计算式以上述表达式能够保证p
water
≥0，在接近空气密度的物质被扫描时，不存在水和碘的吸收正负相消的情况。
[0059]
需要说明的是，计算式(7)和(8)中的幂指数a、b以及拟合n、m的阶数由拟合偏差最小的前提下幂指数a、b以及阶数n,m最小的原则定出。
[0060]
需要说明的是，上述步骤s500中，根据p
water
和p
iodine
重建基物质图时，利用基于fbp和/或fdk的近似重建算法或解析重建算法，或者通过迭代重建算法算出基物质图。
[0061]
在实现时，具体如下：
[0062]
通过理论模拟方法以及计算式(7)和(8)进行能谱分解计算式的系数计算得到f
nm
和g
nm
；
[0063]
通过双层探测器和两个多半圈的ct扫描方式中的任意一种方式得到两个能谱下的吸收p
kv1
和p
kv2
；
[0064]
通过计算式(5)和(6)进行射束硬化校正得到p1和p2；
[0065]
通过计算式(7)和(8)进行能谱分解后得到p
water
和p
iodine
；
[0066]
利用基于fbp，fdk的近似重建算法或者解析重建算法，或者迭代重建算法重建出基物质图，并根据基物质图计算出所需的能谱图像。
[0067]
应当注意的是，本发明的实施例有较佳的实施性，且并非对本发明作任何形式的限制，任何熟悉该领域的技术人员可能利用上述揭示的技术内容变更或修饰为等同的有效实施例，但凡未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何修改或等同变化及修饰，均仍属于本发明技术方案的范围内。