专利名称：：多阶段数字超声波束形成器的制作方法
技术领域：
：0001本发明总体而言涉及诊断超声成像，具体而言涉及数字波束形成。
背景技术：
：0002诊断超声成像的重要性受到广泛地认可，并且随着成像分辨率以及可获得用途和特征的范围稳步地增加而发展。曾经只有在设施最完备的医院才可获得的昂贵奢侈品，而现在诊断超声成像即便在一些私人诊所也是普通且几乎是常规提供的手段。0003医学诊断超声机器是民用领域最复杂的信号处理系统之一，并且是少数能迫使现代处理器达到其性能极限中的一个。因此一幅超声图像的质量直接受到许多因素的影响，这些因素包括用于产生被发射进入身体所需模式的超声信号并用于接收和处理这些超声信号的回波返回所涉及的硬件以及所使用的例程两者。0004图1示出了一种标准超声成像系统的主要部件。用户将各种常规的扫描参数录入输入/输出(I/O)单元100，该元件可以包括常规的硬件，系统软件，以及诸如键盘、旋钮、鼠标和域按钮的设备。将输入元件连接到处理系统200，所述处理系统200通常是一组电连接且协调工作的、带有必要系统软件(例如操作系统)的、诸如微处理器和数字信号处理器的处理器210。因此，1/O单元100可以不需要本身拥有任何处理能力，而是简单的充当主处理系统200的外围设备。0005正如在已知系统中那样，处理系统200设置、调节和监测常规发射控制电路220的操作参数。这个控制电路220产生电控制和驱动信号并将其施加于超声探头，即，换能器300，其包括由电声(通常是压电的，例如PZT)元件阵列310。如本领域所熟知的，在向压电元件施加适当频率的电信号时该压电元件产生超声波。0006为了产生被发射进入身体("询问(interrogation)区域"或"感兴趣区域"400)的超声信号，以给定频率来电激励这些元件的全部或一些子组并将其相对于每个元件在阵列中的位置以如下方式单独地调整相位和时间移位，即使所述信号的同相部分形成相干发射(Tx)波束:选择时间延迟以便相干地叠加来自期望角度方向的回波信号，而那些来自其他方向的信号并不进行相干叠加并因此趋向于相互抵消。通常可以改变(变迹(apodize))向阵列元件所施加的激励信号的振幅，以便控制得到的波束宽度并抑制旁瓣。同样可以改变阵列的孔径，也就是说，可以改变激励元件(通常并不是所有元件被同时激励)的"宽度"，如同可以改变照相机的光圈以便改变焦深并提供一致的横向分辨率。0007根据特定的信号包络激励元件以使对给定区域的成像最佳的各种己知方法一般称之为发射波束形成。通过改变发射信号的振幅和相位调整，发射波束不仅能够聚焦于特定深度，而且能够进行偏轴转向，即发射波束不与那些激励元件的中心元件垂直。例如，在图1中，发散波束己经以相对于阵列310表面的一角度进行转向。各个发射信号的异相部分将趋向于相互抵消。发射波束通常在焦深处会聚，发射波束超过焦深后再次发散。发射波束在方位角順向方向LAT和高度方向EL上进行转向，并在深度/轴向方向AX上进行聚焦，以便将波束的超声能量集中在询问区域400内的期望点上，例如结构410上。0008在询问区域中的任意反射结构被声穿透之后，通常通过切换电路230将换能器从发散模式切换到接收模式，以便从成像区域内任意组织410反射回到阵列表面上的超声将使各元件产生相应的电信号。然而，任意特定元件感测到的超声将是从成像区域内任意位置反射回来的所有超声能量的合并，而不管能量源自那些元件。在一般情况下，每个元件都将接收所有受激励形成发射波束的元件所产生的超声能量中的一些。在接收(Rx)模式中，每个元件将因此产生时变信号，而该时变信号实际上将不同于所有其他元件的接收信号。此外，仍然是反射能量的同相部分占优势。问题是如何最好地确定这些相干部分的本来面目，也就是，在不同元件上的哪些感测值代表来自期望焦点的返回信号，而带有适当时间延迟的返回信号对应于在其他阵列元件上的返回信号的部分。这一过程称之为接收波束形成，并且接收波束形成器的首要任务就是确定由所有有效Rx阵列元件所接收的来自靶的相干信号并对所述相干信号进行求和。0009众所周知，压电阵列元件310将由回波信号引起的微小机械振动转换成相应的射频(RF)电信号。然后由接收控制器240向返回信号施加放大和其他常规的信号调整。这种处理根据识别回波信号和形成接收波束的需要包括多个步骤，所述回波信号对应于在询问区域中的每个扫描元件。接收控制器240同样将超声射频(RF)返回信号(典型为几兆赫到几十兆赫的数量级)转换成用于处理的低频范围，通常将所有或部分的接收控制器240集成到处理系统200自身中。这是超声成像领域所熟知的。接收控制电路包括接收波束形成器，下文将更加详细地描述该接收波束形成器的各方面。需要指出的是，所述发射和接收控制器(就此而言，或者为发射和接收波束形成器)因此不一定实现为分立元件，而是可以相对于两者的硬件和软件整体或部分地进行组合。0010在常规的B模式扫描中，询问区域内的每个点随后表示为强度(亮度)值。询问区域因此能够表示为亮度值或信号强度值的离散模式(矩阵)，所述亮度值或信号强度值作为帧数据存储在存储器250中。其他扫描模式，例如对于多普勒成像，产生并存储相应的值。0011询问区域通常具有与用户期望看到的所显示内容不相同的形状；即便这时，构成帧数据的数字声学强度值通常并不具有适于直接驱动常规灰阶或彩色显示的形式。然后通常向常规扫描转换器260施加所选子集(扫描平面)的声学强度值，所述扫描转换器260将数字声学值转换成适于在驱动显示装置120过程中使用的显示强度值，所述显示装置120可以是I/O单元100的一部分。显示装置120通常包括被分成像元或"像素"的X-Y(或极)矩阵或模式的屏幕(例如，LED或CRT)，所述像元或"像素"构成用户可以观看和解释的图像。0012图2示出了接收波束形成器500的主要部件。所示的波束形成器500出于图示目的进行了简化，并且大小和距离并未按比例绘制。例如，即便是一维线性阵列通常具有数十、甚至数百个阵列元件，而不是图2中所示的八个。而且，通常不是同时激励所有的超声元件，而是根据某一模式激励所有超声元件中的预定子组。在图2中，假设发射波束己经聚焦于焦点FP上；从点500回到阵列元件的线代表声返回信号，也就是，发射超声的回波。0013在图示的情况下，仅仅为了易于图示，焦点FP位于垂直于最上面阵列元件的表面。假设超声通过至少大体上声学均匀的组织进行传播,从焦点FP反射回来的发射波束的超声能量将在到达其他阵列元件前到达最上面的阵列元件。虚线520示出了朝向阵列元件310返回的相干超声的行程；也就是，线520示出了发射波束的反射的相干部分。0014正如波形线530图示的，反射的超声将使每个元件产生相应的模拟电信号。在波形线上的点代表与来自阵列元件中相应元件的接收数据值对应的瞬时值。通常在通过一组放大器540放大后，来自电声阵列元件310的模拟信号被模数转换器(ADC)550转换成数字值，所述模数转换器550的采样率受到时钟电路(CLK)555的控制。0015通常接着将时间序列的经转换数字值缓存在先进先出(FIFO)锁存器560中(为了速度，虽然它们可以被存储在其他类型的存储器设备中)。在图2中，每个存储的数字值表示为锁存器560内信号通道上的小方块。在锁存器内的实心(填满的)方块表示与来自同一靶的回波相对应的值。0016然后在处理模块(570)中对所述数字值进行处理，所述处理过程通常涉及运行计算机可读软件代码的主处理器210(图1)，所述计算机可读软件代码不仅定义了所实现的波束形成例程，而且定义了这样的程序，所述程序涉及更大处理器的专用处理器部分，以便形成代表分配给焦点的值的输出信号，例如亮度值。所述处理过程通常包括如下这样的步骤-时间选通、增益补偿、衍射补偿、动态接收延迟、动态接收孔径大小控制和动态接收变迹。0017由于其在超声成像中的重要性，因此在产业上一直努力改善接收波束形成器。事实上，有许多关于该主题的科技论文、甚至授权的专利。这些授权的专利包括如下授权的美国专利-美国专利号题目7,199,7386,937,1766,716,1716,676,6026,635，0196,602,1946,592,5246,565,5106,500,1206,491,6346,447,4526,436,0486,436,0446,363,0336,245,0166,172,9396,155細6,146,3286,123,6716,113,5456,110,116SigmadeltabeamformerandmethodwithreducedartifactUltrasonicsignalacquisitioninthedigitalbeamformerSystemandmethodforinterfacinganultrasoundtransducerwithacomputingdeviceperformingbeamformingprocessingTwodimensionalarrayswitchingforbeamforminginavolumeScanheadassemblyforultrasonicimaginghavinganintegralbeamformeranddemountablearrayDualbeamformerultrasoundsystemfor2Dand3DimagingTransmitbeamformerdelayarchitectureandmethodfordiagnosticmedicalultrasoundMethodandapparatusforservicingremoteultrasoundbeamformerfromcentralservicefacilityBeamformingsystemusinganalograndomaccessmemorySub-beamformingapparatusandmethodforaportableultrasoundimagingsystemSystemandmethodforreductionofparallelbeamformingartifactsUltrasonicdiagnosticimagingsystemwithscanheadelevationbeamformingSystemandmethodforadaptivebeamformerar)odizationMethodandapparatusfortransmitbeamformersystemUltrasoundimagingsystemhavingpost-beamformersignalprocessingusingdeconvolutionalgorithmMethodandapparatusfortransmitbeamformersystemUltrasonicimagingsystemwithbeamformingusingunipolarorbipolarcodedexcitationMethodandapparatusforultrasonicbeamformingusinggolay-codedexcitationMethodandapparatusfordistributed,agilecalculationofbeamformingtimedelaysandapodizationvaluesUltrasonicbeamformingwithimprovedsignal-to-noiseratiousingorthogonalcomplementarysetsMethodandapparatusforreceivebeamformersystem6,108,2736,104,6736,042,5476,029,1165,995,4505,984,8695,928,1525,921,9325.882.3075,856，9555.851.1875,844,1395.827.1885,797,8475,720,2895,696,7375.685.3085,676,1475,675,5545,653,2365,617,8625細扁TransmitbeamformerwithfrequencydependentfocusMethodandapparatusfortransmitbeamformersystemMethodandapparatusforreceivebeamformersystemMethodandapparatusforabasebandprocessorofareceivebeamformersystemMethodandapparatusfortransmitbeamformersystemMethodandapparatusforultrasonicbeamformingusinggolay-codedexcitationMethodandapparatusforabasebandprocessorofareceivebeamformersystemMethodandapparatusforabasebandprocessorofareceivebeamformersystemMethodandapparatusforreceivebeamformersystemMethodandapparatusfortransmitbeamformersystemMethodandapparatusforultrasonicbeamformingwithspatiallyencodedtransmitsMethodandapparatusforprovidingdynamicallyvariabletimedelaysforultrasoundbeamformerMethodandapparatusforreceivebeamformersystemMethodandapparatusforcomplexbandpassfilteringanddecimationinultrasoundbeamformerMethodandapparatusforageometricaberrationtransforminanadaptivefocusingultrasoundbeamformersystemTransmitbeamformerwithfrequencydependentfocusMethodandapparatusforreceivebeamformersystemUltrasonicreceivebeamformerwithphasedsub-arravsMethodandapparatusfortransmitbeamformerApparatusforreal-timedistributedcomputationofbeamformingdelaysinultrasoundimagingsystemMethodandapparatusforbeamformersystemwithvariableapertureTransmitbeamformerwithfrequencydependentfocus5,581,517Methodandapparatusforfocuscontroloftransmitandreceivebeamformersystems5,573,001Ultrasonicreceivebeamformerwithphasedsub-arrays5,570,691Methodandapparatusforreal-time,concurrentadaptivefocusinginanultrasoundbeamformerimagingsystem5,566,675Beamformerforphaseaberrationcorrection5,555,534Methodandapparatusfordopplerreceivebeamformersystem5,551,433Methodandapparatusforageometricaberrationtransforminanadaptivefocusingultrasoundbeamformersystem5,544,128Multi-beamdigitalbeamformingmethodandapparatus5,517,537Integratedacousticleakdetectionbeamformingsystem5,388,461Beamformingtimedelaycorrectionforamulti-elementarrayultrasonicscannerusingbeamsum-channelcorrelation5,388,079Partialbeamforming5,369,624Digitalbeamformerhavingmulti-phaseparallelprocessing4,644,795Highresolutionmultilineultrasonicbeamformer4,542,653Apparatusandmethodforbeamforminginanultrasonictransducerarray0018常规的数字超声波束形成器通常产生一个发射波束，而接收多个返回的接收波束，这是因为单个发射波束可以被不同结构反射回来。常规的单级波束形成器通常能够通过使用动态接收延迟、动态接收孔径尺寸控制和动态接收变迹来处理这些多个波束。即便如此，现有技术的波束形成器不能考虑从多次发射上的多个接收波束获得的特定信息。
发明内容0019本发明提供一种超声成像系统，包括主处理系统和超声探头，所述主处理系统包括至少一个处理器和至少一个存储器，而所述超声探头包括电声元件阵列。所述阵列可以是一维的(直线或非直线，即弯曲的)或具有更高维，例如二维的。波束形成系统包括发射(Tx)波束形成器和接收(Rx)波束形成器。将发射波束形成器设置成激励阵列中的至少一个。接收波束形成器包括第一阶段，对于每个发射波束，将所述第一阶段设置成确定多个一次(primary)接收波束。接收波束形成器还包括第二阶段(其可与第一阶段一起实现为集成的波束形成器，或者分别地实现)，将所述第二阶段设置成根据对至少两个一次接收波束的时间对准合并来计算至少一个二次接收波束。0020根据相关的操作方法，激励电声元件阵列中的至少一个元件以产生一系列声穿透一个感兴趣区域的发射波束。然后对于至少一个所述发射波束，确定多个一次接收波束。然后，波束形成器的第二阶段根据对至少两个所述一次接收波束的时间对准合并来计算至少一个二次接收波束。然后根据至少一个二次接收波束生成图像。0021波束形成器的第二阶段可以时间对准至少两个一次接收波束，所述至少两个一次接收波束对应于来自不同的发射波束的超声返回信号、但却在共同的接收波束线(通过焦点和在阵列上相应的波束起点的超声传播线)上进行接收，并且然后根据所述至少两个时间对准的一次接收波束的函数来计算每个二次接收波束。一个适合函数的示例是可对至少两个一次接收波束的加权、延迟求和。一个时间对准的示例包括计算至少一个一次接收波束相对于另一个一次波束的时间延迟。0022一个计算时间延迟方法的示例包括模拟与第一原始波束的预期发射条件相对应的条件；对于指定时间的所述条件，确定第一原始波束上的发射脉冲的第一瞬时发射脉冲压力场(700)，其对应于指定深度；确定沿着从第一原始波束的第一波束起点(A)的半径方向的发射脉冲的中心；确定第一发射起点(A)与指定深度下的瞬时发射脉冲压力场和第一原始波束的交叉点(B)之间的第一发射时间(RO);确定第二发射起点(C)与对于第二个一次接收波束的第二脉冲压力场的交叉点(D)之间的第二发射时间(Rl);根据预先确定的关于第一发射时间和第二发射时间的函数确定时间差分；以及在计算至少一个二次接收波束的计算过程中，向第二个一次接收波束施加根据预先确定的关于时间差分的函数所计算出的延迟。0023图1示出了现有技术已有的一种超声成像系统的主要部件；0024图2示出了现有技术已有的一种发射和一次接收波束形成器的操作；0025图3示出了一维阵列中发射波束的分布、具体而言是用于说明本发明原理的编号规则；0026图4a-4j示出了如何根据一系列发射波束形成一次接收波束的二次序列的一个示例；0027图5和图6是如何对准来自不同发射的一次接收波束以便形成用于产生次级接收波束的两个不同表示；0028图7示出了用于形成次级接收波束的延迟值的计算；0029图8示出了一种包括两级波束形成的超声成像系统的各部件；0030图9a-9c示出了用于二维阵列的两级波束形成。具体实施例方式0031根据本发明的多级波束形成器及其新颖的操作方法主要参考一维阵列(其可以是直线或弯曲的，或者是更高维阵列的各元件中的一维子组，等)进行描述，但是本发明的原理并不限制于此，并且事实上图示出2D示例。本发明可以用于具有任意数目元件的任何阵列，只要阵列中有足够的元件允许接收元件/波束进行如下文解释的"重叠"。0032为了更好地理解本发明的原理，假设本发明将使用N个元件的线阵来实现，如图3中所示。发射(Tx)波束可以源于某一元件的中心或者两个元件之间。仅仅为了说明目的，假定将起点在阵列元件(m-l)和m之间的Tx波束编号为(2m-2);将其起点在元件m中心的波束编号为(2m-l)而将其起点在阵列元件m和(m+l)之间的波束被编号为2m。然后将起点在阵列的左右边缘的波束(如在图3中所示)分别标记为0和2N。这样，其起点基本上在元件4和5之间的发射波束为波束8;其起点基本上在元件2的中心的发射波束为Tx3;等等。因此，对于N个元件的线阵，通常会有(2N+1)个可能的波束起点，例如257个波束起点对应128个元件的阵列。0033需要指出的是，"起点"并不意味着发射波束必须从该点垂直于阵列表面进行转向；因而，许多发射波束可以源于相同的点，但可能转向到不同的角度和深度。同样，波束的起点位于特定点并不一定限定这些个阵列元件用于产生所述波束。例如，Tx5可能是通过单独激励阵列元件3产生的，或者通过激励元件2-4、或仅是2和4，或者元件l-5、或元件l、3和5，等等产生，所有这些都表示对称的元件组合，即这些元件组合以元件3的中心作为所述组合的中心。然是，甚至这种对称性也不是必须的使用已知方法同样能够计算延迟以便特定的Tx波束具有并不位于激励元件组的中心的起点，艮卩，所述Tx波束来自于非对称模式的元件。例如，当使用元件1-15来形成发射波束时，同样能够确定Tx8。0034如何产生特定的发射波束通常将取决于像期望的孔径和变迹这样的考虑。本发明不要求任何特定的孔径、变迹等，只要能够确定出Tx平均起点。另一种考虑波束起点的方法是，起点就是系统计算孔径中所有元件的延迟的那个位置。波束是发自波束起点并穿过焦点FP的线。这对于发射和接收波束都是正确的，尽管声音传播的方向显然会被反转(与从焦点出来的方向相反)。在发射和接收波束形成中重要的计算是延迟的计算，因而波束起点是延迟为零的点，这意味着相对于它来测量其他延迟。然后可以将延迟长度定义为在元件到焦点间长度和波束起点到焦点间长度之间的差。0035在图3中所示的编号规则同样用于对接收波束进行编号。需要指出的是，相同的元件、甚至是相同编号的元件并不一定用于接收产生发射波束所使用的回波信号。这又再次成为应用哪个孔径等的问题。作为一个示例，元件1-4可能受到激励以形成Tx4，但是然后可能使用在元件3-5处接收到的信号来形成在线7上具有零延迟或"起点"的接收波束；这样，系统可能于Tx4上发射，但在Rx7上"收听(listen)"。0036图4a-4j示出了如何从两个不同的发射波束Tx8和Txl2形成不同的一次接收波束。仅仅通过图示来看，发射波束被显示为由四个相邻元件产生元件3-6对应于Tx8而元件5-8对应于Txl2。类似地，仅仅是举例，假设在接收模式中使用来自三个相邻元件的信号。为简单起见，同样假设发射线和接收线位于它们各自有效元件组的中心。如上面所解释的，这些假定仅仅反映出若干可能选择之一。发射中所使用的元件用粗体轮廓线示出。用于接收的元件通过其上的括号、还有对其下面起点编号的识别进行指示。这样，在图4a中，Tx8由元件3-6产生并且在元件5-7上进行接收，以形成Rxll。0037在本发明中，由来自单次发射的回波信号所形成的波束被称为一次接收波束。符号PB(i,J)在此用于指示从其起点为i的发射波束产生的并且由其接收波束起点在位置j处的任意选定元件组接收(即，在线j上进行接收)的一次接收波束。0038在图4a-4e中，发射波束在每种情形中都是Tx8，但是用于从同一次发射Tx8形成五个并行接收波束的元件，从用于图4b所示形成PB(8，3)的元件l-3(Rx3)，移位到用于图4c所示形成PB(8，5)的元件2-4(Rx5)，移位到用于图4d所示形成PB(8，7)的元件3-5(Rx7)，移位到图4e所示形成PB(8，9)的元件4-6(Rx9)，移位到图4a所示形成PB(8，11)的元件5-7(Rxll)。这样，在接收模式中的移位是每个接收波束一个元件(两个波束号)。0039在图4f-4j中，发射波束是Txl2。相对于在图4a-4e中所示的情况，发射波束的起点因而已经向右移位了四个波束号(两个元件)。需要指出的是，即使对于相同的耙，由于不同的Tx波束，来自该相同靶的回波将位于原始波束的不同深度，例如，在线PB(8，5)和PB(12，5)上。根据本发明的第二级波束形成器(后面将对此做出详细描述)消除或至少基本上减小该差。在图4f中，对于Txl2在接收模式中使用的元件是用于形成PB(l2，5)的元件2-4;在图4g中，是形成PB(lW)的元件3-5(Rx7);在图4h中，是形成PB(l2，9)的元件4-6(Rx9);在图4i中，是形成PB(12，11)的元件5-7(Rxll);在图4j中，是形成PB(l2J3)的元件5-7(Rx13)。这样，在图4f-4j中随后的每一个一次接收波束都是通过相对前一个一次接收波束向右移位两个波束号(一个元件)产生的。0040图5示出了部分超声扫描的许多可能示例之一。在图5中，编号的圆圈代表一次波束。在该所示的情况中，对于Tx波束计算(以任何已知方式)八个一次波束，然后将孔径向右(如在图5中所示)"移位"四个波束号，形成新的发射波束并且计算新的接收波束，然后将孔径再向右移位四个波束号，等等。图5示出了对于五个发射波束的接收模式中的激励模式。再次仅仅是举例，假设五个发射波束Tx8、Txl2、Txl6、Tx20和Tx24。在图5中，深色圆圈代表对于相应每一次发射所计算的一次接收波束。这样，第一发射波束Tx8产生一次接收波束PB(8，1)、PB(8，3)、PB(8，5)、PB(8，7)、PB(8，9)、PB(8，11)、PB(8，13)和PB(8，15)。以表格的形式，针于五个发射波束中每一个的一次接收波束是-Tx棘号接收波束号81、3、5、7、9、11、13、15125、7、9、11、13、15、17、19169、11、13、15、17、19、21、232013、15、17、19、21、23、25、272417、19、21、23、25、27、29、31,.■表l0041这同样在图6中示出，图6示出并标记出每个发射波束的各个接收波束，其中发射波束在时间上间隔一量At，At将是选定帧率的函数。如图5、图6和表1所示，若干接收波束包括与相同声穿透结构有关的、但是来自不同发射波束的信息。因而，在这一简单示例中，一次接收波束13、15、17和19有四个不同的值。例如，对于一次接收波束13，有值PB(8，13)、PB(12，13)、PB(16，13)和PB(16，13)，而对于接收波束19，有值PB(12，19)、PB(16，19)、PB(2(U9)和PB(24，19)。0042依据波束号的定义，如果在接收孔径中使用与n个波束号对应的阵列元件，并且所述孔径在每次发射时移位m个波束，那么除了在边缘附近，每个波束位置将有n/m个一次波束。n和m的值通常是代表分辨率与速度的权衡的设计选择。需要指出的是，能够依据波束号或元件号来定义移位操作，并且编号操作在任意一种情况中都仅是规则的选择问题。0043一次波束的形成在此被称为第一或一次波束形成阶段，其在图6中标记为第一阶段。本发明具有使用了常规波束形成器忽视的信息的二汰波束形成。该二次波束形成阶段提供对相同线上各一次波束(如图所示，具有相同的波束号)的动态波束形成。0044本发明根据线p上各一次接收波束的函数而计算(使用存储在或在任何适合的媒介，例如存储器250内的可执行编码)的二皿收波束SB(p)。一个适合的函数是对一次接收波束PB(c，p)的时延加权求和，即对其中，所述和是对p值的c(全部或选定子集)上获得。在此，PB(c,p,t-《)是一次接收波束PB(c,p)时延量《，并用值W。进行加权(可能是简单地统一加权，意味着根本不进行微分加权)。换句话说，将对于不同次发射、但在相同束上接收的至少一些一次接收波束的值进行合并。需要指出的是，所述求和是在所述波束的整个长度上进行的，也就是，在选定区间内的所有时间t上进行求和。而且，W。和《同样可以是时间(深度)相关的，并且可能对于不同波束具有不同的值。0045因而，返回到上述示例，可以根据对PB(8，13)、PB(12，13)、PB(16，13)和PB(16，13)的时延加权求和来计算SB(13)。在实践中，这意味着每个二次接收波束包括来自若干发射波束的能量(信息)。因此，对孔径、移位、发射定时(所有这些可以视具体情况使用已知的运算方法而确定)、和本发明的应用的正确选择应当得到合成的二次接收波束，其具有更少的噪声，并因此与单独处理的各一次接收波束中的任意一个相比都更清楚。0046作为另一示例，如前面一样假设，有八个相对间隔一个元件(两个波束位置)的并行Rx波束，但是在这种情况下，假设每次发射的移位仅仅是两个波束位置(一个元件)而不是四个波束位置(两个元件)。表2示出了在发射波束号和接收波束号间的对应关系<table>tableseeoriginaldocumentpage17</column></row><table><table>tableseeoriginaldocumentpage18</column></row><table>表0047在表2中可以看到从波束号15并向上，对于每个波束位置有八个一次波束。即，为了形成最终结果，二次波束号15、SB(15)由一次波束PB(8，15)、PB(10，15)、PB(12，15)、PB(14，15)、PB(16，15)、PB(18，15)、PB(20,15)和PB(22，15)形成。如图前面一样，所有这八个一次波束将迸行时延并加权求和，以便与假设的发射起点和接收起点都位于波束15的波束起点的波束对准。0048在表2中，同样观察到，对于波束号9而言，每个波束位置有5个一次波束。为了形成最后的波束号SB(9)，可以对一次波束PB(8，9)、PB(10，9)、PB(12，9)、PB(14，9)和PB(16，9)进行时延和加权求和。换句话说，很可能对于不同的二次波束在求和中有不同的项数；这通常发生在接近阵列边缘时。同样可能简单到在阵列边缘附近不进行二次波束形成，只要这样做优于增加额外的编程负担。一旦计算完二次波束，可以将其如同在常规波速形成器中的一次波束一样作进一步的处理和显示。0049可以使用已知的理论和实验方法来确定权重W。，假设恰好都进行加权。例如，人们可以选择权重以便削弱离轴更远的一次波束。作为另一示例，人们可以选择削弱在时间上相对更靠后发生的一次波束，具体而言其中例如对于心脏成像而言，与那些更静态的对象、例如腹部结构相反，图像区域是自身运动的。这些考虑事项同样在确定有效阵列元件的移位程度时应当受到重视——此外，累加的一次波束越多，二次波束有可能越好，但是这将占用的时间越长，并且对于诸如心脏的动态区域通常优选更快的扫描。在实践中，一次激励的元件数目和移位程度将是要权衡的一部分期望的分辨率越高，进行处理将占用的时间越长。0050二次接收波束的计算中的一部分包括确定动态时延《，所述动态时延《用于在SB(p)计算中时间上正确对准各一次波束，以便确保它们真实地代表相同的波束线。在常规发射波束形成器的延迟计算中，通常根据几何结构的函数计算发射时间一Tx脉冲从Tx元件到Tx焦点的传播时间，其等于从Tx元件到Tx焦点的直线距离除以声速。然而，这要假设，声波脉冲前沿是球形的，并且Tx元件位于中心；或者，这些系统假设Tx元件是点源。0051在根据本发明的两阶段波束形成器的第二阶段中，所述计算不是基于来自个别元件或到达个别元件的信号，而是基于所述原始发射波束。在这些延迟计算中，因此，人们应当同样考虑发射时间，所述发射时间为Tx脉冲从原始发射波束的Tx起点到Tx焦点的传播时间。然而，这一传播时间并不一定正比于从所述发射波束的Tx起点到Tx焦点的直线距离除以声速，这是因为所述原始发射波束的波前不再必须是球形的。一个特殊的情况是其中所有原始发射波束的发射孔径是单个探头元件；在这一特殊情况下，单独的几何计算将足够用来计算发射延迟。但是，在一般情况下，原始发射波束的发射孔径是一组具有延迟分布和变迹的探头元件；这样的孔径不能被看作是一个点源。0052一种计算一次波束的Tx延迟《的方法是通过在给出有关计算将要用于哪种类型的身体区域(心脏、腹部、胎儿、隔塞性脂肪组织等)，即超声将通过什么样的介质的知识的情况下，模拟发射波束的瞬时Tx脉冲压力场。这种模拟的一个示例的主要步骤包括0053l.设定与一次波束的实际Tx条件相匹配的模拟条件(通过计算或者通过使用物理模型)。所述条件应当包括探头的几何形状、Tx脉冲、Tx孔径大小、Tx变迹、探头元件方向性、组织衰减、Tx焦点位置和其他。00542.设定察看的时间(深度)，例如，在Tx后50pec;计算(或测量)这一指定时间Tx条件设定的瞬时Tx脉冲压力场。00553.计算从Tx起点沿着半径方向的脉冲的中心。在2D图像情形中，绘出(精确地或计算地)在2DTx压力场上沿着脉冲中心的曲线；或是在3D图像情形中，绘出在3DTx压力场上沿着脉冲中心的表面。00564.找出脉冲中心曲线和发射波束的交叉点。假设Ro是从Tx起点到交叉点的距离(在时间上)。00575.找出一次波束和脉冲中心曲线的交叉点。假设R!是从一次波束的波束起点到交叉点的距离(在时间上)。对于这一Tx波束在这一深度上的第二阶段BF而言，一次波束延迟等于Ro与Ri之间的差。00586.计算所有一次波束的延迟。这产生在一个深度上的延迟分布。00597.对于沿着整个场的每一深度重复步骤2-6，以便确定来自一个原始发射波束的各一次波束的动态Tx延迟分布。可以通过在计算出的各深度间进行插值来减少计算。00608.对于所有原始发射波束重复步骤1-7，以便确定来自所有原始发射波束的所有一次波束确定Tx延迟分布。0061作为一个简化的示例，参见图7。假设发射波束Txl2和接收波束起源于波束起点7，从而给出PB(12，7)。曲线700表示在发射后Tx脉冲中心为49ps。在图7中，Txl2波束起点被标记为A，而曲线700和Txl2的交叉点被标己为B。假定Ro是从A到B的距离。在该示例中，Ro=49piS。PB(12，7)的波束起点被标记为C，而曲线700和PB(12，7)的交叉点被标记为D。假定R!是从C到D的距离，并且作为举例假设R,二50网。曲线701示出了假设的来自C的脉冲中心50^isec。现在，如果系统将使用起点位于位置7的波束来査看位于D的耙，则它将从C射出发射波束，所述脉冲将传播50psec到达D，然后将被反射回C用于接收。这将对应于一次接收波束PB(7，7)，因而其的往返时间将是5(^s+50^is-100^is。可是，使用两阶段波束形成器，系统能射出Tx12，其的脉冲起源于A，传播49p并到达位于D的靶。然后来自D的回波从D到C用了50^is，以供位置7处波束进行接收。在这种情况下，总的往返时间将是49hs+50^is=99|lis。在波束形成器的第二阶段中，系统因此可以将lp延迟施加给位于该深度的PB(12，7)，以校正Tx波束的偏移量，从而产生看起来像PB(7，7)的波束。0062根据本发明的两阶段波束形成器可以以多种方式进行灵活的编程，以便显示各种特征。通常，如果对阵列和发射控制器进行选择和编程以便能够形成来自一次发射的Np个并行波束，并且如果使用Nt次发射形成图像，并且如果使用k个一次波束形成一个最终的二次波束，那么最终的二次波束的总数将是Nt*Np/k，因为每次发射将给出Np/k个二次波束。0063如果帧率在给定的应用中更重要，则人们例如可以配置系统使得1^2以便每次发射具有Np/2个最终波束。这样，如果Np-16，从一次发射将有八个最终的二次波束，并且由于两个一次波束的第二阶段Tx波束形成器而很少有并行伪影。甚至可以设定k=l以便形成每次发射16个波束(作为举例假设其是选定孔径)，以使帧率最大化，尽管以一些并行伪影为代价。0064如果图像质量最重要，则在Np-16的情况下，人们可以设定k=16以便每次发射具有单个最终的波束。然后从16个一次波束动态地形成每个二次波束；因此，最终的波束将沿着深度具有良好Tx聚焦。0065为了更深的穿透，再次使用Np-16作为示例，人们同样可以设定k-16，从而每次发射产生单个最终的二次波束。然后可以由16个一次波束动态地进行延迟求和得到每个最终波束，以便使得每个最终波束将包括16个原始发射波束的声功率。这应当使2.5MHz的探头增加大致12dB的穿透，这等价于明显增加大约4.8cm的深度。0066在图8中为了清楚起见去除了图1中的一些部件，图8示出了Rx控制器840如何包括第一阶段波束形成器842和第二阶段波束形成器844。同样示出了Tx控制器820的两个公共部件，即孔径控制器822和变迹控制器824。0067并不一定将本发明的两阶段波束形成器实iJl为两个硬件阶段。例如，它实践上可以是将第二阶段的动态Tx延迟合并到第一阶段的动态Rx延迟。然而，第一阶段的动态Rx延迟可以如现有技术一样根据几何结构的函数进行计算；这可以使用硬件计算器来实现以便减少传输更多控制数据的需要。但是，对于第二阶段的《的计算，并不像上面讨论的那样简单。代替的是，可以根据存储在主系统存储器中或者一些其他存储设备中的表860，在实时成像期间对这些值预先计算、预先存储和加载；而且，所述值可以进行实时、或脱机、或者通过这两种方式的组合的计算。0068可以使用存储器/累加器，在第二阶段波束形成器中完成SB(p)值的计算，所述存储器/累加器带有存储器单元(对于速度专用的，或者其他存储器的一部分，例如系统存储器)，它的大小至少要像包含定义所需波束数量的所有值需要的那么大，并且至少等于每次发射的最大并行波束数量。一旦在第二接收波束形成器阶段中计算完SB(p)值，然后将这些值进行后面的常规处理，以便形成感兴趣区域的显示基础。0069本发明的原理同样可以应用到更高维阵列。例如，图9示出了2D阵列900的8x16元件部分，即，具有i=8"行"和j=16"列"的元件。假定E(i=a，b;j=x，y)表示包括a到b行和x到y列的元件的阵列部分。0070与一维情况一样，同样可以定义波束号以便标记2D阵列中各发射波束的可能的不同起点。图9b中示出了用于2D阵列的3乘3区域(如图9a中所示，左上区域)的发射波束起点的一种可能编号方案，其中圆圈指示波束起点而括号中相邻的二维标记指示波束号。这样，波束号(3,3)(和因而相应的起源于所指示点的波束Tx(3，3))位于元件E(2，2)的中心，波束号(4,5)位于E(2，3)和E(3，3)间边缘的中间等等。例如，其他的编号方案可以仅包括位于元件的角和中心的任何起点位置；编号方案仅为在此用于说明本发明原理的规则，尽管一些编号方案通常将需要对软件进行编程以便更容易地控制发射和接收波束形成器。0071如上面提到的，在2D阵列中可以激励部分行或列、或者全部行或列的元件，以便它们将在效果上形成1D阵列。这总体而言对于诊断超声成像，或者具体而言对于实现本发明都不是必须的；相反，本发明的原理可以在比1D更高的维度中实施。作为举例，假设第一发射波束Tx(4，12)由带阴影的元件4x4子组E(i^,4;〗=5，8)形成，并且波束起点对应于显示为空心小圆圈的波束号(4,12)。作为简化假设，用选定的接收孔径进行扫描，以便使得第一一次接收波束(由实心点指示的)对应于具有波束号(4,4)的波束起点，然后在j方向(向右，如图9a和9b所示)将所述孔径移位一个元件(两条波束线)，对于波束号(4，6)计算新的一次波束，然后再次将所述接收孔径向右移位一个元件，对于波束号(4,8)计算另一个一次波束，然后再次将所述接收孔径向右移位一个元件，并且对于波束号(4,10)计算另一个一次波束。0072PB[(ti，tj);(ra,rb)]然后可以用于表示由具有波束号(ti,tj)所指示的起点的发射波束产生的、并在具有(使用同一方案编号的)号(ra，rb)的波束线上所接收的一次接收波束。这样，在先前段落给出的示例中，对于发射波束Tx(4，12)，以任何常规方式计算一次波束PB[(4,12);(4，4)]、PB[(4，12)(4，6)]、PB[(4,12>(4,8)]和PB[(4,12);(4,10)]。0073现在，可以将发射孔径移位例如两个元件(四个波束号)，以便产生来自位于波束号(4，16)的起点的发射波束Tx(4,16)，在图9c中由位于阴影孔径区域元件E(i^，4;〗=7,10)中心的空心圆圈来指示。对于如前所示的接收孔径(再次用实点指示的)使用相同的模式，计算四个一次接收波束为PB[(4，16);(4,8)]、PB[(4,16);(4，10)]、PB[(4,16);(4，12)]和PB[(4,16);(4,14)]。0074在这一简单的示例中，一次接收波束PB[(4，12);(4,8)](来自第一发射)和PB[(4，16M4,8)](来自第二发射)二者都是在相同波束线、即(4，8)上的一次接收波束。同样PB[(《1^(W0)]和PB[(^W)^，10)]都在波束线(4，10)上。然后根据加权(再次，甚至可以为统一加权)时延求和将这些对的每一对合并成单个二次接收波束SB(4，10)。这些"重叠的"接收波束在图9c中被标记为用虚线矩形围绕。0075对于使发射孔径或接收孔径将要单独在方位角方向或仰角方向上进行移位并不是必须的。即使对一个或另一个、或者两者的移位包括"对角"分量，只要对于相同的波束线形成一个以上的一次接收波束，那么这些波束可以合并成二次波束。在连续的发射和接收上进行移位的程度通常将确定"重叠"的程度，因而确定合并多少项(一次接收波束)来形成单个、相应的二次接收波束。0076本发明提供一种多阶段波束形成器(具有关于发射和接收波束形成两者的特征)，其解决了现有技术中的若干缺点它能用亮光检査(candle)来自单次发射的大量并行接收波束。在第二阶段，它能够进行动态发射聚焦、动态发射孔径大小控制和动态发射变迹。根据对诸如孔径、重叠程度等的上述参数的一些选择，这一波束形成器同样允许增加帧率。它还可以提高发射波束形成器的质量和更充分地使用声学可用不功率。权利要求1、一种超声成像方法，其特征在于激励电声元件阵列(310)中的至少一个元件，以便产生一系列声穿透感兴趣区域(400)的发射波束(Tx)；对于至少一个所述发射波束，确定多个一次接收波束(PB)；根据将至少两个所述一次接收波束进行时间对准合并来计算至少一个二次接收波束(SB)；以及根据至少一个所述二次接收波束产生图像。2、如权利要求1所述的方法，进一步的特征在于时间对准所述至少两个一次接收波束，所述至少两个一次接收波束对应于来自不同的所述发射波束的超声返回信号但却在共同的波束线上进行接收；以及根据所述至少两个时间对准的一次接收波束的函数来计算每个二次接收波束。3、如权利要求2所述的方法，进一步的特征在于根据对所述至少两个一次接收波束的时延求和来计算每个二次接收波束。4、如权利要求3所述的方法，进一步的特征在于在所述求和中对每个一次接收波束的贡献进行加权。5、如权利要求2所述的方法，其中每条波束线是穿过焦点和所述阵列上的相应波束起点的超声传播线。6、如权利要求l所述的方法，其中时间对准所述至少两个一次接收波束的步骤包括计算至少一个所述一次接收波束相对于另一个所述一次接收波束的时间延迟。7、如权利要求6所述的方法，其中计算所述时间延迟的步骤包括模拟与第一原始波束的预期发射条件相对应的条件；对于指定时间下的所述条件，确定所述第一原始波束上的发射脉冲的第一瞬时发射脉冲压力场(700)，其对应于指定深度；确定沿着从所述第一原始波束的第一波束起点(A)的至少一个半径方向的所述发射脉冲的中心；确定在所述第一发射起点(A)与所述指定深度处的所述瞬时发射脉冲压力场和所述第一原始波束的交叉点(B)之间的第一发射时间(R0);确定在第二波束起点(C)与对于第二个一次接收波束的第二脉冲压力场的交叉点(D)之间的第二发射时间(Rl);根据预先确定的关于所述第一发射时间和所述第二发射时间的函数确定时间差分；以及在计算至少一个所述二次接收波束的过程中，向所述第二个一次接收波束施加根据预先确定的关于所述时间差分的函数而计算出的延迟。8、一种超声成像系统，其特征在于主处理系统(200)，其包括至少一个处理器(210)和至少一个存储器(250);超声探头(300)，其包括电声元件阵列(310);波束形成系统(820)，其包括发射(Tx)波束形成器和接收(Rx)波束形成器；其中将所述发射波束形成器设置成激励所述阵列的至少一个元件以便产生一系列声穿透感兴趣区域(400)的发射波束；所述接收波束形成器包括第一阶段(844)，对于每个发射波束，将所述第一阶段设置成确定多个一次接收波束；所述接收波束形成器包括第二阶段(842)，将所述第二阶段设置成根据对至少两个所述一次接收波束的时间对准合并来计算至少一个二次接收波束；以及将所述主处理系统设置成根据至少一个所述二次接收波束产生图像，并设置成产生用于在显示装置(100)上显示所述图像的相应信号。9、如权利要求8所述的系统，其中将所述第二阶段设置成时间对准所述多个一次接收波束，所述多个一次接收波束对应于来自不同的所述发射波束的超声返回信号但却在共同的接收波束上进行接收；以及根据所述至少两个时间对准的一次接收波束的函数来计算每个二次接收波束。10、如权利要求9所述的系统，其中将所述第二阶段进一步设置成根据对所述包括的各一次接收波束的时延求和来计算每个二次接收波束。11、如权利要求8所述的系统，其中所述阵列是一维阵列。12、如权利要求8所述的系统，其中所述阵列是二维阵列。全文摘要根据发射波束形成器激励超声阵列中的元件(310)，以便产生一系列发射波束(Tx)。对于每个发射波束，接收波束形成器的第一阶段(844)确定多个一次接收波束(Rx；RB)。第二波束形成器阶段(842)然后根据一次接收波束计算二次接收波束，所述一次接收波束对应于从不同的发射波束到共同接收起点的返回信号。例如，可以根据对一次接收波束的加权、时延求和来计算每个二次接收波束。然后从波束形成器输出至少一个二次接收波束以用于产生显示图像。文档编号A61B8/00GK101427929SQ20081017992公开日2009年5月13日申请日期2008年10月23日优先权日2007年10月23日发明者姚林鑫申请人:迈瑞控股(香港)有限公司