发明内容
为解决上述技术问题,本发明的目的在于提供一种基于随机孔径的心脏超声杂波抑制方法及系统。
为了实现上述发明目的之一,本发明一实施方式提供一种基于随机孔径的心脏超声杂波抑制方法,所述方法包括:获取N个阵元经过延时处理后分别在每一成像位置对齐形成的N个通道信号;
随机生成N维的第一组随机矢量组和第二组随机矢量组;所述第一组随机矢量组和所述第二组随机矢量组均包括随机配置的第一元素和第二元素,所述第一元素为0,所述第二元素为相同的非零值,所述第一组随机矢量组和第二组随机矢量组的第一元素和第二元素对调;
以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值;
根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数;
根据所述判断参数的大小确认每一成像位置的属性,所述属性包括:组织信号,杂波信号;
依据每一成像位置的属性确定每一成像位置的加权叠加系数;
以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出。
作为本发明一实施方式的进一步改进,所述方法还包括:配置所述第二元素的值为1。
作为本发明一实施方式的进一步改进,以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值包括:配置第一加权参数值以sum1表示,配置第二加权参数值以sum2表示,则:
根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数包括:配置判断参数以P表示,则:
其中,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表示,则
其中,i=1,2,……N,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
为了实现上述发明目的另一,本发明一实施方式提供一种基于随机孔径的心脏超声杂波抑制系统,所述系统包括:数据获取模块,用于获取N个阵元经过延时处理后分别在每一成像位置对齐形成的N个通道信号;
检验模块,用于随机生成N维的第一组随机矢量组和第二组随机矢量组;所述第一组随机矢量组和所述第二组随机矢量组均包括:随机配置、且数量相同的第一元素和第二元素,所述第一元素为0,所述第二元素为相同的非零值,所述第一组随机矢量组和第二组随机矢量组的第一元素和第二元素对调;
处理输出模块,用于以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值;
根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数;
根据所述判断参数的大小确认每一成像位置的属性,所述属性包括:组织信号,杂波信号;
依据每一成像位置的属性确定每一成像位置的加权叠加系数;
以所述加权叠加系数对每一成像位置所对应的N个通道信号进行叠加后进行输出。
作为本发明一实施方式的进一步改进,所述检验模块还用于:配置所述第二元素的值为1。
作为本发明一实施方式的进一步改进,所述处理模块以第一组随机矢量和第二组随机矢量分别对每一成像位置的N个通道信号进行加权求和形成第一加权参数值和第二加权参数值时,具体用于:
配置第一加权参数值以sum1表示,配置第二加权参数值以sum2表示,则:
所述处理模块根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数时,具体用于:
其中,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表示,则
其中,i=1,2,……N,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
与现有技术相比,本发明的有益效果是:本发明的基于随机孔径的心脏超声杂波抑制方法及系统,能够每一成像位置的属性自适应调整形成该成像位置输出数据的多个通道信号的复合方式,进而抑制成像过程中的杂波噪声,使图像更干净清晰,提升成像效果。
具体实施方式
以下将结合附图所示的具体实施方式对本发明进行详细描述。但这些实施方式并不限制本发明,本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的变换均包含在本发明的保护范围内。
如图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表示,则:
其中,i=1,2,……N,r(f)表示第一组矢量数据组中序号为f的元素所对应的值,q(f)表示第二组矢量数据组中序号为f的元素所对应的值,x(f)表示矢量通道信号组中序号为f的元素所对应的值。对比步骤S4,根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数包括:
对于步骤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表示,则
其中,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表示,则:
进一步的,所述处理模块300根据第一加权参数值和第二加权参数值获取对应每一成像位置的判断参数时,具体用于:
其中,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表示,则
其中,i=1,2,……N,x(f)表示矢量通道信号组中序号为f的元素所对应的值。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统中模块的具体工作过程,可以参考前述方法实施方式中的对应过程,在此不再赘述。
综上所述,本发明的基于随机孔径的心脏超声杂波抑制方法及系统,能够每一成像位置的属性自适应调整形成该成像位置输出数据的多个通道信号的复合方式,进而抑制成像过程中的杂波噪声,使图像更干净清晰,提升成像效果。
以上所描述的系统实施方式仅仅是示意性的,其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件是逻辑模块,即可以位于芯片逻辑中的一个模块中,或者也可以分布到芯片内的多个数据处理模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施方式方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施方式中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。