CN112401866A - 基于布尔运算的电阻抗成像方法 - Google Patents
基于布尔运算的电阻抗成像方法 Download PDFInfo
- Publication number
- CN112401866A CN112401866A CN202011259573.8A CN202011259573A CN112401866A CN 112401866 A CN112401866 A CN 112401866A CN 202011259573 A CN202011259573 A CN 202011259573A CN 112401866 A CN112401866 A CN 112401866A
- Authority
- CN
- China
- Prior art keywords
- electrical impedance
- spline curve
- boolean operation
- impedance imaging
- imaging
- 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.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 claims abstract description 63
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims description 12
- 239000000654 additive Substances 0.000 claims description 9
- 230000000996 additive effect Effects 0.000 claims description 9
- 238000012804 iterative process Methods 0.000 claims description 5
- 239000000126 substance Substances 0.000 claims description 5
- 238000010561 standard procedure Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 description 22
- 238000002593 electrical impedance tomography Methods 0.000 description 21
- 238000010586 diagram Methods 0.000 description 9
- 230000008901 benefit Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 7
- 206010035664 Pneumonia Diseases 0.000 description 5
- 238000002591 computed tomography Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 210000004072 lung Anatomy 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000002059 diagnostic imaging Methods 0.000 description 3
- 239000003814 drug Substances 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 206010001052 Acute respiratory distress syndrome Diseases 0.000 description 1
- 241000711573 Coronaviridae Species 0.000 description 1
- 206010011409 Cross infection Diseases 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 206010029803 Nosocomial infection Diseases 0.000 description 1
- 208000013616 Respiratory Distress Syndrome Diseases 0.000 description 1
- 201000000028 adult respiratory distress syndrome Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000010420 art technique Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 206010015037 epilepsy Diseases 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000003094 perturbing effect Effects 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
- 210000001835 viscera 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
一种基于布尔运算的电阻抗成像方法，包括以下步骤：建立完备电极模型，通过有限元方法求解得到观测模型；建立基于B样条曲线形式的基本单元形状的矩阵表示形式；对B样条曲线进行布尔运算；基于B样条曲线与布尔运算相结合的形状重构算法；执行重构算法，进行图像重构，以实现基于布尔运算的电阻抗成像。本发明所提出的方法可以重建得到多相、内嵌等复杂形状，得到的重建电阻抗图像成像精度、空间分辨率高，并且对电阻抗成像中普遍存在的模型误差、测量噪音等具有很好的鲁棒性。
Description
技术领域
本发明涉及电阻抗断层成像技术领域，是一种通过布尔运算与电阻抗成像技术相结合来重建对象体内阻抗分布图像的方法。
背景技术
电阻抗成像(Electrical impedance tomography，EIT)是一种无创、无辐射、低成本、能够连续监测的新型功能成像方法，此项医学影像技术是基于人体的各个组织具有不同电阻抗值这一特性，在人体外部施加一个比较小的安全激励电流(电压)，通过放置在体表的电极来测量电压(电流)，从而重建出人体内二维或三维的电阻抗分布图像。由于人体各个器官和组织的电导率不同，因此由EIT技术得到的电导率图像不仅包含了丰富的解剖学的信息，而且能反映出器官和组织电导率相应的生理、病理状态和功能信息。EIT技术在肺通气检测，肿瘤、癫痫等多种疾病的早期发现与诊断等方向有着诱人的应用前景。
自2019年12月起，世界人民面临着新型冠状肺炎的巨大考验，能够早期地发现疑似病人以及对确诊病人进行有效的治疗与监护对控制疫情的发展具有重要的意义。在《新型冠状病毒肺炎诊疗方案》(试行第七版)中提到严重感染者可快速进展为急性呼吸窘迫综合征，EIT可对感染者进行监护，并且在对病人采取肺复张治疗前可利用EIT技术进行可复张性评价。同时可通过EIT影像学监测肺炎感染者病情的发展，判断治疗的效果，最大限度地降低检查时往返于隔离室和放射科之间发生交叉感染的可能。EIT成像技术在新冠肺炎的监护、诊断等方面起着重要作用。
电阻抗技术具有无创、无辐射、功能信息丰富等优点，但其本质是一个病态的非线性、不适定的反问题，目前成像无法达到x射线计算机断层成像(CT)，核磁共振成像(MRI)等技术的成像精度，也由于EIT不适定的本质使得其对于多相、内嵌等复杂介质的重建能力有限，而在医学、工业应用中存在着大量复杂介质的重建情况。因此研究开发高精度、性能稳定的EIT系统及算法，提升其对于多相、内嵌等复杂介质的重建能力，探索其在临床医学以及工业领域中的应用，是当前的热点和难点问题。
在实现本发明的过程中，申请人发现了以下几个现有技术及其存在的缺点。
1、现有技术提出将B样条曲线应用于EIT。此算法将EIT图像重建问题转换为形状重建问题，利用B样条曲线通过调整控制点来逼近目标形状，提升了重建算法对于尖锐特征的保持能力，但是该算法并不具备拓扑演化能力，进行重建前需要获得区域内部物体的个数这一先验信息。对于无法知晓区域内部物体个数的情况，该算法则无法进行重建或无法重建得到可靠的重建结果，这在很大程度上限制这一算法的推广和应用。
2、现有技术中针对每一个目标物体给定多个MMC作为初始迭代值，通过优化组件的几何特征参数、改变组件之间的连通性或分布实现迭代过程形状的改变，并且利用高斯牛顿方法迭代改变电导率分布以及由正问题模型计算得到的电压值，从而通过逆问题最小化目标泛函，最终实现形状重建。但该发明仅针对两相的情况进行重建，对于多相以及内嵌等复杂形状无法进行重建。
3、现有技术中存在一种基于粒子群和正则化高斯-牛顿迭代算法的生物电阻抗成像方法。该方法在正问题计算中采用非均匀剖分，以提高成像精度。并且在生物电阻抗成像问题的求解中，采用粒子群算法产生一个接近真实值的初值作为正则化高斯-牛顿迭代算法的一个初始值，再利用正则化高斯-牛顿迭代算法求解逆问题。该方法采用非均匀剖分以及标准的粒子群法为正则化高斯-牛顿迭代算法产生一个初始值，可以提高成像的精确度并且克服牛顿类算法对初始值敏感的问题，虽然此方法充分利用了先验信息，但在一定程度上不仅增加了计算的复杂度，而且该方法重建的图像的空间分辨率，成像精度仍较低。
发明内容
有鉴于此，本发明的主要目的在于提供一种基于布尔运算的电阻抗成像方法，以期部分地解决上述技术问题中的至少之一。
为了实现上述目的，作为本发明的一方面，提供了一种基于布尔运算的电阻抗成像方法，包括以下步骤：
建立完备电极模型，通过有限元方法求解得到观测模型；
建立基于B样条曲线形式的基本单元形状的矩阵表示形式；
对B样条曲线进行布尔运算；
基于B样条曲线与布尔运算相结合的形状重构算法；
执行重构算法，进行图像重构，以实现基于布尔运算的电阻抗成像。
其中，所述方法适用于绝对成像、多相成像和差分成像。
其中，所述完备电极模型如下所示：
其中，σ(x)为电导率分布，x∈Ω为空间坐标，zl为接触阻抗，Ul和Il分别为电极el上的电压和电流，n表示外单位法向。
其中，所述观测模型的一般表达式为：
V＝U(σ)+e；
其中，V为测量电压，U(σ)为使用有限元法求解得到的正问题解，即计算电压，e为附加高斯噪音。
其中，所述建立基于B样条曲线形式的基本单元形状的矩阵表示形式包括：
将每一个基本单元边界均表示为B样条曲线形式；
将B样条曲线形式经过递归处理；
得到B样条曲线形式的矩阵表示形式。
其中，所述布尔运算对基本单元进行的操作包括联合、相交和相减。
其中，所述基于B样条曲线与布尔运算相结合的形状重构算法包括：
各基本单元经过布尔运算组成新的介质边界形状，而后对介质分别定义相应的电导率，电导率分布与B样条曲线的控制点相关，可表示为
其中，在所述电阻抗成像方法适用于绝对成像的情况下，所述观测模型具体表示为：
其中，V为测量电压，U为通过正问题求解所得的电压，e为附加噪声，δ为有限元网格离散程度参数，Ω为求解域。
其中，在所述电阻抗成像方法适用于绝对成像的情况下，基于最小二乘法及正则化技术，利用高斯牛顿方法求解下列表达式所示的最小化问题实现介质边界的形状重构：
式中，Le为观测噪声协方差矩阵
的Cholesky因子，满足
Reg(·)为正则化项。
其中，所述高斯牛顿方法是通过迭代进行求解所述最小化问题的，迭代过程不断修正电导率值以及基本单元的控制点参数；和/或
在迭代过程中，需要求解雅克比矩阵
其中，
能够由标准方法求解得到，
能够通过扰动的方法进行计算；和/或
迭代过程采用线性搜索方法，直到迭代步长小于一个正值或者等于零时，迭代终止，得到最后的重建图像。
基于上述技术方案可知，本发明的基于布尔运算的电阻抗成像方法相对于现有技术至少具有如下有益效果之一或其中的一部分：
(1)利用B样条曲线表示基本单元，迭代过程中通过调整控制点的位置来对基本单元的形状进行局部的灵活调整，使算法具有重建物体细节特征的能力；
(2)将曲线的控制点作为设计变量，从而将图像重建问题转换为形状重建问题，有效降低未知数的数目以及计算时间与成本；
(3)结合布尔运算使算法具备拓扑演化能力，因此在进行重建之前无需知晓区域内所含的物体个数，从而使得算法能够适用于解决更多的实际问题；
(4)利用布尔运算对基本单元进行联合、相交、相减等操作可以构建得到多相、内嵌等复杂形状；
(5)本发明所提出的算法可以重建得到多相、内嵌等复杂形状，得到的重建电阻抗图像成像精度、空间分辨率高，并且对电阻抗成像中普遍存在的模型误差、测量噪音等具有很好的鲁棒性。
附图说明
图1是本发明实施例提供的联合、相减运算相结合的运算过程示意图；
图2是本发明实施例提供的方法流程图；
图3是本发明实施例提供的3次B样条闭合曲线示意图；
图4是本发明实施例提供的布尔运算后各区域电导率分布示意图；
图5是本发明实施例提供的算法流程图。
具体实施方式
随着科技的发展，人们生活水平不断地提高，人们开始更多地关注身体健康，越来越多的人希望得到更好的医疗服务。医学成像技术是各项医疗检查、医疗诊断中最重要的方式之一。目前常用的成像方式主要有MRI、CT、超声等，通过此类技术医生可以清楚的看到人体内部器官组织的成像，从而做出准确的诊断。但众所周知，CT、超声、MRI等也存在着不少的缺陷，不仅费用高昂，并且存在辐射，长期使用对人体有着非常大的伤害。
而电阻抗成像方法具有无创、无辐射、使用方便且设备价格相对低廉等优势，成为了近年来国内外研究的热点，在医学成像领域拥有着巨大的发展潜能，同时EIT成像技术在新冠肺炎的监护、诊断等方面起着重要作用。传统EIT算法，已经发展了30多年，但其成像精度、实时性、稳定性以及复杂形状重建等多方面仍存在不少局限性。而在医学和工业上大量存在重建区域内部复杂，对所含物体个数不明确的情况，因此提出一种具备拓扑演化能力，能够重建出多相、内嵌等复杂形状，且高分辨率和鲁棒性的电阻抗成像算法具有很大的实际需求。
本发明公开了一种基于布尔运算的电阻抗成像技术。该方法充分利用了B样条曲线的局部可修改性以及灵活控制的特性、布尔运算的拓扑演化以及复杂形状构造能力等优势，在原有的成像算法上进一步提升成像的质量，具备多相、内嵌等复杂形状的重建能力，并且对电阻抗成像中普遍存在的模型误差、测量噪音等具有很好的鲁棒性和稳定性。
具体地，本发明公开了一种基于布尔运算的电阻抗成像方法，包括以下步骤：
建立完备电极模型，通过有限元方法求解得到观测模型；
建立基于B样条曲线形式的基本单元形状的矩阵表示形式；
对B样条曲线进行布尔运算；
基于B样条曲线与布尔运算相结合的形状重构算法；
执行重构算法，进行图像重构，以实现基于布尔运算的电阻抗成像。
其中，所述方法适用于绝对成像、多相成像和差分成像。
其中，所述完备电极模型如下所示：
其中，σ(x)为电导率分布，x∈Ω为空间坐标，zl为接触阻抗，Ul和Il分别为电极el上的电压和电流，n表示外单位法向。
其中，所述观测模型的一般表达式为：
V＝U(σ)+e；
其中，V为测量电压，U(σ)为使用有限元法求解得到的正问题解，即计算电压，e为附加高斯噪音。
其中，所述建立基于B样条曲线形式的基本单元形状的矩阵表示形式包括：
将每一个基本单元边界均表示为B样条曲线形式；
将B样条曲线形式经过递归处理；
得到B样条曲线形式的矩阵表示形式。
其中，所述布尔运算对基本单元进行的操作包括联合、相交和相减。
其中，所述基于B样条曲线与布尔运算相结合的形状重构算法包括：
各基本单元经过布尔运算组成新的介质边界形状，而后对介质分别定义相应的电导率，电导率分布与B样条曲线的控制点相关，可表示为
其中，在所述电阻抗成像方法适用于绝对成像的情况下，所述观测模型具体表示为：
其中，V为测量电压，U为通过正问题求解所得的电压，e为附加噪声，δ为有限元网格离散程度参数，Ω为求解域。
其中，在所述电阻抗成像方法适用于绝对成像的情况下，基于最小二乘法及正则化技术，利用高斯牛顿方法求解下列表达式所示的最小化问题实现介质边界的形状重构：
式中，Le为观测噪声协方差矩阵
的Cholesky因子，满足
Reg(·)为正则化项。
其中，所述高斯牛顿方法是通过迭代进行求解所述最小化问题的，迭代过程不断修正电导率值以及基本单元的控制点参数；和/或
在迭代过程中，需要求解雅克比矩阵
其中，
能够由标准方法求解得到，
能够通过扰动的方法进行计算；和/或
迭代过程采用线性搜索方法，直到迭代步长小于一个正值或者等于零时，迭代终止，得到最后的重建图像。
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明作进一步的详细说明。
本发明的目的是提出基于B样条曲线与布尔运算相结合的电阻抗成像方法，对医学、工业中存在的多相、内嵌等复杂介质进行重建，在一定程度上解决由于噪音和误差等所导致的电阻抗重建图像空间分辨率以及成像精度低的问题。
本发明给定多条B样条曲线作为初始迭代曲线，通过变动相应控制点的位置来调整B样条曲线所表示的基本单元的形状，并利用单个或多个组合形式的布尔运算来对基本单元进行操作，实现迭代过程图形形状的改变，利用高斯牛顿方法迭代改变电导率分布以及由正问题模型计算得到的电压值，从而通过逆问题求解最小化目标泛函，最终实现图像重建。
通过对基本单元进行联合、相交、相减等布尔运算操作构建得到复杂介质形状。以带有内嵌物质的肺部形状为例，可以通过多个形状联合运算组成肺部形状，再对其进行相减运算得到内嵌的物体，相关运算过程示意图如图1所示。
本发明提出基于布尔运算的电阻抗成像技术，不仅具有复杂形状的重建能力，而且对电阻抗成像中普遍存在的模型误差、测量噪音等均具有很好的鲁棒性。
本发明采用B样条曲线形式表示基本单元，并与布尔运算相结合来进行EIT形状重建，充分地利用了B样条曲线的局部可调整性以及布尔运算的拓扑演化等能力。该方法所具有的优势在于：①利用B样条曲线表示基本单元，通过调整控制点的位置来对基本单元进行局部的灵活控制，从而具有重建物体细节特征的能力；②将B样条曲线的控制点作为设计变量，很大程度上减少了未知数的维度以及计算成本和时间；③结合布尔运算使得算法具有了拓扑演化能力，重建时无需知晓区域内部物体的个数；④利用布尔运算对基本单元进行联合、相交、相减等操作可以构建得到多相、内嵌等复杂形状等。基于以上优势，产生了将其与电阻抗成像技术相结合的想法。
本发明主要针对EIT技术的逆问题进行了创新，初始给定多条B样条曲线(基本单元)，以此作为高斯牛顿迭代方法的初始值，通过求解目标泛函，对基本单元进行布尔运算操作(联合、相交、相减等)得到迭代形状，从而改变电导率分布以及计算电压值，最后实现形状重建。
本发明的流程示意图如图2所示。
本发明的具体技术方案如下：
1、正问题：对所需要求解区域进行有限元剖分，建立(FEM)有限元模型。
将16个电极均匀等距地放置在感兴趣区域边界上，轮流恒流激励其中的两个电极，引起区域内部的电压，通过建立完备电极模型(CEM)求解正问题，并在其余电极上测量相应的电压值。
CEM模型如下：
其中，σ(x)为电导率分布，x∈Ω为空间坐标，zl为接触阻抗，Ul和Il分别为电极el上的电压和电流，n表示外单位法向。
利用有限元法逼近可以得到正问题的观测模型：
V＝U(σ)+e
其中V为测量电压，U(σ)为使用有限元法求解得到的正问题解，即计算电压，e为加性高斯噪音，均值为e*，协方差为Γe。
2、基于B样条曲线形式的基本单元形状隐式表达：
将每一个基本单元边界均表示为B样条曲线形式。以2维平面上的n+1个顶点作为B样条曲线的控制顶点，相应的均匀3次B样条曲线C可表示为B样条基函数的线性表示
其中，Ni，k(t)表示第i个k次(通常将k设为3)B样条基函数，通过递归的形式定义如下：
其中i＝0，1，…，n，并且
表示均匀节点向量，定义为
B样条基函数与B样条曲线控制顶点是相互独立定义的，因此，B样条曲线可表示为矩阵相乘的形式
C＝NP
此处，N为B样条基函数所组成的矩阵，P为B样条曲线控制点横纵坐标所组成的矩阵。
如图3所示，为3次B样条闭合曲线示意图，图中，黑色的点为B样条曲线的控制点，黑色实线为B样条曲线，圆圈圈出的点表示重复控制点。
3、对B样条曲线进行布尔运算：
通过B样条曲线表达基本单元边界，定义水平集函数f(x)，在B样条曲线内f(x)＞0，在B样条曲线上f(x)＝0，在B样条曲线外f(x)＜0。由此，对于Nc条B样条曲线
可得Nc个相对应的水平集函数{fj，j＝1，2，…，Nc}。初始给定的B样条曲线可能是两两重叠，或是与其余B样条曲线均不重叠，而在优化过程当中，原先重叠的两条B样条曲线可能会不再重叠，原先不重叠的两条B样条曲线可能会变得相互重叠，这也是本发明所提出的算法的一大优势：拓扑演化能力。针对两两重叠的B样条曲线进行布尔运算相当于对水平集函数进行取max或取min计算。对多条B样条曲线可以采用单种布尔运算操作，同时也可以采用多种布尔运算组合形式对相应的水平集函数进行取max和取min的组合操作。
基于B样条曲线表示与布尔运算相结合的形状重构算法：
各基本单元经过布尔运算组成新的介质边界形状，而后对介质分别定义相应的电导率，下图给出了几种简要形状的电导率分布示意图。值得注意的是示意图中给出了单条B样条曲线情况、两条B样条曲线分别进行单种形式布尔运算以及多条B样条曲线进行多种组合形式布尔运算的情况。在实际算法设计当中，还存在着各式各样的形状以及布尔算法组合，这也体现了本发明所提出的电阻抗重建算法的灵活性。如图4所示，为布尔运算后各区域电导率分布示意图。
综上可知，电导率分布与B样条曲线的控制点相关，可表示为
从数学理论上来看，电阻抗成像技术的求解属于求解椭圆型偏微分方程，在绝对成像框架下其观测模型可以表示为：
其中，V为测量电压，U为通过正问题求解所得的电压，e为附加噪声，δ为有限元网格离散程度参数，Ω为求解域。
在差分成像框架下其观测模型可以表示为：
ΔV≈JΔσ+Δe
其中，ΔV为前后两组测量电压的差值，
为在初始背景电导率值假设下电压关于电导率的雅克比矩阵，
为待求解的前后两组电导率差值，Δe为两组附加噪声之差。
为简便表示，此处以绝对成像为例，基于最小二乘法及正则化技术，利用高斯牛顿方法求解下列最小化问题实现介质边界的形状重构：
式中，Le为观测噪声协方差矩阵
的Cholesky因子，满足
Reg(·)为正则化项。
通过高斯牛顿方法迭代求解上述最小化问题，迭代过程不断修正电导率值以及基本单元的控制点参数。
在迭代过程中，需要求解雅克比矩阵
其中，
可由标准方法求解得到，
可通过扰动的方法进行计算。
迭代过程采用线性搜索方法，直到迭代步长小于一个比较小的正值或者等于零时，迭代终止，得到最后的重建图像。如图5所示为算法的具体流程图。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于布尔运算的电阻抗成像方法，其特征在于，包括以下步骤：
建立完备电极模型，通过有限元方法求解得到观测模型；
建立基于B样条曲线形式的基本单元形状的矩阵表示形式；
对B样条曲线进行布尔运算；
基于B样条曲线与布尔运算相结合的形状重构算法；
执行重构算法，进行图像重构，以实现基于布尔运算的电阻抗成像。
2.根据权利要求1所述的电阻抗成像方法，其特征在于，所述方法适用于绝对成像、多相成像和差分成像。
3.根据权利要求1所述的电阻抗成像方法，其特征在于，所述完备电极模型如下所示：
其中，σ(x)为电导率分布，x∈Ω为空间坐标，zl为接触阻抗，Ul和Il分别为电极el上的电压和电流，n表示外单位法向。
4.根据权利要求1所述的电阻抗成像方法，其特征在于，所述观测模型的一般表达式为：
V＝U(σ)+e；
其中，V为测量电压，U(σ)为使用有限元法求解得到的正问题解，即计算电压，e为附加高斯噪音。
5.根据权利要求1所述的电阻抗成像方法，其特征在于，所述建立基于B样条曲线形式的基本单元形状的矩阵表示形式包括：
将每一个基本单元边界均表示为B样条曲线形式；
将B样条曲线形式经过递归处理；
得到B样条曲线形式的矩阵表示形式。
6.根据权利要求1所述的电阻抗成像方法，其特征在于，所述布尔运算对基本单元进行的操作包括联合、相交和相减。
7.根据权利要求1所述的电阻抗成像方法，其特征在于，所述基于B样条曲线与布尔运算相结合的形状重构算法包括：
各基本单元经过布尔运算组成新的介质边界形状，而后对介质分别定义相应的电导率，电导率分布与B样条曲线的控制点相关，可表示为
8.根据权利要求2所述的电阻抗成像方法，其特征在于，在所述电阻抗成像方法适用于绝对成像的情况下，所述观测模型具体表示为：
其中，V为测量电压，U为通过正问题求解所得的电压，e为附加噪声，δ为有限元网格离散程度参数，Ω为求解域。
9.根据权利要求2所述的电阻抗成像方法，其特征在于，在所述电阻抗成像方法适用于绝对成像的情况下，基于最小二乘法及正则化技术，利用高斯牛顿方法求解下列表达式所示的最小化问题实现介质边界的形状重构：
式中，Le为观测噪声协方差矩阵
的Cholesky因子，满足
Reg(·)为正则化项。
10.根据权利要求9所述的电阻抗成像方法，其特征在于，所述高斯牛顿方法是通过迭代进行求解所述最小化问题的，迭代过程不断修正电导率值以及基本单元的控制点参数；和/或
在迭代过程中，需要求解雅克比矩阵
其中，
能够由标准方法求解得到，
能够通过扰动的方法进行计算；和/或
迭代过程采用线性搜索方法，直到迭代步长小于一个正值或者等于零时，迭代终止，得到最后的重建图像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011259573.8A CN112401866A (zh) | 2020-11-11 | 2020-11-11 | 基于布尔运算的电阻抗成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011259573.8A CN112401866A (zh) | 2020-11-11 | 2020-11-11 | 基于布尔运算的电阻抗成像方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112401866A true CN112401866A (zh) | 2021-02-26 |
Family
ID=74831431
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011259573.8A Pending CN112401866A (zh) | 2020-11-11 | 2020-11-11 | 基于布尔运算的电阻抗成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112401866A (zh) |
Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20030003148A (ko) * | 2002-11-02 | 2003-01-09 | 윤석인 | 다각형 부울리안 연산을 위한 선분 데이터 생성 방법 |
| CN1875883A (zh) * | 2006-06-08 | 2006-12-13 | 天津大学 | 基于傅立叶变换的混频生物阻抗测量方法 |
| WO2008011716A1 (en) * | 2006-07-27 | 2008-01-31 | Tasc Ltd. | High definition impedance imaging |
| CN101510225A (zh) * | 2009-03-26 | 2009-08-19 | 山东理工大学 | 产品stl模型布尔运算方法 |
| CN101609565A (zh) * | 2009-05-22 | 2009-12-23 | 南京师范大学 | 基于L-Rep模型的三维实体布尔运算方法 |
| US20100198101A1 (en) * | 2007-09-24 | 2010-08-05 | Xubo Song | Non-invasive location and tracking of tumors and other tissues for radiation therapy |
| CN102486809A (zh) * | 2010-12-03 | 2012-06-06 | 中国科学院沈阳自动化研究所 | 义齿cad模型的重建方法 |
| CN102579043A (zh) * | 2012-02-27 | 2012-07-18 | 天津大学 | 一种开放扫描电极模式的电阻抗成像系统 |
| CN103462605A (zh) * | 2013-09-06 | 2013-12-25 | 南京邮电大学 | 一种生物电阻抗成像方法 |
| US20140221864A1 (en) * | 2013-02-05 | 2014-08-07 | Drager Medical Gmbh | Electric impedance tomography device and method |
| CN105232044A (zh) * | 2015-09-11 | 2016-01-13 | 中国人民解放军第四军医大学 | 一种多频电阻抗断层成像的谱成像方法 |
| CN106037650A (zh) * | 2016-06-13 | 2016-10-26 | 河北工业大学 | 一种混合变差生物电阻抗成像方法 |
| CN106503801A (zh) * | 2016-10-18 | 2017-03-15 | 天津工业大学 | 基于深度学习的电阻抗层析成像方法 |
| CN109584323A (zh) * | 2018-10-18 | 2019-04-05 | 天津大学 | 超声反射信息约束的腹部病变电阻抗图像重建方法 |
| US20190130642A1 (en) * | 2016-05-24 | 2019-05-02 | Technion Research & Development Foundation Limited | Systems and methods for generating volumetric models |
| CN110141495A (zh) * | 2019-05-23 | 2019-08-20 | 广州中康先觉健康科技有限责任公司 | 一种测量人体穴位电阻抗的适形定位方法 |
| CN110251130A (zh) * | 2019-05-24 | 2019-09-20 | 中国科学技术大学 | 基于可移动变形组件的电阻抗成像方法、装置及系统 |
| WO2020094056A1 (zh) * | 2018-11-06 | 2020-05-14 | 苏州艾吉威机器人有限公司 | 基于b样条曲线的agv运动控制方法 |
| CN111311703A (zh) * | 2020-01-21 | 2020-06-19 | 浙江工业大学 | 一种基于深度学习的电阻抗断层图像重构方法 |
-
2020
- 2020-11-11 CN CN202011259573.8A patent/CN112401866A/zh active Pending
Patent Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20030003148A (ko) * | 2002-11-02 | 2003-01-09 | 윤석인 | 다각형 부울리안 연산을 위한 선분 데이터 생성 방법 |
| CN1875883A (zh) * | 2006-06-08 | 2006-12-13 | 天津大学 | 基于傅立叶变换的混频生物阻抗测量方法 |
| WO2008011716A1 (en) * | 2006-07-27 | 2008-01-31 | Tasc Ltd. | High definition impedance imaging |
| US20100198101A1 (en) * | 2007-09-24 | 2010-08-05 | Xubo Song | Non-invasive location and tracking of tumors and other tissues for radiation therapy |
| CN101510225A (zh) * | 2009-03-26 | 2009-08-19 | 山东理工大学 | 产品stl模型布尔运算方法 |
| CN101609565A (zh) * | 2009-05-22 | 2009-12-23 | 南京师范大学 | 基于L-Rep模型的三维实体布尔运算方法 |
| CN102486809A (zh) * | 2010-12-03 | 2012-06-06 | 中国科学院沈阳自动化研究所 | 义齿cad模型的重建方法 |
| CN102579043A (zh) * | 2012-02-27 | 2012-07-18 | 天津大学 | 一种开放扫描电极模式的电阻抗成像系统 |
| US20140221864A1 (en) * | 2013-02-05 | 2014-08-07 | Drager Medical Gmbh | Electric impedance tomography device and method |
| CN103462605A (zh) * | 2013-09-06 | 2013-12-25 | 南京邮电大学 | 一种生物电阻抗成像方法 |
| CN105232044A (zh) * | 2015-09-11 | 2016-01-13 | 中国人民解放军第四军医大学 | 一种多频电阻抗断层成像的谱成像方法 |
| US20190130642A1 (en) * | 2016-05-24 | 2019-05-02 | Technion Research & Development Foundation Limited | Systems and methods for generating volumetric models |
| CN106037650A (zh) * | 2016-06-13 | 2016-10-26 | 河北工业大学 | 一种混合变差生物电阻抗成像方法 |
| CN106503801A (zh) * | 2016-10-18 | 2017-03-15 | 天津工业大学 | 基于深度学习的电阻抗层析成像方法 |
| CN109584323A (zh) * | 2018-10-18 | 2019-04-05 | 天津大学 | 超声反射信息约束的腹部病变电阻抗图像重建方法 |
| WO2020094056A1 (zh) * | 2018-11-06 | 2020-05-14 | 苏州艾吉威机器人有限公司 | 基于b样条曲线的agv运动控制方法 |
| CN110141495A (zh) * | 2019-05-23 | 2019-08-20 | 广州中康先觉健康科技有限责任公司 | 一种测量人体穴位电阻抗的适形定位方法 |
| CN110251130A (zh) * | 2019-05-24 | 2019-09-20 | 中国科学技术大学 | 基于可移动变形组件的电阻抗成像方法、装置及系统 |
| CN111311703A (zh) * | 2020-01-21 | 2020-06-19 | 浙江工业大学 | 一种基于深度学习的电阻抗断层图像重构方法 |
Non-Patent Citations (3)
| Title |
|---|
| 何炳蔚, 林志航: "多个B样条曲面光滑重建技术研究", 计算机应用研究, no. 06, 28 June 2003 (2003-06-28) * |
| 寿华好;胡良臣;: "法向约束下B样条曲线拟合的实编码GA算法", 浙江工业大学学报, no. 04, 25 August 2016 (2016-08-25) * |
| 郁道银, 黄家祥, 陈晓冬, 孙正, 谢洪波: "基于B样条的冠状动脉树骨架三维重建方法", 工程图学学报, no. 02, 30 June 2005 (2005-06-30) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Saulnier et al. | Electrical impedance tomography | |
| Harikumar et al. | Electrical impedance tomography (EIT) and its medical applications: a review | |
| Ren et al. | A statistical shape-constrained reconstruction framework for electrical impedance tomography | |
| Zhang et al. | Supervised descent learning for thoracic electrical impedance tomography | |
| Hinkle et al. | 4D CT image reconstruction with diffeomorphic motion model | |
| WO2010062220A1 (ru) | Способ неинвазивного электро физиологического исследования сердца | |
| Murphy et al. | Effect of domain shape modeling and measurement errors on the 2-D D-bar method for EIT | |
| Wang et al. | Exploring respiratory motion tracking through electrical impedance tomography | |
| CN111616708A (zh) | 一种精准识别脑卒中颅内病变区域的图像重建方法 | |
| CN106373194B (zh) | 一种人体肺部电阻层析成像有限元模型设计方法 | |
| Li et al. | Electrical-impedance-tomography imaging based on a new three-dimensional thorax model for assessing the extent of lung injury | |
| CN110811596B (zh) | 基于低秩与稀疏约束和非局部全变分的无创心脏电位重建方法 | |
| Song et al. | A nonlinear weighted anisotropic total variation regularization for electrical impedance tomography | |
| Wu et al. | Bayesian shape reconstruction using B-spline level set in electrical impedance tomography | |
| Dimas et al. | An efficient point-matching method-of-moments for 2D and 3D electrical impedance tomography using radial basis functions | |
| Zhu et al. | Electrical Impedance Tomography Guided by Digital Twins and Deep Learning for Lung Monitoring | |
| Pengpan et al. | A motion-compensated cone-beam CT using electrical impedance tomography imaging | |
| Murphy et al. | Reconstructions of conductive and insulating targets using the D-bar method on an elliptical domain | |
| CN109946388B (zh) | 基于统计逆的电学/超声双模态内含物边界重建方法 | |
| CN112401866A (zh) | 基于布尔运算的电阻抗成像方法 | |
| Cao et al. | 2D image reconstruction of a human chest by using Calderon's method and the adjacent current pattern | |
| Rosa et al. | Urinary bladder volume monitoring using magnetic induction tomography: A rotational simulation model for anatomical slices within the pelvic region | |
| Zhang et al. | Survey of EIT image reconstruction algorithms | |
| CN114708350A (zh) | 一种用于脑部电阻抗层析成像的电导率可视化方法 | |
| CN117274413B (zh) | 一种基于eit的电导率图像重建方法、系统及设备 |
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 |