CN102697524B - 全聚焦超声成像方法及其在血流成像中的运用 - Google Patents

全聚焦超声成像方法及其在血流成像中的运用 Download PDF

Info

Publication number
CN102697524B
CN102697524B CN201210135025.3A CN201210135025A CN102697524B CN 102697524 B CN102697524 B CN 102697524B CN 201210135025 A CN201210135025 A CN 201210135025A CN 102697524 B CN102697524 B CN 102697524B
Authority
CN
China
Prior art keywords
wave
ultrasonic
blood flow
imaging
array
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
Application number
CN201210135025.3A
Other languages
English (en)
Other versions
CN102697524A (zh
Inventor
�田睦
熊佑全
Original Assignee
CHENGDU YOUTU TECHNOLOGY CO LTD
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by CHENGDU YOUTU TECHNOLOGY CO LTD filed Critical CHENGDU YOUTU TECHNOLOGY CO LTD
Priority to CN201210135025.3A priority Critical patent/CN102697524B/zh
Publication of CN102697524A publication Critical patent/CN102697524A/zh
Application granted granted Critical
Publication of CN102697524B publication Critical patent/CN102697524B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

全聚焦超声成像方法及其在血流成像中的运用,本发明步骤包括由B超的超声波换能探头以点波源发射方式依次发射超声波信号;并由超声波换能探头上多个阵元接收形成的接收波束,将前后各点阵元发射时所产生的接收波束进行叠加、解调、动态滤波、坐标转换、图像平滑、图像滤波数字信号处理,便形成图像。将前后各点阵元发射时所产生的接收波束进行相位校正、多普勒相位检测便获得血流速度,对所包含的多普勒效应的相位在每一条波束上的每一点做多普勒效应投影校正,形成2D血流图。本方案的发射声功率极低、帧频高、适应血流速度范围广、血流速度精准。

Description

全聚焦超声成像方法及其在血流成像中的运用
技术领域
[0001] 本发明涉及一种医疗用诊断仪器,更进一步说明为B超成像仪器及方法。
背景技术
[0002] 一般的B超工作过程为:当探头获得激励脉冲后发射超声波(同时探头受聚焦延迟电路控制,实现声波的声学聚焦),然后经过一段时间延迟后再由探头接受反射回的回声信号,探头接收回来的回声信号经过滤波,对数放大等信号处理。然后由DSC电路进行数字变换形成数字信号,在CPU控制下进一步进行图像处理,再同图表形成电路和测量电路一起合成视频信号送给显示器形成我们所熟悉的B超图像,也称二维黑白超声图像。
[0003] 但是,传统模拟B超成像采用数控延迟对发射进行波束形成,在接收聚焦上采用延迟线进行多段聚焦。由于每次发射波束只能对一点进行聚焦,离聚焦点远的地方就模糊,所以一般都要发射2~4次再合成,以保证图像远近都较为清晰。这种做法图像每秒桢数少,且只能在少数几点保证聚焦。在接收方面,采用FPGA进行逐点聚焦.与模拟B超相比,在同等条件下,这种方式大为提高了图像质量。以下是一个关于各B超特点的对照表:
Figure CN102697524BD00031
[0005] 在超声成像市场上,几乎99%的产品都采用发射多个阵元(通常在32路以上),形成单一方向的发射波束,并且在发射聚焦点上形成最强发射信号。这种发射方式通常能够在发射聚焦点上形成很好的图像,产生高分辨率。但是在远离发射焦点的地方,其分辨率逐步降低。图像不均匀一致。
[0006] 同样,这种发射方式在彩超中对血流速度的测量也带来帧频过低的问题。因为对每一点的血流速度测量,都需要进行8~15左右的重复发射,才能估算出多普勒效应,这导致在血流速度测量前,必须先选择一个小范围做测量对象,达到减少每帧的发射次数。并且,线密度也因为要提高帧频而降低。举个例子,如果是128阵元,探测深度为18cm,线密度为1,重复发射为15,那么要实现全屏彩色血流的话,帧频会在每秒I帧。通常人们会把线密度降低为0.5,最多实现半屏血流,同时实现4波束合成,这样会让帧频升高至每秒16帧。
[0007] 由于帧频问题,同一方位只能采集到最多发射15次的血流信息,这对估算血流速度带来非常大的困难。可以这么说,在彩色血流成像方面,有相当数量的学术文章都是在研究如何在样本少的情况下尽可能地去估算血流速度。
[0008] 还有一个问题就是对水平于探头的血管血流速度无法测量,因为发射声束是垂直于血流方向,从而多普勒效应很小而无法观察。通常人们做法是人为手工改动发射方式为斜射。这也导致操作上的麻烦。
发明内容
[0009] 为了解决上述问题,本发明提供了一种聚焦超声成像方法及其在血流成像中的运用。
[0010] 本发明的全聚焦超声成像方法如下:
[0011] 该方法包括使用超声主机系统、视频输出设备、超声换能器及超声换能器上的发射阵元、超声主机系统用作信号处理的单片机,并按以下方式进行一帧图像的扫描:
[0012] a、由B超的超声波换能探头从一侧至另一侧以球面波发射方式发射超声波信号,对机体组织进行扫描。
[0013] 球面波发射方式根据其设置,每次以1-4个阵元进行发射。
[0014] b、并由超声波换能探头上多个阵元接收每一个球面波周围形成任意接收波束,任意接收波束为I条以上,接 收波束方向为垂直于探头平面和倾斜超声波换能探头平面;并将接收波束信息返馈至超声主机系统。
[0015] C、超声主机系统将前后每次球面波发射时所产生的接收波束进行叠加,扫描完一侦,即可获得一幅图像,其中包含了多条被叠加过多次的接收波束,并将该接收波束用于常规成像处理,最后在视频输出设备上显于。此处常规成像处理包含解调、动态滤波、坐标转换、图像平滑、图像滤波数字信号处理。
[0016] 如上所述的本发明的全聚焦超声成像方法,更进一步说明为,每次以单个阵元发射;或者每次以1-4个阵元方式发射,发射时,则依次相交并叠加方式,计(阵元1,阵元2),(阵元2,阵元3);(阵元3,阵元4),直至最后一阵元。
[0017] 本发明的全聚焦超声成像在血流成像中的运用方法如下:
[0018] 该方法包括使用一台超声主机系统、视频输出设备、超声换能器及超声换能器上的发射阵元、超声主机系统用作信号处理计算的单片机,并按以下方式对血流进行一帧图像的扫描:
[0019] a、由B超的超声波换能探头从一侧至另一侧以球面波发射方式依次向血流发射超声波信号,超声波信号为一组连续的脉冲波;球面波发射方式根据其设置,每次以1-4个阵元发射。
[0020] b、并由超声波换能探头上多个阵元接收每一个球面波周围形成任意接收波束,任意接收波束是I条以上的波束,波束方向为垂直于探头平面和倾斜探头平面。
[0021] C、超声主机系统跟据以上波束信息,作血流分析计算,并在视频输出设备上显于血流图。
[0022] 血流分析计算方式为:
[0023] a、根据发射球面波的位置做相位校正所形成的波束组。
[0024] b、对一帧图像中的每一线位置,以多条波束线对应于不同阵元发射,对这多条波束线做多普勒相位检测获得血流速度。[0025] C、对所包含的多普勒效应的相位在每一条波束上的每一点做多普勒效应投影校正。即把对应发射球面波而产生的多普勒效应投影在垂直向上Y方向或/和水平向右X方向。同时投影在Y方向和X方向就可以获得两个方向的血流速度,进而形成2D血流图。
[0026] d、对所包含的多普勒效应的幅度在每一条波束上的每一点做多普勒效应幅度校正,其校正方式为乘以一个衰减系数,该衰减系数根据球面波发射位置到处理点的距离与该处理点到所在波束的第一点的距离而计算得来。
[0027] 本发明的有益效果包括如下:
[0028] 声功率超低:对比常规128阵元,64路通道的机器,本方案的发射声功率为1/32或1/64。这对于许多特殊运用的产品方向起着革命性的改进。
[0029] 帧频高:在128阵元,18cm的探测深度下,其黑白图像帧频可以达到15帧,全屏2D血流图同样可以达到15帧每秒。
[0030] 适应血流速度范围广:由于产生的波束众多,血流采样频率快,能够适用的血流速度范围宽,从而省去大量医生工作去调整设备设置。
[0031] 血流速度精准:由于同一点的血流速度是通过不同角度大量采样得来,其准确度比常规单一发射方向的波束要准确。
[0032] 2D血流图:虽然有不少研究在学院中进行,但在实际应用方面,属于首创。
[0033] 本项技术发明抛开传统的发射波束,采用点阵元发射,运用球面波,通过孔径合成计算,实现发射接收 双逐点聚焦图像,多方位血流速度测量而不受血管方位限制,并且血流速度采集样本可以高达128线,相比传统彩超的最多15次样本,这对精确估算血流速度起到了革新性意义。不仅如此,在帧频方面,还能够以每秒16帧的速度在不降低线密度的情况下实现全屏血流。
附图说明
[0034] 图1扫描并与接收波示意图;
[0035] 图2血流扫描与接收波示意图。
[0036] 附图编号说明:1-发射阵元;2_接收波束;3_ —组脉冲波束;4_血液方向具体实施方式
[0037] 参照图1理解如下内容:
[0038] 本技术的主要特点在于采用点波源发射,而不是多个阵元同时发射产生方向很强的单一波束。同时,在发射阵元I周围,形成多个接收波束2,根据使用的硬件规模,如果是128阵元,并且阵元间线密度是2的情况,最多可以形成256条接收波束。通常64路接收波束已经足够。如果每次发射形成64线,而每帧发射128次,那么一共产生128x64线,或标记为B(m,n),其中m = O~127,标记该波束在图像中的空间位置,η = O~63,标记该波束被产生的时间顺序。注意在边界情况下有伪线产生将会被扔掉,为了讨论方便,这里暂且不考虑边界情况。
[0039] 跟普通超声成像一样,每一帧扫描都从探头最左边开始,在第i个阵元发射的时候,产生以i阵元为对称的64条接收波束,标记为:
[0040] [B (1-32,63) B (1-31,62) B (1-30,61).........................B (i,32) B (i+1,31)….B(i+31,0)]
[0041] 在第i+1阵元发射的时候,同样产生以i+Ι阵元为对称的64条接收波束。
[0042] [B(1-31,63)B(1-30,62)B(1-29,61)......B (i,33) B (i + 1,32)____B(i+31, I)
B(i+32,0)]
[0043] 请注意,这两组64波束中有62条波束是在几何位置上同一位置所形成的波束,需完全相叠加。而B(i,0)则应该为输出,B (i+1,63)则为存储64条波束中的最新波束,将会被叠加63次。扫描完一帧,则形成128条被叠加合成过63次的接收波束线,这128条接收波束将按照普通超声成像处理来做坐标转换等图像处理工作。
[0044] 最终用来显示B模式黑白图像的波束为:BF(0~127,0~4095) =sum(B(0~127,O~63,O~4095),2);血流成像过程实施例如下:
[0045] 血流分析跟B模式成像过程有不同之处。首先,发射脉冲不是单一脉冲,而是一串脉冲3,比如说,20个连续脉冲。其次,在同一位置上形成过的64条接收波束线,并不用来叠加,而是用来做血流分析。参照附图2 (点声源遇到血流4的情况)。
[0046] 由于是球型波,血流4带来的多普勒效应不再是对某一个单一方向,而是全部都相对于声源方向。也就是说,在不同的接收波束上的每一个接收点,该处血管中血流产生的多普勒效应都各不相同,每点的多普勒效应都体现出该处血流方向对应于声源方向所产生的多普勒效应。与传统血流检测相比,这种现象是最大的不同之处,同时也是最有革新性意义之处。
[0047] 检测血流的方法目前最为流行的是相位检测,该方法对同一检测点上发射8~12次相同方向及强度的脉冲,由此产生8~12次接收波束,对这一组波束进行一系列处理计算,即能获得该波束上每一点对应于发射脉冲方向的血流速度(如果有血流的话)。由于每次发射方式相同,在这些波束之间做相关计算分析获得相位的方法变得有意义。但这对球型波来说,正是极具挑战的地方。在每条波束上的每一点上,获得的多普勒信息不仅方向不同,而且能量不同。所以,在同一位置上获得的这64条接收波束都需要在方向和能量上获得统一校正,使多普勒效应投影在接收波束方向,并且补偿由于不同发射阵元而带来的距离衰减,才能进行下一步的相位检测,从而获得在接收波束方向上的血流流速。同样原理,垂直于接收波束方向上的血流流速也可以获得。这样,快速并且准确地显示ID或2D血流就能够实现。
[0048] 多普勒效应投影计算
[0049] 通过128次发射,每次波束形成64线,每条线4098个点,那么我们将有128x64x4096的矩阵。我们标记128阵元方向为X方向,4098个采样点为Y方向,那每个X位置上形成64次波束则为Z方向。
[0050]在Z方向的64线中,只有第32线是从正上方发射而得来的接收波束,其他线都是从左右相邻的阵元发射而得到的,因此,这些线需要经过投影校正,把多普勒效应投影到从正上方发射而形成的多普勒效应。
[0051] 现在考虑在X方向的第48线位置的Z方向的第一条线,该线应该是在第16阵元发射的时候对第48线位置形成的接收波束。这条线被标记为B(48,l,0~4095)。现在考虑该线上的一点,比如说第1000点,根据采样频率和超声速度,该点的物理位置应该在正对第48阵元的下方77mm。现在假设探头为128阵元线阵,阵元间距为0.5mm,如果把第16阵元标记为P16,那么P16到B(48,l,1000)再到P48的夹角为11.73度。由于多普勒效应可以从该点的IQ分解后的IQ之间相位变化上估算得到,那么,IQ之间的相位变化也可以投影在不同方位,这等价于多普勒效应的投影。于是,投影过程可以表示为:
[0052] BC(48,1,1000)_1 = B(48,1,1000)_I ;
[0053] BC(48,1,1000)_Q = B (48,1,1000)_Q*cos (11.73);
[0054] C(48,1,1000) = cos (11.73);
[0055] 从而得到如下投影结果:
[0056] BC = B.*C
[0057] 其中C为投影矩阵,BC为投影后的波束矩阵。需要注意的是,每一波束上的每一点的投影校正都与该点的几何位置和与之相对应的发射阵元的几何位置相关。
[0058] 如果多普勒效应投影在Y方向,则计算形成的血流图像与传统图像一致。如果同时投影在X方向并且计算X方向的血流图像,那么精确的血流速度与方向可以完全得到,因而可以形成2D血流图像。
[0059] 具体方式如下:
[0060] 能量校正:
[0061] 在估算相位变化的时候,如果Z方向的每次采样的能量都不能恒定,则会影响相位估算的准确度。与常规彩超的发射波束不同的是,在Z方向的每条接收波束,其相对应的发射超声波在到达某点时所传播的距离不同,从而导致能量不一致。距离可以容易估算得至|J,还是以B(48,l,1000)为例:其需要修正的距离差为:
[0062] D (48,1,1000) = sqrt ((32*0.5mm) ~2+77mnT2) _77mm ;
[0063] E(48,I,1000) = 10~(D(48,1,1000)*2*Fo/10); [0064] 其中E为实际修正比率值,Fo为超声波中心频率。修正后的波束矩阵为:
[0065] BE = BC.*E
[0066] 血流速度估算:采用常规的自相关相位估算方法。该方法属于成熟多年的算法,在此不做重复描述。

Claims (6)

1.全聚焦超声成像方法,包括使用超声主机系统、视频输出设备、超声换能器及超声换能器上的发射阵元、超声主机系统用作信号处理的单片机进行超声成像,其特征在于: 按以下方式进行一帧图像的扫描: a、由超声换能器从一侧至另一侧以球面波发射方式发射超声波信号,对机体组织进行扫描;球面波发射方式根据其设置,每次以I~4个阵元进行发射; b、并由超声换能器上多个阵元接收每一个球面波周围形成任意接收波束,接收波束方向为垂直于超声换能器平面和倾斜超声换能器平面;并将接收波束信息返馈至超声主机系统; C、超声主机系统将前后每次球面波发射时所产生的接收波束进行叠加,扫描完一中贞,即可获得一幅图像,其 中包含了多条被叠加过多次的接收波束,并将该多条被叠加过多次的接收波束用于常规成像处理,最后在视频输出设备上显示。
2.如权利要求1所述的全聚焦超声成像方法,其特征是球面波发射方式采用依次相交并叠加方式进行发射,直至最后一阵元。
3.如权利要求1所述的全聚焦超声成像方法,其特征是任意接收波束为I条以上。
4.全聚焦超声成像在血流成像中的运用,包括使用一台超声主机系统、视频输出设备、超声换能器及超声换能器上的发射阵元、超声主机系统用作信号处理计算的单片机进行超声血流成像,其特征在于, 按以下方式对血流进行一帧图像的扫描: a、由超声换能器从一侧至另一侧以球面波发射方式依次向血流发射超声波信号,超声波信号为一组连续的脉冲波;球面波发射方式根据其设置,每次以1-4个阵元发射; b、并由超声换能器上多个阵元接收每一个球面波周围形成任意接收波束,任意接收波束是I条以上的波束,波束方向为垂直于超声换能器平面和倾斜超声换能器平面; C、超声主机系统跟据以上波束信息,作血流分析计算,并在视频输出设备上显示血流图。
5.如权利要求4所述的全聚焦超声成像在血流成像中的运用,其特征在于,血流分析计算方式为: a、根据发射球面波的位置,对所形成的波束组进行相位校正; b、对一帧图像中的每一线位置,以多条波束线对应于不同阵元发射,对这多条波束线做多普勒相位检测获得血流速度; C、对所包含的多普勒效应的相位在每一条波束上的每一点做多普勒效应投影校正,即把对应发射球面波而产生的多普勒效应投影在垂直向上Y方向和水平向右X方向;同时投影在Y方向和X方向就可以获得两个方向的血流速度,进而形成2D血流图; d、对所包含的多普勒效应的幅度在每一条波束上的每一点做多普勒效应幅度校正,其校正方式为乘以一个衰减系数,该衰减系数根据球面波发射位置到处理点的距离与该处理点到所在波束的第一点的距离而计算得来。
6.根据权利要求1所述的全聚焦超声成像方法,其特征在于,所述多条被叠加过多次的接收波束用于常规成像处理,此处成像处理包含解调、动态滤波、坐标转换、图像平滑、图像滤波数字信号处理。
CN201210135025.3A 2012-05-04 2012-05-04 全聚焦超声成像方法及其在血流成像中的运用 Expired - Fee Related CN102697524B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210135025.3A CN102697524B (zh) 2012-05-04 2012-05-04 全聚焦超声成像方法及其在血流成像中的运用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210135025.3A CN102697524B (zh) 2012-05-04 2012-05-04 全聚焦超声成像方法及其在血流成像中的运用

Publications (2)

Publication Number Publication Date
CN102697524A CN102697524A (zh) 2012-10-03
CN102697524B true CN102697524B (zh) 2014-08-13

Family

ID=46890794

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210135025.3A Expired - Fee Related CN102697524B (zh) 2012-05-04 2012-05-04 全聚焦超声成像方法及其在血流成像中的运用

Country Status (1)

Country Link
CN (1) CN102697524B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10548571B1 (en) * 2014-11-21 2020-02-04 Ultrasee Corp Fast 2D blood flow velocity imaging
CN105011965A (zh) * 2015-06-10 2015-11-04 熊佑全 合成孔径全息超声扫描的聚焦实现方法
CN109414245A (zh) * 2016-09-30 2019-03-01 深圳迈瑞生物医疗电子股份有限公司 超声血流运动谱的显示方法及其超声成像系统
CN108186050B (zh) * 2018-01-03 2021-02-23 声泰特(成都)科技有限公司 一种基于超声通道数据的多普勒血流速度成像方法和系统
CN110811688B (zh) * 2019-12-02 2021-10-01 云南大学 多角度平面波重复复合的超快超声多普勒血流估计方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101156788A (zh) * 2007-11-14 2008-04-09 成都西云科技有限公司 B超成像方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5615675A (en) * 1996-04-19 1997-04-01 Regents Of The University Of Michigan Method and system for 3-D acoustic microscopy using short pulse excitation and 3-D acoustic microscope for use therein
US10914826B2 (en) * 2008-06-26 2021-02-09 Verasonics, Inc. High frame rate quantitative doppler flow imaging using unfocused transmit beams

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101156788A (zh) * 2007-11-14 2008-04-09 成都西云科技有限公司 B超成像方法

Also Published As

Publication number Publication date
CN102697524A (zh) 2012-10-03

Similar Documents

Publication Publication Date Title
CN102697524B (zh) 全聚焦超声成像方法及其在血流成像中的运用
US5172343A (en) Aberration correction using beam data from a phased array ultrasonic scanner
CN102047140B (zh) 具有引导efov扫描的扩展视野超声成像
US8469887B2 (en) Method and apparatus for flow parameter imaging
JP4795675B2 (ja) 医療用超音波システム
CN1264025C (zh) 干涉成象的方法和装置
US7828731B2 (en) Ultrasonographic apparatus, ultrasonographic data processing method, and ultrasonographic data processing program
CN101176675B (zh) 超声波诊断装置及超声波图像处理方法
CN104080407A (zh) 任意路径的m模式超声成像
JPH0661333B2 (ja) 流動する反射体の大きさ及び方向の表示を求める方法及びコヒーレント作像システム
CN101683274B (zh) 偏转凸阵数字扫描变换和凸阵空间复合的方法及成像系统
US9474510B2 (en) Ultrasound and system for forming an ultrasound image
CN101449983A (zh) 超声线阵偏转成像的扫描变换插值方法和装置
CN104688271B (zh) 合成聚焦超声成像方法和装置
US20080177182A1 (en) Ultrasonic imaging apparatus and method for acquiring ultrasonic image
CN104777485B (zh) 超声二维面阵的三维宽波束小区域快速空化成像方法
CN106680825B (zh) 一种声学阵列成像系统与方法
JP4381028B2 (ja) 超音波診断装置
CN105074498A (zh) 具有对梯形的扇区的空间复合的超声诊断成像系统
US10548571B1 (en) Fast 2D blood flow velocity imaging
JP5714221B2 (ja) 超音波診断装置及び超音波送受信方法
JPH0790026B2 (ja) 超音波診断装置
WO2019127621A1 (zh) 超声成像方法、系统和设备
US6048313A (en) Method and apparatus for fractal-based enhancement of ultrasound imaging
JPH10165400A (ja) 超音波診断装置

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model
ASS Succession or assignment of patent right

Owner name: TIAN YIPING

Free format text: FORMER OWNER: CHENGDU YOUTU TECHNOLOGY CO., LTD.

Effective date: 20150624

TR01 Transfer of patent right

Effective date of registration: 20150624

Address after: 400030, Shapingba District, Chongqing New Dragon Street, No. 106, 2-2

Patentee after: Tian Yiping

Address before: 5 No. 610041 Sichuan province Chengdu high tech Zone Gaopeng Avenue

Patentee before: Chengdu Youtu Technology Co.,Ltd.

C41 Transfer of patent application or patent right or utility model
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: 20140813

Termination date: 20160504