提升电阻率成像获取图片质量的方法与流程

本发明涉及一种提升电阻率成像获取图片质量的方法。

背景技术：

目前，电阻率成像方法主要包括反投影算法和牛顿迭代法。其中，牛顿迭代法理论比较完善、应用比较广泛、生成图片的清晰度较高，但可能出现不收敛的情况，若强行降低收敛标准，会导致算法精度下滑。而反投影算法的优点主要是简单直观、计算量小、速度快，对信息量的要求不高，利用该算法图像重建能够比较容易实现。但这种成像算法也有明显的不足，主要是成像精度低、比较模糊，很多情况下，无法得到场域内准确的电阻率分布，只能得到物质分布的大体轮廓，更多的只能用于定性的在线观测分析。

近年来，国内外许多研究机构对反投影算法和牛顿迭代法进行了改进完善，提出了诸如滤波反投影算法、单步灵敏度系数算法、基于奇异值分解的直接法等，但这些方法要么精度低、成像质量差，要么处理过程复杂、计算量大、速度慢。随着现代医学的发展，医生在临床对图像清晰度逐步提升，有必要提升传统的电阻率成像算法的成像精度。

技术实现要素：

为解决现有技术的不足，本发明提供了一种提升电阻率成像获取图片质量的方法，利用图像风格迁移技术对电阻率成像算法获得的图像进行二次处理，在不增加数据获取需求的情况下有效地提高生成图像质量。

一种提升电阻率成像获取图片质量的方法，包括以下步骤：

对仿真设定生成数据进行电阻率成像算法计算；

将仿真设定生成数据和经过电阻率成像算法计算生成的算法生成数据作为输入特征进行风格迁移训练，其中，利用正问题有限元法计算得到的待测区域电压值与实测得到的待测区域电压值进行比较判断是否终止迭代过程；

经过多次重复训练获得与实际模型最为相符的雅可比矩阵；

雅可比矩阵求逆即为算法生成数据与仿真设定生成数据的映射函数；

将实际情况下采用电阻率成像算法生成的图像带入风格迁移训练获得的映射函数中进行处理获得新图像。

进一步地，待测区域电压值与实测得到的待测区域电压值满足公式则表示迭代终止；

其中，v⁽ⁿ⁾为待测区域电压值，v0为实测得到的待测区域电压值，ε为设定的系统为保证准确性所允许的最大误差值，||v⁽ⁿ⁾-v0||2则是(v⁽ⁿ⁾-v0)的2范数。

一种提升电阻率成像获取图片质量的方法，包括以下步骤：

对仿真设定生成数据进行电阻率成像算法计算；

将仿真设定生成数据和经过电阻率成像算法计算生成的算法生成数据作为输入特征进行风格迁移训练获得算法生成数据与仿真设定生成数据的映射函数；

将实际情况下采用电阻率成像算法生成的图像带入风格迁移训练获得的映射函数中进行处理获得新图像。

进一步地，在风格迁移训练时，利用正问题有限元法计算得到的待测区域电压值与实测得到的待测区域电压值进行比较判断是否终止迭代过程。

进一步地，测区域电压值与实测得到的待测区域电压值满足公式则表示迭代终止；

其中，v⁽ⁿ⁾为待测区域电压值，v0为实测得到的待测区域电压值，ε为设定的系统为保证准确性所允许的最大误差值，||v⁽ⁿ⁾-v0||2则是(v⁽ⁿ⁾-v0)的2范数。

进一步地，经过多次重复训练获得与实际模型最为相符的雅可比矩阵；雅可比矩阵求逆即为算法生成数据与仿真设定生成数据的映射函数。

本发明的有益之处在于可以在不增加数据获取需求的情况下有效地提高生成图像质量。

附图说明

图1是一种提升电阻率成像获取图片质量的方法的流程图；

图2是电阻率成像图片的示意图；

图3是对图2中电阻率成像图片采用提升电阻率成像获取图片质量的方法进行处理后的示意图。

具体实施方式

以下结合附图和具体实施例对本发明作具体的介绍。

如图1至图3所示，一种提升电阻率成像获取图片质量的方法，包括以下步骤：

对仿真设定生成数据进行电阻率成像算法计算；

将仿真设定生成数据和经过电阻率成像算法计算生成的算法生成数据作为输入特征进行风格迁移训练获得算法生成数据与仿真设定生成数据的映射函数；

将实际情况下采用电阻率成像算法生成的图像带入风格迁移训练获得的映射函数中进行处理获得新图像。

具体而言，一种提升电阻率成像获取图片质量的方法，包括以下步骤：

对仿真设定生成数据进行电阻率成像算法计算；

将仿真设定生成数据和经过电阻率成像算法计算生成的算法生成数据作为输入特征进行风格迁移训练，其中，利用正问题有限元法计算得到的待测区域电压值与实测得到的待测区域电压值进行比较判断是否终止迭代过程；

经过多次重复训练获得与实际模型最为相符的雅可比矩阵；

雅可比矩阵求逆即为算法生成数据与仿真设定生成数据的映射函数；

将实际情况下采用电阻率成像算法生成的图像带入风格迁移训练获得的映射函数中进行处理获得新图像。

作为一种具体的实施方式，待测区域电压值与实测得到的待测区域电压值满足公式则表示迭代终止；其中，v⁽ⁿ⁾为待测区域电压值，v0为实测得到的待测区域电压值，ε为设定的系统为保证准确性所允许的最大误差值，||v⁽ⁿ⁾-v0||2则是(v⁽ⁿ⁾-v0)的2范数。

假设初始情况下待测区域电阻率均匀分布，均为ρ⁽⁰⁾。可以通过电阻率成像正问题有限元法获得边界电势计算值，可以假设第n步迭代的待测区域电阻率分布中间估值为ρ⁽ⁿ⁾，可以获得v⁽ⁿ⁾，可以用公式表示为v⁽ⁿ⁾＝f(ρ⁽ⁿ⁾)。若迭代终止进入算法的ρ⁽ⁿ⁾修正环节。

对ρ⁽ⁿ⁾的修正需要借助雅可比矩阵的求解，首先写出逐步迭代方程：

ρ⁽ⁿ⁺¹⁾＝ρ⁽ⁿ⁾+δρ⁽ⁿ⁺¹⁾。

其中，δρ⁽ⁿ⁺¹⁾为迭代法的修正量，需要借助辅助函数f(ρ)进行计算：

以上作为辅助函数帮助判别ρ⁽ⁿ⁺¹⁾取何值时会获得迭代方法的最小误差，通过其一阶导数的零点进行搜索。

即：f′(ρ⁽ⁿ⁺¹⁾)＝0公式展开可写为：[f′(ρ)]＝[f′(ρ)]^t(f(ρ)-v0)＝0。

其中f′(ρ)即为雅可比矩阵，具体为：

通过多次重复仿真可以获得与实际模型最为相符的雅可比矩阵。其中，k为待测区域内的三角有限元个数，m为测量到的电势总组数。对雅可比矩阵的求解需要借助补偿定理法，常用于相邻电流的激励模式。

以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解，上述实施例不以任何形式限制本发明，凡采用等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。