CN102753104B - 增强波束的系统和方法 - Google Patents
增强波束的系统和方法 Download PDFInfo
- Publication number
- CN102753104B CN102753104B CN201080060936.3A CN201080060936A CN102753104B CN 102753104 B CN102753104 B CN 102753104B CN 201080060936 A CN201080060936 A CN 201080060936A CN 102753104 B CN102753104 B CN 102753104B
- Authority
- CN
- China
- Prior art keywords
- wave beam
- main lobe
- signal
- lobe
- wave
- 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 description 223
- 230000009467 reduction Effects 0.000 claims abstract description 15
- 238000007493 shaping process Methods 0.000 claims description 45
- 238000003384 imaging method Methods 0.000 abstract description 33
- 238000005516 engineering process Methods 0.000 abstract description 20
- 239000000523 sample Substances 0.000 description 208
- 230000008569 process Effects 0.000 description 175
- 230000015572 biosynthetic process Effects 0.000 description 121
- 238000003786 synthesis reaction Methods 0.000 description 121
- 230000006870 function Effects 0.000 description 64
- 238000000465 moulding Methods 0.000 description 34
- 230000011218 segmentation Effects 0.000 description 34
- 239000012071 phase Substances 0.000 description 25
- 230000003321 amplification Effects 0.000 description 21
- 238000003199 nucleic acid amplification method Methods 0.000 description 21
- 238000012545 processing Methods 0.000 description 20
- 230000002194 synthesizing effect Effects 0.000 description 19
- 230000008859 change Effects 0.000 description 16
- 238000001228 spectrum Methods 0.000 description 16
- 238000009826 distribution Methods 0.000 description 11
- 238000000354 decomposition reaction Methods 0.000 description 7
- 230000000877 morphologic effect Effects 0.000 description 7
- 208000030984 MIRAGE syndrome Diseases 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 6
- TVLSRXXIMLFWEO-UHFFFAOYSA-N prochloraz Chemical compound C1=CN=CN1C(=O)N(CCC)CCOC1=C(Cl)C=C(Cl)C=C1Cl TVLSRXXIMLFWEO-UHFFFAOYSA-N 0.000 description 6
- 230000015556 catabolic process Effects 0.000 description 5
- 238000006731 degradation reaction Methods 0.000 description 5
- 230000001965 increasing effect Effects 0.000 description 5
- 230000001629 suppression Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000003379 elimination reaction Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000002156 mixing Methods 0.000 description 4
- 239000011148 porous material Substances 0.000 description 4
- 230000002776 aggregation Effects 0.000 description 3
- 238000004220 aggregation Methods 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 3
- 230000008030 elimination Effects 0.000 description 3
- 239000003623 enhancer Substances 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000007620 mathematical function Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000002902 bimodal effect Effects 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 238000005304 joining Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000011112 process operation Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000012285 ultrasound imaging Methods 0.000 description 2
- 230000004304 visual acuity Effects 0.000 description 2
- 206010021403 Illusion Diseases 0.000 description 1
- 206010034719 Personality change Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000029142 excretion Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000008385 outer phase Substances 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000004540 process dynamic Methods 0.000 description 1
- 238000012958 reprocessing Methods 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000007920 subcutaneous administration Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 238000002604 ultrasonography 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
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/18—Methods or devices for transmitting, conducting or directing sound
- G10K11/26—Sound-focusing or directing, e.g. scanning
- G10K11/34—Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
- G10K11/341—Circuits therefor
-
- 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/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52046—Techniques for image enhancement involving transmitter or receiver
- G01S7/52047—Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
-
- 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/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52079—Constructional features
- G01S7/5208—Constructional features with integration of processing functions inside probe or scanhead
Abstract
本发明展示通过旁瓣降低和/或主瓣锐化来增强波束。各个实施方案利用动态分辨率、改善型动态分辨率和/或增强型动态分辨率技术来合成具有所需属性的波束,例如用于超声波成像的超声波束。各个实施方案同时针对每一样本形成第一样本波束和第二或辅助样本波束,以合成增强型扫描波束。根据本文的动态分辨率技术,可以根据所述两个样本波束的总和来形成新波束。各个实施方案的所合成动态分辨率波束具有降低的旁瓣和相对极小或相对不展开的主瓣。可以应用增强型动态分辨率波束锐化函数来提供进一步增强的波束,以便进一步窄化所述主瓣。
Description
相关申请的交叉引用
本申请要求以下美国临时专利申请的权益:2009年11月9日提交的、标题为“SYSTEMSANDMETHODSFORSCANBEAMSIDELOBEREDUCTIONWHILEREDUCINGTHEMAINLOBEUSINGDYNAMICRESOLUTION”的美国临时专利申请第61/259,346号;以及2009年11月10日提交的、标题为“SYSTEMSANDMETHODSFORSCANBEAMSIDELOBEREDUCTIONWHILEREDUCINGTHEMAINLOBEUSINGDYNAMICRESOLUTION”的美国临时专利申请第61/259,938号;所述美国临时专利申请的全部公开内容以引用的方式并入本文中。
技术领域
本发明涉及波束成形,且更明确地说,涉及例如通过旁瓣降低和/或主瓣塑形来增强波束的系统和方法。
发明背景
在声谱系统中,声音信号由扫描头发射到身体或其它对象中,且反射信号由扫描头接收以进行图像处理。反射信号由声谱系统用以形成身体结构(例如,患者的组织)或其它相关对象的图像。用于这种声谱成像中的扫描头通常是手持型封闭体,其含有一个或多个独立的换能器,且可能含有其它电子器件。
声谱系统扫描头的换能器在发射时将电能转换为机械能(声能)而从其表面辐射出去,且在接收时将撞击其表面上的机械能(声能)转换为电能。换能材料的个别部分称为一个元件,通常制造成例如长方形等特定几何形状。通常,这些换能器元件布置成规则的图案(阵列),其中中心布置为一条直线以形成线性阵列或相控阵列、沿着弧形成曲线阵列,或按照网格形成2D阵列。通常,换能器元件的这种规则图案在元件中心之间具有测量出的重复的间隔,称为间距。在声谱成像操作中,换能器元件通常成组地使用。这样一组换能器元件在一个维度上的总范围就是那一维度上的孔径。举例来说,对于线性阵列,一个维度是换能器元件的高度,而另一个维度是所用换能器元件数量乘以间距。
无论是在发射操作中还是在接收操作中,都可以通过适当地使用上述成组换能器元件来形成超声波束。举例来说,通过以下步骤来形成接收波束:调整换能器元件信号的一个或多个属性(例如,延迟和/或加权以提供与具有所选孔径的换能器元件相对应的换能器元件波束成形信号),并对这些换能器元件波束成形信号求和以提供具有与特定点(所述特定点是“焦点”)相对应的最大信号响应的波束成形信号。以上换能器元件信号属性在本文称为波束成形参数。这些波束成形参数通常用以形成波束以拒斥从不需要的区域(例如,不同于所需“观看方向”的方向)接收到的杂波(例如,不需要的反射信号,等等)。
确切地说,对来自换能器元件组的换能器元件信号施加延迟,使得如果从焦点发出窄脉冲,那么已被如此延迟的信号将在同一时间到达求和装置,且因此将得出最大值。来自焦点外的任何其它点的这种相同的窄脉冲将不会同时到达求和器,因此不会加和成像信号一样大。可以通过使用适当波束成形参数来形成具有特定形状(例如,宽度、长度、方向,等等)的波束。举例来说,可以形成“指向”所需“观看”方向的主瓣。
不依赖于施加延迟来创建波束,可以对孔径进行切趾。切趾(apodization)是在对换能器元件信号求和之前,将潜在唯一的增益值施加(加权)到这些换能器元件信号的过程。可以将特定切趾函数应用到孔径以产生具有所需属性(例如降低旁瓣)的波束,从而进一步拒斥杂波。存在许多标准加权函数可以应用到孔径,但具体举例说明其中三种函数。这三种函数是均匀加权(也称为矩形、棚车(boxcar)、辛克(sinc)或非切趾)、汉宁(Hanning,也称为Hann)加权和余弦加权。汉宁加权(1+cos(x))与余弦加权(cos(x))彼此相关之处在于,汉宁加权是升余弦函数。在超声成像系统中可用于孔径切趾的其它数学函数为汉明(Hamming)函数、布莱克曼-哈里斯(Blackman-Harris)函数或其它特定应用窗函数。
由均匀加权的孔径形成的波束称作辛克波束,由汉宁加权的孔径形成的波束称作汉宁波束,且由余弦函数加权的孔径形成的波束称作余弦切趾波束。通过依序变换超声波束(例如,辛克波束、汉宁波束或余弦切趾波束)来扫描对象,以形成图像。取决于实施方案,可以通过辛克波束或汉宁波束或其它类型的波束来形成超声图像。
当连续调整波束成形参数(例如,延迟),以使焦点沿着特定方向移动时,创建出动态聚焦波束。在提供波束扫描以进行声谱成像时,这些动态波束通常形成为使得焦点对于线性阵列来说沿循笛卡尔空间(Cartesianspace)中的直线,或在相控或曲线阵列中从顶点起沿着单一角。举例来说,通过依序调整换能器元件信号的波束成形参数,可以形成一系列波束以扫描相关体积(例如,可以扫描患者体内的特定面积或深度)。可以聚集来自多个这种扫描波束的信息,以产生相关扫描体积的图像(例如,患者皮下部分的超声图像)。举例来说,在超声B模式操作中,从不同观看方向的多个超声波束(例如,在不同观看方向扫描的波束)接收到的多行回声数据产生图像。这种从扫描波束产生图像在本文称为扫描体积成像。
已知,由汉宁波束获取的单色信号(monochromaticsignal)在数学上等于对获自辛克波束的信号与获自两个空间变换相邻辛克波束的信号均值进行求和,条件是这些波束根据尼奎斯特定理(Nyquisttheorem)隔开。也就是说,左侧辛克波束的第一无讯区(null)和右侧辛克波束的第一无讯区必须与中心辛克波束的峰值对准。基于这些特性,通过处理获自三个邻近的根据尼奎斯特定理隔开的辛克波束的信号,已提出一项技术来改善雷达应用的性能。然而,在超声成像中,根据多个系统参数来选择线密度以获得最佳图像质量,因此,通常无法满足根据尼奎斯特准则来设置波束间或扫描线间的取样间隔。此外,在超声成像中,通常结合可变孔径来实施动态波束成形(如同可能用于提供上述扫描样本波束)。换句话说,不同孔径大小用于在不同深度形成波束。因此,通常无法将例如可在雷达中实施的基于处理辛克波束的杂波减少技术用在超声扫描体积成像中。
图1A说明上述扫描体积成像。具体来说,图1A所示的换能器11(具有换能器元件E1到EN)可以操作以提供这种扫描体积成像。在操作中,对换能器元件E1到EN的换能器元件信号进行处理,以形成指向正被成像的体积15内的特定区域的接收波束。可以形成这些波束来收集关于正被成像的体积15内的对象(也称为相关对象)的信息,所述对象例如对象12(例如,流体填充区域)和表面16(例如,皮肤表面)下存在的对象13(例如,组织结构)。
应了解,用于扫描体积成像的信号(例如,波束成形信号)中的信号与杂波比越高,所产生图像中的对比度分辨率将越高(例如,组织区别越好)。信号杂波的一个来源是上述旁瓣,所述旁瓣通常伴随有所产生波束的主瓣。与所需主瓣相关联的不需要的旁瓣的存在可以从图1A的说明中看到。具体来说,图1A中所说明的主瓣各自具有与其相关联的旁瓣(例如,与主瓣ML5相关联的旁瓣SL5,主瓣与旁瓣的组合在虚线部分中示出以帮助区分这些瓣与复合表示)。旁瓣的数量和水平以及其结构决定有多少不需要的离轴回声整合到所得波束成形信号中,从而使相关对象的所需回声中产生杂波。降低旁瓣的能力改善相关对象(例如图像中的组织)的对比度分辨率或可区分性。
图像降级的另一个来源是用以收集图像信息的主瓣的宽度。举例来说,主瓣的宽度决定正被成像的体积内的对象被波束展开的方式。因此,主瓣的宽度通常与图像的细节分辨率有关。因此,通常需要形成用于上述扫描的波束具有窄焦距,以使得能够清晰地确定所产生图像中的相关对象。
从上述内容可以了解到,主瓣的宽度、旁瓣的水平和旁瓣的结构(例如,旁瓣从主瓣外转的速度)对图像质量有很大影响。举例来说,非常清晰的波束能够实现较高分辨率图像。
使用图1A的换能器11处理信号以产生图像可以包括通过对具有所选孔径的换能器元件所接收到的换能器元件信号适当地实施波束成形参数(例如,延迟和/或权重),来使用所选孔径(例如,一组所选换能器元件,例如换能器元件E11-E15)形成波束。举例来说,可以选择波束成形参数的延迟以提供具有所需焦点的主瓣ML11-ML15(例如,施加适当延迟提供波束,从而扫描正被成像的体积15的特定深度)。另外,波束成形过程可能涉及对从具有所选孔径的换能器元件接收到的信号施加适当权重(切趾过程),以便降低与主瓣相关联的旁瓣。因此,用于产生波束的波束成形参数可以包括复值(complexvalue),使得可以在量值和相位两方面修改从换能器元件接收到的信号。
尽管通常降低波束的旁瓣,但使用孔径切趾过程展开主瓣。与使用上述典型波束成形(使用切趾过程)相关联的不良结果由图1B到图1D示出。图1B示出组织模拟幻象150(大体表示图1A所示的正被成像的体积15),其由左侧的流体填充区域A(例如,可以与图1A正被成像的体积的包括对象12的部分相对应)和右侧的组织区域B(例如,可以与图1A正被成像的体积的包括对象13的部分)组成。假定组织区域B包括具有相等散射横截面的点散射体(例如,点散射体14)的群集。当正被成像的体积(此处由组织模拟幻象150表示)被元件E1-EN的线性阵列所形成的一连串超声波束穿透时,形成图像。由于将从流体填充区域A接收到小散射密度,所以所得图像(在理想情况下)不会有针对流体填充区域A显示的那种灰度,而组织区域B将显示与所述模拟幻象中所示情形具有类似密度的点分布。
如上所述,在常规超声成像系统中,用确定性数学函数切趾阵列的孔径,以部分地抑制旁瓣(进而加宽主瓣),从而改善图像对比度,或者不对孔径进行切趾以维持较窄主瓣,进而得出具有增多的杂波的较小成像点大小。每种情况都导致图像质量降级,且会显示失真的图像。
图1C示出上述两种不同波束配置以说明所述问题。波束BU是未切趾波束(例如,使用均匀加权函数确定波束成形加权分布所形成的辛克波束),其提供具有相对较高水平旁瓣的较窄主瓣。波束BH是切趾波束(例如,使用升余弦加权函数确定波束成形加权分布所形成的汉宁波束),其提供具有校对较低水平旁瓣的较宽(展开)主瓣。图1C中所说明的波束的量值成对数压缩,且旁瓣有圆齿,并逐渐外转。
假定使用波束BU和BH来对相同区域成像,具体来说,图1B的组织模拟幻象150的组织区域B的一部分。波束BU产生源自于由主瓣接收到的反射信号(例如,由点散射体14反射)的相关对象表示101(如可以用于聚集成通过在由组织模拟幻象150表示的区域内在不同观看方向上扫描多个波束BU而产生的图像)。波束BU进一步产生源自于由旁瓣接收到的反射信号的假影101-1到101-8(也可能会作为不需要的杂波聚集到所产生图像中)。同样,波束BH产生源自于由主瓣接收到的反射信号(例如,由点散射体14反射)的相关对象表示100(如可以用以聚集成通过在由组织模拟幻象150表示的区域内在不同观看方向上扫描多个波束BH而产生的图像)。波束BH进一步产生源自于由旁瓣接收到的反射信号的假影100-1到100-4(也可能会作为不需要的杂波聚集到所产生图像中)。如图1C中可见,尽管是对对象的相同区域成像,但与波束BU提供的相关对象表示101相比,波束BH提供的相关对象表示100展开。同样如图1C中可见,波束BU产生的假影(假影101-1到101-8)比波束BH产生的假影(假影100-1到100-4)多(尽管较小)。
可以通过扫描多个波束BU或波束BH以穿透正被成像的体积来产生声谱图像。举例来说,可以聚集通过在由组织模拟幻象150表示的区域中扫描波束BU和BH中的各别波束而创建的表示,以形成相关对象的图像。然而,如可以从图1C的说明了解,当使用波束BU时,所产生图像中的相关对象可能会相对尖锐,这是因为相关对象表示(例如,相关对象表示101)相对较小,但假影(例如,假影101-1到101-8)的数量由于较突出的旁瓣而较高。与使用波束BU相关联的假影也从相关对象表示中的相应表示延伸较长的距离,从而进一步使所产生图像降级。也如从图1C的说明可以了解,当使用波束BH时,所产生图像中的相关对象较不尖锐,这是因为相关对象表示(例如,相关对象表示100)相对较大,但相关联假影(例如,假影100-1到100-n)的数量由于较不突出的旁瓣而较低。此外,与使用波束BH相关联的假影从相关对象表示延伸较短距离。因此,每个上述波束成形技术都会导致所产生的图像通常比所需质量差。如从上述内容可以了解,实现清晰波束而无显著旁瓣以提供高质量成像已证明是不现实的。
发明概述
本发明是针对例如通过使用动态分辨率(DR)波束合成技术来提供波束旁瓣降低的系统和方法。通过从多个波束(称作样本波束)合成DR波束,本发明实施方案的动态分辨率波束成形技术除了提供波束旁瓣降低之外,还提供增强的波束主瓣属性。
各个实施方案通过对于正被成像的体积的每一扫描区域从样本波束获取多个波束成形信号来实施DR波束合成技术。举例来说,针对正被成像的体积(例如,组织区域)的每一扫描区域(例如,每一观看方向)形成第一样本波束(例如,可以使用辛克函数确定波束成形加权分布所形成的未切趾波束)和第二样本波束(例如,可以使用余弦函数确定波束成形加权分布所形成的切趾波束)两者。所得样本波束信号(例如,使用未切趾函数的波束成形信号和使用切趾函数的波束成形信号)用以通过本文的DR波束成形技术的操作来合成对应于来自高分辨率、低旁瓣波束的信号的波束成形信号。根据本发明的优选实施方案,可以在DR波束成形技术中对样本波束进行加权和组合,以得出最小的总功率。所得DR波束优选地具有降低的旁瓣和相对较小或相对未展开的主瓣。
在实施方案的改善型动态分辨率(IDR)波束合成技术中,对DR波束进行分段以合成具有所需属性的IDR波束。举例来说,可以例如使用样本波束(例如,上述第二样本波束)将DR波束分段成其主瓣分量和其旁瓣分量。这些波束分量优选地是独立地操纵或以其它方式进行处理,以便更改其一个或多个属性(例如,应用不同权重)。本文的IDR波束合成技术操作以从被操纵的分段波束分量合成IDR波束(例如,主瓣分量应用较大权重,而旁瓣分量应用较小权重),具体是通过重新组合这些波束分量以合成IDR波束。
如果需要,那么可以将锐化函数应用到DR/IDR波束,以提供更进一步增强的波束。被应用锐化函数的DR/IDR波束在本文称为增强型动态分辨率(XDR)波束。被应用锐化函数的XDR波束提供比相应DR/IDR波束窄的主瓣。另外,这种XDR波束的旁瓣可以在一定程度上得到进一步抑制以实现较高图像质量。
尽管IDR和XDR波束处理的实施方案可以利用如上所述的DR波束处理的波束合成,但本文概念的应用并不限于使用DR波束合成。举例来说,可以关于辛克波束和余弦切趾波束应用IDR和/或XDR波束锐化处理的实施方案,而不使用DR波束(例如,通过处理辛克和余弦切趾波束得到的最小功率波束)。
本发明实施方案的一个特征在于优化了每一波束成形样本在所产生图像的样本空间中的聚焦性能。本发明实施方案的另一特征在于最小化了频谱泄漏且改善了脉冲波(“PW”)、连续波(“CW”)和彩色血流处理中的频谱分辨率,具有或不具有编码激励和代码模式。本发明的又一特征在于实施方案可以轻易地调适成用于许多类型的系统,例如多线波束成形、合成孔径波束成形和高帧速率波束成形。
本文概念的实施方案可以应用到超声成像,以提供波束旁瓣降低。然而,本文概念并不限于关于超声成像的应用。可以对于可见光、红外线、射频和其它成像技术来应用各个实施方案。
上文已相当广泛地概述了本发明的特征和技术优点,以便可以更好地理解本发明的以下详细描述。下文将描述形成本发明权利要求书主旨的本发明的额外特征和优点。所属领域的技术人员应了解,所公开的概念和特定实施方案可以轻易地用作修改或设计用于执行本发明的相同目的的其它结构的基础。所属领域的技术人员还应认识到,这种等效构造并不脱离如由附属权利要求书所阐述的本发明的精神和范围。可以结合附图从以下描述更好地理解认为是本发明特征的新颖特征的组织和操作方法,以及其它目标和优点。然而,应明确理解,每一幅图都是仅为了说明和描述目的而提供的,并不旨在定义对本发明的限制。
附图简述
为了更完整地理解本发明,现在结合附图参考以下描述,在附图中:
图1A到图1C示出旁瓣降低的必要性和试图以常规方法降低旁瓣时所固有的问题;
图2A和图2B示出根据本发明实施方案的适于提供动态分辨率处理、改善型动态分辨率处理和/或增强型动态分辨率处理的系统;
图3A示出根据本发明实施方案的可以由动态分辨率波束合成技术使用的第一样本波束和第二或辅助样本波束的实施例;
图3B示出根据本发明实施方案的可以由图3A的样本波束合成的示例性动态分辨率波束;
图4示出使用本发明实施方案的DR波束产生的图像的表示;
图5A和图5B示出根据本发明实施方案的适于合成例如图3B的动态分辨率波束的动态分辨率波束的图2系统的实施方案的细节;
图6示出根据本发明实施方案的图5A和图5B的DR-XDR处理器的示例性操作;
图6A到图6C示出根据实施方案,根据图6的处理操作来组合第一样本波束信号和第二或辅助样本波束信号的信号,以进行动态分辨率波束合成;
图6D到图6I(3)示出根据实施方案,根据图6的处理操作来隔离主瓣信号分量与第一样本波束信号,以进行改善的动态分辨率波束合成;
图6J到图6L(4)示出根据实施方案,根据图6的处理操作来在增强型动态分辨率波束合成的迭代操作中使用波束塑形函数;
图7A到图7C为根据本发明实施方案的示例性样本波束和从其合成的增强型动态分辨率波束的曲线图;
图8示出应用到一维处理的本发明概念的一个实施例;
图9示出根据本发明实施方案使用不同参数集合处理的各种波束;
图10示出根据本发明实施方案处理的各个波束之间的关系;
图11示出将大放大因子应用到本文实施方案的辅助样本波束时,两个主瓣波束信号的量值之间的差异;
图12示出根据据本发明实施方案而利用的波束塑形函数的一个实施例;
图13A和图13B示出分裂以降低残余主瓣、动态分辨率波束的旁瓣以及所得波束的主瓣;
图14示出本发明实施方案的在多次迭代之后的增强型动态分辨率波束的实施例;
图15A到图15G示出根据本发明实施方案而产生具有不同波束特性的新分量信号;以及
图16A到图16E示出根据本发明实施方案的用于波束分解和合成的波束。
发明详述
图2A示出根据本发明实施方案调适的超声成像系统的实施方案。应了解,参考超声成像描述了示例性实施方案,以便提供更具体实施例来辅助理解本文概念。然而,本发明的概念并不限于关于超声成像的应用。因此,本文概念可以应用于其中使用发射信号的反射的许多技术,例如可见光、红外线和射频成像技术。
超声成像系统200被展示为包括与扫描头220通信的系统单元210。实施方案的系统单元210包括处理器式系统,所述处理器式系统可操作以控制扫描头220的换能器(例如,图1A所示的换能器11)使用扫描波束221发射和接收超声信号,以提供扫描体积成像。因此,实施方案的系统单元210的处理器式系统处理所接收到的超声信号以产生图像211、显示在显示器212上,从而表示正被成像的体积201的一部分。关于可以根据本发明概念加以调适的成像系统的细节提供于同在申请中且共同让渡的标题为“ModularApparatusforDiagnosticUltrasound”的美国专利申请第12/467,899号中,所述美国专利申请案的公开内容在此以引用的方式并入本文中。
关于超声成像系统200的实施方案的更多细节展示于图2B的高级功能方块图中。如图2B所示,扫描头220的换能器可以包括与发射/接收电路221(例如,可以包括放大器、缓冲器、多路复用器,等等)通信的超声元件(例如,具有换能器元件E1到EN的图1A的换能器11)阵列,且可操作以可控制地发射和接收超声信号。
图2B的系统单元210包括波束成形器213、DR/IDR/XDR波束合成处理器214、超声图像处理电路215以及显示器212。波束成形器213操作以对于提供给换能器11/从换能器11提供的信号提供波束成形。DR/IDR/XDR波束合成处理器214操作以提供如本文所述的动态分辨率波束合成处理。超声图像处理电路215操作以使用由DR/IDR/XDR波束合成处理器214合成的动态分辨率(例如,DR、IDR和/或XDR)波束信号来形成超声图像(例如,B模式、M模式、多普勒模式、3维、4维,等等),例如以供显示于显示器212上。
应了解,可以根据本发明实施方案利用图2B所说明功能方块的额外和/或替代功能方块。举例来说,例如在通过超声成像系统200实施数字波束成形或数字信号处理的情况下,可以利用一个或多个模拟数字转换器(ADC)和/或数字模拟转换器(DAC)。此外,功能方块的分布可以与图2A中所示不同。举例来说,例如在扫描头220与系统单元210之间需要“细线”连接的情况下,可以将波束成形器213安置在扫描头220中。
可以基于自适应性波束成形过程来优化与询问组织位置相关联的信号中的信号与杂波比,如可以由波束成形器213来实施。在自适应性波束成形过程中,针对每一样本位置使用对来自所有换能器元件(或换能器元件的某一选定子集)的信号的量值和相位的一系列矩阵操作。矩阵的维度通常与阵列孔径的维度成比例。当孔径如在常规超声成像系统中为较大(例如,为32、64或128)时,实施方案所需的处理功率非常高,且可能对于一些系统应用而言为过高的。
波束成形器213的样本波束波束成形器(samplebeambeamformer)213a和213b利用延迟与求和波束成形过程来根据实施方案产生扫描体积内询问组织位置的波束成形信号。延迟与求和波束成形是通过在补偿换能器元件之间的到达时间差异之后整合从换能器元件(所有换能器元件或换能器元件的某一所选子集)接收到的信号来完成。波束成形器的输出是取决于组织的散射横截面的分布而改变量值和相位的波束成形信号。
在根据实施方案的操作中,波束成形参数的第一集合(例如,延迟和/或权重的第一集合)提供第一样本波束,同时波束成形参数的第二集合(例如,延迟和/或权重的第二集合)提供第二样本波束。举例来说,样本波束波束成形器213a可以实施波束成形参数的第一集合以形成第一样本波束(例如,未切趾波束),而样本波束波束成形器213b实施波束成形参数的第二集合以形成第二或辅助样本波束(例如,切趾波束)。因此,样本波束波束成形器213a和213b同时提供两个不同的波束成形信号以用于根据本发明实施方案的动态分辨率处理。
从上述内容应了解,在波束成形信号中,来自多个换能器元件(所有换能器元件或换能器元件的某一选定子集)的信号被整合成单一波束成形信号。因此,在个别换能器元件处接收到的信号(例如,回声信号)的相位和量值丢失。因此,在波束成形之后改善波束成形信号是一项技术难题。然而,使用本文的DR、IDR和/或XDR波束合成概念,有可能改善波束成形后的波束性能。
超声成像系统200的实施方案实施本文描述的动态分辨率(DR)波束合成、改善型动态分辨率(IDR)波束合成和/或扩展动态分辨率(XDR)波束合成技术。举例来说,DR波束合成技术可以由超声成像系统200的DR/IDR/XDR波束合成处理器214实施,从而获取由波束成形器213的样本波束波束成形器213a和213b针对正被成像的对象(例如,组织)的每一扫描区域提供的多个波束成形信号以合成DR波束。所合成的DR波束的信号可以由超声成像系统200的DR/IDR/XDR波束合成处理器214进一步处理(例如使用本文描述的分段技术),以提供IDR和/或XDR波束合成。
根据本发明实施方案的DR波束合成技术从询问组织位置同时获取两个波束成形信号。举例来说,使用波束成形器213针对正被成像的体积201的每一扫描区域(例如,每一观看方向)中的每一样本点形成第一样本波束(例如,图3A的波束BU)和第二或辅助样本波束(例如,图3A的波束A)两者。一种使所合成的DR波束(例如,图3B的波束BO)中的旁瓣最小化的方式是形成相位、形状和量值在组合样本波束以合成DR波束时使旁瓣最小化的这种样本波束。因此,如可以从图3A的说明了解,实施方案的样本波束适于提供峰值和无讯区,所述峰值和无讯区协作以合成在组合时具有所需属性的DR波束。也就是说,与第一样本波束Bu(θ)相关联的旁瓣可以使用第二或辅助样本波束A(θ)来降低,以产生质量更好的DR波束Bo(θ)。
本发明的用于使旁瓣最小化的DR波束合成技术的实施方案可以用未被切趾而形成的第一样本波束和通过使用余弦函数(例如,cos(θ))对孔径进行切趾而形成的第二样本波束或辅助样本波束来实施。举例来说,来自阵列的未切趾波束图案(图3A的波束BU)可以由辛克函数辛克(θ)来描述。函数辛克(θ)在θ=±nπ或辛克(±π)=0而辛克(0)=1的零点交叉为振动的。如果来自阵列的波束图案是使用余弦函数来切趾(例如,图3A的波束A),那么余弦切趾波束对称地分裂成两个几何变换分量辛克波束,其两个峰值与未切趾波束的第一旁瓣对准,而无讯区位于原点A(0)=0处;A是第二样本波束或余弦切趾波束。因此,当应用于使用辛克函数形成的未切趾样本波束时,余弦切趾波束可以用作辅助样本波束以用于旁瓣降低。应了解,尽管本文中提及辛克波束和余弦切趾波束来演示DR波束合成过程,但根据本发明的概念,可以使用样本波束的其它组合来处理信号以用于重构具有改善的图像质量的图像。
由样本波束形成的样本波束成形信号由图2B的DR/IDR/XDR波束合成处理器214用来合成与本发明的DR波束(例如,图3B的BO)相对应的信号。举例来说,使用样本波束的几何和/或形态特性,对来自样本波束的信号进行组合以合成具有所需特征的DR波束。在样本波束提供组合起来会消除不良属性(例如,旁瓣)的几何和/或形态特性的情况下,实施方案的DR波束合成可以操作以对样本波束求和(例如,BO=BU+αA)。然而,在样本波束提供组合起来会增加不良属性(例如,旁瓣)的几何和/或形态特性的情况下,实施方案的DR波束合成可以操作以减去样本波束(例如,BO=BU-αA)。因此,应了解,本文关于使用样本波束合成DR波束而提供的数学关系可以取决于组合/消除所用特定样本波束的特征而实施正负号改变。
如从下文描述将更好地理解,上述实施例中的α定义从第二样本波束接收到的少量信号,用于消除第一样本波束中不需要的部分。根据实施方案,使用参数α来在DR波束合成中提供旁瓣减缓与主瓣展开之间的平衡。因此,通过本文动态分辨率波束成形技术操作合成的DR波束优选地具有降低的旁瓣和相对较小或未展开的主瓣。
提供以下论述以辅助更好地理解本发明实施方案的DR波束合成过程。当使用第一样本波束Bu(θ)来扫描正被成像的体积时,通过扫描对象O(θ)上的第一样本波束Bu(θ)而接收到的所得波束成形信号Iu(θ)可以被描述为Iu(θ)=∫O(φ-θ)*Bu(φ-θ)dφ。假定第一样本波束Bu(θ)可以分解为两个分量BuM(θ)和BuS(θ),其中波束分量BuM(θ)是所需波束分量(例如,主瓣),且波束分量BuS(θ)是不需要的波束(例如,旁瓣),且其中Bu(θ)=BuM(θ)+BuS(θ)。因此,Iu(θ)=∫O(φ-θ)*Bu(φ-θ)dφ=IuM(θ)+IuS(θ)。
第二样本波束A(θ)可以是辅助波束,其“观看”或指向与第一样本波束Bu(θ)的方向相同的方向(θ)。假定第二样本波束A(θ)也可以分解为两个分量AM(θ)和AS(θ),其中A(θ)=AM(θ)+AS(θ)。因此,IA(θ)=∫O(φ-θ)*A(φ-θ)dφ=IAM(θ)+IAS(θ)。从上文可以看出,可以根据Bo(θ)=Bu(θ)+αA(θ)=BuM(θ)+αAM(θ)+BuS(θ)+αAS(θ)形成DR波束Bo(θ)。因此,可以通过确定α以使‖BuS(θ)+αAS(θ)‖2的差异最小化或来降低旁瓣信号BuS(θ)。这样在效果上等同于
根据本发明实施方案,为了不使DR波束的主瓣降级,需要得自第二或辅助样本波束αIAM(θ)的主瓣的波束信号分量尽可能小。也就是说,非零αIAM(θ)通常导致所合成的DR波束Bo(θ)中主瓣BuM(θ)至少在某种程度上展开。为了在与第二样本波束A(θ)组合时,不更改样本波束Bu(θ)的观看或指向方向,实施方案的第二样本波束A(θ)包括与第一样本波束的主瓣BuM(θ)对应放置的无讯区。换句话说,如果A(0)=0,且Bo(0)=BuM(0),那么当来自旁瓣的杂波能量由于对消过程Bu(θ)+αA(θ)而得以最小化时,BuM(θ)的观看或指向方向不会更改,只是主瓣可能会由于BuM(θ)+αAM(θ)处理而稍稍展开。
可以通过预定最小值和最大值来限定参数阿尔法α。举例来说,实施方案的参数α优选地限定在0与1之间,以避免由数据获取或其它数字处理导致的错误。应了解,可能存在以下情形:通过使用来自第二样本波束(例如,上述余弦切趾波束)的信号来消除来自样本波束(例如,上述辛克波束)的杂波信号并不现实。举例来说,当杂波水平相当高时,或来自不需要的方向的信号过高以至于可能需要比用第二样本波束能够容易实现的波束更好的更高性能旁瓣消除波束时。在这种情形中,可以根据本发明实施方案,将参数α的值设置为预定最大可接受值(例如,1)。
α为零表示以下情形:旁瓣在第一样本波束(例如,辛克波束)中是低的,且不需要来自第二或辅助样本波束(例如,余弦切趾波束)的信号在DR波束合成中用于旁瓣消除。在对比度方面,当α等于1时,需要因为旁瓣而从第二或辅助样本波束(例如,余弦切趾波束)接收到的信号来用于在DR波束合成中消除第一样本波束(例如,辛克波束)的旁瓣。因此,在实施方案的辛克波束/余弦切趾波束DR波束合成实施例中,当α=0时,所合成的DR波束是辛克波束(BO=BU+αA=BU+0A=BU),而当α=1时,所合成的DR波束是升余弦或汉宁波束(BO=BU+αA=BU+1A)。在0与1之间的α值提供宽度在辛克波束主瓣与升余弦切趾波束主瓣之间改变的所合成DR波束主瓣,而旁瓣根据从余弦切趾波束旁瓣接收到的少量信号而降低。
可以实施优化过程以计算消除第一样本波束(例如,辛克波束)的不需要部分所需的来自第二或辅助样本波束(例如,余弦切趾波束)的信号量。然而在本发明实施方案的DR波束合成过程中,在这种优化过程中可以使用目标函数的不同选择。可以选择基于使波束中的杂波信号功率最小化的准则来计算上述参数α。也就是说,本发明实施方案的参数α有效地确定将与第一样本波束信号组合以合成DR波束的少量第二或辅助样本波束信号,从而可以根据实施方案选择用于选择参数α的目标函数,使得当组合样本波束时,杂波信号的功率得以最小化。
在本发明实施方案的DR波束合成过程中,可以在样本之间(例如,观看方向之间)动态地改变参数α。由于实施方案的参数α在旁瓣消除过程中是通过基于来自样本波束(例如,辛克波束和余弦切趾波束)的测量使波束的总杂波功率最小化来加以选择,所以可以对于DR处理图像中的每一样本实现细节分辨率与对比度分辨率的良好平衡。举例来说,在旁瓣高的样本位置,可以将参数α设置为高或上限(例如,在上述实施例中为1),且利用第二样本波束来合成具有低旁瓣的DR波束(其外转速率快)。然而,所合成的DR波束的主瓣宽度可能会受到第二样本波束的影响而展开,从而使得图像分辨率可能会较差。在从邻近位置接收到较少杂波的样本位置,可以将参数α设置为低或下限(例如,在上述实施例中为0),且在合成DR波束时实质上仍不使用第二样本波束。由于在这个实施方案中没有利用切趾,所以所合成的DR波束的主瓣较窄。因此,图像中的旁瓣消除参数在上述实施例中从0变化到1。因此,所述DR波束合成过程有效地在图像中的细节分辨率与对比度分辨率之间做出了尽可能最好的折衷。
为了说明上述概念,图3B的波束BO示出从样本波束BU和A的样本波束信号合成的DR波束Bo(θ)的表示(例如,Bo(θ)=Bu(θ)+αA(θ),其中α=1)。波束BO的旁瓣已大大降低。然而,如由图3A的W1和图3B的W2所示,主瓣与辛克波束的主瓣相比来说相对加宽。尽管这种展开的主瓣通常会降低所得图像的分辨率,但可以由这个实施方案的DR波束通过平衡相对轻微的主瓣展开与显著的旁瓣降低来提供改善的图像质量。
图4示出使用根据本发明实施方案合成的DR波束产生的图像的表示。具体地说,图4示出图像460,其展示正被成像的体积内的相关对象(例如,聚集以表示相关对象的相关对象表示400),其是由扫描样本波束(例如,图3A的样本波束BU和A)产生以合成DR波束(例如,图3B的DR波束BO)。如图所示,在图3B中,图像460的动态显示范围对应于在下部截止点GL与上部截止点GH之间的DR波束信号中的信息。图像460中没有旁瓣假影可见,这是因为所合成的DR波束的旁瓣已快速下降,且旁瓣创建的图像低于下部截止点GL。对应于点散射体14的所产生图像在图像460中仅稍微展开,这是因为DR波束的主瓣宽度减小了。因此,很好地保留了图像的纹理,且极少有或没有假影存在于图像460的区域A或B中。
图5A和图5B示出适于使用多个样本波束(例如,图3A的样本波束BU和A)合成DR波束(例如,图3B的DR波束BO)的超声成像系统200的实施方案的额外细节。图5A和图5B的示例性系统利用波束成形器213来形成两个样本波束信号以用于DR波束合成。所示实施方案的波束成形器213包括可在控制处理器(例如DR/IDR/XDR波束合成处理器214)的控制下操作的各个信号处理、加权和组合电路,以形成如本文所述的样本波束信号。
在所示实施方案中,从换能器元件(例如,图1A所示的换能器元件E1-EN)提供的信号由放大器51-1到51-N(例如,可以包括低噪声放大器)处理,从而提供换能器信号放大以用于随后的波束成形处理。元件52-1到52-N为每一换能器元件信号提供信号相位调整(延迟),且元件53-1到53-N为每一换能器元件信号提供信号振幅调整加权)。因此,元件52-1到52-N和元件53-1到53-N都将波束成形参数应用到换能器元件信号。组合器54和55提供换能器元件信号的组合(例如,求和),以形成所得波束信号。因此,元件52-1到52-N与组合器55协作以提供本发明实施方案的样本波束波束成形器213a,且元件52-1到52-N、元件53-1到53-N与组合器54协作以提供实施方案的样本波束波束成形器213b。
在图5A的实施方案中,将被适当延迟的换能器元件信号在组合器55进行组合,以提供辛克波束信号作为上述第一样本波束信号。将被适当延迟并加权的换能器元件信号在组合器54进行组合,以提供余弦切趾波束信号作为上述第二或辅助样本波束信号。在根据本发明实施方案的操作中,辛克波束与余弦切趾波束是同时形成(例如,使用同一换能器元件信号集合)。
如果需要,那么可以根据实施方案在DR波束合成处理之前、同时或之后提供波束信号调节和/或处理。举例来说,将所示实施方案的余弦切趾波束信号提供给正交带通滤波器56-1,而将辛克波束信号提供给正交带通滤波器56-2以促进向量空间中的DR波束合成处理。额外或替代样本波束信号调节可以包括模拟数字转换(例如,本文使用数字信号处理器(DSP)来合成DR波束)、放大、噪声消除,等等。
将辛克波束信号和余弦切趾波束信号提供给所示实施方案的DR/IDR/XDR波束合成处理器214以用于DR波束信号合成。图5A中所示的DR/IDR/XDR波束合成处理器214的实施方案包含多个波束处理电路,此处展示为DR处理511、IDR处理512和XDR处理513,以提供如本文所述的动态分辨率处理。然而,应了解,本发明实施方案可以不实施所示实施方案的所有波束处理电路。举例来说,如果需要,那么实施方案可以仅实施如本文所述的DR处理(例如,DR处理511)或DR处理与IDR处理的组合(例如,DR处理511和IDR处理512)。DR/IDR/XDR波束合成处理器214可以包括可在指令集合控制下操作以提供如本文所述操作的通用处理器、配置成提供如本文所述波束处理电路的DSP、特定应用集成电路(ASIC)、可编程门阵列(PGA),等等。在下文参考图6的处理更完整地描述这种DR/IDR/XDR波束合成处理器的操作。
应了解,根据本发明实施方案合成的波束信号用于图像产生,例如上述的超声图像产生。因此,实施方案的DR/IDR/XDR波束合成处理器214的输出被提供给用于这种图像产生的电路。举例来说,所示实施方案的所合成的动态分辨率波束信号被提供给检测和压缩电路59,以便去除信号的相位并提供信号量值的映射以进行扫描转换。实施方案的扫描转换器501将信号量值从获取空间变换到显示空间,以通过显示器212呈现给用户。
实施余弦切趾波束的问题之一在于组合换能器元件信号(例如,在图5A的组合器54处)。余弦函数在-1与+1之间振荡。因此,当执行加权与求和时,一个换能器元件信号(通道)可以是-1,而另一个换能器元件信号(通道)可以是+1,由此存在一个通道会与其它通道对消的可能性。与使用余弦函数相关联的每一通道的正负值可能在波束成形中引起延迟与求和过程的对消,从而导致动态范围限制。此外,由于余弦函数是从-1到+1且因此穿过零,所以某些换能器元件将会接收到极小的信号。为了克服这个问题,图像处理系统的前端电路(例如,前端电路的A/D转换器)可能需要具有极宽动态范围。
与使用余弦切趾相关联的上述问题在图5B的实施方案中得以避免,在图5B中,通过波束成形器213形成升余弦切趾波束(例如,汉宁波束)。因此,并不是使用波束成形器来直接形成余弦切趾波束,图5B的实施方案操作以除了形成用于DR波束合成的辛克波束之外,还形成升余弦切趾波束(例如,汉宁波束)。此后,可以通过取第一样本波束(此处为辛克波束)与升余弦切趾波束的信号之间的差来计算用于上述DR波束合成的第二或辅助样本波束信号。因此,在图5B的实施方案中,将被适当延迟的换能器元件信号在组合器55进行组合,以提供辛克波束信号作为上述第一样本波束信号。将在升余弦函数中加权的换能器元件信号在组合器54进行组合,以提供升余弦切趾波束信号。在前端形成这种升余弦切趾波束的优点在于组合器54仅需要处理正信号,从而进一步抑制噪声。
在超声成像中,通常结合可变孔径实施动态聚焦。换句话说,不同孔径大小用于在不同深度形成波束。一般来说,优选地随着深度增加而增加孔径大小,以维持图像中各个深度处的分辨率。可以使用具有随深度而变化的不同大小的孔径形成辛克波束和余弦切趾波束。由于在波束成形器中实施的通道数量通常有限,所以用于波束成形的孔径大小在某一深度停止增大,此时利用了可用于波束成形的所有通道。从这个特定深度开始,所有接收通道用于以恒定孔径形成波束。
在产生例如升余弦波束信号时,图5B中所示实施方案的放大器51-I到51-N的输入表示来自换能器的特定换能器元件的信号。对于N通道波束成形器来说,最大孔径包括N个换能器元件以形成具有N个通道的波束。当孔径随深度而变化时,某些计算可能仅需要使用换能器元件的子集(例如,可以针对选定换能器元件或选定孔径计算用于各个深度的升余弦函数)。在优选实施方案中,对于每一升余弦函数计算来说,实施方案的系统寻求最好的加权集合(例如,用于元件51-I到52-N和元件53-I到53-N中的适当元件的设置),以使用两个样本波束实现最好的DR波束。
在信号由图5B的组合器54和55组合之后,所得波束信号自身由组合器57组合,从而在所示实施方案中提供减法组合,以提供第二或辅助样本波束(此处为余弦切趾波束)以供用于根据实施方案的DR波束合成中。从以上论述应了解,如最初由组合器54产生的升余弦或汉宁波束可以是来自特定情形(例如,在实施方案中,参数α=1的情况)下的DR波束合成处理的所得波束。因此,并非从第一样本波束和第二或辅助样本波束重新产生升余弦切趾波束,而是实施方案可以利用最初由组合器54产生的升余弦切趾波束。因此,组合器54的输出被示为连接到DR/IDR/XDR波束合成处理器214,以除了第一样本波束(辛克波束)和第二样本波束(余弦切趾波束)之外,还将升余弦切趾波束信号提供给DR/IDR/XDR波束合成处理器214。
辛克波束信号和余弦切趾波束信号提供给所示实施方案的DR/IDR/XDR波束合成处理器214,以用于如本文所述的动态分辨率波束信号合成。因此,图5B实施方案的DR/IDR/XDR波束合成处理器214可以如上文关于图5A所述来配置。
如同上述图5A的实施方案,根据本发明实施方案,如果需要,那么可以在DR波束合成处理之前、同时或之后提供波束信号调节和/或处理。举例来说,在图5B的实施方案中,DR/IDRXDR波束合成处理器214在RF范围内操作(也就是说,信号被组合为RF信号),且所示实施方案的所合成的动态分辨率波束信号被提供给正交带通滤波器56以用于信号调节。额外或替代波束信号调节可以包括模拟数字转换(例如,本文使用数字信号处理器(DSP)来合成DR波束)、放大、噪声消除,等等。
如同图5A的实施方案,由图5B实施方案的DR-XDR处理器214合成的DR波束信号被提供给提供图像产生的电路。具体来说,所示实施方案的所合成的DR波束信号被提供给检测和压缩电路59,以便去除信号的相位并提供信号量值的映射以用于扫描转换。所示实施方案的扫描转换器501将信号的量值从获取空间变换到显示空间,以通过显示器212呈现给用户。
图6示出动态分辨率波束合成过程的实施方案的细节,所述过程可以由图5A和图5B的DR/IDR/XDR波束合成处理器214提供,以用于实现根据本文概念的扫描波束的旁瓣降低。具体地说,在图6所示的实施方案中,展示于上部点线上方的处理提供对应于图5A的DR处理511的操作的DR波束合成,展示于上部点线与下部点线之间的处理提供对应于图5A的IDR处理512的操作的IDR波束合成,且展示于下部点线下方的处理提供对应于图5A的XDR处理513的操作的XDR波束合成。
在实施方案的DR/IDR/XDR波束合成处理器214的DR处理511的操作中,第一样本波束信号和第二或辅助样本波束信号优选地使用参数α进行加权,并被组合以合成DR波束信号。因此,在处理601处,获取第一样本波束(例如,辛克波束),其具有针对正扫描对象(例如,组织区域)的扫描区域(例如,每一观看方向)中的每一样本点形成的主瓣601-1和旁瓣601-4。另外,获取第二或辅助样本波束(例如,余弦切趾波束),其具有针对在正被成像的体积的扫描区域中的每一样本点形成的主瓣601-2和旁瓣601-3。在根据实施方案的处理601的操作中,针对每一样本使用上述辛克和余弦切趾波束获取第一样本波束信号(包含主瓣信号分量和旁瓣信号分量)和第二或辅助波束信号
所示实施方案的第一样本波束信号和第二或辅助波束信号输入到处理62中,处理62计算参数α以将其用于对样本波束进行加权组合。如上所述,参数α是实施方案的DR波束合成处理中的旁瓣消除参数,且用于从样本波束(例如,辛克和余弦切趾波束)合成DR波束。在替代实施方案中,可以忽略参数α的计算,且因此,如由虚线610所示,处理流程可以直接进行到处理602。
以下向量分析有助于理解可以根据本发明实施方案使用的处理62中的参数α的计算。上述第一样本波束(例如,辛克波束)信号可以分解成两个分量,使得单位向量定义主瓣的观看方向,其中且单位向量定义旁瓣的观看方向,其中上述第二或辅助样本波束(例如,余弦切趾波束)信号可以由表示。当第二或辅助样本波束(例如,余弦切趾波束)与第一样本波束(例如,辛克波束)的旁瓣对准时,第二或辅助(例如,余弦切趾)波束的单位向量将与具有相反相位的第一样本波束(例如,辛克波束)的旁瓣分量的单位向量对准。也就是说,其中因此,第一样本波束信号的旁瓣分量沿着为 其中
应了解,本文的DR波束合成处理可以通过针对位于每一扫描波束的不同深度处的每一样本计算参数来应用到图像的所有样本。举例来说,和是在深度zn从第n个未处理第一样本波束(例如,上述实施例中的未切趾辛克波束)和第二或辅助样本波束(例如,上述实施例中的切趾余弦切趾波束)获取的信号。因此, 如上所述,根据本发明实施方案,需要参数α(n,zn)的值限定在0与1之间,或0≤α(n,zn)≤1。举例来说,在来自不需要方向的杂波在样本位置(n,zn)处非常大的情况下,那么在DR波束合成处理中将α(n,zn)设置成1,以最大化旁瓣抑制的量。在这个实施例中,将从第二或辅助样本波束(例如,余弦切趾波束)获取的信号100%相加到此时来自第一样本波束(例如,辛克波束)的信号。然而,这样会使主瓣由于求和过程而展开。也就是说,当α(n,zn)≠0时,一定程度的主瓣展开出现在样本(n,zn)中。因此,在实施方案的DR合成过程中,α(n,zn)取决于样本位置(n,zn)附近的杂波功率而在样本之间动态地改变。
情况α=0表示以下情形:从辛克波束接收到的杂波相对较小,因此不需要来自余弦切趾波束的信号来抑制旁瓣,从而拒斥杂波。在这种情况下,正被成像的对象由辛克波束的主瓣描绘,所述辛克波束的宽度根据衍射有限分辨率而展开。
情况α=1表示另一种情形:从辛克波束接收到的杂波非常大,以至于从余弦切趾波束接收到的100%信号都用以抑制旁瓣。这样导致汉宁波束的主瓣展开得比用于对象描绘的衍射有限分辨率大得多。用不同分辨率描绘正被成像的对象可能会引起用于图像解释的对象的感觉失真。另外,在杂波过大以至于汉宁波束(当α=1时)的旁瓣不足以抑制(或外转)杂波的极端情形中,样本会被杂波污染,从而使图像质量降级。
归于波束主瓣或旁瓣的信号分量可以在本发明实施方案的IDR和XDR处理中进行分段。通过操纵从DR波束和余弦切趾波束获取的信号,可以产生对应于具有不同形状、不同几何特性的波束的一类信号分量。可以合成主瓣更为尖锐且旁瓣较低的所需信号分量,从而等化主瓣分辨率与旁瓣水平,进而改善图像质量。
参数α(n,zn)指示杂波的量、其强度以及其在样本点附近的分布方式,且可以在本发明实施方案的IDR和XDR处理中用以控制其它波束成形参数,以锐化主瓣并弱化旁瓣。所示实施方案的参数表63包括查找表,用以使用α来映射所需波束成形参数。这些处理参数可以针对多种深度和条件进行预先确定,可以基于各个操作条件和参数来动态地计算,等等。在一个实施方案中,处理参数表63的参数是针对每一扫描波束信号而逐个孔径地设置成与深度和观看方向一致。
从上述内容可以看出,根据实施方案,通过以下步骤来形成实施方案的DR波束:组合来自辛克波束的信号与来自余弦切趾波束的信号的百分之α以降低旁瓣,其中α可以根据最小功率准则来确定,从而导致 换句话说,根据 合成实施方案的DR波束信号因此,本发明实施方案的处理602操作以从辛克波束信号与余弦切趾波束信号的加权总和形成DR波束信号举例来说,在处理602处,对应于所合成的DR波束Bo(θ)的DR波束信号可以从辛克波束Bu(θ)、信号与加权余弦切趾波束A(θ)、信号的总和来形成(例如,Bo(θ)=Bu(θ)+αA(θ))。
对应于在处理602形成的DR波束信号的所合成的DR波束Bo(θ)具有主瓣602-1以及降低或最小化的旁瓣602-3。图6A到图6C示出使用波束成形器213提供的样本波束、根据处理602的上述操作从辛克波束Bu(θ)、信号和加权余弦切趾波束A(θ)、信号的总和来合成对应于所合成的DR波束Bo(θ)的DR波束信号(例如,Bo(θ)=Bu(θ)+αA(θ))。具体来说,图6A示出从辛克波束(第一样本波束)接收到的辛克波束信号且图6B示出从余弦切趾波束(第二或辅助样本波束)接收到的余弦切趾波束信号图6C示出在α=1的情况下从辛克波束信号去除旁瓣分量所得到的合成的DR波束。
如上所述,一种用于合成DR波束信号的技术是组合辛克波束与余弦切趾波束信号与(例如,0≤α≤1)。然而,替代形成余弦切趾波束来合成DR波束,例如上述图5B实施方案的实施方案形成升余弦波束(例如,汉宁波束)以供用于合成DR波束。举例来说,处理602的操作使用如由图5B的波束成形器214提供的样本波束来提供DR波束信号合成,包括取升余弦波束的信号与辛克波束的信号之间的差,(0≤α≤1)。将设置为升余弦波束,则(0≤α≤1)。这个实施方案可以操作以从升余弦波束和辛克波束形成余弦切趾波束,从而使余弦切趾波束与辛克波束的旁瓣对准。在DR波束合成中使用这种余弦切趾波束实质上抑制了辛克波束信号的旁瓣,同时避免了辛克波束信号的主瓣的展开。
在处理602处合成的DR波束信号(提供受抑制的旁瓣和最小化的主瓣展开)可以由超声成像系统200用于产生高质量图像。然而,本发明实施方案提供额外动态分辨率波束合成处理,以进一步改善合成波束的特征。因此,根据所示实施方案的处理继续进行到处理603,以进行额外动态分辨率波束合成。
在实施方案的DR/IDR/XDR波束合成处理器214的IDR处理512的操作中,利用不同波束(例如,一个或多个样本波束和/或合成的DR波束)的某些几何和形态特性来进一步处理动态分辨率波束。具体地说,波束的信号分解为两个分量:对应于主瓣的分量和对应于旁瓣的另一分量。随后将这些分量信号重新组合(加权总和),以创建对应于具有极窄主瓣和极低旁瓣的新波束(本发明实施方案的IDR波束)的新信号。
图6所示的实施方案的处理603提供对所合成的DR波束信号和一个或多个样本波束信号(例如,第二或辅助样本波束信号)的处理,以对分量波束进行分段。随后将来自这些分量波束的信号用于构成新信号,就像它是从高性能波束(IDR波束)接收的一样。根据本发明实施方案,信号分解和重构处理在每一个样本点处、在每一观看方向上操作,从而优化整个图像的细节分辨率和对比度分辨率。
图6D到图6I中示出根据实施方案的通过处理603的操作提供IDR波束合成的波束分段、操纵和重新组合。如上所述,在参数α>0时,根据实施方案的DR波束合成除了导致旁瓣抑制之外,还将导致主瓣展开。如由处理603提供且示于图6D到图6I(对于α=1.)中的IDR波束处理操作以通过使用所合成的DR波束对来自第一样本波束(未展开,与第二或辅助样本波束组合)的主瓣进行分段来避免这种主瓣展开。为了对来自第一样本波束的主瓣进行分段,实施方案的处理603计算DR波束与第二或辅助波束之间的最小值,(其中提供相位对准),以给出如图6D所示的主瓣展开分量。这个主瓣展开分量实质上包括所合成的DR波束中并不是第一样本波束的主要波束的部分。因此,如图6E所示,可以从DR波束减去主瓣展开分量以提供第一样本波束主瓣分量,应了解,图6E的第一样本波束主瓣分量提供的主瓣比通过上述样本波束组合合成的DR波束的主瓣窄,从而可以用于提供高质量图像产生。
图6D和6E中所表示的示例性实施方案将第一样本波束(例如,辛克波束)的信号有效地分裂成两个分量:主瓣分量和旁瓣分量旁瓣分量信号的量值表示由样本点附近的第一样本波束的旁瓣所接收到的杂波信号。主瓣分量的量值表示在样本点处根据主瓣宽度而展开的第一样本波束的主瓣所穿透的区域中接收到的信号。信号与杂波比取决于正穿透的对象而在点之间有所不同。图像常常因此在不同程度上因杂波而有瑕疵。一种说法认为可以通过仅保留来自第一样本波束的主瓣的分量信号并从其旁瓣消除所述分量来重构清晰、无瑕疵的图像,虽然这个说法看起来似乎真实合理,但未必如此。已发现,对于来自旁瓣的信号量值相对于来自主瓣的信号量值来说较小的成像区域,去除旁瓣分量信号可以导致改善图像质量;但对于来自旁瓣的信号量值相对于来自主瓣的信号量值来说较大的其它区域,消除图像中来自旁瓣的所有信号可能会产生多个多斑点“暗”区域,从而可能会引起图像解释的问题。这是因为波束的旁瓣整合了从样本点附近的扩展区域接收到的回声信号。在从主瓣接收到的信号远低于从其旁瓣接收到的信号的样本点处,杂波分量信号在功能上插补了对象的形状,这样有助于图像解释。
因此,从通过对主瓣分量信号与少量旁瓣分量信号进行加权求和而形成的新波束导出的信号因此对于本发明实施方案来说是优选的。在本发明实施方案的处理603中,主瓣信号和旁瓣信号的权重是取决于来自主瓣的信号的量值、来自旁瓣的杂波信号的量值、旁瓣信号与主瓣信号比和/或每一样本点处的DR参数α来进行编程和选择,以优化图像质量。
根据实施方案的处理603的进一步动态分辨率处理由图6F和图6I表示,其中将主瓣和可以例如提供为适当加权的旁瓣分量的旁瓣外转相加到分段的第一样本波束主瓣分量以合成IDR波束。图6F中所示的旁瓣分量信号从第一样本波束(例如,辛克波束)的信号分段。已知,辛克波束的旁瓣水平以约6db/八度的速率(也就是说,与α>0时的DR波束相比以相对较低的速率)成角度(或距离)地外转。具体来说,对于汉宁波束(α=1时的DR波束),旁瓣以约18db/八度的速率(也就是说,与辛克波束(α=0时的DR波束)相比以快得多的速率)成角度地外转。可以通过取DR波束的信号与从辛克波束分段的主瓣信号之间的差来提取分量信号或
图6G示出对应于波束的信号所述波束的旁瓣按照DR波束(例如,α=1时的汉宁波束)而几何分布,由于减法过程而具有DR波束的小残余主瓣。如图6H和图6I所示,可以将信号减小参数γ(其中γ≤1),接着将其与信号整合以创建新信号 表示对应于波束(IDR波束)的信号,所述波束的瓣如图6I所示成角度(或距离)几何分布。IDR波束的主瓣的中心与辛克波束的主瓣一样窄,或比辛克波束的主瓣更窄。随后将主瓣的低振幅裙缘逐渐展开,并将其与以DR波束的速率外转的旁瓣混合。由于根据本发明实施方案将γ选择成小于1,所以旁瓣将减小20logγdb。例如,当γ=0.125时,旁瓣减小18db。可以在每一样本点处适应于其波束成形参数α、不同类型分量信号的量值或其比率来选择参数γ,以获得最佳图像质量。
应了解,先前实施例示出辛克波束的主瓣可以在IDR处理中从DR波束(当α=1时)分段。比辛克波束的主瓣窄的主瓣也可以从其它DR波束(当0<α<1时)分段。举例来说,当α=0时,如同在IDR处理中,对应于窄主瓣分量波束的信号可以通过计算 来分段。图6I(1)示出分段后的分量波束其比辛克波束的主瓣窄。另外,分量信号可以通过计算来分段,其中是对应于DR波束(当α=0时,如图6I(2)所示)的旁瓣的分量信号。可以通过计算形成新IDR波束。是对应于新波束的信号,所述新波束的旁瓣减小了18db;然而,所述旁瓣的外转速率与DR波束(对于α=0)的外转速率相同,为6db/八度。图6I(3)示出在对数标度下分别在α=0和α=1时的主瓣、旁瓣、DR波束的旁瓣外转速率和IDR波束。
应了解,可以根据本发明实施方案类似地分段对应于其它DR波束的主瓣的信号(当0<α<1时)。DR波束的分段主瓣的宽度随着α增加,且对于0≤α≤1,对于0≤α≤1来说,分段IDR波束的旁瓣减小20logγdb,且对于0≤α≤1来说,旁瓣外转速率随着α从6db/八度增加到18db/八度。
应了解,所合成的DR波束的旁瓣分量实质上受到抑制,但主瓣存在一定程度展开。使用在IDR处理中适当按比例调整的分段分量信号,使得DR波束的主瓣去除展开,且进一步改善了其旁瓣外转速率。选择IDR处理参数γ,以构成对应于所合成的波束的信号,所述波束的主瓣宽度比衍射有限波束窄或与衍射有限波束一样窄。IDR主瓣的量值逐渐外转,且与旁瓣(其水平由参数γ确定)混合。根据实施方案,IDR波束的旁瓣的峰值首先减小,接着继续沿着DR波束的旁瓣而以某个速率外转。这种IDR波束可以用以形成具有改善的感觉的图像对比度和细节图像分辨率的图像。
在处理603处合成的IDR波束信号提供抑制的旁瓣而没有主瓣展开,所述信号可以由超声成像系统200用于产生高质量图像。然而,本发明优选实施方案提供额外动态分辨率波束合成处理,以进一步改善所合成的波束特征。因此,根据所示实施方案的处理继续进行到处理604以进行额外动态分辨率波束合成。
在实施方案的DR/IDR/XDR波束合成处理器214的XDR处理513的操作中,所合成的波束(例如,上述IDR波束)的主瓣宽度通过本发明优选实施方案的扩展DR(XDR)波束合成过程来进一步锐化。根据本发明实施方案,波束通过XDR过程来渐进地或迭代地塑形,以获得更好的波束控制。XDR波束合成过程的波束锐化函数可以塑形被处理的波束的主瓣和旁瓣,以获得与从波束接收到的、具有极窄主瓣和极低旁瓣的信号类似的信号。举例来说,实施方案的XDR波束合成所得到的主瓣具有对应于换能器孔径的宽度,所述宽度实质上大于实际用于形成供合成XDR波束的样本波束的宽度。
图6所示的XDR波束合成处理的处理604创建锐化函数604-1,以用于合成具有锐化主瓣的XDR波束信号。在处理605处将这种锐化函数应用(优选地,迭代地应用)到主瓣分量(例如,主瓣605-1)以使波束变窄。如同上述IDR波束合成,将适当加权的旁瓣分量(例如,旁瓣605-3)添加到主瓣分量,以合成具有所需外转的XDR波束(例如,XDR波束605’)。所示实施方案的迭代控制606与处理604的锐化函数产生和处理605的锐化函数应用协作,以迭代地锐化所合成的XDR波束的主瓣。举例来说,在所示实施方案中,n是XDR处理实现所得的合成XDR波束信号(在处理607处输出)所需的迭代数量。图6所示XDR波束处理的每一迭代提供增强的新的合成波束,由所述迭代的XDR波束605’表示。
在锐化(窄化)时,操纵提供给实施方案的XDR波束合成处理的波束主瓣的属性(例如,第一未切趾样本波束的主瓣分量,其也对应于所合成的IDR波束的实施方案的主瓣),以提供波束塑形函数。在理解实施方案的波束塑形函数时,应了解,以上实施例(通过组合第一样本波束与第二或辅助样本波束(例如,Bo(θ)=Bu(θ)+αA(θ))合成DR波束)中的主瓣展开来源于以下情形:第二或辅助样本波束的量值在由第一样本波束的主瓣确定的角间隔中并非零。随着将来自第二样本波束的信号的百分之α加和到第一取样波束的信号以降低第一样本波束的旁瓣分量信号,来自第一样本波束的主瓣的信号分量也将以相同百分比与来自第二样本波束的信号分量相加。举例来说,如图3A所绘示,余弦切趾波束(所示实施方案中的第二或辅助波束)的双峰值主瓣由两个几何变换的辛克波束构成,一个位于辛克波束主瓣左侧,另一个位于辛克波束主瓣右侧。DR波束的主瓣在观看方向±θi(其中|θi|≤θα)处将截断余弦切趾波束的双峰值主瓣。当DR波束的旁瓣消除因子是α时,±θα是主瓣的零点交叉,其界定达2θα的波束分辨率。由于在实施方案的DR波束中,α被限定在0与1之间,所以当α=0时,|θ0|=π,且当α=1时,|θ1|=2π。因此,对获取自余弦切趾波束的信号与辛克波束的信号进行求和将会减少在|θ|≥θα区域中从辛克波束接收到的杂波。辛克波束的主瓣也会由于对-θα≤θ≤θα区域中的波束信号进行加和而展开。
然而,如果从余弦切趾波束接收到的信号被放大κ倍,那么余弦切趾波束与DR波束的交点θi的位置将取决于增益因子κ而沿着DR波束的主瓣上下移动。增益κ越大,θi越靠近主DR波束的中心。当κ接近零时,|θi|=θα。如果来自DR波束的分段主瓣的信号的量值与来自余弦切趾波束的放大后的信号相当且κ为放大因子,那么与的交点将位于±θDm,其中|θDm|≤π。放大因子κ越高,±θDm越靠近原点。注意,2|θDm|是主瓣的展开量,且值|θDm|越小,分辨率越好。举例来说,升余弦或汉宁波束(例如,所合成的DR波束,其中第一样本波束是辛克波束,第二或辅助样本波束是余弦切趾波束,且参数α=1)是与余弦切趾波束组合的辛克波束。因此,从通过取增益升高(例如,κ>1)余弦切趾波束与汉宁波束之间的最小值而计算出的汉宁波束减去分量信号,给出具有锐化主瓣的波束。所述过程有效地创建波束塑形函数以锐化汉宁波束的主瓣。随着汉宁波束的主瓣被锐化,汉宁波束的旁瓣结构也被修改,从而产生具有极窄主瓣和极低旁瓣的新波束,进而改善图像的细节分辨率和对比度分辨率。
图6J到图6L示出使用实施方案的波束塑形函数的XDR波束合成。在提供实施方案的波束塑形函数时,用来提供XDR波束合成处理的波束(例如,第一样本波束、所合成的DR波束、所合成的IDR波束或先前迭代合成的XDR波束)的主瓣分量通过图6J所示的处理604的操作来分裂成两个分量(例如,窄主瓣分量和残余主瓣分量)。举例来说,使用波束塑形函数(例如,图6K的波束塑形函数ψ(θ)将分裂成分量(窄主瓣分量)和(残余主瓣分量),其中且 可以通过如图6L中所示对分量信号进行加权求和来在处理605处合成XDR波束。举例来说,可以通过组合加权的(μ)窄主瓣分量加权的(ρ)残余主瓣分量和加权的(γ)旁瓣分量来形成XDR波束其中可以据此选择加权因子(μ、ρ和γ)。
随着DR波束的峰值与余弦切趾波束的无讯区对准,波束的峰值将不会被处理改变。
因此,信号的峰值在理论上等于信号的峰值;或在实施方案中,切趾处理中的数字错误可能会引起第一样本波束与第二样本波束的微小不对准,从而可能导致处理波束振幅降低。引入比例因子μ来进行增益等化,使得在第一迭代处,对应于第一样本波束(辛克波束)的主瓣的信号分解成两个分量和由于是来自比与对准的窄的波束的信号,且是远离波束中心的信号分量,所以可以通过组合信号与来形成新主瓣或 其中ρ≤1。当ρ=1时, DR波束的主瓣不改变。当ρ<1时,形成新主瓣所述瓣的中心窄,而所述瓣的外侧裙缘是通过减小的分量信号而形成。由于新主瓣是从分解自的分量和创建,所以瓣和以相同角间隔展开。信号的波束展开,连续连接到DR波束的旁瓣。
从上文可以了解,可以通过计算 来形成XDR波束。可以通过馈入在处理器处形成的XDR波束n次来使XDR处理迭代,直到在n次迭代的XDR波束 中产生令人满意的主瓣。每一次迭代的比例参数μn′s可以通过使参数ρn标准化而不损失一般性来设置为1。参数γn′s在每一次迭代中可以设置为不同值;然而,γn′s也可以设置成相等,以简化所有迭代中的计算。当XDR波束的主瓣在每一次迭代中分裂时,新主瓣裙缘处的信号量值根据参数ρi而减小,以将裙缘水平设置为比新主瓣低20logρidb。
图6L(1)到图6L(4)演示如何根据本发明实施方案在由处理604和605提供的XDR处理的一次迭代中形成XDR波束。用于处理实施例的参数设置如下:α=1;κ=2;μ=1;ρ=±0.125;且γ=0.015625。对应于第一样本波束(辛克波束)的主瓣的信号在IDR处理中如图6L(1)所示而分段。通过设置κ=2来升高第二样本波束(余弦切趾波束)的增益,使信号分解为两个分量:和其中这两个信号分量接着被组合成所述信号对应于由的波束构建的窄主瓣和对应于的各瓣的减小裙缘构成的新主瓣。瓣减少系数ρ被设置为0.125,这意味着新主瓣下降18db而过渡性地外转,其中裙缘开始与旁瓣混合。ρ也可以设置成负值,从而迫使主瓣如图6L(2)所示穿过零。已发现,负第一旁瓣增强对象的边界,从而改善一些成像应用的组织区分。波束的远侧旁瓣与减小了因子γ=0.015625(或减小36db)的DR波束的分段旁瓣混合。XDR波束、第一样本波束(辛克波束或DR波束,对于α=0,)和DR波束(对于α=1)在图6L(3)中重叠。图6L(4)以对数标度绘制了XDR波束。如所演示,在所述实施例中,XDR波束的主瓣比辛克波束窄得多,且旁瓣水平实质上低于以18db/八度外转的DR波束。可以重复处理604和605的操作以渐进地窄化所合成的XDR波束主瓣。
如同上文合成的IDR波束,通过以上XDR处理的操作分段的主瓣分量提供窄波束。然而,本发明实施方案操作以组合加权残余主瓣和旁瓣分量(例如,0<ρ<1且0<γ<1),以关于所合成的XDR波束提供所需外转和旁瓣水平。
应了解,也可以通过在数学上组合基于本文原理分段的不同分量信号来产生具有不同波束特性的新分量信号。举例来说,如图6D所示的分段分量信号可以与根据创建的分量信号组合,以获得新分量信号其中是对应于展开程度小于汉宁主瓣但大于其它分段主瓣(例如对于所有α′s,和来自任何波束的主瓣的信号)的主瓣的分量信号。这种处理用图形绘示于图15A到图15C中。三个分量信号(各自从具有不同波束宽度的分段主瓣获取,其中在图15D中重叠。取信号与的差以获得信号其表示从图15E所示的旁瓣分量波束获取的信号。类似地,可以通过取信号与之间的差来分段信号其表示从图15E所示的旁瓣分量波束获取的信号。旁瓣分量信号可以按比例调整,且与组合以获得新信号,所述新信号等于从被塑形主瓣获取的信号。举例来说,对应于信号的波束示出于图15G中。被塑形主瓣可以与具有所需外转的其它旁瓣分量信号组合,以获得具有所需波束的信号。
本发明实施方案的DR、IDR和XDR信号分段技术是基于分量波束根据窄带几何和形态特性在辛克波束、余弦切趾波束、汉宁波束以及其它波束间的合成和分解。然而,可以看出,当宽带信号被波束成形和处理(如同一般成像应用一样)时,本文概念同样有效。图7A到图7C示出样本波束和根据上述概念从其合成的XDR波束的示例性空间时间轮廓的曲线图,以说明本发明波束合成的结果。具有32个换能器元件的阵列由居中于3.5MHz、具有2MHz带宽的高斯脉冲(Gaussianpulse)激发。将适当时间延迟和不同权重应用到每一元件所接收到的信号以聚焦60mm处的点,所述点的展开程度取决于波束特性而不同。图7A的样本波束71包括未切趾第一样本波束(例如,辛克波束),图7B的第二或辅助样本波束72包括切趾波束(例如,汉宁波束或升余弦切趾波束),且XDR波束73是根据本发明的概念使用第一样本波束与第二或辅助样本波束而合成的。图7C示出XDR波束73的所得主瓣比图7A的样本波束71的主瓣窄得多,且具有极少或没有旁瓣。同样,图7C的XDR波束73的主瓣比图7B的样本波束72(切趾波束)的主瓣窄得多。
图8进一步示出适用于一维处理的本发明概念的实施例。曲线图800示出根据本文教示的概念产生的辛克波束801、汉宁波束802和所合成的波束803。具体来说,曲线图800是傅立叶频谱(Fourierspectrum),其示出信号的频谱分量包括十个正弦波a0到a9,这些波以标准化的频率从-0.35到0.25展开。这些信号的动态范围是140dB,其中信号a3处于0dB,且信号a7处于负140dB。负133dB处不均匀分布的噪声被相加到测试信号,且低于此噪声水平。如可以看出,如果对测试信号进行傅立叶变换而不使用窗函数,那么仅分辨出强信号分量a1、a2、a3、a4和a5。其它信号a6、a7、a8、a9和a0被这些强分量的复合旁瓣破坏。每一频谱分量的频谱峰值表示辛克函数的主瓣,其旁瓣以6dB/八度的速率外转。强分量的旁瓣干扰了频率,且使得不能检测到a0、a9、a8、a6和a7。当对信号进行汉宁切趾时,总旁瓣水平较低,也就是说,以较快速率或18dB/八度外转。因此,分辨出额外信号a6。然而,如由汉宁波束802所示,主波束展开,且频谱分量a7、a8、a9和a0仍无法分辨。
余弦切趾傅立叶频谱可以通过用两个最邻近的相邻频谱分量的平均值来替换每一频谱分量而获得。使用未切趾分量和余弦切趾分量,可以计算出DR/IDR/XDR频谱。对于傅立叶频谱中的每一频谱分量,计算杂波消除参数α。每一频谱分量的旁瓣接着受到抑制而不能进行DR处理,随后在实施方案的XDR处理中设置μ=1,γ=1,ρ=1,及κ=2。如所合成的波束803所示,分辨出所有频谱分量,包含水平嵌入于噪声中的分量:a7、a8和a0。
图9示出使用不同参数集合处理的辛克波束905、汉宁波束906、DR波束901、XDR1波束902、XDR2波束903和XDR3波束904。辛克波束905未被切趾,且汉宁波束906被切趾。DR波束901是通过仅使用如上所述的DR处理而产生。通过使用本文所论述的XDR处理,可以产生XDR1波束902、XDR2波束903和XDR3波束904。
图10示出辛克波束1001、汉宁波束1002和余弦切趾波束1003。注意,余弦切趾波束用作导引波束来控制功率水平。余弦切趾波束用以上下移动增益,以便检测系统应切割辛克波束的位置。辛克波束与余弦切趾波束中的最小值将确定主波束穿过零轴的位置。
从上述内容可以了解,XDR波束合成处理的操作例如通过应用将DR波束信号分解成不同分量信号的非线性和线性信号处理,来进一步窄化超声波束的主瓣,和/或进一步降低旁瓣。这些分量信号可以接着用以合成对应于虚拟波束的新波束信号(XDR波束信号),所述虚拟波束的主瓣窄且旁瓣低。或者,可以通过实质上不形成上述DR波束,而直接进行到图6的XDR处理,例如通过维持α恒定(例如,1或0.5),来合成实施方案的XDR波束。
应了解,上述DR/IDR/XDR波束合成技术可以实施于许多不同成像技术中。举例来说,DR/IDR/XDR波束合成技术可以关于线性扫描转换、空间混合等等来实施。对于多波束空间混合,从空间上移位的波束检测到的信号可以被整合以减少连贯斑点。为了增强分辨率,可以同时形成未切趾波束和汉宁波束,以在发生混合处理之前从每一观看方向获得DR/IDR/XDR信号。
尽管上文已参考一维换能器阵列描述了实施方案,但应了解,本发明的概念应用于多维换能器阵列。举例来说,DR/IDR/XDR波束合成概念可以直接应用到二维(2D)波束成形。
上文已描述根据本发明实施方案提供DR、IDR和XDR波束合成,下文将提供关于DR、IDR和XDR波束合成的各个函数和实施方案的细节。应了解,下文所述函数和实施方案可以用于上述系统和方法中,以提供根据本发明实施方案的DR、IDR和/或XDR波束合成。
波束分解和DR/IDR/XDR波束合成
在根据本发明实施方案的操作中,可以通过使所形成的波束在每一观看方向上的功率最小化来计算在特定位置从DR波束获取的信号。由于波束的功率在每一波束位置处改变,所以根据本发明实施方案的功率最小化准则来优化DR波束的主瓣宽度和旁瓣水平。基于DR波束、第一样本波束和第二样本波束的几何和形态特性,可以对与不同分量主瓣和旁瓣波束相对应的信号进行分段。这些分量波束进一步分解,接着通过算数上操纵这些分量波束来合成具有所需特性的新波束。根据本发明实施方案的DR、IDR和XDR波束分解和合成处理,计算对应于每一新波束的信号。根据本发明实施方案,在图6的处理601、602、603、604和605中实施使用分解分量波束构建对应于具有所需特性的所合成波束的信号的处理。这些处理用图形详细绘示,图中使用了与第一样本波束、第二样本波束、未切趾波束(α=0)、汉宁切趾波束(α=1)、如图3A、图6A-6L(4)和图15A到图15G中所示针对其它α(0<α<1)的DR波束组合的信号。
本发明实施方案的波束分解和合成处理可以在软件和/或硬件配置中实施。如图5A和图5B中所示,这种波束分解和合成处理可以在QBP(正交带通)滤波之前或在QBP滤波之后实施。基于图5B的方法是优选的,因为所述处理在从第一波束和第二波束获取的信号由QBP滤波器滤波之后仍存在噪声的情况下更为稳健。
如本文所述,波束分解可以用于根据本发明实施方案的动态分辨率波束合成中。在根据本发明实施方案的样本波束分解和合成的另一个更详细实施例中,将Iu和Ic设置为通过使用以有角度地展开的辛克波束Bu和余弦切趾波束Bc接收到的两个获取的样本信号(例如,来自散射体分布的回声)。由于所接收到的信号的振幅是来自根据每一波束在每一有角方向上的振幅加权的被声穿透的散射体的所有回声的整合,所以这些样本信号可以表示为和因此,可以通过根据以下方程式对信号Iu与Ic求和来计算来自新波束Bα的信号:其中Bα=Bu+αBc。
如上所述,可以将从余弦切趾波束接收到的信号与从辛克波束接收到的信号进行组合以获得新信号。这个新信号有效地从旁瓣水平降低而其主瓣加宽的新波束接收。举例来说,设置关于信号其相应波束Bα的几何和形态特性随着参数α而改变。一般来说,当α大于或等于零而小于或等于1(0≤α≤1)时,主瓣随着参数α而单调加宽,而旁瓣水平单调降低,且旁瓣外转速率单调增加。
使用图16A和图16B绘示针对不同α值(0,0.25,0.5,1)的波束特性。当其中是汉宁波束(升余弦切趾波束)。当其中是辛克波束。在汉宁波束中,第一旁瓣比辛克波束的第一旁瓣低-18db。汉宁旁瓣以18db/八度的较快速率外转,而辛克旁瓣以6db/八度的速率外转。然而,汉宁波束的主瓣将比辛克波束的主瓣宽50%。
应了解,如图16C所示,当θ>±π时,辛克波束和余弦切趾波束以相反相位振荡。因此,在Bα=Bu+αBc处理中,在大于±π的角位置处,辛克波束旁瓣的量值消除出现于>±π的所有角位置处。然而,在±π<θ≤±2π角位置处(辛克波束的第一旁瓣位于此处),旁瓣消除处理导致变换零点交叉点且展开Bα波束的主瓣。
同样如图16C所示,辛克波束的第一次零点交叉位于±π处,因此当根据Bα=Bu+αBc形成波束Bα以获得信号时,由于位置±π处的辛克波束对Bα波束形成没有贡献,所以±π处的波束的量值将等于以α按比例调整的余弦切趾波束的量值,或在波束指向方向上,Bα波束的增益最大,而余弦切趾波束Bc的增益最小,或对于任何α,Bα(0)=Bαmax,Bc(0)=0。在角位置θ<±π处,对于任何α,波束Bα的增益始终大于αBc的增益。当0≤α≤1时,|αBc|≤|Bc|,因此,因此,如果取与之间的信号最小值,那么
其中Mn_α是从通过取两个波束Bα与Bc之间的所有角θ处的绝对最小增益而形成的最小波束接收到的信号;θα是波束Bc与Bα相交的角度,且θα≤π。图16C中所示的C0,C0.25,C0.5,C1是波束Bc与Bα之间的相交点;Z0,Z0.25,Z0.5,Z1是波束Bα分别对于α=0,α=0.25,α=0.5,和α=1的零点交叉点。注意,随着α从零增加到1,零点交叉渐进地从±π移动到±2π。随着参数α增加,波束Bα的主瓣展开,而旁瓣水平降低。
将信号的相位设置为与从波束Bα接收到的信号的相位相同,或那么
其中是信号的相位。
由于波束Bα与Bc在θα处相交,因此Bα(θα)=Bc(θα)。此外,在波束指向方向上,Bα(0)=Bαmax,且Bc(0)=0。从减去获得新信号或
对于0≤α≤1。
是振幅在θ=0处最大的信号;的量值在±θα从其峰值朝向零对称地展开。操作意味着信号分解成两个分量:和分量信号表示对应于分量波束Ψam_n_α的信号,其限定于-θα≤θ≤θα区域。信号的峰值与的峰值对准。分量信号是的残余信号,其对应于在区域θ<-θα和θ>θα中保留波束Bα的旁瓣结构的分量波束Ψαs_n_α。
辛克波束的主瓣角度限定在区域-π≤θ≤π中。由于θα≤|π|,所以分量波束的波束宽度比辛克波束的波束宽度窄。另外,Ψαm_n_0(0)=Bu(0)、波束Ψαm_n_0的最大增益和其指向方向也不会改变。
辛克波束BuM和波束Ψαm_n_α=0=Ψαm_n_0的主瓣在图16D中重叠。新分量波束Ψαs_n_0可以根据Ψαs_n_0=Bum-Ψαm_n_0从辛克波束Bum的主瓣分解。Ψαs_n_0是残余主瓣波束,其由峰值与分量波束Ψαm_n_0的无讯区对准的双瓣组成。通过最小化由分量波束Ψαm_n_0接收到的信号的振幅或功率,分辨率将渐近地趋向于波束Ψαm_n_0。
应了解,辛克波束的主瓣可以根据以下处理来分段。图16E示出波束Bu和Bc被求和以形成波束Bα=1(汉宁波束)。由于Bu(π)=0,且Bα=1(π)=Bc(π)=0。因此,可以通过取波束Bα=1与Bc的振幅之间的最小值来创建波束通过将波束的相位指派为与波束Bα=1的相位相同,当θ≤π时波束的量值将与Bc的量值相同。在区域θ>π中,的量值将与波束Bα=1的量值相同。因此,辛克波束Bum的主瓣可以接着通过来分段,且辛克波束的旁瓣可以通过根据Bus=Bu-Bum=Bu-Ψαm_n_1从辛克波束Bu去除主瓣Bum来分段。
可以根据以下处理来有效地计算来自这些波束的信号。将参数α设置为1:α=1,那么 由于 因此, (对于|θ|=π)。因此
通过令来将Mn_1的相位设置为与向量相同。当来自波束的信号量在减法处理中从信号去除时,无讯区将位于θ=±π处。由于辛克波束的主瓣宽度由-π与+π之间的有角区域确定,所以辛克波束的主瓣因此被分段:
因此,可以通过获得辛克波束的旁瓣;且可以通过获得的旁瓣。由于DRi信号是使用不同α以最小杂波功率在不同样本位置从不同波束接收到的信号,或给定和本发明实施方案的DR波束可以通过对主瓣与不同量的求和使得来合成;其中γ≤1。
结果示出,对于图像中的所有信号,与DR波束信号相关联的所合成的DR波束的波束宽度与衍射有限辛克波束的波束宽度相同,而与DR波束信号相关联的所合成的DR波束的旁瓣低于辛克波束和最小功率波束的旁瓣。
来自DR波束的信号包括两个分量信号:来自旁瓣的分量信号或降低的旁瓣和从主瓣接收到的分量信号。保持来自旁瓣的信号不变,来自分段的主波束的分量信号可以进一步分裂成分量信号,所述分量信号包含:对应于新主瓣的至少一个分量,其波束宽度远窄于来自的波束宽度;以及从残余信号分解成新旁瓣的其它分量,所述新旁瓣是从主瓣分裂而来。来自所需方向的波束形成主瓣,在所述主瓣处,接收到信号在同一方向上,因为余弦切趾波束中无讯区的位置,使得实际上没能从接收到任何信号。换句话说,实质上没有来自所需方向(例如,观看方向)的信号存在于中。
当将放大因子κ应用到时,仅放大来自不需要的位置的信号,而极少信号是来自所需方向。因此,两个信号与之间的量值差表示来自当κ较大时主瓣有效地较窄的波束的信号。此特性示出于图11中。
令分段主瓣为BuM(θ),且由BuM(θ)接收到的信号为 表示当根据权重分布BuM(θ)(从θ=-π到θ=π)将对象进行加权求和时所接收到的信号;其中BuM(θ≤-π)=0,BuM(θ≥-π)=0。在对应于信号的波束的观看方向,bm为波束BuM(θ)的最大增益,而在相同观看方向上,来自余弦切趾波束的信号
用放大因子κ放大从余弦切趾波束获取的信号在效果上等同于将增益κ应用到余弦切趾波束。如果将从余弦切趾波束接收到的信号从辛克波束的信号中减去,那么产生主瓣由波束塑形函数塑形的新波束。这种波束塑形函数是增益κ的函数。这种波束塑形函数不仅塑形主瓣,而且修改旁瓣的结构,因此同时降低其水平。由于余弦切趾波束将与辛克波束在相交点BuM(θm)=Bc(θm)处在θ=±θm相交。从上文可见,通过计算与的量值之间的最小值提取的分量信号具有以下特性:
其中;
其中是信号的相位;且
其中是信号的相位。
因此,分量信号和取决于放大的余弦切趾波束与DR波束的主瓣和旁瓣的互动方式以及这些分量波束在波束放大因子κ下的形态而改变其特性。
假设κ选为足够大,使得的情况。如果从减去信号那么将去除从旁瓣接收到的所有信号。同样,由于所述减法处理,使得在±θm处创建了新无讯区,其中波束与相交(BuM(θm)=Bc(θm),如图11所示)。这种情况在效果上等同于将来自DR波束(在区域-π≥θ≥π中,波束的主瓣限定于其中)的主瓣的信号分裂成两个新信号分量;一个来自新形成的窄化主瓣其由θ-m<θ<θm和θm≤π限定;另一个来自新形成的两个旁瓣,其处于区域π≥θ>θm和-π≤θ<-θm中。
使用从较窄主瓣的波束获取的信号构建图像改善图像质量。然而,在减法处理中完全消除来自分段旁瓣的分量信号可能会在图像中引入空洞,从而使图像质量降级。为了更好地控制从DR波束减去的分量信号量,根据实施方案引入参数η≤1,其中表示从DR波束去除以合成所需高质量波束的一部分信号
当将一定量的信号从接收自DR波束的信号减去时,新信号在效果上等同于从主瓣由波束塑形函数塑形的波束接收。这种受到处理波束的主瓣被塑形以产生较窄主瓣,且可以通过这种函数降低旁瓣。举例来说,
注意,上文示出主瓣由塑形函数ψs(θ)塑形,从而产生处于区域|θ|<θm中的新主瓣这个新主瓣可以表示为
由于当|θ|<θm时,辛克波束BuM(θ)的主瓣的增益不为零,所以
其中BuM_m(θ)=ψ(θ)BuM(θ),且这样示出主波束是根据余弦切趾波束的增益与辛克波束的增益的比率在界定波束相交点θm且|θ|<θm的区域中来塑形。被塑形波束的形态将取决于参数η和κ。
上文示出当|θ|<θm时,DR波束的主瓣被波束塑形函数ψ(θ)(以图形示于图12中)修改。实施方案的塑形函数ψ(θ)具有以下特性:
1.从Bc(0)=0开始,ψ(0)=1;换句话说,在所需观看方向,波束增益最大化,且在信号处理中不改变。
2.当η和κ都是零时,波束的形态与DR波束相同,且对于所有θ,ψ(θ)=1。
3.由于Bc(θ)小于主瓣区域中的BuM(θ),那么当|θ|≤θm时,ψ(θ)<1。因此,波束BuM_m(θ)始终比波束BuM(θ)窄。
4.由于Bc(θ)和BuM(θ)在波束观看方向上都是对称的,所以波束塑形函数是对称函数。
波束塑形函数ψ(θ)在区域|θ|≤θm中改变波束的增益。塑形波束的量和形态将取决于参数η和κ。当η和κ都为零时,由于Bc(θ)和BuM(θ)在波束观看方向上是对称的,所以波束塑形函数ψ(θ)是对称函数。由于|θm|≤π,所以BuM(θ)的波束宽度始终窄于BuM(θ)的波束宽度。
上述分量信号和取决于放大的余弦切趾波束与DR波束的主瓣和旁瓣的互动方式以及这些分量波束在放大因子κ下的形态而改变。当放大因子κ小时,位置±θm较靠近于±π。因此,被处理波束的主瓣将相对宽于用大放大因子κ处理的波束。当放大因子κ小时,余弦切趾波束与DR波束的旁瓣之间的互动较为复杂。被处理波束的旁瓣形态将取决于应用放大因子之后的DR波束的旁瓣的相对振幅和余弦切趾波束的旁瓣的量值,因为来自分量的信号优选地设置为与DR波束同相,且其量值表示两个信号和中的较小者。因此,当从DR波束减去时,旁瓣的信号将始终较小,从而导致旁瓣抑制。
用于对分量信号进行分段的减法处理提示,减和被减信号的相位应保持相同。对于仅关注信号量值的成像应用来说,实施仅量值操作同样有效。在这种情况下,需要减和被减信号的信号正负号保持相同。换句话说,在DR/XDR处理中用替换的任何操作。举例来说,当在将波束成形信号正交带通且抽取成实数据流和虚数据流之后,实施算法以处理所述波束成形信号时,用 替换在这种情境中,代表信号的实部,且是信号的虚部。如此一来,DR/XDR处理可以分别对带通数据的实部和虚部进行操作。被处理的信号分量也将包括两个部分:实部和虚部接着,将被处理信号的实部和虚部重新组合成接着检测信号对其进行压缩,并将其扫描转换为所得图像视频。
然而,对于仅2D灰度成像应用来说,在DR/IDR/XDR处理中可以忽略信号的相位。在这种情况下,在DR/IDR/XDR处理中仅需要信号量值。此外,在XDR处理中,替代处理信号量值,可以使用信号功率来减少计算。由于划分和检测处理被并入到实施方案的DR处理中,所以计算成本可能相对较高。对于低成本实施方案,可以跳过DR处理,从而为了降低成本而损失性能或速度。
也可以在对信号进行正交带通之前实施算法。接着,在信号减法处理中,用sign(IRF)替换以进行有效的分量信号分段。
如前所述,从辛克波束主瓣接收到的信号可以是通过处理信号及(当κ=1时)而分段的随后,可以通过处理信号和(通过设置κ>1)来进一步分裂。当将κ设置为小于1的值时,例如,设置κ=0.22225,也可以形成主瓣大小大致等于汉宁波束的波束。这个较宽主瓣可以使用如上所述的处理来进一步分裂。接着可以用不同分量信号合成新波束。
从获自DR波束的信号去除上述分量信号有效地创建波束锐化函数。这种函数锐化主瓣,并降低被处理波束的旁瓣。图13A示出经过以下处理的主瓣:使用κ=2分裂,设置ρ=0.0625以降低残余主瓣,并设置γ=0.125以降低DR波束的旁瓣(在这种情况下为的旁瓣)。所得信号和其相应的波束绘示于图13B中。值得注意的是,由于主瓣窄化和旁瓣降低而减少了假影,所以过取样相关对象,从而实现比当前可能情况更清晰的显示图像将是有利的。
应了解,可以反复应用取DR波束与在不同放大因子下放大的余弦切趾波束之间的最小信号,并将所述最小信号的相位设置为与DR波束同相的处理,以获得对应于具有所需窄宽度主瓣和所需低旁瓣水平的波束的信号。也就是说,将信号分段为分量并使用不同放大因子κ、旁瓣降低参数γ和波束塑形参数η塑形波束的XDR波束合成处理可以反复应用到被处理波束,以获得具有所需主瓣和旁瓣特性的新的被处理波束。可以根据成像参数定义这些特性,例如成像过程中的细节分辨率、对比度分辨率和动态范围。用不同处理参数进行10次迭代后的波束示出于图14中。
如果需要,那么可以通过减少(例如,乘以小于1的减少因子)主瓣的被分段残余旁瓣来实现更高旁瓣水平控制。举例来说,DR和XDR处理中的增益和减少因子也可以设置成α的函数,以适应于每一样本处旁瓣的功率。
上述分析是基于辛克波束、余弦切趾波束、汉宁波束和其它波束的窄带形成中的几何特性。然而,可以通过实验证明,当宽带信号被波束成形和处理(如同一般成像应用一样)时,本文概念同样有效。此外,所述概念可以直接应用于频谱分析、二维阵列波束成形、多波束空间混合和多波束平行波束成形。
注意,可以使用本文DR/IDR/XDR技术实施方案的变化来实现不同结果。举例来说,向量格式可以由实部和虚部表示,且DR、IDR和XDR可以在实部和虚部中单独进行处理。保持相位不变等效于保持实部的正负号或虚部的正负号不变。可以通过保持信号的正负号不变而非保持相位不变来进行RF域中的仅量值处理。
如果需要,那么可以实现应用增益和减少因子或将增益或减少表示为数学函数的不同序列。对于恒定α(任何α≤1)形成波束,随后进行XDR处理以获得分辨率增强和旁瓣抑制。使用两个波束,实施方案可以应用最小值来分段分量波束或分量信号,且随后混合分量信号以合成新波束,使用增益和减少因子以构建新波束。
实施方案可以用多个增益和减少因子递归地形成新波束,以达到新的高性能波束。因此,概念适用于较高维度波束成形。举例来说:
令κj为来自余弦切趾波束的信号的放大因子,且0<κj≤κmax;
令μ为被处理主瓣的增强因子;μj≥1;
Fρ为被处理旁瓣的减少因子;|ρj|≤1;
开始
令将来自被处理波束的第一被处理信号设置为
j=1
当j<n时
计算
方法1
方法2
j=j+1;
结束
输出和
示例性DR/IDR/XDR波束合成算法
为了辅助理解上述本发明的概念,下文提供可以由图2A、图2B、图5A和图5B的DR/IDR/XDR波束合成处理器214实施的示例性DR/IDR/XDR波束合成算法。应了解,所述算法仅为可操作以根据本发明的概念提供DR/IDR/XDR波束合成的算法的实施例。
令为在第x个波束和第y个样本处来自余弦切趾波束的信号;
令为在第x个波束和第y个样本处来自辛克波束的信号;
通过如下计算逐个样本提取信号
若α>1,设α=1,α<0,设α=0,
则计算
如果需要IDR波束合成,那么开始如下IDR处理:
分段成分量信号
将信号与混合成新信号(合成的IDR波束信号);
如果需要XDR波束合成,那么开始如下XDR处理:
形成新旁瓣分量
形成XDR波束
为了使用根据本发明实施方案的迭代方法形成XDR波束:
方法1
令κj为来自余弦切趾波束的信号的放大因子,且
0<κj≤κmax;
令μ为被处理主瓣的增强因子,μ≥1;
令ρj为被处理旁瓣的减少因子;
ρj≤1;κ=[k1k2…kn];ρ=[ρ1ρ2…ρn];
开始
将来自被处理波束的第一被处理信号设置为并将IDR主瓣分量保存在缓冲器中。
令 且
令XDR波束仅包括旁瓣分量;
令
j=1
当j<n时
计算
j=j+1;
结束
方法2
令κj为用于来自余弦切趾波束的信号的放大因子,且0<κj≤κmax;
令μj为被处理主瓣的增强因子;μj≥1;
令ρ为被处理旁瓣的减少因子;
ρj≤1;
κ=[k1k2…kn];μ=[μ1μ2…μn];ρ=[ρ1ρ2…ρn];
开始
令将来自被处理波束的第一被处理信号设置为
j=1
当j<n时
计算
j=j+1;
结束
输出
再次参考图6,在实施上述示例性方法时,在601中获取在深度z处来自第n个辛克波束和余弦切趾波束的信号和接着根据使的幂次最小化的准则来在62中计算波束成形参数α,并在603中形成DR波束
接着通过查找表63传递参数α,以获得可用于适用于成像应用的IDR和XDR波束分解和合成的波束成形参数μi、ρi、κi、γm和γS。由于参数α指示样本位置(n,z)处的所需信号与不需要的杂波的相对量,所以α小时杂波小,且α大时杂波大,可用于IDR和XDR波束分解和合成的参数μi、ρi、κi、γm和γS可以是α的函数。
所合成的波束获得最好图像质量所需的特性取决于多个系统参数而改变,例如:扫描头的元件间距、可用于波束成形器中的通道数量、图像中的线密度、所使用的平行波束的数量、声穿透信号的频率和带宽,等等;因此,参数μi(α)、ρi(α)、κi(α)、γm(α)和γS(α)的不同集合对于不同成像应用是优选的。可以基于所需图像特征、正被成像的对象、成像系统的配置等等来用实验方法确定特定函数关系。
DR波束和参数μi、ρi、κi、γm、γS在处理器603中用以形成所需IDR波束。DR波束、IDR波束和波束成形参数也可以用于在处理器604和605中迭代地形成XDR波束。
图6L(1)到图6L(4)示出使用参数κ=2,ρ=0.125,γ=0.015625、n=1分解α=1的DR波束,以通过使用上述方法1计算合成XDR波束的实施例。3dbFWHM是0.7122π,其比辛克波束的FWHM窄21.6%。由于对于α=1,DR波束是峰值旁瓣水平为-31db的汉宁波束。在XDR波束合成过程中,使用参数γ来减小汉宁旁瓣20*log10(γ)=-36.12db,从而导致-31+20*log10(γ)=-67.12,如图6L(4)所示。所合成的XDR波束的旁瓣外转速率与汉宁波束相同,为-18db/八度。
根据上述动态分辨率波束合成技术,实施方案的DR/IDR/XDR波束合成在合成XDR动态分辨率波束信号时执行以下操作,其中在深度z从第n个辛克波束获取到的信号:Iu(n,z)=Iur(n,z)+jIui(n,z)和在深度z从第n个余弦切趾波束获取到的信号:Ic(n,z)=Icr(n,z)+jIci(n,z):
对于‖Ic(n,z)‖≠0;
计算
计算IDR(n,z)=Iu(n,z)+αIc(n,z);
LetIDR_r(n,z)=real(IDR(n,z));IDR_i(n,z)=imag(IDR(n,z));
LetIc_r(n,z)=real(Ic(n,z));IDR_i(n,z)=imag(Ic(n,z));
令sIr=sign(IDR_r);sIi=sign(IDR_i)
计算mIr=sIr*min(|IDR_r(n,z)|,κ|Ic_r(n,z)|);
mIi=sIi*min(|IDR_r(n,z)|,κ|Ic_r(n,z)|);
设IXDR(n,z)=IDR(n,z);
若sign(IDR_r(n,z))*sign(Icr(n,z))>0&mIr≠|IDR_r(n,z)|
IXDR_r(n,z)=IDR_r(n,z)+(ρ-1)mIr(n,z)
结束
若sign(IDR_i(n,z))*sign(Ici(n,z))>0&mIi≠|IDR_i(n,z)|;
IXDR_i(n,z)=IDR_i(n,z)+(ρ-1)mIi(n,z);
结束。
尽管已详细描述本发明以及其优点,但应理解,可以在不脱离如由所附权利要求书界定的本发明的精神和范围的情况下,在本文进行各种改变、替代和更改。此外,本发明的范围不希望受限于本说明书中描述的处理、机器、制造、物质组成、构件、方法和步骤的特定实施方案。所属领域的普通技术人员将从本发明的公开内容容易了解,可以根据本发明利用当前存在或稍后将开发的、与本文描述的相应实施方案执行实质上相同功能或实现实质上相同结果的处理、机器、制造、物质组成、构件、方法或步骤。因此,所附权利要求书希望将这些处理、机器、制造、物质组成、构件、方法或步骤包含于其范围内。
Claims (20)
1.一种操作超声波系统的方法,包括:
从一个或多个超声波换能器元件接收换能器元件信号;
加权并且求和各所述换能器元件信号以产生第一波束和第二波束,其中所述第一波束包括具有宽度的第一主瓣和与所述第一主瓣相邻的一对第一旁瓣,并且其中所述第二波束包括具有与所述一对第一旁瓣大致对准的一对峰值并且具有与所述第一旁瓣的相位大体相反的相位的第二主瓣;以及
组合所述第一波束和第二波束以生成第三波束,所述第三波束包括宽度小于所述第一主瓣的所述宽度的第三主瓣。
2.根据权利要求1所述的方法,其中同时获得所述第一波束与所述第二波束。
3.根据权利要求1所述的方法,其中所述第一波束包括未切趾波束,并且所述第二波束包括切趾波束。
4.根据权利要求3所述的方法,其中所述第一波束包括辛克波束(sincbeam),并且所述第二波束包括余弦切趾波束。
5.根据权利要求1所述的方法,进一步包括:
选择加权因子,其中组合所述第一波束和第二波束包括将所述第二波束和所述加权因子之积与所述第一波束求和。
6.根据权利要求5所述的方法,其中选择所述加权因子进一步包括:
所述第二波束和所述加权因子之积确定所述第二波束的至少一个分数部分,当与所述第一波束组合时生成与所述第一波束和第二波束相比具有一个或多个减小的旁瓣的所述第三波束。
7.根据权利要求5所述的方法,其中选择所述加权因子包括:
将所述加权因子的值的选择限制在0与1之间包含0和1的值范围。
8.根据权利要求7所述的方法,其中所述选择所述加权因子进一步包括:
当所述第一波束与所述第二波束中的一个波束内的杂波高时,为所述加权因子选择所述范围中的高值。
9.根据权利要求7所述的方法,其中所述选择所述加权因子包括:
当所述第一波束与所述第二波束中的一个波束内的杂波低时,为所述加权因子选择所述范围中的低值。
10.根据权利要求5所述的方法,进一步包括:
将所述所述第一波束和所述第二波束的生成重复多个循环;以及
在所述多个循环的各循环间动态地改变所述加权因子;以及
将所述第二波束和所述加权因子之积与所述第一波束求和重复多个循环,以生成多个第三波束迭代。
11.根据权利要求10所述的方法,进一步包括确定如下的所述第三波束迭代中的一个或多个:所述第三波束迭代具有的主瓣的宽度小于所述第一主瓣的所述宽度并且与所述第一波束和所述第二波束相比具有减小的旁瓣。
12.一种操作超声波系统的方法,所述方法包括:
从一个或多个超声波换能器元件接收换能器元件信号;
加权并且求和各所述换能器元件信号以产生第一波束和第二波束,其中所述第一波束包括具有宽度的第一主瓣和与所述第一主瓣相邻的一对第一旁瓣,并且其中所述第二波束包括具有与所述一对第一旁瓣大致对准的一对峰值并且具有与所述第一旁瓣的相位大体相反的相位的第二主瓣;以及
生成第三波束,所述第三波束包括宽度小于所述第一主瓣的所述宽度的第三主瓣,其中生成所述第三波束包括将所述第二波束和加权因子之积与所述第一波束组合。
13.根据权利要求12所述的方法,进一步包括:
选择所述加权因子,其中选择所述加权因子包括:所述第二波束和所述加权因子之积确定所述第二波束的至少一个分数部分,当与所述第一波束组合时生成与所述第一波束和第二波束相比具有一个或多个减小的旁瓣的所述第三波束。
14.根据权利要求13所述的方法,进一步包括:
将所述第一波束和所述第二波束的生成重复多个循环;以及
在所述多个循环的各循环间动态地改变所述加权因子;以及
重复所述第三波束的生成,以获取多个第三波束迭代。
15.根据权利要求14所述的方法,进一步包括确定如下的所述第三波束迭代中的一个或多个:所述第三波束迭代具有的主瓣的宽度小于所述第一主瓣的所述宽度并且与所述第一波束和所述第二波束相比具有减小的旁瓣。
16.根据权利要求12所述的方法,进一步包括:
计算所述第一波束的主瓣分量,其中所述主瓣分量包括所述第一主瓣的至少一部分;
将所述主瓣分量分段成第一主瓣部分和第二主瓣部分;
计算所述第三波束的旁瓣分量;
利用所述第一主瓣部分、所述第二主瓣部分和所述第三波束的所述旁瓣分量的加权和来计算第四波束。
17.根据权利要求16所述的方法,其中所述第四波束包括第四主瓣,其具有的宽度小于所述第一主瓣和所述第三主瓣的宽度。
18.一种操作超声波系统的方法,所述方法包括:
从一个或多个超声波换能器元件接收换能器元件信号;
加权并且求和各所述换能器元件信号以产生第一波束和第二波束,其中所述第一波束包括第一主瓣和与所述第一主瓣相邻的一对第一旁瓣,并且其中所述第二波束包括具有与所述一对第一旁瓣大致对准的一对峰值并且具有与所述第一旁瓣的相位大体相反的相位的第二主瓣;以及
组合所述第一波束和所述第二波束以生成第三波束,所述第三波束具有第三主瓣,其中组合所述第一波束和第二波束包括:
选择加权因子,其中选择所述加权因子包括:所述第二波束和所述加权因子之积确定所述第二波束的至少一个分数部分,当与所述第一波束组合时生成与所述第一波束和第二波束相比具有一个或多个减小的旁瓣的波束;以及
将所述第二波束和所述加权因子之积与所述第一波束求和;
将所述第一主瓣分段成第一主瓣分量和第二主瓣分量;以及
利用所述第一主瓣分量、所述第二主瓣分量和所述第三波束的旁瓣分量的加权和来生成第四波束,其中所述第四波束包括如下主瓣,该主瓣具有的宽度小于所述第一主瓣和所述第三主瓣的宽度。
19.根据权利要求18所述的方法,其中生成所述第四波束包括将所述第一主瓣分量、所述第二主瓣分量和所述第三波束的旁瓣分量分别乘以第一、第二和第三波束成形参数。
20.根据权利要求19所述的方法,其中所述第一、第二和第三波束成形参数基于所述加权因子获得。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US25934609P | 2009-11-09 | 2009-11-09 | |
US61/259,346 | 2009-11-09 | ||
US25993809P | 2009-11-10 | 2009-11-10 | |
US61/259,938 | 2009-11-10 | ||
PCT/US2010/056000 WO2011057252A1 (en) | 2009-11-09 | 2010-11-09 | Systems and methods for beam enhancement |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102753104A CN102753104A (zh) | 2012-10-24 |
CN102753104B true CN102753104B (zh) | 2016-03-09 |
Family
ID=43970420
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201080060936.3A Expired - Fee Related CN102753104B (zh) | 2009-11-09 | 2010-11-09 | 增强波束的系统和方法 |
Country Status (6)
Country | Link |
---|---|
US (3) | US8876719B2 (zh) |
EP (1) | EP2498683B1 (zh) |
JP (2) | JP5973349B2 (zh) |
CN (1) | CN102753104B (zh) |
BR (1) | BR112012010958B1 (zh) |
WO (1) | WO2011057252A1 (zh) |
Families Citing this family (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9282945B2 (en) * | 2009-04-14 | 2016-03-15 | Maui Imaging, Inc. | Calibration of ultrasound probes |
KR101659723B1 (ko) | 2009-04-14 | 2016-09-26 | 마우이 이미징, 인코포레이티드 | 복수 개구 초음파 어레이 정렬 설비 |
WO2011057252A1 (en) * | 2009-11-09 | 2011-05-12 | Sonosite, Inc. | Systems and methods for beam enhancement |
JP6274724B2 (ja) | 2010-02-18 | 2018-02-07 | マウイ イマギング,インコーポレーテッド | 多開口超音波撮像を用いた点音源送信及び音速補正 |
US8903012B2 (en) * | 2010-05-04 | 2014-12-02 | Hughes Network Systems, Llc | Continuous phase modulation for satellite communications |
EP2627257B1 (en) | 2010-10-13 | 2019-04-17 | Maui Imaging, Inc. | Concave ultrasound transducers and 3d arrays |
WO2012051305A2 (en) | 2010-10-13 | 2012-04-19 | Mau Imaging, Inc. | Multiple aperture probe internal apparatus and cable assemblies |
EP2638860B1 (en) * | 2010-11-09 | 2020-02-26 | Konica Minolta, Inc. | Beam-forming method, ultrasonic diagnosis device, program, and integrated circuit |
JP6049371B2 (ja) * | 2011-11-09 | 2016-12-21 | 東芝メディカルシステムズ株式会社 | 超音波診断システム |
EP2785253B1 (en) | 2011-12-01 | 2023-11-15 | Maui Imaging, Inc. | Motion detection using ping-based and multiple aperture doppler ultrasound |
CN104080407B (zh) | 2011-12-29 | 2017-03-01 | 毛伊图像公司 | 任意路径的m模式超声成像 |
JP6438769B2 (ja) | 2012-02-21 | 2018-12-19 | マウイ イマギング,インコーポレーテッド | 多数開口超音波を用いた物質の硬度の決定 |
CN104203110B (zh) | 2012-03-26 | 2017-06-06 | 毛伊图像公司 | 用于通过应用加权因子来提高超声图像品质的系统和方法 |
JP5489144B1 (ja) * | 2012-04-27 | 2014-05-14 | コニカミノルタ株式会社 | ビームフォーミング方法、及び超音波診断装置 |
JP5412604B1 (ja) | 2012-06-04 | 2014-02-12 | パナソニック株式会社 | 超音波診断装置、およびビームフォーミング方法 |
EP2883079B1 (en) | 2012-08-10 | 2017-09-27 | Maui Imaging, Inc. | Calibration of multiple aperture ultrasound probes |
EP3893022A1 (en) | 2012-08-21 | 2021-10-13 | Maui Imaging, Inc. | Ultrasound imaging system memory architecture |
JP6039305B2 (ja) * | 2012-08-23 | 2016-12-07 | キヤノン株式会社 | 被検体情報取得装置、情報処理装置および被検体情報取得方法 |
JP6103857B2 (ja) * | 2012-08-28 | 2017-03-29 | キヤノン株式会社 | 被検体情報取得装置、表示方法、及びプログラム |
US9510806B2 (en) | 2013-03-13 | 2016-12-06 | Maui Imaging, Inc. | Alignment of ultrasound transducer arrays and multiple aperture probe assembly |
US9211110B2 (en) | 2013-03-15 | 2015-12-15 | The Regents Of The University Of Michigan | Lung ventillation measurements using ultrasound |
JP2015071028A (ja) * | 2013-09-05 | 2015-04-16 | セイコーエプソン株式会社 | 超音波測定装置、超音波画像装置及び超音波測定方法 |
JP6197505B2 (ja) * | 2013-09-05 | 2017-09-20 | セイコーエプソン株式会社 | 超音波測定装置、超音波画像装置及び超音波測定方法 |
US9883848B2 (en) | 2013-09-13 | 2018-02-06 | Maui Imaging, Inc. | Ultrasound imaging using apparent point-source transmit transducer |
US10401493B2 (en) | 2014-08-18 | 2019-09-03 | Maui Imaging, Inc. | Network-based ultrasound imaging system |
US10238298B2 (en) * | 2015-01-08 | 2019-03-26 | Nec Corporation | Method and apparatus for photoacoustic tomography using optical orbital angular momentum (OAM) |
US10430688B2 (en) * | 2015-05-27 | 2019-10-01 | Siemens Medical Solutions Usa, Inc. | Knowledge-based ultrasound image enhancement |
US9824442B2 (en) * | 2015-08-20 | 2017-11-21 | Siemens Medical Solutions Usa, Inc. | View direction adaptive volume ultrasound imaging |
CN106341165B (zh) * | 2015-11-17 | 2020-02-04 | 北京智谷睿拓技术服务有限公司 | 波束呈现方法和设备 |
EP3408037A4 (en) | 2016-01-27 | 2019-10-23 | Maui Imaging, Inc. | ULTRASONIC IMAGING WITH DISTRIBUTED NETWORK PROBES |
JP6289785B2 (ja) | 2016-01-29 | 2018-03-07 | 三菱電機株式会社 | アンテナ装置及びグレーティングローブの低減方法 |
JP6467386B2 (ja) | 2016-09-14 | 2019-02-13 | 矢崎総業株式会社 | 車両用表示装置 |
JP6819184B2 (ja) * | 2016-09-26 | 2021-01-27 | セイコーエプソン株式会社 | 超音波測定装置および超音波測定装置の制御方法 |
KR101858137B1 (ko) * | 2017-03-08 | 2018-05-16 | 대진대학교 산학협력단 | 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법 |
EP3382423A1 (en) | 2017-03-27 | 2018-10-03 | Koninklijke Philips N.V. | Methods and systems for filtering ultrasound image clutter |
WO2019057592A1 (en) * | 2017-09-22 | 2019-03-28 | Koninklijke Philips N.V. | METHODS AND SYSTEMS FOR ULTRASONIC CONTRAST STRENGTHENING |
EP3598950A1 (en) | 2018-07-24 | 2020-01-29 | Koninklijke Philips N.V. | Ultrasound controller unit and method |
US11125867B2 (en) * | 2018-10-10 | 2021-09-21 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Adaptive weighting for adaptive ultrasound imaging |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1708698A (zh) * | 2002-11-08 | 2005-12-14 | 皇家飞利浦电子股份有限公司 | 用于自动设置超声换能器阵列的发射孔径和旁瓣缩减的方法和设备 |
Family Cites Families (49)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5601083A (en) * | 1986-05-02 | 1997-02-11 | Anderson; Forrest | Real time 3D imaging device using filtered ellipsoidal backprojection |
US5570691A (en) * | 1994-08-05 | 1996-11-05 | Acuson Corporation | Method and apparatus for real-time, concurrent adaptive focusing in an ultrasound beamformer imaging system |
JPH10142334A (ja) * | 1996-09-12 | 1998-05-29 | Secom Co Ltd | 超音波計測装置 |
US5935068A (en) * | 1996-11-01 | 1999-08-10 | The Trustees Of The University Of Pennsylvania | System and method for improving ultrasound image contrast by amplitude compression of ultrasonic wavefront signals |
US7104956B1 (en) * | 1996-11-08 | 2006-09-12 | Research Corporation Technologies, Inc. | Finite amplitude distortion-based inhomogeneous pulse echo ultrasonic imaging |
US6023977A (en) * | 1997-08-01 | 2000-02-15 | Acuson Corporation | Ultrasonic imaging aberration correction system and method |
US6312379B1 (en) * | 1997-08-15 | 2001-11-06 | Acuson Corporation | Ultrasonic harmonic imaging system and method using waveform pre-distortion |
US6111535A (en) * | 1997-08-25 | 2000-08-29 | Veridian Erim International, Inc. | Method of minimizing leakage energy in a synthetic aperture radar process |
US5911692A (en) * | 1998-01-20 | 1999-06-15 | General Electric Company | Sparse two-dimensional wideband ultrasound transducer arrays |
US6312384B1 (en) * | 1998-03-31 | 2001-11-06 | General Electric Company | Method and apparatus for flow imaging using golay codes |
US5902243A (en) * | 1998-04-15 | 1999-05-11 | Acuson Corporation | Ultrasonic imaging method with multiple pulse cancellation |
US6014897A (en) * | 1998-09-02 | 2000-01-18 | Mo; Larry Y. L. | Method and apparatus for improving sidelobe performance of sparse array using harmonic imaging |
US6142946A (en) * | 1998-11-20 | 2000-11-07 | Atl Ultrasound, Inc. | Ultrasonic diagnostic imaging system with cordless scanheads |
US6135963A (en) * | 1998-12-07 | 2000-10-24 | General Electric Company | Imaging system with transmit apodization using pulse width variation |
US6432056B1 (en) * | 1999-10-08 | 2002-08-13 | Atl Ultrasound | Ultrasonic diagnostic imaging system with high frame rate synthetic transmit focus |
US6315723B1 (en) * | 1999-10-08 | 2001-11-13 | Atl Ultrasound | Ultrasonic diagnostic imaging system with synthesized transmit focus |
US6282963B1 (en) * | 1999-10-12 | 2001-09-04 | General Electric Company | Numerical optimization of ultrasound beam path |
WO2001057514A1 (en) * | 2000-01-31 | 2001-08-09 | Angelsen Bjoern A J | Correction of phasefront aberrations and pulse reverberations in medical ultrasound imaging |
US6658141B1 (en) * | 2000-02-25 | 2003-12-02 | Medison Co., Ltd. | Filtering method and apparatus for improving resolution of ultrasound image |
JP2001238883A (ja) * | 2000-02-29 | 2001-09-04 | Aloka Co Ltd | 超音波診断装置 |
US6487433B2 (en) * | 2001-01-08 | 2002-11-26 | General Electric Company | Method and apparatus using golay-coded excitation for echocardiology |
KR100393370B1 (ko) * | 2001-04-25 | 2003-07-31 | 주식회사 메디슨 | 직교 골레이 코드를 이용하는 초음파 영상 형성 방법 및장치 |
KR100406097B1 (ko) * | 2001-12-26 | 2003-11-14 | 주식회사 메디슨 | 가중 쳐프 신호를 이용한 초음파 영상 형성 장치 및 방법 |
KR100406098B1 (ko) * | 2001-12-26 | 2003-11-14 | 주식회사 메디슨 | 가중된 직교 쳐프 신호를 이용한 동시 다중 송신 집속기반의 초음파 영상 형성 장치 및 방법 |
KR100419806B1 (ko) * | 2001-12-31 | 2004-02-21 | 주식회사 메디슨 | 평면파를 이용하는 초음파 영상의 합성 구경 집속 방법 |
US6582369B1 (en) * | 2002-01-02 | 2003-06-24 | Computed Ultrasound Global Corporation | Method for dynamic focus control |
US6905465B2 (en) * | 2002-04-05 | 2005-06-14 | Angelsen Bjoern A. J. | Corrections for pulse reverberations and phasefront aberrations in ultrasound imaging |
CA2513446C (en) * | 2003-01-14 | 2014-08-26 | University Of Virginia Patent Foundation | Ultrasound imaging beam-former apparatus and method |
JP4235006B2 (ja) * | 2003-02-12 | 2009-03-04 | アロカ株式会社 | 超音波診断装置 |
US8038616B2 (en) * | 2003-05-30 | 2011-10-18 | Surf Technology As | Acoustic imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties |
US20050277835A1 (en) * | 2003-05-30 | 2005-12-15 | Angelsen Bjorn A | Ultrasound imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties |
WO2005008280A1 (en) * | 2003-07-17 | 2005-01-27 | Angelsen Bjoern A J | Corrections for wavefront aberrations in ultrasound imaging |
US7889787B2 (en) * | 2003-08-04 | 2011-02-15 | Supertex, Inc. | Ultrasound transmit beamformer integrated circuit and method |
US9116229B2 (en) * | 2003-08-04 | 2015-08-25 | Microchip Technology Inc. | Ultrasound transmit beamformer integrated circuit and method |
US20050096544A1 (en) * | 2003-10-30 | 2005-05-05 | Xiaohui Hao | Method and apparatus for single transmission Golay coded excitation |
US20060074320A1 (en) * | 2004-08-27 | 2006-04-06 | Yoo Yang M | Home ultrasound system |
US9451932B2 (en) * | 2004-12-30 | 2016-09-27 | Crystalview Medical Imaging Limited | Clutter suppression in ultrasonic imaging systems |
JP5341352B2 (ja) * | 2004-12-30 | 2013-11-13 | クリスタルビュー メディカル イメージング リミテッド | 超音波画像化システムのクラッター抑制本出願は、2004年12月30日に出願された米国仮特許出願No.60/640,368に基づき優先権の利益を主張する。本出願は、2004年1月7日に出願した米国仮特許出願No.60/534,390に関連し、明細はここに参照のため組み込まれている。 |
ATE410752T1 (de) * | 2005-03-17 | 2008-10-15 | Koninkl Philips Electronics Nv | Verfahren und einrichtung zur iterativen rekonstruktion von herzbildern |
WO2007007414A1 (ja) * | 2005-07-14 | 2007-01-18 | Rion Co., Ltd. | 遅延和型センサアレイ装置 |
WO2007027703A2 (en) * | 2005-08-29 | 2007-03-08 | University Of Toledo | System for extended high frame rate imaging with limited-diffraction beams |
US8761477B2 (en) * | 2005-09-19 | 2014-06-24 | University Of Virginia Patent Foundation | Systems and method for adaptive beamforming for image reconstruction and/or target/source localization |
US7324042B2 (en) * | 2005-11-15 | 2008-01-29 | The Boeing Company | Monostatic radar beam optimization |
US7876973B2 (en) * | 2006-01-12 | 2011-01-25 | Integrity Applications Incorporated | Edge ringing artifact suppression methods and apparatuses |
US7874991B2 (en) * | 2006-06-23 | 2011-01-25 | Teratech Corporation | Ultrasound 3D imaging system |
US8254654B2 (en) * | 2007-10-31 | 2012-08-28 | University Of Southern California | Sidelobe suppression in ultrasound imaging using dual apodization with cross-correlation |
US8744155B2 (en) * | 2008-02-16 | 2014-06-03 | University Of Virginia Patent Foundation | Imaging or communications system utilizing multisample apodization and method |
US10080544B2 (en) * | 2008-09-15 | 2018-09-25 | Teratech Corporation | Ultrasound 3D imaging system |
WO2011057252A1 (en) | 2009-11-09 | 2011-05-12 | Sonosite, Inc. | Systems and methods for beam enhancement |
-
2010
- 2010-11-09 WO PCT/US2010/056000 patent/WO2011057252A1/en active Application Filing
- 2010-11-09 JP JP2012538082A patent/JP5973349B2/ja active Active
- 2010-11-09 BR BR112012010958-8A patent/BR112012010958B1/pt not_active IP Right Cessation
- 2010-11-09 CN CN201080060936.3A patent/CN102753104B/zh not_active Expired - Fee Related
- 2010-11-09 US US13/508,966 patent/US8876719B2/en active Active
- 2010-11-09 EP EP10829276.4A patent/EP2498683B1/en not_active Not-in-force
-
2014
- 2014-08-22 US US14/466,902 patent/US9420998B2/en active Active
-
2016
- 2016-07-14 JP JP2016139832A patent/JP6227724B2/ja active Active
- 2016-08-22 US US15/243,878 patent/US20160354062A1/en not_active Abandoned
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1708698A (zh) * | 2002-11-08 | 2005-12-14 | 皇家飞利浦电子股份有限公司 | 用于自动设置超声换能器阵列的发射孔径和旁瓣缩减的方法和设备 |
Also Published As
Publication number | Publication date |
---|---|
JP2013509971A (ja) | 2013-03-21 |
US9420998B2 (en) | 2016-08-23 |
BR112012010958A2 (pt) | 2017-11-07 |
US8876719B2 (en) | 2014-11-04 |
BR112012010958B1 (pt) | 2021-09-08 |
JP2016179313A (ja) | 2016-10-13 |
US20120289835A1 (en) | 2012-11-15 |
US20160354062A1 (en) | 2016-12-08 |
EP2498683A4 (en) | 2014-10-22 |
US20140363068A1 (en) | 2014-12-11 |
CN102753104A (zh) | 2012-10-24 |
WO2011057252A1 (en) | 2011-05-12 |
JP5973349B2 (ja) | 2016-08-23 |
JP6227724B2 (ja) | 2017-11-08 |
EP2498683B1 (en) | 2017-04-12 |
EP2498683A1 (en) | 2012-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102753104B (zh) | 增强波束的系统和方法 | |
US9198636B2 (en) | Continuous transmit focusing method and apparatus for ultrasound imaging system | |
JP2807113B2 (ja) | 音響走査方法及び装置 | |
US8690781B2 (en) | Coherent image formation for dynamic transmit beamformation | |
US7722541B2 (en) | Multiple receive beams for rapid acquisition | |
JP2006204923A (ja) | コヒーレンス度適応超音波撮像 | |
EP1782094A2 (en) | Ultrasound imaging using non-linear manipulation of forward propagation properties of a pulse | |
WO1996003918A1 (en) | Method and apparatus for adjustable frequency scanning in ultrasound imaging | |
US6368276B1 (en) | Deep penetration beamformed television | |
Chang et al. | A new synthetic aperture focusing method to suppress the diffraction of ultrasound | |
KR101911470B1 (ko) | 평면파를 이용하여 3차원 초음파 영상을 생성하는 장치 및 방법 | |
EP3596493B1 (en) | Methods and systems for filtering ultrasound image clutter | |
Kortbek et al. | P2b-1 synthetic aperture focusing applied to imaging using a rotating single element transducer | |
Schou et al. | Synthetic aperture sequential beamforming using spatial matched filtering | |
Kim et al. | New synthetic aperture imaging technique with dynamic apodization window | |
Austeng et al. | Synthetic aperture composite focusing for high framerate medical ultrasound imaging | |
Zhang et al. | A Spatial Composite Denoising Algorithm Based on Synthetic Aperture Continuous Beamforming |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160309 |