CN108024791B - 将肺滑动与外部运动区分 - Google Patents

将肺滑动与外部运动区分 Download PDF

Info

Publication number
CN108024791B
CN108024791B CN201680053990.2A CN201680053990A CN108024791B CN 108024791 B CN108024791 B CN 108024791B CN 201680053990 A CN201680053990 A CN 201680053990A CN 108024791 B CN108024791 B CN 108024791B
Authority
CN
China
Prior art keywords
sub
motion
image
imaging
lung
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
Application number
CN201680053990.2A
Other languages
English (en)
Other versions
CN108024791A (zh
Inventor
B·I·拉朱
徐泾平
王守罡
周诗未
A·M·加德斯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of CN108024791A publication Critical patent/CN108024791A/zh
Application granted granted Critical
Publication of CN108024791B publication Critical patent/CN108024791B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54Control of the diagnostic device
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • A61B8/085Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42Details of probe positioning or probe attachment to the patient
    • A61B8/4245Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient
    • A61B8/4254Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient using sensors mounted on the probe
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • A61B8/5276Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G01S15/892Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array the array being curvilinear
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52077Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference

Abstract

相对于被成像的医学对象的体外运动(130)通过成像或者根据成像探头上的运动检测器而被检测,并且从医学图像收回使得可以确定肺滑动是否存在,或者被测量以确定肺滑动检测是否要由于过度体外运动而暂停。对应于图像的相应子区域的图像子区域(164、168)被选择用于图像到图像比较,使得选定的子区域仅包含相对于对图像的采集中的成像深度比图像内的解剖标志更浅的身体组织。基于比较的结果,可以初始化涉及检查比标志更深的图像数据的肺滑动检测。运动传感器可以检测体外运动,并且基于其输出执行以下操作:将图像逐对共配准(170)到收回体外运动的效应的程度和/或确定是否暂停关于肺滑动是否在呼吸期间在所述对象中发生的决定。

Description

将肺滑动与外部运动区分
技术领域
本发明涉及使用成像来检测气胸,并且更具体地,涉及在存在外部运动的情况下的这样的检测。
背景技术
在世界范围内,胸部创伤占所有创伤的20%。早期诊断和适当的处置的及时选择是针对最佳结果的关键要素。医学成像在决定制定过程中扮演重要角色。在钝胸部创伤情况之中,气胸(PTX)表示在肋骨骨折之后第二最常见的损伤。PTX能够作为钝或穿透创伤的结果发生。假如患者不需要机械通气或者空气运输,小PTX通常无症状并且能够在没有胸腔引流的情况下安全地管理。在这些情况下,任何尺寸的PTX应当利用胸廓切开管处置,以避免张力PTX的创建。大的PTX能够导致呼吸窘迫,并且张力PTX能够导致心肺衰竭。在张力PTX中,空气在没有逃逸路线的情况下泄露到胸膜腔中,考虑到单向阀效应并且是威胁生命的条件。
尽管计算机断层摄影(CT)被认为是其体积的PTX检测和测量中的黄金标准,但是CT的主要问题在于,其不能在床旁在不稳定的患者上执行。在床旁完成的胸部X射线(具体而言前后仰卧位胸部放射线照相)能够丢失PTX情况的达到30%。一项研究发现,甚至当PTX体积为可调尺寸(350ml到500ml)时,用于诊断PTX的胸部X射线的灵敏度是低的。
超声成像(特别地定点照护超声(POC-US))是变得更常用于PTX的床旁检测的成像模态。然而,关于超声检查的主要问题是对于训练的需要和检查结果的高操作者相关性。贡献因子是新一类用户,其通常未特别地训练以成为超声医师。一项研究示出,与PTX检测中的超声成像的有经验的用户的90.9%的平均灵敏度相比较,无经验的用户仅能够实现57%的平均灵敏度。
在正常肺中,高频率线性阵列换能器或者弯曲阵列被定位在胸部的最前面点处以识别胸膜线。该线表现为近似地被定位为比肋骨的最浅范围深半厘米的回声水平线。胸膜线包括紧密相对的脏层胸膜和壁层胸膜。在正常肺中,脏层胸膜可以被看到抵着壁层胸膜来回地滑动,从而产生被称为“肺滑动”的生理现象,具有随着对象呼吸的闪耀或者闪烁外观。(一个或多个)B线伪影(还被称为“彗星尾”伪影)是垂直强回声线,其可以从相对胸膜在后部延伸到屏幕的底部或者最大成像深度。通常,上缘地并且下方地识别肋骨两者是可能的,并且从那里将胸膜线在水平方向上定位在两个肋骨之间。肺滑动和(一条或多条)垂直B线的存在排除PTX。
在具有大自由空气的PTX的情况下,空气积聚在壁层胸膜与脏层胸膜的层之间并且将两个层分开。在空气存在的情况下,空气的积聚防止超声波束传播通过肺。相应地,不存在肺滑动或(一条或多条)B线。胸膜线将仅包括壁层,其被看作单条固定线。
在具有中等或较小量的自由空气的PTX的情况下,存在部分正常肺,以及具有PTX的部分异常部分。在这种情况下,在两个部分之间可以存在肺点。此处,肺滑动可以横向地但非前面地注意,因为自由空气将优选地在胸部腔内上缘地并且前面地采集。可以通过查看若干肋间空间中的胸膜线来检测肺滑动的区和滑动的缺少之间的过渡的肺点。
Burlina等人的美国专利No.8914097(“Burlina专利”)(通过引用将其整体并入本文)通过计算在不同的时间处采集的图像的图像间运动的向量来检测肺滑动,每幅图像包括胸膜界面。向量形成跨每幅图像的全部的场。可以根据关于方向和幅度的向量中间的一致来检测探头平移。根据Blurlina专利,所识别的探头运动被防止影响PTX检测。Burlina专利提到运动观察可以出于更大的鲁棒性被组合,但是未以其他方式建议或者指定可以如何防止对PTX检测的不利影响。
发明内容
下面在本文中提出的事物涉及解决以上所讨论的问题中的一个或多个。
在某些方面中,本发明包括被配置为检测肺滑动的成像系统(100)。所述成像系统可以包括超声探头,所述超声探头具有透镜表面;以及处理器,其被配置为执行以下步骤:选择第一超声图像中的第一子区域(164)和第二超声图像中的第二子区域(182),所述第一子区域和所述第二子区域包括表示对象中比第一超声图像和第二超声图像内的解剖标志更靠近于所述透镜表面的组织的图像数据;通过比较第一子区域与第二子区域之间的相对对准来估计对象中的体外(extracorporeal)运动的效应;并且基于体外运动的效应的幅度来确定是否暂停对所述对象中的肺滑动的成像。所述系统可以包括被附接到所述探头的运动传感器。在一些方面中,所述的比较相对对准可以包括比较来自所述运动传感器的针对所述第一子区域和所述第二子区域的位置数据。
在一些方面中,体外运动涉及在采集所述第一超声图像和第二超声图像所跨的时间段(143)内所述透镜表面相对于所述对象的皮肤表面的部分(128)的运动。所述系统可以包括显示器(115),所述显示器被配置为显示所述第一子区域与第二子区域之间的相似度的空间图。所述相似度的空间图可以基于所述第一子区域与第二子区域之间的图像数据的互相关。所述处理器可以被配置为基于所述第一子区域与第二子区域之间的相对对准的第一运动阈值来确定是否暂停成像。所述处理器还可以被配置为发布指示是否暂停对肺滑动的成像的用户通知。所述用户通知可以例如包括二元指示符(S256)。在一些方面中,所述第一子区域和第二子区域(164、182)被预定义为比对应于远离所述透镜表面的距离的深度阈值(177)更浅。深度阈值(177)可以不大于1.8厘米。在某些方面中,所述解剖标志包括胸膜线(152)。所述空间图可以包括多个局部相关系数。在一些方面中,如果所述第一运动阈值被超过,则所述处理器被配置为暂停对肺滑动的成像。所述显示器还可以被配置为基于所述第一子区域与第二子区域之间的相对对准来显示具有指示肺滑动的可能性的长度的条。
在一些方面中,本发明包括成像扫描器被配置用于检测肺滑动,并且更具体地用于在呼吸期间的不同的时间处从对象采集图像。所述扫描器可以包括:
a)电路,其被配置用于为了检测相对于所述对象的体外运动,选择用于图像到图像比较的所述图像的相应那些的图像子区域,使得选定的子区域仅包含相对于所述采集中的成像深度比所述图像内的解剖标志更浅的身体组织;并且用于基于所述比较的结果针对所述肺滑动检测初始化比所述标志更深的图像数据的检查;和/或
b)检测所述体外运的运动传感器以及被配置用于基于所述传感器的输出执行以下操作的电路两者:
i)估计所述体外运动对一对图像之间的相对对准的效应;并且在收回所估计的效应的程度上,逐对共同配准以用于基于共配准对的两者构成之间的比较来决定肺滑动是否在呼吸期间的所述对象中发生;和/或
ii)确定是否暂停关于肺滑动是否在所述呼吸期间在所述对象中发生的决定。
下面在本文中还提出了对应的方法和对应的计算机可读介质。
在未按比例绘制的以下附图的帮助下,下面进一步阐述了具有肺滑动检测初始化能力的创造性扫描器的细节。
附图说明
图1是具有肺滑动检测初始化能力的示范性创造性成像系统的示意图,该示图还图示了本发明的背景下的相关概念;并且
图2是示出根据本发明的示范性操作的流程图。
具体实施方式
图1通过说明性和非限制性范例示出了具有肺滑动检测初始化能力的成像扫描器或成像系统100。通常,扫描器100也将具有完整的肺滑动检测能力。扫描器100包括图像分析器102和成像探头104。成像模态可以是超声,并且在下面在本文中所描述的实施例中被假定。然而,代替地可以使用其他医学成像模态。例如,救护车中的便携式磁共振成像(MRI)设备可以产生受由救护车运动引起的患者到线圈运动不利地影响的图像。
图像分析器102在数据和电源总线106上包括处理电路108、设备接口110、用户接口112和PTX检测器114。用户接口112包括显示器115和用户输入控制(未示出)。
处理电路108可以包括存储器116和处理器118。如在本文中进一步描述的,处理器可以被配置为执行以下步骤:选择第一超声图像中的第一子区域和第二超声图像中的第二子区域,第一子区域和第二子区域包括表示比第一超声图像和第二超声图像内的解剖标志更靠近于透镜表面的对象中的组织的图像数据;通过比较第一子区域与第二子区域之间的相对对准来估计对象中的体外运动的效应;并且基于体外运动的效应来确定是否暂停所述对象中的肺滑动的成像。
探头104可以包括一个或多个运动传感器120。后者被并入在探头104内或者以其他方式物理地附接到探头。结合Augustine等人(“Augustine”)的共同受让的美国专利No.5529070中的图5和图5a描述了并入在探头104内的实施方式,通过引用将其整体公开内容并入本文。
然而,根据关于体外运动检测或者校正在本文中提出的事物,代替于任何这样的传感器120或者除此之外,可以使用成像自身。
探头104可以抵着对象122(诸如人类或者动物患者)的皮肤表面手动地保持。探头104包括换能器124,诸如具有面126的一维或二维换能器阵列。跨在作为正交于方位角方向和海拔方向的方向的探头104的范围方向128上的面126执行成像。声凝胶层常常是被用于将面126与皮肤表面的部分128分离的介质。通过相对于皮肤表面部分128(或更一般地相对于对象122)的偶然的手运动,成像可以刚性地移动,如由图1中的双头箭头130所表示的。由救护车运动引起的手或患者移动可以创建成像中的刚性移动。探头104还可以在垂直于图纸的方向上移动,在图纸的平面中振动,围绕图纸的平面中的轴震动或者围绕平面内的正交轴扭转。探头104与图像分析器102之间或者图像分析器的部件108-114中的任何中间的通信中的任何可以被实施为有线或无线连接。
处理电路108、设备接口110、用户接口112和PTX检测器可利用以下项中所描述的肺滑动检测功能实施:a)Jinping Xu等人题的为“Device and method for automaticpneumothorax detection”的共同所有的PCT申请PCT/EP2015/071625(“Xu申请”),其中,所有发明人是具有2014年9月25日的有效递交日期的本专利申请中的发明人,通过引用将其整体内容并入本文;或者b)Burlina专利。此外,其他电路被包括在例如处理电路108和PTX检测器114中以用于实现肺滑动检测初始化能力。在本文中所提到的电路可以被描述为被配置为执行功能。配置被理解为例如要么通过硬连线要么借助于表示指令或者数据的所存储的电容性电荷。例如,实现用于滑动运动分析的功能性的指令的运行对处理电路进行配置。
扫描器100经由接收到的超声对对象122进行成像,并且所述超声在脉冲回波上背景下还由探头104发射。所采集的图形数据132在从探头104的范围(或者“轴向”)方向128上的相应的成像深度134处。成像在由吸气136和呼气138的周期组成的对象122的呼吸122期间发生。
自动化PTX检测的一些现有方法如下。
在Xu申请中,成像可以被采集为图像数据132的帧140、142。
在图1中的第一示范性帧140中,上肋骨144和下肋骨146针对正常肺被看到。在较大的成像深度134中,阴影148、150产生于肋骨144、146的超声吸收。作为解剖标志的胸膜线152在正常肺中包括紧密相对的壁层胸膜154和脏层胸膜156。
在PTX或不完整PTX的情况下,两个胸膜154、156针对对应地整个长度或者长度的部分通过纵向地被设置在肋骨144、146之间的空气远离彼此更远地分离。在后者情况中,肺点存在,即,其中,分离变得异常。因此,当PTX存在时,胸膜线152包括壁层胸膜。当不完整的PTX存在时,胸膜线152包括其中空气迫使两个胸膜154、156分开的仅壁层胸膜以及其中其紧密相对的两者胸膜。
在不同的随后的时间处采集第二帧142。
有问题地,影响成像的相对于对象122的体外运动130能够已经在跨被彼此比较的帧140、142的采集时间tj、tj+k的时间段143内发生,该比较旨在评估沿着胸膜线152的可能运动,其指示肺滑动的存在或者缺少。
Xu申请定义覆盖肋骨144、146之间的胸膜线152的部分的第一子区域155。定义第二子区域156,其仅包含其成像深度134比胸膜线152更浅的背景软组织。针对被比较的帧140、142,测量两个子区域155、156之间的相对运动。这基于以下关键洞察:当不存在相对于对象122的体外运动130(诸如偶然的手运动)时,不管可能肺滑动和其他内部运动,胸膜线152上方的区域将在成像中被看作固定的;然而,使成像相对于对象122位移的这样的运动将使得该区域表现为在成像中移动。因此,相对运动被评估以在存在的情况下检测肺滑动,从而降低体外运动130的不利影响。相对运动可以根据两个互相关157a、157b计算,一个使用第一子区域155并且另一个使用第二子区域156。应注意到,第一子区域155包含具有比胸膜线152更深的成像深度134的组织158。
如以上在本文中所提到的,Burlina专利形成在图像的全部被比较期间的运动向量161的场160。根据Blurlina专利,所识别的探头运动被防止影响PTX检测。Burlina专利提到运动观察可以出于更大的鲁棒性被组合,但是未以其他方式建议或者指定可以如何防止对PTX检测产生不利影响。此处同样应注意到,Burlina专利包括与具有比胸膜线152更深的成像深度134的组织162有关的运动向量。当信噪比(SNR)在胸膜线下弱时,运动向量估计对于Burlina技术而言是有问题的。
有利地根据在本文中提出的事物,关于PTX是否存在的自动化决定要么在体外运动130的检测时段期间推迟要么这样的运动对所采集的图像之间的相对对准的效应163被收回为决定过程的初始化。
具体地并且通过范例,选择与Xu申请的第二子区域156类似的帧166的“较早”子区域164。在一些实施例中,成像系统的处理器被配置为选择第一超声图像中的第一子区域(164)和第二超声图像中的第二子区域(182),第一子区域和第二子区域包括表示对象中相比于第一超声图像和第二超声图像内的解剖标志更靠近于透镜表面的组织的图像数据。例如,该对的帧166、168使用子区域164来互相关169。使用(一个或多个)最佳滞后用作的作为校正的通过例如将它们中的一个相对于另一个移动造成的帧166、168的逐对对准(或“共配准”)170充当Xu方案和/或Burlina方案的初始化171。初始化171改进针对任一方案的SNR。
作为另一范例,足够浅度的帧166、168的整个上部分上的局部互相关(LCC)的空间图中的值可以建议推迟Xu或Burlina方案的应用直到体外运动130足够地减弱。
对于一维换能器阵列而言,相对于对象122的体外运动130能够导致横向和轴向方向上的无意识运动,其效应163可以通过图像tj、tj+k中的一幅与另一幅的对准170收回,从而影响初始化171。
在二维换能器阵列的情况下,此外,无意识运动可以通过使用三维互相关在正交(即,海拔)方向上测量。
互相关系数可以是皮尔逊积矩互相关系数。其可以是标准化相关系数。在超声成像的领域中公知的,标准化互相关函数指示作为时间以及序列中的对应点的相对位移的函数的时间帧的序列中的对应点处的超声数据的互相关。互相关系数还可以简单地基于平方差的和。
操作性地,如在由图2中所提供的范例中看到,在呼吸时间ti处经由探头104采集图像(步骤S202)。当图像采集继续(步骤S204)时,指数i随着每幅图像被采集而增大(步骤S206)。
同时地开始PTX检测调节过程210和PTX检测预校正过程220。
还与图像采集同时,探头104可以抵着对象122的皮肤表面部分128手动地保持就位(步骤是S208)。因此,偶然的手运动230能够影响成像。即使探头机器地保持就位,便携式应用中的救护车运动240可以同样地引入相对于对象122的体外运动130。
在PTX检测预校正过程220中,图像tj和tj+k被选择(步骤S210)。此处,k可以是1,这表示所采集的下一图像,或k可以大于1;因此,在对象122的呼吸期间的不同的时间处采集两幅图像。电路108在形成对172的构成图像tj和tj+k之间的相似度的局部指示的空间图174中比较图像tj、tj+k的对172。在一些方面中,相似度的空间图基于超声图像的第一子区域与第二子区域之间的图像数据的互相关。在该背景下,局部指示可以是由使小局部区(诸如8个像素)互相关而导出的局部互相关系数(LCCC)176。图174的每个条目对应于被用于该条目的8个的中心像素。因此,LCCC 176是一种类型的相似度-指示符-幅度度量。经受比较(例如,互相关)的图像数据132是射频(RF)数据。然而,也可以使用其他数据,诸如,包络检测数据或者医学数字成像和通信(DICOM)图像数据。
如果校正是基于图像的(步骤S214),而不是基于运动传感器的,则处理路径取决于图像tj和tj+k的对172的浅部分要如何被识别。例如,比较相对对准可以包括比较来自运动传感器的针对超声图像中的第一子区域和第二子区域的位置数据。
如果校正基于图像并且基于成像深度134的预定义截止要被使用(步骤S216),则截止必须满足浅度准则,即,浅度的最小程度(步骤S218)。胸膜线152通常距成年人中的皮肤表面大约2厘米(cm)。因此,可以检查高于1.5cm或者甚至例如1.8cm的帧数据。截止因此在例如1.8cm或更少的患者独立和时间成像帧独立阈值TPFI 177处发生。任选地,能够未假定胸膜线152在特定深度处。
因此,如果预定义截止未被使用(步骤S216),则胸膜线152要被定位。处理取决于从图像tj和tj+k的对172的像素亮度检测是否被采用(步骤S220)。
如果像素亮度检测要被采用(步骤S220),则可以检测特征,诸如大约2cm标记或更低的高亮度线(步骤S222)。Burlina专利例如公开了将霍夫变换与一种形式的亮度阈值化组合。然而,其他解剖结构(诸如肌纤维)可以表现为亮线。
Xu申请中使用的另一方法基于图像tj和tj+k的图174和对172的过滤形式。如果该方法被使用,则处理遵循在Xu申请中所公开的技术(步骤S224)。为了降低处理,该方法可以以至少在一些情况下通过检测局部空间相似度指示符(诸如互相关系数176)中的急剧变化找到胸膜线152的简化形式来实施。利用对应于范围方向128的图174的一个维度,扫描在针对急剧变化的该方向上行进,并且过程通过在垂直于范围方向的方向上移动扫描来迭代地完成。急剧变化的集合横向地对应于胸膜线152。
在浅度的必要程度从刚才所讨论的方法中的一个或多个确定时,将LCCC 176或者更一般地使用来自针对对172的浅部分的相应的一个或多个相似度指示符的一个或多个幅度中间的幅度导出的至少一个相似度指示符幅度度量与基于图像的体外运动阈值TIBEM180进行比较178(步骤S226)。部分足够浅,使得其成像深度134处处小于胸膜线152的成像深度。部分通过选择较早帧166中的较早子区域164来定义。较早的子区域164仅包含比解剖标志(即,胸膜线)更浅的身体组织。阈值TIBEM 180可以具有靠近于一的值。
如果例如LCCC 176的预定义百分比或更多小于阈值TIBEM 180(步骤S228),则这指示实质性体外运动130存在。在这种情况下,选择较早帧166的先前描述的较早子区域164,并且选择较晚帧168的较晚子区域182,其大于并且涵盖所选择的较早子区域164(步骤S230)。两个帧166、168基于其相应的子区域164、182互相关(步骤S232)。由于子区域164、182比胸膜线152位于更浅,其中,身体组织保持无运动,互相关估计体外运动130对被相互比较的对172的构成图像之间的相对对准的效应163。应注意,备选地,较大的子区域可以从较早帧166选择,两个帧166、168对应地互相关。根据互相关找到一个或多个最小滞后(步骤S234)。两个子区域164、182然后被共对准170,从而在收回体外运动130的效应的程度上将两个帧166、168对准,并且因此将图像tj和tj+k的对172对准(步骤S236)。
如果另一方面校正是基于运动传感器的(步骤S214),则将传感器输出与基于传感器的体外运动阈值TSBEM 184进行相比较(步骤S238)。
如果传感器输出落在TSBEM 184之下(步骤S240),则对172的构成图像tj和tj+k共对准(或者“共配准”)170到收回经由(一个或多个)运动传感器120检测的运动的效应的程度(步骤S242)。
关于PTX检测调节过程210,当步骤S234中的一个或多个最小滞后或者步骤S238中的运动传感器输出变得可用时(步骤S244),其可以在其变得可用时实时与相应的运动阈值进行比较(步骤S246)。因此,任何基于图像或者传感器的阈值化可以利用例如与超声B模式扫描的帧速率相同的帧速率来执行。如果相应的第一运动阈值被超过(步骤S248),则肺滑动检测被暂停,或者任何这样的暂停被维持(步骤S250)。如果相应的第二运动阈值未达到(步骤S252),则肺滑动检测被启用,或者如果当前被暂停,要么实现要么考虑暂停的移除(步骤S254)。可以考虑体外运动130对于至少最小时间段是否已经足够低。第一运动阈值和第二运动阈值可以或可以不相同。在显示器115上,可以示出肺滑动检测的启用或者暂停的二元指示(步骤S256)。备选地或者额外地,其他用户反馈是可能的,诸如听觉反馈。反馈可以帮助用户处理探头104以保持其稳定。可以被显示的另一条目是其长度指示肺滑动的可能性的条(步骤S258)。可以根据隐藏马尔可夫模型确定肺滑动的概率,如在Burlina专利中所描述的。而且,例如被实施以显示LCCC 176的图174可以被示出在屏幕上并且实时更新(步骤S260)。步骤S256-260中的用户反馈中的任何可以实时更新,例如,以与超声B模式扫描的帧速率相同的帧速率。如果处理要继续(步骤S261),则返回步骤S244的分支被实现以重新迭代PTX检测调节过程210。
另一方面,如果基于图像的阈值TIBEM 180被满足或者被超过(步骤S228),如果传感器输出满足或者超过基于传感器的TSBEM 184(步骤S240),如果小于LCCC 176的预定义百分比小于阈值TIBEM 180,或者如果图像tj和tj+k的对172要收回体外运动130相对于对象122的效应163(基于成像(步骤S236)或者传感器输出(步骤S242)被逐对共配准170),则处理取决于肺滑动检测是否当前被启用(步骤S262)。其初始地默认被启用,但是经受根据步骤S248-S254的改变。如果其当前被启用(步骤S262),则已经共同配准到收回体外运动130相对于对象120的效应163的程度的图像的对172根据例如要么Xu申请要么Burlina专利经受对的构成图像之间的自动化比较,以自动地决定是否肺滑动在呼吸136、138期间在对象中发生。如果处理要继续(步骤S266),则返回步骤S210的分支被实现以重新迭代PTX检测预校正过程220。
PTX检测预校正过程220可以包括基于成像的预校正(步骤S214)、基于传感器的预校正或者两者。在三种情况中的任一种中,可以包括PTX检测调节过程210。
扫描器100被配置用于在除在步骤S208中保持探头104之外不需要用户交互的情况下自动地操作。
如从图2看到的,扫描器100可以被配置为实时操作。然而,如在Burlina专利中,所采集的成像能够被存储在存储器116中以用于之后由处理器118使用,此处在初始化肺滑动检测以及检测自身中。
在本文中提出的事物增加了医学紧急情况中的超声的无经验用户的自信并且提供鲁棒的PTX检测。
相对于被成像的医学对象的体外运动通过成像或者根据成像探头上的运动检测器被检测,并且从医学图像收回使得可以确定肺滑动是否存在,或者被测量以确定肺滑动检测是否要由于过度体外运动而暂停。对应于图像的相应子区域的图像子区域被选择用于图像到图像比较,使得选定的子区域仅包含相对于图像的采集中的成像深度比图像内的解剖标志更浅的身体组织。基于比较的结果,可以初始化涉及检查比标志更深的图像数据的肺滑动检测。运动传感器可以检测体外运动,并且基于其输出:逐对共配准图像到收回体外运动的效应的程度和/或确定是否暂停关于肺滑动是否在呼吸期间在对象中发生的决定。
应当注意,上文所提到的实施例图示而不是限制本发明,并且本领域的技术人员将能够在不脱离权利要求的范围的情况下设计许多备选实施例。
例如,帮助用户保持探头104稳定的听觉反馈可以包括语音指令。
在权利要求中,被放置在括号之间的任何附图标记不应被解释为对权利要求的限制。词语“包括”和其词形变化的使用不排除除权利要求中记载的那些外的元件或步骤的存在。在元件前面的量词“一”或“一个”不排除多个这样的元件的存在。本发明可以借助于包括若干不同元件的硬件并且借助于具有计算机可读存储介质的适合地编程的计算机和/或借助于具有机器可访问存储介质的集成电路来实施。计算机程序可以随时地、暂时地或者在较长时间段内被存储在适合的计算机可读介质(诸如光学存储介质或者固态介质)上。这样的介质仅在不是暂态传播信号的意义上是非暂态的,但是包括其他形式的计算机可读介质,诸如寄存器存储器、处理器高速缓存和RAM。尽管在互不相同的从属权利要求中记载了特定措施,但是这并不指示不能有利地使用这些措施的组合。

Claims (14)

1.一种被配置为检测肺滑动的成像系统(100),包括:
超声探头,其具有透镜表面;以及
处理器,其被配置为执行以下步骤:
选择第一超声图像中的第一子区域(164)和第二超声图像中的第二子区域(182),所述第一超声图像和所述第二超声图像包括表示胸膜线(152)的图像数据,并且所述第一子区域和所述第二子区域中的每个包括表示对象中相比于所述第一超声图像和所述第二超声图像内的所述胸膜线更靠近于所述透镜表面的组织的图像数据;
通过比较所述第一子区域与第二子区域之间的相对对准来估计对象中的体外运动的效应;并且
基于体外运动的所述效应的幅度来确定是否暂停对所述对象中的所述胸膜线处的肺滑动的成像。
2.根据权利要求1所述的系统,还包括被附接到所述探头的运动传感器。
3.根据权利要求2所述的系统,其中,所述的比较相对对准包括比较来自所述运动传感器的针对所述第一子区域和所述第二子区域的位置数据。
4.根据权利要求1所述的系统,其中,所述体外运动涉及在采集所述第一超声图像和所述第二超声图像所跨的时间段(143)内所述透镜表面相对于所述对象的皮肤表面的部分(128)的移动。
5.根据权利要求1所述的系统,还包括显示器(115),所述显示器被配置为显示所述第一子区域与所述第二子区域之间的相似度的空间图。
6.根据权利要求5所述的系统,其中,所述相似度的空间图基于所述第一子区域与所述第二子区域之间的图像数据的互相关。
7.根据权利要求6所述的系统,其中,所述处理器被配置为基于所述第一子区域与所述第二子区域之间的相对对准的第一运动阈值来确定是否暂停成像。
8.根据权利要求1所述的系统,还包括电路,所述电路被配置为发布指示是否由于体外运动而暂停对肺滑动的成像的用户通知。
9.根据权利要求8所述的系统,其中,所述用户通知包括二元指示符(S256)。
10.根据权利要求1所述的系统,其中,所述第一子区域(164)和所述第二子区域(182)被预定义为比与远离所述透镜表面的距离对应的深度阈值(177)更浅。
11.根据权利要求10所述的系统,其中,所述深度阈值(177)不大于1.8厘米。
12.根据权利要求5所述的系统,其中,所述空间图包括多个局部相关系数。
13.根据权利要求7所述的系统,其中,如果超过所述第一运动阈值,则所述处理器被配置为暂停对肺滑动的成像。
14.根据权利要求13所述的系统,其中,所述显示器被配置为基于所述第一子区域与所述第二子区域之间的所述相对对准来显示具有指示肺滑动的可能性的长度的条。
CN201680053990.2A 2015-09-17 2016-09-13 将肺滑动与外部运动区分 Active CN108024791B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201562219686P 2015-09-17 2015-09-17
US62/219,686 2015-09-17
PCT/IB2016/055430 WO2017046692A1 (en) 2015-09-17 2016-09-13 Distinguishing lung sliding from external motion

Publications (2)

Publication Number Publication Date
CN108024791A CN108024791A (zh) 2018-05-11
CN108024791B true CN108024791B (zh) 2021-09-07

Family

ID=56940117

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680053990.2A Active CN108024791B (zh) 2015-09-17 2016-09-13 将肺滑动与外部运动区分

Country Status (5)

Country Link
US (1) US11185311B2 (zh)
EP (1) EP3349664B1 (zh)
JP (1) JP6732896B2 (zh)
CN (1) CN108024791B (zh)
WO (1) WO2017046692A1 (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11134916B2 (en) * 2015-12-30 2021-10-05 Koninklijke Philips N.V. Ultrasound system and method for detecting pneumothorax
RU2740257C2 (ru) 2016-03-24 2021-01-12 Конинклейке Филипс Н.В. Ультразвуковая система и способ обнаружения скольжения легкого
US11446004B2 (en) * 2017-06-09 2022-09-20 Tokitae Llc Ultrasound systems and methods of identifying fluids in body regions using the same
WO2019034743A1 (en) * 2017-08-16 2019-02-21 Koninklijke Philips N.V. ULTRASONIC DETERMINATION OF A DYNAMIC AERONIC BRONCHOGRAM AND ASSOCIATED DEVICES, SYSTEMS AND ASSOCIATED METHODS
CN111542742B (zh) * 2017-11-27 2022-01-11 徕卡生物系统成像股份有限公司 载片架确定系统
EP3775986A4 (en) * 2018-04-09 2022-01-05 Butterfly Network, Inc. ULTRASONIC SYSTEM CONFIGURATION METHODS AND APPARATUS USING MULTIPLE IMAGING PARAMETER VALUES
WO2020020920A1 (en) * 2018-07-27 2020-01-30 Koninklijke Philips N.V. Devices, systems, and methods for lung pulse detection in ultrasound
US11627932B2 (en) * 2018-08-07 2023-04-18 Bfly Operations, Inc. Methods and apparatuses for ultrasound imaging of lungs
CN109597071B (zh) * 2018-12-29 2023-04-25 内蒙古工业大学 滑动式弧形阵列微波成像方法
JP7260884B2 (ja) * 2019-02-21 2023-04-19 ケース ウエスタン リザーブ ユニバーシティ 三次元(3d)超音波画像の処理
CN114007513A (zh) * 2019-07-10 2022-02-01 深圳迈瑞生物医疗电子股份有限公司 超声成像设备及检测b线的方法、装置、存储介质
WO2021055289A1 (en) * 2019-09-20 2021-03-25 Bard Access Systems, Inc. Automatic vessel detection tools and methods

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9025431D0 (en) 1990-11-22 1991-01-09 Advanced Tech Lab Three dimensional ultrasonic imaging
US5782766A (en) * 1995-03-31 1998-07-21 Siemens Medical Systems, Inc. Method and apparatus for generating and displaying panoramic ultrasound images
US6511426B1 (en) * 1998-06-02 2003-01-28 Acuson Corporation Medical diagnostic ultrasound system and method for versatile processing
US6589176B2 (en) 2001-12-05 2003-07-08 Koninklijke Philips Electronics N.V. Ultrasonic image stabilization system and method
WO2008073560A2 (en) 2006-10-06 2008-06-19 Verathon Inc. Systems and methods for lung imaging, pneumothorax detection and endotracheal tube insertion
US8798347B2 (en) * 2010-03-15 2014-08-05 Siemens Aktiengesellschaft System and method for image-based respiratory motion compensation for fluoroscopic coronary roadmapping
EP2723241B1 (en) * 2011-06-27 2014-11-19 Koninklijke Philips N.V. Ultrasound-image-guide system and volume-motion-base calibration method
EP3363365B1 (en) 2011-12-12 2022-02-09 Koninklijke Philips N.V. Automatic imaging plane selection for echocardiography
US9743857B2 (en) * 2011-12-16 2017-08-29 The Regents Of The University Of California Digital manometry finger-mountable sensor device
US20130184584A1 (en) * 2012-01-17 2013-07-18 Richard E. Berkey Systems and methods for computerized ultrasound image interpretation and labeling
US8914097B2 (en) * 2012-01-30 2014-12-16 The Johns Hopkins University Automated pneumothorax detection
US9554778B2 (en) * 2013-11-08 2017-01-31 Siemens Medical Solutions Usa, Inc. Responsive power saving in ultrasound
JP6557328B2 (ja) 2014-09-25 2019-08-07 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 自動気胸検出のための装置および方法

Also Published As

Publication number Publication date
EP3349664B1 (en) 2020-03-11
US20180344293A1 (en) 2018-12-06
JP2018527112A (ja) 2018-09-20
WO2017046692A1 (en) 2017-03-23
US11185311B2 (en) 2021-11-30
CN108024791A (zh) 2018-05-11
EP3349664A1 (en) 2018-07-25
JP6732896B2 (ja) 2020-07-29

Similar Documents

Publication Publication Date Title
CN108024791B (zh) 将肺滑动与外部运动区分
US20200222027A1 (en) Device and method for automatic pneumothorax detection
US10929951B2 (en) Anatomically intelligent echocardiography for point-of-care
JP7268087B2 (ja) モデルベースセグメンテーションを用いた像撮像誘導
US11191518B2 (en) Ultrasound system and method for detecting lung sliding
JP5498299B2 (ja) 2次元超音波映像に対応する2次元ct映像を提供するシステムおよび方法
US10347035B2 (en) Diagnostic image generation apparatus and diagnostic image generation method
KR20140096919A (ko) 의료 영상 정합 방법 및 장치
JP2015500083A (ja) 心エコー検査のための自動的な画像化平面選択
Kattea et al. Differentiating pneumothorax from the common radiographic skinfold artifact
CN111557692B (zh) 目标器官组织的自动测量方法、超声测量装置及介质
JP5885234B2 (ja) 疾患判定装置および超音波画像形成装置における画像解析方法
CN112801957A (zh) 一种基于超声应变成像的气胸自动检测系统
KR20110102105A (ko) 인체 내 오브젝트를 스캔하는 3차원 초음파 검사기 및 3차원 초음파 검사기 동작 방법
CN114521912A (zh) 用于增强胸膜线的可视化的方法和系统

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