CN103391749B - 通过超声波成像分析二尖瓣返流 - Google Patents
通过超声波成像分析二尖瓣返流 Download PDFInfo
- Publication number
- CN103391749B CN103391749B CN201180068167.6A CN201180068167A CN103391749B CN 103391749 B CN103391749 B CN 103391749B CN 201180068167 A CN201180068167 A CN 201180068167A CN 103391749 B CN103391749 B CN 103391749B
- Authority
- CN
- China
- Prior art keywords
- backflowing
- model
- aperture
- ultrasonic
- doppler
- 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.)
- Expired - Fee Related
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- 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/0883—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- 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/8979—Combined Doppler and pulse-echo imaging systems
-
- 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/8993—Three dimensional imaging systems
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Hematology (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Cardiology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
描述了一种超声诊断成像系统,其量化通过二尖瓣的返流。所述超声系统中的流动量化处理器(34)产生邻近返流孔口的流速场的数学模型。速度场模型产生指向所述返流孔口的速度矢量的值。针对超声系统的超声物理性质和操作的影响修改这些建模值以产生预期速度值。将所述预期速度值与所述超声系统测量的实际多普勒速度比较,并且将差异累加成均方误差,所述均方误差用于调节所述模型的参数,例如孔口位置和流速。当这种迭代处理利用期望的比较收敛时,从最终调节的模型导出的参数用于计算真实孔口位置、流率和体积流量。
Description
相关申请的交叉引用
本申请要求享有2010年12月23日提交的美国临时专利申请61/426669的权益。
技术领域
本发明涉及医学诊断超声系统,具体而言,涉及使用诊断超声成像系统分析二尖瓣返流。
背景技术
返流是一种严重的医学状况,其需要分析和适当的处置。刚好在左心室收缩以向身体中泵送血液之前,二尖瓣必须完全闭合,从而使得收缩将所有血流喷射到主动脉中。如果该瓣膜未完全闭合,左心室中的一些血液将通过不完全密封的瓣膜中的开口喷射回左心房。这种血液的回流,通常为通过不完全闭合的瓣膜小叶向回喷出的血流的小的瞬间的喷射,减小了从心脏的血液流出,并且因此降低了每次心脏收缩的效率。那么心脏必须要更快速地泵送,以便为身体供给其必要的滋养血流供给。由于其效率低,心脏过度工作,从而导致心力衰竭。
临床医生多年来已经使用超声成像来尝试检测返流血流。最初是通过查看心脏左侧的超声图像中血液的上述喷射来完成瓣膜返流的超声检测的。在过去二十年期间,通过二维(2D)彩色流动多普勒成像已经方便了对喷射的观察,在二维(2D)彩色流动多普勒成像中,通过仔细搜索泄漏心脏瓣膜附近这些异常局部流速来检测血液小喷射的高速和湍流。不过这样的图像平面的采集给这种主观方式提出了挑战:在所述图像平面中,喷射是最普遍存在的,并且喷射与心脏、瓣膜运动和二尖瓣附近的血流湍流以及喷射的瞬间发生耦合。近年来,在能够通过超声方式观察喷射位置的情况中,临床医生使用称为PISA(邻近等速表面积法的缩写)的技术,以尝试对返流血流进行量化。在这种方法中,通过彩色流动多普勒成像对嫌疑瓣膜和LV心室内部以及邻近瓣膜的区域成像。在发生喷射的时候,在邻近区域中形成流动会聚区域(FCR),因为该区域中的血流速度会瞬间向返流孔口加速。这种流动模式在彩色流动图像中导致混叠,因为流速瞬间超过用于彩色流动图像的速度范围。此刻捕获并在显示屏上定格彩色流动图像。之后测量FCR第一混叠线处的速度v,并测量从该混叠线到瓣膜孔口的假定中心的距离r。之后使用这两个测量结果,使用表达式Qt=2πr2v计算通过该孔口的流率。
当执行这种流程时会出现若干困难。一个困难是,当在彩色流动图像中捕获喷射时在喷射最大峰值处获得最大的准确度。心脏周期期间喷射的持续时间能够仅为300-450毫秒,然而,典型的彩色流动帧速率可能在10-20帧每秒的范围中。因此,有可能彩色流动图像帧之一的采集时间与当喷射恰好处于其峰值时的时刻不相同。临床医生能够针对额外的心动周期重复彩色流动采集序列,或者能够满足于由在恰好喷射峰值之外处进行测量而造成的不准确性。
另一个问题是,在彩色流动图像中不容易界定瓣膜孔口的中心。瓣膜组织产生大的超声反射,并且当进行扫描时快速运动,因而在图像中能够呈现为大块的模糊或不清楚的块体。因此,测量结果r的准确度能够因为不能估计孔口的精确位置而受到影响。
第三个问题是,基本PISA技术仅仅是单个的一维测量。仅进行一次速度测量,并且在计算中仅使用单个到孔口的半径r。该方法假设FCR中血流的其余部分与单个测量的部分表现相同。显然,进行单个测量中的任何不准确都将产生不准确的结果。
试图消除这些不准确性的基本PISA技术的扩展是在二维图像中描绘FCR外边界的弧周围进行多次速度测量。测量从每个速度测量点到孔口的距离r,并且使用多个测量结果计算流量Qt。尽管多个测量可以防止对单个不准确测量的依赖,但出现了另一个问题。一般通过如下方式完成单个测量(1D)技术:通过心脏顶端并直接通过假定的返流孔口瞄准超声探头的中心射束,并沿这条射束线测量速度v和距离r。在射束线因而与返流的这种矢量对准时,测量的速度将是准确的,并且不受非零多普勒角度的影响。如所知的,超声多普勒速度测量受到流动方向和超声射束方向之间角度的影响。与射束方向直接在一条线上的流动将得以准确测量,而与射束方向成非零角度的流动将由角度的余弦减小。与射束方向正交(90°)的流动将不会产生多普勒响应。因此,必须要根据流动矢量与射束方向之间的角度来校正超声测量的多普勒速度,以便产生更准确的速度测量结果。在2DPISA技术的情况中,由于在进行速度测量的每个点处的其流动矢量和射束方向之间的角度变化,沿FCR弧的额外速度测量将是不准确的。因此,尽管多个测量能够补偿单个速度测量中出现的误差,但合计的测量将由于每个测量点处多普勒角度的变化而不充分测量流率。
PISA技术中固有的又一个问题是,在瓣膜的返流孔口的精确位置的获知中的不准确。如上文提到的,由于血流向着孔口并通过孔口突然加速,所以在瓣膜泄漏附近,普遍存在彩色流动图像中的混叠。因此,彩色流动图像因孔口附近多普勒图像中混叠颜色的闪现而混乱。此外,在心脏收缩期间二尖瓣平面是运动的。由于PISA技术的准确度依赖于为了测量距离r的对孔口位置的获知,对于精确获知孔口位置的这些障碍能够导致不准确的r测量,并因此导致流率的不准确计算。实验证明,孔口位置小到一个毫米的误差都能够导致在计算通过漏洞的流率中的显著误差。此外,已知许多返流瓣膜没有单个的针孔漏洞,而是沿不精确瓣膜闭合的缝隙的漏洞。因此,假设漏洞是由通过瓣膜的单个针孔造成可能不是有效的。
发明内容
因此,期望提供一种用于量化二尖瓣返流的流率和体积流量的方法和装置,其克服了PISA技术的以上限制。本发明的目的是准确识别瓣膜返流的孔口的精确位置。本发明的另一目的是在每个测量位置进行不受多普勒角度影响的一系列测量。本发明的另一目的是能够识别和量化不仅由于单个针孔漏洞,而且由于沿瓣膜闭合的缝隙导致的二尖瓣返流。
根据本发明的原理,描述了一种用于量化返流血流的超声诊断成像系统和方法。在返流瓣膜的超声图像中二尖瓣漏洞的假设位置附近描绘了血流的弧形(二维)或半球(三维)区域。所述区域的内边界处在与漏洞的孔口相邻的混叠区域处或刚好超过与漏洞的孔口相邻的混叠区域,并且所述区域的外边界从所述内边界向外并且与所述内边界同心。对返流做出初始假设,包括假设返流孔口的位置。使用这些假设来计算返流的速度场模型。根据超声系统的超声物理性质和/或设置来调节模型的参数,以确定在所描绘区域中的待由超声系统观察的预期速度矢量,VOBS。由超声系统在所描绘区域之内进行速度测量,并且将测量的速度与预期的速度相比较。计算测量值与预期值之间的差异,并通过非线性曲线拟合,对模型值中的一个或多个进行调节。迭代地重复所述处理,直到预期值和测量值能接受地收敛。能够使用经迭代调节的模型参数来产生流率或瓣膜中泄露孔口位置的量化测量。
根据本发明的另一方面,通过以上过程自动更新所述超声图像中所述返流孔口的初始假设位置,以在所述图像中指示所述孔口的实际位置。因此,尽管在二尖瓣附近存在图像混乱,也向临床医生显示出所述返流孔口的精确位置。
根据本发明的又一方面,在沿着二尖瓣的多个空间位置处重复以上技术。在所述空间位置处测量的合计返流提供了由是二尖瓣闭合中的缝隙而不是单个的泄露针孔口的返流漏洞造成的返流的测量。
根据本发明的另一方面,描述了若干壁滤波器特性,所述壁滤波器特性对于根据本发明的返流测量的实施方式是优选的,其包括这样的壁滤波器:所述壁滤波器具有在中等流速下到达峰值的特性,中等流速常常在邻近返流孔口的流动区域中发现。
附图说明
在附图中:
图1以框图形式图示了根据本发明的原理构建的超声诊断成像系统。
图2为图示了图1的流动量化处理器的功能的框图。
图3图示了返流喷射的超声多普勒图像。
图4图示了根据本发明的超声图像,所述超声图像具有一开始假设的孔口位置和自动指示的孔口位置。
图5图示了根据本发明的用于量化沿瓣膜的缝隙的返流的技术。
图6图示了适用于本发明的超声系统中的壁滤波器的响应特性。
具体实施方式
首先参考图1,以框图形式示出了根据本发明的原理构建的超声诊断成像系统。在图1中,在超声探头10中提供了换能器阵列10',以用于发射超声波并接收回波信息。换能器阵列10'优选为能够在三维中,例如在关于二尖瓣位置的高度和方位角中,进行扫描的换能器元件的二维阵列,以用于进行3D成像。所述换能器阵列耦合到所述探头中的微射束形成器12,微射束形成器12控制阵列元件的信号发射和接收。如美国专利5997479(Savord等人)、6013032(Savord)和6623432(Powers等人)中描述的,微射束形成器能够对换能器元件的组或“片”接收的信号进行至少部分射束形成。由探头线缆将所述微射束形成器耦合到发射/接收(T/R)开关16,发射/接收(T/R)开关16在发射和接收之间切换,并保护主射束形成器20以避免高能发射信号。由耦合到所述T/R开关和射束形成器20的发射控制器18指导在微射束形成器12控制下的来自换能器阵列10的超声射束的发射,发射控制器18从用户对用户接口或控制面板38的操作接收输入。由所述发射控制器控制的功能之一是引导射束的方向。可以从(正交于)换能器阵列的正前方,或在针对更宽视场的不同角度下引导射束。
将微射束形成器12产生的部分射束形成的信号耦合到主射束形成器20,在那里将来自元件各个片的部分射束形成的信号被组合成完全射束形成的信号。例如,主射束形成器20可以具有128个通道,其每个都从12个换能器元件的片接收部分射束形成的信号。以这种方式,由二维阵列的超过1500个换能器元件接收的信号能够有效贡献于单个射束形成的信号。
将射束形成的信号耦合到信号处理器22。信号处理器22能够通过多种方式处理所接收的回波信号,所述多种方式诸如为带通滤波、抽选、I和Q分量分离以及谐波信号分离,谐波信号分离用于分离线性和非线性信号,以便能够识别从组织和微泡返回的非线性回波信号。所述信号处理器还可以执行额外的信号增强,诸如斑点减少、信号混合和噪声消除。
将所处理的信号耦合到B模式处理器26和多普勒处理器28。B模式处理器26采用幅度检测以对身体中的结构成像,所述结构诸如是心壁组织、二尖瓣以及血细胞。可以以谐波模式或基本模式,或者两者的组合,来形成身体结构的B模式图像,如美国专利6283919(Roundhill等人)和美国专利6458083(Jago等人)中所描述的。多普勒处理器28处理来自组织和血流的时间上有区别的信号,以用于检测图像场中物质的运动,诸如血细胞的流动。多普勒处理器通常包括壁滤波器,所述壁滤波器具有这样的参数:可以设置所述参数以通过和/或拒绝从身体中选定类型的材料返回的回波。例如,能够将所述壁滤波器设置为具有通带特性,所述通带特性通过来自较高速度材料的相对低幅度的信号,而拒绝来自较低或零速度材料的相对强的信号。这种通带特性将通过来自流动血液的信号,而拒绝来自附近静止或缓慢运动对象的信号,诸如来自心壁的信号。反转特性将通过来自心脏运动组织的信号,而拒绝血流信号,以用于所谓的组织运动的组织多普勒成像、检测和描绘。所述多普勒处理器接收并处理来自图像场中不同点的时间离散回波信号的序列,来自特定点的回波序列称为系综(ensemble)。能够利用多普勒频率与指示血流速度的速度的对应性,使用以在相对短间隔上的快速相继方式接收的回波系综来估计流动血液的多普勒频移。使用在更长时段上接收的回波系综来估计更慢流动血液或缓慢运动组织的速度。针对迅速发生喷射的二尖瓣返流评估,一般采用短系综长度(更少样本),从而能够实现高的采集帧速率。可以通过形式的方程估计多普勒频移Δf,其中f0为发射频率,c是超声传播速度,v是速度,并且θ是射束方向和血流方向之间的角度。
将由B模式和多普勒处理器产生的结构和运动信号耦合到扫描转换器32和多平面重新格式化器44。所述扫描转换器以这样的空间关系布置回波信号:从所述空间关系,以期望的图像格式接收所述回波信号。例如,所述扫描转换器可以将回波信号布置成二维(2D)扇形格式或锥形三维(3D)图像。所述扫描转换器能够在B模式结构图像上叠加颜色以产生彩色多普勒图像,所述颜色对应于图像场中在与它们多普勒估计速度对应的点处的运动,所述彩色多普勒图像描绘了图像场中组织和血流的运动。所述多平面重新格式化器将从身体的体积区域中的公共平面中的点接收的回波转换成该平面的超声图像,如美国专利6443896(Detmer)中所描述的。体积绘制器42将3D数据集的回波信号转换成如从给定参考点处观看的投影3D图像,如美国专利6530885(Entrekin等人)中所描述的。从扫描转换器32、多平面重新格式化器44和体积绘制器44将2D或3D图像耦合到图像处理器30以用于进一步增强、缓冲和暂时存储,以在图像显示器40上显示。
根据本发明的原理,将由多普勒处理器28产生的血流速度值耦合到流动量化处理器34。所述流动量化处理器如下所述运行,以产生通过返流孔口的流率、通过所述孔口的体积流量以及所述孔口的空间位置的测量。所述流动量化处理器可以从用户控制面板38接收输入,所述输入例如为如下描述的孔口位置的初始估计。将来自所述流动量化处理器的输出数据耦合到图形处理器36,以用于利用显示器40上的图像再现来自所述处理器的输出数据。图形处理器36还能够生成图形叠加,以用于和超声图像一起显示。这些图形叠加能够包含标准识别信息,例如患者姓名、图像的日期与时间、成像参数等。出于这些目的,所述图形处理器从用户接口38接收输入,所述输入例如为键入的患者姓名。所述用户接口还耦合到发射控制器18以控制来自换能器阵列10'的超声信号生成,并因此控制由换能器阵列和超声系统产生的图像。所述用户接口还耦合到多平面重新格式化器44,以用于根据如下描述的本发明,选择并控制多个多平面重新格式化的(MPR)图像的显示,所述显示可以用于量化MPR图像中的返流。
图2描述了图1的流动量化处理器的运行。所述处理器基于如下描述的在包含区域之内的采样点处的流速场的数学模型。在2D图像实施例中,所述包含区域优选为由以返流孔口位置{x0,y0,z0}为中心的两个弧描绘的弧形条带。在3D图像实施例中,所述弧形条带为半球形壳体的形式,如下文结合图3描述的。在优选实施例中,该模型为如下形式的流速场的矢量速度模型:
或
其中针对流动矢量 并且r是从所述包含区域中三维坐标为(x,y,z)的点到返流孔口的距离,并且所述孔口位于三维空间中的坐标{x0,y0,z0}处。这意味着,当用矢量方式表达以指示到达所述孔口的方向时,距离r的形式为
在图2的方框50中,以这种方式对所述包含区域中每个点的流速进行建模。一开始,该模型能够以针对未知参数的假定的或估计的值开始,所述未知参数包括流量和返流孔口的位置{x0,y0,z0}。例如,该模型能够以作为返流的公称特性的参数开始。用户能够输入参数,例如通过在超声图像中指示返流孔口的假设位置。或者能够使用诸如PISA的已知技术来计算用作针对流速场的起始参数的值。
流速场模型对速度矢量进行近似,如果已知真实的生理速度矢量,这会是准确的。在方框50的输出处将由该模型近似的速度矢量指示为~V真实。方框52之后由于诸如所用超声系统的超声物理性质和操作参数的实际因素,对~V真实施加一些限制和调节。考虑到由超声系统观察到的实际速度值,之后调节或缩放流速场模型以将这些实际因素考虑进来。能够对其进行调节的实际因素之一是多普勒角度。如上文提到的,仅当流动方向与超声射束方向一致时,即多普勒角度为零时,超声系统执行的多普勒测量才是精确准确的。针对射束方向和流动方向之间的所有其他角度,对速度测量不充分。在上文提出的多普勒频移方程中,在结果中由cosθ项对多普勒角度加权,其中θ为多普勒角度。超声系统的另一个实际因素是壁滤波器偏置。多普勒壁滤波器通常呈现出非线性特性,这对于血流检测而言,将在DC(无运动)处具有零响应,并在选定的±f高频处上升到最大响应。或者,如下文讨论的,所述壁滤波器能够在低于由系综采样率的尼奎斯特极限指定的频率处呈现出最大响应。在其处测量多普勒频移的身体中的样本体积将不是身体中的单个点,而是将具有有限的尺寸,从而导致多普勒信号的返回指示速度的范围。所述壁滤波器的非均匀响应能够导致宽范围的速度以在响应特性的不同位置经历更大的增益,例如,更大幅度地强调较高速度,而非较低速度。这种非均匀响应能够产生速度分布的感觉中心的偏移,这称为壁滤波器偏置。也能够通过调节模型考虑壁滤波器偏置效应。能够考虑的另一相关因素是谱分布,即由从样本体积到阵列换能器的活动孔径的每个接收元件的不同路径和角度导致的多普勒谱变宽效应。参见美国专利5606972(Routh)。能够考虑的又一因素是混叠效应,即当血流运动速度超过能够由回波系综的采样率的尼奎斯特极限清晰检测的速度时,多普勒频率和速度的误报。在方框52中,通过对流速场模型进行点积调节而考虑诸如这些因素的因素,针对前面描述的模型,所述流速场模型可以为如下形式:
其中是正在针对其进行调节的物理和超声系统因素的矢量表达。图2中所示的结果是VOBS,即由预期效应以及在超声系统测量的信号中的预期内容所调节的模型。
比较器54将来自流速场模型的预期速度值VOBS与多普勒处理器28产生的来自场中的点(采样体积)的实际速度测量结果VMEAS进行比较。在方框56中,将预期值和接收值之间的差求平方,以产生针对每个点的误差项。在整个包含区域上对所述误差项积分以从整个区域产生均方误差项,所述整个包含区域可以是一维(例如,1D线)、二维(例如,2D线弧形区域)或三维(例如,3D半球形壳体)包含区域。之后使用所述误差项调节场模型的参数,例如r和流率,以导致后续测量的误差项的减小。优选的调节技术是使用非线性曲线拟合以向着误差减小的方向调整所述模型。可以使用的一种这样的非线性曲线拟合技术是Levengerg-Marquardt算法,该算法细化返流孔口位置的坐标以及朝向或通过所述孔口的流量或流率(Qt)。
迭代地重复图2的循环以减小所述误差项。每次通过循环,更精确调整的场模型就得以调节并与来自接受区域的超声速度测量结果进行比较,从而迭代地减小所述误差项。当已经将所述误差项减小到能接受的低值时,该模型呈现出孔口位置和通过孔口的血流的期望实际测量结果。这些项以及如期望的其他项产生为流动量化处理器34的输出并在显示器上将其呈现给用户。
图3是概略的2D超声图像,其图示了本发明的实施以测量二尖瓣返流。线100表示存在通过孔口O的返流泄露的二尖瓣平面。当左心室收缩时,血液喷射102逃逸回右心房。在PISA方法中,将由延伸通过孔口O的多普勒射束110探寻这种喷射。美国专利6719697(Li)提出了对PISA技术的改进,其中从多普勒射束110的位置在心脏周期上产生彩色M模式显示。设置彩色多普勒显示处理器以在图像右侧显示由彩色条120指示的极限,彩色条120以相对方向从零速度中心参考(BK=黑色)延伸到+V和-V的最大速度。根据速度的这个范围显示颜色的范围(Y=黄色;DR-暗红;等)。当测量的速度超过+V和-V的尼奎斯特极限时,产生混叠。与孔口紧邻的是流动会聚区域(FCR)104。在流动会聚区域之外的是由本发明的实施方式使用的接受区域112(S)。
PISA技术和本发明两者都以如下假设为前提:在左心室中孔口附近返流将会聚并向孔口位置O流动。这由接受区域S中的流动矢量和指示。但如多普勒射束110的方向图示的,即使当使用相控阵列扇形探头的变化射束角时,实际上并非所有流动矢量与射束方向对准。因此,针对不同的流动矢量将有不同大小的多普勒角度,如上文描述的,在图2的方框52中对流速场模型进行多普勒角度调节考虑了这一点。
FCR104的彩色多普勒图像将是混乱而不稳定的。这是因为,在这个区域中,血细胞的空间取向和血流速度非常迅速地改变,因为在心脏收缩期间血流会改变方向并瞬间向孔口O加速。在这个间隔期间采集的系综样本常常彼此不相关,这消除了通常多普勒估计处理器所依赖的系综之内的相关性。结果,即使认为FCR之内的流动可以是大致层流的,但彩色流动显示能够是强烈湍流和混叠。优选将接受区域S的内边界108设置到流动会聚区域104之外,以避免利用来自FCR的不稳定速度估算。设置边界108的一种方式是将其设置在速度剪切边界处或超过速度剪切边界。能够参考彩色流动显示完成这一点,或者参考速度剪切阈值自动完成这一点。
可以相对于能够测量的速度来设置接受区域的外边界106。尽管返流的速度在孔口附近相对高,但离所述孔口距离越大,它们逐步变得越低。能够将外边界106设置在距孔口O这样的距离处:在该距离处,仍然能够可靠地测量低多普勒速度。可以相对于尼奎斯特极限的百分比或相对于低流速,例如5mm/sec,确定这一距离。因此能够将外边界106设置在这样的距离:在该距离处,多普勒处理器仍然能够实现对低流动的能接受灵敏度。
在图3范例中的接受区域112被看作具有在孔口O处的旋转中心的二维弧形区域S。可以在一维、二维或三维中使用本发明。一维实施方式和模型可以仅考虑例如弧108和106之间的射束线110的线段。二维实施方式会是对诸如图3中接受区域112的平面接受区域采样的实施方式。三维实施会考虑弧形接受,例如112,但在以所述孔口为中心的整个半球中。对于二维和三维实施方式,预计会有更高准确度。优选的实施方式将使用如图1所示的具有二维阵列换能器的3D成像探头,其中在孔口周围的半球或四分之一球体积壳体中进行采集和建模。能够利用二维成像实现更高帧速率,在这种情况中,可以利用多平面重新格式化器44,从3D探头扫描的体积选择通过所述孔口和喷射的MPR帧。能够以高采集帧速率重复扫描期望的平面,并且利用2D MPR图像中所示的二维接受区域完成速度测量和流动计算。
尽管如图3所示能够存在从那里流动以环流模式向所述孔口会聚的接受区域,但已经发现在一些情况中存在其他非环流模式。换言之,指向所述孔口的流动矢量可以被布置为纯圆弧以外的模式。作为图3圆形的备选,接受区域可以具有抛物线或扁平圆形的形状。在三维实施方式中,接受区域可以是抛物面或扁球壳体。接受区域,即在其中与模型相比较地测量速度的区域或体积,其形状可以在迭代通过图2的处理循环期间动态改变,从而该过程将调整并会聚到接受区域,所述接受区域具有从用于初始模型中的形状改变的形状。
如之前提到的,临床医生常常难以准确定位超声图像上微小返流孔口的位置。当心脏不收缩时,所述孔口不存在,并且在其存在时,其位置被所得湍流,诸如在流动会聚区域附近的那些湍流的成像伪影破坏。根据本发明的另一方面,本发明的实施方式能够自动在超声图像上指示泄漏孔口的位置。图4图示了超声图像,其中前面结合图3描述的元件具有相同的附图标记。这幅图像还包括小框130,小框130由临床医生在开始诊断时放置在图像上,以指示返流孔口的假设位置。由临床医生操纵对诸如跟踪球的控制面板38的控制,以在图像中临床医生认为所述孔口所在的位置放置图标130。使用指示的孔口图标130的坐标引发图2中所示的流动量化处理器的处理,由此将所述坐标用作方框50的数学模型中的所述孔口的初始坐标。处理器34响应于测量速度值VMEAS的接收迭代地细化指向孔口位置的流速矢量的建模值。已经发现孔口位置的垂直放置对模型速度矢量与测量速度值的收敛具有最大影响。在心缩期期间心脏收缩同时,收缩运动导致二尖瓣在图像中向着左心室中心垂直运动。因此,心脏的自然收缩运动能够是返流评估中最大误差的来源。处理器34迭代地细化孔口位置,以减小边界106与108之间的流速场(接受区域S)中测量速度值与估计速度值之间的误差不一致。在流动量化处理器迭代地收敛于所述孔口的真实坐标时,最后确定的坐标由流动量化处理器34和图形处理器36用于自动将孔口图标130移动到图像中计算出的坐标。或者,所述图形处理器将在显示器上在迭代模型调节确定的真实坐标位置处放置另一个(计算出的)孔口图标132。如之前提到的,小到1mm的所述孔口的错误放置就能够显著影响孔口流率和体积流量的计算,体积流量是流率在这样的时间上的积分:在所述时间期间,心脏收缩(大约为心脏周期的1/3)。流量体积将在收缩间隔的大致中间处到达峰值。在图4中可以看出,该系统已经将计算出的孔口图标132放置在图像场中其真实位置处,这与临床医生对孔口位置的初始估计不同。当利用2D成像执行超声研究时,孔口坐标将一般在(x,y,z)笛卡尔坐标中。当使用3D成像时,一般将采用球坐标系。此外,如图4的超声图像图示的,在显示的超声图像中图形描绘了接受区域S,并且在显示器上显示的接受区域的显示形状和图形描绘中并入了对所述接受区域形状的任何修改。
如上文提到的,返流孔口并非始终作为闭合瓣膜中的单一针孔口漏洞而存在。可以由缝隙开口或多个小漏洞产生漏洞。图5是使用本发明实施方式评估瓣膜闭合中的缝隙和多个漏洞的流率和体积流量的图示。如图5图示的,针对沿泄露二尖瓣100'排列的一系列孔口位置执行本发明的过程。能够使用这些离散针孔口位置对闭合瓣膜中的复杂的缝隙状孔口或多个孔口建模。图5示出了三种这种确定的范例,每种确定都具有其自身的接受区域和流动矢量(由小箭头指示),所述接受区域由外部区域边界106、106'和106''指示,所述流动矢量指向沿瓣膜100'的不同孔口位置的会聚。以矢量方式将与每个针孔口相关联的矢量速度场与其他孔口相结合,从而与单个结合的矢量速度场一起定义了单个结合的接受区域。针对多个孔口,结合的矢量速度场现在近似于真实的生理速度矢量(~V真实,图2中方框50的输出)。结合矢量速度场的处理确定了与沿二尖瓣闭合的整个缝隙对应的所有针孔口的流率或体积流量。尽管图5给出了从三个针孔口生成的矢量速度场范例,但可以使用任意数量的孔口对缝隙进行充分建模。可以如图5中所示的,结合用于每次测量的接受区域,因为针对结合孔口确定的流动是指向沿缝隙的唯一孔口或位置的矢量流动。在三维实施方式中,测量线不限于平面中的直线,而可以跟随二尖瓣小叶闭合的非线性路径。
图6图示了可以在本发明的实施方式中采用的两个多普勒壁滤波器的响应特性。图6的制图的横坐标单位是采样数据壁滤波器的尼奎斯特极限,其中+1和-1的极限是对应于血流速度的滤波器的归一化尼奎斯特极限。零(DC或无流量)状况处于中心。纵坐标指示滤波器特性的相对幅度响应。响应曲线140是典型的用于测量血流的壁滤波器响应曲线,并且其可以用于本发明的实施方式中。这种特性在中心具有零响应,从而导致对诸如静止血管壁的静止对象没有响应。可以看出,响应随着曲线从零向外延伸而逐渐增大,最大响应在尼奎斯特极限处,其针对在选定采样间隔处无混叠情况下能检测的最高速度流动。结果,这种响应特性对于最高流速下的流动最灵敏。
响应曲线150是针对本发明实施方式的壁滤波器通常优选的响应曲线。可以看出,这条曲线150的响应在相对高但中等的采样率下达到峰值,刚好高于±0.5尼奎斯特,优选在1/2到2/3尼奎斯特的范围中,并在尼奎斯特极限处下降到零响应。这种滤波器特性被设计为对较低流速更灵敏,在接受区域S的外边界106附近能够预期有较低流速。这有助于补偿对低速度流动的低灵敏度,这源自于通常用于返流彩色流动成像的短系综。能够通过曲线150对图中心附近的较低速度的更大响应看出这一点。尽管可以看出曲线150在尼奎斯特极限处下降到零,但当与对于低速度流动的更大灵敏度交换时,高速度灵敏度的这种损失一般是能接受的。
Claims (17)
1.一种用于测量返流的诊断超声系统,包括:
具有换能器阵列的超声探头,其用于向身体中返流的位置发射超声能量并从所述身体中所述返流的位置接收超声回波;
图像处理器,其响应于所接收的回波,用于产生所述返流的位置的超声图像;
多普勒处理器,其响应于所接收的回波,用于产生所述返流的位置附近的血流速度的多普勒超声测量结果;
流动量化处理器,其被配置为基于包括对所述返流的位置的假设的关于所述返流的初始假设产生所述返流的位置附近的流速场的模型,将从所述模型导出的近似流速值与所述多普勒处理器产生的速度测量结果进行比较,并响应于所述比较调节所述模型;以及
显示装置,其耦合到所述图像处理器和所述流动量化处理器,用于显示所述返流的位置的所述超声图像以及从所述流速场的经调节的模型导出的流量测量结果。
2.根据权利要求1所述的诊断超声系统,其中,所述流动量化处理器还被配置为针对当使用超声系统测量速度时遇到的因素,来调节近似流速值。
3.根据权利要求2所述的诊断超声系统,其中,所述因素包括由于用于测量速度的所述超声系统的超声物理性质和操作参数而产生的影响中的至少一个。
4.根据权利要求3所述的诊断超声系统,其中,所述因素之一包括多普勒角度。
5.根据权利要求3所述的诊断超声系统,其中,所述因素之一包括壁滤波器偏置。
6.根据权利要求3所述的诊断超声系统,其中,所述因素之一包括多普勒频谱加宽。
7.根据权利要求3所述的诊断超声系统,其中,所述因素之一包括混叠效应。
8.根据权利要求1所述的诊断超声系统,其中,所述流速场包括孔口附近的接受区域,所述接受区域呈现二维弧形形状。
9.根据权利要求1所述的诊断超声系统,其中,所述流速场包括孔口附近的接受区域,所述接受区域呈现三维半球形状。
10.根据权利要求1所述的诊断超声系统,其中,响应于所述比较调节所述模型包括调节从所述模型导出的近似流速值。
11.根据权利要求1所述的诊断超声系统,其中,响应于所述比较调节所述模型包括调节所述流速场的形状。
12.根据权利要求11所述的诊断超声系统,其中,经调节的所述流速场的形状包括抛物线形状、扁平圆形状、抛物面或扁球体中的一种。
13.根据权利要求1所述的诊断超声系统,其中,响应于所述比较调节所述模型包括调节返流孔口的位置或者到返流孔口的距离。
14.一种用于利用超声系统评估返流的方法,包括:
利用超声探头和图像处理器对身体中返流的位置进行成像;
基于包括对所述返流的位置的假设的关于所述返流的初始假设产生所述返流的位置附近的流速场的模型;
产生所述返流的位置附近的血流速度的多普勒超声测量结果;
将从所述模型导出的近似流速值与所述血流速度的多普勒超声测量结果进行比较;
响应于所述比较调节所述模型的一个或多个参数;并且
显示所述返流的位置的超声图像以及从经调节的模型的参数导出的流动测量结果。
15.根据权利要求14所述的方法,还包括在所述比较的步骤之前,针对所述超声系统的超声物理性质或操作参数中一个或多个的影响,调节所述流速场的所述模型。
16.根据权利要求14所述的方法,其中,产生还包括产生位于返流孔口附近的弧形流速场的模型。
17.根据权利要求14所述的方法,其中,产生还包括产生位于返流孔口附近的半球流速场的模型。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201061426669P | 2010-12-23 | 2010-12-23 | |
US61/426,669 | 2010-12-23 | ||
US201161466046P | 2011-03-22 | 2011-03-22 | |
US61/466,046 | 2011-03-22 | ||
PCT/IB2011/055739 WO2012085797A1 (en) | 2010-12-23 | 2011-12-16 | Analysis of mitral regurgitation by ultrasonic imaging |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103391749A CN103391749A (zh) | 2013-11-13 |
CN103391749B true CN103391749B (zh) | 2015-08-19 |
Family
ID=45524894
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180068167.6A Expired - Fee Related CN103391749B (zh) | 2010-12-23 | 2011-12-16 | 通过超声波成像分析二尖瓣返流 |
Country Status (7)
Country | Link |
---|---|
US (1) | US10674993B2 (zh) |
EP (1) | EP2654570B1 (zh) |
JP (1) | JP6060089B2 (zh) |
CN (1) | CN103391749B (zh) |
BR (1) | BR112013015817A2 (zh) |
RU (1) | RU2569695C2 (zh) |
WO (1) | WO2012085797A1 (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4119098A1 (en) | 2009-04-15 | 2023-01-18 | Edwards Lifesciences CardiAQ LLC | Vascular implant and delivery system |
WO2013116807A1 (en) * | 2012-02-03 | 2013-08-08 | Los Alamos National Security, Llc | Systems and methods for synthetic aperture ultrasound tomography |
US9681951B2 (en) | 2013-03-14 | 2017-06-20 | Edwards Lifesciences Cardiaq Llc | Prosthesis with outer skirt and anchors |
US20150342569A1 (en) * | 2014-05-30 | 2015-12-03 | Siemens Medical Solutions Usa, Inc. | Transparency control for medical diagnostic ultrasound flow imaging |
US11055899B2 (en) * | 2016-06-10 | 2021-07-06 | Koninklijke Philips N.V. | Systems and methods for generating B-mode images from 3D ultrasound data |
CN109069114B (zh) * | 2017-02-16 | 2021-10-22 | 深圳迈瑞生物医疗电子股份有限公司 | 超声医学检测设备及成像方法、成像系统、显示终端 |
EP3692925A1 (en) * | 2019-02-11 | 2020-08-12 | Koninklijke Philips N.V. | Methods and systems for valve regurgitation assessment |
JP7479935B2 (ja) | 2020-05-26 | 2024-05-09 | キヤノンメディカルシステムズ株式会社 | 体液解析装置、体液解析装置の制御方法、およびプログラム |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4913159A (en) * | 1989-03-17 | 1990-04-03 | Hitachi Medial Corp. | Method for determining blood flow through a narrowed orifice using color doppler echocardiography |
US5062427A (en) * | 1988-05-06 | 1991-11-05 | Kabushiki Kaisha Toshiba | Ultrasonic doppler apparatus |
US5373847A (en) * | 1992-11-12 | 1994-12-20 | University Of Pittsburgh Of The Commonwealth System Of Higher Education | Method of color doppler mapping of blood flow in a patient |
US5846202A (en) * | 1996-07-30 | 1998-12-08 | Acuson Corporation | Ultrasound method and system for imaging |
US6383139B1 (en) * | 1996-06-28 | 2002-05-07 | Sonosite, Inc. | Ultrasonic signal processor for power doppler imaging in a hand held ultrasonic diagnostic instrument |
US6468216B1 (en) * | 2000-08-24 | 2002-10-22 | Kininklijke Philips Electronics N.V. | Ultrasonic diagnostic imaging of the coronary arteries |
US6719697B2 (en) * | 2001-02-27 | 2004-04-13 | Koninklijke Philips Electronics N.V. | Ultrasonic quantification of valvular regurgitant blood flow |
CN101219064A (zh) * | 2006-09-11 | 2008-07-16 | 株式会社东芝 | 超声波诊断装置以及超声波诊断装置的控制方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1178411A1 (ru) * | 1984-02-10 | 1985-09-15 | Всесоюзный научно-исследовательский и испытательный институт медицинской техники | Ультразвуковой импульсно-допплеровский измеритель скорости кровотока |
US5606972A (en) | 1995-08-10 | 1997-03-04 | Advanced Technology Laboratories, Inc. | Ultrasonic doppler measurement of blood flow velocities by array transducers |
US6283919B1 (en) | 1996-11-26 | 2001-09-04 | Atl Ultrasound | Ultrasonic diagnostic imaging with blended tissue harmonic signals |
US6458083B1 (en) | 1996-11-26 | 2002-10-01 | Koninklijke Philips Electronics N.V. | Ultrasonic harmonic imaging with adaptive image formation |
US6013032A (en) | 1998-03-13 | 2000-01-11 | Hewlett-Packard Company | Beamforming methods and apparatus for three-dimensional ultrasound imaging using two-dimensional transducer array |
US5997479A (en) | 1998-05-28 | 1999-12-07 | Hewlett-Packard Company | Phased array acoustic systems with intra-group processors |
US6530885B1 (en) * | 2000-03-17 | 2003-03-11 | Atl Ultrasound, Inc. | Spatially compounded three dimensional ultrasonic images |
US6443896B1 (en) | 2000-08-17 | 2002-09-03 | Koninklijke Philips Electronics N.V. | Method for creating multiplanar ultrasonic images of a three dimensional object |
US6620103B1 (en) * | 2002-06-11 | 2003-09-16 | Koninklijke Philips Electronics N.V. | Ultrasonic diagnostic imaging system for low flow rate contrast agents |
DE102008053073B4 (de) | 2008-10-24 | 2010-08-05 | Tomtec Imaging Systems Gmbh | Dreidimensionale Ableitung einer proximalen isokinetischen Schale einer proximalen Flusskonvergenzzonze sowie dreidimensionale PISA-Flussmessung |
US9204858B2 (en) * | 2010-02-05 | 2015-12-08 | Ultrasonix Medical Corporation | Ultrasound pulse-wave doppler measurement of blood flow velocity and/or turbulence |
-
2011
- 2011-12-16 BR BR112013015817-4A patent/BR112013015817A2/pt not_active Application Discontinuation
- 2011-12-16 US US13/994,794 patent/US10674993B2/en active Active
- 2011-12-16 EP EP11811389.3A patent/EP2654570B1/en not_active Not-in-force
- 2011-12-16 JP JP2013545590A patent/JP6060089B2/ja not_active Expired - Fee Related
- 2011-12-16 WO PCT/IB2011/055739 patent/WO2012085797A1/en active Application Filing
- 2011-12-16 CN CN201180068167.6A patent/CN103391749B/zh not_active Expired - Fee Related
- 2011-12-16 RU RU2013134341/14A patent/RU2569695C2/ru active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5062427A (en) * | 1988-05-06 | 1991-11-05 | Kabushiki Kaisha Toshiba | Ultrasonic doppler apparatus |
US4913159A (en) * | 1989-03-17 | 1990-04-03 | Hitachi Medial Corp. | Method for determining blood flow through a narrowed orifice using color doppler echocardiography |
US5373847A (en) * | 1992-11-12 | 1994-12-20 | University Of Pittsburgh Of The Commonwealth System Of Higher Education | Method of color doppler mapping of blood flow in a patient |
US6383139B1 (en) * | 1996-06-28 | 2002-05-07 | Sonosite, Inc. | Ultrasonic signal processor for power doppler imaging in a hand held ultrasonic diagnostic instrument |
US5846202A (en) * | 1996-07-30 | 1998-12-08 | Acuson Corporation | Ultrasound method and system for imaging |
US6468216B1 (en) * | 2000-08-24 | 2002-10-22 | Kininklijke Philips Electronics N.V. | Ultrasonic diagnostic imaging of the coronary arteries |
US6719697B2 (en) * | 2001-02-27 | 2004-04-13 | Koninklijke Philips Electronics N.V. | Ultrasonic quantification of valvular regurgitant blood flow |
CN101219064A (zh) * | 2006-09-11 | 2008-07-16 | 株式会社东芝 | 超声波诊断装置以及超声波诊断装置的控制方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2014503291A (ja) | 2014-02-13 |
US10674993B2 (en) | 2020-06-09 |
WO2012085797A1 (en) | 2012-06-28 |
BR112013015817A2 (pt) | 2018-05-15 |
EP2654570A1 (en) | 2013-10-30 |
RU2013134341A (ru) | 2015-01-27 |
CN103391749A (zh) | 2013-11-13 |
RU2569695C2 (ru) | 2015-11-27 |
JP6060089B2 (ja) | 2017-01-11 |
US20130274606A1 (en) | 2013-10-17 |
EP2654570B1 (en) | 2017-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103391749B (zh) | 通过超声波成像分析二尖瓣返流 | |
CN103379864B (zh) | 用于超声二尖瓣返流分析的壁滤波器 | |
US5860931A (en) | Ultrasound method and system for measuring perfusion | |
CN104983443B (zh) | 频谱多普勒超声成像中的自动多普勒选通定位 | |
CN110785126B (zh) | 用于超声系统独立衰减系数估计的方法 | |
CN103391748B (zh) | 通过超声成像分析来自缝隙孔的二尖瓣反流 | |
CN103945771A (zh) | 超声波诊断装置以及图像处理装置 | |
US10231693B2 (en) | Automated identification of the location of a regurgitant orifice of a mitral valve in an ultrasound image | |
CN111511286B (zh) | 对血管狭窄的超声测量 | |
Takahashi et al. | Speckle-enhanced cardiac blood flow imaging with high frame rate ultrasound | |
Jensen et al. | P3C-3 In-Vivo Vector Velocity Imaging Using Directional Cross-Correlation | |
Bjastad et al. | Velocity sensitivity mapping in tissue Doppler images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
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: 20150819 Termination date: 20211216 |