CN103298412A - 使用经胸肺部多普勒超声测量肺部血压 - Google Patents

使用经胸肺部多普勒超声测量肺部血压 Download PDF

Info

Publication number
CN103298412A
CN103298412A CN2011800505316A CN201180050531A CN103298412A CN 103298412 A CN103298412 A CN 103298412A CN 2011800505316 A CN2011800505316 A CN 2011800505316A CN 201180050531 A CN201180050531 A CN 201180050531A CN 103298412 A CN103298412 A CN 103298412A
Authority
CN
China
Prior art keywords
patient
lung
pulmonary
power
feature
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.)
Granted
Application number
CN2011800505316A
Other languages
English (en)
Other versions
CN103298412B (zh
Inventor
尤伦·帕提
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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Publication of CN103298412A publication Critical patent/CN103298412A/zh
Application granted granted Critical
Publication of CN103298412B publication Critical patent/CN103298412B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/06Diaphragms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/04Measuring blood pressure
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • 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/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • 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
    • 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/5284Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving retrospective matching to a physiological signal
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Hematology (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

可以通过将超声能量传送到患者的肺内、及检测肺内的血管和围绕该血管且充满空气的肺泡之间的运动边界引起的反射超声的多普勒频移来分析患者心或肺的运作。边界的运动由血管中的压力波引起的,该压力波会导致那些血管的直径发生变化。将探测到的多普勒频移用于在多个不同气压水平中的每个气压水平处获取功率和速度的数据,然后基于该获取到的数据对患者的肺部血压进行估算。

Description

使用经胸肺部多普勒超声测量肺部血压
交叉引用相关的申请
此申请要求2010年10月21日提出的美国临时申请61/405,454的利益,该项美国临时申请通过引用合并于此。
背景技术
将超声多普勒用于频谱测量动脉及静脉中的血液的流速已经是很既定的做法。进行这种测量的一种普遍使用的规程基于三个典型阶段:用超声图像对目标面积进行初始识别标识(其中流量待测);将标记放在该图像适合的位置处;及将回波装置从图像模式切换到频谱多普勒检查模式,以实时显示流速。此规程可以用于,例如,测量肺静脉中的血流量。
另一相对新的规程用于经颅多普勒(TCD)测量及一些外周血管研究。在此规程中,超声束不依赖成像,直接瞄准目标的已知位置。由于人类颅骨的结构和位置及其成分相对固定和已知,仅通过此规程的回波多普勒(即无需成像)即可对大脑基部特定的血管,诸如威利斯环的动脉,进行研究。无需成像即可对流速进行测量的事实使得穿过颅骨的骨骼进行的测量成为可能,其中颅骨的骨骼会使得超声束衰减及散射到不能获得实际图像的程度。
虽然经颅多普勒测量现在常规使用于研究脑内结构,但迄今为止人们却不认为此项技术有可能运用到经胸监测肺部血管。这是因为肺包含大量气泡,而这些气泡相比骨骼,对超声束的衰减和散射要严重得多。鉴于这一点,除了不会被肺组织掩蔽的初始的较大的肺部血管段以外,肺血管内的动脉和静脉流速及肺组织本身都还未能通过多普勒超声进行研究。
发明内容
本发明的一个方面涉及一种估算患者肺部血压的方法。所述方法包含的步骤是:用经胸肺部多普勒超声,在多个不同气压水平下从患者的至少一个肺中连续地获取功率和速度的数据。之后,基于所获取到的数据,估算患者的肺部血压。可选地,可以从患者肺内的至少两个不同位置处获取功率和速度的数据。可选地,该估算步骤包含识别至少一个气压水平,在所示至少一个气压水平下,总功率会接近零或下降至在小于气压水平不升高时获取到的总功率的10%。
本发明的另一方面涉及一种测量患者肺部血压的装置。此装置包含一个压力传感器和一个转换器,压力传感器被配置成用于测量患者至少一个肺内的气压,转换器被配置成用于将超声能量传送进患者至少一个肺内的目标区域、探测从该目标区域反射的超声能量及生成基于该探测到的超声能量的输出。此装置还包含一个多普勒信号处理器,所述处理器被配置成用于处理转换器的输出、在多个不同气压水平中的每个气压水平处从患者的至少一个肺中连续获取功率和速度的数据及基于获取到的功率和速度的数据对患者的肺部血压进行估算。可以在此具体实施方式中执行上文描述的可选项。
本发明的另一方面涉及一种确定患者肺部血压水平的方法。此方法包含的步骤是:将超声能量传送进患者的至少一个肺内、检测由至少一个肺内的血管和围绕所述血管且充满空气的肺泡之间的运动边界引起的反射超声能量的多普勒频移、改变肺内的气压、监测检测到的多普勒频移如何响应于压力改变而发生变化及基于在监测步骤中监测到的变化确定患者的肺部血压水平。
本发明的另一方面涉及一种确定患者是否患肺动脉高压的方法。此方法包含的步骤是:使患者至少一个肺内的气压升高到预计会使未患肺动脉高压的患者血流量下降的水平、在气压升高时从患者的肺中获取至少一组功率和速度的数据、及基于在获取步骤中获取到的功率和速度的数据确定总功率是否高于阈值。
附图说明
图1是经胸肺部多普勒(“TPD”)系统的具体实施方式的框图。
图2绘示了由图1的系统生成的输出的实施例。
图3是对图2中示出的输出中五个特征的示意性描述。
图4A绘示了临床多普勒测量的“经典模型”。
图4B绘示了用TPD收集到的多普勒信号的原始位置。
图5A将正常人的TPD输出和对肺部动脉、静脉内血液流速的描记进行了比较。
图5B到E是正常呼吸和不同呼吸动作过程中的TPD输出。
图6绘示了来自正常人的超过十次心动周期的平均TPD输出。
图7A绘示了之后有扩散房性期前收缩的正常窦性心律的TPD输出。
图7B绘示了呈现非扩散房性期前收缩时的TPD输出。
图8绘示了呈现期前收缩时的TPD输出。
图9绘示了发生房颤时的TPD输出。
图10A到C描述了TPD输出的三个特征的平均峰值正流速和负流速的实验数据。
图11A是正常人和异常人之间的速度差别的图解表示。
图11B是正常人、慢性阻塞性肺病(COPD)和纤维化患者之间功率差异的图解表示。
图12是施行肺部血压测量的系统的框图。
图13绘示了TPD信号如何响应于肺气压增加而进行变化。
图14绘示了TPD信号如何在不同肺气压下发生差异。
图15绘示了TPD信号如何响应于肺气压内的变化而进行变化。
图16绘示了TPD功率水平如何响应于肺气压内的变化而进行变化。
图17A绘示了功率水平在两种不同压力下到达了零。
图17B绘示了正常人的功率读数。
图17C绘示了肺动脉高压患者的功率读数。
图18绘示了由自动特征识别算法确定的特征之间的边界。
具体实施方式
发明者已经认知在大量心肺病理生理条件和疾病中,可以对肺循环和肺部光散射性质进行重大修改,而且这样的信息对于诊断和治疗相当重要。此处描述的具体实施方式被设计成用多普勒超声监测肺内动脉和静脉的功能性,及监测围绕肺内动脉和静脉的肺组织的完整性和功能性。本文中将其称为“经胸肺部多普勒”或“TPD”。
图1是一幅此类具体实施方式的框图。与探针11(包含超声转换器)连在一起的多普勒超声机12被用于以常规方法随着时间的推移,在每个相关速度下确定实验对象10的目标区域内的功率。这可以通过以下步骤完成,这些步骤有:生成脉冲超声束、收集反射能量、计算多普勒频移及处理由此获得的数据,从而提供超声反射器的功率及对应的速度的矩阵。相配的多普勒超声机12的一实施例是Sonara/tek脉冲经颅多普勒仪(可从美国威斯康辛州麦迪森的Viasys公司得到),这是一种脉冲多普勒系统。多普勒超声机12将其捕捉到的数据传送到个人计算机13,个人计算机13载有生成常规多普勒超声显示的软件,而在常规的多普勒超声显示(如与计算机13相连的监控器)中,x轴表示时间、y轴表示速度,功率用色彩表示。控制超声参数的适合的软件也可以从Viasys处得到。应注意到,在可替换的具体实施方式中,多普勒超声机12和个人计算机13的功能可以合并入单台设备。
较佳地,还提供ECG系统14。ECG系统14与常规的ECG导联15连接,并以任何常规方式生成输出。该输出较佳的是与多普勒超声机12在时间上同步,以便ECG和超声两者的显示内容可以显示在同一时标上。ECG系统14的输出以任何常规方式被提供至个人计算机13。在可替换的具体实施方式中,个人计算机13可以被换成与多普勒超声机12组合在一起。
探针11可以是标准TCD探针,所述标准TCD探针诸如直径21mm、焦距为4cm的2MHz的传感器。适合的探针可从Viasys处得到,并与Viasys的Sonara/tek机一起使用。还可以使用用来对外周或心脏血管进行多普勒超声测量的常规探针。然而这些应用通常使用窄束以提供较高的空间分辨率,而较高的空间分辨率有助于制作相对小的目标的几何表征,其中所述窄束经常由相控阵换能器定形。虽然这些窄束在TPD的背景下可以产生有用的结果,但一些较佳的可替换具体实施方式使用相对较宽的束,例如有效横截面至少为1/2cm2(例如,在1/2到3cm2之间)的束。用较小的转换器、及用单组元的转换器即可使用相对较宽的束,而无需使用其他解剖应用中普遍使用的相控阵换能器。肺包含有不特定几何形状的、由血管(动脉和静脉两者)及其周围的肺组织组成的相对较大的集合体这一事实,当使用较宽束时可为系统所利用。
应注意到,由于散射会使得无法用超声对肺进行成像,因此人们不得不在仅有已知的解剖学构造而无其他指引的情况下扫描目标。还应注意到,散射会削弱通过相控阵或通过机械手段进行扫描的优势。而且,由于肺的整体深度范围都会引起散射,因此对于在肺部的应用而言,CW(连续波)超声比PW(脉冲波)多普勒超声的效果要差。因此,一些较佳具体实施方式运用带有宽束的PW超声。可选地,此类具体实施方式可以采用设置在本体表面上的复式传感器。
可选地,可以使用特别选取或设计的超声探针和/或适合的束功率控制,该束功率控制包含动态可调的束形状和尺寸,以便能对变化的组织容量进行测量。应注意到,与将多普勒用于其他组织目标的情形相比,此处来自相对较大容量的平均信号或完整信号包含了非常有价值的信息。
除了用于从多普勒信号生成显示的标准软件外,个人计算机13较佳地包含用于激活TPD及用于选择期望的操作模式、显示模式和存储模式的软件。该个人计算机13还包含或具有通往适当的数据存储资源(例如,本地或远程的硬盘)的路径。该个人计算机13较佳地用一种或多种降噪算法(NR)处理速度和功率与时间之比的数据,该一种或多种降噪算法能最优化地使肺组织引起的信号散射和衰减产生的噪音减到最小。
降噪的一种较佳途径含有两个阶段——平均化和边缘检测。在第一阶段,通过平均化N特征信号的功率/速度来获取来自多个心动周期的平均的信号,该处每个N信号较佳地代表单个心动周期。N较佳的是4到20之间的整数(例如,10)。较佳地,每个信号在各端都以R-波为界,尽管在可替换的具体实施方式中可以将心动周期上的其他点用作时间参考点。计算出来的平均信号被假定具有患者的声谱图特性,并作为之后要确定的相关特征的基础。应注意到,虽然施行此平均化阶段是优选地,但在可替换的具体实施方式中,可跳过此阶段,并可以用来自单个心动周期的数据施行后续步骤。
第二阶段是边缘检测和包络计算。在此阶段,考虑到振幅和时间,我们线性的描绘功率及速度信号描记与时间之比,并由此将代表血管运动(即,信号)的分段与代表噪音的分段隔开。在这个阶段中可以使用一种或多种降噪算法。在一较佳实施例中,将两种特定边缘检测算法应用于数据,此处称为算法A和算法B。算法A和算法B都被应用到平均化的信号上,并计算经平均化的图像中的信号和噪音之间的边缘(即,包络)。
算法A是本地的、一维的方法,其中在给定时间处信号和噪音之间的边缘(eA)仅根据此时间近端数据的统计情况进行定义。此算法包含两个步骤:第一步骤中,在任何给定的时间(ti),我们通过在ti近端内寻找最低能量的区域来定义每个功率频谱A(ti)的阈值'thr(ti)'。接着,我们将thr(ti)设置成与此区域内的最高能量水平相等。然后,我们将thr(ti)应用到A(ti)上,并将超过thr(ti)的所有A(ti)的部分视作相当于运动区域,及将所有其他的部分视作相当于噪音。
在算法A的第二个步骤中,我们用噪音的统计情况改进流量和噪音之间的初始区别:在此步骤中,我们假定向下估算(流量包含在噪音区域内);调节包络检测以从噪音区域中排除流量像素点;并通过流像素点的相对较高的值在噪音区域内对它们进行识别。象征性地,该步骤可以用以下三个步骤表示:
(a)对于每个t={l,2,..,N},计算P(t)={噪音区域内的A(t)的平均值}
(b)定义阈值'thr2',该阈值基于{P(l),P(2),...P(N)}的平均值和标准差
(c)在P(t')>thr2的每个t'处,通过升高上包络或降低下包络来减小P(t'),直到P(t')<=thr2。
为了取得更好的效果,较佳地,将步骤(a)-(c)重复多次(例如,10次)。
算法B是一种边缘检测算法,该算法将数据看作二维图像。在此方法中,信号被看成一种物体,该物体被由其分割出的噪音所围绕,并且据此计算出边缘(eB)。此分割方法是Chan-Vese算法的一种实施。(参见Chan T.F.,Vese L.A.,Active contours without edges(无边缘的活动轮廓)Image Processing IEEE,刊登于卷10,期刊2:266-277(2001年2月)的,该内容通过引用合并于此。)
之后,将通过算法A(eA=[eA(tl),eA(t2),...])计算出来的边缘与通过算法B(eB=[eB(tl),eB(t2),...])计算出来的边缘组合。组合这二种边缘的一种适合方法是假定目标边缘在被发现的两种边缘之间穿过。该组合可以通过多种方法实现。一种方法是取每个点上的从算法A和算法B中获得的结果的简单平均值。组合这两种边缘的另一种方法是设计一个重量(w=[w(tl),w(t2),...])的阵列,如下所示:(1)使间隙处图像的功率水平沿着时间整合;(2)将结果线性转化成具有最大值'1'及最小值'0';及(3)随后,在一个时间点ti上的用于边缘的输出通过下述等式定义:e(ti)=w(ti)*eA(ti)+(1-w(ti))*eB(ti)。
较佳地,结果性输出经由一维中滤波器(例如,等级3)进行平滑化及显示,及图2绘示了结果性输出的一个实施例。应注意到,在可替换的具体实施方式中,可以仅使用一种算法(即,算法A或算法B或一种不同的牛顿拉夫森算法),或者单独采用或者与其他牛顿拉夫森算法结合采用。
图2绘示了用2MHz的多普勒超声系统获得的正常人右肺中的超声反射器的速度22,该多普勒超声系统的探针放置在胸骨右侧约3cm、剑状骨(约第四个肋间隙)尖端水平线上方7cm处。超声束大致垂直于胸部表面。图2中,较暗的区域对应于较高的功率。图2的底部还可以较佳地显示出常规的ECG24。从达14cm的深度(出入口)处的记录中和从未被心脏支配的左肺区域中获得了类似的记录。在表面以下8-9cm深的地方记录到整个右肺的最大信号强度。
在常规TCD系统(即,3-10kHz)中使用的同样的脉冲重复频率(PRF)也可以用于TPD系统。然而,TPD声波图22包含多个中速信号,该多个中速信号与心动周期具有同样的周期性,并且通常只会达到约30cm/sec的值。鉴于这些相对较低的峰值速度(与较大动脉中多普勒流动测量相比),可以将TPD REF的值设定得低于标准脉动多普勒系统的值。通过将REF降到1到3kHz之间,可以增加有效束穿透深度。这是非常重要的,原因在于,肺内的超声速度比脂肪、肌肉等中的超声速度要低30-50%,从而会降低有效束穿透深度。较佳地,软件被配置成能算入此较低的速度。可以通过识别肺信号出现的最浅的点(即,带有很大回返的信号)来探测到肺内信号发起的临界点。应注意到,对不同肺深度的测量会产生非常类似的描记,且对于其他明显是正常人的描记也具有大体上相似的特性。
可以看出,在各极性处(正极或负极),通常可以用相对较高的能量和大致三角形的形状来标识出五个重要特征。这五个特征在图3中被示意性地描绘出,并编号为#1-5。这些特征中的每个特征都包含正分量(即,正速度,指示出流动的方向朝向探针),及相应的负分量(即,负速度,指示出流动的方向远离探针),并带有很大程度的正/负对称性。由此,这些特征中的每个特征都指示出相反方向上的同时发生的运动。如图3中所示,这些特征与心动周期(注意ECG24中的R波26)是同步的。
操作理论
从整个肺记录到的上述信号似乎具有不一样的来源。众所周知,肺由许多肺泡管、肺泡囊和肺泡组成,这许多的肺泡管、肺泡囊和肺泡可以被认为是被非常薄的膜密封起来的微小气体容积。肺泡能被假定成由球状体合理地代表,该肺泡的尺寸范围是50到150μ。当暴露在超声波中时,这些天然的肺组成部分在许多方面都非常类似于用在声纳成像中的超声造影剂。(超声造影剂是充满气的微泡并具有高度的回声反射性,即,物体反射超声波的能力。)肺泡和软组织之间的回声反射性区别很大,因此大多数能量会被反射。
尽管散射会使得无法获取肺结构的超声图像,但其确实有助于检测软组织和肺泡之间的高反射边界的运动。呼吸作用甚至心脏收缩及行进在血管和肺部血管中的机械脉冲波引起了此边界运动。众所周知,肺部血管具有很高的顺从性(即,比体循环的血管的顺从性大得多),且围绕该血管的肺泡组织中充盈的空气是高度压缩的。由此,肺部动脉和静脉中的压力波会造成它们的直径发生巨大的变化。这些变化反过来会使高反射边界移动,从而使在其附近的肺泡、肺泡囊等压缩并移动。由于在组织内和在空气中的超声传播速度差别很大,所以在组织和空气的边界处就会出现机械耦合不匹配,该机械耦合不匹配会导致高的回声反射性和很强的超声反射,而这种情况下是来自引起多普勒频移的运动反射器。这些反射通常会在噪音水平以上约100dB(与从动脉中流动的血液中测得的通常强度相比,后者的在噪音水平以上约30到40dB范围内)。因为这些信号如此之强,所以即使回返信号中有部分被固定肺组织层掩蔽,多普勒系统还是可以收集到回返信号,其中所述掩蔽会使超声能量衰减约40dB/cm。
图4A和图4B绘示了常规多普勒信号和用TPD手段通过胸壁收集到的信号之间的差别。图4A描述了临床多普勒测量的“传统模型”,其中设备测量了动脉和静脉血流42产生的,或更具体的,红细胞43经由那些血管44进行运动(会反射超声波)而产生的多普勒频移。
图4B描述了用TPD收集到的多普勒信号的原始位置。此处,压力中的变化会导致血管直径变化,其原因在于,随着心跳产生会驱使血液32穿过血管的压力脉冲,血管壁34会暂时性地向外膨胀,并压缩环绕其周围的充满空气的肺泡、肺囊泡等35。由运动的血管-肺泡边界引起的反射超声的多普勒频移通过TPD系统被转化成功率-和-速度与时间之比的图,并通过TPD系统被显示。可以预期的是,这些信号中的大部分由小尺寸及中等尺寸的动脉和静脉产生。以此种模式产生的信号的独特特征(与身体其他部分中的血液流动产生的信号相比)在于其双向性。此现象可能是因为肺实质从各侧包围血管,使得无论相关束方向在哪个方向上,较近的边界都会朝向束源运动而远处的边界都会远离束源运动。结果,产生了相反极性的类似信号。有时候,如图2所描绘的,信号看起来几乎完全对称。这样的对称在非肺部血流的记录中非常罕见。
值得注意的是,对流经血管的血流进行常规的多普勒测量时,其中,运动的是血流本身,探针要被放置成能使超声束尽可能地与流动轴线平行以获取最大速度。相比之下,此处所述的实现TPD测量的运动是垂直于血流的方向的,因此最优化的位置是正交于流动轴线且平行于血管半径。但是由于肺内有众多的血管,所以在TPD的背景下,放置就不那么关键(与对经由血管的血流进行常规多普勒测量相比)。
由于图2中的特征总是具有对应于ECG24的R-R间隔的重复周期,所以我们可以得出结论:它们必定源自会反射超声能量而又与心跳同步运动的结构。这些实体可以是心脏本身、肺部血管中流动的血液、搏动的血管或其与肺泡、肺泡囊、空气等的连接处。
记录到的信号将被称为肺多普勒速度信号(LDVS)。图5A将正常人的典型LDVS52与肺部动脉、静脉内血压流速的描记53、54进行了比较,比较了单个心动周期的情况,且心动周期的持续时间被均一化到了同样的时标(注意ECG24的R-波26)。呈现出了显著的关联性。图5B-E将正常呼吸(图5B)的LDVS56与不同呼吸动作过程中记录到的LDVS56进行了比较,比较了多个心动周期。例如,在FRC(功能余气量)处屏气的过程中,特征57具有正常的形状和速度,但衰减了的强度。在瓦尔萨尔瓦动作(图5D)过程中,其中胸腔压力大幅提升,特征58看起来几乎消失了。相比之下,在穆勒动作(图5E)过程中,在胸腔内产生了负压力,LDVS59的速度和信号功率都得到了增强。
五个特征(#1-5)与心跳及与关联的机械事件的同步性表明信号源与心脏和血管产生的脉动是相关的,及呼吸动作对特征的强烈调整(参见图5C-E)表明肺实质的状态会严重影响这些特征的形状。尽管肺实质有很强的机械倾倒性质,但事实是,能排出心脏和较大的血管的参与而对整个肺都记录到了类似的信号。由此,最有可能的是,脉动是通过在肺内沿血管传播,包含沿相对较小的血管传播。
基于上述的操作理论,我们对图2和3中描绘的五个特征作出以下诠释:特征#1,一般会很显著,紧随着R波出现,且与收缩期心室收缩相符。特征#2,具有较低的峰值速度,与ECG的T波、复极化及心室松弛相符。特征#3经常是双峰隆起,并且会延续相对较长时间,看上去主要出现在舒张快速充盈期。特征#4,具有典型的低的峰值速度,对应于心舒张后期,特征#4的后部经常会与可探测的信号无关。特征#5,一般会有较高的峰值速度,与心房收缩相符。
这五个特征的相对振幅、上升时间及下落时间、持续时间等能提供关于如下方面的信息:血流血液动力学、多种心血管系统成分的被动机械性质及活动(收缩)力。此外,显示内容提供的信息主要与肺系统相关。
为了证实回返是由运动的组织-空气边界产生的理论,用模型制作了多普勒声波图,其中混合微型空气气泡(小于0.5mm)的假血(多普勒试验流体707,美国康涅狄格州ATS实验室公司)在适当的管中流动。在声波图中,气泡表现为明亮的“光点”。流动的假血和气泡的功率图谱揭示了如下现象:由运动的空气气泡产生的峰值功率比在类似情况下记录到的流动的假血和冠脉血流的峰值功率高约40dB。这些结果与上述的理论一致。
对10名正常的、年龄在27-72岁的志愿者的右肺或左肺进行了测量,测量方法是将超声传感器放置在坐着或者仰卧着的志愿者的胸壁上方。通过标准的超声耦合剂,将焦距为4cm的21mm、2MHz的传感器与胸壁阻抗匹配。用脉冲TCD设备(美国威斯康辛州麦迪森的Viasys公司)在脉冲重复率(PRF)为3kHz时,在胸壁上方的不同进行测量。传送的脉冲功率达到允许的最大ISPTA.3(492mW/cm2)的10%。这些对象被连接至标准的三导联ECG(以色列约克尼穆的Norav医药公司),三导联ECG的输出包含在了显示内容内。
对得到的速度-和-功率与时间之比的描记进行的观察能提供当那些描记偏离所期望的正常描记时,大体上在特定位置上的对肺实质的机械性能的诊断信息。这些诊断信息可以包含与组织结构(可能与肺气肿、肺纤维化、肺不张等相关)、脉管系统、或肺泡内或周围存在流体(如患充血性心脏衰竭或肺炎、诸如栓塞及溢血的血管事件)等相关的信息。这些从正常描记偏离的偏差起因于弹性性质、多种组织成分的质量及其空间分布上的变化。这样的变化会导致描记的一些方面发生全局或局部的相应的变化,所述描记的一些方面包括描记的功率图谱轮廓、时间常数、持续时间或振幅(相对的或绝对的)方面。诸如深呼吸、用力呼气、屏气、瓦耳萨耳瓦动作、运动等的生理操作,可以被用于提升诊断能力。应注意到,当反射自任何肺内成分的超声波穿过介于其与胸壁之间的肺实质时,该被反射的超声波会被修改。此组织表现得像具有具体特性的机械过滤器。这些特性取决于相关实质的状态,以至如Gavriely N.,Y.Palti & G.Elroy(Spectral Characteristics of Normal Breath Sounds,J.Appl.Physiol(正常呼吸音的图谱特性,应用生理学杂志.)50:307-314(1981),该项文献通过引用合并于此)所描述的穿过此过滤器的信号的功率谱反映在声音信号的过滤特性上。
可选地,可以使用ECG24的R波26作为参照点,对来自单个实验对象的信号取多个心动周期的平均值。例如图6,描绘了来自正常人的十个心动周期的平均值62,记录的是右肺。能够看到对应于上述特征#1-5的五个特征#61-65。对于其他正常人进行的描记也大致类似。
对心脏功能的探测和表征
此处所述的TPD系统的一项有益运用是作为一种工具,该工具通过对肺的TPD测量间接探知心血管系统的功能。这可以实现的原因在于,心脏活动的产出沿肺部血管从其在心脏内的原始位置扩散到整个肺容量。在这种途经中利用TPD可以探知并表征一些与正常机械心脏活动偏离的临床上的显著偏差,以下将给出一些实施例。
图7A描述了在由房性期前收缩引起的心律不齐情况下,肺信号与正常模式相比发生的变化,房性期前收缩是一种附加的异常心脏收缩,。图7A的左侧描绘了正常窦节律的典型信号,及右侧描绘了扩散到心室的房性期前收缩71(即,由窦房结产生的早期电跳动生成的信号)的出现。这些信号基本上是正常的节律集合体的复制,即这些信号包含额外的心房收缩(特征#5),跟着是额外的心室收缩(特征#1)和心室放松(特征#3)。当它们发生得足够早,心房收缩信号(特征#5)可能会及时地叠置在之前的心室放松(特征#3)上。图7B描述了由房性期前收缩73产生的信号的特性,房性期前收缩73使得心房紧张(特征#5)不会从心房扩散到心室,如异常附加的特征#5*之后无特征#1和#3所表明的。
图8描述了心室内电异常活动82生成的期前收缩紧张(特征#1*)所产生的信号。图9描绘了对应于罹患房颤的患者的心室收缩的原始位置(#1)的信号。因为没有特征#5(代表心室收缩),所以这种情况在图9中是显而易见的。同样可见的是,心室期前收缩的特性与心房期前收缩的特性完全不同,反映出机械活动本质之间的有较大差异。这样记录的描记会有助于内科医生确定异常活动进行扩散的路径。
上述讨论的有关图7A、7B、8和9的任何异常特征的出现,可由此用作患者具有相应问题的指示。这些通过依据视觉即可实现,也就是通过观看显示内容及识别相关特征。在可替换的具体实施方式中,可以使用模式识别软件以自动识别相关特征。
多位置测量
可以在不同的肺深度上进行TPD测量,且这样的测量一般会显示出非常类似的描记,从而指示出信号在肺容量中的大范围传播。也可以在对象身体的不同位置进行测量,诸如在肋间隙(例如,在第二根肋骨和第三根肋骨之间,或在第五根和第六根肋骨之间),及在肋骨上方的位置。当在多个位置进行这样的测量时,一些情况下,每个位置采得的测量结果在信号形状、速度、功率方面会有显著的差异。发明者已经认知到,这些记录全体和具体的记录差异可以用于辅助诊断特定的生理情况。
在一实施例中,在三个不同位置,对两名慢性阻塞性肺病(COPD)患者的右肺进行测量,测量位于每名患者的右肺上:上部区在第2到3根肋骨的水平,中部区在第4根肋骨的水平,下部区在第5到6根肋骨的水平。与正常人在上、中、下位置测得的数据非常相似的情况有所不同,在COPD患者中,上部区的信号要远远小于中部区的信号,中部区的信号接下来会远远小于下部区的信号。而且,在不同区的信号形状(例如,速度和功率的对称度)也有差异。可以用这种偏离于正常状况的偏差来预知COPD的出现。类似地,可以用其他偏离于正常状况的偏差来预知其他异常情况的存在。
在这三个位置(即,上部、中部和下部)的每个位置上,对一组患者(包含正常的患者、COPD患者、肺结节患者、肺纤维化患者)进行测量来得到特征#1,3,和5的平均峰值正速度和负速度。该实验数据在图10A-C中有描述,正、负速度在y轴线上。正常患者的数据位于左侧,在FS和DUL之间的患者患有COPD,在BAD和BUJ之间的患者患有肺结节,在RL和EHOE之间的患者患有肺纤维化。在图10A中,3组柱形图(左、中、右)中的各组代表了每名患者的平均峰值正、负速度(单位是cm/sec)的测量结果,所述平均峰值正、负速度分别在每名患者的上、中、下区域获得的特征#1的平均峰值正、负速度。图10B和10C描绘了作为特征#3和5的相应的数据。应注意到,在图10A-C的每张图中,只有一名患者的数据包含有标签U、M和L(分别表示上、中、下部区),以避免杂乱。
对图10A-C描绘的数据进行检查揭示出,在正常患者中,特征#1的速度在所有的三个区中的大致类似。但是在COPD患者中,上部区的速度要远远低于中部区的速度,而中部区的速度要远远低于下部区的速度。对于特征#5而言,也确实具有同样的情形。因此,特征#1和5的那些相对速度的出现能用来预知COPD的出现。可以将区分正常患者与COPD患者的试验固定下来(例如,如果中部峰值速度的读数至少是上部峰值速度读数的两倍,且下部峰值速度的读数至少是上部峰值速度读数的三倍,那就指示是COPD)。可替换地,可以用以下将进行说明的参数化法来获得阈值水平。由此,我们会发现,可以将不同位置处特征的速度之间的差异用于辅助区分正常人和患不同疾病的患者。
图11A是正常人和COPD患者之间的差异的图解表示,其基础是,这两组患者的平均值,图11A强调了在上、中、下部区的特征#1和#5的峰值速度的区别。
可选地,上述数据可以与“功率声波图”数据组合,如2010年4月30日提出的美国申请12/771,091中所描述的,该项申请通过引用合并于此。然后,个人计算机13(如图1中所示)要被编程成能从超声回返中提取功率数据,如‘091申请中所述的。图11B表明的是,如此获取的功率数据能用来在正常人、患COPD的患者及罹患肺纤维化的患者之间进行区分。后者,能较好传导超声能量的结缔组织替代充满空气的肺泡,并由此会获得更高的总功率值。还应注意到,在肺纤维化的情形下(对比正常人和COPD情形),经常会在上部肺段中记录到最大的功率信号。可以用此来预知肺纤维化的存在。
充血性心力衰竭(CHF)、肺气肿及肺水肿的区别也可以通过他们不同的多普勒标记图被定性。例如,肺水肿患者中,功率会比正常人低,但是在CHF患者中,由于肺内过量的流体(会比正常患者平常是空气的情形下提供更少的信号衰减),功率可能会比正常人高。肺实质和脉管系统内的局部变化可以引起不同肺区之间的功率区别。可以用TPD检测到这些区别,并依靠这些区别来诊断那些情况,或者是从功率-和-速度与时间之比的显示内容中观察到,或者自动地使用适当的模式识别或参数化软件。可以将类似的概念用于其他病理。
对肺部血压的测量
很多情况及诸如CHF的肺部和心脏疾病都会导致肺部血压升高。尽管探测、定性、及对肺动脉高压(PH)的跟进也很重要,但所有的现有技术工艺都存在问题。有时候,用复合超声成像可以进行间接且不精确的估算。但是,唯一可信的测量方法是侵入式的——经由心脏将测量导管导入肺部血管中。相比之下,用TPD测量肺部血压是迅速、简易、有效且无需侵入的。
使用传统的血压计时,围绕外周血管(例如,上臂的,桡骨的)的压力会升高,而基于血管脉动中的变化来监测动脉脉搏及确定最高和最低压。在此框架内,通过血液流动和脉动停止处的压力来确定收缩压。如上文所解释的,TPD记录下的信号会反映出肺部血管内的脉动。这些血管被肺实质围绕,而肺实质由多个其压力可控的气室组成。因为这一点,就可以实现通过升高肺部气压和监测TPD信号来确定肺部血压,从而在不同压力情况下确定血流和经由肺内血管的血管脉动。
图12是施行此测量方法的系统的框图。在整个步骤中,TPD探针/传感器225要放置在患者的胸部226上,及肺部信号被TPD224进行记录、显示的处理。为了获取读数,将肺部气压升高之后回复正常。改变肺部气压的一种方式是让患者201使其肺膨胀到预定的程度,之后用力地向连接至储气器212的管200内吹气(例如,经由一次性接口管210)。这种情况下,主要由患者的横膈膜增加压力。压力会较佳地显示在显示器218(或压力计,未示出)上让患者看到,且患者会被告知用吹的行为使得压力保持在要求的水平。患者还会被告知要保持其声门打开,以使得整个系统内的压力均衡。如果使用这种方法,可以省略泵215和关联的硬件和软件。改变肺部气压的另一种方法是,用泵215来升高肺222内的压力,其中,泵215由控制器216和处理器217控制以驱使肺压达到期望的水平。当使用泵时,反馈较佳地通过压力传感器214获得。应注意到,在任何情形下,期望的压力水平都可以随着时间改变以遵循期望的曲线(例如,通过先增加压力,之后再通过渐进或分步的方法使压力缓慢降低)。
图13描绘了TPD信号如何响应于肺气压逐渐升高而进行变化,及在血管性质中产生的变化。当压力升高时,血管最终将会在外压等于或超过血压的节点上崩塌(完全或局部),这个节点在图13中发生在11和11.5秒之间(由箭头132指示)。这种现象类似于当压在上臂或桡骨动脉上的常规血压袖带施加的压力超过特定水平时,血流停止的方式。
图14描绘了当肺气压升高并维持在高水平时发生在患者身上的变化的实施例。右侧指示出的压力水平处的不同特征(上述的#1-5)在幅度和特性上的变化携带了信息,该信息是关于在相关血管中的多种血压水平的信息。应注意到,五个特征#1-5中的每个特征的变化可以发生在不同的压力下。例如,信号#1的正分量会在约16mm Hg的压力下消失,而负信号1*几乎完整无缺。在约10mm Hg的压力下,信号#3(3*)的负分量已经衰减,而正分量只会在更高的压力下衰减。在10mm Hg的压力下,信号#4也几乎消除了。
应注意到,正常的肺部血压(通过侵入方式将压力传感器导入相应血管测量所得)的肺动脉舒张压一般在10-15mm Hg的范围内,肺动脉收缩压力一般在25-30mm Hg的范围内,及肺部循环的静脉侧(肺静脉)的压力约为8-10mm Hg。但是由于这些值是对于压力转换器能导入的主要的大血管而言的,所以基于TPD的测量方法中的较低压力水平就很有意义,原因在于,较小血管中的压力水平很可能较低(尽管还未被证实)。因此可以将用TPD测得的压力与肺部循环的适当成分关联起来。
可以逐渐地升高肺气压,以在单个增压过程中记录特征#1-5在一个压力范围内的变化。图15给出了这种测量方法的实施例,其中压力缓慢增加,之后维持在约2kPa的高位平稳段,如中部面板154所描绘的。可以在顶部面板152的右侧区段看到,经由较小肺部血管的血流响应与压力减少地而进行回复。应注意到,此实施例中压力的升高涉及到肺膨胀,直至总肺容量达到3L,如通过肺活量测定法进行测定,如下侧面板156所描绘的。
内科医生可以基于不同的TPD特征#1-5何时萎缩或消失来作出对上述变化的解析及肺循环压的确定。可替换地,可以用适合的模式识别软件来自动检测相关的变化。
图16描绘了肺压升高时,TPD记录的信号的功率水平。可以通过点162来确定肺血管床压力,点162处功率幅度接近零(或降至最大值的10%以下)。图17A描绘了压力升至不同水平并在该水平维持相对较长时间(例如,10-20秒)的情形。信号如所述地衰减,并且,在压力水平与静脉循环的压力水平(描述的实施例中是12mm Hg)一致时,信号接近零。在新一轮的压力升高过程中,例如图17A中的升到15mm Hg,血流和脉动停止。然而,随着血流停止,沿回路的压力降归零,以使得整个系统逐渐达到高收缩压,且所有的血管再次膨胀,因此血流和脉动随着时间(由脉管系统的容量确定)再次出现。这可以从图17A中相应的测得的功率点中观察到。对这样的脉动将进行记录直到应用的和维持的压力值升高到等于或超过收缩动脉压的压力值。没有任何脉动的压力与肺动脉收缩压一致。由此,在两个点处曲线会接近零功率水平(或降至最大值的10%以下)。我们相信,曲线接近零功率水平(即,对此实验对象的压力读数为约12mm Hg)的第一点172相当于静脉侧的肺压。我们相信,曲线接近零功率水平(即,压力读数为约20mm Hg)的第二点174相当于动脉侧的肺压。
图17B和17C分别对正常人的总功率读数(图17B)和肺动脉高压患者的总功率读数(图17C)进行了比较。高血压患者很明显有更高的压力读数。可以通过在已知的时间间隔内,从功率和速度的数据(即,包含上述所有的特征#1-#5)中计算每个相关速度处的功率的总和来获取总功率。
由此,就可以通过用经胸的肺部多普勒超声连续获取的功率和速度的数据来实现对患者肺部血压的估算,其中功率和速度的数据来自至少一名患者的肺在多个不同气压水平中的每个气压水平下的数据。然后可以基于获得的数据估算患者的肺部血压。
可以通过监测总功率随着气压变化的情况来确定患者的肺部血压水平。可以通过提供不同种类血压的数值估算来确定此水平,如上所述。在可替换的具体实施方式中,可以提供对肺部血压水平的二进制指示,其中一个二进制状态指示出正常的肺部血压,另一个二进制状态指示出PHT(肺动脉高压),以下将要进行说明。
生成关于患者是否患有PHT的二进制指示的一种方式是,将患者至少一个肺内的气压升高至一个水平,该水平处预计健康患者(即,不患肺动脉高压的患者)的血流会发生下降。一旦这样做了,气压升高后,就能获得来自患者至少一个肺的功率和速度的数据。如果总功率(从获得的功率和速度数据计算得到)超过阈值(例如,到了预期的总功率的20%后观察患者肺内的气压是否没有升高),那么我们就可以预知该患者患有PHT。完成此试验需要的压力升高量是例如10、15或20mmHg。在更高的压力下,该试验会更加可靠。
关于患者是否患有PHT的二进制指示的另一种方式是,使用分类算法。这种方法依赖于从TPD获得的功率和速度的数据中对分类特征的抽取。这种分类特征的实例包含:在胸壁上的多个位置(例如,肋间隙("ICS")#2,#4和#6)及从表面挑选的多个距离处,对应于不同特征(例如,特征#1,#3和#5)的速度(峰值、平均值、中间值等)和功率整数值。
较佳的分类算法的一实施例使用了以下的4个分类特征:
A=在ICS4和ICS6之间的特征#3中的峰值速度的比率。
B=在ICS4和ICS2之间的特征#3中的功率整数值的比率。
C=在ICS2和ICS6之间的特征#1中的峰值速度的比率。
D=在ICS4和ICS6之间的特征#5中的功率整数值的比率。
将这四个特征均一化到[0-1]的范围,然后将其应用到Fisher线性判别,在该判别中会将选取出来的特征线性组合成一个判别特征。基于33个正常人和20个PHT患者的分类算法得出了以下方程式,用以指明患者是正常的或患有PHT:X=4.8499A+6.3762B-3.3423C-4.6710D。在此实施例中,最优化的判定阈值为0,如果X>0那么表明患者患有PHT,而如果X<0,则表明患者是正常人。Fisher线性判别在Ronald Fisher(1936)撰写的《The use of Multiple Measurements inTaxonomic Problems》(《分类学问题中的多测量方法的使用》)中进行了描述,该文刊于Annals of Eugenics(优生学的编年史)7,179-188页,该文通过引用合并于此。
较佳分类算法的另一实施例使用了与上述同样的4个分类特征A-D及带有径向基函数(RBF)的部分的支持向量机(SVM)。SVM在Chih-Chung Chang和Chih-Jen Lin.2011.LIBSVM中描述的:是支持向量机的数据库。SVM还在报刊上被描述,W.H.等人(2007)“Section16.5Support Vector Machines”Numerical Recipes:The Art of ScientificComputing(“段16.5,支持向量机”数值秘诀:科学计算的艺术)(第三版)纽约:剑桥大学出版社。这些文献都通过引用合并于此。
在五重交互效度分析中,将实验对象随机地划分成五个小组。在5个小组中,保留一个小组作为试验该模型的核实数据,剩下的4个小组为培训数据。然后将交互效度分析处理进行5次,将5个子样本中的每个样本都作为核实数据使用正好一次。最终结果取5次重复的平均值。使用上文标识的分类特征A-D,五重交互效度分析结果是90.5%(48/53)的真实分类。
在另一较佳实施例中,用下文的四个分类特征E-H来替代上文标识的分类特征A-D:
E=在ICS4中和ICS6中的特征#3的峰值速度之间的比率。
F=在ICS2中和ICS6中的特征#1的峰值速度之间的比率。
G=在ICS4中和ICS6中的特征#5的峰值速度之间的比率。
H=在ICS4中和ICS2中的特征#3的功率之间的比率。
可以在其他可替换实施例中使用的分类特征的实例包含:特征#1,#3,#5的持续时间;特征#1,#3,#5中的峰值速度;特征#1,#3,,#5中的峰值速度时间;特征#1,#3,#5中的功率数值;特征#1,#3,#5中的峰值功率整数;特征#1,#3,#5中的峰值功率时间;特征#1,#3,#中的峰值速度时间中的功率;特征#1,#3,#5中的峰值功率时间(正-负)之间的延迟;特征#1,#3,#5中的正负峰值功率数值之比;正/负特征#1,#3,#5中的速度和功率数值之间的相关性;特征#1,#3,#5中的正负峰值速度数值之间的相关性;特征#1,#3,#5中加权功率的峰值速度;特征#1和特征#5的上升斜率;及特征#1的下降斜率。还可以使用在不同的ICS中及距表面不同距离时的上述所有分类特征的线性或非线性组合。
自动特征识别
上文讨论中频繁提到了特征#1-5。可选地,能识别那些特征中各特征之间的线性的软件可被提供在个人计算机13(图1示出)中。自动特征识别("AFR")可以实施在与图6有关的上述经平均化的信号上,实施在单个信号(例如,如图2中所描述的)上,或在包含在NR内的平均化操作(即,上常规降噪的第一阶段)后被实施。图18是基于后者的自动特征识别的实施例。在图18中,特征#1-5中的每个特征被计算得出的包络上的两个局部最小点所划定界限,并根据其峰值速度(即,最大点)和经平均化的信号的ECG波形的相对位置被定义。图18中,这些局部最小点在不同的特征之间限定了过渡区181-185,并用虚线表示。在规则的心搏节律中,限定特征与ECG信号24的关系如下:#1--在第一R波26后的第一速度峰值段;#2--在特征#1后但在ECG的T波前的第一速度峰值段;#3--在T波结束后的第一速度峰值段;#4--被特征#3和特征#5所界定的段;#5—紧挨下一个R波和下个特征#1之前的速度峰值段。
AFR非常有用的原因在于,计算出来的绝对和相对参数可用于区分和诊断病理及其位置,其中绝对和相对参数是以这些段为特性的。这些参数对于依赖参数化法而对多种情况进行自动识别而言非常有益,讨论如下。
参数化法
参数化法可用于表征多种特征,以诊断及估算多种病理的程度,这多种病理诸如COPD、结节病、纤维性哮喘、肺气肿、肺动脉高压、肺栓塞、肿块、肺部血管的动脉硬化、肺不张、心脏收缩功能障碍及心律不齐等。多种参数的量化可以在具体段及多种参数的关系上完成,也可以在原始图谱中信号的变化性上完成(即,在它被平均化前)。参数化法可以用2010年2月5日提出的美国申请12/700,828("'828申请")中描述的方法来实施,该项申请通过引用合并于此。
其中的一些数据源自由多普勒测量提供的功率图谱本身。这些功率图谱的特征也可以被参数化,例如,具体速度处的功率、曲线的平均斜率、正负特征处的不同斜率的数量等。参数还可以源自速度和功率与时间之比的描记。以下表格包含可用于对TPD结果进行参数化的参数的实施例,及其定义:
Figure BDA00003072933800191
Figure BDA00003072933800192
Figure BDA00003072933800201
其他特征
年龄
体重
性别
身高
可以使用这些参数,用如'828申请中所描述地实施学习和将步骤分类。
总结
对于以下组织和结构的多普勒标记图会根据病理而发生变化:肺气肿、肺栓塞、肺动脉高压、肺血管狭窄和畸形、与肺纤维化相关的情况、肺炎、肺不张、气胸、充血性心力衰竭、肺部实体瘤、肺部血管中显现的多种心功能不全、肿块及异物等。由此,用TPD收集到的肺多普勒信号可以被用于提供关于肺实质和脉管系统的结构和完整性的见解和具有潜在价值的诊断信息。因此,TPD可以作为用于诊断肺疾病及肺功能的一种新的非侵入式且无损的工具。其也可以实现对衰竭的肺部或心血管系统的状态的连续监测,及有助于确定功效从而为最优化的治疗实现计量定标。
TPD的一项附增的独特诊断能力是,确定肺部血管树组分的顺从性(倒电容),该顺从性在动脉硬化及其他的血管情况下会变化。可以基于脉搏在血管中的扩散速度来测量血管的顺从性,原因在于,血管越僵硬,扩散得越快。就肺而言,扩散速度可以通过沿扩散路径的不同位置处的任何肺信号出现的时间(或其峰值等)之间的延迟来确定。在不同的肺位置或单个位置以下的不同深度处的不同记录内,可以通过手动或适当的软件作出这样的延迟测量。
虽然本发明在披露过程中提到了某些实施例、数值修改、变更,但也可能有不背离本发明所说明的领域和范围的实施例的变化,如附带的权利要求书中所限定的。据此,应预期到的是,本发明不局限于所述实施例,而具有由以下的权利要求书中的语言限定的整个范围,及其等同物。

Claims (28)

1.一种估算患者的肺部血压的方法,包括步骤:
用经胸肺部多普勒超声,在多个不同气压水平中的每一个下从所述患者的至少一个肺中连续地获取功率和速度的数据;及
基于在所述获取步骤中获取到的所述数据,估算所述患者的所述肺部血压。
2.根据权利要求1所述的方法,其中所述获取步骤包括从所述患者的所述至少一个肺内的至少两个不同位置处获取功率和速度的数据。
3.根据权利要求1所述的方法,其中所述获取步骤包括传送有效横截面至少为1/2cm2的脉冲波超声束。
4.根据权利要求3所述的方法,其中所述获取步骤包括传送脉冲重复频率在1到3kHz之间的脉冲波超声束。
5.根据权利要求1所述的方法,其中所述估算步骤包括对至少一个所述气压水平进行识别,在所述至少一个所述气压水平下的总功率接近零。
6.根据权利要求1所述的方法,其中所述估算步骤包括对至少一个所述气压水平进行识别,在所述至少一个所述气压水平下,总功率小于在所述患者的所述至少一个肺内的所述气压水平不升高时所获取的总功率的10%。
7.根据权利要求1所述的方法,其中所述获取步骤包括从所述患者的所述至少一个肺内的至少两个不同位置处获取功率和速度的数据,并且所述估算步骤包括将分类算法应用到在所述获取步骤中获取的所述数据。
8.根据权利要求1所述的方法,其中所述获取步骤包括从所述患者的所述至少一个肺内的ICS2、ICS4和ICS6处获取功率和速度的数据,并且所述估算步骤包括将分类算法应用到(a)在ICS4和ICS6之间的特征#3中的峰值速度的比率,(b)在ICS4和ICS2之间的特征#3中的功率整数值的比率,(c)在ICS2和ICS6之间的特征#1中的峰值速度的比率,比及(d)在ICS4和ICS6之间的特征#5中的功率整数值的比率,其中特征#1紧随着ECG的R波出现,且与收缩期的心室收缩相一致,特征#3与舒张期的快速充盈阶段相一致,及特征#5与心房收缩相一致。
9.一种用于测量患者的肺部血压的装置,所述装置包括:
压力传感器,所述压力传感器被配置成测量所述患者的至少一个肺内的气压;
转换器,所述转换器被配置成将超声能量传送到所述患者的至少一个肺内的目标区域中、检测从所述目标区域反射的超声能量、及生成基于所述检测到的超声能量的输出;及
多普勒信号处理器,所述多普勒信号处理器被配置成处理所述转换器的所述输出、在多个不同气压水平中的每一个下从所述患者的至少一个肺内连续获取功率和速度的数据、及基于所述获取到的功率和速度的数据对所述患者的所述肺部血压进行估算。
10.根据权利要求9所述的装置,其中所述功率和速度的数据是从所述患者的所述至少一个肺内的至少两个不同位置处获取。
11.根据权利要求9所述的装置,其中所述转换器传送有效横截面至少为1/2cm2的脉冲波超声束。
12.根据权利要求11所述的装置,其中所述转换器传送脉冲重复频率在1到3kHz之间的脉冲波超声束。
13.根据权利要求9所述的装置,其中所述处理器通过对至少一个所述气压水平进行识别来估算所述肺部血压,在所述至少一个所述气压水平下的总功率接近零。
14.根据权利要求9所述的装置,其中所述处理器通过对至少一个所述气压水平进行识别来估算所述肺部血压,在所述至少一个所述气压水平下的总功率小于当所述患者的所述至少一个肺内的所述气压水平不升高时所获取的总功率的10%。
15.根据权利要求9所述的装置,其中从所述患者的所述至少一个肺内的至少两个不同位置处获取所述功率和速度的数据,及所述处理器通过将分类算法应用到在所述获取步骤中所获取的所述数据来估算所述肺部血压。
16.根据权利要求9所述的装置,其中从所述患者的所述至少一个肺内的ICS2、ICS4和ICS6处获取所述功率和速度的数据,并且所述处理器通过应用分类算法来估算所述肺部血压,所述应用算法被应用至:(a)在ICS4和ICS6之间的特征#3中的峰值速度的比率,(b)在ICS4和ICS2之间的特征#3中的功率整数值的比率,(c)在ICS2和ICS6之间的特征#1中的峰值速度的比率,及(d)在ICS4和ICS6之间的特征#5中的功率整数值的比率,其中特征#1紧随着ECG的R波出现,且与收缩期的心室收缩相一致,特征#3与舒张期的快速充盈阶段相一致,及特征#5与心房收缩相一致。
17.一种确定患者的肺部血压水平的方法,所述方法包括步骤:
将超声能量传送进所述患者的至少一个肺内;
检测由所述至少一个肺内的血管和围绕所述血管且充满空气的肺泡之间的运动边界引起的反射超声能量的多普勒频移;
改变所述至少一个肺内的气压;
监测所述检测到的多普勒频移如何响应于所述压力的改变而发生变化;及
基于所述监测步骤中监测到的所述变化确定所述患者的肺部血压的水平。
18.根据权利要求17所述的方法,其中所述传送步骤包括传送有效横截面至少为1/2cm2的脉冲波超声束。
19.根据权利要求18所述的方法,其中所述传送步骤包括传送脉冲重复频率在1到3kHz之间的脉冲波超声束。
20.根据权利要求17所述的方法,其中所述确定步骤包括确定所述患者是否患有肺动脉高压。
21.根据权利要求17所述的方法,其中所述确定步骤包括确定所述患者的所述肺部血压的至少一个数值。
22.根据权利要求17所述的方法,其中所述传送步骤包括:
在相当于至少一个心动周期的时间段内,将超声能量传送到所述患者的至少一个肺内的第一位置;
在相当于至少一个心动周期的时间段内,接着将超声能量传送到所述患者的至少一个肺内的第二位置。
23.一种确定患者是否患肺动脉高压的方法,包括步骤:
使所述患者的至少一个肺内的气压升高到一个水平,在所述水平下,预计未患肺动脉高压的患者的血流量会下降;
在所述气压升高时,从所述患者的至少一个肺中获取至少一组功率和速度的数据;及
基于在所述获取步骤中获取到的所述功率和速度的数据,确定总功率是否高于阈值。
24.根据权利要求23所述的方法,进一步包括,在所述确定步骤中,当确定所述总功率高于所述阈值时,将所述患者分类为患有肺动脉高压。
25.根据权利要求23所述的方法,其中所述阈值是所述总功率的20%,预期在这种情形下正常患者肺内的气压不会升高。
26.根据权利要求25所述的方法,进一步包括将所述患者的至少一个肺内的所述气压升高至少10mm Hg。
27.根据权利要求25所述的方法,进一步包括将所述患者的至少一个肺内的所述气压升高至少15mm Hg。
28.根据权利要求25所述的方法,进一步包括将所述患者的至少一个肺内的所述气压升高至少20mm Hg。
CN201180050531.6A 2010-10-21 2011-10-19 使用经胸肺部多普勒超声测量肺部血压 Expired - Fee Related CN103298412B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US40545410P 2010-10-21 2010-10-21
US61/405,454 2010-10-21
PCT/IB2011/002493 WO2012052824A1 (en) 2010-10-21 2011-10-19 Measuring pulmonary blood pressure using transthoracic pulmonary doppler ultrasound

Publications (2)

Publication Number Publication Date
CN103298412A true CN103298412A (zh) 2013-09-11
CN103298412B CN103298412B (zh) 2017-05-31

Family

ID=44999811

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180050531.6A Expired - Fee Related CN103298412B (zh) 2010-10-21 2011-10-19 使用经胸肺部多普勒超声测量肺部血压

Country Status (6)

Country Link
US (1) US8968203B2 (zh)
EP (1) EP2629673A1 (zh)
JP (1) JP6148620B2 (zh)
CN (1) CN103298412B (zh)
CA (1) CA2836278A1 (zh)
WO (1) WO2012052824A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106175832A (zh) * 2016-06-27 2016-12-07 联想(北京)有限公司 一种检测血压的方法及移动终端
CN108926360A (zh) * 2018-05-30 2018-12-04 飞依诺科技(苏州)有限公司 一种查找目标频谱线的峰值速度点的方法和装置
CN114518470A (zh) * 2022-03-17 2022-05-20 国网河南省电力公司电力科学研究院 一种变压器内部绝缘油不均匀流速场声学成像检测方法

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012064769A2 (en) 2010-11-08 2012-05-18 Vasonova, Inc. Endovascular navigation system
CA2862511A1 (en) * 2012-01-26 2013-08-01 Yoram Palti Diagnosing lung disease using transthoracic pulmonary doppler ultrasound during lung vibration
US11020095B2 (en) 2015-01-14 2021-06-01 Echosense Jersey Limited Data compression to facilitate remote medical analysis and diagnosis
WO2016205891A1 (en) 2015-06-26 2016-12-29 Resmed Limited Diagnosis and monitoring of cardio-respiratory disorders
JP2019523119A (ja) 2016-08-03 2019-08-22 ピーアイ−ハーベスト ホールディング アクチェンゲゼルシャフト 血管内の血圧を含む身体内の圧力を非侵襲的に測定するシステムおよび方法
EP3524164A1 (en) 2016-08-03 2019-08-14 PI-Harvest Holding AG A system, method and software for non-invasive measurement of intravascular, in particular intracardiac blood pressure
US11877885B2 (en) 2017-07-31 2024-01-23 University of Pittsburgh—of the Commonwealth System of Higher Education System for automatic and early detection flap failure
CN113558659B (zh) * 2021-07-30 2023-07-04 重庆安酷科技有限公司 一种高精度超声波肺功能检测仪及其检测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006102489A (ja) * 2004-09-07 2006-04-20 Toshiba Corp 超音波ドプラ診断装置及び診断パラメータ計測方法
EP1870035A1 (en) * 2006-06-21 2007-12-26 IPRM Intellectual Property Rights Management AG Apparatus and computer program for determining a patient's volemic status represented by cardiopulmonary blood volume
WO2008138141A1 (en) * 2007-05-11 2008-11-20 Institut De Cardiologie De Montréal Labelled adrenomedullin derivatives and their use for imaging and therapy.
WO2009025667A1 (en) * 2007-08-21 2009-02-26 Massachusetts Institute Of Technology Determining cardiac output and left ventricular pressure
US20090131765A1 (en) * 2007-11-16 2009-05-21 Broncus Technologies, Inc. Method and system for measuring pulmonary artery circulation information

Family Cites Families (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3314991B2 (ja) * 1993-08-30 2002-08-19 株式会社日立メディコ 超音波ドプラ装置
US5546949A (en) * 1994-04-26 1996-08-20 Frazin; Leon Method and apparatus of logicalizing and determining orientation of an insertion end of a probe within a biotic structure
US5615682A (en) * 1995-10-26 1997-04-01 Hewlett-Packard Company Ultrasound transducer cable management system
WO1997045058A1 (en) * 1996-05-28 1997-12-04 Robin Medical Technologies, Ltd. Method and apparatus for cardiologic echo-doppler image enchancement by gated adaptive filtering in time domain
US6520911B1 (en) * 1996-07-03 2003-02-18 The United States Of America As Represented By The Department Of Health And Human Services Ultrasound-hall effect imaging system and method
US5836882A (en) * 1997-03-17 1998-11-17 Frazin; Leon J. Method and apparatus of localizing an insertion end of a probe within a biotic structure
US6210168B1 (en) * 1998-03-16 2001-04-03 Medsim Ltd. Doppler ultrasound simulator
EP1224436A1 (en) * 1999-10-29 2002-07-24 Universität Zürich Method of volumetric blood flow measurement
JP3532809B2 (ja) * 1999-12-03 2004-05-31 忠和 鄭 超音波診断装置
US7120489B2 (en) * 2000-05-08 2006-10-10 Brainsgate, Ltd. Method and apparatus for stimulating the sphenopalatine ganglion to modify properties of the BBB and cerebral blood flow
US7640062B2 (en) * 2000-05-08 2009-12-29 Brainsgate Ltd. Methods and systems for management of alzheimer's disease
US7117033B2 (en) * 2000-05-08 2006-10-03 Brainsgate, Ltd. Stimulation for acute conditions
US7146209B2 (en) * 2000-05-08 2006-12-05 Brainsgate, Ltd. Stimulation for treating eye pathologies
WO2002069805A1 (en) * 2001-03-02 2002-09-12 Palti Yoram Prof Method and apparatus for detecting arterial stenosis
US20020161304A1 (en) * 2001-04-30 2002-10-31 Eide Per Kristian Monitoring pressure in a body cavity
AU2002331850A1 (en) * 2001-09-11 2003-03-24 Pulmonx Methods of endobronchial diagnosis using imaging
US7186219B2 (en) * 2001-10-11 2007-03-06 Osypka Medical Gmbh Calibration of a doppler velocimeter for stroke volume determination
JP4245428B2 (ja) * 2003-07-29 2009-03-25 東芝医用システムエンジニアリング株式会社 超音波ドプラ診断装置及び画像データ生成方法
WO2005074809A1 (de) 2004-02-05 2005-08-18 Rueter Dirk Vorrichtung zur messung der funktion einer lunge
US20060058854A1 (en) * 2004-09-14 2006-03-16 Scimed Life Systems, Inc. Method for stimulating neural tissue in response to a sensed physiological event
ATE515978T1 (de) 2006-10-26 2011-07-15 Cardiogal Ltd Nicht-invasive herzparameter-messung
US8050476B2 (en) * 2006-11-22 2011-11-01 General Electric Company Heart rate demodulation of periodic movement in ultrasound data for object identification
JP4945300B2 (ja) * 2007-04-25 2012-06-06 株式会社東芝 超音波診断装置
DE102008040266A1 (de) * 2008-07-09 2010-01-14 Biotronik Crm Patent Ag Implantierbare Messanordnung
US20100274133A1 (en) 2009-02-05 2010-10-28 Yoram Palti Detecting a stenosis in a blood vessel
JP5642770B2 (ja) 2009-05-01 2014-12-17 パルティ、ヨーラム ドップラーに基づく流量測定
US20110054322A1 (en) * 2009-09-01 2011-03-03 Paolo Zanatta Microembolic signals detection during cordiopulmonary bypass
WO2011051787A2 (en) * 2009-10-27 2011-05-05 Palti Yoram Prof Transthoracic pulmonary doppler ultrasound

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006102489A (ja) * 2004-09-07 2006-04-20 Toshiba Corp 超音波ドプラ診断装置及び診断パラメータ計測方法
EP1870035A1 (en) * 2006-06-21 2007-12-26 IPRM Intellectual Property Rights Management AG Apparatus and computer program for determining a patient's volemic status represented by cardiopulmonary blood volume
WO2008138141A1 (en) * 2007-05-11 2008-11-20 Institut De Cardiologie De Montréal Labelled adrenomedullin derivatives and their use for imaging and therapy.
WO2009025667A1 (en) * 2007-08-21 2009-02-26 Massachusetts Institute Of Technology Determining cardiac output and left ventricular pressure
US20090131765A1 (en) * 2007-11-16 2009-05-21 Broncus Technologies, Inc. Method and system for measuring pulmonary artery circulation information

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106175832A (zh) * 2016-06-27 2016-12-07 联想(北京)有限公司 一种检测血压的方法及移动终端
CN106175832B (zh) * 2016-06-27 2019-07-26 联想(北京)有限公司 一种检测血压的方法及移动终端
CN108926360A (zh) * 2018-05-30 2018-12-04 飞依诺科技(苏州)有限公司 一种查找目标频谱线的峰值速度点的方法和装置
CN108926360B (zh) * 2018-05-30 2021-03-19 飞依诺科技(苏州)有限公司 一种查找目标频谱线的峰值速度点的方法和装置
CN114518470A (zh) * 2022-03-17 2022-05-20 国网河南省电力公司电力科学研究院 一种变压器内部绝缘油不均匀流速场声学成像检测方法
CN114518470B (zh) * 2022-03-17 2023-08-08 国网河南省电力公司电力科学研究院 一种变压器内部绝缘油不均匀流速场声学成像检测方法

Also Published As

Publication number Publication date
US20120101381A1 (en) 2012-04-26
JP6148620B2 (ja) 2017-06-14
EP2629673A1 (en) 2013-08-28
CN103298412B (zh) 2017-05-31
US8968203B2 (en) 2015-03-03
WO2012052824A1 (en) 2012-04-26
JP2013543740A (ja) 2013-12-09
CA2836278A1 (en) 2012-04-26

Similar Documents

Publication Publication Date Title
CN103298412A (zh) 使用经胸肺部多普勒超声测量肺部血压
CN103607956A (zh) 经胸的心肺监测器
Fukuda et al. Utility of right ventricular free wall speckle-tracking strain for evaluation of right ventricular performance in patients with pulmonary hypertension
Mor-Avi et al. Current and evolving echocardiographic techniques for the quantitative evaluation of cardiac mechanics: ASE/EAE consensus statement on methodology and indications endorsed by the Japanese Society of Echocardiography
WO2014024758A1 (ja) 医用画像診断装置および医用画像分析方法
JP6443517B2 (ja) 生体情報測定装置および生体情報測定装置もしくは生体情報測定システムの作動方法
JP6054743B2 (ja) 経胸壁肺ドップラー超音波
JP6382036B2 (ja) 超音波診断装置及び画像処理装置
US20120165674A1 (en) Ultrasound diagnosis apparatus, image processing apparatus, and image processing method
WO2003077765A1 (fr) Systeme d&#39;echographie
CN104254283A (zh) 在肺振动时用经胸肺多普勒超声诊断肺疾病
JPH06503483A (ja) 自動デジタルイメージ定量化システム
US11730381B1 (en) Systems, apparatuses, and methods for locating blood flow turbulence in the cardiovascular system
Jin et al. Evaluation of myocardial stiffness in cardiac amyloidosis using acoustic radiation force impulse and natural shear wave imaging
Trambaiolo et al. Assessment of left atrial appendage wall velocities by transesophageal tissue Doppler echocardiography: a clinical study in patients with sinus rhythm
CN114271850B (zh) 超声检测数据的处理方法及超声检测数据的处理装置
US20220370140A1 (en) Methods and systems for modeling a cardiac system
Taha et al. Pre-ejection period by radial artery tonometry supplements echo Doppler findings during biventricular pacemaker optimization
CN104323803B (zh) 基于平面波超声成像的声带振动成像与测量系统及方法
Tavoni Development of New Techniques for Clinical Applications of Jugular Venous Pulse with Ultrasound Devices
Raat Accuracy and precision of stroke volume and left ventricular ejection fraction quantification with handheld echocardiography
Andersen Contraction Fronts Of The Left Cardiac Ventricle: A Case For High Frame Rate Ultrasound
Mohajery et al. Automated 3D velocity estimation of natural mechanical wave propagation in the myocardium
Duncan Non-invasive cardiac imaging for the quantification of ventricular function: potential and future applications
Ge Real-time 3D echocardiography for congenital heart disease: From fetus to adults

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170531

Termination date: 20211019