CN113662638B - 一种用于超声穿刺引导的超声成像方法及装置 - Google Patents
一种用于超声穿刺引导的超声成像方法及装置 Download PDFInfo
- Publication number
- CN113662638B CN113662638B CN202110962138.XA CN202110962138A CN113662638B CN 113662638 B CN113662638 B CN 113662638B CN 202110962138 A CN202110962138 A CN 202110962138A CN 113662638 B CN113662638 B CN 113662638B
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- scanning
- image data
- deflection angle
- optimal
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title abstract description 31
- 230000005540 biological transmission Effects 0.000 claims abstract description 48
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000002604 ultrasonography Methods 0.000 claims abstract description 20
- 238000013329 compounding Methods 0.000 claims abstract description 18
- 230000006835 compression Effects 0.000 claims abstract description 15
- 238000007906 compression Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims description 33
- 239000000523 sample Substances 0.000 claims description 19
- 230000005570 vertical transmission Effects 0.000 claims description 12
- 238000012285 ultrasound imaging Methods 0.000 claims description 10
- 238000013507 mapping Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- 230000001351 cycling effect Effects 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- 230000006872 improvement Effects 0.000 abstract description 4
- 238000004422 calculation algorithm Methods 0.000 description 6
- 150000001875 compounds Chemical class 0.000 description 4
- 230000002238 attenuated effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000002787 reinforcement Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
- A61B17/34—Trocars; Puncturing needles
- A61B17/3403—Needle locating or guiding means
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
- A61B8/0841—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating instruments
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
- A61B17/34—Trocars; Puncturing needles
- A61B17/3403—Needle locating or guiding means
- A61B2017/3413—Needle locating or guiding means guided by ultrasound
Abstract
本发明公开一种用于超声穿刺引导的超声成像方法及装置,具体步骤:预设扫描表和延时序列;按照扫描表和延时序列完成超声的发射和接收,生成最优偏转角度范围扫描表和对应的延时序列;按照扫描表和最优偏转角度范围对应的延时序列完成超声的发射和接收,生成最优偏转角度扫描表和对应的延时序列;按照扫描表和最优偏转角度对应延时序列完成超声的发射和接收,获取偏转原始图像数据和垂直原始图像数据,并对所获取的原始图像数据进行空间复合,经对数压缩后得到空间复合后的显示图像,实现穿刺针的增强显示。本发明在不改变传统超声穿刺引导装置方法的基础上,仅修正扫描表和延时序列并更改空间复合策略,来实现成像质量的提升和计算效率的提升。
Description
技术领域
本发明涉及到超声成像与临床应用交叉技术领域,具体涉及一种用于超声穿刺引导的超声成像方法及装置。
背景技术
超声引导穿刺在临床检查、治疗等领域应用广泛,而在超声引导穿刺中超声成像质量对于穿刺引导具有非同寻常的重要意义。但是超声波会在人体内反射,折射和衰减,所以探头发射的超声波只有一部分能够反射回超声探头并经过波束合成形成B型超声图像。而超声穿刺引导所用的B型超声图像特别要求在图像中穿刺针能够突出显示,但如果反射回来的声波太少,强度太弱,在超声图像对比度降低,图像中的穿刺针和周围组织难以区分,从而则势必造成超声穿刺过程中对穿刺针位置的判断失误,并进而对穿刺的准确性造成影响。
针对超声声波衰减造成的组织成像所需信号随传输距离逐渐衰减这一固有问题,时间增益补偿(TGC)被引入到了超声成像数据处理当中,即在数据处理时,根据声学衰减特性,对远端回波信号进行预设倍率的放大补偿。但在超声引导穿刺这一特定应用中,穿刺针往往贯穿整个扫描深度且具有以反射为主的、不同于周围组织的声学特性,另外,穿刺针的角度随穿刺部位、医生手法经验等多重因素影响很难做到标准、可重复的穿刺针位置放置。故利用常规的TGC方法很难获得穿刺针灰度均一,且灰度高于周围组织的理想成像状态。
传统的超声探头更多的是实现超声信号的垂直发射和垂直接收,但对于某些与超声阵元排列方向倾斜的声波强反射物体,采用常规成像方法则有部分信号丢失,即无法接收到部分组织反射信号的风险,尽管超声波的偏转发射是可行的,但是,偏转发射获得的图像实际上的成像范围在探头方向偏转了一定角度,并不符合用户通常的观察习惯,需要对垂直发射和偏转发射得到的图像进行合成,而在穿刺流程中,穿刺针与人体表面所成的角度受多方面因素影响,难以量化为一个有效的固定角度甚至是一个固定角度范围。也就是说,超声穿刺因缺乏对偏转目标的角度信息,难以采用经典的空间复合方法来对偏转方向上的穿刺针增强显示。
图1给出了一种用于穿刺引导的空间复合超声成像探头超声发射方式,即在传统空间复合方法基础上,使用多个偏转角度完成空间复合。从理论上来讲,这种方法是可行的。但是,该方法由于需要对扫描面进行多次超声发射接收、再运用多帧复合算法合成超声图像,则将大大降低超声成像帧率,且超声成像帧率随采用的感兴趣角度个数增加而降低。考虑到超声穿刺引导对成像帧率的要求,有部分厂商仅选取2、3个可能的角度进行偏转发射,并优化算法计算流程,采用专有计算硬件进行数据处理,能够勉强获得满足超声穿刺扫描应用需求的超声产品。但是如果穿刺角度与机器内置的几个偏转发射角度相差较大,这样的方法获取的超声穿刺图像仍有无法看清穿刺针的风险。
在超声穿刺引导领域,便携式超声往往更受临床医师青睐:一方面是便携式超声设备价格低廉,易于临床科室开展使用;另一方面是相较于大型设备,便携式超声设备小巧易用,功能专一性强。但便携式超声的本身特性决定了其所携带的计算硬件规模无法与大型设备相比。因此,在便携式超声设备中若要保证超声穿刺引导的图像清晰还需要设计在理论上计算量更小的优化图像处理方案,同时该方案也要解决现有技术中部分特殊情况下穿刺针不可见的问题。
发明内容
本发明针对现有技术或增加计算代价或降低成像帧率,且存在图像质量降低风险的固有缺陷,提出一种用于超声穿刺引导的超声成像方法及装置,在保证成像质量的同时,降低成像计算代价。
为实现上述技术目的,本发明采用的技术方案是:
一种用于超声穿刺引导的超声成像方法,具体步骤为:
步骤1:预设扫描表和延时序列,所述扫描表的左侧部分用于形成可视化图像用于显示,右侧部分的扫描表用于实现多个固定位置的多角度扫描线发射;
步骤2:超声设备按照预设的扫描表和延时序列完成超声的发射和接收,获取第一帧显示图像和最优的偏转发射角度范围,生成最优偏转角度范围扫描表和对应的延时序列;
步骤3:超声设备按照预设的扫描表和最优偏转角度范围对应的延时序列控制通道完成超声的发射和接收,获取第二帧显示图像和最优偏转角度,生成最优偏转角度扫描表和对应的延时序列;
步骤4:超声设备按照扫描表的左侧部分和最优偏转角度对应延时序列控制通道完成超声的发射和接收,获取偏转原始图像数据和垂直原始图像数据,并对所获取的原始图像数据进行空间复合,并经对数压缩后得到空间复合后的显示图像,实现穿刺针的增强显示。
进一步地,在步骤2中,具体步骤为:
步骤2.1:按照扫描表的左侧部分对扫描线实施正常垂直发射、不加延时进行超声发射/接收,再按照扫描表的右侧部分的固定位置阵元扫描线根据预先设置的扫描延时序列进行超声偏转发射/接收;
步骤2.2:对正常发射接收的扫描线进行超声图像数据处理、对数压缩,形成第一帧显示图像;
步骤2.3:对偏转发射接收的扫描线进行超声图像数据处理,按照不同偏转角度分组,进行最优偏转角度范围分析;若无法找到最优偏转角度范围,则返回步骤2.1;否则生成最优偏转角度范围扫描表和对应的延时序列。
进一步地,在步骤2.3中,分析最优偏转角度范围的方法为:
步骤2.3.1:以多条不同位置但延时设置相同的扫描线为一组对扫描线原始数据进行分组;
步骤2.3.2:对每个扫描组中每一条扫描线像素求灰度最大值;
步骤2.3.3:按照扫描序列对应的实际空间位置将每一扫描组中灰度最大值对应坐标点映射到真实坐标系中,
映射公式为:
x=x0+nr·sinθ
y=nr·cosθ
其中x,y表示真实世界的横纵坐标,x0表示当前扫描线原点在真实世界中的位置,n表示最大值点在当前扫描线序列的位置,r表示为扫描线的纵向分辨率,θ为扫描偏转角度;
步骤2.3.4:根据不同的分组,分别检验该组中多个位置坐标是否共线,统计共线的数据分组;若所有分组中没有共线的位置坐标,则不存在最优角度范围,重复步骤2.1,否则将进入下一步骤;
步骤2.3.5:统计共线的角度分组,并计算共线的倾斜角角度;
步骤2.3.6:去除像素灰度较低、倾斜角过大过小的扫描线组;
步骤2.3.7:将各组倾斜角映射到固定范围内;
步骤2.3.8:多组数据进行最优倾斜角度投票,得到最优偏转角度范围。
进一步地,在步骤3中,具体步骤为:
步骤3.1:按照扫描表左侧部分对扫描线实施正常垂直发射、不加延时进行超声发射/接收,再按照扫描表的右侧部分的固定位置阵元扫描线在最优偏转角度范围附近根据步骤2所获得的对应延时进行超声偏转扫描发射/接收;
步骤3.2:对正常发射接收的扫描线进行超声图像数据处理、对数压缩,形成第二帧显示图像;
步骤3.3:对偏转发射接收的扫描线进行超声图像数据处理,按照不同偏转角度分组,进行最优偏转角度计算,生成最优偏转角度扫描表和对应的延时序列。
进一步地,在步骤3.3中,最优偏转角度的计算方法为:
3.3.1:将扫描线原始数据按角度分组;
3.3.2:获取每一分组中每一条扫描线的最大灰度像素值;
3.3.3:按组求最大灰度像素值的均值,取最大灰度均值序列中的极大值对应扫描角度,即为最优偏转角度。
进一步地,在步骤4中,具体步骤为:
步骤4.1:按照扫描表左侧部分和最优偏转角度对应的延时序列对扫描线实施超声发射/接收,获取偏转原始图像数据;
步骤4.2:按照扫描表左侧部分,对扫描线实施正常垂直发射、不加延时进行超声发射/接收,获取垂直原始图像数据;
步骤4.3:对步骤4.1和步骤4.2步骤获取的原始图像数据进行空间复合,实现穿刺针的增强显示;
步骤4.4:循环步骤4.1到步骤4.3,直至扫描结束。
进一步地,在步骤4.4中,若出现穿刺针显示不清晰,则返回到步骤2重复进行。
一种用于超声穿刺引导的超声成像装置,所述装置包括:
超声探头,用于穿刺时,获取垂直原始图像数据和偏转原始图像数据;
控制模块,用于控制超声探头发射/接收超声波,并控制超声探头各阵元通道关闭/开启、从而控制是否相对延时的超声发射接收策略;
数据处理模块,用于对所获取的偏转原始图像数据进行分析、计算,得到最优偏转角度;并对所述垂直原始图像数据和偏转原始图像数据进行空间复合;
显示模块,与数据处理模块通信连接,用于显示复合后的帧图像
与现有技术相比,本发明的有益效果有:
本发明的用于超声穿刺引导的超声成像方法及超声成像装置,在不改变传统超声穿刺引导装置方法的基础上,仅仅修正扫描表和延时序列并更改空间复合策略,即可实现成像质量的提升和计算效率的提升。
本发明的方法和装置针对便携式超声的应用具有良好的前景和价值,可以让超声穿刺引导应用场景中的超声设备针显影功能的强化和设备小型化成为可能。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本申请的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1是现有技术的一种用于穿刺引导的空间复合超声成像探头超声发射方式示意图;
图2是本发明的用于超声穿刺引导的超声成像方法的总体流程图;
图3是本发明的改进的超声探头发射和接收扫描表示意图;
图4是本发明中的扫描线设置示意图;
图5是本发明中的超声穿刺引导最优扫描角度范围分析算法流程图;
图6是本发明中的超声穿刺引导最优扫描角度计算算法流程图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围,而是仅仅表示本申请的选定实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
一种用于超声穿刺引导的超声成像方法,需要说明,本发明中称经超声数据处理后、对数压缩前的数据为原始图像数据,本发明中实现超声穿刺引导扫描在扫描开始之后先按照扫描表控制扫描获得两帧图像以及部分穿刺角度附加信息,然后再利用该两帧得到的穿刺针角度信息开展后续的穿刺成像,最后成像采用垂直阵元方向的发射扫描线和垂直于穿刺针的最优角度发射扫描线进行空间复合成像。图2给出本发明中用于超声穿刺引导的超声成像方法的总体流程图。下面分别对这该总体流程进行介绍。
当超声引导穿刺开始之后,超声设备开始工作。具体步骤为:
步骤1:预设扫描表和延时序列。
所述的扫描表如图3所示。扫描表在左侧部分是传统经典扫描表,扫描表横向坐标轴编号代表了成像扫描线的编号,而纵向坐标轴编号代表可控制的发射/接收阵元标号,扫描表中的0、1表示该发射阵元在完成某一条特定扫描线成像数据获取时是关闭还是开启,在该表中可以看出扫描线数目为M,可控发射/接收阵元的为N,系统同时可以控制的阵元数目(一般称为通道数)设为5。根据扫描表设置,并控制各个通道发射/接收不存在相对延时的超声发射接收策略,可实现超声信号的垂直发射和垂直接收,根据超声数据处理和对数压缩形成可视化图像用于显示;而右侧部分的扫描表功能是实现多个固定位置的多角度扫描线发射,如图4:其中用到的固定位置为位于阵元排列方向第0.25N、0.5N和0.75N的3个阵元(设阵元数目为N),扫描表相应的序列应该设置为该固定位置附近有效;而多个偏转角度则是不同扫描线成像所需的多个发射/接收通道实施差异性延时形成,假定偏转角度数目为k。因此,整个扫描表的宽度为M+3k。
步骤2:超声设备按照预设的扫描表和延时序列完成超声的发射和接收,获取第一帧显示图像和最优的偏转发射角度范围,生成最优偏转角度范围扫描表和对应的延时序列。
具体的,步骤2.1:按照扫描表的左侧部分对扫描线实施正常垂直发射、不加延时进行超声发射/接收,再按照扫描表的右侧部分的固定位置阵元扫描线根据预先设置的扫描延时序列进行超声偏转发射/接收。
步骤2.2:对正常发射接收的扫描线进行超声图像数据处理、对数压缩,形成超声图像。
步骤2.3:对偏转发射接收的扫描线进行超声图像数据处理,按照不同偏转角度分组,进行最优偏转角度范围分析。若无法找到最优偏转角度范围,则返回步骤2.1;否则生成最优偏转角度范围扫描表和对应的延时序列。
具体的,图5给出了超声穿刺引导最优偏转角度范围分析算法流程图,按照扫描表中发射、接收之后获取并进行了基本数据处理后的扫描线数据,以多条不同位置但延时设置相同的扫描线为一组,可分为k个扫描组。由于金属穿刺针对于超声反射系数较大,在某些偏转角度下获取的超声扫描线会呈现某点的高灰度值,因此,接下来可对k个扫描组中的每一条扫描线求取灰度最大值,并可按照扫描序列对应的实际空间位置将这些最大值坐标点映射到真实世界坐标系当中,映射公式为:
x=x0+nr·sinθ
y=nr·cosθ
其中x,y表示真实世界的横纵坐标,x0表示当前扫描线原点在真实世界中的位置,n表示最大值点在当前扫描线序列的位置,r表示为扫描线的纵向分辨率,θ为扫描偏转角度。
当完成坐标映射之后,根据不同的分组,分别检验该组中多个位置坐标是否共线。若所有这些分组中都没有共线的位置坐标,则重新按照该扫描表发射超声成像;否则,统计这些共线的角度分组,并计算共线的倾斜角角度。
考虑到图像中噪声干扰,还需去除一些最大值低于阈值的角度组;而对于人体中有些能够呈现高反射、且表面平整的组织,相比于穿刺针其走向更趋于平行于探头表面,根据这一特征可以通过设置角度阈值将其去除。
至此,得到某些候选的角度分组,每一个角度分组对应一个计算出的共线倾斜角。
根据计算出的共线倾斜角,统计其落入的角度区间,最终根据有效落入不同区间的角度分组数目投票找到最优的偏转角度范围。
若由于阈值处理导致此时无法计算出合理的穿刺角度范围,则重新按照该扫描表发射超声成像。
步骤3:超声设备按照扫描表和最优偏转角度范围对应的延时序列控制通道完成超声的发射和接收,获取第二帧显示图像和最优偏转角度,生成最优偏转角度扫描表和对应的延时序列。
按照步骤2所生成的最优偏转角度范围扫描表和对应延时序列,进行下一帧回波数据采集。类似第一帧数据采集,扫描表左侧数据用于形成帧图像用于显示,右侧数据用于计算最优扫描角度。图6给出了超声穿刺引导最优扫描角度计算算法流程图:扫描线数据首先按角度分组,然后获取每一组中每一条扫描线的最大灰度像素值,并按组求最大灰度像素值的均值,最后找到均值最大组对应的超声扫描线发射角度,这个角度就认为是最优偏转角度。
具体的:
步骤3.1:按照扫描表左侧部分对扫描线实施正常垂直发射、不加延时进行超声发射/接收,再按照扫描表的右侧部分的固定位置阵元扫描线在最优偏转角度范围附近根据步骤2所获得的对应延时进行超声偏转扫描发射/接收。
步骤3.2:对正常发射接收的扫描线进行超声图像数据处理、对数压缩,形成超声图像。
步骤3.3:对偏转发射接收的扫描线进行超声图像数据处理,按照不同偏转角度分组,进行最优偏转角度计算,生成最优偏转角度扫描表和对应的延时序列。
步骤4:超声设备按照扫描表左侧部分和最优偏转角度对应延时序列控制通道完成超声的发射和接收,获取偏转原始图像数据和垂直原始图像数据,并对所获取的原始图像数据进行空间复合,实现穿刺针的增强显示。
步骤4.1:按照扫描表左侧部分和最优偏转角度对应的延时序列对扫描线实施超声发射/接收,获取偏转原始图像数据;
步骤4.2:按照扫描表左侧部分,对扫描线实施正常垂直发射、不加延时进行超声发射/接收,获取垂直原始图像数据;
步骤4.3:对步骤4.1和步骤4.2步骤获取的原始图像数据进行空间复合,实现穿刺针的增强显示;
步骤4.4:循环步骤4.1到步骤4.3,直至扫描结束。若在其中若有穿刺针显示不清晰,则返回到步骤2进行。
为进一步详细说明本发明的具体内容和步骤,以下以假设第一帧、第二帧设置的扫描角度个数k为10,阵元N数目为256,扫描线M数目为256为例具体说明。本发明的步骤如下:
S1,预设扫描表和延时序列,针对不同的发射角度,以10度为分度计算各个角度所需的偏转延时序列;根据图3的扫描表,设置(256+3*10)*256大小的第一帧、第二帧扫描表以及256*256大小的后续帧扫描表。
S2.1,按照(256+3*10)*256扫描表,首先对左侧部分256个扫描线实施正常垂直发射、不加延时进行超声发射/接收,再对右侧部分分列3个阵元位置的30个扫描线以10度为分度,从0度到90度根据预先设置的扫描延时进行超声扫描发射/接收。
S2.2,对正常发射接收的扫描线进行超声图像数据处理、对数压缩,形成超声图像。
S2.3,对偏转发射接收的扫描线进行超声图像数据处理,按照不同偏转角度三三分组,以30度为一个范围区间进行最优偏转角度范围分析。若无法找到最优偏转角度范围,则返回S2.1;否则生成最优偏转角度范围扫描表和对应的延时序列,继续进行下述步骤。
S3.1,按照(256+3*10)*256扫描表,首先对左侧部分256个扫描线实施正常垂直发射、不加延时进行超声发射/接收,再对右侧部分分列3个阵元位置的10组共计30个扫描线以3度为分度,在最优偏转角度范围附近根据S2.3生成最优偏转角度范围对应的延时序列进行超声扫描发射/接收。
S3.2,对正常发射接收的扫描线进行超声图像数据处理、对数压缩,形成超声图像。
S3.3,对偏转发射接收的扫描线进行超声图像数据处理,按照不同偏转角度三三分组,进行最优偏转角度计算,生成最优偏转角度扫描表和对应的延时序列。
S4.1,以256*256扫描表和最优偏转角度对应的延时序列,对扫描线实施超声发射/接收,实现偏转原始图像数据获取。
S4.2,以256*256扫描表,对扫描线实施正常垂直发射、不加延时进行超声发射/接收,实现垂直原始图像数据获取。
S4.3,对S4.1和S4.2步骤获取的原始图像数据进行空间复合,对数压缩,得到空间复合后的显示图像,实现穿刺针的增强显示。
S4.4,循环S4.1到S4.3步骤,直至扫描结束。若在其中若有穿刺针显示不清晰,则返回到S2步骤进行。
本发明还公开一种用于超声穿刺引导的超声成像装置,所述装置包括:
超声探头,用于穿刺时,获取垂直原始图像数据和偏转原始图像数据;
控制模块,用于控制超声探头发射/接收超声波,并控制超声探头各阵元通道关闭/开启、从而控制是否相对延时的超声发射接收策略,以得到垂直原始图像数据和多位置多角度扫描线偏转原始数据;
数据处理模块,用于对所获取的偏转原始图像数据进行分析、计算,得到最优偏转角度;并对所述垂直原始图像数据和偏转原始图像数据进行空间复合;
显示模块,与数据处理模块通信连接,用于显示复合后的帧图像。
在穿刺时,本装置的超声探头根据控制模块输出的信号向人体组织发射/接收超声波,超声探头的各阵元可关闭/开启,可设置相对延时,以得到垂直原始图像数据和多位置多角度的偏转原始图像数据,所述数据处理模块对所获得的多位置多角度的偏转原始图像数据进行最优偏转角度范围/最优偏转角度分析计算,生成最优偏转角度范围/最优偏转角度扫描表和对应的延时序列,所述控制模块根据数据处理模块的计算结果,实时更新对应的延时序列,以得到最优偏转角度的偏转原始图像数据,然后所述数据处理模块根据垂直原始图像数据和最优偏转角度的偏转原始图像数据进行空间复合、对数压缩,得到空间复合后的帧图像,实时显示于所述显示模块,引导穿刺。
本发明的方法和装置在不改变传统超声穿刺引导装置方法的基础上,仅仅修正扫描表和延时序列并更改空间复合策略,来实现成像质量的提升和计算效率的提升。本发明的方法和装置针对便携式超声的应用具有良好的前景和价值,可以让超声穿刺引导应用场景中的超声设备针显影功能的强化和设备小型化成为可能。
当然,本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (4)
1.一种用于超声穿刺引导的超声成像方法,其特征在于:具体步骤为:
步骤1:预设扫描表和延时序列,所述扫描表的左侧部分用于形成可视化图像用于显示,右侧部分的扫描表用于实现多个固定位置的多角度扫描线发射;
步骤2:超声设备按照预设的扫描表和延时序列完成超声的发射和接收,获取第一帧显示图像和最优的偏转发射角度范围,生成最优偏转角度范围扫描表和对应的延时序列;具体步骤为:
步骤2.1:按照扫描表的左侧部分对扫描线实施正常垂直发射、不加延时进行超声发射/接收,再按照扫描表的右侧部分的固定位置阵元扫描线根据预先设置的扫描延时序列进行超声偏转发射/接收;
步骤2.2:对正常发射接收的扫描线进行超声图像数据处理、对数压缩,形成第一帧显示图像;
步骤2.3:对偏转发射接收的扫描线进行超声图像数据处理,按照不同偏转角度分组,进行最优偏转角度范围分析;若无法找到最优偏转角度范围,则返回步骤2.1;否则生成最优偏转角度范围扫描表和对应的延时序列;分析最优偏转角度范围的方法为:
步骤2.3.1:以多条不同位置但延时设置相同的扫描线为一组对扫描线原始数据进行分组;
步骤2.3.2:对每个扫描组中每一条扫描线像素求灰度最大值;
步骤2.3.3:按照扫描序列对应的实际空间位置将每一扫描组中灰度最大值对应坐标点映射到真实坐标系中,
映射公式为:
x=x0+nr·sinθ
y=nr·cosθ
其中x,y表示真实世界的横纵坐标,x0表示当前扫描线原点在真实世界中的位置,n表示最大值点在当前扫描线序列的位置,r表示为扫描线的纵向分辨率,θ为扫描偏转角度;
步骤2.3.4:根据不同的分组,分别检验该组中多个位置坐标是否共线,统计共线的数据分组;若所有分组中没有共线的位置坐标,则不存在最优角度范围,重复步骤2.1,否则将进入下一步骤;
步骤2.3.5:统计共线的角度分组,并计算共线的倾斜角角度;
步骤2.3.6:去除像素灰度较低、倾斜角过大过小的扫描线组;
步骤2.3.7:将各组倾斜角映射到固定范围内;
步骤2.3.8:多组数据进行最优倾斜角度投票,得到最优偏转角度范围;
步骤3:超声设备按照预设的扫描表和最优偏转角度范围对应的延时序列控制通道完成超声的发射和接收,获取第二帧显示图像和最优偏转角度,生成最优偏转角度扫描表和对应的延时序列;具体步骤为:
步骤3.1:按照扫描表左侧部分对扫描线实施正常垂直发射、不加延时进行超声发射/接收,再按照扫描表的右侧部分的固定位置阵元扫描线在最优偏转角度范围附近根据步骤2所获得的对应延时进行超声偏转扫描发射/接收;
步骤3.2:对正常发射接收的扫描线进行超声图像数据处理、对数压缩,形成第二帧显示图像;
步骤3.3:对偏转发射接收的扫描线进行超声图像数据处理,按照不同偏转角度分组,进行最优偏转角度计算,生成最优偏转角度扫描表和对应的延时序列;最优偏转角度的计算方法为:
3.3.1:将扫描线原始数据按角度分组;
3.3.2:获取每一分组中每一条扫描线的最大灰度像素值;
3.3.3:按组求最大灰度像素值的均值,取最大灰度均值序列中的极大值对应扫描角度,即为最优偏转角度;
步骤4:超声设备按照扫描表的左侧部分和最优偏转角度对应延时序列控制通道完成超声的发射和接收,获取偏转原始图像数据和垂直原始图像数据,并对所获取的原始图像数据进行空间复合,并经对数压缩后得到空间复合后的显示图像。
2.根据权利要求1所述的用于超声穿刺引导的超声成像方法,其特征在于:在步骤4中,具体步骤为:
步骤4.1:按照扫描表左侧部分和最优偏转角度对应的延时序列对扫描线实施超声发射/接收,获取偏转原始图像数据;
步骤4.2:按照扫描表左侧部分,对扫描线实施正常垂直发射、不加延时进行超声发射/接收,获取垂直原始图像数据;
步骤4.3:对步骤4.1和步骤4.2步骤获取的原始图像数据进行空间复合,实现穿刺针的增强显示;
步骤4.4:循环步骤4.1到步骤4.3,直至扫描结束。
3.根据权利要求2所述的用于超声穿刺引导的超声成像方法,其特征在于:在步骤4.4中,若出现穿刺针显示不清晰,则返回到步骤2重复进行。
4.一种包含上述权利要求1-3任一所述的用于超声穿刺引导的超声成像方法的超声成像装置,其特征在于:所述装置包括:
超声探头,用于穿刺时,获取垂直原始图像数据和偏转原始图像数据;
控制模块,用于控制超声探头发射/接收超声波,并控制超声探头各阵元通道关闭/开启、及控制是否相对延时的超声发射接收策略;
数据处理模块,用于对所获取的偏转原始图像数据进行分析、计算,得到最优偏转角度;并对所述垂直原始图像数据和偏转原始图像数据进行空间复合;
显示模块,与数据处理模块通信连接,用于显示复合后的帧图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110962138.XA CN113662638B (zh) | 2021-08-20 | 2021-08-20 | 一种用于超声穿刺引导的超声成像方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110962138.XA CN113662638B (zh) | 2021-08-20 | 2021-08-20 | 一种用于超声穿刺引导的超声成像方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113662638A CN113662638A (zh) | 2021-11-19 |
CN113662638B true CN113662638B (zh) | 2024-04-19 |
Family
ID=78544789
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110962138.XA Active CN113662638B (zh) | 2021-08-20 | 2021-08-20 | 一种用于超声穿刺引导的超声成像方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113662638B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105496515A (zh) * | 2015-12-04 | 2016-04-20 | 深圳华声医疗技术有限公司 | 穿刺增强方法及系统 |
WO2018053739A1 (zh) * | 2016-09-20 | 2018-03-29 | 深圳华声医疗技术有限公司 | 穿刺增强方法、装置及系统 |
WO2018053776A1 (zh) * | 2016-09-21 | 2018-03-29 | 深圳华声医疗技术有限公司 | 超声图像标识的显示方法及装置 |
CN109953771A (zh) * | 2016-06-27 | 2019-07-02 | 中国科学院苏州生物医学工程技术研究所 | 超声成像方法、超声弹性成像方法及微型超声装置 |
CN110101411A (zh) * | 2019-05-28 | 2019-08-09 | 飞依诺科技(苏州)有限公司 | 超声成像空间复合方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010077632A2 (en) * | 2008-12-08 | 2010-07-08 | Silicon Valley Medical Instruments, Inc. | System and catheter for image guidance and methods thereof |
-
2021
- 2021-08-20 CN CN202110962138.XA patent/CN113662638B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105496515A (zh) * | 2015-12-04 | 2016-04-20 | 深圳华声医疗技术有限公司 | 穿刺增强方法及系统 |
CN109953771A (zh) * | 2016-06-27 | 2019-07-02 | 中国科学院苏州生物医学工程技术研究所 | 超声成像方法、超声弹性成像方法及微型超声装置 |
WO2018053739A1 (zh) * | 2016-09-20 | 2018-03-29 | 深圳华声医疗技术有限公司 | 穿刺增强方法、装置及系统 |
WO2018053776A1 (zh) * | 2016-09-21 | 2018-03-29 | 深圳华声医疗技术有限公司 | 超声图像标识的显示方法及装置 |
CN110101411A (zh) * | 2019-05-28 | 2019-08-09 | 飞依诺科技(苏州)有限公司 | 超声成像空间复合方法及系统 |
Non-Patent Citations (2)
Title |
---|
《Puncture procedures utilizing transvaginal ultrasonic guidance》;I.E Timor-Tritsch;D.B Peisner;A.Monteagudo;《Ultrasound in Obstetrics & Gynecology》;第1卷(第2期);全文 * |
《医用自动穿刺针偏转系统的设计与实现》;朱 勒;《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113662638A (zh) | 2021-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11471131B2 (en) | Ultrasound imaging system and method for displaying an acquisition quality level | |
US8834374B2 (en) | Setting an optimal image parameter in an ultrasound system | |
US11331076B2 (en) | Method and system for displaying ultrasonic elastic measurement | |
US20160192907A1 (en) | Clutter suppression in ultrasonic imaging systems | |
US20060058651A1 (en) | Method and apparatus for extending an ultrasound image field of view | |
US6572549B1 (en) | High frame rate extended field of view ultrasound imaging system and method | |
US20080262352A1 (en) | Clutter Suppression in Ultrasonic Imaging Systems | |
JP4903271B2 (ja) | 超音波撮像システム | |
US20160074017A1 (en) | Ultrasound optimization method and ultrasonic medical device therefor | |
US9622724B2 (en) | Ultrasound imaging system and method for tracking a specular reflector | |
CN105496515A (zh) | 穿刺增强方法及系统 | |
US8724880B2 (en) | Ultrasonic diagnostic apparatus and medical image processing apparatus | |
US20200138403A1 (en) | Ultrasound imaging system and method for detecting position and orientation of a coherent reflector | |
US9833218B2 (en) | Apparatus for processing ultrasonic image and method thereof | |
US20180000458A1 (en) | Ultrasonic imaging device and method for controlling same | |
CN109310397A (zh) | 超声成像设备、超声图像增强方法及引导穿刺显示方法 | |
CN110881997A (zh) | 超声波诊断装置以及体数据取入方法 | |
CN113662638B (zh) | 一种用于超声穿刺引导的超声成像方法及装置 | |
US20230320698A1 (en) | Ultrasonic diagnostic apparatus, method for controlling ultrasonic diagnostic apparatus, and control program for ultrasonic diagnostic apparatus | |
US20230061869A1 (en) | System and methods for beamforming sound speed selection | |
US11559280B2 (en) | Ultrasound imaging system and method for determining acoustic contact | |
CN111855824B (zh) | 超声装置及其控制方法 | |
LU500526B1 (en) | Image Optimization Method for Neonatal Lung Ultrasonography | |
CN111436971B (zh) | 一种超声混响伪像抑制方法 | |
US20240061109A1 (en) | Ultrasound diagnosis apparatus, image processing apparatus, and computer program product |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |