CN105982694B - 抑制超声噪声的信号处理方法 - Google Patents
抑制超声噪声的信号处理方法 Download PDFInfo
- Publication number
- CN105982694B CN105982694B CN201510042287.9A CN201510042287A CN105982694B CN 105982694 B CN105982694 B CN 105982694B CN 201510042287 A CN201510042287 A CN 201510042287A CN 105982694 B CN105982694 B CN 105982694B
- Authority
- CN
- China
- Prior art keywords
- dynamic range
- signal
- echo
- noise
- ultrasonic
- 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.)
- Active
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开一种抑制超声噪声的信号处理方法,将接收到的回波信号经匹配滤波和包络检测后,根据发射频率及探查对象确定动态范围调整参数,从调整起点开始对信号随深度变化调整动态范围,最后随经对数压缩输出灰阶值。其中将探测深度按调整变化率分为若干区域,根据超声在人体衰减所引起的幅度变化及噪声特征,将各区域动态范围调整为最佳设定。若当前成像效果未能满足需求,可通过更改动态范围调整参数进行重新调整。并且对相邻两个动态范围对应的灰阶输出做了平滑处理,减弱图像在视觉上的突变。本发明的方法,调整了动态范围会使深部的S/N比及可用性得到改善,使分辨力及图像细微度得以改善,图像的总体质量得到提高。
Description
技术领域
本发明属于超声成像领域,具体涉及一种抑制超声噪声的信号处理方法。
背景技术
超声成像中,信噪比(signal-to-noise ratio,SNR)即回波信号中有用信号与干扰信号的功率谱之比,但通常功率谱难以计算,通常使用信号与噪声的方差之比来近似估计信噪比,信噪比越大,说明回波信号中的有效信息越多,图像质量越好。由于超声成像是从大量背景噪声中提取有用信息,成像信息容易受到各种噪声的干扰,成像质量总的来说并不理想。且如图1所示,由于超声回波信号在人体组织中会随着传播深度而衰减, S/N随深度的增加而变小,降低了设备的实用性。
在超声脉冲回波中,回波信息的一大特点即其幅度不一,在接收器中首先出现的是强发射信号,随后到来的是人体各界面和结构的回波信号,回波信号的幅度分布在一个很广的范围中,定义以分贝(db)数表示的最大信号幅度与最小信号幅度之比为信号的幅度动态范围。在临床上,当超声频率一定时,动态范围因为探测对象不同而略有不同。由于探测深度范围内超声传播衰减的存在,有用回波信号幅度的变化范围也会因此受到影响。为了能够获得接收回波信号中所含的所有信息,成像设备必须具有较大的输入动态范围,但由于噪声的存在,成像过程中需要实时调整动态范围来去除不必要的噪声信号。合适的动态范围既要保证病灶内部低弱回声信号的显示,又要保证病灶边界、强回声的突出。
传统系统中,通过系统预设或者用户实时调节来控制动态范围参数。由于人体的差异性,所获得的信号与噪声均具有较大的差异,预设参数不能满足所有临床需求,而用户实时调节需要一定的操作时间,往往不能快速方便地调节得到最优的显示效果,而且对用户的操作有比较高的要求。
在现有技术中,有一种自动优化动态范围的方法,实时采集一段时间内的多普勒线数据,估计多普勒线数据中所包含的平均噪声水平,根据该平均噪声水平计算噪声特征量;从多普勒线数据中提取信号特征量,并根据所述噪声特征量和信号特征量计算信号的实际动态范围,对系统动态范围参数进行调整,将其限制在与所述信号的实际动态范围相关的范围内。该方法虽然能够达到对参数进行优化的目的,但上述一段时间为一个或多个心动周期,一秒或其它时间长度,动态范围参数调整精度会受到限制。
还有一种自动优化动态范围的方法,先在已有各种参数下获取多普勒信号,去除其中的干扰信号,再结合单独方式获取的噪声信号,对增益与动态范围进行优化,动态范围设置为高出平均噪声的有效信号的实际动态范围,增益则设置为在不丢失信号的基础上对噪声进行最小化。但其噪声信号在不发射超声能量的情况下获得,不能准确估计噪声。
发明内容
本发明的目的是克服现有技术的不足,提出一种抑制超声噪声的信号处理方法,达到改善成像效果的目的。
为了实现这一目的,本发明所采取的技术方案包括以下步骤:
1)超声回波经波束合成后形成回波数字信号,然后经过匹配滤波和包络检测处理环节;
2)经对滤波后信号的分析,根据超声诊断应用需求设置不同的动态范围调整参数,其中包括:动态范围调整起点,动态范围调整率和动态范围调整深度;根据动态范围调整参数,分别设置对应的系统增益参数;
3)对数压缩模块对输入波束进行限制输入动态范围的对数压缩。
具体的,步骤1中,动态滤波器根据不同位置选择不同的中心频率,随着深度的变化,滤波器系数也随之变化;从近场到远场不同的衰减分成若干段不同的滤波系数,分时滤波;经过动态滤波后,再提出回波的包络,得到最终的图像信息。
步骤2中,根据当前图像特征调整动态范围调整参数,并根据当前调整需求与系统预先设定的增益列表,进行增益参数优化。所确定的系统接收回波信号时,当回波信号到达系统下发动态范围调整起点时,开始进行动态范围调整;根据回波信号特征及噪声特征量确定每个时间点的动态范围,随深度一级一级进行调整,相邻两个动态范围调整点之间的距离差固定,由动态范围调整率决定。
在步骤3对数压缩之前对图像进行低通平滑滤波,可减小图像的方差。为了减弱上述处理引起图像视觉上的突变,还可以步骤3之后对回波信号平滑处理。
步骤3中,对数压缩可以采用查找表方式:根据动态范围调整要求随着图像深度不断进行调整,每个断层的动态范围对应一个动态范围编码;将动态范围和输入信号作为查找表的地址输入端,该地址对应存储的内容即为所对应的动态压缩输出;不同动态范围对应不同的输出曲线。
本发明的优点是:将接收到的回波信号经匹配滤波和包络检测后,根据发射频率及探查对象确定动态范围调整参数,从调整起点开始对信号随深度变化调整动态范围,最后随经对数压缩输出灰阶值。其中将探测深度按调整变化率分为若干区域,根据超声在人体衰减所引起的幅度变化及噪声特征,将各区域动态范围调整为最佳设定。若当前成像效果未能满足需求,可通过更改动态范围调整参数进行重新调整。并且对相邻两个动态范围对应的灰阶输出做了平滑处理,减弱图像在视觉上的突变。 本发明的方法调整了动态范围会使深部的S/N比及可用性得到改善,使分辨力及图像细微度得以改善,图像的总体质量得到提高。
附图说明
图1是信号幅度随探测深度变化示意图。
图2是本发明的总体流程图。
图3是输入输出对应关系图。
图4是动态范围随深度一级一级调整示意图。
图5是本发明进行动态范围调整和对数压缩获得图像灰阶信号的示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
为了提高超声图像噪声抑制的效果,本发明采用如下方法进行处理,如图2所示:
第一步:对回波超声信号进行动态匹配滤波。其作用在于选择携带有用信息的频率分量,滤除近场的低频噪声信号和远场的高频噪声,有利于提高近场分辨率和远场信噪比。
超声回波经波束合成后形成回波数字信号,然后经过匹配滤波和包络检测处理环节。动态滤波器根据不同位置选择不同的中心频率,随着深度的变化,滤波器系数也随之变化。从近场到远场不同的衰减分成若干段不同的滤波系数,分时滤波。经过动态滤波后,提升了回波信号的信噪比,但我们所需要的组织反射信息仍作为调制信号存在于回波中。我们需要提出回波的包络,得到最终的图像信息。回波检测可通过数字正交包络检波得到信号幅度波形。也可以通过绝对值包络检波得到。
第二步:经对滤波后信号的分析,根据超声诊断应用需求设置不同的动态范围调整参数,具体包括:动态范围调整起点,动态范围调整率和动态范围调整深度。
本发明实施例中根据人体差异需求最高设置每探测深度内进行多档动态范围调整,用户可根据当前探查对象实际体型特征选择不同动态范围调整参数。
根据动态范围调整要求,为了能接收回波信号的所有信息,系统必须具有大的输入动态范围,对系统增益参数进行调整。根据用户选择的动态范围调整参数,分别设置对应的系统增益参数。
第三步:在对数压缩前对图像进行低通平滑滤波,减小图像的方差。
第四步:根据当前图像特征调整动态范围调整参数,其中包括动态范围调整起点,动态范围调整率和动态范围调整深度。并根据当前调整需求与系统预先设定的增益列表,进行增益参数优化。利用增益对信号线性放大,不会对信噪比产生影响。对于小信号而言,幅值比较低,在后面对数压缩转换时会出现问题,偏离真正的对数值,因此需根据动态范围调整参数设置相应的增益值。
所确定的系统接收回波信号时,当回波信号到达系统下发动态范围调整起点时,开始进行动态范围调整。根据回波信号特征及噪声特征量确定每个时间点的动态范围,随深度一级一级进行调整。如图4所示,P0为动态范围起始点,P0至Pn之间的距离即为调整深度,相邻两个动态范围调整点之间的距离差固定,由动态范围调整率决定。其中n最大值可以根据实际系统的性能去决定。在实际使用中,用户可根据探查对象的区别选择探测深度内动态范围更改的次数。其中动态范围调整参数依据对超声回波在人体中的衰减特性,探头发射频率以及噪声特征来决定。
为了能更清晰的显示断层图像,接收通道又必须具有小的动态范围,系统一般采用对数压缩来实现该功能。对数压缩模块对输入波束进行限制输入动态范围的对数压缩。如图5所示,本发明应用实例中,对数压缩采用查找表方式。根据动态范围调整要求随着图像深度不断进行调整,每个断层的动态范围对应一个动态范围编码。将动态范围和输入信号作为查找表的地址输入端,该地址对应存储的内容即为所对应的动态压缩输出。不同动态范围对应不同的输出曲线,如图3所示。
第五步:为了减弱上述处理引起图像视觉上的突变,还可以对进行噪声抑制动态压缩的回波信号平滑处理,此平滑处理可以参见现有技术中的相关方法。
本实例采用加权拼接方法,在图像相邻区域内设定一个平滑变化的权重因子进行加权处理,消除相邻动态范围对数压缩处理后图像的拼缝及模糊线性,提高融合质量。拼接完成之后对输出值进行选择。在超出动态范围调整深度部分的回波信号的对数值按照Pn最终调整动态范围下取当前信号的对数值,在调整深度范围内的值若超出了灰阶上限,则取上限值,其余取加权值。
Claims (6)
1.抑制超声噪声的信号处理方法,其特征是,包括以下步骤:
1)超声回波经波束合成后形成回波数字信号,然后经过匹配滤波和包络检测处理环节;
2)经对滤波后信号的分析,根据超声诊断应用需求设置不同的动态范围调整参数,其中包括:动态范围调整起点,动态范围调整率和动态范围调整深度;根据动态范围调整参数,分别设置对应的系统增益参数;
3)对数压缩模块对输入波束进行限制输入动态范围的对数压缩;
其中,步骤1中,动态滤波器根据不同位置选择不同的中心频率,随着深度的变化,滤波器系数也随之变化;从近场到远场不同的衰减分成若干段不同的滤波系数,分时滤波;经过动态滤波后,再提出回波的包络,得到最终的图像信息。
2.如权利要求1所述的抑制超声噪声的信号处理方法,其特征是,在步骤3对数压缩之前对图像进行低通平滑滤波。
3.如权利要求1所述的抑制超声噪声的信号处理方法,其特征是,步骤3之后再对回波信号平滑处理。
4.如权利要求1所述的抑制超声噪声的信号处理方法,其特征是,步骤2中,根据当前图像特征调整动态范围调整参数,并根据当前调整需求与系统预先设定的增益列表,进行增益参数优化。
5.如权利要求1所述的抑制超声噪声的信号处理方法,其特征是,步骤2中,所确定的系统接收回波信号时,当回波信号到达系统下发动态范围调整起点时,开始进行动态范围调整;根据回波信号特征及噪声特征量确定每个时间点的动态范围,随深度一级一级进行调整,相邻两个动态范围调整点之间的距离差固定,由动态范围调整率决定。
6.如权利要求1所述的抑制超声噪声的信号处理方法,其特征是,步骤3中,对数压缩采用查找表方式:根据动态范围调整要求随着图像深度不断进行调整,每个断层的动态范围对应一个动态范围编码;将动态范围和输入信号作为查找表的地址输入端,该地址对应存储的内容即为所对应的动态压缩输出;不同动态范围对应不同的输出曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510042287.9A CN105982694B (zh) | 2015-01-27 | 2015-01-27 | 抑制超声噪声的信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510042287.9A CN105982694B (zh) | 2015-01-27 | 2015-01-27 | 抑制超声噪声的信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105982694A CN105982694A (zh) | 2016-10-05 |
CN105982694B true CN105982694B (zh) | 2019-03-19 |
Family
ID=57034914
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510042287.9A Active CN105982694B (zh) | 2015-01-27 | 2015-01-27 | 抑制超声噪声的信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105982694B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10975686B2 (en) * | 2017-04-20 | 2021-04-13 | General Electric Company | Detection system including sensor and method of operating such |
CN109602444A (zh) * | 2018-12-06 | 2019-04-12 | 余姚市华耀工具科技有限公司 | 肾脏缺失验证机构 |
CN111281423A (zh) * | 2018-12-07 | 2020-06-16 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声图像优化方法和超声成像设备 |
CN111429432B (zh) * | 2020-03-24 | 2024-05-03 | 聚融医疗科技(杭州)有限公司 | 基于射频处理和模糊聚类的热消融区域监测方法及系统 |
CN113679418B (zh) * | 2021-08-24 | 2024-03-22 | 合肥工业大学 | 一种基于压缩信噪比系数加权的超声成像方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101273907A (zh) * | 2007-03-28 | 2008-10-01 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置以及处理方法 |
CN101766497A (zh) * | 2008-12-31 | 2010-07-07 | 深圳迈瑞生物医疗电子股份有限公司 | 一种声谱图像的信号处理方法和系统 |
CN101836868A (zh) * | 2009-03-16 | 2010-09-22 | 深圳迈瑞生物医疗电子股份有限公司 | 自适应的动态范围调整方法和设备 |
CN101879076A (zh) * | 2009-05-08 | 2010-11-10 | 深圳迈瑞生物医疗电子股份有限公司 | 多普勒超声成像的自动优化方法及其装置 |
CN101919714A (zh) * | 2010-08-18 | 2010-12-22 | 秦皇岛市康泰医学系统有限公司 | 嵌入式b型超声诊断设备及其信号处理方法 |
CN101972152A (zh) * | 2010-12-01 | 2011-02-16 | 深圳市蓝韵实业有限公司 | 一种b型图像增益参数自动优化方法 |
CN102133089A (zh) * | 2011-03-28 | 2011-07-27 | 李鹏 | 血流速度估计方法和利用该方法的彩色血流图生成方法 |
-
2015
- 2015-01-27 CN CN201510042287.9A patent/CN105982694B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101273907A (zh) * | 2007-03-28 | 2008-10-01 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置以及处理方法 |
CN101766497A (zh) * | 2008-12-31 | 2010-07-07 | 深圳迈瑞生物医疗电子股份有限公司 | 一种声谱图像的信号处理方法和系统 |
CN101836868A (zh) * | 2009-03-16 | 2010-09-22 | 深圳迈瑞生物医疗电子股份有限公司 | 自适应的动态范围调整方法和设备 |
CN101879076A (zh) * | 2009-05-08 | 2010-11-10 | 深圳迈瑞生物医疗电子股份有限公司 | 多普勒超声成像的自动优化方法及其装置 |
CN101919714A (zh) * | 2010-08-18 | 2010-12-22 | 秦皇岛市康泰医学系统有限公司 | 嵌入式b型超声诊断设备及其信号处理方法 |
CN101972152A (zh) * | 2010-12-01 | 2011-02-16 | 深圳市蓝韵实业有限公司 | 一种b型图像增益参数自动优化方法 |
CN102133089A (zh) * | 2011-03-28 | 2011-07-27 | 李鹏 | 血流速度估计方法和利用该方法的彩色血流图生成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105982694A (zh) | 2016-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105982694B (zh) | 抑制超声噪声的信号处理方法 | |
US6120446A (en) | Diagnostic medical ultrasonic imaging system and method with adaptive gain | |
US7787680B2 (en) | System and method for processing an image | |
JP4597491B2 (ja) | 医用診断イメージングシステム及びイメージの利得を適応的に制御するための方法 | |
CN101766497B (zh) | 一种声谱图像的信号处理方法和系统 | |
JP3432204B2 (ja) | 超音波診断装置 | |
JP4891038B2 (ja) | 映像処理システム及び方法 | |
KR101140934B1 (ko) | 초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치 | |
JP5324589B2 (ja) | 超音波画像診断方法及びシステム | |
EP2995259A1 (en) | Ultrasound optimization method and ultrasonic medical device therefor | |
US20150359507A1 (en) | Ultrasound diagnosis apparatus and ultrasound image processing method | |
US6899681B1 (en) | Automated power level for contrast agent imaging | |
JP2009066420A (ja) | 医療診断超音波イメージングの体積画像のゲイン最適化 | |
JP2006043457A (ja) | 超音波システムを制御する方法及びシステム | |
JP4996247B2 (ja) | 超音波診断装置 | |
US20110054317A1 (en) | Tracking and optimizing gain and contrast in real-time for ultrasound imaging | |
US20100069755A1 (en) | Ultrasonic diagnostic apparatus | |
US20060173312A1 (en) | Ultrasonic imaging system | |
CN104757999A (zh) | 一种基于超声基波和谐波的非线性成像方法和系统 | |
CN101879076B (zh) | 多普勒超声成像的自动优化方法及其装置 | |
CN102697521A (zh) | 一种自适应的医疗超声系统成像中的声学参数优化的方法 | |
WO2018000359A1 (zh) | 一种增强超声造影图像的方法、系统及超声造影成像设备 | |
US8313435B2 (en) | Clutter signal filtering in an ultrasound system | |
KR102606973B1 (ko) | 초음파 신호 데이터의 깊이 의존 감쇠를 보상하는 방법 및 시스템 | |
JP4808373B2 (ja) | Bモード画像のバンディング抑制に関連する応用のための方法及び装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 214028 Changjiang Road in Jiangsu province Wuxi new area, industrial park five, 51 No. 53 No. 228 block Applicant after: Wuxi CHISON medical Polytron Technologies Inc Address before: 214028 Changjiang Road in Jiangsu province Wuxi new area, industrial park five, 51 No. 53 No. 228 block Applicant before: Xiangsheng Medical Image Co., Ltd., Wuxi |
|
GR01 | Patent grant | ||
GR01 | Patent grant |