基于电阻抗成像的三维血液灌注图像产生方法与装置与流程

本发明涉及电阻抗成像技术，更具体涉及一种基于电阻抗成像的三维血液灌注图像产生方法与装置。

背景技术：

电阻抗成像(electricalimpedancetomography，eit)是一种无创的、以人体或其他生物体内部的电阻率分布为目标的重建体内组织图像的技术。人体是一个大的生物电导体，各组织、器官均有一定的阻抗，当人体的局部器官发生病变时，局部部位的阻抗必然与其他部位不同，因而可以通过阻抗的测量来对人体器官的病变进行诊断。

现有技术的电阻抗成像方法只能重建二维的血液灌注图像。这个二维的图像反映的是待测人体区域某个断面内由于血液灌注引起的电阻抗变化。二维的图像难以反映三维空间中某个体积区域内的血液灌注情况。

因此，希望能够提供一种基于电阻抗成像的三维血液灌注图像产生方法与装置，该方法与装置能够基于电阻抗测量结果，重建关于血液灌注的三维图像。

技术实现要素：

如前所述，为了解决现有技术中存在的问题，本发明的目的是提供一种三维血液灌注图像的生成方法与装置。所述方法和装置通过对待测人体区域进行电阻抗测量，从测量信号中提取血液灌注信号，最后重建三维血液灌注图像；或者先根据测量信号重建三维差分图像，然后从三维差分图像中提取测量信号中的血液灌注信号所反映的三维血液灌注图像。利用所述方法和装置产生三维血液灌注图像可以通过显示设备进行显示。

根据本发明的实施例，本发明的第一方面提供了一种基于电阻抗成像的三维血液灌注图像产生方法。所述方法可以包括如下步骤：利用分布在三维空间内的电极阵列，对待测人体区域进行电阻抗测量，得到电阻抗测量信号；基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像。

根据本发明的第一方面的方法，优选地，可以从电阻抗测量信号中提取血液灌注信号，然后利用所提取的血液灌注信号，通过图像重建算法重建三维血液灌注图像。

在此情况下，所述的从电阻抗测量信号中提取血液灌注信号可以进一步包括利用信号的时频特性从电阻抗测量信号中提取血液灌注信号。

进一步，优选地，所述的利用信号的时频特性从电阻抗测量信号中提取血液灌注信号可以包括使用带通滤波器从电阻抗测量信号中分离特定频率范围的信号，以形成血液灌注信号。

根据本发明的第一方面的方法，另一方面，也可以根据电阻抗测量信号，通过图像重建算法重建三维差分图像，然后从三维差分图像中提取电阻抗测量信号中的血液灌注信号所反映的三维血液灌注图像。

在此情况下，所述的从三维差分图像中提取电阻抗测量信号中的血液灌注信号所反映的三维血液灌注图像可以进一步包括利用三维差分图像中像素的时频特性提取三维血液灌注图像。

进一步，优选地，所述的利用三维差分图像中像素的时频特性提取三维血液灌注图像可以包括对三维差分图像中每个像素的时域信号，使用带通滤波器从中分离出特定频率范围的信号，以形成三维血液灌注图像相应像素的时域信号。

根据本发明的第一方面的方法，优选地，所述电极阵列可以设置在一条或多条阻抗带上、电极背心上或电极头罩上，使得电极呈三维分布。

根据本发明的实施例，本发明的第二方面提供了一种基于电阻抗成像的三维血液灌注图像产生装置。所述装置可以包括：分布在三维空间内的电极阵列，用于对待测人体区域进行电阻抗测量，得到电阻抗测量信号；图像重建处理器，其执行存储于存储器中的程序，以便使得：基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像。

根据本发明的第二方面的装置，优选地，所述图像重建处理器可以进一步执行存储于存储器中的程序，以便使得：从电阻抗测量信号中提取血液灌注信号；以及利用所提取的血液灌注信号，通过图像重建算法重建三维血液灌注图像。

根据本发明的第二方面的装置，另一方面，所述图像重建处理器也可以进一步执行存储于存储器中的程序，以便使得：通过图像重建算法重建三维差分图像；以及从三维差分图像中提取电阻抗测量信号中的血液灌注信号所反映的三维血液灌注图像。

根据本发明的第二方面的装置，优选地，所述电极阵列可以设置在一条或多条阻抗带上、电极背心上或电极头罩上，使得电极呈三维分布。

根据本发明实施例的血液灌注图像产生方法和装置可以产生血液灌注引起的电阻抗变化的三维图像，相比于现有技术的二维图像，可以更加直观地反映出人体区域三维空间内一个体积区域上的血液灌注情况，有利于图像分析对比、疾病检测和诊断。

附图说明

下面参考附图结合实施例说明本发明。

图1是根据本发明实施例的基于电阻抗成像的三维血液灌注图像产生方法的示意流程图。

图1a和图1b分别是根据本发明实施例的基于电阻抗成像的三维血液灌注图像产生方法的两种优选实施例的示意流程图。

图2是根据图1a的优选实施例图示说明人体胸腔测量数据实例，其中，(a)是时域信号，(b)是频域信号。

图3是根据图1a的优选实施例图示说明通过滤波提取出的灌注相关信号，其中，(a)为时域信号，(b)为频域信号。

图4是根据图1b的优选实施例图示说明三维差分图像以及像素点的时频信号示例，其中，(a)为三维差分图像，(b)为示例像素的时域信号，(c)为示例像素的频域信号。

图5是根据图1b的优选实施例图示说明通过滤波提取出的示例像素的灌注相关信号，其中，(a)为时域信号，(b)为频域信号。

图6是根据本发明的优选实施例的人体肺部三维血液灌注图像示意图。

图7是根据本发明实施例的基于电阻抗成像的三维血液灌注图像产生装置的示意框图。

具体实施方式

附图仅用于示例说明，不能理解为对本发明的限制。下面结合附图和实施例对本发明的技术方案做进一步的说明。

图1是根据本发明实施例的基于电阻抗成像的三维血液灌注图像产生方法100的示意流程图。

图1的方法100开始于步骤110，在此步骤，测量人体的电阻抗信号。具体地说，利用分布在三维空间内的电极阵列，对待测人体区域进行电阻抗测量，得到电阻抗测量信号。

电阻抗测量首先需要在待测人体区域周围固定电极阵列。所述电极阵列包含若干个分布在三维空间内的电极。然后，通过电极阵列对待测人体区域进行激励并测量由此产生的响应，如：轮流对电极施加电流激励，并依次在其他电极上测量由此产生的电压信号。

具体地说，在一个实施例中，将包含电极的传感模块固定在人体被测量部位，如胸腔、脑部、腹部或四肢的周围，采用阻抗带、电极背心、电极头罩等电极阵列的形式。在有些实施例中，电极可以采用体内电极的形式。所谓的体内电极是指，将该电极置入人体食道、气管等人体内位置。

根据本发明的优选实施例，可以利用一条或多条阻抗带、电极背心、电极头罩等形式的电极阵列进行信号测量。也就是说，所述电极阵列设置在一条或多条阻抗带上、电极背心上或电极头罩上，使得电极呈三维分布。具体地说，为了使重建图像具有三维分辨率，所述电极阵列通常分布在三维空间内，而不是分布在一个二维平面或近似的二维平面内。为了使电极阵列分布在三维空间内，可以采用多条阻抗带或采用电极呈三维分布的电极背心、电极头罩等方案。

测量得到的电阻抗测量信号可以是电压信号，具体可以是复电压信号。复电压信号可以用幅度和相位的形式表达，或者可以用实部和虚部的形式表达。

接下来，在步骤120，基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像。

步骤120可以通过两种方式来实现。图1a和图1b分别是根据本发明实施例的基于电阻抗成像的三维血液灌注图像产生方法的两种优选实施例100a和100b的示意流程图。

如图1a中所示，根据本发明的一种优选实施例100a，在步骤110之后，在步骤121a，从电阻抗测量信号中提取血液灌注信号。

在此步骤中，需要从上一步骤获得的电阻抗测量信号中提取血液灌注信号。根据本发明的优选实施例，可以利用信号的时频特性从上一步骤获得的电阻抗测量信号中提取血液灌注信号。具体地，利用滤波器从所测量的电阻抗信号中分离血液灌注信号。

下面以一个人体胸腔测量信号的例子来说明。

图2是根据图1a的优选实施例图示说明人体胸腔测量数据实例。图2(a)显示了某个测量数据的时域图形。图中曲线代表特定电极激励时，在特定电极上测量得到的电压信号。其他激励-测量情况得到的数据与之类似。需要说明的是，图中纵坐标为从数字电压表直接读取的数值，仍未将其转换为电压值。图2(b)显示了该测量信号的频域图形。具体地，可以由图2(a)中信号经过傅里叶变换得到。从图2(b)可以分辨出由呼吸引起的信号分量和由灌注引起的信号分量。为了提取灌注信号分量，可以设计一个带通滤波器。

滤波后的信号如图3所示。图3是根据图1a的优选实施例图示说明通过滤波提取出的灌注相关信号。其中，图3(a)对应滤波后信号的时域图形，图3(b)对应滤波后信号的频域图形。

在上述例子中，使用带通滤波器从电阻抗测量信号中分离特定频率范围的信号，以形成血液灌注信号。

在步骤122a，利用所提取的血液灌注信号，通过图像重建算法重建三维血液灌注图像。

一般来说，如果以肺部血液灌注图像的重建为例，其重建过程为：首先从测量数据中提取灌注信号，然后利用不同时刻灌注信号的差进行图像重建。所述三维血液灌注图像反映由于血液灌注引起的待测人体区域内的电阻抗变化，例如电导率的变化。由此，相应地展示出不同时刻肺部血液含量的变化。

在本发明的一个优选实施例中，所述图像重建算法是线性差分成像算法。下面就以一个线性差分成像算法的例子来说明。

假设上一步骤所提取出的灌注信号的时域形式为u(t)，其中t为时间变量。eit差分重建可以表述为如下最小二乘问题：

minδσ||j·δσ-δu||²+α||r·δσ||²，

其中，j为雅可比矩阵，δu＝u(t2)-u(t1)为信号在时刻t2相对于时刻t1的变化，δσ为上述两个时刻由于血液灌注引起的电导率变化，r为正则化矩阵，α为正则化参数。δσ定义在离散化的三维模型中，如四面体网格或体素网格。

上述问题的解为δσ*＝(j^t·j+αr^t·r)^-1·j^t·δu.

令d＝(j^t·j+αr^t·r)^-1·j^t，则上述公式可以表示为：

δσ*＝d·δu.

上述δσ*即为计算所得的血液灌注图像。

尽管在上述例子中，具体使用了线性差分成像算法作为图像重建算法来计算并重建血液灌注的三维图像，但是本领域普通技术人员应该理解，本发明可以利用的图像重建算法可以包含多种重建算法：线性的或非线性的、迭代的或非迭代的、随机的或确定性的图像重建算法等。

回到图1以及图1b。如图1b所示，图1的步骤120还包括另外一种实施方式100b，即在图1的步骤110之后，执行步骤121b，即根据电阻抗测量信号，通过图像重建算法重建三维差分图像。

在此步骤，所述图像重建算法可以使用与本发明的第一种优选实施例100a中相同的图像重建算法。

然后，在步骤122b，从三维差分图像中提取电阻抗测量信号中的血液灌注信号所反映的三维血液灌注图像。根据本发明的优选实施例，可以利用图像信号的时频特性从三维差分图像中提取三维血液灌注图像。具体地，利用滤波器从三维差分图像中分离三维血液灌注图像。

下面仍以上述人体胸腔测量信号的例子来说明。

图4是根据图1b的优选实施例图示说明重建出的三维差分图像。图4(a)是三维差分图像的示意图。图4(b)是从三维差分图像上选取了一个示例像素点，所绘制出的该像素点的时域信号。图4(c)是该像素点的频域信号。具体地，图4(c)中信号可以由图4(b)中信号经过傅里叶变换得到。由图4(c)可以分辨出由呼吸引起的信号分量和由灌注引起的信号分量。为了从中提取灌注信号分量，可以设计一个带通滤波器。

带通滤波器对上述三维差分图像中的每个像素的时域信号执行滤波操作。对上述示例像素点滤波后的信号如图5所示。其中，图5(a)是滤波后的时域信号。图5(b)是滤波后的频域信号。

对上述三维差分图像中的每个像素滤波之后，即可得到三维血液灌注图像。所述三维血液灌注图像反映由于血液灌注引起的待测人体区域内的电阻抗变化，例如电导率的变化。由此，相应地展示出不同时刻肺部血液含量的变化。图6显示了利用包括实施例100a和实施例100b在内的上述方法所产生的人体肺部的三维血液灌注图像的一个示意图。

由此可以看出，在图1的方法100的步骤120的两种实现方式之间，区别在于：先从电阻抗测量信号中提取血液灌注信号再进行图像重建，还是先根据测量信号重建三维差分图像再从图像中提取专门的灌注图像。本发明的保护范围希望包括这两种实现方式。

图7是根据本发明实施例的基于电阻抗成像的三维血液灌注图像产生装置700的示意框图。

根据本发明的实施例，参见图7，基于电阻抗成像的三维血液灌注图像产生装置700可以包括：分布在三维空间内的电极阵列710，用于对待测人体区域进行电阻抗测量，得到电阻抗测量信号；图像重建处理器720，其执行存储于存储器中的程序，以便使得：基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像。

根据本发明的一个优选实施例，所述电极阵列710可以设置在一条或多条阻抗带上、电极背心上或电极头罩上，使得电极呈三维分布。

根据本发明的优选实施方式，所述图像重建处理器720进一步执行存储于存储器中的程序，以便使得：从电阻抗测量信号中提取血液灌注信号；以及利用所提取的血液灌注信号，通过图像重建算法重建三维血液灌注图像。

具体地说，所述图像重建处理器720可以利用信号的时频特性从电阻抗测量信号中提取血液灌注信号。更具体地说，所述图像重建处理器720可以使用带通滤波器从电阻抗测量信号中分离特定频率范围的信号，以形成血液灌注信号。

在本发明的一个具体实施例中，所述图像重建算法是线性差分成像算法。然而，本领域技术人员应该理解，本发明可以利用的图像重建算法可以包含多种重建算法：线性的或非线性的、迭代的或非迭代的、随机的或确定性的图像重建算法等。

另一方面，根据本发明的另一种实现方式，所述图像重建处理器720进一步执行存储于存储器中的程序，以便使得：通过图像重建算法重建三维差分图像；以及从三维差分图像中提取电阻抗测量信号中的血液灌注信号所反映的三维血液灌注图像。

具体地说，所述图像重建处理器720在重建三维差分图像时，可以使用与第一种实现方式中相同的图像重建算法。所述图像重建处理器720可以利用所述三维差分图像中像素的时频特性从中提取三维血液灌注图像。更具体地说，所述图像重建处理器720可以使用带通滤波器从所述三维差分图像中分离三维血液灌注图像。

此外，尽管图7中未示出，本领域技术人员应该能够认识到，装置700还可以包括显示单元，例如lcd显示器，用于呈现图像重建后的三维血液灌注图像，其显示的示意例如如图6中所示。由于本发明更关注于三维血液灌注图像的产生，因此显示单元并非本发明的方法或装置所必不可少的元素。本发明的范围应以权利要求书为准，而不受到本发明说明书中记载或未记载的任何实施例的限制。

此外，本领域普通技术人员应该认识到，本发明的方法可以实现为计算机程序。如上结合图1、1a、1b所述，通过一个或多个程序执行上述实施例的方法，包括指令来使得计算机或处理器执行结合附图所述的算法。这些程序可以使用各种类型的非瞬时计算机可读介质存储并提供给计算机或处理器。非瞬时计算机可读介质包括各种类型的有形存贮介质。非瞬时计算机可读介质的示例包括磁性记录介质(诸如软盘、磁带和硬盘驱动器)、磁光记录介质(诸如磁光盘)、cd-rom(紧凑盘只读存储器)、cd-r、cd-r/w以及半导体存储器(诸如rom、prom(可编程rom)、eprom(可擦写prom)、闪存rom和ram(随机存取存储器))。进一步，这些程序可以通过使用各种类型的瞬时计算机可读介质而提供给计算机。瞬时计算机可读介质的示例包括电信号、光信号和电磁波。瞬时计算机可读介质可以用于通过诸如电线和光纤的有线通信路径或无线通信路径提供程序给计算机。

因此，根据本发明，还可以提议一种计算机程序或一种计算机可读介质，用于记录可由处理器执行的指令，所述指令在被处理器执行时，使得处理器执行一种基于电阻抗成像的三维血液灌注图像产生方法，包括如下操作：利用分布在三维空间内的电极阵列，对待测人体区域进行电阻抗测量，得到电阻抗测量信号；基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像。

上面已经描述了本发明的各种实施例和实施情形。但是，本发明的精神和范围不限于此。本领域技术人员将能够根据本发明的教导而做出更多的应用，而这些应用都在本发明的范围之内。

也就是说，本发明的上述实施例仅仅是为清楚说明本发明所做的举例，而非对本发明实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其他不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、替换或改进等，均应包含在本发明权利要求的保护范围之内。