CN113749636A - 一种三维通气图像产生方法、控制器及装置 - Google Patents

一种三维通气图像产生方法、控制器及装置 Download PDF

Info

Publication number
CN113749636A
CN113749636A CN202110111098.8A CN202110111098A CN113749636A CN 113749636 A CN113749636 A CN 113749636A CN 202110111098 A CN202110111098 A CN 202110111098A CN 113749636 A CN113749636 A CN 113749636A
Authority
CN
China
Prior art keywords
ventilation
dimensional
image
electrical impedance
measured
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
CN202110111098.8A
Other languages
English (en)
Other versions
CN113749636B (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 CN202110111098.8A priority Critical patent/CN113749636B/zh
Priority to EP21922489.6A priority patent/EP4285816A1/en
Priority to JP2023544419A priority patent/JP2024505852A/ja
Priority to US18/268,891 priority patent/US20240057887A1/en
Priority to PCT/CN2021/132763 priority patent/WO2022160899A1/zh
Publication of CN113749636A publication Critical patent/CN113749636A/zh
Application granted granted Critical
Publication of CN113749636B publication Critical patent/CN113749636B/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/08—Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/0809—Detecting, measuring or recording devices for evaluating the respiratory organs by impedance pneumography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/085—Measuring impedance of respiratory organs or lung elasticity
    • 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
    • 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/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • 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/7271—Specific aspects of physiological measurement analysis
    • A61B5/7278—Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00—Medical imaging apparatus involving image processing or analysis
    • A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for 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/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

Landscapes

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

Abstract

本发明公开了一种三维通气图像产生方法、控制器及装置，所述方法包括：通过信号提取算法和图像重建算法，根据对待测目标区域进行电阻抗测量得到的电阻抗信号，生成三维通气图像，其中，对待测目标区域进行电阻抗测量利用在待测目标区域的外围呈三维分布的电极阵列实现。本发明提供了能够反映人体胸腔内由于人体通气引起的电阻抗变化的三维通气图像，从而反映人体胸腔在三维空间中各个体积内的通气情况。

Description

一种三维通气图像产生方法、控制器及装置
技术领域
本发明属于电阻抗成像应用技术领域，具体涉及一种三维通气图像产生方法、控制器及装置。
背景技术
EIT(Electrical Impedance Tomography，电阻抗成像)技术是一种无创的、以人体或其他生物体内部的电阻率分布为目标的重建体内组织图像的技术。人体是一个大的生物电导体，各组织、器官均有一定的阻抗，当人体的局部器官发生病变时，局部部位的阻抗必然与其他部位不同，因而可以通过阻抗的测量来对人体器官的病变进行诊断。
现有技术只能生成二维的通气图像，这个二维的图像反映的是待测人体胸腔区域某个断面内由于气体含量变化引起的电阻抗变化。然而，二维的图像难以反映人体胸腔在三维空间中某个体积内的通气情况。
现在亟须一种三维通气图像产生方法、控制器及装置。
发明内容
本发明所要解决的技术问题是如何生成三维通气图像，从而反映人体胸腔在三维空间中各个体积内的通气情况。
针对上述问题，本发明提供了一种三维通气图像产生方法、控制器及装置。
第一方面，本发明提供了一种三维通气图像产生方法，包括以下步骤：
通过信号提取算法和图像重建算法，根据对待测目标区域进行电阻抗测量得到的电阻抗信号，生成三维通气图像，其中，对待测目标区域进行电阻抗测量利用在待测目标区域的外围呈三维分布的电极阵列实现。
在本发明的一些实施例中，通过信号提取算法和图像重建算法，根据对待测目标区域进行电阻抗测量得到的电阻抗信号，生成三维通气图像，包括以下步骤：
通过信号提取算法从对待测目标区域进行电阻抗测量得到的电阻抗信号中提取通气相关信号；
通过图像重建算法根据所述通气相关信号重建三维通气图像。
在本发明的一些实施例中，所述电阻抗信号包括通气相关信号和血液灌注相关信号，通过信号提取算法从对待测目标区域进行电阻抗测量得到的电阻抗信号中提取通气相关信号，包括以下步骤：
利用低通滤波器从对待测目标区域进行电阻抗测量得到的电阻抗信号中提取通气相关信号，其中，所述低通滤波器的截止频率大于通气相关信号的二次谐波频率，且小于血液灌注相关信号的基波频率。
在本发明的一些实施例中，通过信号提取算法和图像重建算法，根据对待测目标区域进行电阻抗测量得到的电阻抗信号，生成三维通气图像，包括以下步骤：
通过图像重建算法根据对待测目标区域进行电阻抗测量得到的电阻抗信号重建三维图像；
通过信号提取算法从所述三维图像中提取三维通气图像。
在本发明的一些实施例中，通过信号提取算法从所述三维图像中提取三维通气图像，包括以下步骤：
根据多个时刻下的三维图像数据列出三维图像中每个像素的时间序列，其中，每个像素的时间序列由每个像素在不同时刻下的值组成；
从三维图像中每个像素的时间序列中提取通气相关像素的时间序列；
根据通气相关像素的时间序列构建三维通气图像。
在本发明的一些实施例中，从三维图像中每个像素的时间序列中提取通气相关像素的时间序列通过频域滤波法、主成分分析法和神经网络法中的任意一种算法实现。
在本发明的一些实施例中，所述信号提取算法为频域滤波法、主成分分析法和神经网络法中的任意一种算法。
在本发明的一些实施例中，所述图像重建算法为线性差分重建算法或基于神经网络的图像重建算法。
第二方面，本发明提供了一种三维通气图像产生控制器，其包括存储器和处理器，该存储器上存储有计算机程序，该计算机程序被处理器执行时实现上述方法的步骤。
第三方面，本发明提供了一种三维通气图像产生装置，包括：
在待测目标区域的外围呈三维分布的电极阵列，其用于对待测目标区域进行电阻抗测量，并将测量得到的电阻抗发送至三维通气图像产生控制器；和
上述三维通气图像产生控制器。
与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
应用本发明的三维通气图像产生方法，通过信号提取算法和图像重建算法，根据对待测目标区域进行电阻抗测量得到的电阻抗信号，生成三维通气图像，其中，对待测目标区域进行电阻抗测量利用在待测目标区域的外围呈三维分布的电极阵列实现，能够提供三维通气图像，从而反映人体胸腔在三维空间中各个体积内的通气情况。
本发明的其它特征和优点将在随后的说明书中阐述，并且部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例共同用于解释本发明，并不构成对本发明的限制。在附图中：
图1示出了本发明实施例一三维通气图像产生方法的流程图；
图2示出了本发明实施例一三维通气图像产生方法的另一流程图；
图3(a)示出了本发明实施例二三维通气图像产生方法的流程图；
图3(b)示出了本发明实施例二三维通气图像产生方法的另一流程图；
图4(a)示出了本发明实施例二人体胸腔测量数据的时域信号示意图；
图4(b)示出了本发明实施例二人体胸腔测量数据的频域信号示意图；
图5(a)示出了本发明实施例二人体胸腔测量数据滤波后的通气相关信号的时域信号示意图；
图5(b)示出了本发明实施例二人体胸腔测量数据滤波后的通气相关信号的频域信号示意图；
图6示出了本发明实施例二利用图3(a)中示出的三维通气图像产生方法生成的人体胸腔三维通气图像示意图；
图7示出了本发明实施例二利用图3(b)中示出的三维通气图像产生方法生成的人体胸腔三维差分图像示意图；
图8(a)示出了图7中示例像素点的时域信号示意图；
图8(b)示出了图7中示例像素点的频域信号示意图；
图9(a)示出了图7中示例像素点数据滤波后的通气相关信号的时域信号示意图；
图9(b)示出了图7中示例像素点数据滤波后的通气相关信号的频域信号示意图；
图10示出了本发明实施例二利用图3(b)中示出的三维通气图像产生方法生成的人体胸腔三维通气图像示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。需要说明的是，只要不构成冲突，本发明中的各个实施例以及各实施例中的各个特征可以相互结合，所形成的技术方案均在本发明的保护范围之内。
实施例一
为解决现有技术中存在的上述技术问题，本发明实施例提供了一种三维通气图像产生方法，其中，本实施例的三维通气图像产生方法通过两种方式实现，具体如图1和图2所示。
参照图1，本实施例的三维通气图像产生方法，包括以下步骤：
S110，通过信号提取算法从对待测目标区域进行电阻抗测量得到的电阻抗信号中提取通气相关信号，其中，对待测目标区域进行电阻抗测量利用在待测目标区域的外围呈三维分布的电极阵列实现，电极阵列可以采用多条阻抗带或者是采用电极呈三维分布的电极背心；
S120，通过图像重建算法根据所述通气相关信号重建三维通气图像。
在一个具体实施方式中，所述电阻抗信号包括通气相关信号和血液灌注相关信号，通过信号提取算法从对待测目标区域进行电阻抗测量得到的电阻抗信号中提取通气相关信号，包括以下步骤：
利用低通滤波器从对待测目标区域进行电阻抗测量得到的电阻抗信号中提取通气相关信号，其中，所述低通滤波器的截止频率大于通气相关信号的二次谐波频率，且小于血液灌注相关信号的基波频率。
在步骤S110中，所述信号提取算法为频域滤波法、主成分分析法和神经网络法中的任意一种算法。
在步骤S120中，所述图像重建算法为线性差分重建算法或基于神经网络的图像重建算法。
参照图2，本实施例的三维通气图像产生方法，包括以下步骤：
S210，通过图像重建算法根据对待测目标区域进行电阻抗测量得到的电阻抗信号重建三维图像，其中，对待测目标区域进行电阻抗测量利用在待测目标区域的外围呈三维分布的电极阵列实现；
S220，根据多个时刻下的三维图像数据列出三维图像中每个像素的时间序列，其中，每个像素的时间序列由每个像素在不同时刻下的值组成；
S230，从三维图像中每个像素的时间序列中提取通气相关像素的时间序列；
S240，根据通气相关像素的时间序列构建三维通气图像。
在步骤S230中，从三维图像中每个像素的时间序列中提取通气相关像素的时间序列通过频域滤波法、主成分分析法和神经网络法中的任意一种算法实现。
在步骤S210中，所述图像重建算法为线性差分重建算法或基于神经网络的图像重建算法。
实施例二
为解决现有技术中存在的上述技术问题，本发明实施例基于实施例一提供了一种应用于人体胸腔的三维通气图像产生方法，其中，本实施例的三维通气图像产生方法通过两种方式实现，具体如图3(a)和图3(b)所示。
如图3(a)所示，本实施例的三维通气图像产生方法包括以下步骤：首先，对待测人体胸腔区域进行电阻抗测量；然后，从测量信号中提取通气相关信号；最后，重建三维通气图像。具体过程如下：
第一步，对待测人体胸腔区域进行电阻抗测量。在所述电阻抗测量中，首先，需要在待测人体胸腔周围固定电极阵列。所述电极阵列包含若干个分布在三维空间内的电极。然后，通过电极阵列对待测人体胸腔进行激励并测量由此产生的响应，即：轮流对电极施加电流激励，并依次在其他电极上测量由此产生的电压信号；
第二步，从上一步测量得到的电阻抗信号中提取通气相关信号。在本步骤的一个实施例中，利用滤波器从所测量的电阻抗信号中提取通气相关信号。滤波器可以使用有限脉冲响应滤波器或无限脉冲响应滤波器等。下面是对人体胸腔进行测量的实施例。图4(a)示出了测量数据的时域信号。图中曲线代表特定电极激励时，在特定电极上测量得到的电压信号。其他激励-测量情况得到的数据与之类似。需要说明的是，图中纵坐标为从数字电压表直接读取的数值，仍未将其转换为电压值。图4(b)示出了测量数据的频域信号。图4(b)所示信号由图4(a)所示信号经过傅里叶变换得到。由图4(b)可以分辨出通气相关信号和血液灌注相关信号。为了提取通气相关信号，设计一个低通滤波器，可以是有限脉冲响应低通数字滤波器，该滤波器的截止频率大于通气相关信号的二次谐波频率，且小于血液灌注相关信号的基波频率。滤波后的信号图形如图5(a)和图5(b)所示，其中，图5(a)为时域信号图形，图5(b)为频域信号图形。
在本步骤的另一个实施例中，使用基于PCA(Principle Component Analysis，主成分分析)的方法提取通气相关信号。具体地，假设测量信号为u。其尺寸为Nt×Nc，其中，Nt为采样点数，Nc为特征数(在这里为测量通道数)。利用主成分分析得到信号的主成分
Figure BDA0002916831200000061
其中，pi(i＝1,2,…,Nc)的尺寸为Nt×1，且其对应的特征值依次减小。将前若干个主成分(如p1，p2)作为模板对信号u进行模板匹配滤波，得到通气相关信号uV。
在本步骤的另一个实施例中，使用基于神经网络的方法提取通气相关信号。具体地，所述基于神经网络的方法分为训练和预测两个阶段。在训练阶段，利用训练数据通过有监督或无监督的方法训练一个通气相关信号提取网络；在预测阶段，利用训练好的通气相关信号提取网络来提取电阻抗测量信号中的通气相关信号。
第三步，利用第二步所提取的通气相关信号，通过图像重建算法重建三维通气图像。所述三维通气图像反映由于呼吸引起的待测人体区域内的电阻抗变化。在本步骤的一个实施例中，所述图像重建算法是线性差分重建算法。下面是线性差分重建算法进行重建三维通气图像的实施例。
假设第二步所提取出的通气相关信号的时域形式为u(t)，其中，t为时间变量。EIT差分重建可以表述为如下最小二乘问题：
minδσ||J·δσ-δu||2+α||R·δσ||2，
其中，J为雅可比矩阵，δu＝u(t)-u(tref)为信号在时刻t相对于参考时刻tref的变化，δσ为上述两个时刻人体内由于通气引起的电导率变化，R为正则化矩阵，α为正则化参数。参考时刻tref既可以设置为在整个图像重建过程中固定不变，也可以设置为随着图像重建过程的进行而动态地更新。δσ定义在离散化的三维模型中，如四面体网格或体素网格。上述问题的解为
δσ*＝(JT·J+αRT·R)-1·JT·δu.
令D＝(JT·J+αRT·R)-1·JT，则上述公式可以表示为：
δσ*＝D·δu.
上述δσ*即为计算所得的三维通气图像。
图6示出了利用上述方法所产生的人体胸腔的三维通气图像的一个示意图。
在本步骤的另一个实施例中，所述图像重建算法为基于机器学习的方法。EIT差分成像可以表示为：
Figure BDA0002916831200000071
其中，
Figure BDA0002916831200000072
为重建算子，δu为不同时刻测量数据的变化，δσ为相应时刻电导率的变化。所述基于机器学习的方法分为训练和预测两个阶段。首先，在训练阶段，给定训练数据{δui,δσi}，可以训练一个机器学习模型
Figure BDA0002916831200000073
来近似算子
Figure BDA0002916831200000074
在预测阶段，给定差分测量信号δu，可以通过
Figure BDA0002916831200000076
来预测相应的电导率变化：
Figure BDA0002916831200000075
除了上述实施例中的图像重建算法外，本步骤还可以使用各种线性的或非线性的、迭代的或非迭代的、随机的或确定性的图像重建算法。
如图3(b)所示，本实施例的三维通气图像产生方法包括以下步骤：首先，通过对待测人体胸腔区域进行电阻抗测量，然后重建三维差分图像，最后从三维差分图像中提取三维通气图像。具体过程如下：
第一步，对待测人体胸腔区域进行电阻抗测量。
第二步，利用上一步测量得到的电阻抗信号，通过图像重建算法重建三维差分图像。所述三维差分图像反映了待测人体胸腔内的电阻抗变化，该电阻抗变化可能是由人体通气或血液灌注引起的。所述图像重建算法可以采用上述图像重建算法。图7示出了利用图4所示数据和线性差分重建算法生成的三维差分图像。
第三步，从上一步得到的三维差分图像中提取通气图像。在本步骤的一个实施例中，利用滤波器从三维差分图像中提取通气图像。假设N个时刻的三维差分图像可排列为一个矩阵A＝{a1,a2,…,aM}T，其中，ai(i＝1,2,…,M)为像素i在N个时刻的值所组成的列向量，M为三维图像中像素的总数。对每个像素i的时间序列ai(i＝1,2,…,M)进行低通滤波可以得到通气图像上对应像素的时间序列。具体地，假设滤波函数为f(·)，则通气图像为AV＝{f(a1),f(a2),…,f(aM)}T。
图8(a)和其对应的频谱图8(b)示出了图7中人体胸腔三维差分图像示例像素点的时间序列。由图8(b)可以分辨出通气相关信号和血液灌注相关信号。为了提取通气相关信号，设计一个低通滤波器，可以是有限脉冲响应低通数字滤波器，该滤波器的截止频率大于通气相关信号的二次谐波频率，且小于血液灌注相关信号的基波频率。
图9(a)和频谱图9(b)示出了对图7中示例像素点进行滤波后的时域信号。对三维差分图像中的每个像素进行上述低通滤波后，即可得到三维通气图像，如图10所示。在本步骤的另外两个实施例中，可以使用基于主成分分析和基于神经网络的方法提取通气图像。
需要说明的是，图8(a)、图8(b)、图9(a)和图9(b)中的A.U.为任意单位。
本实施例的应用于人体胸腔的三维通气图像产生方法提供了能够反映人体胸腔内由于人体通气引起的电阻抗变化的三维通气图像，从而反映人体胸腔在三维空间中各个体积内的通气情况。
实施例三
为解决现有技术中存在的上述技术问题，本发明实施例提供了一种三维通气图像产生控制器。
本实施例的三维通气图像产生控制器.，其包括存储器和处理器，该存储器上存储有计算机程序，该计算机程序被处理器执行时实现实施例一和实施例二所述的方法的步骤。
实施例四
为解决现有技术中存在的上述技术问题，本发明实施例还提供了一种三维通气图像产生装置。
本实施例的三维通气图像产生装置，包括：
在待测目标区域的外围呈三维分布的电极阵列，其用于对待测目标区域进行电阻抗测量，并将测量得到的电阻抗发送至三维通气图像产生控制器；
实施例三所述的三维通气图像产生控制器。
本实施例的三维通气图像产生装置，还包括：图像显示装置，其用于对所述三维通气图像产生控制器产生的三维通气图像并进行显示。
虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (10)

1.一种三维通气图像产生方法，其特征在于，包括以下步骤：
通过信号提取算法和图像重建算法，根据对待测目标区域进行电阻抗测量得到的电阻抗信号，生成三维通气图像，其中，对待测目标区域进行电阻抗测量利用在待测目标区域的外围呈三维分布的电极阵列实现。
2.根据权利要求1所述的方法，其特征在于，通过信号提取算法和图像重建算法，根据对待测目标区域进行电阻抗测量得到的电阻抗信号，生成三维通气图像，包括以下步骤：
通过信号提取算法从对待测目标区域进行电阻抗测量得到的电阻抗信号中提取通气相关信号；
通过图像重建算法根据所述通气相关信号重建三维通气图像。
3.根据权利要求2所述的方法，其特征在于，所述电阻抗信号包括通气相关信号和血液灌注相关信号，通过信号提取算法从对待测目标区域进行电阻抗测量得到的电阻抗信号中提取通气相关信号，包括以下步骤：
利用低通滤波器从对待测目标区域进行电阻抗测量得到的电阻抗信号中提取通气相关信号，其中，所述低通滤波器的截止频率大于通气相关信号的二次谐波频率，且小于血液灌注相关信号的基波频率。
4.根据权利要求1所述的方法，其特征在于，通过信号提取算法和图像重建算法，根据对待测目标区域进行电阻抗测量得到的电阻抗信号，生成三维通气图像，包括以下步骤：
通过图像重建算法根据对待测目标区域进行电阻抗测量得到的电阻抗信号重建三维图像；
通过信号提取算法从所述三维图像中提取三维通气图像。
5.根据权利要求4所述的方法，其特征在于，通过信号提取算法从所述三维图像中提取三维通气图像，包括以下步骤：
根据多个时刻下的三维图像数据列出三维图像中每个像素的时间序列，其中，每个像素的时间序列由每个像素在不同时刻下的值组成；
从三维图像中每个像素的时间序列中提取通气相关像素的时间序列；
根据通气相关像素的时间序列构建三维通气图像。
6.根据权利要求5所述的方法，其特征在于，从三维图像中每个像素的时间序列中提取通气相关像素的时间序列通过频域滤波法、主成分分析法和神经网络法中的任意一种算法实现。
7.根据权利要求2或4所述的方法，其特征在于，所述信号提取算法为频域滤波法、主成分分析法和神经网络法中的任意一种算法。
8.根据权利要求2或4所述的方法，其特征在于，所述图像重建算法为线性差分重建算法或基于神经网络的图像重建算法。
9.一种三维通气图像产生控制器，其包括存储器和处理器，其特征在于，该存储器上存储有计算机程序，该计算机程序被处理器执行时实现如权利要求1至8中任一项所述方法的步骤。
10.一种三维通气图像产生装置，其特征在于，包括：
在待测目标区域的外围呈三维分布的电极阵列，其用于对待测目标区域进行电阻抗测量，并将测量得到的电阻抗发送至三维通气图像产生控制器；和
根据权利要求9所述的三维通气图像产生控制器。
CN202110111098.8A 2021-01-26 2021-01-26 一种三维通气图像产生方法、控制器及装置 Active CN113749636B (zh)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CN202110111098.8A CN113749636B (zh) 2021-01-26 2021-01-26 一种三维通气图像产生方法、控制器及装置
EP21922489.6A EP4285816A1 (en) 2021-01-26 2021-11-24 Three-dimensional ventilation image generation method, and controller and apparatus
JP2023544419A JP2024505852A (ja) 2021-01-26 2021-11-24 三次元換気画像生成方法、コントローラおよび装置
US18/268,891 US20240057887A1 (en) 2021-01-26 2021-11-24 Three-dimensional ventilation image generation method, and controller and apparatus
PCT/CN2021/132763 WO2022160899A1 (zh) 2021-01-26 2021-11-24 一种三维通气图像产生方法、控制器及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110111098.8A CN113749636B (zh) 2021-01-26 2021-01-26 一种三维通气图像产生方法、控制器及装置

Publications (2)

Publication Number Publication Date
CN113749636A true CN113749636A (zh) 2021-12-07
CN113749636B CN113749636B (zh) 2022-06-24

Family

ID=78786469

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110111098.8A Active CN113749636B (zh) 2021-01-26 2021-01-26 一种三维通气图像产生方法、控制器及装置

Country Status (5)

Country Link
US (1) US20240057887A1 (zh)
EP (1) EP4285816A1 (zh)
JP (1) JP2024505852A (zh)
CN (1) CN113749636B (zh)
WO (1) WO2022160899A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115153500A (zh) * 2022-07-04 2022-10-11 北京华睿博视医学影像技术有限公司 用于呈现呼吸阻抗与呼吸阻抗变化率关系的方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102793551A (zh) * 2011-05-24 2012-11-28 柯尼卡美能达医疗印刷器材株式会社 胸部诊断辅助信息生成系统
CN104853677A (zh) * 2012-12-12 2015-08-19 柯尼卡美能达株式会社 图像处理装置以及程序
CN109864712A (zh) * 2019-04-02 2019-06-11 北京华睿博视医学影像技术有限公司 电阻抗成像设备和方法
CN111053556A (zh) * 2019-12-26 2020-04-24 北京华睿博视医学影像技术有限公司 基于监督下降法的电阻抗成像方法与装置
CN111067521A (zh) * 2019-12-31 2020-04-28 北京华睿博视医学影像技术有限公司 基于电阻抗成像的三维血液灌注图像产生方法与装置
US20200221970A1 (en) * 2017-06-28 2020-07-16 Drägerwerk AG & Co. KGaA Device and method for processing and visualizing data relating to cardiac and pulmonary circulation, obtained by means of an electrical impedance tomography device
CN112057073A (zh) * 2020-09-08 2020-12-11 北京华睿博视医学影像技术有限公司 呼吸及血液灌注图像同步生成方法、设备和系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE9803862L (sv) * 1998-11-11 2000-03-20 Siemens Elema Ab Elektriskt impedanstomografisystem
KR101812587B1 (ko) * 2016-11-18 2018-01-30 주식회사 바이랩 피험자의 영상 모니터링 장치 및 그 방법과, 영상 모니터링 시스템
WO2020007991A1 (en) * 2018-07-04 2020-01-09 Navix International Limited System and method for conductivity-based imaging

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102793551A (zh) * 2011-05-24 2012-11-28 柯尼卡美能达医疗印刷器材株式会社 胸部诊断辅助信息生成系统
CN104853677A (zh) * 2012-12-12 2015-08-19 柯尼卡美能达株式会社 图像处理装置以及程序
US20200221970A1 (en) * 2017-06-28 2020-07-16 Drägerwerk AG & Co. KGaA Device and method for processing and visualizing data relating to cardiac and pulmonary circulation, obtained by means of an electrical impedance tomography device
CN109864712A (zh) * 2019-04-02 2019-06-11 北京华睿博视医学影像技术有限公司 电阻抗成像设备和方法
CN111053556A (zh) * 2019-12-26 2020-04-24 北京华睿博视医学影像技术有限公司 基于监督下降法的电阻抗成像方法与装置
CN111067521A (zh) * 2019-12-31 2020-04-28 北京华睿博视医学影像技术有限公司 基于电阻抗成像的三维血液灌注图像产生方法与装置
CN112057073A (zh) * 2020-09-08 2020-12-11 北京华睿博视医学影像技术有限公司 呼吸及血液灌注图像同步生成方法、设备和系统

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115153500A (zh) * 2022-07-04 2022-10-11 北京华睿博视医学影像技术有限公司 用于呈现呼吸阻抗与呼吸阻抗变化率关系的方法及装置

Also Published As

Publication number Publication date
US20240057887A1 (en) 2024-02-22
JP2024505852A (ja) 2024-02-08
CN113749636B (zh) 2022-06-24
EP4285816A1 (en) 2023-12-06
WO2022160899A1 (zh) 2022-08-04

Similar Documents

Publication Publication Date Title
CN111067521B (zh) 基于电阻抗成像的三维血液灌注图像产生方法与装置
Dusek et al. Electrical impedance tomography methods and algorithms processed with a GPU
Hauptmann et al. Multi-scale learned iterative reconstruction
EP0599877A1 (en) ELECTRIC IMPEDANCE TOMOGRAPHY.
CN110251130B (zh) 基于可移动变形组件的电阻抗成像方法、装置及系统
Gong et al. Higher order total variation regularization for EIT reconstruction
CN114270397A (zh) 使用电特性断层成像确定流体和组织体积估计的系统和方法
CN113749636B (zh) 一种三维通气图像产生方法、控制器及装置
CN112057073A (zh) 呼吸及血液灌注图像同步生成方法、设备和系统
CN114758031A (zh) 电容耦合电阻抗层析成像图像重建方法和装置
JP7300811B2 (ja) 医用情報処理装置、医用情報処理方法、およびプログラム
Nahvi et al. Wideband electrical impedance tomography
CN115797493B (zh) 基于一维系统矩阵稀疏采样的磁场自由线磁粒子成像方法
CN114469046A (zh) 电阻抗成像方法、系统、存储介质以及电子设备
Aggarwal et al. Double temporal sparsity based accelerated reconstruction of compressively sensed resting-state fMRI
Ain et al. Dual modality electrical impedance and ultrasound reflection tomography to improve image quality
CN115830156B (zh) 精准电阻抗层析成像方法、装置、系统、介质及设备
Gong et al. EIT imaging regularization based on spectral graph wavelets
Wang et al. Patch-based sparse reconstruction for electrical impedance tomography
Xie et al. Non‐invasive reconstruction of dynamic myocardial transmembrane potential with graph‐based total variation constraints
Zhang et al. HybridDenseU-Net: learning a multi-scale convolution and dense connectivity CNN for inverse imaging problems
Hashemzadeh et al. A hybrid analytical-numerical algorithm for determining the neuronal current via EEG
Vauhkonen et al. Estimation of organ boundaries in electrical impedance tomography
Pan et al. Targeted-ROI imaging in electron paramagnetic resonance imaging
CN110988043B (zh) 基于多频加权频差的多介质分离成像方法

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: 20231102

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: 102629 213, floor 2, building 10, No. 9, Tianfu Road, Daxing District, Beijing

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