CN105044706B - 一种自适应波束形成方法 - Google Patents
一种自适应波束形成方法 Download PDFInfo
- Publication number
- CN105044706B CN105044706B CN201510342472.XA CN201510342472A CN105044706B CN 105044706 B CN105044706 B CN 105044706B CN 201510342472 A CN201510342472 A CN 201510342472A CN 105044706 B CN105044706 B CN 105044706B
- Authority
- CN
- China
- Prior art keywords
- array
- array element
- covariance matrix
- wave beam
- adaptive beamformer
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 35
- 238000003384 imaging method Methods 0.000 claims abstract description 17
- 238000009499 grossing Methods 0.000 claims description 18
- 230000008707 rearrangement Effects 0.000 claims description 6
- 238000011068 loading method Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000002604 ultrasonography Methods 0.000 description 3
- 238000007792 addition Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明涉及一种自适应波束形成方法。在一个实施例中，该方法包括以下步骤：根据发射波束的空间位置重新构造换能器阵元的排序,将靠近波束的阵元置于新阵列的中心、远离波束的阵元置于新阵列的两端；根据重新排列的阵元顺序计算用于最小方差波束形成的协方差矩阵；根据估计得到的所述协方差矩阵计算自适应波束形成的输出结果。本发明实施例通过重新构造换能器阵元排序，使阵列两端的回波信息得到了充分的利用，改善了图像边缘区域和中心区域的分辨率的不一致性，提高了超声成像的质量。
Description
技术领域
本发明涉及超声成像中的波束形成方法，尤其涉及一种基于改进空间平滑方式的自适应波束形成方法。
背景技术
传统医学超声诊断及超声无损检测系统普遍使用延时叠加的波束形成方法，超声成像的质量较差。
最小方差自适应波束形成方法虽然能提高超声成像的质量，但最小方差波束形成方法的核心问题之一协方差矩阵的估计，因超声回波信号具有较强的相干性，需要进行空间平滑的方式处理。传统空间平滑过程，是将整个阵列划分为若干子阵列，每个子阵列包含有相同数目的阵元，并从阵列的一端依次向另一端滑动构造出一组子阵列，然后将这些子阵列的信号进行叠加并取平均，最终再结合时间平滑及对角加载等技术得到用于最小方差波束形成方法的协方差矩阵。
这种空间平滑方式没有充分利用所有阵元所接收到的回波信息，位于原阵列中心位置的阵元得到更多的加成，阵列两端的阵元的有效信息却没有得以充分的利用。导致在B模式成像下，成像结果中图像边缘部分的分辨率相较中心位置有较大损失，即使得整个图像的分辨率由于位置关系而变得不一致，且这种现象随着换能器孔径的增大而愈发明显。即使近年来引入超声成像中的基于前后向空间平滑的自适应波束形成方法能够在一定程度上提高成像的质量，但也主要反映在改善了图像的对比度，而在分辨率方面依然存在上述问题。
因此，现有的空间平滑方式下的自适应波束形成方法对换能器阵列两端的阵元信息利用率低，导致在分辨率，尤其是点散射目标的辨识程度上存在缺陷。急需对现有空间平滑方式加以改进，以全面提高超声图像的整体质量。
发明内容
本发明的目的在于解决现有的空间平滑方式下的自适应波束形成方法存在的上述缺陷。
为实现上述目的，本发明提供了一种自适应波束形成方法，具体包括以下步骤：
根据发射波束的空间位置重新构造换能器阵元的排序,将靠近波束的阵元置于新阵列的中心、远离波束的阵元置于新阵列的两端；根据重新排列的阵元顺序计算用于最小方差波束形成的协方差矩阵；根据估计得到的所述协方差矩阵计算自适应波束形成的输出结果。
进一步的，所述阵元的排序包括：
假设在B模式超声成像下，生成一帧图像需发射M-1个波束，即每条波束均匀的分布在各个阵元之间，则每一根扫描线对应的所有阵元接收的回波数据可以用M×N的矩阵来描述，其中M为阵元数目，N为采样数目，经过时延预处理后每个阵元接收到的信号可以用向量表示为其中l＝1，2，…，M，k＝1，2，…，N。
根据每次发射的波束与阵元的位置关系重新构造阵列的阵元排序，让靠近波束的阵元位于新阵列的中心，与此对应，远离发射波束的阵元位于新阵列的两端。从而得到回波数据以用于最小方差波束形成方法中的协方差矩阵。
进一步的，所述协方差矩阵的计算包括：
根据重新构造的阵列进行空间平滑与时间平滑处理，即将新阵列以L个阵元为一组划分成M-L+1个子阵列，其中，每一个子阵列由L个重新排序的阵元组成，由新阵元序号1开始滑动至M，按下式进行计算，得到协方差矩阵R(k)：
根据对角加载方法对所述协方差矩阵进行处理，以提高协方差矩阵的稳健性：
R(k)＝R(k)+εI
其中，I为单位阵，以及
加载因子ε的计算公式如下：
ε＝Δ·tr(R(k))
其中，Δ一般介于1/10L与1/1000L之间；
进一步的，根据估计得到的所述协方差矩阵计算自适应波束形成的输出结果，包括：
根据所述协方差矩阵，利用最小方差波束形成法计算出一组最优的权向量w(k)，计算公式如下：
其中，R-1(k)是所述协方差矩阵中所得到的协方差矩阵的逆，导向矢量a中的元素均取为“1”；
根据上述方法得到的所述一组最优的权向量计算B模式成像下一根扫描
线对应的波束形成的输出结果，由此进行图像变换最终显示成像：
y(k)＝wH(k)·xd(k)
其中，xd是所述新阵列的回波数据。
本发明的技术优势在于，根据发射波束与阵元的空间位置关系来构造用于空间平滑的阵列，使阵列边缘的阵元回波信息得到多叠加，充分利用了换能器各个阵元所接收到的回波信息，改善了超声成像在横向分辨率方面的不一致性，尤其能提高对点散射目标的辨识程度，提高了超声成像的质量。
附图说明
图1为本发明实施例提供的一种自适应波束形成方法流程示意图；
图2为本发明实施例发射波束与阵元位置关系的示意图；
图3为本发明实施例使用的空间平滑方法的子阵列划分示意图；
图4为本发明实施例对一个由8个阵元组成的阵列在第1次发射波束时所构造的新阵列示意图；
图5为本发明实施例对比使用其他几种空间平滑方式下的自适应波束形成方法的B模式成像仿真图；
图6为图5所示实施例在60mm深度截面的波束图。
具体实施方式
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
图1为本发明实施例提供的一种自适应波束形成方法流程示意图。如图1所示，本发明实施例应用于B模式下的超声成像。等间距的M个阵元的超声换能器阵列工作在B模式成像线性扫描方式下，对所接收到的回波数据进行波束形成处理，每一根扫描线对应着M个阵元的接收信号，并且每个阵元的采样数据长度为N，这些时域上的回波数据已经经过动态聚焦处理，即在每一采样深度都将回波数据进行“对齐”，得到的xd作为下述过程中的默认处理对象。本发明实施例的波束形成处理方法包括步骤101-103：
步骤101，根据发射波束与换能器阵元的空间位置关系对阵元进行重新排序。如图2所示，阵元1和阵元2距离第1条发射波束最近，阵元3、4、…、M依次远离之。如图3所示，利用等距线阵的平移不变性,将阵元均匀划分为相互重叠的M-L+1个子阵列，每个子阵列都包含了数量相同的重新排序的阵元。每一次空间平滑均从新阵元序号1开始滑动至M。
图4为本发明实施例对一个由8个阵元组成的阵列在第1次发射波束(即阵列的最左侧)时所构造的新阵列示意图。如图4所示，对M＝8的阵列，子阵阵列元数目L＝4，对应在第1条发射波束时重新排列的阵元顺序，这种排序方式的原则在于使靠近发射波束的阵元位于新阵列的中心、远离的则位于新阵列的两端。位于阵列中靠近波束发射位置的阵元，经过重新排序后，始终位于新阵列的中心。如图4所示，经过换能器阵元的重新排序，新阵列中阵元依次为：7、5、3、1、2、4，6、8(1至8分别为排序前从左到右的阵列号)。空间平滑以子阵为单位，每个子阵由数量相同的阵元组成，空间平滑时由新阵元序号1(即原第7号)开始滑动至8。即从排序前阵列的阵元7开始依次经过阵元5、3、1、2、4、6、8，第1子阵的四个阵元依次为排序前的阵元7、5、3、1，第2子阵的四个阵元依次为排序前的阵元5、3、1、2…第5子阵的四个阵元依次为排序前的阵元2、4、6、8。
步骤102，根据重新排列的阵元顺序计算用于最小方差波束形成的协方差矩阵，所述步骤102的具体处理方法包括步骤1021-1022：
步骤1021，根据由步骤101中重新构造的阵列进行空间平滑与时间平滑处理，即将新阵列以L个阵元为一组划分成M-L+1个子阵列，其中，每一个子阵列由L个重新排序的阵元组成，由新阵元序号1开始滑动至M，协方差矩阵计算的公式如下：
从公式看，不同子阵的输出向量中多次涉及阵元1及其附近的阵元，而较少涉及边缘的阵元7和8，因此协方差矩阵的计算中对阵元1及其附近的阵元做了更多的加权。
如前文所述，这里的阵元1是根据发射波束与换能器阵元的空间位置关系重新排序的，实际对应于发射波束的阵元。因此，在本发明的实施例中，通过发射波束与换能器阵元的空间位置关系重新排序，协方差矩阵的计算中对发射波束附近的阵元做了相似的加权，无论该发射波束附近的阵元是位于阵列中心还是阵列的边缘。
步骤1022，为提高协方差矩阵的稳健性，使用对角加载技术进行如下的处理：
R(k)＝R(k)+εI
其中，I为单位阵，加载因子ε的计算公式如下：
ε＝Δ·tr(R(k))
参数Δ一般介于1/10L与1/1000L之间；
步骤103，根据得到的协方差矩阵计算自适应波束形成的输出结果，所述步骤103的具体处理方法包括步骤1031-1032：
步骤1031，根据由步骤102中所得到的协方差矩阵，利用最小方差波束形成的基本原理计算出一组最优的权向量，计算公式如下：
其中R-1(k)是由步骤1022中所得到的协方差矩阵的逆，导向矢量a中的元素均取为“1”；
步骤1032，根据步骤1031中所得到的权向量计算B模式成像下一根扫描线对应的波束形成的输出：
y(k)＝wH(k)·xd(k)
y(k)即为用于最终的显示成像的处理结果。
图5为本发明实施例对比使用其他几种空间平滑方式下的自适应波束形成方法的B模式成像仿真图。如图5所示，使用Field II软件进行仿真模拟，并将本发明(BSAMV)对比了其他几种常用的波束形成方法，即延时叠加法(DAS)、标准的最小方差波束形成方法(SMV)以及基于前后向空间平滑的最小方差波束形成方法(FBMV)。其中M＝128，L＝64，K＝2，Δ＝1/100L。
图6为图5所示实施例在60mm深度截面的波束图。如图6所示，从仿真的结果中可以看出本发明能够改善成像的横向分辨率，缩小了位于图像边缘区域和中心区域的目标的分辨率方面的差距，同时对于点散射目标具有更好的辨识能力。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (3)
1.一种自适应波束形成方法，其特征在于，所述方法包括：
根据发射波束和换能器阵元的空间位置重新构造换能器阵元的排序,将靠近波束的阵元置于新阵列的中心、远离波束的阵元置于新阵列的两端；
根据重新排列的阵元顺序计算用于最小方差波束形成的协方差矩阵；
所述协方差矩阵的计算中对发射波束附近的阵元做了相似的加权；所述发射波束附近的阵元位于阵列中心或边缘；
根据估计得到的所述协方差矩阵计算自适应波束形成的输出结果。
2.根据权利要求1中所述的方法，其特征在于，所述根据重新排列的阵元顺序计算用于最小方差波束形成的协方差矩阵，包括：
根据重新构造的阵列进行空间平滑与时间平滑处理，即将新阵列以L个阵元为一组划分成M-L+1个子阵列，其中，每一个子阵列由L个重新排序的阵元组成，由新阵元序号1开始滑动至M，按下式进行计算，得到协方差矩阵
其中，K是选取的时间平滑参数，p是新阵列中子阵列序号，是的共轭转置，是第p个子阵的输出向量；
根据对角加载方法对所述协方差矩阵进行处理，得到对角加载后的协方差矩阵：
其中，I为单位阵；
按下式对加载因子ε进行计算：
其中，Δ一般介于1/10L与1/1000L之间。
3.根据权利要求1中所述的方法，其特征在于，所述根据估计得到的所述协方差矩阵计算自适应波束形成的输出结果，包括：
根据估计得到的所述协方差矩阵，利用最小方差波束形成法计算出一组最优的权向量w(k)：
其中，是所述协方差矩阵中所得到的协方差矩阵的逆，导向矢量a中的元素均取为“1”；
根据上述方法得到的所述一组最优的权向量计算B模式成像下一根扫描线对应的波束形成的输出结果，由此进行图像变换最终显示成像：
y(k)＝wH(k)·xd(k)
其中，xd是所述新阵列的回波数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510342472.XA CN105044706B (zh) | 2015-06-18 | 2015-06-18 | 一种自适应波束形成方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510342472.XA CN105044706B (zh) | 2015-06-18 | 2015-06-18 | 一种自适应波束形成方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105044706A CN105044706A (zh) | 2015-11-11 |
| CN105044706B true CN105044706B (zh) | 2018-06-29 |
Family
ID=54451382
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510342472.XA Expired - Fee Related CN105044706B (zh) | 2015-06-18 | 2015-06-18 | 一种自适应波束形成方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105044706B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3509325B1 (en) * | 2016-05-30 | 2021-01-27 | Oticon A/s | A hearing aid comprising a beam former filtering unit comprising a smoothing unit |
| CN107544059A (zh) * | 2017-07-20 | 2018-01-05 | 天津大学 | 一种基于对角加载技术的稳健自适应波束形成方法 |
| CN111208213A (zh) * | 2020-02-25 | 2020-05-29 | 重庆大学 | 融合交替乘子迭代的谱寻求子带最小方差超声成像算法 |
| CN116421218A (zh) * | 2021-10-22 | 2023-07-14 | 武汉联影医疗科技有限公司 | 一种超声波发射方法和系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101238992A (zh) * | 2008-02-01 | 2008-08-13 | 哈尔滨工业大学 | 基于相关性分析的超声成像系统自适应波束形成器 |
| CN102499712A (zh) * | 2011-09-30 | 2012-06-20 | 重庆大学 | 一种基于特征空间的前后向自适应波束形成方法 |
| CN102895000A (zh) * | 2012-11-06 | 2013-01-30 | 重庆大学 | 一种基于自适应加权的双聚焦波束合成方法 |
-
2015
- 2015-06-18 CN CN201510342472.XA patent/CN105044706B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101238992A (zh) * | 2008-02-01 | 2008-08-13 | 哈尔滨工业大学 | 基于相关性分析的超声成像系统自适应波束形成器 |
| CN102499712A (zh) * | 2011-09-30 | 2012-06-20 | 重庆大学 | 一种基于特征空间的前后向自适应波束形成方法 |
| CN102895000A (zh) * | 2012-11-06 | 2013-01-30 | 重庆大学 | 一种基于自适应加权的双聚焦波束合成方法 |
Non-Patent Citations (3)
| Title |
|---|
| 一种改进的Capon波束形成算法在超声成像中的研究;杜春宁;《中国生物医学工程学报》;20071031;第26卷(第5期);全文 * |
| 利用回波相关性的医学超声成像相位畸变校正方法;夏新源;《声学学报》;20150331;第40卷(第2期);全文 * |
| 最小方差波束形成与广义相干系数融合的医学超声成像方法;吴文焘;《声学学报》;20110131;第36卷(第1期);第67页右栏-68页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105044706A (zh) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105044706B (zh) | 一种自适应波束形成方法 | |
| Kortbek et al. | Sequential beamforming for synthetic aperture imaging | |
| JP5587743B2 (ja) | 超音波画像処理装置 | |
| CN101683274B (zh) | 偏转凸阵数字扫描变换和凸阵空间复合的方法及成像系统 | |
| CN106097427A (zh) | 子容积的连续定向的增强超声成像 | |
| CN105074498B (zh) | 具有对梯形的扇区的空间复合的超声诊断成像系统 | |
| JPH04244145A (ja) | ソナー映像変換方法 | |
| Senouf et al. | High frame-rate cardiac ultrasound imaging with deep learning | |
| JP5600285B2 (ja) | 超音波画像処理装置 | |
| CN109363714A (zh) | 一种超声成像设备及其超声成像方法 | |
| US20170269208A1 (en) | Ultrasound transducer arrays with variable patch geometries | |
| Zeng et al. | Correspondence-beam-domain eigenspace-based minimum variance beamformer for medical ultrasound imaging | |
| CN108107428B (zh) | 用于mimo阵列的相移偏移成像方法及装置 | |
| CN104902825A (zh) | 超声波拍摄装置 | |
| Rabinovich et al. | Multi-line acquisition with minimum variance beamforming in medical ultrasound imaging | |
| Cruza et al. | A new beamforming method and hardware architecture for real time two way dynamic depth focusing | |
| CN108885258A (zh) | 具有一维片块的二维超声阵列换能器 | |
| US8235906B2 (en) | System and method for accelerated focused ultrasound imaging | |
| US20210256700A1 (en) | Systems and methods for ultrasonic imaging | |
| US20210338209A1 (en) | High quality high frame rate ultrasound imaging with diverging transmit beams | |
| Peretz et al. | Deep learning applied to beamforming in synthetic aperture ultrasound | |
| US20190187279A1 (en) | Rapid synthetic focus ultrasonic imaging with large linear arrays | |
| Dziewierz et al. | A design methodology for 2D sparse NDE arrays using an efficient implementation of refracted-ray TFM | |
| CN106037805A (zh) | 超声成像的方法及装置 | |
| CN106997045B (zh) | 基于超声系统点扩散函数测量和压缩感知的超声成像方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180629 |