CN104523291B - 一种超声回波信号的脉冲压缩与正交解调处理方法及装置 - Google Patents

一种超声回波信号的脉冲压缩与正交解调处理方法及装置 Download PDF

Info

Publication number
CN104523291B
CN104523291B CN201410720045.6A CN201410720045A CN104523291B CN 104523291 B CN104523291 B CN 104523291B CN 201410720045 A CN201410720045 A CN 201410720045A CN 104523291 B CN104523291 B CN 104523291B
Authority
CN
China
Prior art keywords
signal
sampled
frequency spectrum
spectrum shift
pulse compression
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
Application number
CN201410720045.6A
Other languages
English (en)
Other versions
CN104523291A (zh
Inventor
王平
蒋辉
于铁柱
唐英勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CHONGQING BOEN FUKE MEDICAL EQUIPMENT Co Ltd
Original Assignee
CHONGQING BOEN FUKE MEDICAL EQUIPMENT Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by CHONGQING BOEN FUKE MEDICAL EQUIPMENT Co Ltd filed Critical CHONGQING BOEN FUKE MEDICAL EQUIPMENT Co Ltd
Priority to CN201410720045.6A priority Critical patent/CN104523291B/zh
Publication of CN104523291A publication Critical patent/CN104523291A/zh
Application granted granted Critical
Publication of CN104523291B publication Critical patent/CN104523291B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明实施例公开了一种超声回波信号的脉冲压缩与正交解调处理方法及装置,该方法包括:对接收的超声回波信号进行波束合成,获取合成信号;对所述合成信号进行频谱搬移,获取频谱搬移信号;对所述频谱搬移信号进行低通滤波,获取低频分量信号;对所述低频分量信号进行K倍降采样,获取降采样后信号,所述K为大于2的整数;对所述降采样信号进行脉冲压缩和取模运算,获取用于成像的扫描线数据;采用本发明的方法及装置,可有效减小脉冲压缩环节的运算量,进而提高整个超声回波信号处理环节的运算效率。

Description

一种超声回波信号的脉冲压缩与正交解调处理方法及装置
技术领域
本发明涉及超声成像技术领域,特别是涉及一种超声回波信号的脉冲压缩与正交解调处理方法及装置。
背景技术
在传统的超声成像中,发射脉冲一般采用短脉冲信号,但是受到医用超声安全剂量的限制,所发射的超声波的能量均较弱,经人体组织衰减反射后,所得到的超声回波信号较弱,成像效果较差。
目前,为了提高超声成像的质量,一般采用编码激励技术。所谓编码激励技术,是指在超声波的发射模式下对发射信号进行编码激励,在接收模块下,对超声回波信号进行处理,以提取有用信号并去除噪声干扰。采用编码激励技术,可有效提高超声成像的信噪比、信号穿透力和纵向分辨率等,从而有效提高超声成像的质量。在现有技术中,一般采用如下方法,对超声编码回波信号进行处理:首先对接收的编码回波信号进行波束合成,然后进行脉冲压缩、频谱搬移和低通滤波,再然后进行K倍降采样和取模,即可获得用于成像的扫描线数据。
但是,由于上述脉冲压缩环节置于K倍降采样环节之前,这将大大增加脉冲压缩环节的运算量,从而增加整个超声回波信号处理环节的运算量,大幅度降低了成像系统的实时性。
发明内容
本发明实施例中提供了一种超声回波信号的脉冲压缩与正交解调处理方法及装置,以减小脉冲压缩环节的运算量,进而减小整个超声回波信号处理环节的工作量。
为了解决上述技术问题,本发明实施例公开了如下技术方案:
一种超声回波信号的脉冲压缩与正交解调处理方法,包括:
对接收的超声回波信号进行波束合成,获取波束合成信号;
对所述波束合成信号进行频谱搬移,获取频谱搬移信号;
对所述频谱搬移信号进行低通滤波,获取低频分量信号;
对所述低频分量信号进行K倍降采样,获取降采样信号,所述K为大于2的整数;
对所述降采样信号进行脉冲压缩,获取脉冲压缩信号;
对所述脉冲压缩信号进行取模运算,获取用于成像的扫描线数据。
优选的,所述频谱搬移信号包括第一频谱搬移信号和第二频谱搬移信号,所述对波束合成信号进行频谱搬移,获取频谱搬移信号,包括:
将所述波束合成信号r[n]分别与cos[2πf0n]和-sin[2πf0n]相乘,以获取所述第一频谱搬移信号和第二频谱搬移信号;所述f0为所述波束合成信号的中心频率,所述n为整数。
优选的,所述低频分量信号包括第一低频分量信号和第二低频分量信号,所述对频谱搬移信号进行低通滤波,获取低频分量信号,包括:
对所述第一频谱搬移信号进行低通滤波,以获取所述第一低频分量信号I[n];
对所述第二频谱搬移信号进行低通滤波,以获取所述第二低频分量信号Q[n]。
优选的,所述降采样信号包括第一降采样信号和第二降采样信号,所述对低频分量信号进行K倍降采样,获取降采样信号,包括:
对所述I[n]进行K倍降采样,以获取所述第一降采样信号IK[n];
对所述Q[n]进行K倍降采样,以获取所述第二降采样信号QK[n]。
优选的,所述对降采样信号进行脉冲压缩,获取脉冲压缩信号,包括:
利用公式I′Kc[n]+jQ′Kc[n]=(IK[n]+jQK[n])*(cBK[n])对所述降采样信号的同相分量IK[n]和正交分量QK[n]进行脉冲压缩,获取脉冲压缩信号的同相分量I′Kc[n]和正交分量Q′Kc[n];所述cBK[n]为脉冲压缩滤波函数,所述所述IBK[n]为脉冲压缩滤波器的同相分量,所述QBK[n]为脉冲压缩滤波器的正交分量,所述LPF表示低通滤波,所述↓K表示K倍降采样,所述c[n]表示脉冲压缩滤波函数。
一种超声回波信号的脉冲压缩与正交解调处理装置,包括:
波束合成模块,用于对接收的超声回波信号进行波束合成,获取波束合成信号;
频谱搬移模块,用于对所述波束合成信号进行频谱搬移,获取频谱搬移信号;
低通滤波模块,用于对所述频谱搬移信号进行低通滤波,获取低频分量信号;
降采样模块,用于对所述低频分量信号进行K倍降采样,获取降采样信号,所述K为大于2的整数;
脉冲压缩模块,用于对所述降采样信号进行脉冲压缩,获取脉冲压缩信号;
取模运算模块,用于对所述脉冲压缩信号进行取模运算,获取成像的扫描线数据。
优选的,所述频谱搬移信号包括第一频谱搬移信号和第二频谱搬移信号,所述频谱搬移模块包括:
第一频谱搬移单元,用于将所述波束合成信号r[n]与cos[2πf0n]相乘,以获取所述第一频谱搬移信号,所述f0为所述波束合成信号的中心频率,所述n为整数;
第二频谱搬移单元,用于将所述波束合成信号r[n]与-sin[2πf0n]相乘,以获取所述第二频谱搬移信号。
优选的,所述低频分量信号包括第一低频分量信号和第二低频分量信号,所述低通滤波模块,包括:
第一低通滤波单元,用于对所述第一频谱搬移信号进行低通滤波,以获取所述第一低频分量信号I[n];
第二低通滤波单元,用于对所述第二频谱搬移信号进行低通滤波,以获取所述第二低频分量信号Q[n]。
优选的,所述降采样信号包括第一降采样信号和第二降采样信号,所述降采样模块,包括:
第一降采样单元,用于对所述I[n]进行K倍降采样,以获取所述第一降采样信号IK[n];
第二降采样单元,用于对所述Q[n]进行K倍降采样,以获取所述第二降采样信号QK[n]。
优选的,所述脉冲压缩模块,包括:
脉冲压缩单元,用于利用公式I′Kc[n]+jQ′Kc[n]=(IK[n]+jQK[n])*(cBK[n])对所述降采样信号的同相分量IK[n]和正交分量QK[n]进行脉冲压缩,获得压缩信号的同相分量I′Kc[n]和正交分量Q′Kc[n];所述cBK[n]为复基带压缩滤波函数,所述所述IBK[n]为压缩滤波器的同相分量,所述QBK[n]为压缩滤波器的正交分量,所述LPF表示低通滤波,所述↓K表示K倍降采样,所述c[n]表示传统脉冲压缩滤波函数。
本发明的有益效果包括:假设采样点个数为N,压缩滤波系数为M,采用现有技术中的超声回波信号的脉冲压缩和正交解调处理方法,脉冲压缩环节的运算量为N·M次乘法运算,而由于本发明的K倍降采样环节置于脉冲压缩环节之前,因此采样点个数变为N/K,压缩滤波系数变为M/K,脉冲压缩环节的运算量变为4N·M/K2次乘法运算(4表示需进行4次卷积运算),K为大于2的整数,因此可见,采用本发明的方法及装置,可减小脉冲压缩环节的运算量,进而减小了整个超声回波信号处理环节的运算量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的超声回波信号的脉冲压缩和正交解调处理方法的一流程示意图;
图2为本发明实施例提供的超声回波信号的脉冲压缩和正交解调处理方法的另一流程示意图;
图3为传统的超声回波信号的脉冲压缩和正交解调处理方法的一流程示意图;
图4为本发明实施例提供的超声回波信号的脉冲压缩与正交解调处理装置的一结构示意图;
图5为本发明实施例提供的超声回波信号的脉冲压缩与正交解调处理装置的另一结构示意图;
图6为本发明实施例提供的超声回波信号的脉冲压缩与正交解调处理装置的又一结构示意图;
图7为本发明实施例提供的超声回波信号的脉冲压缩与正交解调处理装置的另一结构示意图;
图8为本发明实施例提供的超声回波信号的脉冲压缩与正交解调处理装置的又一结构示意图。
具体实施方式
本发明实施例提供一种超声回波信号的脉冲压缩与正交解调处理方法及装置,以减小脉冲压缩环节的运算量,从而减小整个超声回波信号处理环节的运算量。
为了使本技术领域的人员更好地理解本发明中的技术方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
首先对本发明实施例的超声回波信号的脉冲压缩和正交解调处理方法进行说明,如图1所示,至少包括以下步骤:
S11:对接收的超声回波信号进行波束合成,获取波束合成信号;
S12:对所述波束合成信号进行频谱搬移,获取频谱搬移信号;
S13:对所述频谱搬移信号进行低通滤波,获取低频分量信号;
其中,所述低频分量信号的特性与超声波束合成信号的特性一致;
S14:对所述低频分量信号进行K倍降采样,获取降采样信号,所述K为大于2的整数;
S15:对所述降采样信号进行脉冲压缩,获取脉冲压缩信号;
S16:对所述脉冲压缩信号进行取模运算,获取成像的超声扫描线数据。
在本发明实施例中,假设采样点个数为N,压缩滤波系数为M,采用现有技术中的超声回波信号的脉冲压缩和正交解调处理方法,脉冲压缩环节的运算量为N·M次乘法运算,而由于本发明的K倍降采样环节置于脉冲压缩环节之前,因此采样点个数变为N/K,压缩滤波系数变为M/K,脉冲压缩环节的运算量变为4N·M/K2次乘法运算(4表示需进行4次卷积运算),K为大于2的整数,因此可见,采用本发明的处理方法,可减小脉冲压缩环节的运算量,进而减小了整个超声回波信号处理环节的运算量。
在本发明的另一可行实施例中,对超声回波信号的脉冲压缩与正交解调整个处理过程,如图2所示:
其中,所述频谱搬移信号包括第一频谱搬移信号和第二频谱搬移信号,上述所有实施例中的步骤S12可具体为:
将所述波束合成信号r[n]分别与cos[2πf0n]和-sin[2πf0n]相乘,以获取第一频谱搬移信号和第二频谱搬移信号;所述f0表示所述合成信号的中心频率,所述n为整数。
其中,所述低频分量信号包括第一低频分量信号和第二低频分量信号,上述所有实施例中的步骤S13可具体为:
对所述第一频谱搬移信号进行低通滤波,以获取第一低频分量信号I[n];
对所述第二频谱搬移信号进行低通滤波,以获取第二低频分量信号Q[n];
其中,所述降采样信号包括第一降采样信号和第二降采样信号,上述所有实施例中的步骤S14可具体为:
对所述I[n]进行K倍降采样,以获取第一降采样信号IK[n];
对所述Q[n]进行K倍降采样,以获取第二降采样信号QK[n]。
在本发明实施例中,采样频率通常大于4f0
在本发明实施例中,上述所有实施例中的步骤S15可具体包括:
利用公式I′Kc[n]+jQ′Kc[n]=(IK[n]+jQK[n])*(cBK[n])对所述降采样信号的同相分量IK[n]和正交分量QK[n]进行脉冲压缩,获得压缩信号的同相分量I′Kc[n]和正交分量Q′Kc[n];所述cBK[n]为复基带压缩滤波函数,所述所述IBK[n]为压缩滤波器的同相分量,所述QBK[n]为压缩滤波器的正交分量,所述LPF表示低通滤波,所述↓K表示K倍降采样,所述c[n]表示传统脉冲压缩滤波函数。
由于可见,在本发明实施例中,脉冲压缩需要4次卷积运算,而输入信号(IK[n]和QK[n])以及压缩滤波系数cBK[n]已进行了K倍降采样,此时脉冲压缩的运算总量降为原脉冲压缩运算总量的K2/4倍。由于在实际应用中,采样因子K总大于2,采用本发明的方法,脉冲压缩的运算量将小于现有技术中脉冲压缩的运算量,且成平方律变化。但在实际应用中,K不能一味增大,否则成像效果也会逐渐下降。根据采样定理,降采样率最大为fs/2f0=10倍。
在现有技术中,对超声回波信号的整个处理过程,可如图3所示,可见,r[n]代表超声回波信号的数字量,c[n]为传统脉冲压缩滤波函数。
根据离散卷积定义可得
将公式2代入公式1中可得:
I K c + jQ K c = L P F [ ( r [ n ] · e - j 2 πf 0 n ) * ( c [ n ] · e - j 2 πf 0 n ) ] ↓ K ≈ L P F ( r [ n ] · e - j 2 πf 0 n ) ↓ K * L P F ( c [ n ] · e - j 2 πf 0 n ) ↓ K = ( I K + jQ K ) * ( I B K + jQ B K ) = I K c ′ + jQ K c ′
由上可见,本发明的处理方法与传统的处理方法相比,误差只产生在低通滤波环节与卷积环节,其它环节都是等效的。因此,本发明的方法与传统方法所生在的扫描线数据相当,即成像效果相当,但本发明的运算量远低于传统方法。
需要说明的是,本发明的方法不仅完整地保留编码信号各项特性,克服了传统超声成像在图像信噪比、纵向分辨率等方面缺陷,同时极大地降低了编码激励技术的运算复杂度,并且还能通过改变降采样因子K对运算量动态调整。因此,本发明方法可以有效降低滤波器阶数,降低超声硬件系统的成本,对于提高超声系统的实时性具有重要的参考价值。
通过以上的方法实施例的描述,所属领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:只读存储器(ROM)、随机存取存储器(RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
与上述方法相对应的,本发明还提出了一种超声回波信号的脉冲压缩与正交解调处理装置,如图4所示,至少包括:
波束合成模块41,用于对接收的超声回波信号进行波束合成,获取波束合成信号;
频谱搬移模块42,用于对所述波束合成信号进行频谱搬移,获取频谱搬移信号;
低通滤波模块43,用于对所述频谱搬移信号进行低通滤波,获取低频分量信号;
降采样模块44,用于对所述低频分量信号进行K倍降采样,获取降采样信号;
脉冲压缩模块45,用于对所述降采样信号进行脉冲压缩,获取脉冲压缩信号;
取模运算模块46,用于对所述脉冲压缩信号进行取模运算,获取成像的扫描线数据。
假设采样点个数为N,压缩滤波系数为M,采用现有技术中的超声回波信号的脉冲压缩和正交解调处理方法,脉冲压缩环节的运算量为N·M次乘法运算,而由于本发明的K倍降采样环节置于脉冲压缩环节之前,因此采样点个数变为N/K,压缩滤波系数变为M/K,脉冲压缩环节的运算量变为4N·M/K2次乘法运算(4表示需进行4次卷积运算),K为大于2的整数,因此可见,采用本发明的装置,可减小脉冲压缩环节的运算量,进而减小了整个超声回波信号处理环节的运算量。
在本发明的其它可行实施例中,上述所有实施例中的频谱搬移信号包括第一频谱搬移信号和第二频谱搬移信号,如图5所示,频谱搬移模块42包括:
第一频谱搬移单元51,用于将所述合成信号r[n]与cos[2πf0n]相乘,以获取第一频谱搬移信号,所述f0表示所述合成信号的中心频率,所述n为整数;
第二频谱搬移单元52,用于将所述合成信号r[n]与-sin[2πf0n]相乘,以获取第二频谱搬移信号。
在本发明的其它可行实施例中,所述低频分量信号包括第一低频分量信号和第二低频分量信号,如图6所示,低通滤波模块43,包括:
第一低通滤波单元61,用于对所述第一频谱搬移信号进行低通滤波,以获取第一低频分量信号I[n];
第二低通滤波单元62,用于对所述第二频谱搬移信号进行低通滤波,以获取第二低频分量信号Q[n]。
在本发明的其它可行实施例中,所述降采样信号包括第一降采样信号和第二降采样信号,如图7所示,降采样模块44,包括:
第一降采样单元71,用于对所述I[n]进行K倍降采样,以获取第一降采样信号IK[n];
第二降采样单元72,用于对所述Q[n]进行K倍降采样,以获取第二降采样信号QK[n]。
在本发明的其它可行实施例中,如图8所示,脉冲压缩模块45,包括:
脉冲压缩单元81,用于利用公式I′Kc[n]+jQ′Kc[n]=(IK[n]+jQK[n])*(cBK[n])对所述降采样信号的同相分量IK[n]和正交分量QK[n]进行脉冲压缩,获得压缩信号的同相分量I′Kc[n]和正交分量Q′Kc[n];所述cBK[n]为脉冲压缩滤波函数,所述所述IBK[n]为压缩滤波器的同相分量,所述QBK[n]为压缩滤波器的正交分量,所述LPF表示低通滤波,所述↓K表示K倍降采样,所述c[n]表示传统脉冲压缩滤波函数。
需要说明的是,在本文中,诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅是本发明的具体实施方式,使本领域技术人员能够理解或实现本发明。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (10)

1.一种超声回波信号的脉冲压缩和正交解调处理方法,其特征在于,包括:
对接收的超声回波信号进行波束合成,获取波束合成信号;
对所述波束合成信号进行频谱搬移,获取频谱搬移信号;
对所述频谱搬移后的信号进行低通滤波,获取有效的低频分量信号;
对所述低频分量信号进行K倍降采样,获取降采样信号,所述K为大于2的整数;
对所述降采样信号进行脉冲压缩,获取脉冲压缩信号;
对所述脉冲压缩信号进行取模运算,获取用于成像的扫描线数据。
2.根据权利要求1所述的方法,其特征在于,所述频谱搬移信号包括第一频谱搬移信号和第二频谱搬移信号,所述对波束合成信号进行频谱搬移,获取频谱搬移信号,包括:
将所述波束合成信号r[n]分别与cos[2πf0n]和-sin[2πf0n]相乘,以获取所述第一频谱搬移信号和第二频谱搬移信号;所述f0为所述波束合成信号的中心频率,所述n为整数。
3.根据权利要求2所述的方法,其特征在于,所述低频分量信号包括第一低频分量信号和第二低频分量信号,所述对频谱搬移后的信号进行低通滤波,获取有效的低频分量信号,包括:
对所述第一频谱搬移信号进行低通滤波,以获取所述第一低频分量信号I[n];
对所述第二频谱搬移信号进行低通滤波,以获取所述第二低频分量信号Q[n]。
4.根据权利要求3所述的方法,其特征在于,所述降采样信号包括第一降采样信号和第二降采样信号,所述对低频分量信号进行K倍降采样,获取降采样信号,包括:
对所述I[n]进行K倍降采样,以获取所述第一降采样信号IK[n];
对所述Q[n]进行K倍降采样,以获取所述第二降采样信号QK[n]。
5.根据权利要求4所述的方法,其特征在于,所述对降采样信号进行脉冲压缩,获取压缩信号,包括:
利用公式I′Kc[n]+jQ′Kc[n]=(IK[n]+jQK[n])*(cBK[n])对所述降采样信号的同相分量IK[n]和正交分量QK[n]进行脉冲压缩,获取压缩信号的同相分量I′Kc[n]和正交分量Q′Kc[n],“*”表示卷积运算;所述cBK[n]为脉冲压缩滤波函数,所述所述IBK[n]为脉冲压缩滤波器的同相分量,所述QBK[n]为压缩滤波器的正交分量,所述LPF表示低通滤波,所述↓K表示K倍降采样,所述c[n]表示传统脉冲压缩滤波函数。
6.一种超声回波信号的脉冲压缩与正交解调处理装置,其特征在于,包括:
波束合成模块,用于对接收的超声回波信号进行波束合成,获取波束合成信号;
频谱搬移模块,用于对所述波束合成信号进行频谱搬移,获取频谱搬移信号;
低通滤波模块,用于对所述频谱搬移信号进行低通滤波,获取低频分量信号;
降采样模块,用于对所述低频分量信号进行K倍降采样,获取降采样信号,所述K为大于2的整数;
脉冲压缩模块,用于对所述降采样信号进行脉冲压缩,获取脉冲压缩信号;
取模运算模块,用于对所述脉冲压缩信号进行取模运算,获取成像的扫描线数据。
7.根据权利要求6所述的装置,其特征在于,所述频谱搬移信号包括第一频谱搬移信号和第二频谱搬移信号,所述频谱搬移模块包括:
第一频谱搬移单元,用于将所述波束合成信号r[n]与cos[2πf0n]相乘,以获取所述第一频谱搬移信号,所述f0为所述波束合成信号的中心频率,所述n为整数;
第二频谱搬移单元,用于将所述波束合成信号r[n]与-sin[2πf0n]相乘,以获取所述第二频谱搬移信号。
8.根据权利要求7所述的装置,其特征在于,所述低频分量信号包括第一低频分量信号和第二低频分量信号,所述低通滤波模块,包括:
第一低通滤波单元,用于对所述第一频谱搬移信号进行低通滤波,以获取所述第一低频分量信号I[n];
第二低通滤波单元,用于对所述第二频谱搬移信号进行低通滤波,以获取所述第二低频分量信号Q[n]。
9.根据权利要求8所述的装置,其特征在于,所述降采样信号包括第一降采样信号和第二降采样信号,所述降采样模块,包括:
第一降采样单元,用于对所述I[n]进行K倍降采样,以获取所述第一降采样信号IK[n];
第二降采样单元,用于对所述Q[n]进行K倍降采样,以获取所述第二降采样信号QK[n]。
10.根据权利要求9所述的装置,其特征在于,所述脉冲压缩模块,包括:
脉冲压缩单元,用于利用公式I′Kc[n]+jQ′Kc[n]=(IK[n]+jQK[n])*(cBK[n])对所述降采样信号的同相分量IK[n]和正交分量QK[n]进行脉冲压缩,获得压缩信号的同相分量I′Kc[n]和正交分量Q′Kc[n];所述cBK[n]为脉冲压缩滤波函数,所述所述IBK[n]为压缩滤波器的同相分量,所述QBK[n]为压缩滤波器的正交分量,所述LPF表示低通滤波,所述↓K表示K倍降采样,所述c[n]表示传统脉冲压缩滤波函数。
CN201410720045.6A 2014-12-02 2014-12-02 一种超声回波信号的脉冲压缩与正交解调处理方法及装置 Expired - Fee Related CN104523291B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410720045.6A CN104523291B (zh) 2014-12-02 2014-12-02 一种超声回波信号的脉冲压缩与正交解调处理方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410720045.6A CN104523291B (zh) 2014-12-02 2014-12-02 一种超声回波信号的脉冲压缩与正交解调处理方法及装置

Publications (2)

Publication Number Publication Date
CN104523291A CN104523291A (zh) 2015-04-22
CN104523291B true CN104523291B (zh) 2016-09-14

Family

ID=52839062

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410720045.6A Expired - Fee Related CN104523291B (zh) 2014-12-02 2014-12-02 一种超声回波信号的脉冲压缩与正交解调处理方法及装置

Country Status (1)

Country Link
CN (1) CN104523291B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113397591B (zh) * 2021-06-16 2022-11-29 复旦大学 微纳马达实时成像与追踪方法、装置及微纳马达控制方法
CN113892973B (zh) * 2021-12-13 2022-03-01 深圳迈瑞动物医疗科技有限公司 超声成像设备及其超声回波数据降采样因子的确定方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201333056Y (zh) * 2008-12-05 2009-10-28 中国医学科学院生物医学工程研究所 高频超声编码激励与接收系统
KR101999078B1 (ko) * 2010-06-09 2019-07-10 리전츠 오브 더 유니버스티 오브 미네소타 초음파 치료의 전달을 제어하기 위한 이중 모드 초음파 트랜스듀서(dmut) 시스템 및 방법
CN101999909A (zh) * 2010-11-26 2011-04-06 温州医学院眼视光研究院 全数字眼科超声生物显微镜
JP6074299B2 (ja) * 2013-03-22 2017-02-01 富士フイルム株式会社 超音波診断装置、超音波診断装置の信号処理方法およびプログラム

Also Published As

Publication number Publication date
CN104523291A (zh) 2015-04-22

Similar Documents

Publication Publication Date Title
US10371804B2 (en) Ultrasound signal processing circuitry and related apparatus and methods
CN107461611B (zh) 基于小波和emd重构相结合的泄漏检测方法及泄漏检测装置
CN103178853B (zh) 基于压缩感知的稀疏信号欠采样方法
CN204495996U (zh) 宽带声纳接收器
CN101496728B (zh) 超声频率复合成像方法与装置
CN105997148B (zh) 脉冲多普勒超高谱分辨率成像处理方法及处理系统
CN107192553A (zh) 基于盲源分离的齿轮箱复合故障诊断方法
CN104523291B (zh) 一种超声回波信号的脉冲压缩与正交解调处理方法及装置
CN101161205B (zh) 多普勒血流声音信号抗混叠的方法及装置
EP1607044B1 (en) Ultrasonic imaging device
US6704438B1 (en) Apparatus and method for improving the signal to noise ratio on ultrasound images using coded waveforms
CN103126725A (zh) 一种超声成像的方法和装置
Hosman et al. Design and characterization of an MFSK-based transmitter/receiver for ultrasonic communication through metallic structures
Zhang et al. Fault diagnosis for gearbox based on EMD-MOMEDA
CN105212964B (zh) 基于rf数据超声成像处理方法及系统
CN107561586B (zh) 一种气泡压制的方法和装置
Guo et al. Underwater target detection and localization with feature map and CNN-based classification
CN105919625A (zh) 脉冲多普勒壁滤波处理方法及处理系统
CN108266432A (zh) 一种泵源液压系统压力脉动自适应滤波方法及系统
CN101690671A (zh) 一种多普勒超声成像系统动态解调装置
CN103728632B (zh) 卫星数据捕获装置和捕获方法
US11488574B2 (en) Method and system for implementing a modal processor
CN104224232A (zh) 超声谐波成像方法及装置
US10825443B2 (en) Method and system for implementing a modal processor
He et al. Signal-noise separation of sensor signal based on variational mode decomposition

Legal Events

Date Code Title Description
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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160914

Termination date: 20201202

CF01 Termination of patent right due to non-payment of annual fee