CN108024791B - 将肺滑动与外部运动区分 - Google Patents
将肺滑动与外部运动区分 Download PDFInfo
- 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
Links
- 230000002685 pulmonary effect Effects 0.000 title description 3
- 210000004072 lung Anatomy 0.000 claims abstract description 55
- 238000003384 imaging method Methods 0.000 claims abstract description 52
- 239000000523 sample Substances 0.000 claims abstract description 26
- 230000000694 effects Effects 0.000 claims abstract description 21
- 238000002604 ultrasonography Methods 0.000 claims description 30
- 238000000338 in vitro Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 abstract description 39
- 230000029058 respiratory gaseous exchange Effects 0.000 abstract description 8
- 201000003144 pneumothorax Diseases 0.000 description 41
- 238000000034 method Methods 0.000 description 23
- 210000004224 pleura Anatomy 0.000 description 14
- 230000008569 process Effects 0.000 description 14
- 238000012937 correction Methods 0.000 description 11
- 238000012545 processing Methods 0.000 description 9
- 210000001519 tissue Anatomy 0.000 description 9
- 210000000038 chest Anatomy 0.000 description 6
- 239000013598 vector Substances 0.000 description 6
- 239000000470 constituent Substances 0.000 description 5
- 230000002411 adverse Effects 0.000 description 4
- 208000014674 injury Diseases 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000008733 trauma Effects 0.000 description 3
- 238000012285 ultrasound imaging Methods 0.000 description 3
- QLRRUWXMMVXORS-UHFFFAOYSA-N Augustine Natural products C12=CC=3OCOC=3C=C2CN2C3CC(OC)C4OC4C31CC2 QLRRUWXMMVXORS-UHFFFAOYSA-N 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000001936 parietal effect Effects 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 239000000725 suspension Substances 0.000 description 2
- 206010051093 Cardiopulmonary failure Diseases 0.000 description 1
- 206010038687 Respiratory distress Diseases 0.000 description 1
- 208000027790 Rib fracture Diseases 0.000 description 1
- 206010052428 Wound Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000011976 chest X-ray Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000003601 intercostal effect Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005399 mechanical ventilation Methods 0.000 description 1
- 210000001087 myotubule Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 210000003281 pleural cavity Anatomy 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 210000000115 thoracic cavity Anatomy 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 230000009278 visceral effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/54—Control of the diagnostic device
-
- 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/085—Detecting 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/42—Details of probe positioning or probe attachment to the patient
- A61B8/4245—Details 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/4254—Details 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- 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
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
- A61B8/5276—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
- G01S15/892—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array the array being curvilinear
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details 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/52077—Details 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所述的系统,其中,所述显示器被配置为基于所述第一子区域与所述第二子区域之间的所述相对对准来显示具有指示肺滑动的可能性的长度的条。
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)
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)
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. | 自動気胸検出のための装置および方法 |
-
2016
- 2016-09-13 JP JP2018514320A patent/JP6732896B2/ja active Active
- 2016-09-13 US US15/760,244 patent/US11185311B2/en active Active
- 2016-09-13 CN CN201680053990.2A patent/CN108024791B/zh active Active
- 2016-09-13 WO PCT/IB2016/055430 patent/WO2017046692A1/en active Application Filing
- 2016-09-13 EP EP16766646.0A patent/EP3349664B1/en active Active
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 |