CN110097608A - 修正路径追踪描述的连续波超声层析成像重建方法 - Google Patents

修正路径追踪描述的连续波超声层析成像重建方法 Download PDF

Info

Publication number
CN110097608A
CN110097608A CN201910195283.2A CN201910195283A CN110097608A CN 110097608 A CN110097608 A CN 110097608A CN 201910195283 A CN201910195283 A CN 201910195283A CN 110097608 A CN110097608 A CN 110097608A
Authority
CN
China
Prior art keywords
indicate
field domain
path
pixel
indicates
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
CN201910195283.2A
Other languages
English (en)
Other versions
CN110097608B (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201910195283.2A priority Critical patent/CN110097608B/zh
Publication of CN110097608A publication Critical patent/CN110097608A/zh
Application granted granted Critical
Publication of CN110097608B publication Critical patent/CN110097608B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)

Abstract

本发明涉及一种修正路径追踪描述的连续波超声层析成像重建方法，包括如下步骤：获取边界测量值；基于收发探头相对几何位置进行直线线性投影，构建系数矩阵；使用同步代数重建方法进行成像迭代计算求得场域内衰减系数的分布；使用高斯滤波函数对步骤三的结果进行图像滤波；使用直方图均衡化映射对步骤四的结果进行图像增强；使用差分曲率方法对步骤五的结果进行内含物轮廓提取，该方法通过计算能量熵函数判断图像各像素点是否对应内含物轮廓，在内含物轮廓对应像素点，能量熵函数较大；根据步骤六中提取的内含物轮廓求取声速分布；对给定发射超声探头和接收超声探头进行路径追踪计算；构建新的系数矩阵投影衰减测量值τ和系数矩阵R；迭代终止。

Description

修正路径追踪描述的连续波超声层析成像重建方法
技术领域
本发明属于超声层析成像技术领域，涉及实现一种采用修正路径追踪方法描述的连续波超声成像重建 方法，用于实现场域内低对比度两相介质的重建。
背景技术
超声层析成像技术(Ultrasonic Tomography,UT)是一种结构性成像技术，其通过在被测场域外布置超 声传感器阵列并施加一定的激励以得到边界电压测量数据，以此来重建被测场域内部的折射系数、衰减系 数或声阻抗分布情况。相比其他软场成像技术例如电阻抗层析成像(Electrical Impedance Tomography,EIT) 和电磁层析成像(MagneticImpedance Tomography,MIT)，UT具有非侵入、分辨率高的优点，相比精度较 高的硬场成像技术如X射线层析成像(X-ray Computed Tomography,X-CT)及光学层析成像方法(OpticalCoherence Tomography,OCT)，UT使用安全、结构简单、可以实现实时成像。此外UT还有着非接触、方 向性好、成本低等优势，是一种较为理想的过程可视化检测监视手段。UT作为一种层析成像技术手段， 在多相流可视化检测、化工石油输送、航空发动机探查以及生物医学诊断中均有广泛的应用。
完整的UT系统主要包含三个部分：超声换能器设计、制造与安装；信号激励与采集系统；超声成像 图像重建算法。其中超声成像算法通过对从采集系统得到的换能器接收信号进行处理，通过解调提取测量 幅值或渡越时间，得到某个确定激励下的全部换能器的有效测量数据，进一步通过图像重建方法得到场域 内含物介质分布的合理估计。目前，超声成像算法作为一种主要利用硬场特性的成像方法，超声成像方法 严重依赖于场域边界换能器的数量，其逆问题求解具有严重的病态性(对测量值得微小扰动会导致重建结 果的大幅度变化)和欠定性(所需求解的方程数远小于未知量的数目，方程有无穷多解)。为克服这个问 题，学者们提出了许多图像重建算法，其中，基于路径的投影重建算法是一种克服病态性的有效手段。该 方法通过计算激励、接收换能器间的路径，将收发探头间的时间延迟或幅值衰减均匀的分配给所计算路径 上的每一个像素，通过对不同收发探头间的路径进行计算并对同一像素在不同路径上的估计值进行叠加， 得到场域内每个像素值的有效估计，以达到可视化测量和图像重建的目的。典型的投影重建方法有徐立军 等人1998年在《仪器仪表学报(Chinese Journal of Scientific Instrument)》第17卷，1-7页，发表的“气液 两相泡状流体监测用超声层析成像系统的研究(Investigation of UltrasoundTomography System used for Monitoring Bubbly Gas/Liquid Two-phase fluid)”中提到的二值反投影方法、Rahim等人在《Sensors and Actuators(传感器与执行器)》第135卷，337-345页发表的“Non-invasive imaging of liquid/gas flowusing ultrasonictransmission mode tomography(超声对液体、气体的非侵入性成像)”中提到的采用R-L函数的 线性滤波反投影方法、Gordon等人在《Journal of theoretical biology(理论生物学期刊)》第29卷，第3期， 471-481页发表的“Algebraic reconstruction techniques(ART)for three-dimensional electron microscopy and X-rayphotography(用于三维电子显微镜和X射线CT的代数重建技术)”中提出的代数重建方法、苏邦良 等人在《Chemical Engineering Journal(化学工程期刊)》第77卷，37-41页发表的“The use ofsimultaneous iterative reconstruction technique for electrical capacitancetomography(同步迭代重建技术在电容层析成像中 的应用)”中提出的同步迭代重建方法、Anderson等人在《Ultrasound Imaging(超声成像)》第6卷，81-94 页发表的“Simultaneousalgebraic reconstruction technique(SART):a superior implementation of the ARTalgorithm(同步代数重建技术(SART)：ART算法更优越的实现)”中提出的同步代数重建方法等。其中， SART算法以其收敛快速、残差较小的优点被广泛使用在超声成像的研究中。
不同的应用对象对超声成像图像重建方法提出了不同的要求，面向工业两相流动过程尤其是以油水两 相流为代表的弱声阻抗比介质分布，超声成像重建算法需要满足实时性的要求(即成像速度每秒大于24 幅)。此外，针对弱声阻抗比介质，除基于已有线性成像模型设计高精度成像算法外，还应对其正问题中 的非线性现象进行有效刻画。在满足实时性要求方面：鲍勇等人在《Measurement Science and Technology (测量科学和技术)》第28卷，074002页发表的“Real-time temperature fieldmeasurement based onacoustic tomography(基于超声层析成像的实时温度场测量)”中提出了多探头同步激励方法；在正问题的非线性现 象描述方面，Anders等人在《Ultrasound Imaging(超声成像)》第12卷，268-291页发表的“ARay-tracing approach to restoration and resolutionenhancement in experimental ultrasound tomography(实验超声断层扫描 中恢复和分辨率增强的射线追踪方法)”中推导了基于程函方程的超声传播路径追踪模型，瞿晓磊等人在 《Japanese Journal of Applied Physics(日本应用物理学报)》第56卷，07JF14页发表的“Computational cost reduction by avoiding ray-linking iteration in bent-raymethod for sound speed image reconstruction in ultrasound computed tomography(使用超声成像进行声速图像重建时通过避免在弯曲射线方法中进行射线链接迭代来 降低计算成本)”中提出了虚拟接收换能器的概念以提升路径追踪的计算时间，贾佳斌等人在《2017IEEE International Conference on Imaging System and Techniques(2017IEEE国际成像系统与技术大会)》发表的 “Nonlinear temperature field reconstructionusing acoustic tomography(超声层析成像非线性温度场重建)” 中应用非线性路径追踪方法进行温度场分布的有效重建。
在超声成像图像重建算法中，基于路径追踪的正问题非线性描述方法与超声成像实时性的要求间存在 矛盾。首先，基于路径追踪的正问题非线性描述方法的目的在于优化确定发射探头与接收探头的最优传播 路径，其本质为迭代优化过程，若所有发射与接收探头间传播路径均采用路径追踪方法进行计算，则重建 时间会大幅提高，无法满足实时性成像的要求。其次，场域内用于成像的像素数量的多少会直接影响重建 算法的计算量及计算时间，随着场域内像素数量提高，重建算法对计算资源的占用较大，在迭代重建算法 中，过多的像素无法满足超声成像的实时性要求。
除非线性正问题描述及超声成像实时性的矛盾外，探头数量的多少对重建图像精度及重建图像分辨率 有至关重要的影响，及超声成像重建算法与探头间的有效投影路径数目紧密相关：投影路径越多，成像精 度越高，伪影越少。但在UPT技术实际应用的过程中，受限于场域尺寸及信号激励幅值限制，场域边界的 探头数目并不能无限增大；另一方面，超声作为机械波在场域内的传播需要一定的渡越额时间，换能器数 目过度无法满足可视化监测的实时性要求。在UPT技术的实际应用中，超声换能器的数目一般不超过32 个。较高精度的超声成像要求和较快的数据成像速度需要形成了较大的矛盾。因此需要一种在低投影数量 和低像素数量下的超声成像正问题非线性描述方法及相应图像重建算法，以期在实时性成像要求下，重建 结果达到较高的精度和较少的伪影。
发明内容
本发明针对超声层析成像逆问题图像重建，提出一种修正路径追踪描述的连续波超声成像重建方法， 可以获得高精度、高信噪比的边界测量数据，提高重建结果的求解精度和图像分辨率。技术方案如下：
一种修正路径追踪描述的连续波超声层析成像重建方法，包括如下步骤：
步骤一：获取边界测量值，在被测场域外表面均匀的布置一定数量的超声换能器，并使用正弦波电压 对发射探头进行激励。对除发射外的其他超声探头，记录其接收到正弦信号连续50个单峰值的平均值， 记为边界电压测量值。据此，获取重建所需的投影衰减测量值τ，具体计算方式为
其中，fc是激励信号的中心频率，As为背景介质(水)下的边界电压测量值，Ar为存在内含物(水 中的离散油泡)情况下的边界电压测量值，ln表示对数符号。
步骤二：基于收发探头相对几何位置进行直线线性投影，构建系数矩阵R，其元素表示为Ri,j，表 示场域内第i条投影路径穿过场域内第j个像素的相对长度，同时对应系数矩阵中第i行、第j列的元 素。
步骤三：使用同步代数重建方法(SART)进行成像迭代计算求得场域内衰减系数的分布，具体计算 方式为：
a(k+1)＝a(k)+αDp(DrR)T(τ-Ra(k))
其中，α表示迭代步长，k表示迭代次数，a(k)表示第k次迭代时的像素单元衰减系数分布， Dp＝diag(1/R+,1,1/R+,2,····,1/R+,N)，Dr＝diag(1/R1,+,1/R2,+,····,1/RM,+)，R+,1表示对系数矩阵内第一列所有 元素进行求和，R1,+表示对系数矩阵内第一行所有元素进行求和，N表示场域内像素个数，M表示场 域内投影路径数量，diag()表示对角矩阵，()T表示矩阵转置。
步骤四：使用高斯滤波函数对步骤三的结果进行图像滤波。
步骤五：使用直方图均衡化映射对步骤四的结果进行图像增强。
步骤六：使用差分曲率方法对步骤五的结果进行内含物轮廓提取，该方法通过计算能量熵函数判断图 像各像素点是否对应内含物轮廓，在内含物轮廓对应像素点，能量熵函数较大；在非内含物轮廓对应像素 点，能量熵函数约等于零。能量熵函数计算方式表示为：
其中，ax表示成像结果水平方向的一阶导数，ay表示成像结果垂直方向的一阶导数，axx表示对ax求 水平方向的偏导数，ayy表示对ay求垂直方向的偏导数，axy表示对ax求垂直方向的偏导数，|·|表示绝对 值。
步骤七：根据步骤六中提取的内含物轮廓求取声速分布，轮廓内部声速为内含物(油)的声速， 轮廓外部声速为背景介质(水)的声速。
步骤八：根据步骤七中得到的声速分布，对给定发射超声探头和接收超声探头进行路径追踪计算。 路径追踪计算过程的输入变量为起始点坐标和起始角，起始点坐标为发射超声探头坐标，起始角为发射与 接收探头间连线的角度。路径追踪的输出为该条路径上的坐标点向量和路径追踪与场域边界的交点坐标 (终止点)。路径追踪的具体计算方式为：
其中，表示下一个像素路径点的方向矢量，表示当前像素路径点的路径矢量，表 示路径方向矢量对当前像素求一阶偏导数，表示路径方向矢量对当前像素求二阶偏导数，O((Δs)3) 表示以像素尺寸为自变量的高阶无穷小，s表示路径追踪计算当前所在的像素点，Δs表示像素点的尺寸， n表示场域内折射系数分布，表示场域内折射系数分布的梯度分布，≈为约等于号。
步骤九：根据步骤八中的路径追踪终止点坐标和不同路径的坐标点向量构建新的系数矩阵投影衰 减测量值τ和系数矩阵R。投影衰减测量值τ的计算方法为根据路径追踪终止点坐标的位置对边界测量值 进行线性插值，系数矩阵R的计算方法为计算不同路径的坐标点向量构成的近似曲线同场域内不同像素相 交的长度。
步骤十：重复步骤三至步骤九直至残差满足要求其中，Rea(k)＝||R·a(k)-τ||，ε为人为设定的残差阈值，m为步骤三至步骤九迭代过程中的迭代次数。
本发明提供的修正路径追踪描述的连续波超声成像重建方法。所提出修正路径追踪描述的连续波超声 成像重建方法，实质性特点是：采用路径追踪方法精确表征超声层析成像正问题中的非线性特性，采用低 精度迭代、高精度可视化的策略进行逆问题迭代重建。在正问题描述方面,使用连续波超声激励以获得弱 声阻抗比两相介质准确的衰减信息，通过衰减分布的成像结果映射声速分布以进行路径追踪计算；提出使 用虚拟接收换能器进行路径追踪计算以大幅减小路径追踪计算时间；在逆问题重建方面，采用低分辨率成像迭代计算以满足算法实时性要求，通过使用基于周期性协方差矩阵的高斯过程回归方法在保证成像精度 的同时实现高分辨率图像重建。所提出的超声成像重建方法可以获得高精度、高信噪比的边界测量数据， 同时给出了超声传播过程中非线性特性的准确描述，在满足实时性要求的情况下，成像结果可以有效解决 传统反投影类算法在进行油水两相分布逆问题求解时成像精度差、成像分辨率低、成像伪影严重的问题。 所提算法有效扩展了超声透射成像算法的应用，在满足工业流动过程实时性成像要求的基础上，显著提高 重建结果的求解精度和图像分辨率。
附图说明
图1为修正路径追踪描述的连续波超声层析成像重建方法基本流程图；
图2为本发明中采用连续波成像的超声层析成像测量方法示意图；
图3为本发明中提出的利用虚拟接收探头进行路径追踪计算的示意图；
图4为本发明的四个典型仿真模型，并分别给出了相应的总变差正则化(TV)；
图5给出传统算法与本算法的重建指标对比。
具体实施方式
结合附图和实施例对本发明的修正路径追踪描述的连续波超声层析成像重建方法加以说明。
本发明的修正路径追踪描述的连续波超声层析成像重建方法，实施例中针对工业输油管道中油水两相 流的成像这一UPT技术的常见应用形式使用本发明所提算法进行计算。在正问题描述方面，使用连续波超 声激励以获得弱声阻抗比两相介质准确的衰减信息，通过衰减分布的成像结果映射声速分布以进行路径追 踪计算。在逆问题重建计算方面，提出使用虚拟接收换能器进行路径追踪计算以大幅减小路径追踪计算时 间，与此同时，采用低分辨率成像迭代计算以满足算法实时性要求，通过使用基于周期性协方差矩阵的高 斯过程回归方法在保证成像精度的同时实现高分辨率图像重建。所提出的超声成像重建方法可以获得高精 度、高信噪比的边界测量数据，同时给出了超声传播过程中非线性特性的准确描述，在满足实时性要求的 情况下，成像结果可以有效解决传统反投影类算法在进行油水两相分布逆问题求解时成像精度差、成像分 辨率低、成像伪影严重的问题。所提算法有效扩展了超声透射成像算法的应用，在满足工业流动过程实时 性成像要求的基础上，显著提高重建结果的求解精度和图像分辨率。
图1为本发明所提算法的基本流程图，主要分为SART迭代成像、声衰减分布映射声速分布、路径追 踪三部分
图2为超声层析成像系统的基本原理图示意，在对进行测量时，共计16个超声换能器均匀的沿管壁 安装负责激励、接收超声波。采用循环激励、一发全收的测量模式，探头按顺时针方向均匀分布。16个超 声探头按顺序接入峰峰值50V、频率1MHz的方波电压激励，探头通道切换时间间隔2.5ms。四与此同时， 16个通道同步接收稳态时刻的电压正弦信号持续1ms，并通过正交解调得到接收电压有效值。每次测量总 计获得16×15＝240个边界电压测量数据。
图3为使用虚拟接收探头的路径追踪方法示意，路径追踪起始点为发射超声探头位置坐标，路径追踪 起始角为发射探头与接收探头连线角度，将路径追踪终止点作为虚拟接收探头，根据其位置坐标进行边界 测量值线性插值，得到每条路径追踪路线对应的测试值。
图4、图5中分别给出了传统UT成像算法的成像结果与本算法的成像结果与重建指标对比，重建指 标包括相对误差(RE)和图像相关系数(CC)两种，其计算方法表示为：
其中，σ表示重建的像素单元电导率分布，σ*表示真实情况下的电导率分布，σj和σj *表示第 j个像素单元重建的和真实的电导率分布，和表示重建的和真实的电导率分布的平均值。
本算法实施例包括如下具体步骤：
(1)获取边界测量值，在被测场域外表面均匀的布置一定数量的超声换能器，并使用正弦波电压对 发射探头进行激励。对除发射外的其他超声探头，记录其接收到正弦信号连续50个单峰值的平均值，记 为边界电压测量值。据此，获取重建所需的投影衰减测量值τ，具体计算方式为
其中，fc是激励信号的中心频率，As为单一背景介质下(水)的边界电压测量值，Ar为存在内含物 (为离散油泡)情况下的边界电压测量值，ln表示对数符号。
(2)基于收发探头相对几何位置进行直线线性投影，构建系数矩阵R，其元素表示为Ri,j，表示 场域内第i条投影路径穿过场域内第j个像素的相对长度，同时对应系数矩阵中第i行、第j列的元素。
(3)：使用同步代数重建方法(SART)进行成像迭代计算求得场域内衰减系数的分布，具体计算方 式为：
a(k+1)＝a(k)+αDp(DrR)T(τ-Ra(k))
其中，α表示迭代步长，采用单因素变量法进行优化，最优取值为1.75，k表示迭代次数，a(k)表 示第k次迭代时的像素单元衰减系数分布，Dp＝diag(1/R+,1,1/R+,2,····,1/R+,N)， Dr＝diag(1/R1,+,1/R2,+,····,1/RM,+)，R+,1表示对系数矩阵内第一列所有元素进行求和，R1,+表示对系数矩 阵内第一行所有元素进行求和，N表示场域内像素个数，M表示场域内投影路径数量，diag()表 示对角矩阵，()T表示矩阵转置。
(4)：使用高斯滤波函数对(3)的结果进行图像滤波。
(5)：使用直方图均衡化映射对(4)的结果进行图像增强。
(6)：使用差分曲率方法对(5)的结果进行内含物轮廓提取，该方法通过计算能量熵函数判断图像 各像素点是否对应内含物轮廓，在内含物轮廓对应像素点，能量熵函数较大；在非内含物轮廓对应像素点， 能量熵函数约等于零。能量熵函数计算方式表示为：
其中，ax表示成像结果水平方向的一阶导数，ay表示成像结果垂直方向的一阶导数，axx表示对ax求 水平方向的偏导数，ayy表示对ay求垂直方向的偏导数，axy表示对ax求垂直方向的偏导数，|·|表示绝对 值。
(7)：根据(6)中提取的内含物轮廓求取声速分布，轮廓内部声速为内含物(油)的声速，轮 廓外部声速为背景介质(水)的声速。
(8)：根据(7)中得到的声速分布，对给定发射超声探头和接收超声探头进行路径追踪计算。 路径追踪计算过程的输入变量为起始点坐标和起始角，起始点坐标为发射超声探头坐标，起始角为发射与 接收探头间连线的角度。路径追踪的输出为该条路径上的坐标点向量和路径追踪与场域边界的交点坐标 (终止点)。路径追踪的具体计算方式为：
其中，表示下一个像素路径点的方向矢量，表示当前像素路径点的路径矢量，表 示路径方向矢量对当前像素求一阶偏导数，表示路径方向矢量对当前像素求二阶偏导数，O((Δs)3) 表示以像素尺寸为自变量的高阶无穷小，s表示路径追踪计算当前所在的像素点，Δs表示像素点的尺寸， n表示场域内折射系数分布，表示场域内折射系数分布的梯度分布，≈为约等于号。
(9)：根据(8)中的路径追踪终止点坐标和不同路径的坐标点向量构建新的系数矩阵投影衰减 测量值τ和系数矩阵R。投影衰减测量值τ的计算方法为根据路径追踪终止点坐标的位置对边界测量值进 行线性插值，系数矩阵R的计算方法为计算不同路径的坐标点向量构成的近似曲线同场域内不同像素相交 的长度。
(10)：重复(3)～(9)直至残差满足要求其中，Rea(k)＝||R·a(k)-τ||， ε为人为设定的残差阈值，这里选定为0.0001，m为(3)～(9)迭代过程中的迭代次数。
以上所述实施例为本发明的几个示例模型，本发明不局限于该实施例和附图所公开的内容。凡是不脱 离本发明所公开的精神下完成的等效或修改，都在本发明保护的范围。

Claims (2)

1.一种修正路径追踪描述的连续波超声层析成像重建方法，用于超声层析成像，包含以下步骤：
步骤一：获取边界测量值，在被测场域外表面布置一定数量的超声换能器，并使用正弦波电压对发射探头进行激励，对除发射外的其他超声探头，记录其接收到正弦信号连续多个单峰值的平均值，记为边界电压测量值，据此，获取重建所需的投影衰减测量值τ，具体计算方式为
其中，fc是激励信号的中心频率，As为场域内充满单一背景介质时的边界电压测量值，Ar为场域内存在连续的背景介质和离散的内含物介质情况下的边界电压测量值，ln表示对数符号。
步骤二：基于收发探头相对几何位置进行直线线性投影，构建系数矩阵R，其元素表示为Ri,j，表示场域内第i条投影路径穿过场域内第j个像素的相对长度，同时对应系数矩阵中第i行、第j列的元素。
步骤三：使用同步代数重建方法(SART)进行成像迭代计算求得场域内衰减系数的分布，具体计算方式为：
a(k+1)＝a(k)+αDp(DrR)T(τ-Ra(k))
其中，α表示迭代步长，k表示迭代次数，a(k)表示第k次迭代时的像素单元衰减系数分布，Dp＝diag(1/R+,1,1/R+,2,····,1/R+,N)，Dr＝diag(1/R1,+,1/R2,+,····,1/RM,+)，R+,1表示对系数矩阵内第一列所有元素进行求和，R1,+表示对系数矩阵内第一行所有元素进行求和，N表示场域内像素个数，M表示场域内投影路径数量，diag()表示对角矩阵，()T表示矩阵转置。
步骤四：使用高斯滤波函数对步骤三的结果进行图像滤波。
步骤五：使用直方图均衡化映射对步骤四的结果进行图像增强。
步骤六：使用差分曲率方法对步骤五的结果进行内含物轮廓提取，该方法通过计算能量熵函数判断图像各像素点是否对应内含物轮廓，在内含物轮廓对应像素点，能量熵函数较大；在非内含物轮廓对应像素点，能量熵函数约等于零。能量熵函数计算方式表示为：
其中，ax表示成像结果水平方向的一阶导数，ay表示成像结果垂直方向的一阶导数，axx表示对ax求水平方向的偏导数，ayy表示对ay求垂直方向的偏导数，axy表示对ax求垂直方向的偏导数，|·|表示绝对值。
步骤七：根据步骤六中提取的内含物轮廓求取声速分布，轮廓内部声速为内含物的声速，轮廓外部声速为背景介质的声速。
步骤八：根据步骤七中得到的声速分布，对给定发射超声探头和接收超声探头进行路径追踪计算。路径追踪计算过程的输入变量为起始点坐标和起始角，起始点坐标为发射超声探头坐标，起始角为发射与接收探头间连线的角度。路径追踪的输出为该条路径上的坐标点向量和路径追踪与场域边界的交点坐标(终止点)。路径追踪的具体计算方式为：
其中，表示下一个像素路径点的方向矢量，表示当前像素路径点的路径矢量，表示路径方向矢量对当前像素求一阶偏导数，表示路径方向矢量对当前像素求二阶偏导数，O((Δs)3)表示以像素尺寸为自变量的高阶无穷小，s表示路径追踪计算当前所在的像素点，Δs表示像素点的尺寸，n表示场域内折射系数分布，▽n表示场域内折射系数分布的梯度分布，≈为约等于号。
步骤九：根据步骤八中的路径追踪终止点坐标和不同路径的坐标点向量构建新的系数矩阵投影衰减测量值τ和系数矩阵R。投影衰减测量值τ的计算方法为根据路径追踪终止点坐标的位置对边界测量值进行线性插值，系数矩阵R的计算方法为计算不同路径的坐标点向量构成的近似曲线同场域内不同像素相交的长度。
步骤十：重复步骤三-步骤九直至残差满足要求其中，Rea(k)＝||R·a(k)-τ||，ε为人为设定的残差阈值，m为步骤三-步骤九迭代过程中的迭代次数。
2.根据权利要求1所述的超声层析成像重建方法，其特征在于：所述步骤一中边界测量值的获取，是指将被测对象置于超声层析成像被测场域中，场域外均匀分布多个超声换能器，对发射探头进行连续波正弦信号电压激励，采用电压激励电压测量的基本测量模式，循环激励、一发全收的测量策略，投影衰减测量值为单一背景介质存在下的边界测量电压值和背景介质中含有离散相介质的边界测量电压值的比值。
CN201910195283.2A 2019-03-14 2019-03-14 修正路径追踪描述的连续波超声层析成像重建方法 Active CN110097608B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910195283.2A CN110097608B (zh) 2019-03-14 2019-03-14 修正路径追踪描述的连续波超声层析成像重建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910195283.2A CN110097608B (zh) 2019-03-14 2019-03-14 修正路径追踪描述的连续波超声层析成像重建方法

Publications (2)

Publication Number Publication Date
CN110097608A true CN110097608A (zh) 2019-08-06
CN110097608B CN110097608B (zh) 2023-04-07

Family

ID=67443162

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910195283.2A Active CN110097608B (zh) 2019-03-14 2019-03-14 修正路径追踪描述的连续波超声层析成像重建方法

Country Status (1)

Country Link
CN (1) CN110097608B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110706298A (zh) * 2019-09-10 2020-01-17 天津大学 正则化加权最小二乘的透射反射双模式超声成像重建方法
CN112067691A (zh) * 2020-08-31 2020-12-11 天津大学 油气水三相介质衰减谱融合多频超声层析成像方法
CN113092049A (zh) * 2021-03-25 2021-07-09 北京理工大学 一种三维跨界面成像方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DK516883D0 (da) * 1982-11-12 1983-11-11 Schlumberger Ltd Fremgangsmade og anlaeg til diffraktionstomografi
CN102305828A (zh) * 2011-05-13 2012-01-04 中北大学 基于环绕式阵列的超声层析成像检测系统及其方法
CN102629326A (zh) * 2012-03-19 2012-08-08 天津工业大学 一种基于单目视觉的车道线检测方法
CN104101687A (zh) * 2014-07-10 2014-10-15 天津大学 基于超声多普勒与电学多传感器的多相流可视化测试方法
CN106226392A (zh) * 2016-07-05 2016-12-14 天津大学 基于超声衰减机理模型的油水两相流相含率测量方法
CN108287199A (zh) * 2017-12-30 2018-07-17 天津大学 基于传播路径网格剖分的超声透射模式层析成像方法
CN108442420A (zh) * 2018-03-21 2018-08-24 大连理工大学 基于分区超声波层析成像的地下连续墙缺陷检测方法
CN109102552A (zh) * 2018-07-26 2018-12-28 天津大学 一种非均匀形状约束的像素值域滤波超声成像重建方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DK516883D0 (da) * 1982-11-12 1983-11-11 Schlumberger Ltd Fremgangsmade og anlaeg til diffraktionstomografi
CN102305828A (zh) * 2011-05-13 2012-01-04 中北大学 基于环绕式阵列的超声层析成像检测系统及其方法
CN102629326A (zh) * 2012-03-19 2012-08-08 天津工业大学 一种基于单目视觉的车道线检测方法
CN104101687A (zh) * 2014-07-10 2014-10-15 天津大学 基于超声多普勒与电学多传感器的多相流可视化测试方法
CN106226392A (zh) * 2016-07-05 2016-12-14 天津大学 基于超声衰减机理模型的油水两相流相含率测量方法
CN108287199A (zh) * 2017-12-30 2018-07-17 天津大学 基于传播路径网格剖分的超声透射模式层析成像方法
CN108442420A (zh) * 2018-03-21 2018-08-24 大连理工大学 基于分区超声波层析成像的地下连续墙缺陷检测方法
CN109102552A (zh) * 2018-07-26 2018-12-28 天津大学 一种非均匀形状约束的像素值域滤波超声成像重建方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
JIANMIN LI ET AL.: "Learning localGaussianprocessregressionforimagesuper-resolution", 《NEUROCOMPUTING》 *
SASCHA LANGENER ET AL.: "A real-time ultrasound process tomography system using a reflection-mode reconstruction technique" *
XIAOLEI QU ET AL.: "Computational cost reduction by avoiding ray-linking iteration in bent-ray method for sound speed image reconstruction in ultrasound computed tomography", 《JAPANESE JOURNAL OF APPLIED》 *
杨柳: "基于环形阵列的超声层析检测方法" *
贺梦颖: "医学图像的增强技术研究与仿真", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110706298A (zh) * 2019-09-10 2020-01-17 天津大学 正则化加权最小二乘的透射反射双模式超声成像重建方法
CN110706298B (zh) * 2019-09-10 2023-04-07 天津大学 正则化加权最小二乘的透射反射双模式超声成像重建方法
CN112067691A (zh) * 2020-08-31 2020-12-11 天津大学 油气水三相介质衰减谱融合多频超声层析成像方法
CN112067691B (zh) * 2020-08-31 2024-01-30 天津大学 油气水三相介质衰减谱融合多频超声层析成像方法
CN113092049A (zh) * 2021-03-25 2021-07-09 北京理工大学 一种三维跨界面成像方法
CN113092049B (zh) * 2021-03-25 2022-04-15 北京理工大学 一种三维跨界面成像方法

Also Published As

Publication number Publication date
CN110097608B (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
CN112067691B (zh) 油气水三相介质衰减谱融合多频超声层析成像方法
US8614707B2 (en) 3D and real time electrical capacitance volume-tomography sensor design and image reconstruction
CN110097608A (zh) 修正路径追踪描述的连续波超声层析成像重建方法
CN101794453B (zh) 基于回归分析的节点映射图像重构方法
CN106442631B (zh) 基于电学/超声双模态融合的层状界面重建方法
CN109598769B (zh) 总变差正则化约束的超声成像同步代数迭代重建方法
CN104688224B (zh) 一种应用于声学非均匀媒介磁声耦合成像重建方法
CN102608205B (zh) 基于变波速相位迁移的多层物体无损检测超声成像方法
Liu et al. Real-time reconstruction for low contrast ultrasonic tomography using continuous-wave excitation
CN109884183A (zh) 透射反射模态融合的超声层析成像方法
Liu et al. Nonlinear ultrasonic transmissive tomography for low-contrast biphasic medium imaging using continuous-wave excitation
Jiang et al. Contactless electrical impedance and ultrasonic tomography: Correlation, comparison and complementarity study
CN110706298B (zh) 正则化加权最小二乘的透射反射双模式超声成像重建方法
Wang et al. Shape reconstruction for electrical impedance tomography with V 2 D-Net deep convolutional neural network
CN109102552B (zh) 一种非均匀形状约束的像素值域滤波超声成像重建方法
Liang et al. Ultrasound guided electrical impedance tomography for 2D free-interface reconstruction
CN109946388B (zh) 基于统计逆的电学/超声双模态内含物边界重建方法
Li et al. A novel gas/liquid two-phase flow imaging method through electrical resistance tomography with DDELM-AE sparse dictionary
Hou et al. Multi-frequency ultrasound tomography based on modified matrix regularization method and wavelet fusion
Liang et al. Dynamic imaging for time-varying distribution using electrical/ultrasonic dual-modality tomography
Zhang et al. Survey of EIT image reconstruction algorithms
CN113870377A (zh) 基于V-ResNet肺部成像方法
Huang et al. ECT image reconstruction method based on multi-exponential feature extraction
Willemink et al. Imaging of acoustic attenuation and speed of sound maps using photoacoustic measurements
CN108254418A (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