超声中的信号处理ppt课件

超声中的信号处理 一 1 内容提要 概述超声前端波束合成信号处理 2 1 针对两年之内的新员工 2 初步理解超声各种模式的成像原理及过程 3 对超声图像如何得到 有一个明确的概念 概述 目的 3 好莱坞大片的飞机以及军舰上经常能看到上面的图像 雷达系统 发射无线电 通过检测回波的信息来判断有无目标以及目标的距离速度等信息 超声技术80 来源于雷达 脉冲波 连续波 合成孔径 相控阵 超声系统 发射超声波 通过检测从人体组织反射回来的回波信息 以图形化的方式再现人体内部器官的状况 供医生诊断 概述 超声系统与雷达的相似之处 4 线阵颈动脉血流速度图 概述 几幅超声图像例子 肾脏血流能量图 凸阵肝脏血流速度图 心脏黑白图 5 从最开始的扫描 到最后的计算机显示 上面这些图形是如何形成的 首先看看超声信号处理结构 超声处理结构图 前端处理 中端处理 后端处理 概述 超声处理流结构框图 6 前端处理 博大精深 是超声系统的精髓 发射聚焦 逐级驱动产生高压 高压开关 对发射来说是将N个逻辑通道连接到K N个物理阵元中的N个物理阵元 对接收来说是从K N个物理阵元中选出N个物理阵元输出到N个逻辑通道 其中M7系统的N 64 K 2 模拟前端 小信号放大 低噪放 VGA 抗混叠滤波 AD变换波束合成 延时RAM 聚焦计算 变迹 求和 概述 超声前端处理结构图 7 中端处理 噪声中必要的信号检测手段 正交解调 频移 抽取前的抗混叠滤波 抽取 B信号处理 包络检测 对数压缩 C信号处理 解交织 壁滤波 自相关 参量估计 D信号处理 壁滤波 频谱估计 声音信号处理 概述 超声中端处理结构图 8 后端处理 图像的美化以及显示 B信号处理 帧相关 图像增强 扫描变换 C信号处理 帧相关 KEYHOLE去噪 扫描变换 D信号处理 谱后处理 扫描变换 各种模式的融合 概述 超声后端处理结构图 9 一个女孩子整体形象的构成 先天条件 必要条件 后天弥补超声前端 先天条件 相貌 身材 肤色 超声中端 必要妆扮 如刷牙 洗脸 超声后端 弥补妆扮 胭脂 香水 减肥药 概述 超声前端 中端 后端关系 10 先天条件好 不需要中端 后端弥补手段 因此说超声前端是超声系统的精髓 概述 超声前端 中端 后端关系 11 内容提要 概述超声前端波束合成信号处理 12 发射波形 发射波形考虑的因素 频率 级性 脉冲串个数对图像清晰度 分辨率 灵敏度的影响 已超出我的所知 13 发射波形 第一 第二个波形的发射频谱 发射脉宽宽度和频谱带宽成反比 单级性 在零频附近很大的一部分的发射能量被浪费了 14 发射波形 双级性发射无零频 能量利用率高 单个双级性脉冲的发射带宽比多个双级性脉冲的发射要宽 常用于B M 黑白图像更强调的是清晰度 多个双级性脉冲的发射常用于C PW 相位检测图像更着重的是灵敏度 右边第三图是发射加窗后的频谱图 降低了旁瓣 主瓣有所展宽 15 发射波形 加HANNING窗后的线性调频 chirp 波 从频域上看 旁瓣更低 200dB以下 基本无旁瓣 优点1 有利于自然谐波成像 优点2 能量利用率更高 16 发射波形 传统的发射波形和加窗的线性调频波形成像效果对比 17 前端发射框图 18 前端接收框图 19 发射 接收 高压隔离细节 20 前端ATGC控制改善信号的动态范围 21 前端ATGC控制改善信号的动态范围 22 内容提要 概述超声前端波束合成信号处理培训课程新人培养 23 在发射焦点法线两侧各32个阵元 发射脉冲信号 各阵元的发射信号的先后顺序不同 使得64路信号同时到达焦点 形成发射聚焦 很明显 中间阵元距离焦点短 发射最晚 最两侧阵元距离焦点远 发射最早 发射聚焦 24 波束合成 接收线 25 假设发射焦点法线两侧有两条线 可以认为线上的每一点都是一个点光源 向四面八方反射 收集反射波到达不同阵元的等相位点数据 相加合并后为该条接收线的一个样点 上图中圆弧和延时RAM相交的位置 即为该点回波的等相位点在延时RAM中的位置 波束合成 接收线上不同点的等相位点在延时RAM中的位置 26 下一次发射 让聚焦点向右偏移一点距离 法线也向右偏移同样的距离 参与发射的阵元同样向右偏移 波束合成 下一条线的发射聚焦 27 波束合成 再次得到两条接收线的数据 28 发射 接收 如此反复下去 直到收集完毕形成一帧所需的所有线数据 将收集到的线数据经过后续的处理 以一帧图像的方式直观显示出来 波束合成 收集一帧图像所需的线数据 29 多波束的概念 在DC6中 认为焦点法线两侧各有一条线 两条线同一深度点的在延时RAM中的等相位点的位置稍有不同 波束合成后的值也不同 在M5中 认为焦点法线两侧各有两条线 四条线同一深度点的在延时RAM中的等相位点不同 因此波束合成后的值也不同 在M7中 可以做到4 8 12波束兼容 波束合成 多波束的概念 30 多波束接收 同一深度不同波束的点在延时RAM中等相位点的位置稍有不同 上图为双波束接收的例子 波束合成 多波束的概念 31 系统时钟频率Fc即为系统总的资源 假设波束个数为N 即采样率为Fc N DC6 Fc 80M 2波束 因此采样率 40M M5 Fc 120M 4波束 因此采样率 30M M7 上图中对于4波束来说 采样率为160M 4 40M 对于8波束来说 采样率为160M 8 20M 波束合成 M7兼容4 8 12波束 32 反射点到阵元的距离不同 所需的时间不同 回波数据先到达中间阵元 然后到达两边的阵元 要收集所有阵元的等相位回波数据 等最两边的阵元的数据到达后才可以合成 因此需要延时RAM来存储中间先到来的等相位点的数据 接收聚焦的定义 在不同通道的延时RAM中求出等相位数据的过程 波束合成 聚焦的定义 33 为什么要接收多个通道的等相位点的数据 实际每一个通道的数据既有信号也有噪声 信号的等相位叠加是幅度的叠加 而噪声是不相关的 噪声的叠加是能量的叠加 因此当通道数增加一倍时 信号的幅度变为原来的两倍 能量变为原来的四倍 即能量增加了6dB 而噪声能量增加了一倍 即3dB 因此理论上 通道每增加一倍 SNR增加 6 3 3dB 为什么不增加更多的通道 当通道增加到一定值后SNR增益增加有限 但是随之而来的硬件前端成本 波束合成的资源成本却成倍增加 功耗也随之增加 折衷考虑 DC6 64通道 M5 32通道 M7 64通道 波束合成 多通道的相干叠加 34 相干叠加的概念贯穿于超声信号处理的整个过程 波束合成 相干叠加提高信噪比 35 沿着接收线方向有很多个样点 通常是几千个样点 计算每一个样点所有通道在延时RAM中等相位点的位置 计算量非常大 因为相邻样点弧线与延时RAM的交点变化非常小 以前DC6 M5为了减少计算量 认为沿接收线方向相邻的若干个接收点的半径相同 M7聚焦 改进了聚焦的设计方法 沿着接收线的方向 每一个采样点都计算其在延时RAM中的等相位点的位置 称为逐点聚焦 波束合成 逐点聚焦与隔点聚焦 36 实际等相位曲线在延时RAM中的位置 不一定刚好落在整数个样点上 如图中黑色圈圈所示 如果就近取出RAM数据的方法 得到图中所示折线 即用折线来近似代替等相位曲线 聚焦精度等于采样率 波束合成 聚焦精度 37 假设蓝色曲线为采样前的回波曲线 采样后延时RAM中的值为黑色箭头所示 现在等相位曲线和延时RAM相交点为图中红色箭头所示 现在问题转化为由已知的采样点求出原模拟曲线上任一点 由采样定理可知 如果采样频率大于信号最高频率的两倍 可以由采样后的离散点求出原模拟信号上的任何一点处的值 公式为以SINC函数为权重的对所有离散点的加权求和 波束合成 插值计算提高聚焦精度 38 插值计算 也可认为是采样率的变换 后面在将降采样率滤波器的时候 还会有详细的介绍 插值精度的选择 可以认为是几倍的升采样 实际中用到的是4倍的升采样 即可以得到两个样点之间的三个等时间间隔的采样点 上图中红色所示 聚焦精度为1 4chip 波束合成 插值计算提高聚焦精度 39 4倍升采样插值后 插值后的位置如上图黑线所示 和理论上的圆圈点仍然存在允许的误差 波束合成 插值计算提高聚焦精度 40 两种滤波方法 传统的插值滤波 在MATLAB中的命令为 y resample p q x 如下图蓝色所示 大约30阶 对于64个通道来说 阶数太高 计算量太大 信号带宽窄 可以用陷波器的方法去掉三份镜像 如下图红色所示 如此设计的陷波器 阶数非常低 需要的乘法器数目少 稍微牺牲一点性能 来换取过大的计算量 波束合成 插值实现方法 41 波束合成 插值前后频谱 42 波束合成 陷波器的简单实现方法 用极零点的方法来实现陷波器 43 用极零点设计滤波器的方法 得到该FIR陷波器的8个系数 有意思的是 该FIR的8个系数并非对称 但是相频曲线却是线性相位的 波束合成 陷波器的简单实现方法 44 波束合成 插值提高聚焦精度的例子 45 实际上由于探头的原因 不可能收集到所有的等相位点的数据 右图中只收到一部分等相位点的数据 图中黑点所示 等相位曲线上橙色点无法接收到 也可以认为都接收到了 但是权重不同 黑色点的权重为1 橙色点的权重为0 波束合成 变迹 46 将弧线展开为直线 等相位曲线上所有有用点加的窗如上图所示 因为不能收集所有的有用数据 加窗不可避免 时域加窗 即时域和数据相乘 对应于频域卷积 卷积的结果 是使得理想频谱被展宽了 而且出现了旁瓣 在数字信号处理的窗函数的截断效应里 有详细的论述 波束合成 加权的概念 47 波束合成 不同的窗函数 幸好前人专门研究如何降低截断效应引起的旁瓣 开发了各种其它的窗函数 HANNING HAMMING BLACKMAN 48 局部放大幅频特性曲线 可得下列结论 矩形窗的旁瓣的最高值是 13 26dB 而HANNING窗的旁瓣的最高值是 31 45dB 即HANNING窗对旁瓣的抑制比矩形窗好 HANNING窗的旁瓣下降趋势比矩形窗的下降趋势更快 HANNING窗的主瓣比矩形窗的主瓣宽 波束合成 不同窗函数的优缺点 49 变迹 给不同通道的等相位点数据加窗的过程 从右图成像结果可以看出 矩形窗旁瓣的影响 HANNING窗效果好很多 缺点 加HANNNING窗后 亮点略有展宽 任何一个处理过程都有利有弊 关键看利弊的大小 如何取舍 波束合成 变迹的定义 50 右图是变迹曲线随深度变化的示意图 远离中心阵元接收到的等相位点的数据 认为其作用小 离中心阵元越远 给于的权重越小 甚至权重取0 如图中单根红色曲线所示 随着深度增加 同一接收通道的权重逐渐增加 如图中红色曲线和黄线的相交点 图中最窄的变迹曲线是最近场的变迹曲线 波束合成 变迹曲线随深度而变 51 波束合成小结 波束合成 数据延时 聚焦计算 变迹 求和 52 内容提要 概述超声前端波束合成信号处理 53 内容提要 信号处理 检波信号处理 正交解调信号处理 B信号处理信号处理 C信号处理信号处理 多普勒信号处理信号处理 声音处理后端信号处理 54 检波 回波表达式 波束合成将多个通道 32 64等 的回波数据合并为一路数据 其特性和单通道的回波数据类似 回波不仅在幅度上受人体器官的调制 而且相位上由于多普勒效应也受到运动组织或者血流的调制 55 回波表达式 对于B M信号处理 仅关心回波的包络信息A t 对于C CM PW CW不仅关心回波的幅度 而且关心回波的相位信息 t B M恢复包络的方式 1 传统的绝对值后滤波的包络检波 黑白超声 DC6的B处理 2 希尔伯特变换法检波 产品中没有用到 3 正交解调的方法得出包络 M5 M7 检波 回波表达式 56 早期的黑白产品 以及第一款彩超DC6 对于B M信号的包络检测 都是用对信号取绝对值后低通滤波以恢复原信号的包络 检波 包络检波1 57 也可以采用HILBERT变换的方法 HILBERT变换的本质是一个90度的移相器 检波 包络检波2 58 上图为用HILBERT变换的方法求回波的包络 检波 包络检波2 59 内容提要 信号处理 检波信号处理 正交解调信号处理 B信号处理信号处理 C信号处理信号处理 多普勒信号处理信号处理 声音处理后端信号处理 60 从下图很清楚看出 正交解调分为三部分 1 频移 2 去镜像以及抗混叠滤波 3 抽取在实际的逻辑实现过程中 2和3一起完成 正交解调 组成 61 乘以正余弦因子 正交解调 理论推导 低通滤波得到基带的幅度信息以及相位信息 求模 得出相位 用于估算速度 频谱 62 信号处理 正交解调处理过程频谱图 63 第一步 频移 将信号搬迁到零频附近 第二步 低通滤波 过滤掉二倍频成分 同时降采样之前必须实现低通滤波 防止频谱混叠 为什么要降采样 假设对于深度为5cm以及10cm 采样率为40M 有用样点数分别为 0 05 2 1540 40e6 26000 1 2 1540 40e6 5200而最后显示的纵向点数是一定的 大约为 480点 即把2600个样点或5200个样点变换到最终的480个点 因此对于5cm以及10cm 降采样因子分别为 2600 480约为5 5200 480约为10 正交解调 步骤 64 上图是B信号回波的频谱示意图 横轴表示时间即探测深度 纵轴表示回波频谱成份 颜色深的红色表示能量集中信号随深度变化趋势 从图中可以看出 高频成分衰减的快 低频成份衰减的慢 因此 1 回波的中心频率随着深度增加而向低频漂移 2 回波的带宽随深度增加而减小 1 2 是设计正交解调的指导因素 正交解调 回波特点 65 正交解调 频移 66 正交解调 频移 频移对应着时域的调制 即乘以COS SIN旋转因子 对于B来说 频移的归一化角频率随时间而减小 即产生随深度增加而归一化角频率减小的正余弦因子 对于C PW来说 频移值不变 这是由超声不同模式的成像特点决定的 67 正交解调 正余弦因子的产生 几款产品COS SIN旋转因子产生的方法 68 正交解调 不同模式的正余弦因子示意 69 正交解调 分数倍采样率变换 掌握采样率变换的理论以及实现技术升采样前后频谱关系以及镜像抑制滤波 降采样前后频谱关系以及抗混叠滤波 滤波器的多相分解 6个NOBLE恒等式 利用NOBLE恒等式来简化变采样滤波器的设计 DC6 M5 M7的一个重要而必需的步骤 从理论到逻辑实现都比较枯燥 这里简单带过 70 正交解调 分数倍采样率变换 插值 增加抽样率以增加数据的过程 又称为上采样upsample 抽取 降低抽样率以减少数据的过程 又称为下采样downsample 71 分数倍采样率变换 插值 将x n 的抽样率由fs增加L倍 变为Lfs 简单的办法是将x n 每两个点之间补L 1个零 设补零后的信号为v n 则 上式说明 补 L 1 个0后 信号的带宽压缩为原来的1 L 同时产生了L 1个映象 通过一低通滤波器H jw 消除多余的L 1个映象的同时 对补充的零值进行平滑 过程如下 72 分数倍采样率变换 插值 73 分数倍采样率变换 抽取 将x n 的抽样率由fs降低M倍 变为fs M 最简单的办法是将x n 每M个点抽取1个 依次组成一个新的序列y n 则 如果原信号成分包含的角频率超过了pi M 则直接做1 M的抽取后必然会发生频谱混叠 为了防止抽取后的混叠 需要在抽取前对x n 滤波 滤波器的截止角频率为pi M 过程如下 74 分数倍采样率变换 抽取 75 待续 超声中的信号处理 二 超声中的信号处理 三 76 AGLOBALPARTNERYOUCANRELYON END 77