医学超声成像—合成孔径法(Synthetic Aperture Methods)[2]
医学超声成像—合成孔径法(Synthetic Aperture Methods)[2]
- 3. 医学超声成像(Medical Ultrasound Imaging)
- Transducer and Arrays(传感器阵列)
- 信号处理算法
- 成像质量评估参数(Image quality Factors)
Reference:Yao H . Synthetic Aperture Methods forMedical UltrasonicImaging[J]. 1997.
3. 医学超声成像(Medical Ultrasound Imaging)
- Ultrasound imaging: non-invasive nature and apparent safty
- 超声成像系统:transducer、scanning method、scan conversion technique
- 超声成像评价指标:空间分辨率(spatial resolution)、对比度(contrast resolution)、时间分辨率(temporal resolution)、信噪比(signal-to-noise ratio)
Transducer and Arrays(传感器阵列)
- 压电传感器(piezoelectric transducer):将电信号转化为机械信号,将机械信号转化为电信号;可以发射短脉冲、接收回声
- 线性阵列:用于腹部(abdominal imaging)造影成像等
- 环形阵列:用于血管(intravascular imaging)造影等
- 压电传感器材料
衡量压电材料指标是ECC(Electrical-to-mechanical coupling coefficient):
锆钛酸盐(PZT,lead zirconate titanate),ECC较高,但声阻抗大
聚偏氟乙稀(PVDF,polyvinylidene fluoride )
- 单元传感器与阵列传感器
单元传感器(single element):控制简单,不需要beamforming
多元传感器(multi-element):能够提高信噪比,包括线性阵列、环形阵列、2D阵列
环形阵列(annular array):将球形单晶切割成同心环来制造环形阵列。每个环相互独立,有自己的处理电路,一般有2-8个环。为每个环加相应的延时,改变接收时间,得到了一个可变的电子透镜。聚焦平面关于线性阵列对称。
线性阵列(linear phased array):单元沿直线等间距排列,每个都可以发射接收声信号,单元间隔小于λ/2\lambda/2λ/2
曲线阵列(curved linear array):单元排列在一个曲面上,可以更好地与被测曲面接触
2-D阵列:所需要的单元多,线性阵列需要32-128个单元,2-D阵列需要32×32,128×12832\times 32, 128\times 12832×32,128×128。但是有利于3-D成像
信号处理算法
超声成像系统示例图:
Transduer:发射声波信号,接收声波信号
RF-units:负责聚焦、位移等,得到不同波形、方向的超声波,对接收信号做beamforming(波束形成)
过程:将超声波束聚焦到一个方向,接收得到RF数据,聚焦求和得到单线(single line)RF数据,再移到下一个方向
扫描方法:
通常有两种扫描方法:机械扫描、电子扫描。扫描方向包括线性扫描、扇形(sector)扫描
1. 机械扫描:通过传感器搭载的电机移动等,线性扫描通常需要较大的声学窗口(acoustic window),扇形扫描需要较小的
2. 电子扫描:通过传感器阵列实现。对于扇形扫描,通过传感器单元电信号的延时实现不同方向的扫描,电信号延迟函数是关于单元位置的线性函数,为避免栅瓣(grating lobes),内部单元之间的间距需要小于λ/2\lambda /2λ/2;对于线性扫描,逐个或者逐组激发,单元间距可大于λ/2\lambda /2λ/2。
脉冲回波B模式成像(Pulse- Echo B-mode Imaging)
基于反射、散射回波成像
相位调制多普勒方法(phase demodulated Doppler methods):对于移动物体,当物体靠近观测源时,频率会变高,反之,频率变低。
A模式(Amplitude mode):最古老最简单的脉冲回波超声,回声幅度是关于穿透深度的函数
B模式(Brightness mode):将A模式中的幅度信息转化为亮度信息,可以得到2-D扇形扫描区
C模式(Constant depth mode):与B模式类似,区别在于扫描平面固定在距探头特定深度处
M模式(Motion mode):对于运动物体,一个轴表示穿透深度(纵向时间),另一个轴表示时间(物体的空间位移)
扫描转换(Scan Conversion)
在B模式下,短脉冲波按照特定方向像人体传播,在不同组织(具有不同波阻抗)交界面被反射,传感器接收回波,可以重构出被测物体的切面图。得到极坐标图,纵轴为范围(深度)rrr,横轴为扫描角度θ\thetaθ。Scan conversion将数据从极坐标图变为笛卡尔坐标图。
笛卡尔坐标与极坐标的转换关系为:
进一步通过插值得到笛卡尔坐标系的图像。
插值
为降低时间成本,采用最近邻插值或者双线性插值。通常采用双线性插值:
插值公式为:
其中,点S1,S2,S3,S4S_1,S_2,S_3,S_4S1,S2,S3,S4的极坐标分别为:(r1,θ1),(r2,θ2),(r3,θ3),(r4,θ4)(r1,\theta 1), (r2, \theta 2), (r3,\theta 3), (r4, \theta 4)(r1,θ1),(r2,θ2),(r3,θ3),(r4,θ4)。
增益补偿(Gain Compensation)
超声在传播过程中会衰减,在不同深度的相同探测目标会有不同的回波幅值。- TGC(time-gain compensation):增益补偿是关于深度的线性函数。利用超声在介质中的衰减公式计算补偿公式即可。
对数压缩(Logarithmic compression)
降低超声图像的动态范围,降低强弱信号之间的显示差距,使得弱信号更明显。对数压缩公式:
- TGC(time-gain compensation):增益补偿是关于深度的线性函数。利用超声在介质中的衰减公式计算补偿公式即可。
成像质量评估参数(Image quality Factors)
最常用的评估参数:空间分辨率、对比分辨率、时间分辨率、信噪比
空间分辨率(Spatial resolution)
空间分辨率指两个反射点之间可以被区别的最小距离,用于衡量分辨小结构、边缘的能力。对于2-D成像,空间分辨率包括轴向分辨率、横向分辨率。
轴向分辨率(axial resolution)
沿着声波传播轴的分辨率,取决于脉冲的长度TpT_pTp(与传感器带宽BwB_wBw反相关),轴向分辨率为:
c为声波传播速度。传感器频率越高,脉宽越小,轴向分辨率越高。但是高频率会导致声波衰减变大,需要在这两者之间做权衡。
横向分辨率(lateral resolution)
垂直于声波传播轴的分辨率,主要取决于波束(beam)宽度,其次是孔径特征(由单元数量决定)、传感器频率与调焦。对于电子相控阵波束,横向分辨率为:
横向分辨率取决于主瓣角宽度θml\theta_{ml}θml,与传感器表面的距离RRR:
当角宽度较小时,公式可以简化为:
【问题】:角宽度θml\theta_{ml}θml如何确定?
当传感器主频增加时,波束宽度会减小,从而有更好的横向分辨率。
旁瓣问题:由阵列产生,会引入其他方向的信号噪声。栅瓣(grating lobe)是旁瓣的一种,当主瓣大角度转向时会产生栅瓣,会引入人工噪声。决定栅瓣的一个参数是两个相邻传感器单元的间距ppp(pitch),当满足一下条件时,栅瓣的影响可被完全消除:
其中,λmin\lambda_{min}λmin表示脉冲频谱中最高频率分量对应的波长。
对比度分辨率(contrast resolution)
对比度分辨率指最小能分辨的波阻抗。影响对比度分辨率的主要因素时空间分辨率、系统动态范围、旁瓣水平和图像噪声。空间分辨率影响高对比度的目标(波阻抗相差较大),图像斑点(speckle)噪声影响低对比度目标。
对比度分辨率也被称作图像动态范围,定义为最强信号与最弱信号之间的比值,通常以dB表示。
时间分辨率(temporal resolution)
时间分辨率定义为对运动物体的分辨能力,取决于系统成像时间(帧率:每单位时间可获取的图片数量)。获取一帧图像的时间TTT等价于处理一个波束的时间TbT_bTb与波束数量MMM的乘积。对于距离为R的物体,至少需要等待2R/c2R/c2R/c来获取信号,此时才能发射下一个波束。为了确保上一个脉冲已经全部接收完成,我们继续等待T0T_0T0才去发射下一个波束。帧率即周期TTT的倒数:
帧率在实时成像系统中时很重要的指标。
信噪比(Signal-to-noise Ratio)
表示有效信号与干扰噪声之间的比例。对于一个传感器单元,加噪图像满足:
s(t)s(t)s(t) 表示信号成分,n(t)n(t)n(t)表示噪声成分。
SNR满足:
当噪声分量与信号分量不相关时上式成立。噪声主要有两种类型:斑块噪声(speckle)、随机电子噪声。随机电子噪声主要由信号处理电路产生。散斑噪声(Speckle)是由来自许多小散射中心的声波的相干散射产生的。
对于M个单元的传感器阵列,SNR是原来的M倍。
总结
文章后面部分,使用electronic beamforming(电子扫描构成超声波束)和B-mode获得超声图像。由于使用了electronic beamforming,因此采用传感器阵列。
医学超声成像—合成孔径法(Synthetic Aperture Methods)[2]相关推荐
- 计算机在超声的应用,计算机在医学超声成像中应用.pdf
综 述- 计算机在医学超声成像中的应用 李义兵1余大昆1刘晓东2 (1.咸宁学院生物医学工程系湖北省咸宁市437100:2.襄樊市中心医院湖北省襄樊市441021) 摘要介绍了计算机在医学超声成像领域 ...
- 经典文献翻译:合成孔径雷达快速反投影算法(Fast Backprojection Algorithm for Synthetic Aperture Radar)
排版不方便,图片不好上传,Word版翻译及原文下载地址://download.csdn.net/download/u010112513/12576783 Ali F. Yegulalp 翻译:Ivan ...
- 超声成像_人工智能如何帮助转变医学超声成像
超声成像 介绍 (Introduction) Diagnostic ultrasound is a popular imaging method used for a variety of scree ...
- 合成孔径雷达影像变化检测研究进展
遥感影像变化检测技术用于检测同一地点在一段时间内所发生的变化情况,具有重要的应用价值,而基于合成孔径雷达(synthetic aperture radar,SAR)影像的变化检测由于其传感器具有不受时 ...
- 基于深度学习的合成孔径雷达自聚焦
文章目录 引言 什么是合成孔径雷达 什么是自聚焦 经典自聚焦方法 基于机器学习的方法 基于极速学习机的方法 基于深度学习的SAR自聚焦 代码附录 引言 本文全面介绍合成孔径雷达自聚焦概念和方法.想获取 ...
- matlab散斑噪声,合成孔径雷达图像的散斑噪声滤波方法与流程
本发明属于星载合成孔径雷达和光学图像信息处理领域,涉及一种新的合成孔径雷达图像的散斑噪声滤波方法. 背景技术: 合成孔径雷达(synthetic aperture radar,SAR)是一种高分辨的微 ...
- matlab 雷达拼图,(完整版)SAR合成孔径雷达图像点目标仿真报告(附matlab代码)
SAR 图像点目标仿真报告 徐一凡 1 SAR 原理简介 合成孔径雷达(Synthetic Aperture Radar .简称SAR)是一种高分辨率成像雷达技术.它利用脉冲压缩技术获得高的距离向分辨 ...
- 遥感相关专业英语词汇汇总
对网络上英语专业词汇进行汇总整理,如下表: 序号 英语基础词汇 摄影测量与遥感专业英语词汇 英语基础词汇 英语基础词汇-生态环境与遥感 1 Absorbed photosynthetically ac ...
- sar图像去噪matlab,一种基于总曲率的SAR图像变分去噪方法与流程
本发明属于数字图像处理技术领域,具体涉及一种基于总曲率的SAR图像变分去噪方法. 背景技术: : 相干斑噪声是合成孔径雷达(Synthetic Aperture Radar,简称SAR)图像的重要特征 ...
最新文章
- 成都Uber优步司机奖励政策(4月12日)
- android视频分享功能吗,Unity 调用 Android 分享功能(基于ShareRec SDK视频分享)
- 小程序服务器api接口,百度智能小程序API 开放接口
- 十二月无书稿,总个结(2011)
- 【追光者系列】Hikari连接池大小多大合适?(第一弹)
- python dataframe将字符转换为数字_python中如何将华氏温度转换为摄氏温度?
- IIS 配置 PHP 环境搭建:web 文件管理器
- 软件推荐:Picasa
- Springboot毕设项目公共机房的值班管理系统wyz7b(java+VUE+Mybatis+Maven+Mysql)
- 为什么国内抖音没有网页版,原因竟然是这样!
- uv-pv-vv的区别
- Core的.mli文件在哪里
- Mac虚拟机Net模式设置
- hexo添加文章更新时间
- Django分布式路由
- 一位拖延症患者程序员的自我救赎!
- 【社媒营销】四大跨境营销渠道分析
- oracle表给用户授权
- 多线程 - 三种实现
- LinuxC进程编程
热门文章
- oracle11g連不上em,oracle11gem重建失败的几点解决办法.doc
- 所有的环境都配置好了,将新创建的weex项目导入到Android studio中的时候,出现的错误。...
- BP客户主数据信用数据批量修改
- html中图片自动循环滚动代码,实现长图片自动循环滚动效果
- 烤仔的朋友们丨政策暖风吹来,国内公链们的春天来了?
- 算法很美 笔记 4.多维数组与矩阵
- Linux 磁盘划分 LVM 逻辑分区管理步骤
- 电脑怎么写入便签并同步到手机版便签上?
- java实现word模板文件填充
- “不确定性问题”(Nondeterministic Problem,NP)与哥德尔不完全定理“