CN112263278A - 基于随机孔径的心脏超声杂波抑制方法及系统 - Google Patents

基于随机孔径的心脏超声杂波抑制方法及系统 Download PDF

Info

Publication number
CN112263278A
CN112263278A CN202011291164.6A CN202011291164A CN112263278A CN 112263278 A CN112263278 A CN 112263278A CN 202011291164 A CN202011291164 A CN 202011291164A CN 112263278 A CN112263278 A CN 112263278A
Authority
CN
China
Prior art keywords
imaging position
group
weighted
attribute
random
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.)
Granted
Application number
CN202011291164.6A
Other languages
English (en)
Other versions
CN112263278B (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.)
Vinno Technology Suzhou Co Ltd
Original Assignee
Vinno Technology Suzhou 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 Vinno Technology Suzhou Co Ltd filed Critical Vinno Technology Suzhou Co Ltd
Priority to CN202011291164.6A priority Critical patent/CN112263278B/zh
Publication of CN112263278A publication Critical patent/CN112263278A/zh
Application granted granted Critical
Publication of CN112263278B publication Critical patent/CN112263278B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0883Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
    • 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/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02Preprocessing
    • G06F2218/04Denoising

Abstract

本发明提供了一种基于随机孔径的心脏超声杂波抑制方法及系统,所述方法包括:获取N个阵元经过延时处理后分别在每一成像位置对齐形成的N个通道信号;随机生成N维的第一组随机矢量组和第二组随机矢量组;以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值;根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数;根据判断参数的大小确认每一成像位置的属性;依据每一成像位置的属性确定每一成像位置的加权叠加系数;以加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出。本发明抑制成像过程中的杂波噪声,使图像更干净清晰,提升成像效果。

Description

基于随机孔径的心脏超声杂波抑制方法及系统
技术领域
本发明涉及医疗设备成像领域,尤其涉及一种基于随机孔径的心脏超声杂波抑制方法及系统。
背景技术
超声成像技术具有快速、实时等特点,一直是心脏医学领域最普遍采用的检测手段。由于心脏运动剧烈,腔室结构复杂造成反射多,心血管各部分大小、厚度极其不均匀,心脏瓣膜运动状况受到疾病影响不规律等原因使得心脏在进行超声波图像采集过程中,充斥着大量的杂波噪声,这些噪声在无回声的心腔内尤为明显,严重干扰了临床诊断。
由于心超成像帧频较高,复杂度较高的信号、图像处理方法难以适用,目前解决噪声问题的现有技术均存在技术缺陷。
现有技术中,一种心超降噪方法为:提高persistence(平滑)设置,这种方法有一定的降噪效果,但是这种方法也会降低图像帧频,使心脏图像序列发生卡顿和模糊,丢失时间细节信息。另一种降噪方法为:通过降低图像增益来抑制噪声,该种方式会使得心脏结构变暗,图像对比度下降,影响成像质量。
发明内容
为解决上述技术问题,本发明的目的在于提供一种基于随机孔径的心脏超声杂波抑制方法及系统。
为了实现上述发明目的之一,本发明一实施方式提供一种基于随机孔径的心脏超声杂波抑制方法,所述方法包括:获取N个阵元经过延时处理后分别在每一成像位置对齐形成的N个通道信号;
随机生成N维的第一组随机矢量组和第二组随机矢量组;所述第一组随机矢量组和所述第二组随机矢量组均包括随机配置的第一元素和第二元素,所述第一元素为0,所述第二元素为相同的非零值,所述第一组随机矢量组和第二组随机矢量组的第一元素和第二元素对调;
以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值;
根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数;
根据所述判断参数的大小确认每一成像位置的属性,所述属性包括:组织信号,杂波信号;
依据每一成像位置的属性确定每一成像位置的加权叠加系数;
以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出。
作为本发明一实施方式的进一步改进,所述方法还包括:配置所述第二元素的值为1。
作为本发明一实施方式的进一步改进,以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值包括:配置第一加权参数值以sum1表示,配置第二加权参数值以sum2表示,则:
Figure BDA0002783906460000021
根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数包括:配置判断参数以P表示,则:
Figure BDA0002783906460000022
其中,i=1,2,……N,r(f)表示第一组矢量数据组中序号为f的元素所对应的值,q(f)表示第二组矢量数据组中序号为f的元素所对应的值,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
作为本发明一实施方式的进一步改进,根据所述判断参数的大小确认每一成像位置的属性包括:
根据判断参数与预设参数阈值的大小关系,判断每一成像位置的属性,所述预设参数阈值以P1表示;P1∈(0,0.3);
若P<P1,则当前成像位置的属性为组织信号;
若P之P1,则当前成像位置的属性为杂波信号。
作为本发明一实施方式的进一步改进,依据每一成像位置的属性确定每一成像位置的加权叠加系数之前,所述方法还包括:
配置成像位置属性与加权叠加系数的映射关系;
若当前成像位置的属性为组织信号,则配置其加权叠加系数为Q1;
若当前成像位置的属性为杂波信号,则配置其加权叠加系数为Q2;
其中,Q1=1>Q2;
以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出包括:
配置每一成像位置的输出数据以y表示,则
Figure BDA0002783906460000031
其中,i=1,2,……N,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
为了实现上述发明目的另一,本发明一实施方式提供一种基于随机孔径的心脏超声杂波抑制系统,所述系统包括:数据获取模块,用于获取N个阵元经过延时处理后分别在每一成像位置对齐形成的N个通道信号;
检验模块,用于随机生成N维的第一组随机矢量组和第二组随机矢量组;所述第一组随机矢量组和所述第二组随机矢量组均包括:随机配置、且数量相同的第一元素和第二元素,所述第一元素为0,所述第二元素为相同的非零值,所述第一组随机矢量组和第二组随机矢量组的第一元素和第二元素对调;
处理输出模块,用于以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值;
根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数;
根据所述判断参数的大小确认每一成像位置的属性,所述属性包括:组织信号,杂波信号;
依据每一成像位置的属性确定每一成像位置的加权叠加系数;
以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出。
作为本发明一实施方式的进一步改进,所述检验模块还用于:配置所述第二元素的值为1。
作为本发明一实施方式的进一步改进,所述处理模块以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值时,具体用于:
配置第一加权参数值以sum1表示,配置第二加权参数值以sum2表示,则:
Figure BDA0002783906460000041
所述处理模块根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数时,具体用于:
配置判断参数以P表示,则:
Figure BDA0002783906460000042
其中,i=1,2,……N,r(f)表示第一组矢量数据组中序号为f的元素所对应的值,q(f)表示第二组矢量数据组中序号为f的元素所对应的值,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
作为本发明一实施方式的进一步改进,所述处理模块根据所述判断参数的大小确认每一成像位置的属性时,具体用于:
根据判断参数与预设参数阈值的大小关系,判断每一成像位置的属性,所述预设参数阈值以P1表示;P1∈(0,0.3);
若P<P1,则当前成像位置的属性为组织信号;
若P≥P1,则当前成像位置的属性为杂波信号。
作为本发明一实施方式的进一步改进,所述处理模块在依据每一成像位置的属性确定每一成像位置的加权叠加系数之前,还用于:
配置成像位置属性与加权叠加系数的映射关系;
若当前成像位置的属性为组织信号,则配置其加权叠加系数为Q1;
若当前成像位置的属性为杂波信号,则配置其加权叠加系数为Q2;
其中,Q1=1>Q2;
所述处理模块以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出时,具体用于:
配置每一成像位置的输出数据以y表示,则
Figure BDA0002783906460000051
其中,i=1,2,……N,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
与现有技术相比,本发明的有益效果是:本发明的基于随机孔径的心脏超声杂波抑制方法及系统,能够每一成像位置的属性自适应调整形成该成像位置输出数据的多个通道信号的复合方式,进而抑制成像过程中的杂波噪声,使图像更干净清晰,提升成像效果。
附图说明
图1是本发明一实施方式提供的基于随机孔径的心脏超声杂波抑制方法的流程示意图;
图2是本发明一实施方式提供的基于随机孔径的心脏超声杂波抑制系统的模块示意图。
具体实施方式
以下将结合附图所示的具体实施方式对本发明进行详细描述。但这些实施方式并不限制本发明,本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的变换均包含在本发明的保护范围内。
如图1所示,本发明第一实施方式中提供一种基于随机孔径的心脏超声杂波抑制方法,所述方法包括:
S1、获取N个阵元经过延时处理后分别在每一成像位置对齐形成的N个通道信号;
S2、随机生成N维的第一组随机矢量组和第二组随机矢量组;所述第一组随机矢量组和所述第二组随机矢量组均包括:随机配置、且数量相同的第一元素和第二元素,所述第一元素为0,所述第二元素为相同的非零值,所述第一组随机矢量组和第二组随机矢量组的第一元素和第二元素对调;
S3、以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值;
S4、根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数;
S5、根据所述判断参数的大小确认每一成像位置的属性,所述属性包括:组织信号,杂波信号;
S6、依据每一成像位置的属性确定每一成像位置的加权叠加系数;
以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出。
对于步骤S1,超声成像过程中,超声探头的阵元对组织进行扫查,接收到的超声回波经过模数变换、延时处理等操作后,每一阵元在每一成像位置形成一个通道信号;超声阵元数量较多,通常为96、128等。
较佳的,对于步骤S1,分别将对应每一成像位置的N个通道信号依序排列形成矢量通道信号组;
本发明具体示例中,配置N个阵元在任一成像位置对齐形成的通道信号组为[x1,x2,...,xN],即任一成像位置依次形成的N个通道信号分别以x(f)表示,i=1,2,……N,N为阵元的总数量,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
需要说明的是,最终所需输出的为帧图像,每一成像位置多个阵元的复合数据形成输出帧图像中对应位置的输出数据;以下内容中还会继续描述。
对于步骤S2,本发明可实现方式中,由于杂波以信号形式反映在通道内是随机、杂乱、不规则的,如此,将每一成像位置对齐形成的多个通道信号采用不同孔径视角,则杂波噪声水平差异较大。
基于此,本发明选用截然不同的孔径在每一成像位置分别进行成像,并进一步观察两个孔径的成像结果,如果两个孔径成像的结果类似,则说明该点对应组织信号,予以保留;如果两个孔径成像的结果差异较大,则说明该点大概率是杂波,应予以抑制。
本发明具体实施方式中,第一元素设定为0,在实际应用中,相当于开孔完全关闭;所述第二元素为相同的非零值,相当于开孔打开;本发明其他实施方式中,若以0代表开孔关闭,1代表开孔完全打开,则元素的值处于0-1之间时,表示打开孔的面积占开孔完全打开时的比例,即表示开孔大小。
本发明具体示例中,配置第一元素和第二元素的数量相同,且配置第一元素的值为0,配置所述第二元素的值为1。
本发明一具体实例中,配置第一组随机矢量组为r,该第一组矢量数据组中序号为f的元素所对应的值以r(f)表示;配置第二组随机矢量组为q,该第二组矢量数据组中序号为f的元素所对应的值以q(f)表示,i=1,2,……N。
本发明一具体示例中,r=[1,0,0,0,1,...,1],N维;q=[0,1,1,1,0,...,0],N维。
对于步骤S3,以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值;
配置第一加权参数值以sum1表示,配置第二加权参数值以sum2表示,则:
Figure BDA0002783906460000071
其中,i=1,2,……N,r(f)表示第一组矢量数据组中序号为f的元素所对应的值,q(f)表示第二组矢量数据组中序号为f的元素所对应的值,x(f)表示矢量通道信号组中序号为f的元素所对应的值。对比步骤S4,根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数包括:
配置判断参数以P表示,则:
Figure BDA0002783906460000072
对于步骤S5,根据所述判断参数的大小确认每一成像位置的属性包括:
根据判断参数与预设参数阈值的大小关系,判断每一成像位置的属性,所述预设参数阈值以P1表示;P1∈(0,0.3);
若P<P1,则当前成像位置的属性为组织信号;
若P之P1,则当前成像位置的属性为杂波信号。
本发明一具体示例中,配置P1=0.25;在本发明实现过程中,当判断参数小于预设参数阈值时,在当前成像位置,采用不同孔径进行数据处理后,其获得的信号基本一致,则判定当前位置收到组织信号,需要保留;反之,在当前成像位置,采用不同孔径进行数据处理后,其获得的信号差异较大,则判定当前位置为收到杂波信号,需要抑制。
经过属性判断,可以将一帧数据按照成像位置的属性划分为多个自适应区域;进一步的,可以根据区域不同选择不同的复合加权方式对每一成像位置的多个通道信号进行复合处理,以提升成像效果,该加权复合的方式目的是保持组织信号不变,同时,对杂波信号进行抑制,如此,抑制心脏组织和心腔内的杂波噪声,使图像更干净清晰。
较佳的,在步骤S6之前,所述方法还包括:配置成像位置属性与加权叠加系数的映射关系;
若当前成像位置的属性为组织信号,则配置其加权叠加系数为Q1;
若当前成像位置的属性为杂波信号,则配置其加权叠加系数为Q2;
其中,Q1=1>Q2。
本发明具体示例中,对于组织信号,对其进行保持,因此,配置其加权叠加系数为1,而对于杂波信号,需要抑制,因此,可以在0和1之间选择一个数值配置给Q2;本发明较佳实施方式中,在0至1之间,且接近0的位置选择一个数值配置给Q2,以更利于杂波抑制。
较佳的,配置Q2∈(0,0.3);本发明具体示例中,配置Q2=0.25。
进一步的,对于步骤S6,以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出包括:
配置每一成像位置的输出数据以y表示,则
Figure BDA0002783906460000081
其中,i=1,2,……N,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
需要说明的是,在输出数据的表达公式上,没有考虑额外的变迹系数,如果有使用Hamming等变迹系数,则对结果也没有影响。
在该实施方式中,根据每一成像位置的属性调整其复合加权系数,对于信组织信号,保持其为基础系数,即判断参数小于预设参数阈值时,对于当前成像位置直接输出N个通道信号的叠加信号;对于杂波信号,在基础系数基础上降低其加权系数,以对其进行抑制,即判断参数不小于预设参数阈值时,对于当前成像位置的N个通道信号输出按照配置的加权系数进行加权复合后叠加信号;如此,经过该种方式的加权复合,可以抑制心脏组织和心腔内的杂波噪声,使图像更干净清晰。
结合图2所示,本发明一实施方式提供一种基于随机孔径的心脏超声杂波抑制系统,所述系统包括:数据获取模块100,检验模块200以及处理输出模块300。
数据获取模块100用于获取N个阵元经过延时处理后分别在每一成像位置对齐形成的N个通道信号;
检验模块200用于随机生成N维的第一组随机矢量组和第二组随机矢量组;所述第一组随机矢量组和所述第二组随机矢量组均包括:随机配置、且数量相同的第一元素和第二元素,所述第一元素为0,所述第二元素为相同的非零值,所述第一组随机矢量组和第二组随机矢量组的第一元素和第二元素对调;
处理输出模块300用于以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值;
根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数;
根据所述判断参数的大小确认每一成像位置的属性,所述属性包括:组织信号,杂波信号;
依据每一成像位置的属性确定每一成像位置的加权叠加系数;
以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出。
本发明较佳实施方式中,所述检验模块还用于:配置所述第二元素的值为1。
较佳的,所述处理模块300以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值时,具体用于:
配置第一加权参数值以sum1表示,配置第二加权参数值以sum2表示,则:
Figure BDA0002783906460000101
进一步的,所述处理模块300根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数时,具体用于:
配置判断参数以P表示,则:
Figure BDA0002783906460000102
其中,i=1,2,……N,r(f)表示第一组矢量数据组中序号为f的元素所对应的值,q(f)表示第二组矢量数据组中序号为f的元素所对应的值,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
较佳的,所述处理模块300根据所述判断参数的大小确认每一成像位置的属性时,具体用于:
根据判断参数与预设参数阈值的大小关系,判断每一成像位置的属性,所述预设参数阈值以P1表示;P1∈(0,0.3);
若P<P1,则当前成像位置的属性为组织信号;
若P之P1,则当前成像位置的属性为杂波信号。
较佳的,所述处理模块300在依据每一成像位置的属性确定每一成像位置的加权叠加系数之前,还用于:
配置成像位置属性与加权叠加系数的映射关系;
若当前成像位置的属性为组织信号,则配置其加权叠加系数为Q1;
若当前成像位置的属性为杂波信号,则配置其加权叠加系数为Q2;
其中,Q1=1>Q2;
所述处理模块以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出时,具体用于:
配置每一成像位置的输出数据以y表示,则
Figure BDA0002783906460000111
其中,i=1,2,……N,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统中模块的具体工作过程,可以参考前述方法实施方式中的对应过程,在此不再赘述。
综上所述,本发明的基于随机孔径的心脏超声杂波抑制方法及系统,能够每一成像位置的属性自适应调整形成该成像位置输出数据的多个通道信号的复合方式,进而抑制成像过程中的杂波噪声,使图像更干净清晰,提升成像效果。
以上所描述的系统实施方式仅仅是示意性的,其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件是逻辑模块,即可以位于芯片逻辑中的一个模块中,或者也可以分布到芯片内的多个数据处理模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施方式方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施方式中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于随机孔径的心脏超声杂波抑制方法,其特征在于,所述方法包括:
获取N个阵元经过延时处理后分别在每一成像位置对齐形成的N个通道信号;
随机生成N维的第一组随机矢量组和第二组随机矢量组;所述第一组随机矢量组和所述第二组随机矢量组均包括:随机配置、且数量相同的第一元素和第二元素,所述第一元素为0,所述第二元素为相同的非零值,所述第一组随机矢量组和第二组随机矢量组的第一元素和第二元素对调;
以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值;
根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数;
根据所述判断参数的大小确认每一成像位置的属性,所述属性包括:组织信号,杂波信号;
依据每一成像位置的属性确定每一成像位置的加权叠加系数;
以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出。
2.根据权利要求1所述的基于随机孔径的心脏超声杂波抑制方法,其特征在于,所述方法还包括:配置所述第二元素的值为1。
3.根据权利要求1所述的基于随机孔径的心脏超声杂波抑制方法,其特征在于,以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值包括:配置第一加权参数值以sum1表示,配置第二加权参数值以sum2表示,则:
Figure FDA0002783906450000011
根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数包括:配置判断参数以P表示,则:
Figure FDA0002783906450000012
其中,i=1,2,……N,r(i)表示第一组矢量数据组中序号为i的元素所对应的值,q(i)表示第二组矢量数据组中序号为i的元素所对应的值,x(i)表示矢量通道信号组中序号为i的元素所对应的值。
4.根据权利要求1所述的基于随机孔径的心脏超声杂波抑制方法,其特征在于,根据所述判断参数的大小确认每一成像位置的属性包括:
根据判断参数与预设参数阈值的大小关系,判断每一成像位置的属性,所述预设参数阈值以P1表示;P1∈(0,0.3);
若P<P1,则当前成像位置的属性为组织信号;
若P≥P1,则当前成像位置的属性为杂波信号。
5.根据权利要求1所述的基于随机孔径的心脏超声杂波抑制方法,其特征在于,依据每一成像位置的属性确定每一成像位置的加权叠加系数之前,所述方法还包括:
配置成像位置属性与加权叠加系数的映射关系;
若当前成像位置的属性为组织信号,则配置其加权叠加系数为Q1;
若当前成像位置的属性为杂波信号,则配置其加权叠加系数为Q2;
其中,Q1=1>Q2;
以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出包括:
配置每一成像位置的输出数据以y表示,则
Figure FDA0002783906450000021
其中,i=1,2,……N,x(i)表示矢量通道信号组中序号为i的元素所对应的值。
6.一种基于随机孔径的心脏超声杂波抑制系统,其特征在于,所述系统包括:
数据获取模块,用于获取N个阵元经过延时处理后分别在每一成像位置对齐形成的N个通道信号;
检验模块,用于随机生成N维的第一组随机矢量组和第二组随机矢量组;所述第一组随机矢量组和所述第二组随机矢量组均包括随机配置的第一元素和第二元素,所述第一元素为0,所述第二元素为相同的非零值,所述第一组随机矢量组和第二组随机矢量组的第一元素和第二元素对调;
处理输出模块,用于以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值;
根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数;
根据所述判断参数的大小确认每一成像位置的属性,所述属性包括:组织信号,杂波信号;
依据每一成像位置的属性确定每一成像位置的加权叠加系数;
以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出。
7.根据权利要求6所述的基于随机孔径的心脏超声杂波抑制系统,其特征在于,所述检验模块还用于:配置所述第二元素的值为1。
8.根据权利要求6所述的基于随机孔径的心脏超声杂波抑制系统,其特征在于,所述处理模块以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值时,具体用于:
配置第一加权参数值以sum1表示,配置第二加权参数值以sum2表示,则:
Figure FDA0002783906450000031
所述处理模块根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数时,具体用于:
配置判断参数以P表示,则:
Figure FDA0002783906450000032
其中,i=1,2,……N,r(i)表示第一组矢量数据组中序号为i的元素所对应的值,q(i)表示第二组矢量数据组中序号为i的元素所对应的值,x(i)表示矢量通道信号组中序号为i的元素所对应的值。
9.根据权利要求6所述的基于随机孔径的心脏超声杂波抑制系统,其特征在于,所述处理模块根据所述判断参数的大小确认每一成像位置的属性时,具体用于:
根据判断参数与预设参数阈值的大小关系,判断每一成像位置的属性,所述预设参数阈值以P1表示;P1∈(0,0.3);
若P<P1,则当前成像位置的属性为组织信号;
若P≥P1,则当前成像位置的属性为杂波信号。
10.根据权利要求6所述的基于随机孔径的心脏超声杂波抑制系统,其特征在于,所述处理模块在依据每一成像位置的属性确定每一成像位置的加权叠加系数之前,还用于:
配置成像位置属性与加权叠加系数的映射关系;
若当前成像位置的属性为组织信号,则配置其加权叠加系数为Q1;
若当前成像位置的属性为杂波信号,则配置其加权叠加系数为Q2;
其中,Q1=1>Q2;
所述处理模块以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出时,具体用于:
配置每一成像位置的输出数据以y表示,则
Figure FDA0002783906450000041
其中,i=1,2,……N,x(i)表示矢量通道信号组中序号为i的元素所对应的值。
CN202011291164.6A 2020-11-18 2020-11-18 基于随机孔径的心脏超声杂波抑制方法及系统 Active CN112263278B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011291164.6A CN112263278B (zh) 2020-11-18 2020-11-18 基于随机孔径的心脏超声杂波抑制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011291164.6A CN112263278B (zh) 2020-11-18 2020-11-18 基于随机孔径的心脏超声杂波抑制方法及系统

Publications (2)

Publication Number Publication Date
CN112263278A true CN112263278A (zh) 2021-01-26
CN112263278B CN112263278B (zh) 2022-08-19

Family

ID=74339550

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011291164.6A Active CN112263278B (zh) 2020-11-18 2020-11-18 基于随机孔径的心脏超声杂波抑制方法及系统

Country Status (1)

Country Link
CN (1) CN112263278B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5617862A (en) * 1995-05-02 1997-04-08 Acuson Corporation Method and apparatus for beamformer system with variable aperture
US20060058661A1 (en) * 2004-08-18 2006-03-16 Kabushiki Kaisha Toshiba Ultrasonic diagnostic apparatus and ultrasonic diagnostic method
JP2006349477A (ja) * 2005-06-15 2006-12-28 Mitsubishi Electric Corp レーダ装置
US20120224686A1 (en) * 2011-02-17 2012-09-06 Mazurenko Ivan L Stochastic vector based network echo cancellation
US20160143614A1 (en) * 2013-06-28 2016-05-26 Koninklijke Philips N.V. Rib blockage delineation in anatomically intelligent echocardiography

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5617862A (en) * 1995-05-02 1997-04-08 Acuson Corporation Method and apparatus for beamformer system with variable aperture
US20060058661A1 (en) * 2004-08-18 2006-03-16 Kabushiki Kaisha Toshiba Ultrasonic diagnostic apparatus and ultrasonic diagnostic method
JP2006349477A (ja) * 2005-06-15 2006-12-28 Mitsubishi Electric Corp レーダ装置
US20120224686A1 (en) * 2011-02-17 2012-09-06 Mazurenko Ivan L Stochastic vector based network echo cancellation
US20160143614A1 (en) * 2013-06-28 2016-05-26 Koninklijke Philips N.V. Rib blockage delineation in anatomically intelligent echocardiography

Also Published As

Publication number Publication date
CN112263278B (zh) 2022-08-19

Similar Documents

Publication Publication Date Title
Mateo et al. Finding out general tendencies in speckle noise reduction in ultrasound images
US9451932B2 (en) Clutter suppression in ultrasonic imaging systems
US20130343627A1 (en) Suppression of reverberations and/or clutter in ultrasonic imaging systems
WO2013128301A2 (en) Clutter suppression in ultrasonic imaging systems
WO2018000359A1 (zh) 一种增强超声造影图像的方法、系统及超声造影成像设备
Tasnim et al. Study of speckle noise reduction from ultrasound B-mode images using different filtering techniques
JP2022023982A (ja) 超音波画像を処理するための方法及びシステム
CN112263278B (zh) 基于随机孔径的心脏超声杂波抑制方法及系统
CN109674493A (zh) 医用超声自动追踪颈动脉血管的方法、系统及设备
Touil et al. Analysis of motion tracking in echocardiographic image sequences: Influence of system geometry and point-spread function
CN111815527B (zh) 基于Weibull分布的混合高阶变分超声图像去噪方法
CN110731795A (zh) 空间复合成像的处理方法和装置
Hasan et al. Blind deconvolution of ultrasound images using $ l_ {1} $-norm-constrained block-based damped variable step-size multichannel LMS algorithm
Yao et al. Spatial compounding of large numbers of multi-view 3D echocardiography images using feature consistency
Toffali et al. Improving the Quality of Monostatic Synthetic-Aperture Ultrasound Imaging Through Deep-Learning-Based Beamforming
CN110869799B (zh) 用于处理超声图像的方法和系统
CN110431443B (zh) 用于对超声图像杂波进行滤波的方法和系统
CN112150370A (zh) 一种空间复合成像方法和装置
CN111340805A (zh) B超图像处理装置、脂肪肝b超图像处理装置及b超系统
Lin et al. Ultrasound image compounding based on motion compensation
Apostolakis et al. Adaptive Slow-time Singular Value Thresholding (SVT) based on Stein's Unbiased Risk Estimate (SURE) for Ultrasound Image Random Noise Reduction
Kleckler et al. Characterization of Heterogeneous Perfusion in Contrast-Enhanced Ultrasound
Wu et al. Enhancing the image resolution of acoustic microscopy with FISTA
Yao et al. Spatial compounding of large sets of 3D echocardiography images
Matrone et al. Enhanced ultrasound harmonic imaging using the filtered-delay multiply and sum beamformer

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Applicant after: Feiyinuo Technology (Suzhou) Co.,Ltd.

Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Applicant before: VINNO TECHNOLOGY (SUZHOU) Co.,Ltd.

Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Applicant after: Feiyinuo Technology Co.,Ltd.

Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Applicant before: Feiyinuo Technology (Suzhou) Co.,Ltd.

GR01 Patent grant
GR01 Patent grant