CN111067521B - 基于电阻抗成像的三维血液灌注图像产生方法与装置 - Google Patents
基于电阻抗成像的三维血液灌注图像产生方法与装置 Download PDFInfo
- Publication number
- CN111067521B CN111067521B CN201911424085.5A CN201911424085A CN111067521B CN 111067521 B CN111067521 B CN 111067521B CN 201911424085 A CN201911424085 A CN 201911424085A CN 111067521 B CN111067521 B CN 111067521B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- image
- blood perfusion
- electrical impedance
- signal
- 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.)
- Active
Links
- 230000008081 blood perfusion Effects 0.000 title claims abstract description 106
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000003384 imaging method Methods 0.000 title claims abstract description 23
- 238000002847 impedance measurement Methods 0.000 claims abstract description 52
- 210000000746 body region Anatomy 0.000 claims abstract description 11
- 230000015654 memory Effects 0.000 claims description 10
- 238000003491 array Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 abstract description 6
- 201000010099 disease Diseases 0.000 abstract description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 abstract description 2
- 238000003745 diagnosis Methods 0.000 abstract description 2
- 230000002349 favourable effect Effects 0.000 abstract description 2
- 238000010191 image analysis Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 12
- 230000010412 perfusion Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 210000000038 chest Anatomy 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 210000004072 lung Anatomy 0.000 description 4
- 210000000056 organ Anatomy 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 101100521334 Mus musculus Prom1 gene Proteins 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000002593 electrical impedance tomography Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 210000001015 abdomen Anatomy 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 210000003238 esophagus Anatomy 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 210000003414 extremity Anatomy 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 210000003437 trachea Anatomy 0.000 description 1
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)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (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 (6)
1.一种基于电阻抗成像的三维血液灌注图像产生方法，包括如下步骤：
利用分布在三维空间内的电极阵列，对待测人体区域进行电阻抗测量，得到电阻抗测量信号；
根据电阻抗测量信号，通过图像重建算法重建三维差分图像；以及
从三维差分图像中提取电阻抗测量信号中的血液灌注信号所反映的三维血液灌注图像。
2.如权利要求1所述的方法，其中，所述电极阵列设置在一条或多条阻抗带上、电极背心上或电极头罩上，使得电极呈三维分布。
3.如权利要求1所述的方法，其中，所述的从三维差分图像中提取电阻抗测量信号中的血液灌注信号所反映的三维血液灌注图像进一步包括：
利用三维差分图像中像素的时频特性提取三维血液灌注图像。
4.如权利要求3所述的方法，其中，所述的利用三维差分图像中像素的时频特性提取三维血液灌注图像进一步包括：
对三维差分图像中每个像素的时域信号，使用带通滤波器从中分离出特定频率范围的信号，以形成三维血液灌注图像相应像素的时域信号。
5.一种基于电阻抗成像的三维血液灌注图像产生装置，包括：
分布在三维空间内的电极阵列，用于对待测人体区域进行电阻抗测量，得到电阻抗测量信号；
图像重建处理器，其执行存储于存储器中的程序，以便使得：根据电阻抗测量信号，通过图像重建算法重建三维差分图像；以及从三维差分图像中提取电阻抗测量信号中的血液灌注信号所反映的三维血液灌注图像。
6.如权利要求5所述的装置，其中，所述电极阵列设置在一条或多条阻抗带上、电极背心上或电极头罩上，使得电极呈三维分布。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911424085.5A CN111067521B (zh) | 2019-12-31 | 2019-12-31 | 基于电阻抗成像的三维血液灌注图像产生方法与装置 |
| PCT/CN2020/106754 WO2021135211A1 (zh) | 2019-12-31 | 2020-08-04 | 基于电阻抗成像的三维血液灌注图像产生方法与装置 |
| 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 |
| 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 CN111067521A (zh) | 2020-04-28 |
| CN111067521B true 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) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111067521B (zh) * | 2019-12-31 | 2021-10-01 | 北京华睿博视医学影像技术有限公司 | 基于电阻抗成像的三维血液灌注图像产生方法与装置 |
| CN112057073A (zh) * | 2020-09-08 | 2020-12-11 | 北京华睿博视医学影像技术有限公司 | 呼吸及血液灌注图像同步生成方法、设备和系统 |
| CN114469046A (zh) | 2020-10-23 | 2022-05-13 | 北京华睿博视医学影像技术有限公司 | 电阻抗成像方法、系统、存储介质以及电子设备 |
| CN113749636B (zh) * | 2021-01-26 | 2022-06-24 | 北京华睿博视医学影像技术有限公司 | 一种三维通气图像产生方法、控制器及装置 |
| CN113208578B (zh) * | 2021-05-26 | 2022-04-26 | 北京华睿博视医学影像技术有限公司 | 血液灌注信号分离模型训练方法及设备 |
| CN114081468B (zh) * | 2022-01-19 | 2022-05-06 | 北京华睿博视医学影像技术有限公司 | 具有准周期性凹凸持液结构的电极片及电极带 |
| CN115040100B (zh) * | 2022-06-14 | 2023-10-27 | 安影科技（北京）有限公司 | 一种视神经血流灌注数值快速采集方法 |
| CN116524123B (zh) * | 2023-04-20 | 2024-02-13 | 深圳市元甪科技有限公司 | 一种三维电阻抗断层扫描图像重建方法及相关设备 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 | 北京华睿博视医学影像技术有限公司 | 电阻抗成像设备和方法 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006012181A1 (en) * | 2004-06-24 | 2006-02-02 | Ksn Energies, Llc | Electrical impedance tomography to characterize tissue |
| WO2009040742A2 (en) * | 2007-09-27 | 2009-04-02 | Philips Intellectual Property & Standards Gmbh | Methods for imaging the blood perfusion |
| 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 |
| CN109965877B (zh) * | 2019-04-08 | 2023-03-14 | 中国人民解放军陆军军医大学第二附属医院 | 基于腹部电阻抗频谱参数的便携式组织血流灌注监护仪 |
| CN111067521B (zh) * | 2019-12-31 | 2021-10-01 | 北京华睿博视医学影像技术有限公司 | 基于电阻抗成像的三维血液灌注图像产生方法与装置 |
-
2019
- 2019-12-31 CN CN201911424085.5A patent/CN111067521B/zh active Active
-
2020
- 2020-08-04 WO PCT/CN2020/106754 patent/WO2021135211A1/zh unknown
- 2020-08-04 EP EP20909510.8A patent/EP4085828A4/en active Pending
- 2020-08-04 JP JP2022538885A patent/JP7430303B2/ja active Active
- 2020-08-04 US US17/790,414 patent/US20230045401A1/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 | 北京华睿博视医学影像技术有限公司 | 电阻抗成像设备和方法 |
Non-Patent Citations (1)
| Title |
|---|
| 机械通气下肺血流灌注状况的电阻抗断层成像评估新方法研究;曲志华等;《中国医疗设备》;20190131;第34卷(第1期);第7页第2栏第5段 * |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2023508167A (ja) | 2023-03-01 |
| CN111067521A (zh) | 2020-04-28 |
| US20230045401A1 (en) | 2023-02-09 |
| EP4085828A1 (en) | 2022-11-09 |
| JP7430303B2 (ja) | 2024-02-13 |
| EP4085828A4 (en) | 2024-01-24 |
| WO2021135211A1 (zh) | 2021-07-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111067521B (zh) | 基于电阻抗成像的三维血液灌注图像产生方法与装置 | |
| JP7244136B2 (ja) | 電気インピーダンス断層撮影装置及び方法 | |
| US12059260B2 (en) | Channel integrity detection and reconstruction of electrophysiological signals | |
| US9581627B2 (en) | Method and system for tomographic imaging | |
| KR20200141478A (ko) | 서로 다른 반복 시간을 가지는 두 개의 MRI 이미지들로부터 얻어진 1MHz 미만 저주파 교류 전도도 추산치 | |
| US10555682B2 (en) | Signal processing methods | |
| JP2012179352A (ja) | 電流双極子を構築するシステムおよび電流双極子を構築する方法 | |
| Rasch et al. | Dynamic MRI reconstruction from undersampled data with an anatomical prescan | |
| CN114270397A (zh) | 使用电特性断层成像确定流体和组织体积估计的系统和方法 | |
| CN110811596B (zh) | 基于低秩与稀疏约束和非局部全变分的无创心脏电位重建方法 | |
| CN109758149B (zh) | 一种结合时差成像和频差成像的电阻抗成像方法 | |
| CN106821380A (zh) | 基于乘性正则化的生物医学电阻抗成像方法和装置 | |
| JP4414591B2 (ja) | 脳および心臓の一次電流の断層撮影システムおよび方法 | |
| WO2024007877A1 (zh) | 用于呈现呼吸阻抗与呼吸阻抗变化率关系的方法及装置 | |
| CN116824048B (zh) | 一种传感器、雅可比矩阵求解方法、三维成像系统及方法 | |
| Sze | Detection of breast cancer with electrical impedance mammography | |
| CN113749636B (zh) | 一种三维通气图像产生方法、控制器及装置 | |
| Abdi et al. | EIT in breast cancer imaging: application to patient-specific forward model | |
| Cortés et al. | Electrical Impedance Tomography Simulation for Detection of Breast Tumors Based on Tumor Emulators | |
| Zolgharni et al. | Imaging haemorrhagic cerebral stroke by frequency-difference magnetic induction tomography: numerical modelling | |
| Islam | Kalman Algorithm Based Electrical Impedance Tomography Imaging | |
| KR20230078217A (ko) | 자기 공명 영상으로부터 뇌 전기전도도 지도를 재구성하는 방법 및 장치 | |
| Baldewsing et al. | Aalders, MC, see van der Meer, FJ, T-MI Oct 05 1369-1376 Abbey, CK, see Zemp, RJ, T-MI Mar 05 300-310 Abboud, T., see Kybic, J., T-MI Jan 05 12-28 Abolmaesumi, P., see Yao, W., T-MI Aug 05 997-1010 Abramoff, MD, see Niemeijer, M., T-MI May 05 584-592 |
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. |