Slab切块图像的连接方法与流程

本发明涉及一种Slab切块图像的连接方法，属于CT图像重建技术领域。

背景技术：

心脏扫描重建是中高端CT机的必备功能，由于心脏处于跳动状态，低速CT机扫描重建的心脏图像会有严重的运动伪影，图像信息错位，不足以清晰重建出心脏冠脉血管等更精细的组织结构。

现有的心脏图像CT重建技术包括：

1、GE的革命式CT(revolution CT)单周最快扫描时间0.2s，同时宽体探测器可达16cm的z向覆盖率，可基本涵盖成人心脏的大小，一次扫描就可重建出任意心率下整个心脏“运动冻结”状态的图像，由于扫描速度快，硬件要求很高，且探测器比其他厂商宽出几倍，成本大幅提高；

2、西门子的Force开源CT，单周最快扫描时间达0.25s，由于其采用垂直双源扫描结构，所以仅需90度扫描就可重建出心脏图像，时间分辨率可到达66ms，也可实现任意心率下“运动冻结”状态的心脏图像重建，由于该CT机使用小球管，成像算法极为复杂，具有一定局限性，推广应用较为困难；

3、利用ECG心跳信号门控技术实现心脏图像重建，即根据患者的ECG心跳信号，挑选出心脏跳动最平缓的时间段进行图像重建，如图1所示，挑选出Slab1、Slab2、Slab3切块图像，将三个切块图像依序相接进行图像重建。受限于扫描时间、探测器宽度等硬件结构的影响，CT机单次扫描可重建出心脏的部分图像，称为Slab切块图像，每个Slab切块图像是由若干同时重建出来的薄层二维切片图像组成的三维厚片图像，相邻Slab切块图像相接部分的图像不一致是心脏图像重建的普遍问题；即使相邻Slab切块图像所选的生数据部分为平稳区，由于心脏相邻周期路径不一定相同，同样存在图像不一致的问题。

技术实现要素：

鉴于上述原因，本发明的目的在于提供一种Slab切块图像的连接方法，可解决相邻Slab切块图像相接部分图像不一致的问题，实现心脏图像的CT重建。

为实现上述目的，本发明采用以下技术方案：

一种Slab切块图像的连接方法，包括，利用ECG心跳信号门控方法，得到用于CT重建的多个Slab切块图像，

将多个Slab切块图像依序相接，

相邻的两Slab切块图像，计算相接处对应的两两图像的相似度，根据相似度值对两两图像进行插值计算，得到相接处的最终图像。

计算两图像的相似度的方法是，利用相关系数衡量两幅图像中感兴趣区域的相似度，利用余弦相似度衡量两幅图像整体的灰度信息匹配程度；

计算相关系数的公式为：

其中，t和r分别表示两幅图像的像素值，N表示所选感兴趣区域的窗口边长，和分别表示位于窗内的像素均值，对于多个感兴趣区域，计算相应的相关系数，并计算多个相关系数的均值

计算余弦相似度的公式为：

其中，A和B分别是两幅图像的直方图信息。

根据相似度值对两图像进行插值计算的方法是，

w2＝1-w1

插值计算后的图像＝图像1*w1+图像2*w2

其中，w1表示图像1的权重值，w2表示图像2的权重值。

对于相邻的第一Slab切块图像与第二Slab切块图像，第一Slab切块图像的末端图像依序为图像m-1，m，第二Slab切块图像的前端图像依序为图像1，2，图像m-1，m与1，2无重叠，则，分别计算图像m-1与图像1、图像m与图像2的相似度值，然后利用相应的相似度值对相应的两图像进行插值计算，得到相应的最终图像。

对于相邻的第一Slab切块图像与第二Slab切块图像，第一Slab切块图像的末端图像依序为图像m-2，m-1，m，第二Slab切块图像的前端图像依序为图像1，2，3，图像m-1与图像1重叠，图像m与图像2重叠，则，分别计算图像m-1与图像1、图像m与图像2的相似度值，然后利用相应的相似度值对相应的两图像进行插值计算，得到相应的最终图像。

本发明的优点是：

本发明的Slab切块图像的连接方法，可解决相邻Slab切块图像相接部分图像不一致的问题，CT重建的心脏图像更加清晰，心脏的冠状面、矢位面及横截面血管位置更加准确，避免发生血管错位情况。

附图说明

图1是EGG心跳信号时间序列与重建位置的对应关系示意图。

图2是本发明的相邻Slab切块图像相接处无重叠的情况示意图。

图3是本发明的相邻Slab切块图像相接处有部分重叠的情况示意图。

具体实施方式

以下结合附图和实施例对本发明作进一步详细的描述。

本发明公开的Slab切块图像的连接方法，包括：

1、利用ECG心跳信号门控方法，得到用于CT重建的多个Slab切块图像；

2、将多个Slab切块图像依序相接，对于相邻的两Slab切块图像，计算相接处多幅图像的相似度，根据相似度值对图像进行插值计算，得到相接处的最终图像。其中，

1)对于相邻Slab切块图像相接处图像无重叠的情况

如图2所示，Slab3切块图像与Slab4切块图像为相邻的两Slab切块图像，图像21、22为Slab3切块图像的末端图像，图像1、2为Slab4切块图像的前端图像，即，图像21、22、1、2为该两Slab切块图像相接处的图像，该两Slab切块图像相接时，计算图像21与图像1的相似度，得到的相似度值作为线性插值系数对该两幅图像进行插值计算，得到插值计算后的图像代替Slab3切块图像的图像22；相类似的，计算图像22与图像2的相似度，得到的相似度值作为线性插值系数对该两幅图像进行插值计算，得到插值计算后的图像代替Slab4切块图像的图像1；之后，将Slab3切块图像与Slab4切块图像相接即可。

计算两幅图像的相似度的方法是，采用对两种度量相似度的函数进行混合计算，包括，利用相关系数衡量两幅图像中所感兴趣区域内的相似度，利用余弦相似度衡量两幅图像整体的灰度信息匹配程度。

计算相关系数的公式为：

其中，t和r分别表示两幅图像的像素值，N表示所选感兴趣区域的窗口边长，和分别表示位于窗内的像素均值，感兴趣区域可以是多个，计算多个感兴趣区域对应的相关系数，然后计算多个相关系数的均值

计算余弦相似度的公式为：

其中，A和B分别是两幅图像的直方图信息，通过直方图计算余弦值，越接近1图像相似度越高。

线性插值系数的计算公式为：

w2＝1-w1

插值计算后的图像＝图像1*w1+图像2*w2

其中，w1表示图像1的权重值，w2表示图像2的权重值。

2)对于相邻Slab切块图像相接处有部分图像重叠的情况

如图3所示，Slab3切块图像与Slab4切块图像为相邻的两Slab切块图像，图像20、21、22为Slab3切块图像的末端图像，图像1、2、3为Slab4切块图像的前端图像，Slab3切块图像与Slab4切块图像相接处的图像部分重叠，即，图像21与图像1重叠，图像22与图像2重叠。此种情况下，按照上述方法，计算重叠处的图像21与图像1的相似度，计算结果作为线性插值系数，再利用该线性插值系数对该两幅图像进行插值计算，得到插值计算后的图像，代替Slab3切块图像的图21及Slab4切块图像的图像1；同理，计算重叠处的图像22与图像2的相似度，计算结果作为线性插值系数，再利用该线性插值系数对该两幅图像进行插值计算，得到插值计算后的图像，代替Slab3切块图像的图22及Slab4切块图像的图像2；之后，将Slab3切块图像与Slab4切块图像相接，重叠部分的图像择一即可。

以上所述是本发明的较佳实施例及其所运用的技术原理，对于本领域的技术人员来说，在不背离本发明的精神和范围的情况下，任何基于本发明技术方案基础上的等效变换、简单替换等显而易见的改变，均属于本发明保护范围之内。