一种超声成像方法及装置
技术领域
本发明涉及超声成像领域,尤其涉及一种超声成像方法及装置。
背景技术
现有超声成像系统的处理流程一般如图1所示:换能器将电信号转换为超声波,向人体发射声波束;换能器接收经过组织反射和散射的超声波束,转换为电信号;经过波束合成120,信号处理130后,将处理后的数据进行DSC(DigitalScanConversion,数字扫描变换)变换140和图像处理150;通过显示模块160显示处理后的图像。
介入性超声技术是近年来随着超声医学发展而迅速发展起来的分支科学,主要是利用超声波传导束具有良好的指向性、高分辨率的二维图像显示、能够实时监测和无创性等技术优点,将超声作为穿刺定位的工具,在超声引导下(导向),将穿刺针、引流导管、输液管或药物等正确地插入或者注入所要达到的病灶、囊腔、体腔、管道和其他特定部位,间接达到使操作者能够在直视病灶或靶器官所在方位的情况下进行操作的效果,提高穿刺的准确性和安全性,从而根据穿刺目的和临床需要进行诊断或治疗。
以穿刺针为例。在进行穿刺或麻醉时,穿刺针和探头表面呈一定的角度刺入人体组织。由于穿刺针的声阻抗大,超声波无法穿透穿刺针,从而形成一个反射界面。一般地,穿刺针的部位是一个强反射的回声,传统的超声成像系统的B模式沿垂直于探头表面的方向发射、接收。如图2所示,探头10发射的声束遇到穿刺针20后反射,根据反射定律,由于声束方向发生了改变,大部分的能量无法在探头的相同位置接收到。尤其是在穿刺针刺入的偏转角较大时,超声B图像上针的显示会变弱,不利于使用者的观察,从而不利于进行穿刺定位操作。
发明内容
本发明实施方式要解决的技术问题是,提供一种超声成像方法及装置,能够既增强显示如穿刺针等介入性物体,又保证了图像的质量不会变差。
根据本发明实施例的第一方面,提供一种超声成像方法,包括:发射步骤,从至少一个第一角度和至少一个第二角度向待检测对象发射超声波束,所述第一角度为有利于接收倾斜插入待检测对象的介入性物体的超声回波的角度,所述第二角度为有利于接收待检测对象内部组织的超声回波的角度;接收步骤,接收从所述至少一个第一角度和所述至少一个第二角度反馈的回波信号,包括至少一个第一角度回波信号和至少一个第二角度回波信号;信号处理步骤,对接收到的回波信号进行信号处理后输出图像,包括至少一个第一角度图像和至少一个第二角度图像;显示步骤,对输出的图像进行显示。
根据本发明实施例的第二方面,提供一种超声成像方法,包括:第一角度图像获取步骤,从多个角度向插入介入性物体的待检测对象发射超声波束,接收从多个角度返回的回波信号,对回波信号进行处理后得到多个角度的图像,根据图像的特征,从该多个图像中选择出有利于显示倾斜插入待检测对象的介入性物体的至少一个图像,将选择出的图像作为第一角度图像输出;第二角度图像获取步骤,从至少一个第二角度向待检测对象发射超声波束,所述第二角度为有利于接收待检测对象内部组织的超声回波的角度,接收从所述至少一个第二角度反馈的至少一个第二角度回波信号,对接收到的至少一个第二角度回波信号进行信号处理后输出至少一个第二角度图像;显示步骤,对输出的图像进行显示。
根据本发明实施例的第三方面,提供一种超声成像装置,包括:发射模块,用于从至少一个第一角度和至少一个第二角度向待检测对象发射超声波束,所述第一角度为有利于接收倾斜插入待检测对象的介入性物体的超声回波的角度,所述第二角度为有利于接收待检测对象内部组织的超声回波的角度;接收模块,用于接收从所述至少一个第一角度和所述至少一个第二角度反馈的回波信号,包括至少一个第一角度回波信号和至少一个第二角度回波信号;信号处理模块,用于对接收到的回波信号进行信号处理后输出图像,包括至少一个第一角度图像和至少一个第二角度图像;显示模块,用于对输出的图像进行显示。
根据本发明实施例的第四方面,提供一种超声成像装置,包括:第一角度图像获取模块,用于从多个角度向插入介入性物体的待检测对象发射超声波束,接收从多个角度返回的回波信号,对回波信号进行处理后得到多个角度的图像,根据图像的特征,从该多个图像中选择出有利于显示倾斜插入待检测对象的介入性物体的至少一个图像,将选择出的图像作为第一角度图像输出;第二角度图像获取模块,用于从至少一个第二角度向待检测对象发射超声波束,所述第二角度为有利于接收待检测对象内部组织的超声回波的角度,接收从所述至少一个第二角度反馈的至少一个第二角度回波信号,对接收到的至少一个第二角度回波信号进行信号处理后输出至少一个第二角度图像;显示模块,用于对输出的图像进行显示。
本发明实施例的有益效果在于:将从第一角度向介入性物体的发射和从第二角度向待检测对象的发射交替进行,既保证显示出介入性物体形成的图像,又保证待检测对象形成的图像的质量不变差。
附图说明
图1是现有超声成像装置的原理性示意框图;
图2是探头垂直组织发射时声波在穿刺针的界面反射的原理示意图;
图3是探头垂直或基本垂直穿刺针发射时声波在穿刺针界面反射的原理示意图;
图4是本发明超声成像方法(装置)中探头分别向介入性物体和待检测对象发射超声波束的方向示意图;
图5(a)、(b)、(c)、(d)分别是本发明超声成像装置四种实施方式的原理性示意图;
图6(a)、(b)、(c)、(d)分别是本发明超声成像方法四种实施方式的流程示意图;
图7(a)、(b)分别是本发明超声成像方法一种实施方式中空间复合关闭与开启的示意图;
图8是本发明超声成像方法一种实施方式中合成处理步骤的流程示意图;
图9(a)、(b)分别是Hough变换的空间示意图;
图10(a)、(b)分别是本发明超声成像方法实施方式中带权重Hough变换的流程示意图;
图11(a)、(b)分别是穿刺针在偏转发射、正常发射中成像的示意图;
图11(c)是将正常图像和偏转图像叠加后的示意图;
图12是正常图像与偏转图像叠加的算法流程示意图;
图13是合成区域中进行投影的示意图;
图14是在合成图像中加框后显示的示意图;
图15是本发明超声成像方法另一种实施方式中合成处理步骤的流程示意图。
具体实施方式
下面通过具体实施方式结合附图对本发明作进一步详细说明。
首先对下面用到的一些术语或概念进行解释。
(1)倾斜角度:是指介入性物体倾斜插入待检测对象时所形成的倾斜角,即介入性物体和待检测对象表面之间的夹角,如图4所示的
(2)偏转角度:当介入性物体倾斜插入待检测对象时,超声波束的发射方向垂直或基本垂直于介入性物体表面时,如图3所示,介入性物体反馈的回波最强,基本垂直指超声波束的发射方向与介入性物体表面的夹角为90°±φ,φ为一较小的角度,例如5°或10°。为获得较强的介入性物体的超声回波,探头向介入性物体发射的超声波束与垂直于探头表面的超声波束之间具有一定范围的夹角,如图4所示的这里将称为第一角度,或称为偏转角度。
(3)偏转发射:是指探头从第一角度向待检测对象(如病灶等特定部位)发射超声波束。第一角度是有利于接收倾斜插入待检测对象的介入性物体(如穿刺针)的超声回波的角度,也就是说,在偏转发射时,探头向介入性物体发射的超声波束垂直或基本垂直于介入性物体。根据基本的几何知识可知,此时等于或基本等于
(4)偏转图像:是指从第一角度接收到的超声回波经信号处理后所得的第一角度图像。
(5)正常角度:是指探头向待检测对象发射超声波束的方向与垂直于探头表面的方向所形成的夹角如图4所示,该夹角也称为第二角度。应理解,该夹角可以为0°(即探头发射的超声波束垂直于待检测对象),也可以为一锐角。
(6)正常发射:是指探头从第二角度向待检测对象发射超声波束。探头从第二角度发射超声波束是为了获得有利于接收待检测对象内部组织的超声回波。若探头从两个或两个以上第二角度向待检测对象发射超声波束,则,相当于是沿不同角度对待检测对象进行扫描,然后再将这些不同角度的图像对应的像素点进行叠加(即对图像的整个区域进行处理),形成一帧图像,这种成像方法称为空间复合成像。空间复合成像的目的是为了降低噪声,增强图像不同组织间的对比。这里将探头从多个不同的第二角度进行发射以成像的方式称为空间复合开启;若探头仅从一个第二角度进行发射以形成一帧图像,这种方式称为空间复合关闭。
(7)正常图像:是指从第二角度接收到的超声回波经信号处理后所得的第二角度图像,包括空间复合开启所形成的一帧图像,也包括空间复合关闭时所得到的一帧图像。
本发明实施方式的总体思想是:将从第一角度向介入性物体的发射和从第二角度向待检测对象的发射交替进行,从而,既保证显示出介入性物体形成的图像,又保证待检测对象形成的图像的质量不变差。
图5(a)是本发明超声成像装置一种实施方式的原理性示意图,如图所示,包括:发射模块511、接收模块513、信号处理模块515和显示模块519。该发射模块511用于从至少一个第一角度和至少一个第二角度向待检测对象发射超声波束,所述第一角度为有利于接收倾斜插入待检测对象的介入性物体的超声回波的角度,也就是说第一角度用于对插入待检测对象中的介入性物体进行成像,所述第二角度为有利于接收待检测对象内部组织的超声回波的角度,也就是说第二角度用于对待检测对象进行成像。接收模块513用于接收从至少一个第一角度反馈的第一角度回波信号以及接收从至少一个第二角度反馈的至少一个第二角度回波信号。信号处理模块515用于对接收到的至少一个第一角度回波信号进行信号处理后输出至少一个第一角度图像以及对接收到的至少一个第二角度回波信号进行信号处理后输出至少一个第二角度图像。显示模块519用于对输出的图像进行显示。
图5(a)的一种变形实施方式包括:第一角度成像模块、第二角度成像模块和显示模块。该第一角度成像模块用于以第一角度对待检测对象发射超声波束,接收待检测对象反馈的回波信号,并进行信号处理获得第一角度图像并输出,其中,该第一角度根据插入待检测对象中的介入性物体的倾斜角度获得。该第二角度成像模块用于以第二角度对待检测对象发射超声波束,接收待检测对象反馈的回波信号,并进行信号处理获得第二角度图像并输出。显示模块用于对输出的图像进行显示。
图5(b)为在图5(a)所示实施方式基础上产生的变形实施方式,如图所示,该实施方式的超声成像装置包括发射模块521、接收模块523、信号处理模块525、合成处理模块527和显示模块529。所述发射模块521、接收模块523、信号处理模块525分别类似于图5(a)所示实施方式的发射模块511、接收模块513、信号处理模块515,合成处理模块527用于将第二角度图像与第一角度图像进行合成处理以生成新的合成图像并输出。同样地,图5(b)的一种变形实施方式包括:第一角度成像模块、第二角度成像模块、合成处理模块和显示模块,各模块的功能描述同前述,在此不再重述。
图5(c)为本发明超声成像装置另一种实施方式的原理性示意图,如图所示,包括:第一角度图像获取模块531、第二角度图像获取模块532、显示模块539;其中,第一角度图像获取模块531用于从多个角度向插入介入性物体的待检测对象发射超声波束,接收从多个角度返回的回波信号,对回波信号进行处理后得到多个角度的图像,根据图像的特征,从该多个图像中选择出有利于显示倾斜插入待检测对象的介入性物体的至少一个图像,将选择出的图像作为第一角度图像输出;第二角度图像获取模块532用于至少一个第二角度向待检测对象发射超声波束,第二角度为有利于接收待检测对象内部组织的超声回波的角度,接收从至少一个第二角度反馈的至少一个第二角度回波信号,对接收到的至少一个第二角度回波信号进行信号处理后输出至少一个第二角度图像;显示模块539用于对输出的图像进行显示。
图5(d)为在图5(c)所示实施方式基础上产生的变形实施方式,如图所示,该实施方式的超声成像装置包括第一角度图像获取模块541、第二角度图像获取模块542、合成处理模块547和显示模块549。该第一角度图像获取模块541、第二角度图像获取模块542、显示模块549分别类似于图5(c)所示实施方式的第一角度图像获取模块531、第二角度图像获取模块532、显示模块539,合成处理模块547用于将第二角度图像与第一角度图像进行合成处理以生成新的合成图像并输出。
基于图5(a)-图5(d)所述的各实施方式的超声成像装置,可得到对应的各实施方式的超声成像方法,详见后文。
基于上述图5(a)所示超声成像装置的超声成像方法实施方式的流程如图6(a)所示,包括:发射步骤S611、接收步骤S613、信号处理步骤S615、显示步骤S619。其中,发射步骤S611包括:从至少一个第一角度发射超声波束(即偏转发射)和从至少一个第二角度发射超声波束(即正常发射),偏转发射与正常发射交替进行;接收步骤S613包括:接收从至少第一角度反馈的回波信号、和接收从至少一个第二角度反馈的回波信号;信号处理步骤S615包括:对接收到的至少一个第一回波信号进行信号处理后输出第一角度图像(即偏转图像)、和对接收到的至少一个第二回波信号进行信号处理后输出第二角度图像(即正常图像);显示步骤S619用于对输出的图像进行显示。
同理,可以得到基于图5(b)所示超声成像装置的超声成像方法实施方式的流程,如图6(b)所示,包括:发射步骤S621、接收步骤S623、信号处理步骤S625、合成处理步骤S627、显示步骤S629。其中,发射步骤S621、接收步骤S623、信号处理步骤S625分别类似于图6(a)所示实施方式的发射步骤S611、接收步骤S613、信号处理步骤S615、显示步骤S619,合成处理步骤S627用于将正常图像与偏转图像进行合成处理以生成新的合成图像并输出。
为使介入性物体在偏转图像中的显示比较强,探头向介入性物体发射的角度(即第一角度)要垂直于或基本垂直于该介入性物体;由前述,第一角度的角度值为或基本为介入性物体向所述待检测对象插入时的倾斜角度值,因此,只要获知介入性物体向所述待检测对象插入时所形成的倾斜角度值,就可以得到偏转发射时第一角度的角度值(以下称为“偏转角度”)。可以采用以下几种方式之一来获得偏转角度。
获取偏转角度的方式一:在用户界面上提供若干个角度,由用户来选择偏转角度。采用这种偏转角度获取时,用户选择的角度通常也是用户插入介入性物体的倾斜角度,或用户插入介入性物体的倾斜角度接近于其选择的角度。
获取偏转角度的方式二:通过介入性物体所在的支撑架上的传感器来获取角度值。例如,在穿刺架的进针槽上放置一个传感器,通过该传感器获取进针的角度,从而,由该进针的角度即可得知偏转角度。可以理解,方式一和方式二为半自动方式。
获取偏转角度的方式三:探头从多个角度进行发射,根据图像的特征选取从该多个角度中选择出最合适的那个角度的角度值。以所采用图像的特征为灰度值为例,假设多个角度分别为ψa、ψb、ψc、ψd,这些角度对应得到的图像分别为Ia、Ib、Ic、Id,利用现有图像处理技术判断Ia、Ib、Ic、Id中哪个图像中的介入性物体的灰度值最大,对应该灰度值最大的图像(如Ic)所对应的角度(如ψc)为最合适的角度。图像的特征还可以是其他的如图像梯度等,具体选择合适角度的过程与图像的特征为灰度值类似,在此不再重述。其中,这多个角度的发射方式可以是选择每隔一段合适的时间(可凭实验给出该时间间隔)进行发射,或者是每次发射都从多个角度进行发射。可以理解,方式三为自动方式。
回到图6,基于上述图5(c)所示超声成像装置的超声成像方法实施方式的流程如图6(c)所示,包括:第一角度图像获取步骤S631、第二角度图像获取步骤S632、显示步骤S639;其中,第一角度图像获取S631包括:从多个角度向插入介入性物体的待检测对象发射超声波束,接收从多个角度返回的回波信号,对回波信号进行处理后得到多个角度的图像,根据图像的特征,从该多个图像中选择出有利于显示倾斜插入待检测对象的介入性物体的至少一个图像,将选择出的图像作为第一角度图像输出;第二角度图像获取步骤S632用于至少一个第二角度向待检测对象发射超声波束,第二角度为有利于接收待检测对象内部组织的超声回波的角度,接收从至少一个第二角度反馈的至少一个第二角度回波信号,对接收到的至少一个第二角度回波信号进行信号处理后输出至少一个第二角度图像;显示步骤S639用于对输出的图像进行显示。
这里,对于步骤S631中提及的从多个角度进行发射是指,假设多个角度分别为ψa、ψb、ψc、ψd,这些角度对应得到的图像分别为Ia、Ib、Ic、Id,利用现有图像处理技术判断Ia、Ib、Ic、Id中哪个图像中的介入性物体的灰度值较大,对应该灰度值较大的图像(如Ic)所对应的角度(如ψc)为最合适的角度(也就是第一角度,或称偏转角度),而该图像(Ic)即为第一角度图像。其中,这多个角度的发射方式可以是选择每隔一段合适的时间(可凭实验给出该时间间隔)进行发射,或者是每次发射都从多个角度进行发射。
同样可以得到基于图5(d)所示超声成像装置的超声成像方法实施方式的流程,如图6(d)所示,包括:第一角度图像获取步骤S641、第二角度图像获取步骤S642、合成处理步骤S647、显示步骤S649;其中,第一角度图像获取步骤S641、第二角度图像获取步骤S642、显示步骤S649分别类似于图5(c)所示实施方式的第一角度图像获取模块531、第二角度图像获取步骤S632、显示步骤S639,合成处理步骤S647用于将第二角度图像与第一角度图像进行合成处理以生成新的合成图像并输出。
一般地,偏转发射的深度和正常发射的深度相同。实际操作中还可以对偏转发射的发射频率使用低频、物理(或称为透镜)方向弱聚焦,从而得到合适的偏转图像。
各实施方式中,进行偏转发射的目的是使对介入性物体的显示较强,改善垂直扫描时介入性物体在图像中的显示变弱,然而,偏转后的图像仍存在缺陷:一是探头下方的图像的有效区域变小,二是偏转后得到的图像除了针的部分显示增强外,组织的分辨率变差,偏转角越大,图像质量则越差,影响使用者对正常组织和病变组织的区分;因此,各实施方式中还增加了正常发射,正常发射可以得到待检测对象的正常图像,不会因受偏转发射影响而使得图像的有效区域变小。这样,通过将偏转发射和正常发射交替进行,则可达到既增强显示介入性物体形成的图像,又保证了正常图像的有效区域不会变小,图像质量不会变差。
下面给出更为具体的实施例以对本发明超声成像方法与装置进行进一步说明。
实施例一:
本实施例中第一角度只有一个,第二角度也只有一个,也就是说,探头仅从一个第一角度和一个第二角度进行发射,此时空间复合关闭,如图7(a)所示,探头70不仅从第一角度ψ1(如图所示的虚线)进行偏转发射,还从第二角度ψ2(如图所示的实线)进行正常发射)。
本实施例采用图6(a)所示流程,在发射步骤S611中,确定了偏转角度后,偏转发射和正常发射交替进行,可以先进行正常发射再进行偏转发射,或者反过来,即先进行偏转发射再进行正常发射。
接收步骤S613中,接收偏转发射和正常发射反馈的回波信号,包括偏转回波信号(即第一方向回波信号)和正常回波信号(即第二方向回波信号)。
信号处理步骤S615中,对接收到的回波信号进行信号处理后输出偏转图像(即第一角度图像)和正常图像(即第二角度图像),其中,信号处理可采用现有超声成像中涉及的信号处理,如模数转换、波束合成、正交调解、壁滤波、空间平滑等。
在显示步骤S619中,对信号处理步骤S615输出的图像进行显示,即显示正常图像和偏转图像,由用户自行对比两幅图以及感知介入性物体在图像中的位置。本领域技术人员应理解,在显示前可以采用DSC分别对两幅图像处理后再显示。
实施例二:
本实施例中第一角度和第二角度仍分别为一个,类似于实施例一,空间复合仍然关闭。
本实施例采用图6(b)所示流程,其中,发射步骤S621、接收步骤S623、信号处理步骤S625分别类似于实施例一中图6(a)的发射步骤S611、接收步骤S613、信号处理步骤S615,在此不再重述。
在合成处理步骤S627中,将正常图像和偏转图像进行合成处理以生成新的合成图像。具体而言,如图8所示,合成处理步骤S800包括:步骤S830,在偏转图像上检测并分割出对应介入性物体的图像;步骤S850,将分割出的对应图像与正常图像进行叠加,得到新的合成图像。这里先从偏转图像分割出介入性物体图像再与正常图像进行叠加,可以避免偏转图像中的组织图像影响正常图像中的组织图像。具体的叠加方法可参看后文。
本实施例中,介入性物体为穿刺针,穿刺针在偏转图像中呈现为具有一定宽度的高亮的直线段或轻微的曲线段。以下以直线段为例对合成处理步骤S800进行详细描述。
步骤S830中,在偏转图像上检测并分割出对应介入性物体的图像,对于本实施例而言,是要在偏转图像上检测并分割出直线段。检测方法有如下几种:
检测方法一:采用现有图像处理技术中的Hough变换。
Hough变换是一种应用广泛的直线检测方法,基本原理如下:
在图像上建立一个直角坐标系,如图9(a)所示,原点在图像的左上角,x轴正方向向左,y轴正方向向下。过图像上任一个点(x,y)在该坐标系下的直线方程可写为:
xcosθ+ysinθ=ρ(1)
θ和ρ与直线的关系如图9(a)所示,从原点向直线做垂线,θ描述为x轴正方向顺时针旋转至垂线所成的角,-90°≤θ≤90°,ρ的大小是原点到直线的距离。令图像x方向的长度为Iwidth,y方向的长度为Iheight,则
过一个点有无数条直线,即对应无数个θ和ρ,从而可以构造一个新的参数空间,称为θ-ρ空间,或称为Hough空间。可以理解,图像上的一直线经过式(1)变换后对应于θ-ρ空间上的一点。如图9(b)所示。
用Hough变换检测直线的算法如下:
(1)如图9(b)所示,将θ和ρ通过采样,可以将θ-ρ空间量化成许多小格,构成Hough矩阵;
(2)通过阈值处理将灰度图像变换成二值图像,根据上式(1)将二值图像上的值为1的各像素点(x,y)变换到θ-ρ空间,如在θ=-90°~90°内以小格的步长计算各个ρ值,所得值落入某个小格内,便使该小格的累加计数器加1;
(3)当图像中值为1的全部的点都变换后,对小格进行检验,计数值最大的小格,其(θ,ρ)值对应与图像中所求的直线。
检测方法二:采用本实施例提出的带权重的Hough变换方法进行穿刺针的检测,其流程如图10(a)所示,针对图像中的每一像素点:
步骤S1001′,根据偏转角度(偏转角度的获取方法同实施例一)确定θ-ρ空间中的θ值的取值范围,根据θ计算相应的ρ值,然后在θ-ρ空间中找到对应的θ和ρ位置;
步骤S1002′,计算权重。在该位置,根据下式(2)对加权参考值进行计算,将所述加权参考值映射为预定范围之间的权重,这里,加权参考值是指可以用于作为权重的一个参数,包括像素点的灰度值、梯度大小、梯度方向、偏转角度、进针位置中的任一种。本实施例中将图像的灰度根据权重映射函数(S-curve)映射成为0~1之间的权重Wgray,权重代表了该位置是穿刺针的可能性的大小。
本实施例中,权重映射函数如下:
其中,x为加权参考值,x∈[xmin,xmax],使用前可以先归一化到[0,1]区间;yhigh、ylow分别为预设定的映射值y的最大值和最小值;k是阶数;x0∈[0,1],称为breakpoint(指所映射的曲线中最陡峭的点);k和x0均用于控制所映射的曲线的形状。
其他实施例中,权重映射函数还可以如下:
式中各参数的定义与前式中的参数定义相同。
步骤S1003′,将该权重加到θ-ρ空间对应的位置上。
步骤S1004′,当图像中所有像素点都完成上述步骤S1001′-步骤S1003′后,可以得到Hough矩阵的最大值,该值对应的θ、ρ所对应的直线则是具有最可能是穿刺针所在的直线。
检测方法三:是对检测方法二的改进,首先对图像进行平滑,平滑的方法包括但不限于低通滤波、FIR、IIR、均值、中值、基于PDE的图像平滑等现有图像平滑方法。然后再对图像的每一像素点采用方法二进行穿刺针对应的直线段的检测。
检测方法四:是对检测方法二或检测方法三的改进,以对检测方法三改进为例,如图10(b)所示,在针对图像的每一像素点采用检测方法二时:
步骤S1001,对图像进行平滑,平滑方法采用现有图像平滑方法,在此不再赘述;
步骤S1002与检测方法二的步骤S1001′相同;
步骤S1003,与前述步骤S1002′类似,只是,加权参考值采用灰度、梯度大小、梯度方向(或者结构方向),进针方向(即偏转角度),进针位置等,采用式(2)分别进行映射,得到各权重Wgray、Wangle、Wmagnitude、Winject_angle、Winject_position;
步骤S1004,将步骤S1003中所有加权参考值各自映射后得到的所有权重相乘,得到一个总权重。将该权重加到θ-ρ空间对应的位置上。
步骤S1005,当图像中所有像素点完成上述步骤S1001-步骤S1004后,可以得到Hough矩阵的最大值,该值对应的θ、ρ所对应的直线则是具有最可能是穿刺针所在的直线。
检测方法五:本实施例中还可以基于穿刺针形状分割,例如采用基于直线或曲线模型的拟合方法(line-fitting)、radon变换等。具体算法可参考现有图像处理技术中的图像分割算法。
检测方法六:本实施例中还可以采用基于特征的阈值分割方法。算法如下:
1)在图像的特征域上进行阈值分割,变成二值图像,图像的特征可以是图像的灰度或梯度等;
2)对二值图像进行形态学操作,例如填充小空洞、去除小噪声等;
3)根据连通域的形状、灰度、位置等信息,选取最符合穿刺针的特征的连通域,即为分割后的穿刺针。
本方法中所涉及的阈值分割、形态学等算法可参考现有图像处理技术,在此不作描述。
可以理解,在偏转图像中,由于探头是垂直或基本垂直于穿刺针发射的,此时超声对穿刺针的显示是最强的,如图11(a)所示,穿刺针1102在偏转图像1101中特别清晰,图示为实线;而在正常图像中,由于探头是沿垂直于探头表面的方向发射,穿刺针1102′在正常图像1101′中会变弱,如图11(b)所示,图示为虚线。
回到步骤S830,在检测到介入性物体(本实施例中为穿刺针)后,分割出的对应介入性物体的图像(以下称为分割偏转图像),该分割偏转图像包括:介入性物体形成的介入图像(本实施例中穿刺针所形成的直线段)、和/或介入性物体形成的介入图像以及包含该介入图像的在预定范围内的图像。其中,介入性物体形成的介入图像以及包含该介入图像的在预定范围内的图像具体而言是指,利用步骤S830中的检测结果,将穿刺针所在直线标定在偏转图像相应的位置上,在该位置附近确定预定范围的像素区域,如图11(c)所示的区域1103,该预定范围的大小可以根据经验值或者不同所需取得。
回到图8,步骤S850中,将分割偏转图像与正常图像进行叠加,得到新的合成图像。应理解,分割偏转图像与正常图像相叠加的区域即为合成区域,即图11(c)所示的区域1103,其他区域为非合成区域。可以采用以下几种方法之一进行叠加(以下仍以介入性物体为穿刺针为例进行说明)。
叠加方法一:
可以将在分割偏转图像中直线段直接叠加在正常图像上的对应位置,即将新的灰度值(也就是偏转图像中直线段的灰度值)代替正常图像在该位置的原值。此外,为了予以区分,可以在该直线段上覆盖一层透明的伪彩。
叠加方法二:
本方法将分割偏转图像和正常图像上与该分割偏转图像对应的区域进行加权处理,将加权后的结果代替正常图像在该区域的原值。具体而言,针对正常图像中的每一个像素点以及与该像素点一一对应的在分割偏转图像中的每个像素点,如图12所示,包括如下步骤:
步骤S1201,判断正常图像的像素点的灰度值是否大于或等于其对应在偏转图像中的灰度值;如果是,则转步骤S1202,如果否,则转步骤S1203;
步骤S1202,将该正常图像的像素点的灰度值作为在新的合成图像中对应像素处的灰度值;
步骤S1203,利用前述权重映射函数(即式(2))对合成参考值进行计算以得到该合成参考值的合成权重,这里,合成参考值是指可以用于作为权重的一个参数,包括以下任一种:像素点的灰度值、Hough矩阵中的值、分割偏转图像中的像素点到介入图像的距离(本实施例中即指合成区域内的像素点到穿刺针所形成的直线段的距离)、分割偏转图像投影后的灰度最大值,具体权重计算如下:
A、对于合成参考值为像素点的灰度值,使用前述式(2)的权重映射函数来映射灰度值以得到权重值,灰度越大,权重越大,将该权重值作为合成权重;
B、对于合成参考值为Hough矩阵中的值时,根据前述步骤S830得到的Hough矩阵,使用前述式(2)的权重映射函数来映射Hough值,如果Hough值的权重小于预设阈值(根据经验值或统计数据设定),则使用正常图像的像素点的像素值作为在新的合成图像中对应像素处的像素值;否则,将该Hough值的权重值作为合成权重;
C、对于合成参考值为分割偏转图像中的像素点到介入图像时,在分割偏转图像中的每一像素点,计算其到检测出的直线段的距离,同样利用式(2)计算出该距离的权重,一般地,离直线段越近的点,其权重值越大,离直线段越远的点,其权重值越小,这样,在合成和非合成区域就会有平缓的过渡;
D、对于合成参考值为分割偏转图像投影后的灰度最大值时,如图12所示,在合成区域内,将合成区域沿竖直方向(如图13所示的竖直虚线方向)投影,此处为多列竖直方向,投影到每列竖直方向后得到灰度最大值,将该灰度最大值利用前述公式(2)映射为权重,灰度越大,权重越大;或者,另一种可选的方案是利用竖直方向上投影的连续性,判断穿刺针的针尖位置,对针体的位置赋予较大的权重,非针体的位置赋予很小的权重。
步骤S1204,根据合成权重对分割出的对应图像(即分割偏转图像)和正常图像进行合成,合成公式为:
I(x,y)=Isteer(x,y)*weight(x,y)+Inormal(x,y)*(1-weight(x,y))(3)
式中,I(x,y)为新的合成图像在像素点(x,y)的像素值,Isteer(x,y)是所述分割出的对应图像在像素点(x,y)的像素值,Inormal(x,y)是正常图像在像素点(x,y)的像素值,weight是所述合成权重。
可以理解,分割偏转图像为介入性物体形成的介入图像时,该叠加方法二需要做稍微的修改。
叠加方法三:
该方法是对叠加方法二的改进,即针对正常图像中的每一个像素点以及与该像素点一一对应的在偏转图像中的每个像素点,同样执行前述步骤S1201、步骤S1202,但在步骤S1203中,合成参考值包括:像素点的灰度值、Hough矩阵中的值、分割偏转图像中的像素点到介入图像的距离、分割偏转图像投影后的灰度最大值,各合成参考值的权重计算方法与前述方法二相同,然后将所有合成参考值映射后得到的各个权重相乘得到一个总权重,该总权重即为合成权重;然后接着执行步骤S1204。
同样地,可以理解,分割偏转图像为介入性物体形成的介入图像时,该叠加方法三需要做稍微的修改。
回到图6(b),显示步骤S629中,对合成处理步骤S627后得到的新的合成图像进行显示。此外,在显示时,可以在合成图像上增加一个框,该框可以是梯形、三角形等,该框的主要目的是在于指示穿刺针的有效位置,穿刺针如果越过该框将得不到前述的加权处理。如图14所示,所增加的框为梯形,图中1402为图像上显示的穿刺针,1403为穿刺目标(即待检测对象)。该框的显示方案还可以用其他形式替代,如覆盖一层透明的伪彩,显示一个指示穿刺针有效位置的棒形的框,和/或显示一个小框使穿刺目标位于该小框内,等。
实际显示时,可以在屏幕上双图(即正常图像和偏转图像)双实时显示得到的合成图像和未进行合成处理过的正常图像,便于用户对比和发现增强图上的误增强部分。用户也可以控制界面只显示合成图像。本领域技术人员应理解,在显示前可以采用DSC分别对两幅图像处理后再显示。
实施例三:
本实施例中第一角度只有一个,第二角度为多个,也就是说,探头从一个第一角度和多个第二角度进行发射,即此时空间复合开启,如图7(b)所示,其中第一角度ψ1为所示的虚线,多个第二角度ψ21、…ψ2n(n为自然数)为所示的实线。
偏转角度的获取方法同实施例一,在此不再重述。
本实施例仍采用图6(a)所示流程,在发射步骤S611中,确定了偏转角度后,偏转发射和正常发射交替进行。实施例中偏转发射位于所有正常发射之前或之后,也就是说,可以先进行所有的正常发射再进行偏转发射,或者反过来,即先进行偏转发射再进行所有的正常发射。可以理解,发射方式为先进行空间复合成像的发射然后再进行偏转发射(即从垂直或基本垂直于介入性物体的方向发射),或者为先偏转发射再进行空间复合成像的发射。或者,可以在偏转发射中做标记,这样可以将属于空间复合(即针对待检测对象的发射)与属于偏转发射(即针对介入性物体的发射)区分开。
接收步骤S613类似于实施例一中的接收步骤S613,在此不再重述。
信号处理步骤S615中,由于采用空间复合开启,因此,在对从多个第二角度的发射、接收及相应的信号处理后,信号处理步骤S615后还有一个空间复合步骤(图未示出),用于对多个第二角度回波信号进行信号处理后的数据进行复合以形成一帧第二角度图像,具体复合方法可以采用现有的空间复合成像方法,在此不再重述。
显示步骤S619类似于实施例一中显示步骤S619,在此不再重述。
实施例四:
本实施例中第一角度只有一个,第二角度仍为多个,类似于实施例三,即此时空间复合仍然开启。
本实施例仍采用图6(b)所示流程,其中,发射步骤S621、接收步骤S623、信号处理步骤S625分别类似于实施例三中图6(a)的发射步骤S611、接收步骤S613、信号处理步骤S615,在此不再重述。
在合成处理步骤S627中,将正常图像和偏转图像进行合成处理以生成新的合成图像,在本实施例中是指,将所有第二角度图像进行复合以生成复合的第二角度图像,然后将复合的第二角度图像与第一角度图像进行合成以生成新的合成图像。如图15所示,合成处理步骤S1500包括:步骤S1510,复合所有正常图像;步骤S1530,在偏转图像上检测并分割出对应介入性物体的图像;步骤S1550,将分割出的对应图像与复合后的正常图像进行叠加,得到新的合成图像。
本实施例仍以介入性物体为穿刺针为例对合成处理步骤S1500进行描述。步骤S1510中采用现有的空间复合成像方法将本次发射中的所有第二角度图像进行复合以形成一帧第二角度图像(即空间复合后的正常图像,简称为正常图像);步骤S1530和S1550分别与实施例二中的步骤S830和S750相类似,在此不再重述。
显示步骤S629与实施例一中的显示步骤S619相似,在此不再重述。
实施例五:
本实施例采用如图6(c)所示的流程,流程中各步骤所采用的处理过程可参考前述实施例一至实施例四中的相类似步骤的处理过程,在此不再重述。
实施例六:
本实施例采用如图6(d)所示的流程,流程中所采用的处理过程可参考前述实施例一至实施例四中的处理过程,在此不再重述。
实施例七:
本实施例是上述各实施例进行的改进。在获得第一角度后,以获得的第一角度进行发射,并在第一角度的基础上微调一个或多个角度后发射,这些角度对应的图像都视为第一角度图像,然后再从这些第一角度图像中择优选择或者采用类似空间复合的方法将这些第一角度图像进行叠加以复合成最终的一帧第一角度图像。而第二角度图像的获取方式可以采用前述六个实施例中任一方法得到。最后,将最终的第一角度图像和最终的第二角度图像进行合成处理,具体处理方法参见前述各实施例中的描述,在此不再重述。
以上各个实施例中,通过在正常发射的同时,增加一次垂直或基本垂直于穿刺针的偏转发射,偏转角度可由用户选择;接收到的偏转发射信号经过信号处理后,形成一帧偏转图像,将偏转图像和正常图像进行合成,即可得到对穿刺针的局部合成处理的图像。如前述,进行偏转发射的目的是使对介入性物体的显示较强,而正常发射可以得到待检测对象的正常图像,而不会因受偏转发射影响而使得图像的有效区域变小;这样,通过将偏转发射和正常发射交替进行,则可达到既增强显示介入性物体形成的图像,又保证了正常图像的有效区域不会变小,图像质量不会变差。
可以理解,本发明各实施方式及实施例的使用范围不仅限于穿刺针,还可以应用于:麻醉、适用于超声成像中对于扫描的某方向敏感的组织或介入性的物体。
上述实施例只是本发明的举例,尽管为说明目的公开了本发明的最佳实施例和附图,但是本领域的技术人员可以理解:在不脱离本发明及所附的权利要求的精神和范围内,各种替换、变化和修改都是可能的。因此,本发明不应局限于最佳实施例和附图所公开的内容。