CN103018284B - 一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法 - Google Patents
一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法 Download PDFInfo
- Publication number
- CN103018284B CN103018284B CN201210544370.2A CN201210544370A CN103018284B CN 103018284 B CN103018284 B CN 103018284B CN 201210544370 A CN201210544370 A CN 201210544370A CN 103018284 B CN103018284 B CN 103018284B
- Authority
- CN
- China
- Prior art keywords
- terminal
- impedance
- relay
- impedance measurement
- terminal 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.)
- Active
Links
Landscapes
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法,属于电学无损检测技术领域。所述方法包括如下步骤:建立电阻抗层析成像传感器的有限元模型,计算相邻激励、相对激励或者对角线激励模式下的四端子灵敏度矩阵;继电器双T型多通道切换模块选通传感器的电极,进行二端子阻抗测量;将二端子阻抗值转换成相应激励模式下的四端子阻抗值;求解敏感场内电导率差异分布,并进行图像重建。本发明适用于非侵入式电阻抗层析成像技术,尤其适用于四端子阻抗测量方式不满足应用需求或者不具备四端子阻抗测量仪器的情形。本发明可有效提高所重建图像的对比度和分辨率,加快成像速度。
Description
【技术领域】
本发明属于电学无损检测技术领域,尤其涉及一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法。
【背景技术】
近三十余年发展起来的电阻抗层析成像技术(Electrical Impedance Tomography—EIT),以其非侵入性、便携性、价格低廉、响应快速等技术优势,在工业和医学领域具有重要的应用前景。EIT技术实质上是根据敏感场的电导率分布获得物场的媒质分布信息。通过在敏感场边界施加激励电流,当场内电导率分布变化时,导致场内电势分布变化,从而场域边界上的测量电压发生变化,利用获得的电阻抗信息通过一定的图像重建算法,可以重建出场内的电导率分布,从而获得物场的媒质分布。
国家知识产权局在电阻抗层析成像领域只授权了一个发明专利:基于微针电极的电阻抗层析成像仪及其微创式测量方法(ZL200610114600.6),利用微针电极跨越人体皮肤高阻抗的角质层,替代传统表面电极,作为电流激励装置以及测量电压信号的传感器,测量信号经过放大、滤波后输入PC机进行数据处理,最后以灰度图或彩图反映断层上各点的电阻抗分布信息。可以达到降低电流激励与与采集带来的误差,同时使得实测对象区域内的电阻抗分布更趋于一阶连续,从而使得重构所得图像分辨率更高,置信度更高。
在电阻抗测量中,二端子模式和四端子模式是两种常见的测量模式,但是各有其优缺点。二端子阻抗测量的最主要缺点是会受到接触阻抗的影响。J.J.Ackmann于1993年指出指出当进行四端子阻抗测量,频率高于5kHz时,会引入不能接受的相位差(“Complex bioelectric impedance measurement system for thefrequency range from5Hz to 1MHz,”Ann.Biomed.Eng.,vol.21,no.2,pp.135-146,Mar.1993.)。Szczepanik和Rucki于2000年指出当进行四端子阻抗测量,频率高于5kHz时,在响应电极上测得的电压会有较大的波动(Z.Szczepanik,and Z.Rucki,“Frequency analysis of electrical impedance tomography system,”IEEE Trans.Instrum.Meas.,vol.49,no.4,pp.844-851,Aug.2000.)。Dickin于年指出四端子阻抗测量时容易产生电流源饱和现象,使得测量上限被限制(F.Dickin,and M.Wang,“Electrical resistance tomography for processapplications,”Meas.Sci.Technol.,vol.7,no.3,pp.247-260,Mar.1996.)。
在电阻抗层析成像中,二端子模式和四端子模式也是两种常见的成像模式。在二端子阻抗测量模式下,对传感器的每两个电极依次施加激励电流,再测量这两个电极上感应到的电压值,计算得到二端子电阻抗值。四端子模式下,对传感器两个电极施加激励电流,再测量另两个电极上感应到的电压值,计算得到四端子电阻抗值。对于四端子模式,可以选择相邻激励、相对激励或者对角线等多种激励方式。在利用灵敏度理论来进行电阻抗层析成像中,二端子模式下的成像质量和速度不及四端子模式下的成像质量和速度。
加之在实际应用中,经常会遇到四端子模式不能满足应用需求或者没有四端子阻抗测量仪器两种情形。因此极为有必要发明出一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法,来克服上述缺陷,可有效提高成像对比度和分辨力,并能加快成像速度。
【发明内容】
本发明的目的是提供一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法,以能够在四端子模式不能满足应用需求或者不具备四端子阻抗测量仪器的情形,仍能提供高质量、高速度的图像重建。
为实现上述目的,本发明提供的一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法,采用如下技术方案:
一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法,其特征在于,包括如下步骤:
步骤一,建立电阻抗层析成像(EIT)传感器11的有限元分析模型;
步骤二,根据Geselowitz和Lehr提出的灵敏度理论,结合EIT传感器11的有限元分析模型,计算所述EIT传感器11在相邻、相对或对角线等任意激励模式下的四端子灵敏度矩阵S;
步骤三,利用继电器双T型多通道切换模块14选通电极12为测量端子一、测量端子二、接地、浮空四个状态之一,在传感器只含均匀导电介质和被测对象13进入传感器中这两种情形下,分别通过二端子阻抗测量模块15进行二端子阻抗测量,测量频率为10kHz,并通过GPIB-USB接口将数据实时上传至PC机16中;
步骤四,一对电极A-B被电流激励,测得另一对电极C-D的感应电压,计算所得的四端子阻抗值记为ZCD←AB;四对电极A-C,A-D,B-C或者B-D分别被施加电流激励,并测得该对电极之间的感应电压,计算所得的二端子阻抗值记为zAC,zAD,zBC或者zBD;根据本发明提供的的阻抗测量数据转换公式将所测得的二端子阻抗值矩阵可转换成相邻、相对或对角线等不同激励模式下的四端子阻抗值矩阵;在传感器只含均匀导电介质的情形下的四端子阻抗值矩阵记为Z0,被测对象13进入传感器中的情形下的四端子阻抗值矩阵记为Z,根据ΔZ=Z-Z0计算出ΔZ;
步骤五,利用逆问题求解方法求解ΔZ=SΔσ,得到敏感场内电导率差异分布Δσ;
步骤六,结合EIT传感器11的有限元分析模型,利用所得的电导率差异分布进行图像重建,在PC机16中显示被测对象13的形状、位置等信息。
其中,一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法,其特征在于所述方法步骤三中的切换二端子阻抗测量通道的继电器双T型多通道切换模块14采用Omron公司出品的微型双刀双掷继电器开关G6K构成双T型开关结构来消除二端子阻抗测量的接触阻抗,具体结构为:继电器22的转换触点222连接电极12,继电器22的常开触点221连接继电器21的转换触点212,继电器22的常闭触点223连接继电器23的转换触点232;继电器21的常开触点211连接二端子阻抗测量模块15的测量端子一,继电器21的常闭触点213连接地;继电器23的常开触点231连接二端子阻抗测量模块15的测量端子二,继电器23的常闭触点233浮空;所述方法还通过电极镀铂;选择测量频率大于10kHz来消除二端子阻抗测量的接触阻抗。
本发明提供的一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法的优点是:
1.本发明适用于非侵入式电阻抗层析成像技术,尤其适用于四端子阻抗测量方式不满足应用需求或者不具备四端子阻抗测量仪器的情形。
2.本发明的阻抗测量数据转换公式可通过有限次的二端子阻抗测量有效地变换成相邻、相对或对角线等任意激励模式下的四端子阻抗测量,大大提高了数据的利用率。
3.本发明可有效提高所重建图像的对比度和分辨率,加快成像速度。
4.本发明的的切换二端子阻抗测量通道的继电器双T型多通道切换模块采用Omron公司出品的微型双刀双掷继电器开关G6K构成双T型开关结构来消除二端子阻抗测量的接触阻抗,可通过尽量少的切换动作实现选择电极为端子一、端子二、接地、浮空四个状态之一。本发明还通过给电极镀铂,测量频率选在在10kHz以上等措施有效消除二端子阻抗测量的主要缺点接触阻抗的干扰。
【说明书附图】
图1为一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法示意图;
图2为继电器双T型开关模块的示意图;
图3为四端子阻抗测量模式示意图;
图4为二端子阻抗测量模式示意图;
图5为四种被测对象放入传感器中实验示意图;
图6为一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法经共轭梯度法迭代25次重建图像;
图7为二端子电阻抗层析成像方法经共轭梯度法迭代25次重建图像;
图8为二端子电阻抗层析成像方法经共轭梯度法迭代100次重建图像。
【具体实施方式】
下面结合附图对本发明的一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法的实施方式作进一步的说明,给出实验所得重建图像,并与二端子电阻抗层析成像进行重建结果对比。
步骤一,参照图3,建立电阻抗层析成像(EIT)传感器11的有限元分析模型;
步骤二,参照图3,根据Geselowitz和Lehr提出的灵敏度理论,结合EIT传感器11的有限元分析模型,计算所述EIT传感器11在相邻、相对或对角线等任意激励模式下的四端子灵敏度矩阵S;
步骤三,参照图1-2,继电器双T型开关模块工作方式如下:默认情况下,继电器不动作,电极12可置于浮空状态;当信号控制继电器22的常开触点221闭合,继电器21的常开触点211闭合,电极12可连接精密阻抗分析仪15的端子一;当信号控制继电器23的常开触点231闭合,电极12可连接精密阻抗分析仪15的端子二;当信号控制继电器22的常开触点221闭合,电极12可连接地。PC机16通过USB或者RS-232接口控制继电器双T型开关模块14通过尽量少的切换动作选择电极12为端子一、端子二、接地、浮空四个状态之一;参照图5,被测对象13包括1支圆玻璃棒51,2支圆玻璃棒52,3支圆玻璃棒53,4支玻璃棒54;参照图1、图5,在传感器只含均匀导电介质和被测对象13进入传感器中这两种情形下,分别通过二端子阻抗测量模块15进行二端子阻抗测量,测量频率为10kHz,并通过GPIB-USB接口将数据实时上传至PC机16中;
步骤四,参照图3,一对电极A-B被电流激励,测得另一对电极C-D的感应电压,计算所得的四端子阻抗值记为ZCD←AB;参照图4,四对电极A-C,A-D,B-C或者B-D分别被施加电流激励,并测得该对电极之间的感应电压,计算所得的二端子阻抗值记为zAC,zAD,zBC或者zBD;参照图3-4,根据本发明提供的的阻抗测量数据转换公式将所测得的二端子阻抗值矩阵可转换成相邻、相对或对角线等不同激励模式下的四端子阻抗值矩阵;在传感器只含均匀导电介质的情形下的四端子阻抗值矩阵记为Z0,被测对象13进入传感器中的情形下的四端子阻抗值矩阵记为Z,根据ΔZ=Z-Z0计算出ΔZ;
步骤五,利用逆问题求解方法求解ΔZ=SΔσ,得到敏感场内电导率差异分布Δσ;
步骤六,参照图1,结合EIT传感器11的有限元分析模型,利用所得的电导率差异分布进行图像重建,在PC机16中显示被测对象13的形状、位置等信息;参照图6,本发明的一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法经共轭梯度法迭代25次获得1支玻璃棒重建图像61,2支玻璃棒重建图像62,3支玻璃棒重建图像63,4支玻璃棒重建图像64。
为了更客观地说明本发明的一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法的可行性与有效性,二端子电阻抗层析成像方法分别经共轭梯度法迭代25次和迭代100次进行图像重建。参照图7,二端子电阻抗层析成像方法经共轭梯度法迭代25次获得1支玻璃棒重建图像71,2支玻璃棒重建图像72,3支玻璃棒重建图像73,4支玻璃棒重建图像74。参照图8,二端子电阻抗层析成像方法经共轭梯度法迭代100次获得1支玻璃棒重建图像81,2支玻璃棒重建图像82,3支玻璃棒重建图像83,4支玻璃棒重建图像84。
比较图6、图7说明在相同迭代次数下,本发明的方法所重建的图像在对比度、分辨力都远远优于二端子电阻抗层析成像方法。比较图6、图7、图8说明二端子电阻层析成像方法迭代次数增加到4倍,尽管重建图像质量有所提高,但是仍不能完全区分出四支玻璃棒,因此其成像质量仍不及本发明的方法。表1给出二端子电阻抗层析成像方法与本发明的方法经共轭梯度法重建一幅50*50像素的图像所需时间。结合表1、图6-8,说明本发明的方法在采用共轭梯度法来求解逆问题进行图像重建0.8438秒可获得比二端子电阻抗层析成像方法重建2.9958秒更高的成像质量。
表1.二端子电阻抗层析成像方法与本发明的方法重建图像所需时间
因此,本发明适用于非侵入式电阻抗层析成像技术,尤其适用于四端子阻抗测量方式不满足应用需求或者不具备四端子阻抗测量仪器的情形。本发明可有效提高所重建图像的对比度和分辨率,加快成像速度。
以上所述仅为本发明具体实施方法的基本方案,但本发明的保护范围并不局限于此,任何熟悉本技术领域的人员在本发明公开的技术范围内,可想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。所有落入权利要求的等同的含义和范围内的变化都将包括在权利要求的范围之内。
Claims (2)
1.一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法,其特征在于,包括以下步骤:
步骤一,建立电阻抗层析成像传感器(11)的有限元分析模型;
步骤二,在电阻抗层析成像传感器(11)的有限元分析模型中,根据Geselowitz和Lehr提出的灵敏度理论计算所述电阻抗层析成像传感器(11)在相邻或相对激励模式下的四端子灵敏度矩阵S;
步骤三,利用继电器双T型多通道切换模块(14)选通电阻抗层析成像传感器(11)中的任一电极(12)为测量端子一、测量端子二、接地、浮空四个状态之一,在电阻抗层析成像传感器(11)只含均匀导电介质的情形下和被测对象(13)进入该传感器(11)的均匀导电介质中的情形下,分别通过二端子阻抗测量模块(15)进行二端子阻抗测量,测量频率为10kHz,并通过GPIB-USB接口将数据实时上传至PC机(16)中;
步骤四,电极A、B、C和D是电阻抗层析成像传感器(11)在四端子阻抗测量模式中的四个电极,一对电极A-B被电流激励,测得另一对电极C-D的感应电压,计算所得的四端子阻抗值记为ZCD←AB;四对电极A-C,A-D,B-C和B-D分别被施加电流激励,并测得该对电极之间的感应电压,计算所得的二端子阻抗值分别记为zAC,zAD,zBC和zBD;当激励电极A和B处于相邻或相对的位置时,由阻抗测量数据转换公式将所测得的二端子阻抗值矩阵转换成相邻或相对激励模式下的四端子阻抗值矩阵;在电阻抗层析成像传感器(11)只含均匀导电介质的情形下的四端子阻抗值矩阵记为Z0,被测对象(13)进入该传感器(11)的均匀导电介质中的情形下的四端子阻抗值矩阵记为Z,根据ΔZ=Z-Z0计算出ΔZ;
步骤五,利用逆问题求解方法求解ΔZ=SΔσ,得到敏感场内电导率差异分布Δσ;
步骤六,结合电阻抗层析成像传感器(11)的有限元分析模型,利用所得的电导率差异分布进行图像重建,在PC机(16)中显示被测对象(13)的形状、位置信息。
2.根据权利要求1所述的基于二端子阻抗测量模式的四端子电阻抗层析成像方法,其特征在于,所述方法步骤三中所使用的切换二端子阻抗测量通道的继电器双T型多通道切换模块采用Omron公司出品的微型双刀双掷继电器开关G6K构成双T型开关结构来消除二端子阻抗测量的接触阻抗,具体结构为:第二继电器(22)的转换触点(222)连接电阻抗层析成像传感器(11)的任一电极(12),第二继电器(22)的常开触点(221)连接第一继电器(21)的转换触点(212),第二继电器(22)的常闭触点(223)连接第三继电器(23)的转换触点(232);第一继电器(21)的常开触点(211)连接二端子阻抗测量模块(15)的测量端子一,第一继电器(21)的常闭触点(213)连接地;第三继电器(23)的常开触点(231)连接二端子阻抗测量模块(15)的测量端子二,第三继电器(23)的常闭触点(233)浮空;所述方法还通过电极镀铂,选择测量频率大于10kHz来消除二端子阻抗测量的接触阻抗。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210544370.2A CN103018284B (zh) | 2012-12-14 | 2012-12-14 | 一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210544370.2A CN103018284B (zh) | 2012-12-14 | 2012-12-14 | 一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103018284A CN103018284A (zh) | 2013-04-03 |
CN103018284B true CN103018284B (zh) | 2015-05-13 |
Family
ID=47967128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210544370.2A Active CN103018284B (zh) | 2012-12-14 | 2012-12-14 | 一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103018284B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105982671B (zh) * | 2015-03-03 | 2020-07-07 | 思澜科技(成都)有限公司 | 一种优化的生物阻抗测试电路 |
CN107677706A (zh) * | 2017-08-24 | 2018-02-09 | 西安理工大学 | 高速离心机多相介质分布的可视无线检测系统及检测方法 |
JP2020165923A (ja) * | 2019-03-29 | 2020-10-08 | 日本電産リード株式会社 | 検査装置、及び検査方法 |
CN111766276A (zh) * | 2020-06-16 | 2020-10-13 | 中国农业大学 | 植物茎秆冻融测量方法、装置、系统、设备及存储介质 |
CN113533439B (zh) * | 2021-06-25 | 2022-09-02 | 中国农业大学 | 玉米穗水分分布的电阻抗成像方法、系统及电子设备 |
CN114288555B (zh) * | 2022-01-26 | 2023-02-28 | 云南贝泰妮生物科技集团股份有限公司 | 一种基于皮肤电阻抗的射频美容仪自适应调频系统 |
CN117617939B (zh) * | 2022-12-16 | 2024-09-13 | 上海融易迈医疗健康科技有限公司 | 电阻抗成像加速方法、存储介质及电阻抗成像设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101241094A (zh) * | 2008-03-12 | 2008-08-13 | 天津大学 | 非接触式电阻抗传感器及基于该传感器的图像重建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002053029A1 (en) * | 2000-12-30 | 2002-07-11 | University Of Leeds | Electrical impedance tomography |
-
2012
- 2012-12-14 CN CN201210544370.2A patent/CN103018284B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101241094A (zh) * | 2008-03-12 | 2008-08-13 | 天津大学 | 非接触式电阻抗传感器及基于该传感器的图像重建方法 |
Non-Patent Citations (1)
Title |
---|
"用组合算法逐步提高电阻层析成像的分辨率";石天明等;《仪器仪表学报》;20081231;第29卷(第12期);第2587-2593页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103018284A (zh) | 2013-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103018284B (zh) | 一种基于二端子阻抗测量模式的四端子电阻抗层析成像方法 | |
Wang et al. | Split Bregman iterative algorithm for sparse reconstruction of electrical impedance tomography | |
CN102688041B (zh) | 基于电极采用交错位方式排列的三维电阻抗断层成像方法 | |
CN102961136B (zh) | 一种用于电阻抗断层成像系统的定标装置 | |
CN102106731A (zh) | 一种电阻抗断层成像测量方法 | |
CN102138790A (zh) | 一种三维电阻抗断层成像系统性能测试方法与装置 | |
Chen et al. | Four-terminal imaging using a two-terminal electrical impedance tomography system | |
CN108711178A (zh) | 一种基于闭环控制原理的电容层析成像图像重建方法 | |
Beretta et al. | A reconstruction algorithm based on topological gradient for an inverse problem related to a semilinear elliptic boundary value problem | |
Yue et al. | 3D electromagnetic tomography using a single layer sensor array | |
CN109758149B (zh) | 一种结合时差成像和频差成像的电阻抗成像方法 | |
Kumar et al. | Recent prospects of medical imaging and sensing technologies based on electrical impedance data acquisition system | |
CN103065286A (zh) | 一种准静态电阻抗成像中的图像重构方法 | |
Grysiński et al. | Planar sensors for local conductivity measurements in biological objects—Design, modelling, sensitivity maps | |
Pan et al. | Research on defect scanning and spatial locating based on CCIT sensor of single-pair electrode | |
Ye et al. | Evaluation of effect of number of electrodes in ERT sensors on image quality | |
Zhu et al. | A Wideband Contactless Electrical Impedance Tomography System | |
Zhao et al. | The impact of the measurement accuracy and the excitation pattern on EIT image reconstruction | |
Sun et al. | A new method for electrical impedance tomography with incomplete electrode array | |
CN109102552A (zh) | 一种非均匀形状约束的像素值域滤波超声成像重建方法 | |
ChuanLei et al. | New selection methods of regularization parameter for electrical resistance tomography image reconstruction | |
Yang et al. | An image reconstruction algorithm for electrical impedance tomography using measurement estimation of virtual electrodes | |
CN103549954B (zh) | 一种磁共振成像三维相位去折叠方法及其装置 | |
Yan et al. | A Detection Method for Fast Electrical Impedance Imaging of Grounding Grid Based on Optimized Differential-Multigrid-Homotopy Algorithm | |
CN104614411B (zh) | 基于p向量等差下降的电学层析成像Lp正则化重建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |