一种电阻抗成像的稀疏重建方法专利检索-电阻抗成像测量仪器和方法专利检索查询-专利查询网

一种 电阻抗成像的稀疏重建方法

阅读：154发布：2020-07-17

专利汇可以提供一种电阻抗成像的稀疏重建方法专利检索，专利查询，专利分析的服务。并且本发明公开了一种电阻抗成像的稀疏重建方法，根据动态电阻抗图像的电导率分布具有稀疏性的特点，在进行电阻抗成像的重构时，将重构的目标函数的正则项采用L1范数，并将非均匀的电导率的稀疏性引入到正则项中，采用凸优化问题求解方法计算目标函数，从而得到重构的电阻抗成像。本发明通过采用L1范数正则项有效抵抗噪声，将电导率的非均匀性引入正则项中，从而提高电阻抗图像质量和抗噪能力，而且运算速度快。，下面是一种电阻抗成像的稀疏重建方法专利的具体信息内容。

权利要求

1. 一种电阻抗成像的稀疏重建方法，其特征在于，根据动态电阻抗图像的电导率分布具有稀疏性的特点，在进行电阻抗成像的重构时，将重构的目标函数的正则项采用LI范数，并将非均匀的电导率的稀疏性引入到正则项中，采用凸优化问题求解方法计算目标函数，从而得到重构的电阻抗图像。
2.如权利要求1所述的电阻抗成像的稀疏重建方法，其特征在于，所述的动态电阻抗成像问题看做一个线性欠定系统Ax=b，其中A是敏感性矩阵，b是边界测量电压值；求解的目标函数为min||/l.x:-&| +Ilxll，拟合项采用L2范数最小二乘min||Jx-，正则项采用LI范数I !X2-X1I I1, A为正则化参数，目标为动态电阻抗图像的非均匀电导率X = X2-X1,其中X1是前一帧的电导率分布，X2是当前帧的电导率分布，X是差分成像结果；目标函数采用凸优化问题求解方法计算。
3.如权利要求2所述的电阻抗成像的稀疏重建方法，其特征在于，所述的\为10e_5。
4.如权利要求2所述的电阻抗成像的稀疏重建方法，其特征在于，所述目标函数的求解方法为预校验共轭梯度法，计算方法如下：给定：e >0,初始化
循环：（1)计算搜索方向（Ax，Au)作为牛顿系统// 的近似解，其中
是海森矩阵，g = ▽戎是当前迭代（x，u)的梯度； (2)采用回溯法线性搜索计算s，计算方法如下：定义当前迭代步为X, A和V,下一迭代步为X+, X +和V+,即
;关于y+的留数定义为r+ ； ①计算最大正步长：
初始化回溯：s=0. 99s—，用3去乘s，直到f(x+) < 0 ；持续用 3 去乘 s,直到 I |rt(x+，A +, v+) I I a (l_a s) I |rt(x, A , v) | |2, a =0. 01,^ =0. 5 ； (3 )迭代更新
(4)构造对偶可行点 v ：
(5)评估对偶间隙 (6)如果，退出；

说明书全文

—种 电阻抗成像的稀疏重建方法

技术领域

[0001] 本发明属于电阻抗成像技术领域，涉及一种电阻抗成像的稀疏重建方法。

背景技术

[0002] 电阻抗断层成像技术是一种新型医学成像技术，该技术基于不同的生物组织具有不同的电阻率特性、病理生理功能的改变也会显著改变组织电阻抗特性这一特点，通过在人体目标区域贴放一圈体表电极，并通过这些电极向体内注入弱的、对人体完全无创的交流电流，再测量各电极上的响应电压信号，最后经特定的图像重构算法实现某一截面或区域内组织电阻抗特性的二维或三维分布状况的成像。由于该技术对人体无创无害，且相对现有成像设备而言，该技术具有结构简单、操作简便、设备便携、成像速度快、系统造价低等优点，能够满足长时间、动态观察病变部位的演变过程的特殊需求，因而有望弥补现有医学成像技术的不足，是当前相关领域的研究热点。

[0003] 但在实际应用中，由于测量模式的限制，即人体被测者体表只能设置有限数量的电极，以及测量过程中不可避免的误差和噪声引使得图像存在严重伪影。人体内部的某些重要电阻抗变化可能无法通过电阻抗图像得到及时反映，不利于电阻抗成像的临床推广和应用。因此，为了提高电阻抗图像的成像质量，需要一种能够对噪声不敏感的成像方法。

发明内容

[0004] 本发明解决的问题在于提供一种电阻抗成像的稀疏重建方法，该方法对噪声不敏感，速度更快，成像质量更好。

[0005] 本发明是通过以下技术方案来实现：

[0006] 一种电阻抗成像的稀疏重建方法，根据动态电阻抗图像的电导率分布具有稀疏性的特点，在进行电阻抗成像的重构时，将重构的目标函数的正则项采用LI范数，并将非均匀的电导率的稀疏性引入到正则项中，采用凸优化问题求解方法计算目标函数，从而得到重构的电阻抗图像。

[0007] 所述的动态电阻抗成像问题看做一个线性欠定系统Ax=b，其中A是敏感性矩阵，b是边界测量电压值；求解的目标函数为min||AY-/?|g +K，拟合项采用L2范数最小二mmm\\Ax-bf2 ,正则项采用Li范数I |X2-Xl| I1, A为正则化参数，目标为动态电阻抗图像
的非均匀电导率X = X2-X1,其中X1是前一帧的电导率分布，X2是当前帧的电导率分布，X是差分成像结果；目标函数采用凸优化问题求解方法计算。

[0008] 所述的入为10e_5 ;在目标函数求解时，A、b为已知条件。

[0009] 所述目标函数的求解方法为预校验共轭梯度法，计算方法如下：

[0010]给定：e>0，初始化：=丄 x=0，u=l=(l，…，I) G Rn;/I[0011] 循环：（1)计算搜索方向（Ax, Au)作为牛顿系统
的近似解，其中
是海森矩阵，
於是当前迭代（X，U)的梯度；

[0012] (2)采用回溯法线性搜索计算S，计算方法如下：

[0013] 定义当前迭代步为X, X和V,下一迭代步为X+,入+和V+,即X+= x+s A X, A+ =入+s A A , V+ = v+s A V ;关于y+的留数定义为r+ ；

[0014] ①计算最大正步长：

[0015]

[0016] 初始化回溯：s=0. 99s—，用@去乘s，直到f(x+) < 0 ；

[0017]持续用 @ 去乘 s,直到

[0021] (6)如果退出；

[0022] (7)更新

[0023] 与现有技术相比，本发明具有以下有益的技术效果：

[0024] 本发明提供的电阻抗成像的稀疏重建方法，通过分析动态电阻抗图像的电导率分布X = X2-X1,具体的电导率在数值上大部分值为0或者近似为0，得出其非均匀的电导率具有非均匀性的特点，然后进行电阻抗成像的重构时引入稀疏重建技术，将非均匀的电导率的稀疏性引入到目标函数的正则项中。

[0025] 本发明提供的电阻抗成像的稀疏重建方法，通过采用LI范数正则项有效抵抗噪声，将电导率的非均匀性引入正则项中，从而提高电阻抗图像质量和抗噪能力，而且运算速度快。附图说明

[0026] 图1为动态电阻抗图像非均匀电导率差分成像图，突显出来的区域即为成像目标，可将背景电导率的值看为0或者近似为0

[0027] 图2为动态电阻抗图像非均匀电导率具有非均匀性；横坐标为元素的数目，纵坐标为电导率值

[0028] 图3中成像结果对比图，其中左列是给定的成像目标，中列是传统Tikhonov成像方法，右列是稀疏重建方法成像结果图，噪声水平80dB。具体实施方式

[0029] 下面结合具体的实施例对本发明做进一步的详细说明，所述是对本发明的解释而不是限定。

[0030] 在电阻抗动态成像问题中，通过分析动态电阻抗图像的电导率分布X = X2-X1,具体的电导率在数值上大部分值为0或者近似为0。如图1，图2目标为动态电阻抗图像的非均匀电导率X = X2-X1，其中X1是前一帧的电导率分布，X2是当前帧的电导率分布，X是差分成像结果，因此X是两帧之间的变化，因此电导率X具有稀疏性的特点，然后进行电阻抗成像的重构时引入稀疏重建技术。

[0031] 动态电阻抗成像问题可以看做一个线性欠定系统，Ax=b，其中A是敏感性矩阵，b是边界测量电压值。电导率分布是X = X2-X1,其中X1是前一帧的电导率分布，X2是当前帧的电导率分布，X是差分成像结果。

[0032]目标函数的构造采用最小二乘法，拟合项采用L2范数最小二乘min||Ax-b||22，正则项采用LI范数||X2-X1||,目标为动态电阻抗图像的非均匀电导率X = X2-X1

[0033]目标函数为min||Ar-b||22+λ||x||,其求解是一个凸优化问题，可以采用常用的凸优化问题求解方法去计算，其中动态电阻抗图像的非均匀电导率为X = X2-X1, λ为正则化参数取10e-5。

[0034]目标函数的求解方法采用现有的预校验共轭梯度法，具体计算方法如下：

[0035]给定：e >0，初始化：t=1/λ，x=0，u=l=(l,…，I) ∈Rn

[0036] 循环：（1)计算搜索方向（Δx, Δu)作为牛顿系统h[ΔxΔu]=-g的近似解，其中H=▽2Φt(x,u)∈R2n×2n是海森矩阵，g =▽Φt(x,u）∈R2n 是当前迭代（x，u)的梯度；

[0037] (2)采用回溯法线性搜索计算S，具体计算方法如下：

[0038] 定义当前迭代步为X, A和V,下一迭代步为X+, λ +和V+,即X+ = x+s Δ x,λ+= λ +s Δ λ , v+=v+s ΔV。关于y+的留数定义为r+。

[0039] ①计算最大正步长：

[0040]

[0041] 初始化回溯：s=0. 99smax，用β去乘S，直到f (x+) < 0。

[0042]持续用 β去乘 s,直到 || rt (X+，λ +, v+)|| a ≤(1- a s) | | rt (x,λ，v) ||2,这里a =0. 01, β =0. 5。

[0043] (3 )迭代更新：(x, u) = (x, u)+s( Δx, Δu)；

[0044] (4)构造对偶可行点 V :

[0045] (5)评估对偶间隙 η

[0046] (6)如果η/G(v)≤ε,退出；[0047] (7)更新

[0048] 在仿真条件下（EID0RS3. 5 MATLAB2012a)，边界测量电压施加80dB的高斯白噪声条件下，将上述方法与传统方法进行对比，具体的如图3所示，其中左列为给定的成像目标图，中列为传统方法结果图，右列为稀疏重建结果图，且图中数据为各方法的耗时记录结果(耗时记录单位为秒，记为S)。

[0049] 对比成像质量可以看出：

[0050] 传统方法：伪影多，中心区域不能成像，速度慢；

[0051] 本发明的方法：伪影少，中心区域成像，速度快。

标题	发布/更新时间	阅读量
变电站接地网数字化管理系统	2020-05-20	590
一种非接触式电阻抗层析成像装置及方法	2020-05-21	278
一种基于电阻抗成像技术的检测系统和方法	2020-05-13	614
一种红外生命监护穿衣镜及其使用方法	2020-05-23	513
基于空间谱估计算法的电阻抗成像方法	2020-05-18	494
基于电阻抗成像损伤监测的复合材料强度预测方法	2020-05-22	622
基于低秩稀疏分解的电阻抗图像重建方法	2020-05-21	128
一种具有生命监护功能的穿衣镜及其使用方法	2020-05-16	305
一种用于生物电阻抗成像的数据采集装置	2020-05-15	805
一种电阻抗成像的数据采集摄像机用连接机构	2020-05-08	593

一种电阻抗成像的稀疏重建方法

—种电阻抗成像的稀疏重建方法

技术领域

背景技术

发明内容

该功能需要专业版企业版VIP权限，您可以：