超声非线性成像-谐波成像基本知识

2023-11-18

概述

  最近了解了一下非线性超声中的谐波成像,可以作为知识普及了解一下。大部分来自于文献,如有遗漏可以私信咸鱼进行修改。另外有错误和补充的也可以私信咸鱼
  超声成像因其相对较低的价格、无电磁辐射、无创伤、成像方式多样而在医学成像中得到了广泛的应用。早期的超声成像技术都是基于线性声学原理。然而在实际的医学超声成像应用中,声波在生物组织的非线性传播和造影剂的非线性振动能够产生非线性声信号。在接下来的报告中讲解了如何使用非线性信号中的谐波信号进行成像,以及相关的谐波信号提取方法。最后使用一些代表性的应用实例进行验证,得出谐波成像在分辨率和对比度上确实优于传统基波成像的结论。

超声的非线性

  当声波在均匀介质中传播时,基本上按直线传播,无任何反射。但当声波入射组织中时,由于界面的声阻抗的差异而被反射,或者遇到相对波长较小的颗粒被散射,脉冲回波式成像就是利用反射和散射的回波信号强度作为成像参量成像。成像基础是生物组织的声阻抗率的差异,图像反应了声阻抗率的分布,也就可以显示出组织的结构。
  线性超声是建立在小振幅的假设基础上的,该假设要求声压、质点的振动速度、质点的振动位移、密度的起伏分别相对于静压强、声波的传播速度、声波的波长、静密度为小量,也就是常说的远小于。但是在实际中换能器发射的声波幅度较大,频率较高,此时小振幅假设就不再成立,必须考虑非线性效应。同时,由于生物组织的不均匀性引起的非线性,还有造影剂中微气泡的非线性振动产生的非线性信号。
  在线性超声成像中,非线性是无用的,甚至是负面的。九十年代后期,开始利用由生物组织和超声造影剂的非线性效应而产生的谐波信号来成像,即谐波成像技术。通过大量仿真和实验研究表明,谐波成像有更好的空间分辨率和对比度。但是谐波的强度比基波弱很多,信噪比低。根据谐波产生的来源不同,谐波成像可以划分为组织谐波成像(Tissue Harmonic Imaging)和造影剂谐波成像(Contrast Harmonic Imaging)。

组织谐波成像

  组织谐波主要由入射波的非线性传播引起。中心频率为f0的声波发射到组织中,能得到2f0、3f0等谐波信号,分别称为二次谐波、三次谐波等,如下图所示。三次谐波以及后面的谐波统称为高次谐波。
谐波的产生
  当声波在组织中传播时,会在组织中形成压缩区和疏松区,两个区域的传播速度不同。如下图所示,由于各点的传播速度不同,会导致波形发生畸变,从而产生谐波。随着声波在组织中进一步传播,波形畸变越来越大,谐波也越来越丰富。
波形畸变和谐波滋生
  声波在生物组织中传播时,基波的强度会随着深度的增加而越来越小。谐波的强度会随着深度的增加而越来越大,到后面由于衰减,强度会减小,如下图所示。
基波和谐波声压随深度的变化
  组织谐波成像需要滤除回波中的基波信号,并提取组织产生的谐波信号来进行成像。如果换能器发出的声脉冲包括一个宽频带的频率成份,接收到的信号带宽就呈现出发射带宽的倍数,包括基波或入射带宽和组织产生的谐波的带宽。如果不控制基波带宽,那么接收到的基波和谐波的带宽重叠,谐波信号就会被高强度的基波所淹没,造成谐波图像的质量下降。

组织谐波的提取方法

  一般有三种方法进行组织谐波的提取,分别为带通滤波成像、脉冲反转谐波成像、脉冲幅度调制谐波成像。
  带通滤波成像是将发射的基波信号的频带控制的比较窄,那么产生的二次谐波和基波之间会分开处在各自的频带上,再使用带通滤波器滤除基波成份就可以提取出二次谐波成份,如下图所示。由于该方法需要将基波信号的频带控制的比较窄,由傅里叶变换可知,频带越窄,脉宽越宽,这导致轴向分辨率的降低。
基波和谐波频带分开
  脉冲反转谐波成像通过连续发射两次脉冲,而脉冲的频率一致,相位相差180度。将接收到的两个脉冲的回波信号相加。设接收到的两次响应分别为:
公式1
两者相加得:
在这里插入图片描述
  则基波和奇次谐波相互抵消,偶次谐波变为原来的两倍,得到增强,如下图所示。
脉冲反转成像方法
  脉冲反转谐波成像增加了信号强度,提高了信噪比。在发射的时候,可以采用宽带发射,增加了轴向分辨率。但是该方法需要进行两次脉冲的发射和接收,降低了帧频。而且对运动特别敏感,容易产生运动伪像,对心脏等器官成像尤为明显,需要一些应对措施。
  脉冲幅度调制谐波成像通过先后发送两个幅度不同,频率和相位都相同的脉冲。例如第二个脉冲幅度是第一个脉冲幅度的1/4。在接收的时候将第二个脉冲的回波信号乘以4。然后使用接收到第一次脉冲的回波信号与第二次乘以4之后的信号相减。由于第二次的发射幅度比第一次小得多,所以非线性效应较小,谐波分量不明显,而基波分量几乎是一样的。相减之后,基波分量被抵消,谐波分量得到保留。
  以上三种方法中,脉冲反转谐波成像不仅能很好地消除基波,同时还能增强谐波信号,相对有更好的成像效果。

组织谐波成像的优点

  组织谐波成像相对于基波成像有许多优点,重点介绍一下消除近场多重反射的伪像,抑制旁瓣回声以及能检测更深的组织这三个优点。
  如下图所示,当遇到两个平行的强反射界面时,发射的声波可能经过多次反射产生伪像。但是使用组织谐波成像时,由于谐波随着传播深度先逐渐增强,然后逐渐减弱。在探头与皮肤之间非线性效应不明显,几乎没有谐波,所以使用组织谐波成像可以消除近场多重反射造成的伪像。
近场多重反射
  谐波只有在基波的较高声压处才会产生,弱的基波几乎不产生谐波,这使得谐波相对于基波有更窄的主瓣、更低的旁瓣,如下图所示。而更窄的主瓣代表着更高的空间分辨率,较低的旁瓣代表的就是较高的图像对比度,所以组织谐波成像具有更高的分辨率和对比度。
更窄的主瓣,更低的旁瓣
  谐波是随着在组织中的传播距离越深,谐波声压越大,之后才会衰减。这使得谐波声压的最大值比基波声压的最大值的深度要深。同时谐波只有在回波的时候才会有衰减,所以衰减也比基波小,因此组织谐波成像提高了深层组织的清晰度。

造影谐波成像

  超声造影剂大多是包含直径在1~10um微气泡的液体,可以通过静脉注射流到身体的各部分(肠道相关的可以通过口服的方式)。由于微气泡的声阻抗特性和组织存在很大差异,可以增强基波超声成像的诊断能力。如下图所示。
造影剂基波和谐波对比
  声波通过造影剂中的微气泡产生非线性传播,波形畸变,谐波成份明显增加。相比之下,其它组织产生的谐波相对较少。利用微气泡的非线性特性来增加造影剂和组织的对比度,提高信噪比。但是以上都是建立在组织回波都是线性的假设上,现在确定组织也是可以产生源于成像的谐波信号的,这个就会造成很大的干扰。
  当声压达到一定阈值,发射波的频率两倍于微气泡的共振频率时,微气泡除了产生二次谐波和高次谐波之外,还有次谐波(f0/2,f0/3…)和超谐波(3f0/2,5f0/2)。如下图所示。发射频率两倍于微气泡的共振频率这个条件是为了达到比较好的效果,如果频率达不到这个条件,还是会产生次谐波和超谐波,但是效果不明显。而组织不会产生次谐波和超谐波。这样就进一步提高了对比度。同时由于发射频率高,可以提高图像的轴向分辨率。对于接收的低频率的次谐波信号,回波衰减较小,有利于深部组织的检测。
造影剂产生次谐波和超谐波

应用实例

  通过一些应用中的实例对比,看一下谐波成像的效果。
  下图为组织谐波成像增强心内膜清晰度和近场分辨率,图a为基波成像,图b为组织谐波成像。可以看出组织谐波成像的心尖四腔视图更加明显。
心尖四腔成像基波和谐波对比
  下图为为胸骨旁心脏长轴图像, (a)基础成像, (b)谐波成像和©脉冲反转谐波成像。可以看出组织谐波成像减少了伪影,而脉冲反转谐波成像进一步提高对比度和分辨率。
基波成像和谐波成像
脉冲反转成像
  下图是通过二次谐波成像增强造影剂强度。左图为基频静脉造影剂成像,右图为二次谐波造影剂成像。可以看出使用二次谐波造影剂成像的对比度明显增强。
造影剂基波和谐波成像
  下图中使用的是将8mm血管嵌入到模拟组织的仿体进行的实验。图a是造影前的基波图像,图b是造影后的基波图像。图c是造影前的次谐波图像,图d是造影后的次谐波图像。图e是造影前的谐波图像,图f是造影后的谐波图像。可以看出造影之后的图像对比度明显比造影前的高。而使用了造影剂之后,次谐波比谐波的图像对比度要好,而且更加清晰。使用造影剂之后的基波图像对比度要比谐波的差一点,而且有伪影。
造影前后的基波,次谐波,谐波成像对比

总结

  通过原理和应用实例的比较,相比于谐波成像,谐波成像有着更好的图像分辨率和对比度,而且可以检测到更深的组织图像。但是在实际的应用中,由于谐波的信噪比相对于基波低很多。这就要求设备的灵敏度非常高,增加了设备开发实现的难度,提高了成本。
  在报告中提到的伪影和噪声都是与超声和生物组织相关的。其实在提取信号的时候,设备中的噪声也不可忽视。尤其对于信噪比比较低的谐波信号,不是说探头的灵敏度达到就可以的,电源的要求也会更高。在基波超声设备中,由于基波信号较大,那么相应的电源纹波要求会较低;而在谐波超声设备中,由于谐波的信噪比较低,信号强度也比基波弱很多,那么相应的电源纹波要求就会高很多,否则谐波信号会被淹没在电源噪声中。同时内部电路的电磁兼容性要求也会相应的提升,这些都是谐波成像中必须考虑的问题。

  以上就是关于谐波成像的一些了解,有问题可以在下方留言。
更多文章请关注微信公众号:ubug404
公众号

参考文献

[1] Martin K. Introduction to B-mode imaging[J]. Diagnostic Ultrasound: Physics and Equipment, 2010:1.
[2] 杜功焕,朱哲民,龚秀芬,“声学基础”,南京大学出版社,2001.
[3] 钱祖文,“非线性声学”,科学出版社,2009.
[4] 安普岩.基于脉冲反相技术的组织谐波成像研究[D].辽宁:东北大学,2015.
[5] James D. Thomas, David N. Rubin,Tissue harmonic imaging: Why does it work?,Journal of the American Society of Echocardiography,Volume 11, Issue 8,1998,Pages 803-808.
[6] Peter N Burns, David Hope Simpson, Michalakis A Averkiou, Nonlinear imaging, The Future of Biomedical Ultrasound, Volume 26, Supplement 1, S19-S22, May 01, 2000.
[7] F. Forsberg, W.T. Shi, B.B. Goldberg, Subharmonic imaging of contrast agents, Ultrasonics, Volume 38, Issues 1–8, 2000, Pages 93-98.
[8] 章东,龚秀芬,马青玉,非线性超声医学成像的研究进展,百度文库:https://wenku.baidu.com/view/b0842b6a0a4e767f5acfa1c7aa00b52acec79c14.html

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

超声非线性成像-谐波成像基本知识 的相关文章

  • 理解SPDX

    SPDX The Software Package Data Exchange SPDE an open standard for communicating software bill of material information in
  • 基于RFID技术在物流仓储中的解决方案—铨顺宏FUWIT

    一 行业背景 2011年6月8日国务院出台物流行业新国八条 明确指出要推进物流技术创新和应用 加强物流新技术自主研发 加快仓储物流设备研制 制定和推广物流标准 适时启动物联网的应用示范 推进物流信息资源开放共享 物流信息化 指在物流活动中全
  • 基于Opencv的水位识别,液面识别、高度识别

    Update 代码已经上传到github上了 可以点这里 Cutting 一直说这要整理一下Computer Vision课程的大作业 拖了好久 这两天忙着写一个订单处理的第三方库 陷入了僵局 所以换个口味 把大作业整理一下 Require
  • python入门到精通 _7异常、模块与包

    文章目录 1 异常捕获 1 1 捕获常规异常 1 2 捕获指定异常 1 3 捕获多个异常 1 4 异常的finally 2 模块的导入与自定义 2 1 模块的导入 2 2 模块的自定义 3 安装第三方包与自定义 3 1 自定义包 3 2 安
  • 生成3d地图obj(可导入到ppt中使用)

    一 效果 二 使用 1 网址 https 3d mapper com MAP index php 2 注册 3 使用 创建3D地图 正在生成中 有以下设置功能 截图 管理你的地图
  • 按键松手检测 - 检测是否连续按下

    u8 KEY Scan void static u8 keyup 1 防止检测多次 if keyup KEY0 0 KEY1 0 KEY3 0 delay ms 50 去抖 if KEY0 0 KEY1 0 KEY3 0 keyup 0 i
  • 代码质量有哪些评判标准?

    描述代码质量的词 灵活性 flexibility 可扩展性 extensibility 可维护性 maintainability 可读性 readability 可理解性 understandability 易修改性 changeabili
  • Kruskal算法代码实现

    package com kruskal import java util Arrays public class KruskalCase private int edgeNum 边的个数 private char vertexs 顶点个数
  • 树莓派+8x32的MAX7219点阵模块显示字符串

    https blog csdn net zyaiwmy article details 70224250 utm source blogxgwz1 按照上面这篇文章进行配置了驱程 修改了一下示例代码 效果如下 这里记一下项目的github地
  • Entity Framework Core系列教程-2安装EF Core

    安装Entity Framework Core 这里我们将使用EF Core 3 1 因为它是长期支持版本 开发工具使用Visual Studio 2019 编写 NET Core应用程序 EF Core 3 1 支持 NET Standa
  • 华为机试题-字符串最后一个单词的长度

    计算字符串最后一个单词的长度 单词以空格隔开 字符串长度小于5000 输入描述 输入一行 代表要计算的字符串 非空 长度小于5000 输出描述 输出一个整数 表示输入字符串最后一个单词的长度 示例1 输入 hello nowcoder 复制
  • web前端笔记:html5的标签

    在HTML4 01中 lt b gt lt i gt 是视觉要素 presentationl elements 分别表示无意义的加粗 无意义的斜体 表现样式为 font weight bolder 仅仅表示 这里应该用粗体显示 或者 这里应
  • 人脸重建环境配置时的坑

    ERROR Could not install packages due to an EnvironmentError Errno 28 No space left on device 由于是docker环境 所以把 tmp目录下的文件清理
  • 《Apache MINA 2.0 用户指南》第一章:入门

    最近准备将Apache MINA 2 0 用户指南英文文档翻译给大家 但是我偶然一次百度 发现 Defonds 这位大牛已经翻译大部分文档 原文链接 http mina apache org mina project userguide c
  • [C#] 汉字转拼音,支持多音字

    这份代码大概不是严格意义上正确的 但是一般场景用用应该没问题 而且支持dotnet core public static class Pinyin region 拼音对照表 private static string py 一 yi 丁 d
  • Unity 调用系统键盘

    using System Diagnostics using UnityEngine public class KeyboardEvent MonoBehaviour private Process keyboard
  • 第十章 内部类(下)

    第十章 内部类 下 随着后面所讲的内容越来越深入 所以可能理解得比较慢了 同时这里边的文字描述和示例也越来越多 希望大家能够坚持下去 慢慢看完 相信会有所收获 当然如果像我一样一个字一个字的敲一遍 印象会更深的 不骗你 因为有的东西开始真的
  • SpringCloud Alibaba 框架背后的故事

    前言 Spring Cloud Alibaba是Spring Cloud的一个子项目 它是由阿里巴巴公司推出的 用于构建基于微服务架构的分布式应用程序的开源框架 它与Spring Cloud的其他组件 如Netflix OSS 相结合 为开
  • sqli-labs Less5-6(布尔盲注)

    目录 前言 一 Less 5 1 布尔型的注入相比前四关 思路上最大的不同就是通过对错来获取对你来说有用的信息 1 先找到注入点 2 先判断数据库长度 3 判断数据库名中的每一个字母是什么 4 同样得方法测试表名 5 根据第四步得出表名 推

随机推荐

  • 人性的弱点

    附 本作品来自互联网 本人不做任何负责 内容版权归作者所有 人性的弱点 by Dale Carnegie雷吟译 目录 这本书对你有十二种功用 译者序 前言 原著序 如何从这本书里获得最大效益 第一篇 待人的基本技巧 第一章 如欲采蜜 勿蹴蜂
  • QT槽函数的使用

    QT槽函数的使用 例如 在头文件中设置槽函数 public slots void OnClickedButtonEnsure void OnClickedButtonExit cpp文件 构造函数中写入 connect ui pushBut
  • STM32F103ZET6【标准库函数开发】------配置定时器参数的几个常用函数

    TIM TimeBaseInitTypeDef 基本初始化 TIM OCInitTypeDef 比较输出初始化 TIM ICInitTypeDef 输入捕获初始化 TIM BDTRInitTypeDef 断路和死区初始化 TIM TimeB
  • pycharm创建的虚拟环境为什么用conda env list命令查询不到?

    问题描述 pycharm创建的虚拟环境为什么用conda env list命令查询不到 pycharm开发环境可以创建虚拟环境 目的是为隔绝其他环境种库带来的版本干扰 但是发现一个问题 无论是在windows终端 anaconda终端 Py
  • Java课程设计-学籍信息管理系统

    一 系统分析 学生的学籍信息是记录学生的重要信息档案 如何以电子文档形式记录下学生的学籍信息是每个学校必须做的事情 该学生学籍信息管理系统就是为了方便学校记录下每一个学生的基本信息 生成电子数据库 并且能够做到查询 更改 删除 浏览等功能操
  • Unity 移动端触摸屏操作

    Unity 触屏操作 当将Unity游戏运行到IOS或Android设备上时 桌面系统的鼠标左键可以自动变为手机屏幕上的触屏操作 但如多点触屏等操作却是无法利用鼠标操作进行的 Unity的Input类中不仅包含桌面系统的各种输入功能 也包含
  • 关于类和对象

    1 面向对象是一种技术的开发模式 但是最早的时候所使用的模式是面向过程 面向对象就是一种组件化的设计思想 面向过程 指的是针对于某一问题单独提出解决方案及代码开发 面向对象 以一种组件化的形势进行代码的设计 这样开发出来的代码有一个最大的好
  • Python Tkinter实现简单文本编辑器

    下面看看如何利用Tkinter实现一个文本编辑器 语言 python 运行效果 在 Python 代码中使用 Tkinter 的简单文本编辑器应用程序 这 使用Tkinter的简单文本编辑器应用程序 是编码在 python 程序设计语言 该
  • [编写高质量代码:改善java程序的151个建议]建议103 反射访问属性或方法时将Accessible设置为true...

    转载于 https www cnblogs com akingseu p 3506576 html
  • canvas小练习之鼠标粒子特效

  • Qt 笔记5--Qt 拉伸窗体

    Qt 笔记5 Qt 拉伸窗体 正常情况下 Qt主窗体拉伸 但内部子窗体大小及位置并未发生该表 导致出现一部分空白区域 因此实际中为了美观需要实现其动态拉伸和调整效果 以下为笔者通过重新resizeEvent方法实现的一个简单拉伸案例 相应学
  • 小程序如何进行分包详细介绍

    微信小程序开发过程中 随着业务不断迭代 程序包的体积越来越大 使用分包加载是开发者必须面对的问题 正常情况下 小程序首次启动时 会将整个代码包下载下来 所以如果代码包过大 会影响小程序首次启动时间 因此微信官方对小程序代码包做了大小限制 官
  • JetBrains Resharper 2019.1.3的安装、破解、使用

    最近在浏览GitHub时 看到好多项目下面会有如下图标 几经折腾之后 发现这个一个非常好用的VS插件 称为神级也不为过 可以进行编写提示 代码重构 文件搜索等许多操作 当然目前最新版的VS中已经具有很多类似的功能 但是这个插件从实用性上来说
  • 2014年2月7日星期五(7-3,消除了背面的3D线框立方体)

    好久没进行了 看看吧 与DEMO7 2的区别 只看不同的地方 相机类型不同 归一化平面改为2 2 ar 视口为640 480 执行了背面消除 就是判断面元向量与面元到视点向量的点乘 gt 0 则可见 否则 消除 这一步是在物体剔除后和世界坐
  • 百度飞浆零基础实践深度学习入门课,通过PaddlePaddle实现波士顿房价预测

    百度架构师手把手带你零基础实践深度学习 课程已加入学习 对深度学习有了深入浅出对理解 并通过两种方式实现波士顿房价预测的算法 包括Python和Numpy构建深度学习框架 以及使用paddlepaddle快速搭建 快速实现预测 A 使用Py
  • 使用pl.droidsonroids.gif.GifImageView在安卓中显示动图遇到的问题

    在做一款聊天软件 其中聊天界面需要发送表情 而表情都是动图 在安卓中想要显示动图 就要借助第三方框架 我选的是pl droidsonroids gif GifImageView 使用方法如下 你在github上下载了android gif
  • Webpack5 处理图片资源、图标资源、修改输出目录、清空上一次打包的所有文件

    1 处理图片资源 Webpack5 已经将两个 Loader file loader 和 url loader 功能内置到 Webpack 里了 我们只需要简单配置即可处理图片资源 优化 gt 小图片转化为base64 小图片可以转换为ba
  • android分析monkey的anr

    1 根据anr信息和traces txt判断anr是发生在那个进程中 2 在文件夹logs4android x中搜索 ANR 或者在monkeylog中搜索 anr in 找到发生anr的时间 3 在文件夹logs4android x中搜索
  • 基于深度路由的水下传感器网络(翻译、整理)

    基于深度路由的水下传感器网络 翻译 整理 背景与意义 无线传感器网络引起了研究人员的广泛关注 并在该领域取得了很大进展 起初 这些网络仅涵盖地面的实用程序 但后来 由于在地球被 70 水覆盖 所以大家开始研究传感器网络的水下应用 也就是我们
  • 超声非线性成像-谐波成像基本知识

    文章目录 概述 超声的非线性 组织谐波成像 组织谐波的提取方法 组织谐波成像的优点 造影谐波成像 应用实例 总结 参考文献 概述 最近了解了一下非线性超声中的谐波成像 可以作为知识普及了解一下 大部分来自于文献 如有遗漏可以私信咸鱼进行修改