[0013] 3)由L组阵元向颅内成像区域依次发射一束低频线性调频聚焦扫描声束,并接收来自成像区域范围的超声回波信号;
[0014] 4)利用开放式超声成像系统或数据采集卡对每组阵元接收到的超声回波信号进行采样并转化为射频数据;
[0015] 5)针对每组阵元扫描声束采样并转化得到的超声回波信号射频数据ri(t),根据发射波形设计压缩滤波器h(-t),对各超声回波信号进行脉冲压缩,脉冲压缩结果为ri(t)*h(-t),其中,i=1,...,N,*表示卷积;
[0016] 6)改变线性调频激励信号为-h(t)(π相位),然后重复步骤3)-5);
[0017] 7)将脉冲压缩后的不同相位的(0相位和π相位)信号按照各阵元对应位置进行叠加,得到脉冲逆转叠加后的回波信号
[0018] 8)根据扫描声束内目标点位置,计算阵元接收到的回波信号的时间延迟,然后对脉冲逆转叠加后的回波信号进行时间校正处理,得到延迟后的射频信号rdi(t);
[0019] 9)根据波束合成方法对延迟后的射频信号rdi(t)进行变迹,然后加权叠加:
[0020]
[0021] 其中ωi为各阵元对应的加权值,加权值选取对应波束合成方法的迹函数计算,即合成该目标点数据;
[0022] 10)重复步骤8)-9),对每个扫描声束下的所有目标点脉冲逆转叠加后的信号进行波束合成,得到包含L条扫描线的整幅成像数据;
[0023] 11)对经过波束合成之后的整幅成像数据依次进行包络检波、对数压缩,得到输出图像。
[0024] 优选的,所述时间校正处理包括以下步骤:设当前扫描声束下目标点坐标为(x,0,z),接收阵元的回波信号相对于参考阵元的时间延迟为:
[0025]
[0026] 其中,xi对应每个接收阵元的坐标值,i=1,...,N,c为声速,则每个阵元对应的延迟后的脉冲逆转叠加射频信号为:
[0027]
[0028] 优选的,所述超声线阵换能器的中心频率选取范围为1.5MHz~3MHz,带宽为60%~80%。
[0029] 优选的,所述线性扫描线数选取为L=64~256。
[0030] 优选的,所述波束合成方法选取压缩自适应波束合成方法(万明习,柏晨,徐珊珊,敬博文.“压缩自适应波束合成的平面波超声成像和微泡成像的方法与系统”.申请号:201510079932.4)。
[0031] 一种经颅低频超声平面波单一角度或多角度复合线性调频脉冲逆转微泡成像方法,包括以下步骤:
[0032] 1)将低频经颅专用超声线阵换能器(阵元总数为M)放置于一侧颞骨附近;
[0033] 2)所述超声线阵换能器的M个阵元同时向颅内成像区域发射与阵元平面呈倾斜角度为αk的宽波束(规定向正半轴方向偏转角度为正),k=1,...,n,n为奇数,保证倾斜角度在正负半轴对称且α(n+1)/2=0,当n=1是即为单一角度成像,发射所用的线性调频激励信号为h(t)(假设为0相位),h(t)的中心频率为所述超声线阵换能器的中心频率,时宽及带宽满足所述超声线阵换能器频带特性要求;
[0034] 3)利用开放式超声成像系统或数据采集卡对M个阵元接收的超声回波信号进行采样并转化为射频数据ri(t);
[0035] 4)根据发射波形设计压缩滤波器h(-t),对ri(t)进行脉冲压缩,脉冲压缩结果为ri(t)*h(-t),其中,i=1,...,M,*表示卷积;
[0036] 5)在当前倾斜角度下改变激励信号为-h(t)(π相位),然后重复步骤2)-4);
[0037] 6)将脉冲压缩后的不同相位(0相位和π相位)的信号按照各阵元对应位置进行叠加,得到脉冲逆转叠加数据
[0038] 7)根据成像区域内目标点位置及倾斜角度αk计算当前倾斜角度下阵元接收到的回波信号的时间延迟,然后对当前倾斜角度下的脉冲逆转叠加数据进行时间校正处理,得到延迟后的射频信号rdi(t);
[0039] 8)对延迟后的射频信号rdi(t)进行波束合成,求得变迹加权叠加输出:
[0040]
[0041] 其中ωi为各阵元对应的加权值,加权值选取对应波束合成方法的迹函数计算;遍历成像区域内目标点,得到当前倾斜角度下的波束合成后数据;
[0042] 9)改变αk,重复步骤2)-8),直至得到所有倾斜角度下的波束合成后数据,将每个倾斜角度下的波束合成后数据叠加后再平均,获得最终的成像数据;
[0043] 10)对成像数据依次进行包络检波、对数压缩,得到输出图像。
[0044] 优选的,所述宽波束为所有阵元同时发射、同时接收数据。
[0045] 优选的,所述时间校正处理包括以下步骤:设目标点位置为(x,0,z),回波信号接收阵元位置为(xi,0,0),则接收阵元的回波信号相对于参考阵元的时间延迟为:
[0046]
[0047] 其中,i=1,...,M,对应每个接收阵元,c为声速,则每个阵元对应的延迟后的脉冲逆转叠加射频信号为:
[0048]
[0049] 优选的,采用n=3~7个倾斜角度的宽波束进行复合成像。
[0050] 优选的,所述波束合成方法选取压缩自适应波束合成方法(万明习,柏晨,徐珊珊,敬博文.“压缩自适应波束合成的平面波超声成像和微泡成像的方法与系统”.申请号:201510079932.4)。
[0051] 一种经颅超声微泡成像系统,包括低频经颅专用超声线阵换能器以及回波信号接收模块、脉冲逆转叠加模块、波束合成模块和成像模块;其中:
[0052] 所述超声线阵换能器用于发射相位相反的线性调线性扫描声束或宽波束;
[0053] 所述回波信号接收模块用于对所述超声线阵换能器对应阵元接收的超声回波信号进行采样并转化为射频数据;
[0054] 所述脉冲逆转叠加模块用于对射频数据进行脉冲压缩以及将两组脉冲压缩后的相位相反的射频数据进行叠加;
[0055] 所述波束合成模块用于将经过叠加的射频数据经延迟校正后进行压缩自适应波束合成;
[0056] 所述成像模块用于对经过波束合成得到的信号依次进行包络检波及对数压缩。
[0057] 本发明的有益效果体现在:
[0058] 本发明不同于目前主要针对小动物颅脑的高频超声经颅成像方法,提出了利用低频线性调频信号与线性扫描模式或多角度复合成像技术联合,突破了颅骨对超声信号的遮挡和屏蔽,实现了高灵敏度的颅内组织超声成像。本发明解决了超声信号经颅骨后严重衰减,使得经颅超声探测深度有限、难以对颅内进行检测成像的难题。
附图说明
[0059] 图1为实施例中给出的经颅超声微泡成像系统示意图;
[0060] 图2为经颅超声微泡成像实验中所用声学及结构特性匹配的3D打印颅骨模型;
[0061] 图3为线性调频发射激励信号h(t)的时域图;
[0062] 图4为低频超声线扫描模式线性调频各扫描波束声场示意图(Beam指扫描线);
[0063] 图5为低频超声线性扫描成像结果,其中:(a)经颅低频超声常规线性扫描(0相位)成像,(b)经颅低频超声线性扫描模式线性调频(0相位)成像,(c)经颅低频超声常规线性扫描脉冲逆转微泡成像,(d)经颅低频超声线性扫描模式线性调频脉冲逆转微泡成像,图中实/虚线框内信号强度用来计算造影组织比(CTR);
[0064] 图6为各低频超声线性扫描成像方法CTR比较,其中:(a)经颅低频超声常规线性扫描(0相位)成像,(b)经颅低频超声线性扫描模式线性调频(0相位)成像,(c)经颅低频超声常规线性扫描脉冲逆转微泡成像,(d)经颅低频超声线性扫描模式线性调频脉冲逆转微泡成像;
[0065] 图7为5角度下经颅低频超声平面波线性调频发射声场示意图;其中:(a)-10°角度声场示意图,(b)-5°角度声场示意图,(c)0°角度声场示意图,(d)5°角度声场示意图,(e)10°角度声场示意图;
[0066] 图8为超声宽波束改变角度发射后,其回波信号时间延迟计算示意图;
[0067] 图9为低频超声平面波成像结果,其中:(a)经颅低频超声平面波单一角度脉冲逆转微泡成像,(b)经颅低频超声平面波多角度复合脉冲逆转微泡成像,(c)经颅低频超声平面波单一角度线性调频脉冲逆转微泡成像,(d)经颅低频超声平面波多角度复合线性调频脉冲逆转微泡成像,图中实/虚线框内信号强度用来计算造影组织比(CTR);
[0068] 图10为各低频超声平面波成像方法CTR比较,其中:(a)经颅低频超声平面波单一角度脉冲逆转微泡成像,(b)经颅低频超声平面波多角度复合脉冲逆转微泡成像,(c)经颅低频超声平面波单一角度线性调频脉冲逆转微泡成像,(d)经颅低频超声平面波多角度复合线性调频脉冲逆转微泡成像。
具体实施方式
[0069] 下面结合附图和实施例对本发明做进一步详细说明。
[0070] 以低频经颅专用超声线阵换能器为例,超声线阵换能器阵元总数为M=128,各换能器阵元尺寸为0.25mm×10mm,换能器阵元间距为0.05mm,发射中心频率为2MHz,有效带宽为70%。
[0071] 如图1所示,本发明所举例实验中,低频经颅专用超声线阵换能器在开放式超声成像系统主机的控制下发射超声线性扫描波束/单一或多角度平面波、接收回波信号并将回波信号传送至主机;主机将回波信号传输至开放式超声成像系统接收模块;接收模块对收到的回波信号采样、存储为射频采样数据并传送给计算机;计算机对回波射频采样数据进行脉冲压缩、脉冲逆转、波束合成等处理,随后依照标准的成像过程,包括射频信号包络检测以及对数压缩等,显示得到超声图像。造影微泡按照血管仿体的管道以流动泵推动进行持续流动,按照前述方法对流动的造影微泡进行成像。
[0072] 如图2所示,本发明所举例实验中,低频经颅专用超声线阵换能器都将放置于颞骨附近(此区域对超声波的声衰减最弱);所用颅骨为声学及结构特性与真实颅骨匹配的3D打印颅骨模型,其声学衰减特性、尺寸大小与人体真实颅骨几乎一致,但可以用来任意切割以方便实验,具体请参阅文献(C.Bai,M.Ji,J.Zong,et al.“, A 3D-printed Skull Model with Corresponding Acoustic Characteristic of Human Skull for Ultrasound Brain Imaging and Diagnosis”,Proceedings of International Society for Therapeutic Ultrasound(ISTU)17th Annual Symposium,Nanjing,China,Jun.,2017.)。
[0073] 如图3所示,本发明中超声成像线阵激励信号采用上扫频信号,其中心频率2MHz,时宽7us,带宽1.3MHz,满足低频经颅专用超声线阵换能器频带特性。
[0074] (一)经颅低频超声线性扫描模式线性调频脉冲逆转微泡成像方法[0075] 超声换能器的功能是将输入的电功率转换成机械功率(即超声波)再传递出去,而自身消耗很少的一部分功率。由于经颅传播过程中颅骨对超声波有较大的衰减,而介质对超声信号的衰减与信号频率呈反比,为了使颅内微泡成像更加清晰,利用特定的低频超声换能器发射低频超声信号,可以在一定程度上减少颅骨的衰减作用,从而增加回波信号的强度。为了进一步提高经颅超声微泡成像的灵敏度,本发明在经颅低频超声成像的基础上,提出低频超声线性扫描模式线性调频脉冲逆转微泡成像方法,即利用线性调频信号有效提高成像的信噪比及灵敏度。
[0076] 经颅低频超声线性扫描模式线性调频脉冲逆转微泡成像方法的具体步骤为:
[0077] 1)低频经颅专用超声线阵换能器放置于一侧颞骨附近;
[0078] 2)将N个阵元编成一组,共L=128组。线性调频激励信号为图3所示的h(t)(0相位);
[0079] 3)由一组阵元向成像的微泡区域发射一束低频线性调频聚焦扫描声束(图4),聚焦深度为50mm,并接收来自成像区域范围的回波信号;
[0080] 4)利用开放式超声成像系统对接收到的超声回波信号进行采样并转化为射频数据;
[0081] 5)依次由下一组阵元发射相同的0相位声束,并接收、采样、转化为射频数据;
[0082] 6)在所有组阵元发射接收完成后,将得到的所有扫描波束的射频数据存储并传送给计算机;
[0083] 7)针对每组发射波束下的回波信号ri(t),根据发射波形设计压缩滤波器h(-t),对各回波信号进行脉冲压缩,其结果为ri(t)*h(-t),其中i=1,...,N,*表示卷积;
[0084] 8)改变激励信号为-h(t)(π相位),重复步骤3)-7);
[0085] 9)将脉冲压缩后的0相位、π相位信号按照阵元对应位置进行叠加,即脉冲逆转叠加;
[0086] 10)根据目标点位置,计算回波信号的时间延迟,并对脉冲逆转叠加后的信号进行时间校正:设当前扫描波束下目标点坐标为(x,0,z),计算回波信号相对于参考阵元的时间延迟:
[0087]
[0088] 其中,xi对应每个接收阵元的坐标值,其中i=1,...,N,则每个阵元对应的延迟后的脉冲逆转叠加射频信号为:
[0089]
[0090] 11)根据压缩自适应波束合成方法对延迟后的信号rdi(t)进行变迹,加权叠加输出的信号可以表示为:
[0091]
[0092] 其中ωi为各阵元对应的加权值。遍历当前扫描波束下目标点,即合成一条扫描线。
[0093] 12)重复步骤10)-11),对每个扫描声束下的回波信号进行波束合成,得到128条扫描线;
[0094] 13)对经过波束合成之后128条扫描线数据进行包络检波、对数压缩,最后输出一帧图像。
[0095] 如图5所示,通过比较在相同条件下的经颅超声常规线性扫描(0相位)、线性扫描模式下的线性调频(0相位)、常规线性扫描脉冲逆转、线性扫描模式下的线性调频脉冲逆转微泡成像的结果,以上所述经颅低频超声线性扫描模式线性调频脉冲逆转微泡成像方法得到的成像质量有显著提高(图5d),但是受到线性扫描模式的限制,该方法的成像帧率只有约100Hz。
[0096] 如图6所示,通过量化考察各方法成像的造影组织比(CTR),其中以上所述经颅低频超声线性扫描模式线性调频脉冲逆转微泡成像方法的CTR最高(图6d),约为常规线性扫描成像(图6a)的2.50倍,约为常规线性扫描脉冲逆转成像(图6c)的1.32倍。
[0097] (二)经颅低频超声平面波单一角度或多角度复合线性调频脉冲逆转微泡成像方法
[0098] 低频经颅专用超声线阵换能器的所有阵元同时发射一次超声宽波束,可以覆盖成像范围内的颅内微泡目标所在的大部分区域,从而一次性获取整个微泡区域的射频数据。在垂直于声束的方向上,每一部分的图像数据都是同时采集到的,因此极大地避免了线扫描方式中扫面线之间的采样时间差,进而提高了成像速度,保证了成像的实时性,其单一成像帧率可以达到5000Hz,多角度复合成像帧率可以达到1000Hz。
[0099] 所述经颅低频超声平面波单一或多角度复合线性调频脉冲逆转微泡成像方法,具体步骤为:
[0100] 1)低频经颅专用超声线阵换能器放置于一侧颞骨附近;
[0101] 2)如图7所示,128个阵元同时向成像的微泡区域发射倾斜角度为αk的宽波束(规定向正半轴方向偏转角度为正,且k=1,...,5,当仅取倾斜角度αk=0°即为单一角度成像),发射所用的线性调频激励信号为图3所示的h(t)(0相位);
[0102] 3)发射声波遇到成像目标后发生散射,产生与发射方向相反的回波;
[0103] 4)利用开放式超声成像系统对超声回波信号进行采样并转化为射频数据;
[0104] 5)根据发射波形设计压缩滤波器h(-t),对各回波信号ri(t)进行脉冲压缩,其结果为ri(t)*h(-t),其中i=1,...,128,*表示卷积;
[0105] 6)在当前倾斜角度下改变激励信号为-h(t)(π相位),重复步骤2)-5);
[0106] 7)将脉冲压缩后的0相位、π相位信号按照各阵元对应位置进行叠加,即脉冲逆转叠加;
[0107] 8)如图8所示,假设目标点位置为(x,0,z),回波信号接收阵元位置为(xi,0,0),根据倾斜角度αk按照式(4)计算当前角度的时间延迟,对当前角度脉冲逆转叠加数据进行校正处理;
[0108]
[0109] 其中,i=1,...,128,对应每个接收阵元,c为声速。则每个阵元的延迟校正后的回波信号为:
[0110]
[0111] 9)对校正后的射频数据rdi(t),依据式(6)进行压缩自适应波束合成,求得变迹加权叠加输出:
[0112]
[0113] 其中ωi为各阵元对应的加权值;遍历成像区域内目标点,得到当前角度下的波束合成后数据。
[0114] 10)改变αk,重复步骤2)-9),将所有角度的波束合成后数据连续叠加并平均,获得最终的成像数据;
[0115] 11)对成像数据经包络检波、对数压缩后输出一帧图像。
[0116] 如图9所示,通过比较在相同条件下的经颅低频超声平面波单一角度脉冲逆转、平面波多角度复合脉冲逆转、平面波单一角度线性调频脉冲逆转、平面波多角度复合线性调频脉冲逆转超声微泡成像的结果,经颅低频超声平面波多角度复合线性调频脉冲逆转微泡成像方法的成像质量有显著提高(图9d),成像帧率有了极大提升。
[0117] 如图10所示,量化考察各方法成像的造影组织比(CTR),其中经颅低频超声平面波多角度复合线性调频脉冲逆转微泡成像方法的CTR最高(图10d),约为平面波单一角度成像(图10a)的3倍,约为平面波多角度复合脉冲逆转微泡成像(图10b)的2.18倍,约为平面波单一角度线性调频脉冲逆转(图10c)的1.43倍。