CN105232085B - 基于动态孔径控制的超声剪切波弹性成像方法 - Google Patents
基于动态孔径控制的超声剪切波弹性成像方法 Download PDFInfo
- Publication number
- CN105232085B CN105232085B CN201510794661.0A CN201510794661A CN105232085B CN 105232085 B CN105232085 B CN 105232085B CN 201510794661 A CN201510794661 A CN 201510794661A CN 105232085 B CN105232085 B CN 105232085B
- Authority
- CN
- China
- Prior art keywords
- shear wave
- elements
- displacement
- array elements
- acoustic radiation
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 18
- 230000005855 radiation Effects 0.000 claims abstract description 48
- 238000006073 displacement reaction Methods 0.000 claims abstract description 45
- 239000002245 particle Substances 0.000 claims abstract description 10
- 238000000034 method Methods 0.000 claims description 46
- 239000013598 vector Substances 0.000 claims description 17
- 238000002604 ultrasonography Methods 0.000 claims description 16
- 230000001965 increasing effect Effects 0.000 claims description 9
- 238000002099 shear wave elastography Methods 0.000 claims description 7
- 238000004613 tight binding model Methods 0.000 claims description 7
- 230000003247 decreasing effect Effects 0.000 claims description 6
- 230000001052 transient effect Effects 0.000 claims description 5
- 230000001419 dependent effect Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 10
- 239000003550 marker Substances 0.000 abstract description 8
- 238000010008 shearing Methods 0.000 abstract 3
- 239000011148 porous material Substances 0.000 abstract 1
- 210000001519 tissue Anatomy 0.000 description 28
- 230000003068 static effect Effects 0.000 description 13
- 238000001514 detection method Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000002091 elastography Methods 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- AYEKOFBPNLCAJY-UHFFFAOYSA-O thiamine pyrophosphate Chemical compound CC1=C(CCOP(O)(=O)OP(O)(O)=O)SC=[N+]1CC1=CN=C(C)N=C1N AYEKOFBPNLCAJY-UHFFFAOYSA-O 0.000 description 3
- 238000002126 acoustic radiation force impulse imaging Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002113 ultrasound elastography Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明提供了一种基于动态孔径控制的超声剪切波弹性成像方法,对不同的聚焦深度,开启对应的阵元数目,设置孔径大小,能有效消除栅瓣现象和位移场中的伪位移,得到的质点位移‑时间曲线更符合单点激发剪切波随时间衰减的规律,提升了剪切波传播速度测量准确度。本发明解决了现有技术中在浅聚焦深度时,由于声辐射力栅瓣的影响,导致标记点位移—时间曲线发生畸变,降低剪切波传播速度测量准确度的技术问题。
Description
技术领域
本发明涉及超声弹性成像技术领域,具体涉及一种基于动态孔径控制的超声剪切波弹性成像方法。
背景技术
超声弹性成像是一种以检测生物组织弹性信息为目的的新兴技术,其成像原理是,利用聚焦超声产生的声辐射力(acoustic radiation force,ARF)激励生物组织,根据声辐射力施加前后该组织的形变差求解弹性系数,进而弹性成像。根据声辐射力的施加方式,超声弹性成像可分为三类:准静态法,例如诊断超声中的声流法(Acoustic Streamingin Diagnostic Ultrasound,ASIDU);瞬态法,例如声辐射力脉冲成像(AcousticRadiation Force Impulse Imaging,ARFII)、剪切波弹性成像(Shear Wave ElasticityImaging,SWEI)、超声剪切成像(Supersonic Shear Imaging,SSI);谐波法,例如振动声成像(Vibro-Acoustography)、简谐运动成像(Harmonic Motion Imaging,HMI)。
其中,剪切波弹性成像通过检测由声辐射力激励生物组织产生的剪切波的传播速度,反演出剪切模量,再据此弹性成像。由于只需要聚焦超声单点激励,所以成像过程受外部及生物组织内部环境影响较小,且安全可靠,收到广大临床医护和科研人员的高度重视,近些年来,无论在理论研究和临床应用都取得了众多突破,成为当前超声弹性成像技术中的一个研究热点。超声剪切波弹性成像的核心步骤是测出聚焦声辐射力诱发的剪切波的传播速度,而固定孔径下,当聚焦深度较浅时,聚焦超声产生的声辐射力场会在焦点两侧产生栅瓣,引起位移场在焦点两侧产生伪位移场,从而导致标记质点位移-时间曲线出现畸变,且聚焦深度越浅,畸变程度越大,这样大大降低了剪切波传播速度测量精度。
发明内容
本申请通过提供一种基于动态孔径控制的超声剪切波弹性成像方法,对不同的测量深度,动态控制孔径大小,以获得良好的位移-时间曲线,以解决现有技术中浅聚焦深度时剪切波传播速度不准确的技术问题。
为解决上述技术问题,本申请采用以下技术方案予以实现:
一种基于动态孔径控制的超声剪切波弹性成像方法,包括以下步骤:
S1:在不同的聚集深度z下,通过调整超声换能器的孔径宽度D来控制声辐射力场的形状,以消除声辐射力栅瓣,即:聚集深度z变大,超声换能器的孔径宽度D调大,聚集深度z变小,超声换能器的孔径宽度D调小;
S2:根据超声换能器在孔径宽度D时的有效声压Pe计算声辐射力F;
S3:由声辐射力矢量与三维波动方程的格林函数卷积求得该声辐射力引起质点振动的位移场;
S4:针对每个聚焦深度,利用峰值时间法求解剪切力波的传播速度CT;
S5:按照计算剪切模量μ,其中ρ为被测生物组织密度,CT为步骤S4所得的剪切力波的传播速度;
S6:根据杨氏模量E、剪切模量μ与泊松比ν三者关系式: 重构生物组织的弹性模量并实现不同深度的弹性成像,其中λ为Lame常数。
进一步地,孔径宽度D随焦距f的改变而变化,且满足式中,λ0为载波波长;超声换能器的孔径宽度D=N_elements·d=(N_elements-1)(width+kerf),式中,N_elements为活跃物理阵元数,d为相邻两阵元的中心间距,width为阵元宽度,kerf为阵元间隙;通过调整活跃物理阵元数N_elements的数目实现超声换能器的孔径宽度D的调节,对不同焦距f下的活跃物理阵元数N_elements按照进行动态调整,N_elements取整数,且当阵元总数为偶数时,活跃物理阵元数N_elements取偶数,当阵元总数为奇数时,活跃物理阵元数N_elements取奇数。
动态孔径控制即针对不同的聚焦深度,开启相应的活跃物理阵元数目,当聚焦深度较浅时,只开启少数的中心阵元,边缘的其他阵元则处于静默状态,此时的孔径宽度D较小,随着聚焦深度的逐渐增加,相应的边缘阵元也随之开启,增大了孔径宽度D,从而控制声辐射力场的形状,有效消除栅瓣现象和位移场中的伪位移,使得得到的聚焦中心侧向标记点的位移-时间曲线更符合对剪切波传播速度测量的要求。
更进一步地,步骤S2中声辐射力式中,α为生物组织的声衰减系数,c为超声在生物组织中的传播速度,I为一空间位置时间平均声强,且式中,ρ为生物组织的组织密度,Z为生物组织的声阻抗,Pe为当前时刻的有效声压。
实施时,步骤S3中按照求得该声辐射力引起质点振动的位移场,式中,为瞬态作用力矢量在粘弹性介质中引起的位移场,和为方向向量,t和τ为时间,为格林函数分量。
更具体的,步骤S4中的峰值时间法求解剪切力波的传播速度CT的具体步骤如下:
S41:在聚焦位置的横向方向上,取间距相同的多个标记点,记录每个采样时刻对应的振动位移,得到每个标记点的振动位移时间图;
S42:记录每个标记点振动位移主峰值对应的时刻;
S43:利用最小二乘法拟合出以各标记点的位置信息为自变量,对应的主峰值到达的时刻为因变量的一条直线,该条直线的斜率的倒数即为剪切波传播速度值。
使用峰值时间法(Time To Peak displacement,TTP)最大特点在于求解过程中关注的是各个点主峰值到达时刻,与峰值的具体值无关,因此,利用该方法求解剪切波传播速度误差较小。
与现有技术相比,本申请提供的技术方案,具有的技术效果或优点是:通过动态改变超声换能器孔径的大小来控制声辐射力场的形状,有效消除栅瓣现象和位移场中的伪位移,大大提高了浅聚焦深度的剪切波传播速度的测量准确率。
附图说明
图1为动态孔径控制示意图;
图2为基于动态孔径控制的聚焦超声发射示意图;
图3为聚焦超声诱发生物组织产生剪切波示意图;
图4为某一标记点离焦点不同侧向位置的位移-时间曲线图;
图5为焦距5mm静态孔径控制下标志点归一化位移-时间曲线图;
图6为焦距5mm动态孔径控制下标志点归一化位移-时间曲线图;
图7为焦距7mm静态孔径控制下标志点归一化位移-时间曲线图;
图8为焦距7mm动态孔径控制下标志点归一化位移-时间曲线图;
图9为焦距9mm静态孔径控制下标志点归一化位移-时间曲线图;
图10为焦距9mm动态孔径控制下标志点归一化位移-时间曲线图。
具体实施方式
本申请通过提供一种基于动态孔径控制的超声剪切波弹性成像方法,对不同的测量深度,动态控制孔径大小,以获得良好的位移-时间曲线,以解决现有技术中浅聚焦深度时剪切波传播速度不准确的技术问题。
为了更好的理解上述技术方案,下面将结合说明书附图以及具体的实施方式,对上述技术方案进行详细的说明。
实施例
一种基于动态孔径控制的超声剪切波弹性成像方法,包括以下步骤:
S1:在不同的聚集深度z下,通过调整超声换能器的孔径宽度D来控制声辐射力场的形状,以消除声辐射力栅瓣,即:聚集深度z变大,超声换能器的孔径宽度D调大,聚集深度z变小,超声换能器的孔径宽度D调小;
动态孔径控制是接收回波信号环节中的一个常用技术,而在常规超声成像系统中,发射与接收电子聚焦的物理学和数学原理是相同的。因此,本发明利用类似常规超声成像系统接收模式下的动态孔径控制来实现不同焦距下合理有效激励剪切波产生过程中的动态聚焦。动态孔径控制过程如图1所示。发射聚焦时,通常采用延迟顺序激励阵元的方法,使各阵元按照设计的延时依次发射超声波,在焦点处同相叠加干涉增强,焦点外异相叠加干涉减弱。要求各阵元发射超声时间的关系为:
式(1)中,t0为中心阵元的发射时间,tfn为相对于中心阵元第n阵元发射时间,f为焦距,c为超声波在被测组织中的传播速度,D为孔径大小,从上式可知,最大延迟时间tfnmax,即中心阵元与边缘阵元发射时间的差值t0-tfn随焦距减小而增大。
在聚焦孔径中,对任意两阵元element A和element B,假设element A相对于element B的延迟时间为tfAB,对于空间某一点,element A相对于element B的声程差为Δl,若满足式(2)关系,该点则会出现干涉加强。
结合式(1)、(2)可知,在被测组织浅部检测剪切波传播速度时,大孔径聚焦则会引起过大的tfnmax,从而带来大概率的焦点外的二次干涉加强,造成声辐射力场出现栅瓣,引起标记点位移-时间曲线发生畸变。因此,在浅聚焦位置,需要控制孔径大小,才能保证较高的剪切波传播速度检测准确度。限定最大延迟时间小于一个载波周期,则能促使所有阵元发射的超声波仅在焦点处干涉加强,焦点以外区域干涉减弱。这需要孔径大小D必须随焦距f的改变而变化,且要求:
式(3)中,λ0为载波波长。孔径大小D由阵元宽度width与阵元间隙kerf决定,表示为:
D=N_elements·d=(N_elements-1)(width+kerf) (4)
式(4)中,N_elements为阵元数,d为相邻两阵元的中心间距,width为阵元宽度,kerf为阵元间隙。在设计孔径时,可通过参数N_elements、width和kerf来设置孔径的大小,但在实际应用中,对于某一运行中的孔径,阵元宽度width和阵元间隙kerf往往不易实时改变。因此,针对不同焦距f下的聚焦,采用动态孔径控制的方法,通过开/关电路控制阵元的数量来设置适宜的孔径大小。同时为了保证对剪切波的识别度,即足够大的剪切波振幅,本发明对不同焦距f下阵元数量按照以下关系式进行动态调控:
基于动态孔径控制的聚焦超声发射过程如图2所示。相比于静态孔径的聚焦超声发射模式,动态孔径控制添加了一个调节阵元数量的物理阵元开/关控制模块,以针对不同的焦距,开启对应的适宜的阵元数目。当焦距较小时,只开启少数位于中心的阵元,其它阵元则被关闭。开启的阵元处于活跃状态,正常发射超声波,未被开启的阵元处于静默状态,不执行任何命令。随着焦距的增大,开启的阵元数量逐渐增多,直至所有阵元都被开启。考虑到电子聚焦超声的对称性,若阵元总数为偶数,对不同焦距下的阵元数量取满足关系式(5)的偶数个,若阵元总数为奇数,则对不同焦距下的阵元数量取满足关系式(5)的奇数个。
S2:根据超声换能器在孔径宽度D时的有效声压Pe计算声辐射力F,单位(kg(s2cm2));
其中,声辐射力式中,α为生物组织的声衰减系数(Np/m),c为超声在生物组织中的传播速度,通常设置为1540m/s,I为一空间位置时间平均声强(W/cm2),且式中,ρ为生物组织的组织密度(kg/m3),Z为生物组织的声阻抗(N·s/m3),Pe为当前时刻的有效声压(Pa),可以由计算,其中,TP为脉冲周期,P(t)为脉冲实时声压。
S3:由声辐射力矢量与三维波动方程的格林函数卷积求得该声辐射力引起质点振动的位移场;
S4:针对每个聚焦深度,利用峰值时间法求解剪切力波的传播速度CT;
S5:按照计算剪切模量μ,其中ρ为被测生物组织密度,CT为步骤S4所得的剪切力波的传播速度;
S6:根据杨氏模量E、剪切模量μ与泊松比ν三者关系式: 重构生物组织的弹性模量并实现不同深度的弹性成像,其中λ为Lame常数。
实施时,步骤S3中按照求得该声辐射力引起质点振动的位移场,式中,为瞬态作用力矢量在粘弹性介质中引起的位移场,和为方向向量,t和τ为时间,为格林函数分量。
步骤S4中的峰值时间法求解剪切力波的传播速度CT的具体步骤如下:
S41:在聚焦位置的横向方向上,取间距相同的多个标记点,记录每个采样时刻对应的振动位移,得到每个标记点的振动位移时间图;
S42:记录每个标记点振动位移主峰值对应的时刻;
S43:利用最小二乘法拟合出以各标记点的位置信息为自变量,对应的主峰值到达的时刻为因变量的一条直线,该条直线的斜率的倒数即为剪切波传播速度值。
超声剪切波弹性成像(Shear Wave Elasticity Imaging,SWEI)是指利用超声传感器产生聚焦声辐射力诱发生物软组织产生剪切波,通过检测剪切波的传播状态,重构该组织的弹性模量。超声传感器发射聚焦超声诱发生物组织产生剪切波的原理示意图如图3所示,并定义声辐射力在z轴上聚焦。
剪切波传播速度与剪切模量的关系式为假设生物组织的密度ρ是恒定的,则CT与成正比关系。因此,对生物组织剪切模量的计算,其关键是准确测量出剪切波传播速度。
欲求出剪切波传播速度,需首先求解出施加到生物组织上的声辐射力场。在给定的空间位置,声辐射力F可以表示为
求出声辐射力场后,可以用辐射力矢量与三维波动方程的格林函数卷积求得该辐射力引起质点振动的位移场。
定义为i方向的作用力在j方向的格林函数分量,可以表示为体积分量剪切分量以及它们的耦合分量的和,可以表示为:
式中,为方向向量,t为时间,ai为源在i方向上的分量;
其中 和分别表示如下:
式中,为方向向量,r为方向向量的模量,γi为指向的单位向量i方向的分量,γj为指向的单位向量j方向的分量,ρ为介质的密度,υp为运动体积粘度,υs为运动剪切粘度,cp为沿着超声波传播方向(纵波)的波速,cs为沿着垂直于超声波传播方向(横波/剪切波)的波速,δij为kronecker符号
空间坐标系的三个基本轴表示为i=1,2,3,基于这个坐标系的三个分量为(x1,x2,x3),模量指向的单位向量表示为其中γi=xi/r。
因此,耦合分量更为清楚的表达式:
其中,是误差函数。这样,粘弹性格林函数的耦合项可以简化为:
其中,
利用叠加原理,对于任意瞬态作用力矢量在粘弹性介质中引起的位移场可以通过与以上推导的三维波动方程的格林函数进行卷积求得,表示为:
根据上述方法,本实施例利用超声仿真平台FieldⅡ仿真出聚焦声辐射力激励下质点的振动位移场。图4显示了某一标记点离焦点不同侧向位置的位移-时间曲线图(焦距20mm),从图4中可以看出,在声辐射力激励后,质点的振动位移先增大,然后逐渐衰减,整个侧向位置上的质点振动所形成的剪切波逐渐向两端传播。
对于每个聚焦位置,利用峰值时间法结合最小二乘法求解剪切波的传播速度。该方法首选根据焦点侧向标记点的位移-时间曲线图记录下各个标记点振动位移主峰值到达的时刻,再用最小二乘法拟合出以各标记点的位置信息为自变量,相对应的主峰值到达的时刻为因变量的一条直线,该直线的斜率的倒数便是剪切波传播速度值。可见,标记点位移-时间曲线图的走势直接影响剪切波传播速度检测结果,走势良好的标记点位移-时间曲线图符合单点激发剪切波随时间衰减的规律,呈现出明显主峰。利用峰值时间法结合最小二乘法计算剪切波传播速度的最大优势在于,该方法关注的是峰值达到时刻,而与峰值具体值无关,对微位移的识别具有鲁棒性,计算误差低。
在同一聚焦深度,可以通过改变参数f_number来控制声辐射力场的形状,f_number为声辐射力聚焦深度与孔径宽度的比值,用代数式表示为:
式中,z为声辐射力聚焦深度(mm),D为活跃孔径宽度(mm)。
为了进一步验证本发明的有益效果,本实施例通过超声仿真平台FieldⅡ对多个小焦距下分别基于静态孔径和动态孔径聚焦超声发射方式的声辐射力场及标记点位移-时间曲线图进行仿真测试,利用峰值时间法结合最小二乘法计算剪切波传播速度,选取具有代表性的3个小焦距5mm、7mm和9mm下的仿真结果进行展示。仿真中,剪切波在介质中的传播速度理论值设置为2m/s,对孔径参数的设置如表1。
表1孔径参数设置
实验结果表明,在3个焦距下,基于静态孔径的声辐射力场均出现有栅瓣,且焦距为5mm时的栅瓣幅值明显高于7mm与9mm;基于动态孔径控制的聚焦超声发射方式则有效改善了小焦距下的声辐射力,在对应的每个焦距下的声辐射力场均未出现栅瓣。
图5、6、7、8、9、10分别为5mm、7mm和9mm焦距下,基于静态孔径和动态孔径在焦点侧向2mm、3mm、4mm和5mm这4个位置的标记点的位移-时间曲线图。从图5、6中可以看出,当焦距为5mm时,基于静态孔径控制的位移-时间曲线图出现畸变,无明显主峰,基于动态孔径控制的位移-时间曲线图则走势良好,符合单点激发剪切波随时间衰减的规律,呈现出明显的主峰。当焦距为7mm和9mm时,基于静态孔径控制的位移-时间曲线图均出现不同程度的畸变,相比之下,基于动态孔径控制的位移-时间曲线图则走势良好。
考虑到临床上受功率限制,诊断超声诱发组织的振动位移通常较小,检测系统对剪切波的识别度不高以及最小二乘法对标记点个数的要求等因素,此次仿真取焦点侧向0mm、1mm、2mm……9mm这10个标记点的位置信息与主峰值到达时刻来计算剪切波传播速度。表2为5mm、7mm和9mm焦距下,分别基于静态孔径和动态孔径聚焦超声发射方式的剪切波传播速度的检测结果。由表2可见,在三个小焦距下,基于动态孔径的剪切波传播速度的检测结果较传统基于静态孔径的检测结果更接近理论值,与理论值的相对误差更小。例如焦距为7mm时,与理论值相对误差下降了16.585%;焦距为9mm时,相对误差下降了15.205%。另外,在焦距为5mm时,基于静态孔径的标记点位移-时间曲线图由于主峰值到达时刻出现超前现象(距焦点远的标记点的主峰值到达时刻先于距焦点近的标记点的主峰值到达时刻),因此没有剪切波传播速度检测结果。综上可见,基于动态孔径的剪切波传播速度检测准确度高于传统基于静态孔径的检测准确度。
表2不同聚焦深度下剪切波传播速度检测准确度对比
本申请的上述实施例中,通过提供一种基于动态孔径控制的超声剪切波弹性成像方法,对不同的聚焦深度,开启对应的阵元数目,设置适宜的孔径大小,通过超声仿真平台FieldⅡ在浅聚焦位置,减少孔径大小能有效消除栅瓣现象和位移场中的伪位移,得到的质点位移—时间曲线更符合单点激发剪切波随时间衰减的规律,提升了剪切波传播速度测量准确度。本发明解决了现有技术中在浅聚焦深度时,由于声辐射力栅瓣的影响,导致标记点位移—时间曲线发生畸变,降低剪切波传播速度测量准确度的技术问题。
应当指出的是,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的普通技术人员在本发明的实质范围内所做出的变化、改性、添加或替换,也应属于本发明的保护范围。
Claims (5)
1.一种基于动态孔径控制的超声剪切波弹性成像方法,其特征在于,包括以下步骤:
S1:在不同的聚集深度z下,通过调整超声换能器的孔径宽度D来控制声辐射力场的形状,以消除声辐射力栅瓣,即:聚集深度z变大,超声换能器的孔径宽度D调大,聚集深度z变小,超声换能器的孔径宽度D调小;
S2:根据超声换能器在孔径宽度D时的有效声压Pe计算声辐射力F;
S3:由声辐射力矢量与三维波动方程的格林函数卷积求得该声辐射力引起质点振动的位移场;
S4:针对每个聚焦深度,利用峰值时间法求解剪切力波的传播速度CT;
S5:按照计算剪切模量μ,其中ρ为被测生物组织密度,CT为步骤S4所得的剪切力波的传播速度;
S6:根据杨氏模量E、剪切模量μ与泊松比ν三者关系式: 重构生物组织的弹性模量并实现不同深度的弹性成像,其中λ为Lame常数。
2.根据权利要求1所述的基于动态孔径控制的超声剪切波弹性成像方法,其特征在于,孔径宽度D随焦距f的改变而变化,且满足式中,λ0为载波波长;超声换能器的孔径宽度D=N_elements·d=(N_elements-1)(width+kerf),式中,N_elements为活跃物理阵元数,d为相邻两阵元的中心间距,width为阵元宽度,kerf为阵元间隙;通过调整活跃物理阵元数N_elements的数目实现超声换能器的孔径宽度D的调节,对不同焦距f下的活跃物理阵元数N_elements按照进行动态调整,N_elements取整数,且当阵元总数为偶数时,活跃物理阵元数N_elements取偶数,当阵元总数为奇数时,活跃物理阵元数N_elements取奇数。
3.根据权利要求1所述的基于动态孔径控制的超声剪切波弹性成像方法,其特征在于,步骤S2中声辐射力式中,α为生物组织的声衰减系数,c为超声在生物组织中的传播速度,I为一空间位置时间平均声强,且式中,ρ为生物组织的组织密度,Z为生物组织的声阻抗,Pe为当前时刻的有效声压。
4.根据权利要求1所述的基于动态孔径控制的超声剪切波弹性成像方法,其特征在于,步骤S3中按照求得该声辐射力引起质点振动的位移场,式中,为瞬态作用力矢量在粘弹性介质中引起的位移场,和为方向向量,t和τ为时间,为格林函数分量。
5.根据权利要求1所述的基于动态孔径控制的超声剪切波弹性成像方法,其特征在于,步骤S4中的峰值时间法求解剪切力波的传播速度CT的具体步骤如下:
S41:在聚焦位置的横向方向上,取间距相同的多个标记点,记录每个采样时刻对应的振动位移,得到每个标记点的振动位移时间图;
S42:记录每个标记点振动位移主峰值对应的时刻;
S43:利用最小二乘法拟合出以各标记点的位置信息为自变量,对应的主峰值到达的时刻为因变量的一条直线,该条直线的斜率的倒数即为剪切波传播速度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510794661.0A CN105232085B (zh) | 2015-11-18 | 2015-11-18 | 基于动态孔径控制的超声剪切波弹性成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510794661.0A CN105232085B (zh) | 2015-11-18 | 2015-11-18 | 基于动态孔径控制的超声剪切波弹性成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105232085A CN105232085A (zh) | 2016-01-13 |
CN105232085B true CN105232085B (zh) | 2017-10-31 |
Family
ID=55030113
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510794661.0A Expired - Fee Related CN105232085B (zh) | 2015-11-18 | 2015-11-18 | 基于动态孔径控制的超声剪切波弹性成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105232085B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104640506B (zh) * | 2013-09-18 | 2017-06-30 | 东芝医疗系统株式会社 | 超声波诊断装置、医用图像处理装置和医用图像处理方法 |
CN108135568B (zh) * | 2015-10-08 | 2021-11-26 | 梅约医学教育与研究基金会 | 用于利用持续换能器振动进行超声弹性成像的系统和方法 |
CN106383996A (zh) * | 2016-09-05 | 2017-02-08 | 四川大学 | 肝脏组织弹粘性力学新模型的建立以及基于新模型的应力应变关系和剪切波速度的获得方法 |
CN106618638B (zh) * | 2016-11-04 | 2019-02-26 | 声泰特(成都)科技有限公司 | 一种定量剪切波弹性成像系统 |
CN107440740B (zh) * | 2017-07-21 | 2021-06-25 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性定量方法及装置 |
CN107505233A (zh) * | 2017-07-21 | 2017-12-22 | 无锡海斯凯尔医学技术有限公司 | 介质粘弹性定量方法及装置 |
JP6996035B2 (ja) * | 2017-11-02 | 2022-01-17 | 富士フイルムヘルスケア株式会社 | 超音波診断装置、および、生体組織の物性評価方法 |
WO2020037674A1 (zh) * | 2018-08-24 | 2020-02-27 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像装置和方法、超声弹性检测装置和方法 |
CN109730722B (zh) * | 2019-02-27 | 2021-04-27 | 天津大学 | 基于聚焦超声声振信号的弹性成像方法 |
CN109745077B (zh) * | 2019-02-27 | 2021-04-27 | 天津大学 | 基于聚焦超声声振信号的弹性特性检测方法 |
CN112034427B (zh) * | 2020-09-09 | 2023-10-27 | 京东方科技集团股份有限公司 | 一种相控阵多焦点测量的方法、装置及存储介质 |
WO2022141052A1 (zh) * | 2020-12-29 | 2022-07-07 | 无锡祥生医疗科技股份有限公司 | 剪切波弹性成像方法及系统 |
CN113827278B (zh) * | 2021-10-27 | 2023-09-15 | 青岛海信医疗设备股份有限公司 | 剪切波传播速度的确定方法及装置 |
CN114403920A (zh) * | 2022-01-24 | 2022-04-29 | 陕西师范大学 | 一种分时多点激励剪切波融合成像方法及系统 |
CN114947963A (zh) * | 2022-06-17 | 2022-08-30 | 中国医学科学院生物医学工程研究所 | 一种眼轴测量方法、装置、存储介质和计算机设备 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5584157A (en) * | 1978-12-18 | 1980-06-25 | Picker Corp | Automatic aperture control for ultrasonic wave |
DE3311040A1 (de) * | 1983-03-25 | 1984-09-27 | Siemens AG, 1000 Berlin und 8000 München | Ultraschall-sende/empfangs-system |
US7662114B2 (en) * | 2004-03-02 | 2010-02-16 | Focus Surgery, Inc. | Ultrasound phased arrays |
CN101690677B (zh) * | 2009-10-13 | 2011-09-21 | 西安交通大学 | 聚焦超声治疗组合阵元相控阵及多焦点剪切波成像系统 |
CN101893705B (zh) * | 2010-06-30 | 2013-02-27 | 重庆大学 | 基于超声成像系统的动态孔径控制方法 |
-
2015
- 2015-11-18 CN CN201510794661.0A patent/CN105232085B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105232085A (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105232085B (zh) | 基于动态孔径控制的超声剪切波弹性成像方法 | |
CN102641137B (zh) | 使用幅度-相位调制超声波的粘弹性测量 | |
AU2016347143B2 (en) | Elasticity detection method and device | |
JP5771758B2 (ja) | 超音波診断装置 | |
CN103462643B (zh) | 一种剪切波速度测量方法、装置和系统 | |
JP2010069295A (ja) | 剛性が変化した領域を検出するためのシステム及び方法 | |
KR101922522B1 (ko) | 전단파들을 사용하는 사운드 스피드 이미징 | |
JP2014004375A5 (zh) | ||
Radulescu et al. | Hydrophone spatial averaging corrections from 1 to 40 MHz | |
JP2013140029A (ja) | 液面レベル計測装置、方法及びプログラム | |
WO2008135659A3 (fr) | Dispositif pour mesurer des propriétés viscoélastiques de tissus biologiques et procédé utilisant ce dispositif. | |
US20120250454A1 (en) | Method and system for shaping a cmut membrane | |
WO2015128656A1 (en) | Methods and systems for measuring properties with ultrasound | |
CN102539530A (zh) | 基于时间反转聚焦的水声无源材料回声降低/反射系数测量方法 | |
US9642600B2 (en) | Shear wave attenuation from k-space analysis system | |
CN104605890A (zh) | 一种剪切波峰值波形校正的方法、装置、系统及其应用 | |
CN107085041B (zh) | 一种聚焦声场中声偏振方向的控制方法 | |
Varslot et al. | Aberration in nonlinear acoustic wave propagation | |
JP6420574B2 (ja) | 超音波診断装置及びプログラム | |
JP5059344B2 (ja) | 板厚測定装置および測定方法 | |
JP5895571B2 (ja) | 超音波診断装置 | |
KR101117547B1 (ko) | 평면파를 이용하여 탄성영상을 형성하는 초음파 진단 시스템 및 방법 | |
JP5949090B2 (ja) | 超音波送信ビーム形成制御に用いられる方法、装置、およびソフトウェアプログラム | |
EP0605720A1 (en) | Inspection apparatus | |
JP5393568B2 (ja) | 超音波撮像方法および超音波撮像装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171031 Termination date: 20201118 |
|
CF01 | Termination of patent right due to non-payment of annual fee |