CN114129185B - 波束合成方法、超声成像方法、装置及设备 - Google Patents
波束合成方法、超声成像方法、装置及设备 Download PDFInfo
- Publication number
- CN114129185B CN114129185B CN202111470011.2A CN202111470011A CN114129185B CN 114129185 B CN114129185 B CN 114129185B CN 202111470011 A CN202111470011 A CN 202111470011A CN 114129185 B CN114129185 B CN 114129185B
- Authority
- CN
- China
- Prior art keywords
- signal
- signals
- distance
- line
- depth value
- 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
- 238000000034 method Methods 0.000 title claims abstract description 89
- 238000003384 imaging method Methods 0.000 title claims abstract description 20
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 20
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 20
- 238000012937 correction Methods 0.000 claims description 28
- 238000012545 processing Methods 0.000 claims description 25
- 230000002194 synthesizing effect Effects 0.000 claims description 19
- 238000000926 separation method Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 8
- 238000006243 chemical reaction Methods 0.000 description 18
- 230000005540 biological transmission Effects 0.000 description 14
- 238000004364 calculation method Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 9
- 230000011218 segmentation Effects 0.000 description 9
- 238000001914 filtration Methods 0.000 description 8
- 210000005075 mammary gland Anatomy 0.000 description 4
- 210000000481 breast Anatomy 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 101100188603 Arabidopsis thaliana OFP1 gene Proteins 0.000 description 2
- 101100188604 Arabidopsis thaliana OFP2 gene Proteins 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本申请提供了一种波束合成方法、超声成像方法、装置及设备,包括如下步骤:通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z;根据校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将BF1信号集中的各BF1信号进行指定叠加得到BF2信号。通过实现多波束的信号对准,能够较好地抑制旁瓣,提高图像分辨率及信噪比;通过多波束叠加极大的提高线密度,提高了图像的空间分辨率;装置通过在波束合成模块中增加了校准模块和合成模块,使波束合成后数据更精确;并且通过声场判断,对有效声场数据叠加,相较于传统方法实现了较好的聚焦效果,提高了图像质量。
Description
本申请为申请日为2019年03月18日、申请号为201910204036.4、名称为“波束合成方法、超声成像方法、装置及设备”的发明专利申请的分案申请。
技术领域
本申请涉及计算机技术领域,特别是涉及一种波束合成方法、超声成像方法、装置及设备。
背景技术
超声成像因其具有安全、实时、便携、无创及成本低等优势,被广泛应用于临床医学诊断。而波束合成是超声成像系统中处于核心位置,对成像质量起着决定性的作用。波束合成中主瓣宽度和旁瓣幅度用来判断所形成波束质量的高低,通常主瓣的宽度越窄,那么成像的横向分辨率越高;旁瓣的幅度越小,那么成像的对比度就越大,并且伪影噪声越少。
另外波束形成接收分为单波束和多波束接收,在运动器官的实时诊断中,多波束接收可以极大的提高帧率,然而每次发射接收的多波束之间相位不对齐,这样叠加得到的RF信号同样具有较高的伪影噪声。
传统方法通过进行延时叠加波束,图像质量差,旁瓣等级较高,空间分辨率低,虽然可以通过动态聚焦、幅度变迹和动态孔径三种方法来控制波束主瓣宽度与旁瓣幅度,从而抑制旁瓣改善超声成像的质量,但是同时也抑制主瓣,对图像分辨率有影响;而现有的其他方法亦是以降低旁瓣幅度,提高分辨率,提高信噪比等方面改善聚焦效果,但提高帧率,改善空间分辨率的方法仍然是现有的技术盲区。
发明内容
鉴于上述问题,提出了本发明实施例以便提供一种克服上述问题或者至少部分地解决上述问题的波束合成的方法及装置。
为了解决上述问题,本发明实施例公开了一种波束合成方法,包括如下步骤:
通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z;
根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号。
进一步地,上述根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号的步骤之前,还包括如下步骤:
对上述BF1信号集进行有效信号过滤,并区分标记出上述BF1信号集中的有效BF1信号和无效BF1信号。
进一步地,上述通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z的步骤,包括如下步骤:
通过阵元数N_elements、相邻阵元间距pitch、发射线数nlines和移动步长step换算出线阵中各波束线的间距dx;
通过线阵中各波束线的间距dx、每次发射或接收的波束数量beam、发射线上样本点到发射点的距离z_sample1和焦距F分别换算出中心发射线上的样本点到各条波束线上的虚拟焦点的横向距离Dxibeam和纵向距离Dyibeam;
根据中心发射线上的样本点到各条波束线上的虚拟焦点的横向距离Dxibeam和纵向距离Dyibeam换算出中心发射线上的样本点到各波束线虚拟焦点的距离dr1ibeam;
根据中心发射线上的样本点到各波束线虚拟焦点的距离dr1ibeam和焦距F换算出各波束线对应的校准后的聚焦深度值Z2。
进一步地,上述对上述BF1信号集进行有效信号过滤,并区分标记出上述BF1信号集中的有效BF1信号和无效BF1信号的步骤,包括如下步骤:
通过超声波发射器的孔径A、通道数N_channel、相邻阵元间距pitch和焦距F换算出发射声场的角度θ1;
计算出各波束线上各样本点偏离发射线的角度angle1ibeam;
判断上述角度angle1ibeam是否大于θ1/2;
若否,则判定与上述角度angle1ibeam对应的样本点处于有效声场范围内。
进一步地,上述根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号的步骤,包括如下步骤:
通过接收线的坐标ref1、阵元坐标xe1和各波束线的校准后的聚焦深度值Z2换算出各阵元到聚焦点的距离dn1ibeam;
通过接收线的坐标ref1、阵元坐标xe1、各波束线的校准后的聚焦深度值Z2和各阵元到聚焦点的距离dn1ibeam换算出各阵元的延时τ1;
通过各阵元的延时τ1将与之对应的波束线进行相位对齐,并将信号对准后的波束线的BF1信号集进行叠加得到上述BF2信号。
进一步地,上述通过各阵元的延时τ1将与之对应的波束线进行相位对齐,并将信号对准后的波束线的BF1信号集进行叠加得到上述BF2信号的步骤,包括如下步骤:
根据各阵元的延时τ1和处于有效声场范围内的上述角度angle1ibeam换算出对应样本点的回波信号;
将上述对应样本点的回波信号汇集形成上述信号对准后的BF1信号集。
进一步地,上述通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z的步骤,包括如下步骤:
根据阵元数N_elements、相邻阵元间距pitch和凸阵圆心半径R换算出凸阵圆心角beta;
根据发射线数nlines、移动步长step和上述凸阵圆心角beta换算出接收线间角dx_beta;
根据每次发射或接收的波束数量beam和上述接收线间角dx_beta换算出接收波束到发射波束的夹角Angle_beam;
根据焦距F、上述凸阵圆心半径R和接收波束到发射波束的夹角Angle_beam分别换算出样本点到接收波束虚拟焦点的横向距离Hxibeam和纵向距离Hyibeam,并根据上述上述凸阵圆心半径R和纵向距离Hyibeam换算出校正焦距FA;
根据发射线上样本点到圆心的距离z_sample2、上述凸阵圆心半径R、横向距离Hxibeam和纵向距离Hyibeam换算出中心发射线上的样本点到各波束线虚拟焦点的距离dr2ibeam;
根据上述中心发射线上的样本点到各波束线虚拟焦点的距离dr2ibeam、焦距F和校正焦距FA换算出各波束线对应的校准后的聚焦深度值Z2。
进一步地,上述对上述BF1信号集进行有效信号过滤,并区分标记出上述BF1信号集中的有效BF1信号和无效BF1信号的步骤,包括如下步骤:
根据凸阵圆心半径R、通道数N_channel、焦距F和相邻阵元间距pitch换算出发射声场的角度θ2;
获取各波束线上样本点偏离发射线的角度angle2ibeam;
判断上述角度angle2ibeam是否大于θ2/2
若否,则判定与上述角度angle2ibeam对应的样本点处于有效声场范围内。
进一步地,上述根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号的步骤,包括如下步骤:
根据接收线的起始位置坐标M、阵元的位置坐标N、上述校准后的聚焦深度值Z2换算出阵元到聚焦点的距离dn2ibeam;
根据上述校准后的聚焦深度值Z2和阵元到聚焦点的距离dn2ibeam换算出各阵元的延时τ2;
通过各阵元的延时τ2将与之对应的波束线进行相位对齐,并将信号对准后的波束线的BF1信号集进行叠加得到上述BF2信号。
进一步地,上述通过各阵元的延时τ2将与之对应的波束线进行相位对齐,并将信号对准后的波束线的BF1信号集进行叠加得到上述BF2信号的步骤,包括如下步骤:
根据各阵元的延时τ2和处于有效声场范围内的上述角度angle2ibeam换算出对应样本点的回波信号;
将上述对应样本点的回波信号汇集形成上述信号对准后的BF1信号集。
为了解决上述问题,本发明实施例公开了一种超声成像方法,包括如下步骤,
获取超声波数据并转化为相应的数字信号;
将上述数字信号通过指定的校正合成为射频信号;
将上述射频信号通过指定的信号处理分离出载波信号;
将上述载波信号通过指定的图像处理获得超声图像;
上述将上述数字信号通过指定的校正合成为射频信号的步骤,包括上述任一实施例所述的波束合成方法。
为了解决上述问题,本发明实施例公开了一种波束合成装置,包括如下具体模块:
校准模块,用于通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z;
合成模块,用于根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号。
为了解决上述问题,本发明实施例公开了一种超声成像装置,包括如下具体模块:
获取模块,用于获取超声波数据并转化为相应的数字信号;
波束合成模块,用于将上述数字信号通过指定的校正合成为射频信号;
分离模块,用于将上述射频信号通过指定的信号处理分离出载波信号;
成像模块,用于将上述载波信号通过指定的图像处理获得超声图像,
上述波束合成模块,包括上述任一实施例所述的波束合成装置。
为了解决上述问题,本发明实施例公开了一种超声成像设备,包括产生超声波信号辐射到被测试组织中和吸收反射声波信号的信号收发器,将接收的声波转换为电信号的振动传感器,将接收信号采样和数字化的模数转换器(A/D),补偿由于深度导致的超声波幅值衰减的时间增益补偿器(TGC),将数字信号转化为RF信号的接收波束合成器,将RF信号进行包络提取和解调处理分离出载波信号的信号处理器,将载波信号进行扫描转换和后端图像处理的数字扫描转换器DSC以及用于最终显示的图像的显示器。
其中,上述接收波束合成器,包括相位校正模块、用于合成BF1信号的第一次波束波束合成模块、有效声场判断模块和用于合成BF2信号的第二次波束合成器,且上述相位校正模块、用于合成BF1信号的第一次波束波束合成模块、有效声场判断模块和用于合成BF2信号的第二次波束合成器依次相连。
为了解决上述问题,本发明实施例公开了一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现本发明实施例中任意一项所述的方法。
与现有技术相比,本申请包括以下优点:
本发明实施例中,通过实现多波束的空间对准和相位对齐,能够较好地抑制旁瓣,提高图像分辨率及信噪比;通过多波束叠加极大的提高线密度,提高了图像的空间分辨率;本方法实现算法简单,在消耗少量处理资源的情况下即可实现多帧率,高信噪比,高分辨率的效果;装置通过在波束合成模块中增加了校准模块和合成模块,使波束合成后更数据更精确;并且通过声场判断,对有效声场数据叠加,相较于传统方法实现了较好的聚焦效果,提高了图像质量。
附图说明
图1是本发明一实施例的波束合成方法的步骤流程示意图;
图2是本发明一实施例的波束合成方法的步骤流程示意图;
图3是本发明一实施例的波束合成方法的步骤流程示意图;
图4是本发明一实施例的波束合成方法的步骤流程示意图;
图5是本发明一实施例的线阵聚焦深度校正计算示意图;
图6是本发明一实施例的线阵有效声场计算示意图;
图7是本发明一实施例的线阵延时距离计算示意图;
图8是本发明一实施例的凸阵聚焦深度校正计算示意图;
图9是本发明一实施例的凸阵有效声场计算示意图;
图10是本发明一实施例的凸阵延时距离计算示意图;
图11是本发明一实施例的波束合成装置的模块结构示意图;
图12是本发明一实施例的一种计算机设备的结构示意图。
1、校准模块;2、合成模块;12、计算机设备;14、外部设备;16、处理单元;18、总线;20、网络适配器;22、(I/O)接口;24、显示器;28、系统存储器;30、随机存取存储器(RAM);32、高速缓存存储器;34、存储系统;40、程序/实用工具;42、程序模块。
本发明目的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
为使本申请的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本申请作进一步详细的说明。
显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
另外,各个实施例之间的技术方案可以相互结合,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在,也不在本发明要求的保护范围之内。
最后,本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。
需要说明的是,在本发明任一实施例中,信号对准包含空间对准和相位对齐;空间对准:指接收线位置相同;相位对齐:是指不同接收线(空间对准后的接收线)采用不同的延时时间,使信号能够同时到达到空间位置相同的一点(即聚焦到一点),通过校正使信号能够聚集到同一位置。
需要说明的是,为方便公式的表达,在下列任意实施例中所提及的计算参数“Xibeam”表示为参数任一波束线i对应的“X”,Xibeam(j)表示为任一波束线i中的任一点j对应的“X”,如:将任一波束线对应的校准后的聚焦深度值Z2以Z2ibeam在公式中进行表示。
参照图1,示出了本申请一种波束合成方法,包括如下步骤:
S1、通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z;
S2、根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号。
如上述步骤S1所述,通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z,需要说明的是,对所述聚焦深度值进行校准的方式在本发明实施例中包括两种,分别为根据线阵校正法和凸阵校正法,其中,线阵校正法通过阵元数N_elements、相邻阵元间距pitch、发射线数nlines、移动步长step、每次发射或接收的波束数量beam、发射线上样本点到发射点的距离z_sample1和焦距F对所述聚焦深度值Z进行校正;凸阵校正法通过根据阵元数N_elements、相邻阵元间距pitch、凸阵圆心半径R、发射线数nlines、移动步长step、每次发射或接收的波束数量beam、焦距F、发射线上样本点到圆心的距离z_sample2对所述聚焦深度值Z进行校正。
如上述步骤S2所述,根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号,需要说明的是,步骤S2的实现方法包括两种,分别为根据线阵校准法和凸阵校准法,其中,线阵校准法通过接收线的坐标ref1、阵元坐标xe1、校准后的聚焦深度值Z2将对应的波束线进行相位对齐,并将信号对准后的BF1信号集进行叠加得到所述BF2信号;凸阵校准发通过接收线的起始位置坐标M、阵元的位置坐标N、所述校准后的聚焦深度值Z2将对应的波束线进行相位对齐,并将信号对准后的BF1信号集进行叠加得到所述BF2信号。在通过校准后的聚焦深度值Z2将对应的波束线进行相位对齐前,波束线已经空间对准完毕,通过上述校准后的聚焦深度值Z2对波束进行延时。得到BF1,判断数据的有效性,再叠加得到BF2由于BF1信号的相位已经对齐,所以只需同一空间数据相叠加。
参照图2-4,示出了本申请一种波束合成方法实施例2的步骤流程图,具体可以包括如下步骤:
在本发明实施例中,上述根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号的步骤之前,还包括如下步骤:
S3、对上述BF1信号集进行有效信号过滤,并区分标记出上述BF1信号集中的有效BF1信号和无效BF1信号。
如上述步骤S3所述,对上述BF1信号集进行有效信号过滤,并区分标记出上述BF1信号集中的有效BF1信号和无效BF1信号需要说明的是,步骤S3的实现方法包括两种,分别为根据线阵校准法和凸阵校准法,其中,线阵校准法通过通过超声波发射器的孔径A、通道数N_channel、相邻阵元间距pitch、焦距F和各波束线上各样本点偏离发射线的角度angle1ibeam计算并筛选出有效声场范围内的BF1信号;凸阵校准法根据凸阵圆心半径R、通道数N_channel、焦距F、相邻阵元间距pitch和各波束线上样本点偏离发射线的角度angle2ibeam计算并筛选出有效声场范围内的BF1信号。
在本发明实施例中,上述通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z的步骤,包括如下步骤:
S111、通过阵元数N_elements、相邻阵元间距pitch、发射线数nlines和移动步长step换算出线阵中各波束线的间距dx;
S112、通过线阵中各波束线的间距dx、每次发射或接收的波束数量beam、发射线上样本点到发射点的距离z_sample1和焦距F分别换算出中心发射线上的样本点到各条波束线上的虚拟焦点的横向距离Dxibeam和纵向距离Dyibeam;
S113、根据中心发射线上的样本点到各条波束线上的虚拟焦点的横向距离Dxibeam和纵向距离Dyibeam换算出中心发射线上的样本点到各波束线虚拟焦点的距离dr1ibeam;
S114、根据中心发射线上的样本点到各波束线虚拟焦点的距离dr1ibeam和焦距F换算出各波束线对应的校准后的聚焦深度值Z2。
如上述步骤S111所述,通过阵元数N_elements、相邻阵元间距pitch、发射线数nlines和移动步长step换算出线阵中各波束线的间距dx,需要说明的是,上述计算过程所采用的公式在本发明实施例中优选为:
其中,当移动步长step为1时,接收线间距等于发射线间距。step与最终的线密度相关,当step越大线密度越高,但是多波束叠加的数量就减少,具体关系如下:
波束叠加数×step=beam
式中,beam为每次发射或接收的波束数量。
如上述步骤S112所述,通过线阵中各波束线的间距dx、每次发射或接收的波束数量beam、发射线上样本点到发射点的距离z_sample1和焦距F分别换算出中心发射线上的样本点到各条波束线上的虚拟焦点的横向距离Dxibeam和纵向距离Dyibeam,需要说明的是,上述计算过程所采用的公式在本发明实施例中优选为:
Dx=abs(-(beam-1)×dx×1/2+(ibeam-1)*dx)ibeam=1,2.....beam公式(2)
Dy=abs(z_sample-F)
式中,abs代表为式中的计算结果取绝对值。
如上述步骤S113所述,根据中心发射线上的样本点到各条波束线上的虚拟焦点的横向距离Dxibeam和纵向距离Dyibeam换算出中心发射线上的样本点到各波束线虚拟焦点的距离dr1ibeam,需要说明的是,上述计算过程所采用的公式在本发明实施例中优选为:
以图5为例:图中P1,Pn为发射线上近场和远场的两个样本点,以样本点P1为例,接收波束线以bm2为例进行说明,S为bm2上与发射焦点同位置的虚拟焦点。S到发射线上每个样本点的距离为drbm2,Dx1为发射线上任意样本点P到接收虚拟焦点S的横向距离;Dy1为纵向深度。通过公式(3),即可计算得到发射线上样本点P1到波束线bm2虚拟焦点S的距离,即dr1。
如上述步骤S114所述,根据中心发射线上的样本点到各波束线虚拟焦点的距离dr1ibeam和焦距F换算出各波束线对应的校准后的聚焦深度值Z2,需要说明的是,上述计算过程所采用的公式在本发明实施例中优选为:
Z2ibeam(j)=F±dribeam(j) 公式(4)
需要说明的是,式中,当z_sample1小于F使用’-’;当z_sample1大于F使用’+’;当z_sample1等于F时,即焦点位置。此时焦点处声场集中,受到其他发射声场的影响较小,因此,当z_sample1等于F时,仅采用当前线焦点处的数据。
通过上述公式(4)即可得到任一条将要接收的波束线ibeam上任一个点j的校正聚焦深度值Z2,延时计算时,使用对应的Z2值,即可接收到相位一致的多波束BF1。
在本发明实施例中,上述对上述BF1信号集进行有效信号过滤,并区分标记出上述BF1信号集中的有效BF1信号和无效BF1信号的步骤,包括如下步骤:
S311、通过超声波发射器的孔径A、通道数N_channel、相邻阵元间距pitch和焦距F换算出发射声场的角度θ1;
S312、计算出各波束线上各样本点偏离发射线的角度angle1ibeam;
S313、判断上述角度angle1ibeam是否大于θ1/2;
S314、若否,则判定与上述角度angle1ibeam对应的样本点处于有效声场范围内。
如上述步骤S311所述,通过超声波发射器的孔径A、通道数N_channel、相邻阵元间距pitch和焦距F换算出发射声场的角度θ1,需要说明的是,线阵的声场形状类似沙漏状,如图6所示。以孔径大小为A发射超声波,进行多波束接收,其中P2,P3为指定的一条接收波束线上的两个样本点,从P2,P3到中心发射线做垂线,垂点为N1,N2。
因此发射声场的角度优选通过下列公式计算:
如上述步骤S312所述,计算出各波束线上各样本点偏离发射线的角度angle1ibeam,需要说明的是,上述计算过程所采用的公式在本发明实施例中优选为:
如上述步骤S313所述,判断上述角度angle1ibeam是否大于θ1/2,如图6所示,可知P2不在发射声场之内,其与中心发射线的夹角angle11大于θ1/2;样本点P3在发射声场之内,其与中心发射线的夹角angle12小于θ1/2。因此可以通过判断波束线上各点与中心发射线焦点的夹角判断各点是否在有效声场之内。
如上述步骤S314所述,若否,则判定与上述角度angle1ibeam对应的样本点处于有效声场范围内,具体地,判断后对对应的样本点进行标记,具体标记为:
在本发明实施例中,上述根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号的步骤,包括如下步骤:
S211、通过接收线的坐标ref1、阵元坐标xe1和各波束线的校准后的聚焦深度值Z2换算出各阵元到聚焦点的距离dn1ibeam;
S212、通过接收线的坐标ref1、阵元坐标xe1、各波束线的校准后的聚焦深度值Z2和各阵元到聚焦点的距离dn1ibeam换算出各阵元的延时τ1;
S213、通过各阵元的延时τ1将与之对应的波束线进行相位对齐,并将信号对准后的波束线的BF1信号集进行叠加得到上述BF2信号。
如上述步骤S211所述,通过接收线的坐标ref1、阵元坐标xe1和各波束线的校准后的聚焦深度值Z2换算出各阵元到聚焦点的距离dn1ibeam,需要说明的是,延时距离的计算方法如图7所示,样本点P0在波束线上移动,dn为各阵元到聚焦点距离,因此,上述换算过程所采用的公式在本发明实施例中优选为:
如上述步骤S212所述,通过接收线的坐标ref1、阵元坐标xe1、各波束线的校准后的聚焦深度值Z2和各阵元到聚焦点的距离dn1ibeam换算出各阵元的延时τ1,需要说明的是,上述换算过程所采用的公式在本发明实施例中优选为:
式中,c为组织中的平均声速,c=1540m/s。
如上述步骤S213所述,通过各阵元的延时τ1将与之对应的波束线进行相位对齐,并将信号对准后的波束线的BF1信号集进行叠加得到上述BF2信号,需要说明的是,在BF1信号进行相位对齐的步骤具体为如下过程:
S213-1、根据各阵元的延时τ1和处于有效声场范围内的上述角度angle1ibeam换算出对应样本点的回波信号,需要说明的是,上述换算过程所采用的公式在本发明实施例中优选为:
式中,S1ibeam(j)为第ibeam条接收波束线上第j个样本点的接收回波信号,s1j(τ1)为子孔径中阵元j接收到的回波,Wibeam(j)为信号幅度变迹系数。
通过上述公式(10)即可得到相位对齐的多波束回波信号,即,BF1信号。其中,BF1信号的线数与beam有关,发射1次接收beam次,因此发射nlines次,共接收到得到nlines×beam条波束。
需要说明的是,在将信号对准后的BF1信号集进行叠加得到上述BF2信号之前,还包括如下步骤:
S213-2、将上述对应样本点的回波信号汇集形成上述信号对准后的BF1信号集,将nlines×beam条波束的BF1信号,汇集为BF1信号集,该BF1信号集即为本次波束合成过程中所有有效声场中的BF1信号的总集。
当汇集完BF1信号集后对信号集中的BF1信号进行叠加最后形成BF2信号,由于BF1信号的相位已经对齐,所以只需同一空间数据相叠加即可获得BF2信号。
在本发明实施例中,上述通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z的步骤,包括如下步骤:
S121、根据阵元数N_elements、相邻阵元间距pitch和凸阵圆心半径R换算出凸阵圆心角beta;
S122、根据发射线数nlines、移动步长step和上述凸阵圆心角beta换算出接收线间角dx_beta;
S123、根据每次发射或接收的波束数量beam和上述接收线间角dx_beta换算出接收波束到发射波束的夹角Angle_beam;
S124、根据焦距F、上述凸阵圆心半径R和接收波束到发射波束的夹角Angle_beam分别换算出样本点到接收波束虚拟焦点的横向距离Hxibeam和纵向距离Hyibeam,并根据上述上述凸阵圆心半径R和纵向距离Hyibeam换算出校正焦距FA;
S125、根据发射线上样本点到圆心的距离z_sample2、上述凸阵圆心半径R、横向距离Hxibeam和纵向距离Hyibeam换算出中心发射线上的样本点到各波束线虚拟焦点的距离dr2ibeam;
S126、根据上述中心发射线上的样本点到各波束线虚拟焦点的距离dr2ibeam、焦距F和校正焦距FA换算出各波束线对应的校准后的聚焦深度值Z2。
如上述步骤S121所述,根据阵元数N_elements、相邻阵元间距pitch和凸阵圆心半径R换算出凸阵圆心角beta,需要说明的是,上述换算过程所采用的公式在本发明实施例中优选为:
如上述步骤S122所述,根据发射线数nlines、移动步长step和上述凸阵圆心角beta换算出接收线间角dx_beta,需要说明的是,上述换算过程所采用的公式在本发明实施例中优选为:
需要说明的是,同线阵,当step=1时,发射线间角等于接收线间角,step越大,接收线间角越小,则线密度越高,在beam一定时,则叠加的波束数越少。一般线密度越高,叠加次数越多聚焦效果越好,但是step与叠加次数呈反比关系,因此实际应用中一般会根据实际情况对step、beam、叠加数目进行权衡,以得到较好的聚焦效果。
如上述步骤S123所述,根据每次发射或接收的波束数量beam和上述接收线间角dx_beta换算出接收波束到发射波束的夹角Angle_beam,需要说明的是,上述换算过程所采用的公式在本发明实施例中优选为:
(ibeam=1,2.....beam)。
如上述步骤S124所述,根据焦距F、上述凸阵圆心半径R和接收波束到发射波束的夹角Angle_beam分别换算出样本点到接收波束虚拟焦点的横向距离Hxibeam和纵向距离Hyibeam,并根据上述上述凸阵圆心半径R和纵向距离Hyibeam换算出校正焦距FA,需要说明的是,换算出样本点到接收波束虚拟焦点的横向Hxibeam、纵向距离Hyibeam所采用的公式在本发明实施例中优选为:
Hxibeam=|(F+R)×cos(Angle_beamibeam)|
Hyibeam=|(F+R)×sin(Angle_beamibeam)| 公式(14)。
如图8所示,图中以bm4为例,P4,P5为发射波束线上近场和远场的两个样本点,接收波束虚拟焦点为S,F_curve为焦距圆弧,从F_curve的虚拟焦点S上向中心发射波束线做垂线,垂点为F_Al,样本点到S点的横纵距离如图8所示;由于相对于线阵而言Hy恒定不变为焦距,凸阵虚拟焦点垂直中心发射线,垂点为校正的焦点,则校正焦距为FA,因此
FAibeam=Hyibeam-R公式(15)。
如上述步骤S125所述,根据发射线上样本点到圆心的距离z_sample2、上述凸阵圆心半径R、横向距离Hxibeam和纵向距离Hyibeam换算出中心发射线上的样本点到各波束线虚拟焦点的距离dr2ibeam,需要说明的是,上述换算过程所采用的公式在本发明实施例中优选为:
如上述步骤S126所述,根据上述中心发射线上的样本点到各波束线虚拟焦点的距离dr2ibeam、焦距F和校正焦距FA换算出各波束线对应的校准后的聚焦深度值Z2,需要说明的是,上述换算过程所采用的公式在本发明实施例中优选为:
在本发明实施例中,上述对上述BF1信号集进行有效信号过滤,并区分标记出上述BF1信号集中的有效BF1信号和无效BF1信号的步骤,包括如下步骤:
S321、根据凸阵圆心半径R、通道数N_channel、焦距F和相邻阵元间距pitch换算出发射声场的角度θ2;
S322、获取各波束线上样本点偏离发射线的角度angle2ibeam;
S323、判断上述角度angle2ibeam是否大于θ2/2;
S324、若否,则判定与上述角度angle2ibeam对应的样本点处于有效声场范围内。
如上述步骤S321所述,根据凸阵圆心半径R、通道数N_channel、焦距F和相邻阵元间距pitch换算出发射声场的角度θ2,需要说明的是,通过超声波发射器的孔径A、通道数N_channel、相邻阵元间距pitch和焦距F换算出发射声场的角度θ2,需要说明的是,凸阵的声场如线阵为沙漏状,如图9所示,以孔径A发射超声波,多波束接收。其中,从焦点F到中心发射线做垂线,垂点为N3。
因此发射声场的角度优选通过下列公式计算:
如上述步骤S322所述,获取各波束线上样本点偏离发射线的角度angle2ibeam,需要说明的是,波束线上样本点偏离发射线的角度angle2优选通过以下方法计算获取,通过例子对以下具体步骤进行解释,具体为:
①当样本点位于垂点N3之前
即计算三角形OFP1中,∠OFP1为angle21,其中令OF=a1,FP1=b1,OP1=c1
a1=F+R
b1=dr(1:markibeam,ibeam)ibeam=1,2,3.......beam
c1=z_sample+R
mark为各波束线上校正焦距的点坐标。
②样本点位于垂点N3之后
即计算三角形OFP2中,∠OFP2为angle22,其中令OF=a2,FP2=b2,OP2=c2
a2=F+R
b2=dr(markibeam:end,ibeam)ibeam=1,2,3.......beam
c2=z_sample+R
如上述步骤S323所述,判断上述角度angle2ibeam是否大于θ2/2,如图9所示,以接收的波束线bm4为例,判断波束线上每个样本点是否在有效声场之内。P6,P7为bm4上近场和远场的两个样本点,由图可知P6不在发射声场之内,其与中心发射线的夹角angle21大于θ2/2;样本点P7在发射声场之内,其与中心发射线的夹角angle22小于θ2/2。因此可以通过判断波束线上各点与中心发射线焦点的夹角判断各点是否在有效声场之内。
如上述步骤S324所述,若否,则判定与上述角度angle2ibeam对应的样本点处于有效声场范围内,具体地,判断后对对应的样本点进行标记,具体标记为:
在本发明实施例中,上述根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号的步骤,包括如下步骤:
S221、根据接收线的起始位置坐标M、阵元的位置坐标N、上述校准后的聚焦深度值Z2换算出阵元到聚焦点的距离dn2ibeam;
S222、根据上述校准后的聚焦深度值Z2和阵元到聚焦点的距离dn2ibeam换算出各阵元的延时τ2;
S223、通过各阵元的延时τ2将与之对应的波束线进行相位对齐,并将信号对准后的波束线的BF1信号集进行叠加得到上述BF2信号。
如上述步骤S221所述,根据接收线的起始位置坐标M、阵元的位置坐标N、上述校准后的聚焦深度值Z2换算出阵元到聚焦点的距离dn2ibeam,需要说明的是,阵元到聚焦点的距离dn2ibeam的计算公式通过以下过程获得:
如图10所示,M为接收线的起始位置,N点为阵元的位置,聚焦点Q在波束线上移动,dn为阵元到聚焦点距离,由图可知:
式中,(ref-xe)为阵元到接收线起始位置的距离。
由于阵元到聚焦点的距离为dn
式中,|MP|=Z2ibeam(j),∠NMP=δ
因此可得:
如上述步骤S222所述,根据上述校准后的聚焦深度值Z2和阵元到聚焦点的距离dn2ibeam换算出各阵元的延时τ2,需要说明的是,上述换算过程所采用的公式在本发明实施例中优选为:
τ2j=(Z2ibeam(j)+dnibeam(j))/c 公式(24)
如上述步骤S223所述,通过各阵元的延时τ2将与之对应的波束线进行相位对齐,并将信号对准后的波束线的BF1信号集进行叠加得到上述BF2信号,需要说明的是,在BF1信号进行相位对齐的步骤具体为如下过程:
S223-1、根据各阵元的延时τ2和处于有效声场范围内的上述角度angle2ibeam换算出对应样本点的回波信号,需要说明的是,上述换算过程所采用的公式在本发明实施例中优选为:
式中,S2ibeam(j)为第ibeam条接收波束线上第j个样本点的接收回波信号,s2j(τ2)为子孔径中阵元j接收到的回波。
通过上述公式(25)即可得到相位对齐的多波束回波信号,即,BF1信号。其中,BF1信号的线数与beam有关,发射1次接收beam次,因此发射nlines次,共接收到得到nlines×beam条波束。
需要说明的是,在将信号对准后的BF1信号集进行叠加得到上述BF2信号之前,还包括如下步骤:
S223-2、将上述对应样本点的回波信号汇集形成上述信号对准后的BF1信号集,将nlines×beam条波束的BF1信号,汇集为BF1信号集,该BF1信号集即为本次波束合成过程中所有有效声场中的BF1信号的总集。
当汇集完BF1信号集后对信号集中的BF1信号进行叠加最后形成BF2信号,由于BF1信号的相位已经对齐,所以只需同一空间数据相叠加即可获得BF2信号。
在实际应用中,如图7所示,图中为接收16波束,step=4的情况下的多波束叠加示意图,其他情况,如波束32,step=4等,可以同理得到。
为了解决上述问题,本发明实施例公开了一种超声成像方法,包括如下步骤,
获取超声波数据并转化为相应的数字信号;
将上述数字信号通过指定的校正合成为射频信号;
将上述射频信号通过指定的信号处理分离出载波信号;
将上述载波信号通过指定的图像处理获得超声图像,
上述将上述数字信号通过指定的校正合成为射频信号的步骤,包括上述任一实施例所述的波束合成方法。
对于装置实施例而言,由于其与方法实施例基本相似,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
参照图11,示出了本发明的一种波束合成装置,包括如下具体模块:
校准模块1,用于通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z;
合成模块2,用于根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号。
上述校准模块1,一般用于通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z,需要说明的是,对所述聚焦深度值进行校准的方式在本发明实施例中包括两种,分别为根据线阵校正法和凸阵校正法,其中,线阵校正法通过阵元数N_elements、相邻阵元间距pitch、发射线数nlines、移动步长step、每次发射或接收的波束数量beam、发射线上样本点到发射点的距离z_sample1和焦距F对所述聚焦深度值Z进行校正;凸阵校正法通过根据阵元数N_elements、相邻阵元间距pitch、凸阵圆心半径R、发射线数nlines、移动步长step、每次发射或接收的波束数量beam、焦距F、发射线上样本点到圆心的距离z_sample2对所述聚焦深度值Z进行校正。
上述合成模块2,一般用于根据上述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,再将上述BF1信号集中的各BF1信号进行指定叠加得到BF2信号,需要说明的是,步骤S2的实现方法包括两种,分别为根据线阵校准法和凸阵校准法,其中,线阵校准法通过接收线的坐标ref1、阵元坐标xe1、校准后的聚焦深度值Z2将对应的波束线进行相位对齐,并将信号对准后的BF1信号集进行叠加得到所述BF2信号;凸阵校准发通过接收线的起始位置坐标M、阵元的位置坐标N、所述校准后的聚焦深度值Z2将对应的波束线进行相位对齐,并将信号对准后的BF1信号集进行叠加得到所述BF2信号。在通过校准后的聚焦深度值Z2将对应的波束线进行相位对齐前,波束线已经空间对准完毕,通过上述校准后的聚焦深度值Z2对波束进行延时。得到BF1,判断数据的有效性,再叠加得到BF2由于BF1信号的相位已经对齐,所以只需同一空间数据相叠加。
示出了本发明的一种超声成像装置,包括如下具体模块:
获取模块,用于获取超声波数据并转化为相应的数字信号;
波束合成模块,用于将上述数字信号通过指定的校正合成为射频信号;
分离模块,用于将上述射频信号通过指定的信号处理分离出载波信号;
成像模块,用于将上述载波信号通过指定的图像处理获得超声图像;
上述波束合成模块,包括上述任一实施例所述的波束合成装置。
示出了本发明的一种超声成像设备,包括产生超声波信号辐射到被测试组织中和吸收反射声波信号的信号收发器,将接收的声波转换为电信号的振动传感器,将接收信号采样和数字化的模数转换器(A/D),补偿由于深度导致的超声波幅值衰减的时间增益补偿器(TGC),将数字信号转化为RF信号的接收波束合成器,将RF信号进行包络提取和解调处理分离出载波信号的信号处理器,将载波信号进行扫描转换和后端图像处理的数字扫描转换器DSC以及用于最终显示的图像的显示器。
其中,上述接收波束合成器,包括相位校正模块、用于合成BF1信号的第一次波束波束合成模块、有效声场判断模块和用于合成BF2信号的第二次波束合成器,且上述相位校正模块、用于合成BF1信号的第一次波束波束合成模块、有效声场判断模块和用于合成BF2信号的第二次波束合成器依次相连。
参照图12,示出了本发明的一种实现所述波束合成方法的计算机设备,具体可以包括如下:
上述计算机设备12以通用计算设备的形式表现,计算机设备12的组件可以包括但不限于:一个或者多个处理器或者处理单元16,系统存储器28,连接不同系统组件(包括系统存储器28和处理单元16)的总线18。
总线18表示几类总线18结构中的一种或多种,包括存储器总线18或者存储器控制器,外围总线18,图形加速端口,处理器或者使用多种总线18结构中的任意总线18结构的局域总线18。举例来说,这些体系结构包括但不限于工业标准体系结构(ISA)总线18,微通道体系结构(MAC)总线18,增强型ISA总线18、音视频电子标准协会(VESA)局域总线18以及外围组件互连(PCI)总线18。
计算机设备12典型地包括多种计算机系统可读介质。这些介质可以是任何能够被计算机设备12访问的可用介质,包括易失性和非易失性介质,可移动的和不可移动的介质。
系统存储器28可以包括易失性存储器形式的计算机系统可读介质,例如随机存取存储器(RAM)30和/或高速缓存存储器32。计算机设备12可以进一步包括其他移动/不可移动的、易失性/非易失性计算机体统存储介质。仅作为举例,存储系统34可以用于读写不可移动的、非易失性磁介质(通常称为“硬盘驱动器”)。尽管图12中未示出,可以提供用于对可移动非易失性磁盘(如“软盘”)读写的磁盘驱动器,以及对可移动非易失性光盘(例如CD~ROM,DVD~ROM或者其他光介质)读写的光盘驱动器。在这些情况下,每个驱动器可以通过一个或者多个数据介质界面与总线18相连。存储器可以包括至少一个程序产品,该程序产品具有一组(例如至少一个)程序模块42,这些程序模块42被配置以执行本发明各实施例的功能。
具有一组(至少一个)程序模块42的程序/实用工具40,可以存储在例如存储器中,这样的程序模块42包括——但不限于——操作系统、一个或者多个应用程序、其他程序模块42以及程序数据,这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块42通常执行本发明所描述的实施例中的功能和/或方法。
计算机设备12也可以与一个或多个外部设备14(例如键盘、指向设备、显示器24、摄像头等)通信,还可与一个或者多个使得用户能与该计算机设备12交互的设备通信,和/或与使得该计算机设备12能与一个或多个其他计算设备进行通信的任何设备(例如网卡,调制解调器等等)通信。这种通信可以通过输入/输出(I/O)界面22进行。并且,计算机设备12还可以通过网络适配器20与一个或者多个网络(例如局域网(LAN)),广域网(WAN)和/或公共网络(例如因特网)通信。如图所示,网络适配器20通过总线18与计算机设备12的其他模块通信。应当明白,尽管图12中未示出,可以结合计算机设备12使用其他硬件和/或软件模块,包括但不限于:微代码、设备驱动器、冗余处理单元16、外部磁盘驱动阵列、RAID系统、磁带驱动器以及数据备份存储系统34等。
处理单元16通过运行存储在系统存储器28中的程序,从而执行各种功能应用以及数据处理,例如实现本发明实施例所提供的波束合成方法。
也即,上述处理单元16执行上述程序时实现:通过对原始乳腺图像进行显著性检测,调整上述原始乳腺图像的目标区域和非目标区域之间的对比度,获得显著性图像;对上述显著性图像进行乳腺轮廓分割处理,确定上述显著性图像中的有效乳腺区域;将确定有效乳腺区域后的上述显著性图像进行区域粗分割,保留满足指定灰度阈值的粗分割区域,获得粗分割图像;采用K-means聚类方法对上述粗分割图像中的指定粗分割区域进行细分割,获得细分割图像;根据指定形态学特征对上述细分割图像进行假阳性区域滤除,获得真阳性区域图像;将真阳性区域图像中的真阳性区域标记在上述对显著性图像进行乳腺轮廓分割处理后的图像中,获得标记出了病变区域的乳腺图像。
在上述任一发明实施例中,通过实现多波束的空间对准和相位对齐,能够较好地抑制旁瓣,提高图像分辨率及信噪比;通过多波束叠加极大的提高线密度,提高了图像的空间分辨率;本方法实现算法简单,在消耗少量处理资源的情况下即可实现多帧率,高信噪比,高分辨率的效果;装置通过在波束合成模块中增加了校准模块和合成模块,使波束合成后更数据更精确;并且通过声场判断,对有效声场数据叠加,相较于传统方法实现了较好的聚焦效果,提高了图像质量。
以上对本申请所提供的一种波束合成方法及装置,进行了详细介绍,本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本申请的方法及其核心思想;同时,对于本领域的一般技术人员,依据本申请的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本申请的限制。
Claims (6)
1.一种波束合成方法,其特征在于,包括如下步骤:
通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z;具体地,根据阵元数N_elements、相邻阵元间距pitch和凸阵圆心半径R换算出凸阵圆心角beta;根据发射线数nlines、移动步长step和所述凸阵圆心角beta换算出接收线间角dx_beta;根据每次发射或接收的波束数量beam和所述接收线间角dx_beta换算出接收波束到发射波束的夹角Angle_beam;根据焦距F、所述凸阵圆心半径R和接收波束到发射波束的夹角Angle_beam分别换算出样本点到接收波束虚拟焦点的横向距离Hxibeam和纵向距离Hyibeam,并根据所述所述凸阵圆心半径R和纵向距离Hyibeam换算出校正焦距FA;根据发射线上样本点到圆心的距离z_sample2、所述凸阵圆心半径R、横向距离Hxibeam和纵向距离Hyibeam换算出中心发射线上的样本点到各波束线虚拟焦点的距离dr2ibeam;根据所述中心发射线上的样本点到各波束线虚拟焦点的距离dr2ibeam、焦距F和校正焦距FA换算出各波束线对应的校准后的聚焦深度值Z2;
根据所述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,根据凸阵圆心半径R、通道数N_channel、焦距F和相邻阵元间距pitch换算出发射声场的角度θ2;获取各波束线上样本点偏离发射线的角度angle2ibeam;判断所述角度angle2ibeam是否大于θ2/2;若否,则判定与所述角度angle2ibeam对应的样本点处于有效声场范围内;根据接收线的起始位置坐标M、阵元的位置坐标N、所述校准后的聚焦深度值Z2换算出阵元到聚焦点的距离dn2ibeam;根据所述校准后的聚焦深度值Z2和阵元到聚焦点的距离dn2ibeam换算出各阵元的延时τ2;通过各阵元的延时τ2将与之对应的BF1信号进行相位对齐,并将信号对准后的BF1信号集进行叠加得到BF2信号。
2.根据权利要求1所述的方法,其特征在于,所述通过各阵元的延时τ2将与之对应的BF1信号进行相位对齐,并将信号对准后的BF1信号集进行叠加得到BF2信号的步骤,包括如下步骤:
根据各阵元的延时τ2和处于有效声场范围内的所述角度angle2ibeam换算出对应样本点的回波信号;
将所述对应样本点的回波信号汇集形成所述信号对准后的BF1信号集。
3.一种超声成像方法,包括如下步骤,
获取超声波数据并转化为相应的数字信号;
将所述数字信号通过指定的校正合成为射频信号;
将所述射频信号通过指定的信号处理分离出载波信号;
将所述载波信号通过指定的图像处理获得超声图像,
其特征在于,所述将所述数字信号通过指定的校正合成为射频信号的步骤,包括权利要求1-2任意一项所述的波束合成方法。
4.一种波束合成装置,其特征在于,包括如下具体模块:
校准模块,用于通过指定步骤校准每条波束对应的聚焦深度值Z,并将校准后的聚焦深度值Z2替换相对应波束原有的聚焦深度值Z;具体地,根据阵元数N_elements、相邻阵元间距pitch和凸阵圆心半径R换算出凸阵圆心角beta;根据发射线数nlines、移动步长step和所述凸阵圆心角beta换算出接收线间角dx_beta;根据每次发射或接收的波束数量beam和所述接收线间角dx_beta换算出接收波束到发射波束的夹角Angle_beam;根据焦距F、所述凸阵圆心半径R和接收波束到发射波束的夹角Angle_beam分别换算出样本点到接收波束虚拟焦点的横向距离Hxibeam和纵向距离Hyibeam,并根据所述所述凸阵圆心半径R和纵向距离Hyibeam换算出校正焦距FA;根据发射线上样本点到圆心的距离z_sample2、所述凸阵圆心半径R、横向距离Hxibeam和纵向距离Hyibeam换算出中心发射线上的样本点到各波束线虚拟焦点的距离dr2ibeam;根据所述中心发射线上的样本点到各波束线虚拟焦点的距离dr2ibeam、焦距F和校正焦距FA换算出各波束线对应的校准后的聚焦深度值Z2;
合成模块,用于根据所述校准后的聚焦深度值Z2进行信号对准,并获取信号对准后的波束集中的BF1信号集,根据凸阵圆心半径R、通道数N_channel、焦距F和相邻阵元间距pitch换算出发射声场的角度θ2;获取各波束线上样本点偏离发射线的角度angle2ibeam;判断所述角度angle2ibeam是否大于θ2/2;若否,则判定与所述角度angle2ibeam对应的样本点处于有效声场范围内;根据接收线的起始位置坐标M、阵元的位置坐标N、所述校准后的聚焦深度值Z2换算出阵元到聚焦点的距离dn2ibeam;根据所述校准后的聚焦深度值Z2和阵元到聚焦点的距离dn2ibeam换算出各阵元的延时τ2;通过各阵元的延时τ2将与之对应的BF1信号进行相位对齐,并将信号对准后的BF1信号集进行叠加得到BF2信号。
5.一种超声成像装置,包括如下具体模块:
获取模块,用于获取超声波数据并转化为相应的数字信号;
波束合成模块,用于将所述数字信号通过指定的校正合成为射频信号;
分离模块,用于将所述射频信号通过指定的信号处理分离出载波信号;
成像模块,用于将所述载波信号通过指定的图像处理获得超声图像;
其特征在于,所述波束合成模块,包括权利要求4所述的波束合成装置。
6.一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1~2中任意一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111470011.2A CN114129185B (zh) | 2019-03-18 | 2019-03-18 | 波束合成方法、超声成像方法、装置及设备 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910204036.4A CN110101409B (zh) | 2019-03-18 | 2019-03-18 | 波束合成方法、超声成像方法、装置及设备 |
CN202111470011.2A CN114129185B (zh) | 2019-03-18 | 2019-03-18 | 波束合成方法、超声成像方法、装置及设备 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910204036.4A Division CN110101409B (zh) | 2019-03-18 | 2019-03-18 | 波束合成方法、超声成像方法、装置及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114129185A CN114129185A (zh) | 2022-03-04 |
CN114129185B true CN114129185B (zh) | 2023-11-07 |
Family
ID=67484350
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111470011.2A Active CN114129185B (zh) | 2019-03-18 | 2019-03-18 | 波束合成方法、超声成像方法、装置及设备 |
CN201910204036.4A Active CN110101409B (zh) | 2019-03-18 | 2019-03-18 | 波束合成方法、超声成像方法、装置及设备 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910204036.4A Active CN110101409B (zh) | 2019-03-18 | 2019-03-18 | 波束合成方法、超声成像方法、装置及设备 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN114129185B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110974304B (zh) * | 2019-12-13 | 2021-07-30 | 山东大学齐鲁医院 | 基于穿戴式柔性超声换能器的超声波束合成系统及方法 |
US20240111046A1 (en) * | 2019-11-12 | 2024-04-04 | GE Precision Healthcare LLC | Method and system for flow processing on channel data for application of nonlinear beamforming |
WO2021114071A1 (zh) * | 2019-12-10 | 2021-06-17 | 深圳大学 | 一种经颅超声成像方法、装置及计算机可读存储介质 |
CN111772679B (zh) * | 2020-08-04 | 2023-03-31 | 深圳蓝影医学科技股份有限公司 | 一种超声成像的方法、装置、计算机设备及计算机可读介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6432056B1 (en) * | 1999-10-08 | 2002-08-13 | Atl Ultrasound | Ultrasonic diagnostic imaging system with high frame rate synthetic transmit focus |
CN1788685A (zh) * | 2004-12-15 | 2006-06-21 | 深圳迈瑞生物医疗电子股份有限公司 | 基于双波束及合成孔径的接收方法及其装置 |
CN102727256A (zh) * | 2012-07-23 | 2012-10-17 | 重庆博恩富克医疗设备有限公司 | 基于虚拟阵元的双聚焦波束合成方法及装置 |
CN102895000A (zh) * | 2012-11-06 | 2013-01-30 | 重庆大学 | 一种基于自适应加权的双聚焦波束合成方法 |
CN105997137A (zh) * | 2015-03-27 | 2016-10-12 | 柯尼卡美能达株式会社 | 超声波信号处理装置、以及超声波诊断装置 |
CN109363714A (zh) * | 2018-09-17 | 2019-02-22 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像设备及其超声成像方法 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4989143A (en) * | 1987-12-11 | 1991-01-29 | General Electric Company | Adaptive coherent energy beam formation using iterative phase conjugation |
CN101097256B (zh) * | 2006-06-28 | 2011-03-23 | 深圳迈瑞生物医疗电子股份有限公司 | 超声波束合成聚焦参数实时计算方法及装置 |
CN101664321B (zh) * | 2009-09-07 | 2011-06-15 | 无锡祥生科技有限公司 | 组织声速实时可调的超声诊断设备及其波束合成方法 |
CN101865998A (zh) * | 2010-05-19 | 2010-10-20 | 深圳市蓝韵实业有限公司 | 一种延时参数实时计算和流水线的多波束合成方法 |
CN101915916B (zh) * | 2010-08-04 | 2013-01-30 | 深圳市蓝韵实业有限公司 | 一种超声发射延时处理的方法 |
CN102846336B (zh) * | 2011-06-27 | 2014-10-22 | 深圳市蓝韵实业有限公司 | 超声成像延时控制方法及装置 |
CN102670250B (zh) * | 2012-05-24 | 2014-07-09 | 深圳市开立科技有限公司 | 超声波换能器的阵列校准系统及方法 |
CN103654732B (zh) * | 2013-12-31 | 2015-08-26 | 南京大学 | 一种基于线性延时补偿的光声图像优化方法 |
US9763646B2 (en) * | 2014-06-12 | 2017-09-19 | General Electric Company | Method and systems for adjusting a pulse generated for ultrasound multi-line transmit |
WO2016125509A1 (ja) * | 2015-02-04 | 2016-08-11 | 日立アロカメディカル株式会社 | 超音波撮像装置、および、超音波信号の処理方法 |
CN106971055B (zh) * | 2016-01-12 | 2019-10-25 | 深圳迈瑞生物医疗电子股份有限公司 | 超声造影成像方法和系统 |
CN105832366B (zh) * | 2016-03-18 | 2019-11-15 | 深圳安盛生物医疗技术有限公司 | 一种用于波束合成过程中的延时实时计算方法 |
JP6747108B2 (ja) * | 2016-07-05 | 2020-08-26 | コニカミノルタ株式会社 | 超音波信号処理装置、超音波信号処理方法、及び、超音波診断装置 |
-
2019
- 2019-03-18 CN CN202111470011.2A patent/CN114129185B/zh active Active
- 2019-03-18 CN CN201910204036.4A patent/CN110101409B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6432056B1 (en) * | 1999-10-08 | 2002-08-13 | Atl Ultrasound | Ultrasonic diagnostic imaging system with high frame rate synthetic transmit focus |
CN1788685A (zh) * | 2004-12-15 | 2006-06-21 | 深圳迈瑞生物医疗电子股份有限公司 | 基于双波束及合成孔径的接收方法及其装置 |
CN102727256A (zh) * | 2012-07-23 | 2012-10-17 | 重庆博恩富克医疗设备有限公司 | 基于虚拟阵元的双聚焦波束合成方法及装置 |
CN102895000A (zh) * | 2012-11-06 | 2013-01-30 | 重庆大学 | 一种基于自适应加权的双聚焦波束合成方法 |
CN105997137A (zh) * | 2015-03-27 | 2016-10-12 | 柯尼卡美能达株式会社 | 超声波信号处理装置、以及超声波诊断装置 |
CN109363714A (zh) * | 2018-09-17 | 2019-02-22 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像设备及其超声成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110101409B (zh) | 2022-03-22 |
CN110101409A (zh) | 2019-08-09 |
CN114129185A (zh) | 2022-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114129185B (zh) | 波束合成方法、超声成像方法、装置及设备 | |
US7740583B2 (en) | Time delay estimation method and system for use in ultrasound imaging | |
US8968199B2 (en) | Spatial compound imaging in an ultrasound system | |
US8068647B2 (en) | Method and apparatus for real-time motion correction for ultrasound spatial compound imaging | |
US20070167802A1 (en) | Accurate time delay estimation method and system for use in ultrasound imaging | |
JP5394694B2 (ja) | サイドローブの影響を除去する方法 | |
CN110101411B (zh) | 超声成像空间复合方法及系统 | |
US20130258805A1 (en) | Methods and systems for producing compounded ultrasound images | |
CN107789008B (zh) | 一种基于通道数据的自适应超声波束合成方法和系统 | |
JP2006204923A (ja) | コヒーレンス度適応超音波撮像 | |
JPH10506800A (ja) | 超音波画像における調節可能な周波数走査方法および装置 | |
KR20080039446A (ko) | 실시간 공간 합성을 이용하는 흐름 영상처리를 위한 초음파영상처리 시스템 및 방법 | |
JP5069022B2 (ja) | 超音波撮像で使用するための正確な時間遅延推定の方法及びシステム | |
CN102727255A (zh) | 一种超声图像空间复合成像的方法及装置 | |
WO2019127621A1 (zh) | 超声成像方法、系统和设备 | |
US5249577A (en) | Ultrasonotomography | |
CN107997783B (zh) | 一种基于超声波方向性的自适应超声波束合成方法和系统 | |
US20040054285A1 (en) | Elevation compounding for ultrasound imaging | |
KR101348770B1 (ko) | 도플러 스펙트럼 영상의 합성 영상을 제공하는 초음파 시스템 및 방법 | |
US20180038947A1 (en) | Clutter suppression in ultrasonic imaging systems | |
US20110245676A1 (en) | Method and apparatus for ultrasound signal acquisition and processing | |
CN110731795A (zh) | 空间复合成像的处理方法和装置 | |
KR101388333B1 (ko) | 빔형성 시간 지연 교정 방법과 빔형성 시간 지연을추정하는 초음파 시스템 | |
EP2466330B1 (en) | Ultrasound system and method for processing beam-forming based on sampling data | |
KR20090042184A (ko) | 영상 움직임 추정 및 보상 장치와 그 방법 |
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 |
Address after: No.103, baguang District Service Center, No.2 BaiShaWan Road, baguang community, Kuiyong street, Dapeng New District, Shenzhen, Guangdong 518000 Applicant after: Shenzhen Lanying Medical Technology Co.,Ltd. Address before: 518000 baguang area service center 103, No. 2 BaiShaWan Road, baguang community, Kuiyong street, Dapeng new area, Longgang District, Shenzhen, Guangdong Province Applicant before: SHENZHEN LANYUN MEDICAL IMAGE CO.,LTD. |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |