CN105476664A - 生成多普勒频谱图的方法和设备 - Google Patents

生成多普勒频谱图的方法和设备 Download PDF

Info

Publication number
CN105476664A
CN105476664A CN201610022328.2A CN201610022328A CN105476664A CN 105476664 A CN105476664 A CN 105476664A CN 201610022328 A CN201610022328 A CN 201610022328A CN 105476664 A CN105476664 A CN 105476664A
Authority
CN
China
Prior art keywords
depth
degree
doppler frequency
blood flow
frequency spectrum
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
CN201610022328.2A
Other languages
English (en)
Other versions
CN105476664B (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.)
BEIJING YUEQICHUANGTONG TECHNOLOGY Co Ltd
Original Assignee
BEIJING YUEQICHUANGTONG TECHNOLOGY Co Ltd
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 BEIJING YUEQICHUANGTONG TECHNOLOGY Co Ltd filed Critical BEIJING YUEQICHUANGTONG TECHNOLOGY Co Ltd
Priority to CN201610022328.2A priority Critical patent/CN105476664B/zh
Publication of CN105476664A publication Critical patent/CN105476664A/zh
Application granted granted Critical
Publication of CN105476664B publication Critical patent/CN105476664B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • 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/0808Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the brain
    • 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/0891Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/486Diagnostic techniques involving arbitrary m-mode
    • 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/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • A61B2576/02Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/026Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Hematology (AREA)
  • Vascular Medicine (AREA)
  • Neurology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明提供一种生成多普勒频谱图的方法和设备。该方法包括:获取M波特征数据;识别所述M波特征数据中的血流数据;根据所述血流数据确定多普勒频谱图应针对的深度;以及基于所确定的深度生成多普勒频谱图。上述生成多普勒频谱图的方法和设备可以自动生成理想的多普勒频谱图,而无需用户干预,大大降低用户工作强度,提高了工作效率。

Description

生成多普勒频谱图的方法和设备
技术领域
本发明涉及医疗器械领域,具体地,涉及一种生成多普勒频谱图的方法和设备。
背景技术
多普勒效应是指物体辐射的波的波长由于该物体和观测者的相对运动而产生变化。在运动的波源前面,波被压缩,波长变得较短,频率变得较高;在运动的波源后面时,会发生相反的现象,波长变得较长,频率变得较低。波源的运动速度越高,所产生的多普勒效应越显著。因此,根据波频率改变的程度,可以计算出波源循着观测方向运动的速度。
多普勒超声血流分析是利用多普勒效应、通过非侵入性的检查评价不同生理学特征的一种方法。经颅多普勒超声分析仪是一种定制化的超声设备,专门用于经颅骨的超声检查。经颅多普勒超声分析仪使用体外超声探头中的晶片经颅骨的缝隙或“窗口”向脑血管发射超声波(简称发射波),血流的存在将导致多普勒效应(多普勒频移)的产生,最后超声波被反射回到探头(简称回波),经同一晶片接收,由分析仪进行数据处理得出相应的血流信息。由于采用深度选通(或距离选通)技术,可进行定点血流测定,因而具有很高的距离分辨力,也可对某点血流的性质做出准确的分析。具体地,如果某一检测深度不存在血液流动,那么不产生多普勒效应,与发射波相比,回波的中心频率不会发生改变;而如果某一检测深度存在血液流动,则会产生多普勒效应,与发射波相比,回波的中心频率会发生偏移。如果使用滤波器把中心频率滤掉,只保留多普勒频偏成分,则对于不存在血液流动的深度,频谱信号中将只剩下背景噪声信号(无用信号),而对于存在血液流动的深度,频谱信号中将包括多普勒信号(有效信号)和无用信号。
正常人体中血管深度相对固定,一个典型范围为40毫米(mm)到75mm,中心值约为50mm。由于体型差异,不同被测者的实际血管深度存在一定差异。此外,用户的操作手法也会影响实际的经颅多普勒超声分析仪的采集深度。为了获得效果更好的多普勒频谱图,现有的经颅多普勒设备需要用户手动调整多普勒频谱图的检测深度。这要求用户进行额外操作,增加了用户具体工作量。此外,手动调整多普勒频谱图还要求用户具有一定的经验,以尽快将其调整合适。
发明内容
为了至少部分地解决现有技术中存在的问题,根据本发明的一个方面,提供了一种生成多普勒频谱图的方法,包括:
获取M波(motionmodeDoppler,或称运动模式多普勒)特征数据;
识别所述M波特征数据中的血流数据;
根据所述血流数据确定多普勒频谱图应针对的深度;以及
基于所确定的深度生成多普勒频谱图。
根据本发明的另一方面,还提供了一种生成多普勒频谱图的设备,包括数据采集装置、识别装置、确定装置和制图装置。数据采集装置用于获取M波特征数据。识别装置用于识别所述M波特征数据中的血流数据。确定装置用于根据所述血流数据确定多普勒频谱图应针对的深度。制图装置用于基于所确定的深度生成多普勒频谱图。
上述生成多普勒频谱图的方法和设备可以自动获得理想的多普勒频谱图,无需用户干预,大大降低用户工作强度,提高了工作效率。
在发明内容中引入了一系列简化的概念,这些概念将在具体实施方式部分中进一步详细说明。本发明内容部分并不意味着要试图限定所要求保护的技术方案的关键特征和必要技术特征,更不意味着试图确定所要求保护的技术方案的保护范围。
以下结合附图,详细说明本发明的优点和特征。
附图说明
本发明的下列附图在此作为本发明的一部分用于理解本发明。附图中示出了本发明的实施方式及其描述,用来解释本发明的原理。在附图中,
图1示出根据本发明一个实施例的经颅多普勒超声分析仪的示意性框图;
图2示出根据本发明一个实施例的多普勒频谱图;
图3示出根据本发明一个实施例的生成多普勒频谱图的方法的示意性流程图;
图4示出根据本发明一个实施例的M波图像;
图5示出根据本发明另一个实施例的生成多普勒频谱图的方法的示意性流程图;
图6示出根据本发明另一个实施例的M波图像;
图7示出根据本发明再一个实施例的M波图像;以及
图8示出根据本发明一个实施例的生成多普勒频谱图的设备的示意性框图。
具体实施方式
在下文的描述中,提供了大量的细节以便能够彻底地理解本发明。然而,本领域技术人员可以了解,如下描述仅涉及本发明的较佳实施例,本发明可以无需一个或多个这样的细节而得以实施。此外,为了避免与本发明发生混淆,对于本领域公知的一些技术特征未进行描述。
图1示出了根据本发明一个实施例的经颅多普勒超声分析仪1000的示意性框图。如图1所示,该经颅多普勒超声分析仪可以包括主机1100和探头1200。主机1100可以包括发射装置1110、探头座1120和处理电路1130。探头1200可以是多个。
发射装置1110用于提供发射信号。发射装置1110可以包括发射控制器1111和驱动电路1112。发射控制器1111可以提供特定频率的脉冲序列作为发射信号。驱动电路1112可以用于将该发射信号转换为高压信号,以驱动探头1200。发射控制器1111的每次发射信号相当于在时间轴上进行了一次采样。该发射信号可以表示为随时间变化的一维信号。
探头1200可以经由探头座1120接收驱动电路1112所转换的高压信号,并将其进行电-声转换,以发射超声波。每次发射后,超声波在人体组织内按照一定速度传播,人体软组织中的平均速度约为1540米/秒(m/s)。超声波发送到人体组织及骨骼后,一部分能量会返回到探头。探头1200还能够接收经被测者反射后的超声波,并进行声-电转换,把包含声速信息的回波转变为电信号,以生成接收信号。
处理电路1130可以包括放大器1131、AD采样电路1132和信号处理模块1133。经探头1200声-电转换生成的接收信号通常较为微弱,因此可选地,处理电路1130可以包括放大器1131,以将微弱的电信号转化为较强的电信号。此外,探头1200产生的电信号为模拟信号,因此可选地,处理电路1130可以包括AD采样电路1132,其通过以采样频率Fs1对模拟信号进行采样,将其转变为数字信号。采样频率Fs1可以称为系统采样率。AD采样电路1132可以按照非常高的时间精度进行采样,这样就可以判断出每个采样点的确切时刻。采样点的时刻和发射时刻的差值为超声波在人体组织中的传播时间,由超声波在人体组织中的平均速度即可得到传播距离。由于探头1200位置固定,超声波的传播距离为来回双程的距离。实际采样点所在的深度为传播距离的一半。由于AD采样电路1132的连续采样,可以得到一系列代表不同深度的采样点,这样就可以得到整个检测深度内所有多普勒信息。多普勒频谱图是根据对于同一深度的采样点进行采样所获得的多普勒信息生成的。另外,处理电路1130中还可以包括信号处理模块1133,其用于对数字信号进行处理,以生成多普勒频谱图。多普勒信号本质上为非平稳信号,随时间变化其频率也会随之变化。自相关是一种常用的用于处理信号f(t)的信号处理方法。它的思想是选择一个时频局部化的分析窗函数g(t),假定分析窗函数g(t)在一个短时间间隔内是平稳(伪平稳)的,移动该分析窗函数g(t),使f(t)g(t)在不同的有限时间宽度内是平稳信号。从而,计算出信号f(t)在各个不同时刻的能量和相位。根据多普勒效应理论,相位信息可以转化为速度信息,因此,可以计算出血流速度。在这里,原始信号f(t)是AD采样电路1132进行模数转换所生成的数字信号。
可选地,主机1100可以连接上位机1300,以由上位机1300显示根据频谱信号生成的多普勒频谱图。多普勒频谱图是三维图像,其中的每条显示线对应于一条功率谱。根据多普勒效应原理,频谱信号中的频率正比于血流速度。所以,多普勒频谱图的横坐标为时间轴,纵坐标为当前检测深度的血流速度。图2示出了根据本发明一个实施例的多普勒频谱图。仅当检测深度对准血管的轴向中心时,多普勒频谱图较好地反映血管的血流情况,其是能够有效帮助用户的理想图。
目前市场上的多普勒分析仪需要用户手动调整多普勒频谱图的检测深度,以获得理想的多普勒频谱图,这给用户带来了操作上的不便。
根据本发明一个方面,提供一种生成多普勒频谱图的方法。图3示出根据本发明一个实施例的生成多普勒频谱图的方法300的示意性流程图。如图3所示,该方法300包括步骤S310、步骤S330、步骤S350和步骤S370。
在步骤S310中,获取M波特征数据。如前所述,经颅多普勒超声分析仪可以获得整个检测深度内的所有多普勒信息。该信息可以称为M波特征数据。M波特征数据不仅包括检测时间信息和血管深度信息,还可以包括以下项中的至少一项:血流速度信息和血流能量信息等。
M波特征数据可以表达为M波图像。M波图像的横轴可以表示时间。纵轴可以表示人体内的深度,例如最上面的像素行表示人体表皮,即深度为0mm,越靠近图像的下部,则表示深度越深。M波图像中的像素的亮度可以表示对应采样点的速度信息或能量信息。换言之,某特定时刻的多个深度的速度信息(能量信息)可以构成M波图像中的一列像素。一个时间段的、多个深度的速度信息(能量信息)可以构成一个M波图像。图4示出了根据本发明一个实施例的M波图像。
M波特征数据可以用于定位血管深度,以便生成多普勒频谱图的时候,可以直接从最佳深度获取信号。从图4中可以看到,整个检测深度范围内,血管只占其中很小的部分。由于个体存在差异,血管实际深度会有明显不同。因此,现有技术中,用户需要频繁手动调节检测深度,才能定位到最佳深度。
在步骤S330中,识别M波特征数据中的血流数据。
血流数据是M波特征数据中与血管部分对应的数据。具体地,可以根据M波特征数据中的血流速度信息或血流能量信息等识别M波特征数据中的血流数据。如图4中间的灰色带状部分所示,M波图像中的血管的部分,即有血流的部分,血流速度较快(能量较高)。所以,可以根据血流速度信息或血流能量信息确定血流数据。例如,如果某采样点的血流速度高于速度阈值,可以认为该采样点为有血流部分,则从该采样点获得的所有数据为血流数据。
在步骤S350中,根据血流数据确定多普勒频谱图应针对的深度。
当经颅多普勒超声分析仪定位到血管中心位置时,可以获得最佳质量的多普勒频谱图。而血流数据是表征了人体中血管部分的特征信息的数据。所以,可以根据血流数据确定多普勒频谱图应针对的深度。
血流数据可以包括血管深度信息和血流速度信息。一般而言,血管中越靠近中心的位置,血流速度越快,所以可以根据血管深度信息和血流速度信息确定期望的深度。例如,选择某时刻的血流速度最大的深度作为多普勒频谱图应针对的深度。对于M波图像来说,即选择某列像素中灰度值最小的像素(颜色最深)所对应的深度作为期望的深度。
血流数据还可能包括血流能量信息。与血流速度信息类似地,该血流能量信息也可以用来确定多普勒频谱图应针对的深度。例如,选择某时刻的血流能量最大的深度作为多普勒频谱图应针对的深度。
根据血流速度信息和血流深度信息能够准确地确定多普勒频谱图应针对的深度,避免噪声干扰,并且方法简单,实现容易。
此外,还可以直接根据血管数据中的血管深度信息确定期望的深度。因为期望检测深度为血管的中心位置所在的深度,所以可以根据血管深度信息确定血管的中心位置。例如,血管深度信息表明当前血管处于50mm至53mm之间的位置,则多普勒频谱图应针对的深度为51.5mm处。从M波图像上说,期望的深度为在垂直于血管方向上的血管的中心。
根据血管深度信息确定多普勒频谱图应针对的深度,方法简单,实现容易。
在该步骤S350中,还可以进一步包括:对M波特征数据进行心动周期分析,以根据心动周期将M波特征数据划分为分别与一个心动周期对应的周期M波特征数据;根据多个周期M波特征数据确定多普勒频谱图应针对的深度。
人体中血液流动的速度一般都是按照心动周期进行变化。心脏在收缩期,其向外射血,血管内血流速度加快;心脏在舒张期,血管内血流速度减慢。因为每个心动周期的血液流动的速度大体相同,所以以心动周期为单位进行M波特征数据的分析将有效保证分析的准确性。
可选的,在根据多个周期M波特征数据确定多普勒频谱图应针对的深度时,首先,根据每个周期M波特征数据确定该周期M波特征数据对应的周期深度。根据上面描述可以理解,该确定操作可以根据血管深度信息或者根据血管深度信息与血流速度信息和/或血流能量信息的组合。然后,计算所有周期深度的平均值,以作为多普勒频谱图应针对的深度。该方法确定的期望深度更接近真实值。
可选的,在根据多个周期M波特征数据确定多普勒频谱图应针对的深度时,首先,根据每个周期M波特征数据确定该周期M波特征数据对应的周期深度。然后,取所有周期深度中的中间值作为多普勒频谱图应针对的深度。假设,一共计算获得了5个周期深度:51mm、51.5mm、52mm、52.3mm和54mm。则认为多普勒频谱图应针对的深度为52mm。该方法确定的期望深度更简单、容易。
可选的,在根据多个周期M波特征数据确定多普勒频谱图应针对的深度时,仅根据正常周期M波特征数据,即在正常的心动周期所获得的周期M波特征数据,确定多普勒频谱图应针对的深度。这种方式直接忽略了多个心动周期中的非正常的心动周期,避免了无用干扰。
总之,对于心率不齐的患者,每个心动周期的M波特征数据可能有较大差异。上述方法可以综合考虑多个心动周期的数据,得到了最可靠的多普勒频谱图应针对的深度。
在步骤S370中,基于所确定的深度生成多普勒频谱图。在此步骤中,根据上述步骤S350所确定的深度生成多普勒频谱图。该多普勒频谱图的采样点较好地对准血管的中心位置。由此该图能够较理想地显示血管的频谱信号,为用户提供有效医疗信息。
上述生成多普勒频谱图的方法300可以在无需用户干预的情况下生成更理想的多普勒频谱图,大大降低用户工作强度,提高了工作效率。
图5示出了根据本发明另一个实施例的生成多普勒频谱图的方法500的流程图。如图5所示,该方法500包括步骤S510、步骤S530、步骤S540、步骤S551、步骤S552、步骤S561、步骤S562和步骤S570。其中步骤S510、步骤S530和步骤S570分别与上述方法300中的对应步骤类似。为了简洁,在此不再详细赘述。步骤S551和步骤S552共同构成步骤S550(未示出),其与上述方法300中的步骤S350类似。
在步骤S550确定多普勒频谱图应针对的深度之前,方法500还包括步骤S540:去除血流数据中的噪声。
在数据采集过程中,被测者和用户不可能长时间保持静止。因此,偶然会出现小的身体动作(比如被测者咳嗽)可能会引发干扰。除此之外,外界各种物理条件也可能产生干扰,比如强电磁干扰。这种干扰信号来源不属于被测者,所以会对步骤S550确定多普勒频谱图应针对的深度造成负面影响,使得方法500中自动生成的多普勒频谱图不理想。图6示出了根据本发明一个实施例的M波图像,如图6所示,该M波图像中包括干扰信号对应的部分,即图中细长部分。
一般来说,干扰信号和正常血流信号所生成的图像有明显差异,例如在以下方面:能量范围和分布,持续时间,形态,周期性规律等。例如,干扰信号会使得M波图像中出现形态高且尖的图形。由此,噪声可以是满足以下条件的信号:在小于时间阈值的持续时间内能量增加超出特定阈值并且所涉及的深度范围明显超过正常血管宽度的信号。例如:在小于100ms的时间段内能量增加超过6dB并且所涉及的深度大于4厘米的信号,则可认为是短时干扰。在正确识别并且去除噪声后,可以保证后续步骤所依据的数据全部为有效信号,为准确确定多普勒频谱图应针对的深度提供有力保证,进而保证所生成的多普勒频谱图满足要求。
M波特征数据采集过程中,经常会出现同时采集了多根血管的信号的情况。如图7所示,该M波图像中,包括两根血管。其中只有一根血管是用户希望观察的,这可以通过血管名称来确定。同一名称的血管,其血流方向是固定的,并且深度在一定范围内。比如大脑中动脉,其血流方向为正向,深度范围在40mm到75mm之间。由深度范围和血流方向两个条件,可以判定出哪一根血管为用户需要。目标血管识别后,可以据此确定期望深度。
可选地,方法500可以包括步骤S551和步骤S552。在步骤S551中,根据目标血管的血流方向和深度范围确定血流数据中的目标血流数据。目标血管名称可以是通过输入装置接收的。基于目标血管名称可以确定该目标血管中的血流方向和深度范围。每个血管的血流方向是固定的。而在检测时,如果用户遵守检测规范,那么在M波特征数据中的目标血管的血流方向就是固定的。目标血管的深度范围也是固定的。由此,基于目标血管的血流方向和深度范围即可确定目标血流数据,即目标血管的相关数据。在步骤S552中,根据目标血流数据确定多普勒频谱图应针对的深度。上述步骤S551和步骤S552有效保证了目标深度准确性。由此,可以满足用户要求,提高用户体验。
方法500还包括步骤S561和步骤S562。
在步骤S561中,根据血流数据计算打图质量参数。探头所对的方向与血管中心偏离越小,图越打的正,获得的图像质量越好。
一个有效参数为血管宽度。一般来说,血管截面打得越正,血管宽度越大。血管宽度和取样容积有一定关系,在打正情况下,血管宽度要高于取样容积。另一个有效参数是平均血流能量,平均血流能量越大,血管截面打的越正。由于血管中心流速较高,因此打正情况下,可以测得血管的最高血流速度,这样最高血流速度也是一个重要判断参数。此外,在打偏的情况下,血管充盈度可能较差,因此,血管充盈度也是一个重要参数。
可以为上述每个参数可以设置一个参考值,实际数值和参考值的比值作为一个特征值。比如取样容积为16mm,实际血管宽度为12mm,其比值为0.75。可以设置参数取值范围在0至1之间,0表示质量最差,1表示质量最好。这样,每个参数可以转为一个特征。
优选地,将数据参数转化为用户可以接受的参数指标,将该参数指标反馈给用户更有价值。可选地,将所有特征值组成一个特征向量,求取其二范数。求取二范数时,可以对每个特征进行加权。当然,也可以将每个特征的权重设为相同。这个二范数的取值范围也是0至1。0可以代表图像质量最差;1可以代表图像质量最佳。可以理解,当取值0至0.3时,可以视图像质量为不可接受;当取值0.4至0.6时,可以视图像质量为中等;当取值0.7至1.0时,可以视图像质量为较好。这样就可以帮助用户判断操作手法是否合适,提高其工作效率。
在步骤S562中,提供上述打图质量参数。可以理解,可以直接提供上述血管宽度、平均血流能量、最高血流速度等参数的数值。还可以提供根据前述参数值所计算的参数指标,以利于用户理解。在这里,自动反馈了血管打正信息,帮助用户判断手法是否合适,而无需用户依据经验判断,利于初级用户操作。
本领域普通人员可以理解,在上述方法500中,以步骤S510、步骤S530、步骤S540、步骤S551、步骤S552、步骤S561、步骤S562和步骤S570的顺序进行了描述。但该实现方式仅是为了说明本发明的实施例的示例,其不对本发明造成限制。
根据本发明另一方面,还提供一种生成多普勒频谱图的设备。图8示出了根据本发明一个实施例的生成多普勒频谱图的设备800。如图8所示,设备800包括数据采集装置810、识别装置830、确定装置850和制图装置870。其中,数据采集装置810用于获取M波特征数据。识别装置830用于识别M波特征数据中的血流数据。确定装置850用于根据血流数据确定多普勒频谱图应针对的深度。制图装置870用于基于所确定的深度生成多普勒频谱图。
上述血流数据可以包括血管深度信息和血流速度信息。上述确定装置850可以根据血管深度信息和血流速度信息确定多普勒频谱图应针对的深度。
上述血流数据可以包括血管深度信息和血流能量信息。上述确定装置850可以根据血管深度信息和血流能量信息确定多普勒频谱图应针对的深度。
上述血流数据包括血管深度信息。上述确定装置850可以仅根据血管深度信息确定多普勒频谱图应针对的深度。
可选地,上述确定装置850进一步包括目标确定模块和第一深度确定模块。目标确定模块用于根据目标血管的血流方向和深度范围确定血流数据中的目标血流数据。第一深度确定模块用于根据目标血流数据确定多普勒频谱图应针对的深度。
可选地,上述确定装置850进一步包括周期划分模块和第二深度确定模块。周期划分模块用于对M波特征数据进行心动周期分析,以根据心动周期将M波特征数据划分为分别与一个心动周期对应的周期M波特征数据。第二深度确定模块用于根据多个周期M波特征数据确定多普勒频谱图应针对的深度。
可选地,上述第二深度确定模块进一步包括周期深度确定单元,用于根据每个周期M波特征数据确定对应的周期深度。第二深度确定模块还包括平均单元,用于计算所有周期深度的平均值,以作为多普勒频谱图应针对的深度。
可选地,设备800进一步包括去噪装置,用于去除血流数据中的噪声。
可选地,设备800进一步包括质量参数计算装置和质量参数提供装置。其中质量参数计算装置用于根据血流数据计算打图质量参数。质量参数提供装置用于提供打图质量参数。
可选地,质量参数计算装置包括:特征值计算模块,用于根据血流数据计算包括以下项中的至少一项的特征值:血管宽度、血管充盈度、最高血流速度和平均血流能量。质量参数计算装置还包括质量参数获得模块,用于根据由特征值构成的特征向量计算打图质量参数。
本领域普通技术人员通过阅读上文关于生成多普勒频谱图的方法的详细描述,能够理解上述生成多普勒频谱图的设备的结构、实现以及优点,因此这里不再赘述。
本发明已经通过上述实施例进行了说明,但应当理解的是,上述实施例只是用于举例和说明的目的,而非意在将本发明限制于所描述的实施例范围内。此外本领域技术人员可以理解的是,本发明并不局限于上述实施例,根据本发明的教导还可以做出更多种的变型和修改,这些变型和修改均落在本发明所要求保护的范围以内。本发明的保护范围由附属的权利要求书及其等效范围所界定。

Claims (11)

1.一种生成多普勒频谱图的方法,包括:
获取M波特征数据;
识别所述M波特征数据中的血流数据;
根据所述血流数据确定多普勒频谱图应针对的深度;以及
基于所确定的深度生成多普勒频谱图。
2.根据权利要求1所述的方法,其特征在于,所述确定多普勒频谱图应针对的深度进一步包括:
根据目标血管的血流方向和深度范围确定所述血流数据中的目标血流数据;以及
根据所述目标血流数据确定所述多普勒频谱图应针对的深度。
3.根据权利要求1所述的方法,其特征在于,所述确定多普勒频谱图应针对的深度进一步包括:
对所述M波特征数据进行心动周期分析,以根据心动周期将所述M波特征数据划分为分别与一个心动周期对应的周期M波特征数据;
根据多个周期M波特征数据确定所述多普勒频谱图应针对的深度。
4.根据权利要求3所述的方法,其特征在于,其中,所述根据多个周期M波特征数据确定所述多普勒频谱图应针对的深度进一步包括:
根据每个周期M波特征数据确定对应的周期深度;以及
计算所有周期深度的平均值,以作为所述多普勒频谱图应针对的深度。
5.根据权利要求1所述的方法,其特征在于,其中,所述血流数据包括血管深度信息和血流速度信息,所述确定多普勒频谱图应针对的深度进一步包括:
根据所述血管深度信息和所述血流速度信息确定所述多普勒频谱图应针对的深度。
6.根据权利要求1所述的方法,其特征在于,其中,所述血流数据包括血管深度信息和血流能量信息,所述确定多普勒频谱图应针对的深度进一步包括:
根据所述血管深度信息和所述血流能量信息确定所述多普勒频谱图应针对的深度。
7.根据权利要求1所述的方法,其特征在于,其中,所述血流数据包括血管深度信息,所述确定多普勒频谱图应针对的深度进一步包括:
根据所述血管深度信息确定所述多普勒频谱图应针对的深度。
8.根据权利要求1所述的方法,其特征在于,在所述确定多普勒频谱图应针对的深度之前,所述方法进一步包括:
去除所述血流数据中的噪声。
9.根据权利要求1所述的方法,其特征在于,在所述识别所述M波特征数据中的血流数据之后,所述方法进一步包括:
根据所述血流数据计算打图质量参数;以及
提供所述打图质量参数。
10.根据权利要求9所述的方法,其特征在于,所述根据所述血流数据计算打图质量参数进一步包括:
根据所述血流数据计算包括以下项中的至少一项的特征值:血管宽度、血管充盈度、最高血流速度和平均血流能量;以及
根据由所述特征值构成的特征向量计算所述打图质量参数。
11.一种生成多普勒频谱图的设备,包括:
数据采集装置,用于获取M波特征数据;
识别装置,用于识别所述M波特征数据中的血流数据;
确定装置,用于根据所述血流数据确定多普勒频谱图应针对的深度;以及
制图装置,用于基于所确定的深度生成多普勒频谱图。
CN201610022328.2A 2016-01-13 2016-01-13 生成多普勒频谱图的方法和设备 Active CN105476664B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610022328.2A CN105476664B (zh) 2016-01-13 2016-01-13 生成多普勒频谱图的方法和设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610022328.2A CN105476664B (zh) 2016-01-13 2016-01-13 生成多普勒频谱图的方法和设备

Publications (2)

Publication Number Publication Date
CN105476664A true CN105476664A (zh) 2016-04-13
CN105476664B CN105476664B (zh) 2018-05-01

Family

ID=55664192

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610022328.2A Active CN105476664B (zh) 2016-01-13 2016-01-13 生成多普勒频谱图的方法和设备

Country Status (1)

Country Link
CN (1) CN105476664B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111110278A (zh) * 2019-12-30 2020-05-08 深圳市德力凯医疗设备股份有限公司 一种采集参数的配置方法、存储介质及超声设备
WO2020133332A1 (zh) * 2018-12-29 2020-07-02 深圳迈瑞生物医疗电子股份有限公司 确定超声波发射角度的方法以及超声设备
CN116831537A (zh) * 2023-07-13 2023-10-03 逸超医疗科技(北京)有限公司 一种脉搏波传导速度成像与超声多普勒频谱自动同步方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6423006B1 (en) * 2000-01-21 2002-07-23 Siemens Medical Solutions Usa, Inc. Method and apparatus for automatic vessel tracking in ultrasound systems
CN1738575A (zh) * 2003-02-18 2006-02-22 松下电器产业株式会社 超声波多普勒血流测量装置
US7534209B2 (en) * 2000-05-26 2009-05-19 Physiosonics, Inc. Device and method for mapping and tracking blood flow and determining parameters of blood flow
JP2009213593A (ja) * 2008-03-10 2009-09-24 Panasonic Corp 超音波診断装置
CN102038522A (zh) * 2009-10-21 2011-05-04 株式会社东芝 超声波诊断装置以及超声波数据取得方法
CN103169503A (zh) * 2013-03-22 2013-06-26 北京汇影互联科技有限公司 一种脑部血流测量系统及方法
CN104622505A (zh) * 2013-11-13 2015-05-20 深圳迈瑞生物医疗电子股份有限公司 一种超声颅内血流检测系统和方法
CN105030278A (zh) * 2015-05-21 2015-11-11 深圳市德力凯电子有限公司 一种用于自动扫描颅内脑血管的方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6423006B1 (en) * 2000-01-21 2002-07-23 Siemens Medical Solutions Usa, Inc. Method and apparatus for automatic vessel tracking in ultrasound systems
US7534209B2 (en) * 2000-05-26 2009-05-19 Physiosonics, Inc. Device and method for mapping and tracking blood flow and determining parameters of blood flow
CN1738575A (zh) * 2003-02-18 2006-02-22 松下电器产业株式会社 超声波多普勒血流测量装置
JP2009213593A (ja) * 2008-03-10 2009-09-24 Panasonic Corp 超音波診断装置
CN102038522A (zh) * 2009-10-21 2011-05-04 株式会社东芝 超声波诊断装置以及超声波数据取得方法
CN103169503A (zh) * 2013-03-22 2013-06-26 北京汇影互联科技有限公司 一种脑部血流测量系统及方法
CN104622505A (zh) * 2013-11-13 2015-05-20 深圳迈瑞生物医疗电子股份有限公司 一种超声颅内血流检测系统和方法
CN105030278A (zh) * 2015-05-21 2015-11-11 深圳市德力凯电子有限公司 一种用于自动扫描颅内脑血管的方法及系统

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020133332A1 (zh) * 2018-12-29 2020-07-02 深圳迈瑞生物医疗电子股份有限公司 确定超声波发射角度的方法以及超声设备
CN112672694A (zh) * 2018-12-29 2021-04-16 深圳迈瑞生物医疗电子股份有限公司 确定超声波发射角度的方法以及超声设备
CN111110278A (zh) * 2019-12-30 2020-05-08 深圳市德力凯医疗设备股份有限公司 一种采集参数的配置方法、存储介质及超声设备
CN116831537A (zh) * 2023-07-13 2023-10-03 逸超医疗科技(北京)有限公司 一种脉搏波传导速度成像与超声多普勒频谱自动同步方法

Also Published As

Publication number Publication date
CN105476664B (zh) 2018-05-01

Similar Documents

Publication Publication Date Title
KR101501479B1 (ko) 초음파 최적화 방법과 그를 위한 초음파 의료 장치
JP4667394B2 (ja) 超音波診断装置
US8864671B2 (en) Methods and systems for color flow imaging
KR20170085516A (ko) 전단파 탄성초음파 영상화를 수행하기 위한 영상화 방법 및 장치
JP2013523325A (ja) 超音波撮像の方法および装置
US20100138191A1 (en) Method and system for acquiring and transforming ultrasound data
WO2009118798A1 (ja) 超音波診断装置
US20110092819A1 (en) Ultrasonic diagnosis apparatus and ultrasoinc data acquisition method
CN105662473A (zh) 组织参数检测方法和系统
KR20160073168A (ko) 초음파 진단장치 및 그에 따른 초음파 진단 장치의 동작 방법
CN105476664A (zh) 生成多普勒频谱图的方法和设备
KR101652727B1 (ko) 초음파 도플러 영상 장치 및 그 제어 방법
US11504087B2 (en) Ultrasonic diagnostic device and control method therefor
CN105476661A (zh) 调整多普勒频谱图的方法和设备
CN117064441A (zh) 一种超声成像方法以及超声成像系统
CN112168210B (zh) 医学图像处理终端、超声诊断设备和胎儿图像的处理方法
US20160296201A1 (en) Imaging and measuring system of vocal cord vibration based on plane wave ultrasonography, and method thereof
CN105662472A (zh) 生成多普勒频谱图的方法和设备
KR101051567B1 (ko) 표준 단면 정보를 제공하는 초음파 시스템 및 방법
KR20190070883A (ko) Arfi 이미징을 위한 교정
CN114025672B (zh) 一种超声成像设备和子宫内膜蠕动的检测方法
CN105476663A (zh) 调整多普勒频谱图的方法和设备
KR102223058B1 (ko) 초음파 시스템 및 이의 클러터 필터링 방법
US11602328B2 (en) Ultrasound flow imaging
WO2007119609A1 (ja) 塞栓子観測方法および装置ならびにこれらを用いた超音波診断装置

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