CN100410925C - 超声信号的数字信号处理方法 - Google Patents
超声信号的数字信号处理方法 Download PDFInfo
- Publication number
- CN100410925C CN100410925C CNB2004100114032A CN200410011403A CN100410925C CN 100410925 C CN100410925 C CN 100410925C CN B2004100114032 A CNB2004100114032 A CN B2004100114032A CN 200410011403 A CN200410011403 A CN 200410011403A CN 100410925 C CN100410925 C CN 100410925C
- Authority
- CN
- China
- Prior art keywords
- signal
- ultrasonic
- frequency
- time
- signals
- 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.)
- Expired - Fee Related
Links
Images
Abstract
本发明涉及利用数字信号处理超声信号的方法。将数字超声信号通过提升小波包变换滤波、希尔伯特(Hilbert)变换、正交解调等处理,在合成超声成像的同时进行联合时频分析和特征提取。本发明提出的集成式的数字信号处理构架,可以实时地分析射频信号、解析信号以及正交解调信号。特别地,可以使用联合时频分析技术提取出较传统的快速傅立叶变换更加丰富的有关被测材料完整性与几何特征的重要信息。作为虚拟仪器开发的软件平台,Labwindows/CVI被用于构筑了基于计算机的上层测试与控制组件。这一系统可以实现基于小波的动态滤波、降噪、下采样、联合时频分析,以及专业的测试与控制。本发明可以应用于材料的无损检测与评估、声纳、距离与速度测量以及医学应用等领域。
Description
技术领域:
本发明属于超声信号处理技术,涉及利用数字信号处理(DSP)超声信号的方法。
背景技术:
图1为超声信号处理的传统模型。由于AD转换器的采样速度较低,所以在AD转换器前端必须加入抗混叠滤波器。通常这一抗混叠滤波器的设计非常烦琐,它的电路在印刷电路板上占据了相当大的面积,而且设计实现后无法改变,无法适应不同类型的超声应用。特征提取与噪声衰减方面采用傅立叶变换。
图2为目前人们采用的一种改进的模型。虽然AD转换器的采样速度较高,但是为了提高信噪比并降低抗混叠滤波器设计的复杂度,使用4倍甚至更高的过采样率,而提高的采样率带来更大的数据量。为了提取有用信息并降低数据量,人们采用复杂可编程逻辑器件(FPGA)来实现数字滤波和下采样。FPGA虽具有在线可编程的特点,但是它不可以实时修改,配置并不方便,而且数字滤波器的设计也仍然无法实现自动化。特征提取与噪声衰减方面采用傅立叶变换,缺点是单纯地进行时域或频域的分析对于缺陷特征的识别往往不够充分。
发明内容:
为了解决上述背景技术中抗混叠滤波器设计实现后无法改变,无法适应不同类型超声应用的问题;可编程逻辑器件不能实时修改、配置不方便、设计不能自动化,对缺陷特征的识别不充分等问题,本发明将要提出一种集成式的对超声信号的数字信号处理(DSP)方法。
本发明超声信号的数字处理方法步骤如下:
a、首先将数字超声信号通过提升小波包变换将数字超声信号分解到二进制子频带,根据传感器超声频带,保留所需二进制子频带信号能量,屏蔽或压制其它频带信号能量,而后经小波反变换,综合出滤波后超声射频信号;
b、利用希尔伯特(Hilbert)变换将步骤a产生的超声射频信号转换成超声解析信号;
c、将步骤b产生的超声解析信号通过正交解调产生正交解调信号;
d、步骤c产生的正交解调信号再次通过提升小波包变换,提取低频段信号能量,得到压缩的正交解调信号;
e、步骤d产生的压缩的正交解调信号在进行超声成像的同时进行联合时频分析,得到含有时域和频域的时频信息;
f、将步骤b产生的超声解析信号和步骤e得到的时频信息进行特征提取,则完成超声信号的数字处理。
本发明的优点:
由于采用了提升小波包变换(WPT),使用正交与双正交小波从有用子带中提取超声信号来进行滤波.小波分解的级次和有用的子带根据采样频率、传感器的中心频率和带宽自动确定,就不再需要设计滤波器,解决了抗混叠滤波器设计实现后无法改变、无法适应不同类型超声应用的问题。小波包变换的最大优点是信号特征的提取非常容易,它在特征提取与噪声衰减方面的性能在许多领域中都胜过了传统的快速傅立叶变换。
由于本发明采用DSP构架,将超声信号进行处理,具有高度灵活性,可以实时修改,因而实时性佳,解决了采用可编程逻辑器件不能实时修改、配置不方便、设计不能自动化,对于缺陷特征的识别不充分等问题。特别地,这种架构使用联合时频分析技术,提取被测材料有关超声信号的完整性与几何特征的重要信息较传统的快速傅立叶变换更加丰富。
本发明可以应用于材料的无损检测与评估、声纳、距离与速度测量以及医学应用等领域。
附图说明:
图1是背景技术传统超声信号处理的结构示意图
图2是背景技术常用的超声信号处理的结构示意图
图3是本发明超声信号处理的结构示意图
图4是本发明超声测试系统中数字信号处理结构
具体实施方式:
实施例:如图3图4所示包括时间增益补偿电路、A/D转换器、数字信号处理电路、计算机。
时间增益补偿电路包括:可编程逻辑器件、DA转换器、可编程增益放大器。可编程逻辑器件采用CPLD可编程逻辑器件;DA转换器采用MAX5110高速DA转换器;可编程增益放大器采用AD604可编程增益放大器。时间增益补偿电路由CPLD经高速DA转换器控制可编程增益放大器来实现。
A/D转换器:可采用40Mhz/60Mhz的A/D转换器,如:MAX1183。
数字信号处理电路:采用TI的DSP开发板:TMS320C6701EVM,结构如图5所示。
计算机分析软件采用美国国家仪器公司的Labwindows/CVI虚拟仪器开发系统。
本发明中:
a.利用提升小波包变换将数字超声信号分解到二进制子频带,根据传感器超声频带,保留所需二进制子频带信号能量,屏蔽或压制其它频带信号能量,而后经小波反变换,综合出滤波后信号,即超声射频信号(1)。
b.以超声射频信号作为实部,超声射频信号经过Hilbert变换后的信号作虚部即构成超声解析信号(2),超声解析信号是复信号,只有正频率成分。
Hilbert变换就是将超声射频信号作90度相移,得到超声射频信号的正交信号。定义如下:
c.超声解析信号(2)的幅度即构成了超声射频信号(1)的复包络,即正交解调信号(3),所以将步骤b产生的超声解析信号(2)取模,即得到超声射频信号(1)的包络,即正交解调信号(3);
d.步骤c产生的正交解调信号(3)再次通过提升小波包变换,提取低频段信号能量即得到压缩的正交解调信号(4);
e.将压缩的正交解调信号(4)的信号幅度以色彩表示,每个扫描信号即构成一条色彩随幅度变换的线条,上千条扫描线即构成了彩色的超声图像。
利用Labwindows/CVI中相关的函数对压缩的正交解调信号(4)进行联合时频分析。联合时频分析将一维的时域信号转变成二维的能量对时间和频率的联合分布。可以将变换结果以3维或称之2.5维的方式表现出来(能力大小以光强或色彩值表现)。联合时频分析对应于一个信号,是信号能量在时间和频率方面的综合表达。它能准确的判定特征信号的出现时间和对应的频率及其强度。联合时频分析展现了信号的频率成分随时间变换的过程。通过时间和频率域的同时分析我们得到了更丰富的信息。
f.利用Labwindows/CVI中相关的函数提取信号的瞬时频率的方法,完成精确的目标定位,将步骤b产生的超声解析信号(2)和步步骤e得到的时频信息进行特征提取,则完成超声信号的数字处理。
Claims (1)
1. 超声信号的数字处理方法,其特征在于其步骤如下:
a、首先将数字超声信号通过提升小波包变换将数字超声信号分解到二进制子频带,根据传感器超声频带,保留所需二进制子频带信号能量,屏蔽或压制其它频带信号能量,而后经小波反变换,综合出滤波后超声射频信号;
b、利用希尔伯特(Hilbert)变换将步骤a产生的超声射频信号转换成超声解析信号;
c、将步骤b产生的超声解析信号通过正交解调产生正交解调信号;
d、步骤c产生的正交解调信号再次通过提升小波包变换,提取低频段信号能量,得到压缩的正交解调信号;
e、步骤d产生的压缩的正交解调信号在进行超声成像的同时进行联合时频分析,得到含有时域和频域的时频信息;
f、将步骤b产生的超声解析信号和步骤e得到的时频信息进行特征提取,则完成超声信号的数字处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100114032A CN100410925C (zh) | 2004-12-30 | 2004-12-30 | 超声信号的数字信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100114032A CN100410925C (zh) | 2004-12-30 | 2004-12-30 | 超声信号的数字信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1731388A CN1731388A (zh) | 2006-02-08 |
CN100410925C true CN100410925C (zh) | 2008-08-13 |
Family
ID=35963742
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004100114032A Expired - Fee Related CN100410925C (zh) | 2004-12-30 | 2004-12-30 | 超声信号的数字信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100410925C (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4909410B2 (ja) * | 2006-06-21 | 2012-04-04 | エヌエックスピー ビー ヴィ | 変調された信号を復調する方法、復調器及び受信機 |
CN100417191C (zh) * | 2006-09-01 | 2008-09-03 | 上海大学 | 一种图像组合降噪方法 |
CN101190137B (zh) * | 2006-11-27 | 2011-06-22 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声成像系统的实时数字正交解调方法与装置 |
CN101291517B (zh) * | 2007-04-17 | 2012-01-04 | 华为技术有限公司 | 一种测量处理方法及装置 |
CN101312529B (zh) * | 2007-05-24 | 2010-07-21 | 华为技术有限公司 | 生成上下采样滤波器及实现编码的方法、系统和装置 |
JP5049773B2 (ja) * | 2007-12-27 | 2012-10-17 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置、超音波画像処理プログラム |
CN102176121A (zh) * | 2011-01-18 | 2011-09-07 | 河海大学 | 数字超声经颅多普勒数字解调和信号处理方法及装置 |
CN106488358B (zh) * | 2015-09-09 | 2019-07-19 | 上海其高电子科技有限公司 | 优化声场成像定位方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1367955A (zh) * | 1999-08-17 | 2002-09-04 | 英国电讯有限公司 | 信号发生器和解码器 |
WO2004048903A1 (en) * | 2002-11-25 | 2004-06-10 | Instromet Ultrasonics B.V. | Ultrasonic signal processing method and applications thereof |
-
2004
- 2004-12-30 CN CNB2004100114032A patent/CN100410925C/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1367955A (zh) * | 1999-08-17 | 2002-09-04 | 英国电讯有限公司 | 信号发生器和解码器 |
WO2004048903A1 (en) * | 2002-11-25 | 2004-06-10 | Instromet Ultrasonics B.V. | Ultrasonic signal processing method and applications thereof |
Also Published As
Publication number | Publication date |
---|---|
CN1731388A (zh) | 2006-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100410925C (zh) | 超声信号的数字信号处理方法 | |
Liu et al. | Automatic mode extraction of ultrasonic guided waves using synchrosqueezed wavelet transform | |
CN110559014B (zh) | 基于探头补偿的分数阶傅里叶变换回波成像方法与系统 | |
CN106645949A (zh) | 一种基于低频检波的外差扫频式频谱分析仪 | |
CN109884192A (zh) | 用于焊缝特征导波缺陷回波信号特性提取的稀疏表示方法 | |
Xu et al. | Advanced methods for time-of-flight estimation with application to Lamb wave structural health monitoring | |
Awwad et al. | Communication network for ultrasonic acoustic water leakage detectors | |
Avanesians et al. | Wave separation: Application for arrival time detection in ultrasonic signals | |
CN104698456A (zh) | 基于预调制aic的太赫兹雷达成像方法 | |
Yang et al. | Parametric study on interply tracking in multilayer composites by analytic-signal technology | |
Gauthier et al. | Lightweight and amplitude-free ultrasonic imaging using single-bit digitization and instantaneous phase coherence | |
CN207352066U (zh) | 一种信号检测装置 | |
CN105611018B (zh) | 一种mipi lp信号测试系统及方法 | |
CN107248869B (zh) | 一种基于吕分布的多分量线性调频信号去噪方法 | |
CN101109731B (zh) | 超声信号表示的Gabor变换自适应窗宽选择方法 | |
Shoupeng et al. | Quadrature demodulation based circuit implementation of pulse stream for ultrasonic signal FRI sparse sampling | |
CN102539542A (zh) | 复小波变换法提取陶瓷基芳纶纤维增强复合板中单模式超声导波信息 | |
Zhou et al. | Anefficient quadrature demodulator for medical ultrasound imaging | |
Gammell | Analogue implementation of analytic signal processing for pulse-echo systems | |
CN206096037U (zh) | 基于时间反转和兰姆波的金属管材缺陷的超声检测装置 | |
CN109799284A (zh) | 一种超声回波信号的多次谐波自适应分离方法 | |
CN107843653B (zh) | 一种双频超声换能器及高次谐波的自发自收式测量方法 | |
CN107101715A (zh) | 一种基于交叉小波的变压器振动信号幅频特征量提取方法 | |
CN107888301A (zh) | 一种实车调频收音干扰的测试系统和方法 | |
CN110110607A (zh) | 基于谐波小波包提取狭窄带宽信号的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080813 Termination date: 20100201 |