CN111067521A - 基于电阻抗成像的三维血液灌注图像产生方法与装置 - Google Patents

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

Info

Publication number
CN111067521A
CN111067521A CN201911424085.5A CN201911424085A CN111067521A CN 111067521 A CN111067521 A CN 111067521A CN 201911424085 A CN201911424085 A CN 201911424085A CN 111067521 A CN111067521 A CN 111067521A
Authority
CN
China
Prior art keywords
blood perfusion
signal
image
dimensional
electrical impedance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911424085.5A
Other languages
English (en)
Other versions
CN111067521B (zh
Inventor
张可
张昕
王谊冰
于洋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Beijing Huarui Boshi Medical Imaging Technology Co Ltd
Original Assignee
Beijing Huarui Boshi Medical Imaging Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Huarui Boshi Medical Imaging Technology Co Ltd filed Critical Beijing Huarui Boshi Medical Imaging Technology Co Ltd
Priority to CN201911424085.5A priority Critical patent/CN111067521B/zh
Publication of CN111067521A publication Critical patent/CN111067521A/zh
Priority to EP20909510.8A priority patent/EP4085828A4/en
Priority to US17/790,414 priority patent/US20230045401A1/en
Priority to PCT/CN2020/106754 priority patent/WO2021135211A1/zh
Priority to JP2022538885A priority patent/JP7430303B2/ja
Application granted granted Critical
Publication of CN111067521B publication Critical patent/CN111067521B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/053—Measuring electrical impedance or conductance of a portion of the body
    • A61B5/0536—Impedance imaging, e.g. by tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6802—Sensor mounted on worn items
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6802—Sensor mounted on worn items
    • A61B5/6804—Garments; Clothes
    • A61B5/6805—Vests
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/04—Arrangements of multiple sensors of the same type
    • A61B2562/046—Arrangements of multiple sensors of the same type in a matrix array

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Physiology (AREA)
  • Hematology (AREA)
  • Cardiology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

提供一种基于电阻抗成像的三维血液灌注图像产生方法与装置。所述方法(100)包括：利用分布在三维空间内的电极阵列，对待测人体区域进行电阻抗测量，得到电阻抗测量信号(110)；基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像(120)。由此，可以产生血液灌注引起的电阻抗变化的三维图像，相比于现有技术的二维图像，可以更加直观地反映出人体区域三维空间内一个体积区域上的血液灌注情况，有利于图像分析对比、疾病检测和诊断。

Description

基于电阻抗成像的三维血液灌注图像产生方法与装置
技术领域
本发明涉及电阻抗成像技术，更具体涉及一种基于电阻抗成像的三维血液灌注图像产生方法与装置。
背景技术
电阻抗成像(Electrical Impedance Tomography，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||2+α||R·δσ||2，
其中，J为雅可比矩阵，δu＝u(t2)-u(t1)为信号在时刻t2相对于时刻t1的变化，δσ为上述两个时刻由于血液灌注引起的电导率变化，R为正则化矩阵，α为正则化参数。δσ定义在离散化的三维模型中，如四面体网格或体素网格。
上述问题的解为δσ*＝(JT·J+αRT·R)-1·JT·δu.
令D＝(JT·J+αRT·R)-1·JT，则上述公式可以表示为：
δσ*＝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(随机存取存储器))。进一步，这些程序可以通过使用各种类型的瞬时计算机可读介质而提供给计算机。瞬时计算机可读介质的示例包括电信号、光信号和电磁波。瞬时计算机可读介质可以用于通过诸如电线和光纤的有线通信路径或无线通信路径提供程序给计算机。
因此，根据本发明，还可以提议一种计算机程序或一种计算机可读介质，用于记录可由处理器执行的指令，所述指令在被处理器执行时，使得处理器执行一种基于电阻抗成像的三维血液灌注图像产生方法，包括如下操作：利用分布在三维空间内的电极阵列，对待测人体区域进行电阻抗测量，得到电阻抗测量信号；基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像。
上面已经描述了本发明的各种实施例和实施情形。但是，本发明的精神和范围不限于此。本领域技术人员将能够根据本发明的教导而做出更多的应用，而这些应用都在本发明的范围之内。
也就是说，本发明的上述实施例仅仅是为清楚说明本发明所做的举例，而非对本发明实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其他不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、替换或改进等，均应包含在本发明权利要求的保护范围之内。

Claims (10)

1.一种基于电阻抗成像的三维血液灌注图像产生方法，包括如下步骤：
利用分布在三维空间内的电极阵列，对待测人体区域进行电阻抗测量，得到电阻抗测量信号；
基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像。
2.如权利要求1所述的方法，其中，所述的基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像进一步包括：
从电阻抗测量信号中提取血液灌注信号；以及
利用所提取的血液灌注信号，通过图像重建算法重建三维血液灌注图像。
3.如权利要求1所述的方法，其中，所述的基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像进一步包括：
根据电阻抗测量信号，通过图像重建算法重建三维差分图像；以及
从三维差分图像中提取电阻抗测量信号中的血液灌注信号所反映的三维血液灌注图像。
4.如权利要求1所述的方法，其中，所述电极阵列设置在一条或多条阻抗带上、电极背心上或电极头罩上，使得电极呈三维分布。
5.如权利要求2所述的方法，其中，所述的从电阻抗测量信号中提取血液灌注信号进一步包括：
利用信号的时频特性从电阻抗测量信号中提取血液灌注信号。
6.如权利要求5所述的方法，其中，所述的利用信号的时频特性从电阻抗测量信号中提取血液灌注信号进一步包括：
使用带通滤波器从电阻抗测量信号中分离特定频率范围的信号，以形成血液灌注信号。
7.如权利要求3所述的方法，其中，所述的从三维差分图像中提取电阻抗测量信号中的血液灌注信号所反映的三维血液灌注图像进一步包括：
利用三维差分图像中像素的时频特性提取三维血液灌注图像。
8.如权利要求7所述的方法，其中，所述的利用三维差分图像中像素的时频特性提取三维血液灌注图像进一步包括：
对三维差分图像中每个像素的时域信号，使用带通滤波器从中分离出特定频率范围的信号，以形成三维血液灌注图像相应像素的时域信号。
9.一种基于电阻抗成像的三维血液灌注图像产生装置，包括：
分布在三维空间内的电极阵列，用于对待测人体区域进行电阻抗测量，得到电阻抗测量信号；
图像重建处理器，其执行存储于存储器中的程序，以便使得：基于电阻抗测量信号中的血液灌注信号，通过图像重建算法重建三维血液灌注图像。
10.如权利要求9所述的装置，其中，所述电极阵列设置在一条或多条阻抗带上、电极背心上或电极头罩上，使得电极呈三维分布。
CN201911424085.5A 2019-12-31 2019-12-31 基于电阻抗成像的三维血液灌注图像产生方法与装置 Active CN111067521B (zh)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CN201911424085.5A CN111067521B (zh) 2019-12-31 2019-12-31 基于电阻抗成像的三维血液灌注图像产生方法与装置
EP20909510.8A EP4085828A4 (en) 2019-12-31 2020-08-04 METHOD AND DEVICE BASED ON ELECTRICAL IMPEDANCE TOMOGRAPHY FOR GENERATING A THREE-DIMENSIONAL IMAGE OF BLOOD PERFUSION
US17/790,414 US20230045401A1 (en) 2019-12-31 2020-08-04 Electrical impedance tomography based method and device for generating three-dimensional blood perfusion image
PCT/CN2020/106754 WO2021135211A1 (zh) 2019-12-31 2020-08-04 基于电阻抗成像的三维血液灌注图像产生方法与装置
JP2022538885A JP7430303B2 (ja) 2019-12-31 2020-08-04 電気インピーダンス・トモグラフィーによる３次元血液灌流画像生成方法及び装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911424085.5A CN111067521B (zh) 2019-12-31 2019-12-31 基于电阻抗成像的三维血液灌注图像产生方法与装置

Publications (2)

Publication Number Publication Date
CN111067521A true CN111067521A (zh) 2020-04-28
CN111067521B CN111067521B (zh) 2021-10-01

Family

ID=70321461

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911424085.5A Active CN111067521B (zh) 2019-12-31 2019-12-31 基于电阻抗成像的三维血液灌注图像产生方法与装置

Country Status (5)

Country Link
US (1) US20230045401A1 (zh)
EP (1) EP4085828A4 (zh)
JP (1) JP7430303B2 (zh)
CN (1) CN111067521B (zh)
WO (1) WO2021135211A1 (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112057073A (zh) * 2020-09-08 2020-12-11 北京华睿博视医学影像技术有限公司 呼吸及血液灌注图像同步生成方法、设备和系统
WO2021135211A1 (zh) * 2019-12-31 2021-07-08 北京华睿博视医学影像技术有限公司 基于电阻抗成像的三维血液灌注图像产生方法与装置
CN113208578A (zh) * 2021-05-26 2021-08-06 北京华睿博视医学影像技术有限公司 血液灌注信号分离模型训练方法及设备
CN113749636A (zh) * 2021-01-26 2021-12-07 北京华睿博视医学影像技术有限公司 一种三维通气图像产生方法、控制器及装置
CN114081468A (zh) * 2022-01-19 2022-02-25 北京华睿博视医学影像技术有限公司 具有准周期性凹凸持液结构的电极片及电极带
WO2022083258A1 (zh) * 2020-10-23 2022-04-28 北京华睿博视医学影像技术有限公司 电阻抗成像方法、系统、存储介质以及电子设备

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115040100B (zh) * 2022-06-14 2023-10-27 安影科技（北京）有限公司 一种视神经血流灌注数值快速采集方法
CN116524123B (zh) * 2023-04-20 2024-02-13 深圳市元甪科技有限公司 一种三维电阻抗断层扫描图像重建方法及相关设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060004301A1 (en) * 2004-06-24 2006-01-05 Kasevich Raymond S Clinical application of electrical impedance tomography to characterize tissue
CN101658428A (zh) * 2008-08-28 2010-03-03 西门子（中国）有限公司 灌注成像技术中对位图进行处理的方法和系统
CN101808576A (zh) * 2007-09-27 2010-08-18 皇家飞利浦电子股份有限公司 用于对血液灌注进行成像的方法
CN103381096A (zh) * 2013-04-19 2013-11-06 西安交通大学 骨表微血管血流灌注分离检测与成像方法
DE102017006107A1 (de) * 2017-06-28 2019-01-03 Drägerwerk AG & Co. KGaA Vorrichtung und Verfahren zur Verarbeitung und Visualisierung von mittels eines Elektro-lmpedanz-Tomographie-Gerätes (EIT) gewonnenen Daten hinsichtlich eines Durchblutungszustandes von Herz und Lunge
CN109864712A (zh) * 2019-04-02 2019-06-11 北京华睿博视医学影像技术有限公司 电阻抗成像设备和方法
CN109965877A (zh) * 2019-04-08 2019-07-05 中国人民解放军陆军军医大学第二附属医院 基于腹部电阻抗频谱参数的便携式组织血流灌注监护仪

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2593006B1 (en) * 2010-07-16 2013-12-04 Csem Sa Method and apparatus for the non-invasive measurement of pulse transit times (ptt)
CN102688041B (zh) * 2012-06-08 2014-02-26 思澜科技（成都）有限公司 基于电极采用交错位方式排列的三维电阻抗断层成像方法
DE102012214786A1 (de) * 2012-08-20 2014-05-15 Dräger Medical GmbH Vorrichtung zur Bestimmung der regionalen Verteilung eines Maßes für die Lungenperfusion
TWI542324B (zh) * 2015-05-18 2016-07-21 國立交通大學 三維電阻抗斷層攝影方法
US11266467B2 (en) * 2016-10-25 2022-03-08 Navix International Limited Systems and methods for registration of intra-body electrical readings with a pre-acquired three dimensional image
CN111067521B (zh) * 2019-12-31 2021-10-01 北京华睿博视医学影像技术有限公司 基于电阻抗成像的三维血液灌注图像产生方法与装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060004301A1 (en) * 2004-06-24 2006-01-05 Kasevich Raymond S Clinical application of electrical impedance tomography to characterize tissue
CN101808576A (zh) * 2007-09-27 2010-08-18 皇家飞利浦电子股份有限公司 用于对血液灌注进行成像的方法
CN101658428A (zh) * 2008-08-28 2010-03-03 西门子（中国）有限公司 灌注成像技术中对位图进行处理的方法和系统
CN103381096A (zh) * 2013-04-19 2013-11-06 西安交通大学 骨表微血管血流灌注分离检测与成像方法
DE102017006107A1 (de) * 2017-06-28 2019-01-03 Drägerwerk AG & Co. KGaA Vorrichtung und Verfahren zur Verarbeitung und Visualisierung von mittels eines Elektro-lmpedanz-Tomographie-Gerätes (EIT) gewonnenen Daten hinsichtlich eines Durchblutungszustandes von Herz und Lunge
CN109864712A (zh) * 2019-04-02 2019-06-11 北京华睿博视医学影像技术有限公司 电阻抗成像设备和方法
CN109965877A (zh) * 2019-04-08 2019-07-05 中国人民解放军陆军军医大学第二附属医院 基于腹部电阻抗频谱参数的便携式组织血流灌注监护仪

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
曲志华等: "机械通气下肺血流灌注状况的电阻抗断层成像评估新方法研究", 《中国医疗设备》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021135211A1 (zh) * 2019-12-31 2021-07-08 北京华睿博视医学影像技术有限公司 基于电阻抗成像的三维血液灌注图像产生方法与装置
CN112057073A (zh) * 2020-09-08 2020-12-11 北京华睿博视医学影像技术有限公司 呼吸及血液灌注图像同步生成方法、设备和系统
WO2022083258A1 (zh) * 2020-10-23 2022-04-28 北京华睿博视医学影像技术有限公司 电阻抗成像方法、系统、存储介质以及电子设备
JP7440985B2 (ja) 2020-10-23 2024-02-29 北京華睿博視医学影像技術有限公司 電気インピーダンスイメージング方法、システム、記憶媒体および電子機器
US11925449B2 (en) 2020-10-23 2024-03-12 Beijing Huarui Boshi Medical Imaging Technology Co., Ltd. Electrical impedance imaging method, system, storage medium, and electronic device
CN113749636A (zh) * 2021-01-26 2021-12-07 北京华睿博视医学影像技术有限公司 一种三维通气图像产生方法、控制器及装置
WO2022160899A1 (zh) * 2021-01-26 2022-08-04 北京华睿博视医学影像技术有限公司 一种三维通气图像产生方法、控制器及装置
CN113208578A (zh) * 2021-05-26 2021-08-06 北京华睿博视医学影像技术有限公司 血液灌注信号分离模型训练方法及设备
CN113208578B (zh) * 2021-05-26 2022-04-26 北京华睿博视医学影像技术有限公司 血液灌注信号分离模型训练方法及设备
CN114081468A (zh) * 2022-01-19 2022-02-25 北京华睿博视医学影像技术有限公司 具有准周期性凹凸持液结构的电极片及电极带

Also Published As

Publication number Publication date
EP4085828A1 (en) 2022-11-09
JP2023508167A (ja) 2023-03-01
EP4085828A4 (en) 2024-01-24
US20230045401A1 (en) 2023-02-09
CN111067521B (zh) 2021-10-01
WO2021135211A1 (zh) 2021-07-08
JP7430303B2 (ja) 2024-02-13

Similar Documents

Publication Publication Date Title
CN111067521B (zh) 基于电阻抗成像的三维血液灌注图像产生方法与装置
JP7244136B2 (ja) 電気インピーダンス断層撮影装置及び方法
US20210100466A1 (en) Channel integrity detection and reconstruction of electrophysiological signals
US9581627B2 (en) Method and system for tomographic imaging
US10555682B2 (en) Signal processing methods
JP2012179352A (ja) 電流双極子を構築するシステムおよび電流双極子を構築する方法
CN103083016B (zh) 用于在软场层析成像中进行波束形成的系统和方法
Rasch et al. Dynamic MRI reconstruction from undersampled data with an anatomical prescan
CN114270397A (zh) 使用电特性断层成像确定流体和组织体积估计的系统和方法
WO2024007877A1 (zh) 用于呈现呼吸阻抗与呼吸阻抗变化率关系的方法及装置
CN116824048B (zh) 一种传感器、雅可比矩阵求解方法、三维成像系统及方法
JP4414591B2 (ja) 脳および心臓の一次電流の断層撮影システムおよび方法
CN113749636B (zh) 一种三维通气图像产生方法、控制器及装置
CN110393522A (zh) 一种基于图总变分约束的无创心脏电生理反演方法
Rajagopal et al. Nonlinear electrocardiographic imaging using polynomial approximation networks
US20120268128A1 (en) Method for Principal Frequency Magnetic Resonance Elastography Inversion
Hartinger et al. Eit system and reconstruction algorithm adapted for skin cancer imaging
US11779310B2 (en) Systems and methods for magnetic resonance elastography with unconstrained optimization inversion
Abdi et al. EIT in breast cancer imaging: application to patient-specific forward model
Rué-Queralt et al. Connectome spectrum electromagnetic tomography: a method to reconstruct electrical brain source-networks at high-spatial resolution
Ibrahim et al. Harmonic phase versus sine-wave modulation for measuring regional heart function from tagged MRI images
CN111615355B (zh) 通道完整性检测和电生理信号的重构
Cortés et al. Electrical Impedance Tomography Simulation for Detection of Breast Tumors Based on Tumor Emulators
Islam Kalman Algorithm Based Electrical Impedance Tomography Imaging
KR20230078217A (ko) 자기 공명 영상으로부터 뇌 전기전도도 지도를 재구성하는 방법 및 장치

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231107

Address after: 102629 213, floor 2, building 10, No. 9, Tianfu Road, Daxing District, Beijing

Patentee after: BEIJING HUARUI BOSHI MEDICAL IMAGE TECHNOLOGY Co.,Ltd.

Patentee after: TSINGHUA University

Address before: Room 431-5, building 3, No. 538, yongfengtun, Haidian District, Beijing 100094

Patentee before: BEIJING HUARUI BOSHI MEDICAL IMAGE TECHNOLOGY Co.,Ltd.