计算机断层扫描图像重建方法及装置与流程

【技术领域】

本发明涉及计算机断层扫描(computedtomography,ct)技术领域，尤其涉及一种计算机断层扫描图像重建方法及装置。

背景技术：

计算机断层扫描(computedtomography，简称ct)是用x射线对人体的特定部位按一定厚度的层面进行扫描，由于不同的人体组织对x射线的吸收能力不同，可以用计算机重建出断层面的影像。

在进行心脏ct扫描时，为了消除心脏运动对图像的影响，减少重建的图像中的运动伪影，会在扫描的同时监控心电(ecg)信号，并通过一条与所述心电信号相关的门控曲线(加权函数曲线)对扫描数据进行加权，对心脏运动幅度较小时采集的扫描数据采用较大的权重参与图像重建，而对心脏运动幅度较大时采集的扫描数据采用较小的权重参与图像重建，以得到比较清晰的心脏图像。

使用门控曲线对于扫描数据进行加权时，不同角度(view)采样的扫描数据有一个对应的权重值。而该门控曲线的产生是跟重建视野(fov)相关的，不同重建视野会产生不同的门控曲线，对应不同的时间分辨率。一般地，当重建视野较小时，其图像的像素对应的时间分辨率较高；而当重建视野较大时，为了保证大视野图像边缘处的像素有足够的投影数据，需要使得其对应的相位处加权函数宽度增加，即在该相位处用到了更长时间的投影数据参与图像重建，从而降低了图像的时间分辨率。由于心脏冠脉重建对时间分辨率要求较高，往 往要求重建视野较小。但小视野图像重建不能满足所有心脏扫描的临床实验要求，对于一些扫描协议，如胸痛三联征的扫查，需要大视野重建以观察除心脏周围的一些肺部组织，同时希望位于图像中心的心脏冠脉血管成像清晰。而现有的技术如果按照大视野图像重建往往会减低图像中心区域的冠脉血管的时间分辨率；而如果按照小视野图像重建则会导致图像边缘肺部区域因为数据不足而发黑，或者得不到肺部区域的图像。

因此，需要提出一种新的计算机断层扫描图像重建方法，在重建的图像具有较大重建视野的同时，在图像的中心区域具有较高的时间分辨率。

技术实现要素：

本发明解决的是现有的计算机断层扫描图像重建方法，特别是心脏扫描图像重建的视野过大时，图像分辨率下降的问题。为解决上述问题，本发明提出了一种计算机断层扫描图像重建方法，包括以下步骤：

进行扫描，获得扫描数据，以及生理信号；

获取第一视野图像重建范围，产生第一加权函数，利用所述第一加权函数对所述扫描数据进行加权，对经过第一加权函数加权后的扫描数据进行第一视野图像重建，获得第一图像；

获取第二视野图像重建范围，产生第二加权函数，利用所述第二加权函数对所述扫描数据进行加权，对经过第二加权函数加权后的扫描数据进行第二视野图像重建，获得第二图像；

合并第一图像和第二图像，获得第三图像。

可选地，所述生理信号为心电信号或呼吸信号。

可选地，还包括对所述扫描数据进行预处理的步骤。

可选地，所述合并第一图像和第二图像，包括：对第一图像和第二图像进 行加权合并，得到第三图像。

可选地，所述第一视野小于第二视野，选取第一图像中的一个圆形区域参与合并。

可选地，所述第一视野小于第二视野，将所述第一图像按比例缩小到与第二图像相同的尺度后参与合并。

可选地，所述对第一图像和第二图像进行加权合并，包括：对所述第一图像和第二图像分别进行线性加权，将线性加权后的第一图像和第二图像相加，得到第三图像。

可选地，所述第一视野小于第二视野，选取第一图像中的一个圆形区域参与合并，所述圆形区域的半径为r1，对第一图像和第二图像分别乘以第三加权函数w3(r)和第四加权函数w4(r)后再相加，获得第三图像；

所述第三加权函数w3(r)为

所述第四加权函数w4(r)为：

其中，r为半径，且r0<r1。

本发明还提出了一种计算机断层扫描图像重建装置，包括：

扫描单元，用于进行扫描，获得扫描数据，以及生理信号；

第一图像获得单元，用于获得第一视野图像重建范围，产生第一加权函数，利用所述第一加权函数对所述扫描数据进行加权，对经过第一加权函数加权后 的扫描数据进行第一视野图像重建，获得第一图像；

第二图像获得单元，用于获取第二视野图像重建范围，产生第二加权函数，利用所述第二加权函数对所述扫描数据进行加权，对经过第二加权函数加权后的扫描数据进行第二视野图像重建，获得第二图像；

合并单元，用于合并第一图像和第二图像，获得第三图像。

可选地，还包括预处理单元，用于对所述扫描数据进行预处理。

本发明对比现有技术有如下的有益效果：

采用本发明的图像重建方法，可以兼顾图像的时间分辨率和数据的完备性，从而使得图像中心处有较高的时间分辨率，而图像整体不会出现数据不足的情况，提高图像质量。

进一步地，在本发明的可选实施方式中，对第一图像和第一图像进行加权合并，从而可以避免的在第一图像和第二图像的边界线处产生伪影。

【附图说明】

图1是本发明的计算机断层扫描系统的结构图；

图2是本发明的一实施例的计算机断层扫描图像重建方法的流程图；

图3是本发明的心电信号和与其对应的加权函数曲线的示意图；

图4是本发明的一实施例的图像合并方法的流程图；

图5是本发明的一实施例的图像合并方法的示意图；

图6是本发明的一实施例的计算机断层扫描图像重建装置的结构图；

图7是小视野图像重建、大视野图像重建以及使用本发明的图像重建方法获得的图像的对比图。

【具体实施方式】

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对 本发明的具体实施方式做详细的说明。

图1是一种计算机断层扫描系统的结构示意图，如图1所示，计算机断层扫描系统100包括机架110，机架110具有围绕系统轴线旋转的可旋转的部分130。可旋转的部分130具有相对设置的x射线源131和x射线探测器132的x射线系统。还具有检查床120，在进行检查时，受检者在该检查床120上可以沿着z轴方向被推入到扫描腔体133中。x射线源131绕z轴旋转，探测器132相对于x射线源131一起运动，以采集投影测量数据，这些数据在之后被用于重建图像。还可以进行螺旋扫描，在螺旋扫描期间，通过受检者沿着z轴的连续运动和x射线源131的同时旋转，x射线源131相对于受检者产生螺旋轨迹。高压发生单元134与射线源131相连，以提供电源。

生理信号监控单元141用于监控扫描受检者的生理信号，例如：心电信号或呼吸信号。处理单元142连接探测器132以获得受检者的投影测量数据，供后续处理，例如：图像重建及处理。控制单元140连接高压发生单元134以控制射线源131的扫描过程。控制台及显示器143用以呈现界面、数据和影像给使用者。控制单元140还连接处理单元142和控制台及显示器143以控制该部件的运作。

图2是本发明的一实施例的计算机断层扫描图像重建方法的流程图，如图2所示，本发明的计算机断层扫描图像重建方法包括以下步骤：

步骤201，进行扫描，获得扫描数据，以及生理信号。

优选地，所述生理信号为心电(ecg)信号，所述扫描数据为投影数据。

在本发明的其他实施例中，所述生理信号也可以是呼吸信号。

步骤202，对扫描数据进行预处理。

在一些实施例中，在步骤202中的扫描数据预处理可以是对扫描数据采集 过程中因x射线源、探测器、机械、特殊的数据采集或病人等原因引起的不完善进行补偿的方法，例如，空气校准，中心校准，探测器增益校准，扫描数据降噪处理，或以上方法的任意组合。

在本发明的一些实施例中，所述步骤202的扫描数据预处理过程也可以省略，而直接执行步骤203。

步骤203，获得小视野(fov)图像重建范围，产生第一加权函数，利用所述第一加权函数对所述扫描数据进行加权，对经过第一加权函数加权后的扫描数据进行小视野图像重建，获得小视野图像。

优选地，所述小视野图像重建范围为200mm×200mm的矩形(fov＝200)。

图3是本发明的心电信号和与其对应的加权函数曲线的示意图。图3上半部分为步骤201中获得的心电(ecg)信号的曲线，其横坐标为时间t，纵坐标为心跳幅度，其代表了心跳幅度随时间的变化；图3的下半部分为与所述心电信号对应的加权函数曲线。

通常来说，对心脏运动幅度较小时采集的扫描数据采用较大的权重，从而使得在心脏运动幅度较小时采集的扫描数据更多地参与图像重建；而对心脏运动幅度较大时采集的扫描数据采用较小的权重，从而使得在心脏运动幅度较大时采集的扫描数据较少的参与图像重建。以上方法可以减少心脏运动对于重建图像的影像，减少图像中的运动伪影。

如图3所示，在本实施例中，对于心脏运动幅度较大时对应的权重系数为0，对于心脏运动幅度较大时对应的权重系数为1。

本实施例以心脏的运动为例进行了说明，但也可以适用于其他情况的扫描对象的运动，例如肺部(呼吸)或腹部的运动。

在本发明的一些实施例中，所述第一加权函数为w1(z),其中z为心脏的相 位，z与时间t和进行ct扫描的视角(view)分别对应。将所述第一加权函数w1(z)与对应的扫描数据相乘，将相乘后的数据进行小视野图像重建，获得小视野图像。

步骤204，获得大视野(fov)图像重建范围，产生第二加权函数，利用所述第二加权函数对所述扫描数据进行加权，对经过第二加权函数加权后的扫描数据进行大视野图像重建，获得大视野图像。

优选地，所述大视野图像重建范围为350mm×350mm的矩形(fov＝350)。

在本发明的一些实施例中，所述第二加权函数为w2(z),其中z为心脏的相位，其与时间和进行ct扫描的视角(view)分别对应。将所述加权函数w2(z)与对应的扫描数据相乘，将相乘后的数据进行大视野图像重建，获得大视野图像。

关于产生加权函数的方法，在上文中已具体说明，这里不再敷述。

需要说明的是，第一或第二加权函数不仅和心脏信号相关，在一些实施例中，也和视野大小相关，基于不同的视野大小会产生不同的加权函数，简单来说，当视野增大时，为了保证大视野图像边缘处的像素有足够的扫描数据参与重建，需要使得其对应的相位处加权函数的窗口宽度(如图3所示)增加，即在该相位处用到了更长时间的扫描数据参与图像重建。

如图3所示，将一个心动周期对应的加权函数曲线中加权系数不为零的部分视为一个矩形(也可以称为“选通窗口”)。在不考虑其他因素的情况下，心脏信号(ecg)决定了该矩形的位置，通常将心跳最平缓处的相位点作为该矩形的中心线位置；而重建的视野大小决定了该矩形的宽度w，重建视野越大则该矩形的宽度也越大，反之亦然。即在本实施例中，第一加权函数曲线对应的选通窗口的中心线和第二加权函数曲线对应的选通窗口的中心线重合，而第一加权函数曲线对应的选通窗口的宽度小于第二加权函数曲线对应的选通窗口 的宽度。

在本发明的一些实施例中，上述步骤203或步骤204中获取小视野(fov)图像重建范围，或者获取大视野(fov)图像重建范围的方法可以是预先设定重建视野，也可以由操作者根据实际需求在通过如图1所示的控制台及显示器143上设置重建视野而获得。

上述步骤203或步骤204中的对加权后的扫描数据进行图像重建的方法可以是但不限于，滤波反投影方法，迭代重建方法等。

步骤205，合并小视野图像和大视野图像，得到最终的图像。

图4是本发明的一实施例的图像合并方法的流程图，图5是本发明的一实施例的图像合并方法的示意图。

请参照图4和图5，所述步骤205中的合并小视野图像和大视野图像的方法具体可以包括以下步骤：

步骤401，选取小视野图像中的一个圆形区域。

所述圆形区域的半径为r1。

优选地，选取小视野图像中能够选取的最大的一个圆形区域，例如，所述小视野图像为r1×r1大小的矩形，则选择半径为r1的圆形区域。

步骤402，将所述小视野图像的圆形区域图像按比例缩小到与大视野图像相同的尺度。

在本发明的一些实施例中，可以先执行步骤402再执行步骤401，即先将小视野图像按比例缩小到与大视野图像相同的尺度，再选择小视野图像中的一个圆形区域。

步骤403，将缩小后的小视野图像与大视野图像进行加权合并。

优选地，所述加权合并的权值是按照图像的半径方向进行，将图像分为半 径大于r1的区域，半径大于等于r0小于等于r1(r0<r1)的区域，以及半径小于r0的区域。

对小视野图像和大视野图像(的像素值)分别乘以线性加权函数后再相加，获得最终的合并后的图像。

具体地，对小视野图像乘以第三加权函数w3(r)，对于大视野图像乘以第四加权函数w4(r)，所述第三加权函数和第四加权函数为线性函数。

第三加权函数w3(r)为：

第四加权函数w4(r)为：

需要注意的是，上面所描述的流程图是为插图的目的而提供的，而不是用来限制本发明的范围。对于本领域技术人员，在本公开的范围中，可以进行各种各样的变更和修改，这些变化和修改应该不会偏离本公开的保护范围。例如，所述步骤203获得小视野图像的过程可以在所述步骤204获得大视野图像的过程之前或同时执行，所述步骤201、202、203或步骤204可以分解为多个子步骤,等等。

优选地，本发明的图像重建方法适用于心脏ct扫描，特别是胸痛三联症扫描。但是本领域技术人员应当理解，本发明的图像重建方法还可以用于其他和生理信号相关或可能由于受检者本身的运动而产生运动伪影的ct扫描的图像重建。

本发明还提出了一种计算机断层扫描图像重建装置，如图6所示，所述计 算机断层扫描图像重建装置600包括：

扫描单元601，用于进行扫描，获得扫描数据，以及与所述扫描数据同步的生理信号。

所述生理信号可以是但不限于心电信号或呼吸信号。

第一图像获得单元602，所述第一图像获得单元602与所述扫描单元601相连，用于获得第一视野图像重建范围，产生第一加权函数，利用所述第一加权函数对所述扫描数据进行加权，对经过第一加权函数加权后的扫描数据进行第一视野图像重建，获得第一图像；

第二图像获得单元603，所述第二图像获得单元603与所述扫描单元601相连，用于获取第二视野图像重建范围，产生第二加权函数，利用所述第二加权函数对所述扫描数据进行加权，对经过第二加权函数加权后的扫描数据进行第二视野图像重建，获得第二图像；

合并单元604，所述合并单元604与所述第一图像获得单元602和第二图像获得单元603分别相连，用于合并第一图像和第二图像，获得第三图像。

在一些实施方式中，所述合并单元604可用于执行如下步骤：

选取小视野图像中的一个圆形区域；将所述小视野图像的圆形区域图像按比例缩小到与大视野图像相同的尺度；将缩小后的小视野图像与大视野图像进行加权合并，获得最终的图像。

优选地，所述加权合并的权值是按照图像的半径方向进行，对于小视野图像和大视野图像分别进行线性加权，将线性加权后的小视野图像和大视野图像相加，获得最终的图像。

在本发明的一些实施例中，所述计算机断层扫描图像重建装置600还可以包括预处理单元605，用于对所述扫描数据进行预处理；所述预处理单元605 与所述扫描单元601、第一图像获得单元602以及第二图像获得单元603分别相连。

在本发明的一些实施方式中，所述第一图像获得单元602，第二图像获得单元603，合并单元604以及预处理单元605具体可以通过如图1所示的处理单元142实现。

本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于计算机可读存储介质中，存储介质可以包括但不限于：软盘、光盘、磁光盘、rom(只读存储器)、ram(随机存取存储器)、eprom(可擦除可编程只读存储器)、eeprom(电可擦除可编程只读存储器)、磁卡或光卡、闪存、或适于存储机器可执行指令的其他类型的介质/机器可读介质。

图7是小视野图像重建、大视野图像重建以及使用本发明的图像重建方法获得的图像的对比图，其中，图7(a)是使用大视野(fov＝350，即350mm×350mm的矩形)图像重建获得的图像；图7(b)是使用小视野(fov＝200，即200mm×200mm的矩形)图像重建获得的图像；图7(c)是使用本发明的图像重建方法获得的图像。如图7所示，使用本发明的图像重建方法获得的图像相较于大视野图像重建获得的图像具有更高的分辨率，例如：对于图7(a)和图7(c)中箭头所指出的位于图像中心的组织，图7(c)的成像质量明显好于图7(a)。并且，使用本发明的图像重建方法获得的图像相较于图7(b)所示的小视野图像重建获得的图像具有更大的视野，能够更好的符合实际临床需求。

以上，仅以示例方式阐释了可使用本发明所提供的图像重建方法的计算机断层扫描设备，本领域技术人员应当理解，如使用x射线的c型臂系统等设备， 或组合式医学成像系统(例如：组合式正电子发射断层成像-计算机断层成像，positronemissiontomography-computedtomographytomography,pet-ct)，或使用其它类型射线的断层成像设备等，均可适用本发明所述的计算机断层扫描图像重建方法和装置，本发明对计算机断层扫描设备的类型与结构并不做具体限定。

本发明中，各实施例采用递进式写法，重点描述与前述实施例的不同之处，各实施例中的相同方法或结构参照前述实施例的相同部分。

本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。