CN1086122C - 一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法 - Google Patents
一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法 Download PDFInfo
- Publication number
- CN1086122C CN1086122C CN 95114858 CN95114858A CN1086122C CN 1086122 C CN1086122 C CN 1086122C CN 95114858 CN95114858 CN 95114858 CN 95114858 A CN95114858 A CN 95114858A CN 1086122 C CN1086122 C CN 1086122C
- Authority
- CN
- China
- Prior art keywords
- depth
- value
- degree
- transcranial doppler
- technology
- 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 - Lifetime
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明涉及一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法及仪器,它提供了一种沿声束路径自动搜寻颅内血管位置,动态指示沿声束路径颅内各空间点的血流状态分布情况并定位血管位置的方法,其原理是用多普勒回波信号的能量I2+Q2指示空间点的血流状态,它可帮助医生迅速了解沿某一探头指向的声路径内的颅内血管分布情况和血流状态,使临床医生可以方便准确地进行血管定位和识别,大大降低了诊断难度和诊断时间。
Description
本发明涉及一种经颅多普勒超声脑血流分析仪利用低频脉冲波超声和超声多普勒效应来检测颅内脑血流的状态分布的方法。
在经颅多普勒技术中,超声波探头发射2MHZ的低频超声波脉冲,对于脑内的背向散射回波信号,用下式确定散射源的位置深度:
depth=c·t/2 (1)
其中c为超声波在颅内传播的速度,通常取c=1540m/s,t为发射脉冲与接收的回波信号之间的时延。
目前,常规的经颅多普勒(TCD)仪器,由下列部件:脉冲波超声探头,超声波发射电路,超声波接收前置放大电路,高频放大电路,正交解调电路,带通滤波电路,距离选通电路,A/D转换电路,主控电路,谱分析电路及IBM PC机主机组合构成,主要特点是可以检测分析颅内大血管血流动力学参数,为医生诊断脑血管疾病提供手段。但现有TCD只对采样容积内的回波信号进行分析处理,采样容积的深度根据(1)式设定,通常在仪器开机时设定一个采样容积深度的初始值,医生通过转动探头的方向和通过键盘手动改变采样容积的深度的方法来搜寻血管。由于血管分布于人脑这样一个三维立体空间中,对临床医师来说,搜寻血管成为一项艰苦的工作,具有很大的盲目性。
本发明的目的是克服了现有技术的缺点,提供了一种沿声束路径自动搜寻颅内血管位置,动态指示沿声束路径颅内各空间点的血流状态分布的方法。
为达到上述目的,本发明的技术解决方案是:
一种使用经颅多普勒超声技术检测颅内脑血流状态分布方法,所述方法是在现有的经颅多普勒超声脑血流分析仪的系统程序中增加了血流状态检测并指示,该方法包括下列步骤:
a)设定扫描的起始深度和扫描的终止深度;
b)在上述设定的扫描范围内,以Δ为步距,动态改变扫描深度值;
c)在每个扫描深度处,采样获取多普勒回波信号的同相位成份I、正交相位成份Q值;
d)计算多普勒回波信号能量E;
e)在步骤a)所述的扫描范围内,自动搜索出能量E最大值的位置;
f)在终端设备上显示出能量E随深度的变化,自动锁定能量E最大值的位置。
所述的方法中,步骤d)所述回波信号的能量用E=I2+Q2表示。
所述的方法,还包括下列步骤:
i)在每个扫描深度处重复采样M次I、Q值并计算M组数据,对计算出的M个E值做算术平均,
E作为衡量血流状态的参数。
j)所计算的算术平均值
E做归一化处理。
所述的方法,其中所述参数E的计算可取任何以I2+Q2为根本的变形表达式,I、Q值可以是时域内I、Q信号,也可以是经过付氏变换后的I2+Q2。
所述的方法,其中所述距离选通单元的大小可取为0~(终止深度-初始深度)的任意值,步距Δ的大小可取为0~(终止深度-初始深度)/2的任意值。
所述的方法中,步骤f)包括:
g)在显示器上以横轴为深度,纵轴为E,显示沿声路径不同位置的多普勒回波信号能量,并分别用数字和可动声窗指示出E取最大值的深度点的位置。
h)在输入结束键后,经颅多普勒仪距离选通单元的位置自动放在当前声窗所在的深度上。
本发明采用了上述方法,可以帮助医生迅速了解沿某一探头指向的声路径内的颅内血管分布情况和血流状态,使临床医生可以方便准确地进行血管定位和识别,大大降低了诊断难度和诊断时间。
参照附图对本发明的实施例进行详细的描述。
图1是本发明的所用仪器的系统信号流程图。
图2是本发明动态扫描原理图。
图3为实现动态扫描的硬件电路原理图。
图4是本发明信号I、Q获取原理图。
图5是本发明程序流程图。
图6是本发明的键处理流程图。
图1所示是经颅多普勒超声脑血流分析仪的超声回波信号流程,其原理为:探头接收的回波信号S,经过前级放大、高频放大后,分成两路,分别与同相位(0°)和正交相位(90°)的参考信号相乘,再经过低通滤波得到回波中的多普勒频移成分I1、Q1,I1、Q1滤去血管壁及PRF成分后,在距离选通电路的控制下,对某一深度的多普勒回波信号进行采样,采样数据I、Q送入谱分析仪中,进行复FFT运算,计算其频谱,算出幅度谱
,结果送入IBM PC主机进一步处理及显示。
IBM PC机作为仪器的控制核心,通过写寄存器操作,可以改变控制电路中计数器的初值,从而实现不同的脉冲重复频率、采样容积的大小、扫描深度、探头的发射功率和仪器的接收放大器的增益的控制。
本发明在上述的TCD仪的程序中增加了血流状态声窗程序模块,使距离选通单元在某一深度范围内进行自动扫描,距离选通单元的大小可取为0~(终止深度-初始深度)的任意值。其原理如图2所示(为简化说明,在每个深度处只采样一次):在第1个发射触发脉冲后,令接收延时depth=初始深度;在第二个发射触发脉冲后,延时depth1=depth+Δ,在第3个发射触发脉冲后,延时depth2=depth+2Δ…以此类推,在N个触发脉冲后depthN-1=depth+(N-1)Δ=最大深度,在第N+1个发射触发脉冲后,反向扫描使depthN+1=depth+(N-2)Δ。如此周期性反复扫描,使采样点从初始深度开始到终止深度,以Δ为步距(Δ的大小可取为0~[终止深度-初始深度]/2的任意值)的各个深度点的Dopple回波(I、Q),如图1虚线所示,将采样后的I、Q信号直接送入IBM PC机中,并用下面的公式计算信号的能量:
E=I2+Q2 (2)
计算的结果E做为衡量某深度点的血流状态的参量。
在计算机显示器上,以横轴为深度值,纵轴为E显示沿声路径N个深度位置的血流状态,医生根据每个深度点的E的大小就可以确定该深度是否有血管及血管内的血流状态。
图3为实现动态扫描的硬件电路原理图。颅内血管中的血细胞背向散射的超声回波信号经放大处理后,解调成多普勒回波信号的同相位成份I、正交相位成份Q,经A/D转换为数字量后存于双端口RAM中,双端口RAM的另一端挂在PC总线上,采样后的I、Q数据送到IBM PC起始的内存地址中。图中寄存器74LS273(1,2)联于PC机总线上,PC机CPU以写外设端口的方式对其进行操作,赋给8位计数器(4,5)计数初值,计数器输出通过与非门电路(6,7)产生某一深度的距离选通信号。图4为获取信号I、Q的原理图,I、Q信号经A/D转换为数字量后,结果存于双端口RAM中,由地址计数器产生其写入时序,双端口RAM的另一端挂在PC总线上,占有PC机的内存地址为D800∶0000~D800∶07FF。
在实施时,考虑现有的技术平台,用下面方法实现了上述功能,图5-图6是本发明应用程序流程图,本发明所实现程序(起名为声窗程序)作为TCD仪器系统程序的一个功能模块,以功能键ESC激活调用,在声窗模块内按ESC键,结束声窗模块操作,返回主程序。
进入声窗程序模块后,该模块主要完成下述功能:
1)、获取沿声路径上N个深度位置处的Doppler回波的I、Q值;
在本程序中设定的深度扫描范围在0mm--200mm内任选,深度增加(或减小)的步距Δ为1--100mm内任选;
2)、根据公式(2)计算某一深度点的回波信号,在每个深度点采样M次(取1~50次之内根据需要选择)I、Q值,将M次计算的结果做算术平均:
计算结果
E做为衡量该空间点的血流状态的参数;
3)、在得到了N个深度的
E值后,程序自动找出
E值最大时对应的深度值,退出声窗模块时,该值传给主程序,作为当前的采样容积深度,从而实现了自动查找血管功能;
4)、程序在显示屏幕上以横轴为深度,纵轴为参量
E,显示在每个深度位置的
E值,该值随时间及探头方向的变化而动态改变,程序同时显示一个滑动游标,动态指示
E值最大时对应的深度;
在用纵坐标显示参量E时,如使显示值落于显示屏上窗的高度范围内,对计算的
E值要做归一化处理,处理方法为:
其中EBASE为基底值,ESCALE为比例系数,由于ESCALE可以影响Edisp的变化速度,因此称之为灵敏度系数。
5)、程序在每完成一次沿探头指向的扫描后,查询键盘是否有键按下,如有键按下,则转键处理程序,程序定义了“ESC”键,在声窗程序内按下此键,将退出声窗模块,重新回到系统主程序循环。
程序定义了“↑”/“↓”键,用于改变(4)式中的EBASE值;
程序定义了“←”/“→”键,用于改变(4)式中的ESCALE值;
程序定义了“HOME”/“END”键,用于改变探头的发射功率;
程序定义了“1”/“2”键,用于改变仪器的接收放大器的增益。
Claims (7)
1.一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法,其特征在于:所述方法是在现有的经颅多普勒超声脑血流分析仪的系统程序中增加了血流状态检测并指示,该方法包括下列步骤:
a)设定扫描的起始深度和扫描的终止深度;
b)在上述设定的扫描范围内,以Δ为步距,动态改变扫描深度值;
c)在每个扫描深度处,采样获取多普勒回波信号的同相位成份I、正交相位成份Q值;
d)计算得到所有回波信号能量E;
e)在步骤a)所述的扫描范围内,自动搜索出能量E最大值的位置;
f)在终端设备上显示出能量E随深度的变化,自动锁定能量E最大值的位置。
2.根据权利要求1所述的一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法,其特征在于:步骤d)所述的回波信号的能量用E=I2+Q2表示。
3.根据权利要求1所述的一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法,其特征在于:在每个扫描深度处重复采样M次I、Q值并计算M组数据,对计算出的M个E值做算术平均,
E作为衡量血流状态的参数。
4.根据权利要求3所述的一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法,其特征在于:对所计算的算术平均值
E做归一化处理。
5.根据权利要求2所述的使用经颅多普勒超声技术检测颅内脑血流状态分布的方法,其特征在于:参数E的计算可取任何以I2+Q2为根本的变形表达式,I、Q值可以是时域内I、Q信号,也可以是经过付氏变换后的I2+Q2。
6.根据权利要求1所述的一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法,其特征在于:距离选通单元的大小可取为0~(终止深度-初始深度)的任意值,步距Δ的大小可取为0~(终止深度-初始深度)/2的任意值。
7.根据权利要求1所述的一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法,其特征在于:步骤f)包括:
g)在显示器上以横轴为深度,纵轴为E,显示沿声路径不同位置的多普勒回波信号能量,并分别用数字和可动声窗指示出E取最大值的深度点的位置;
h)在输入结束键后,经颅多普勒仪距离选通单元的位置自动放在当前声窗所在的深度上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 95114858 CN1086122C (zh) | 1995-04-10 | 1995-04-10 | 一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 95114858 CN1086122C (zh) | 1995-04-10 | 1995-04-10 | 一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1133166A CN1133166A (zh) | 1996-10-16 |
CN1086122C true CN1086122C (zh) | 2002-06-12 |
Family
ID=5080343
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 95114858 Expired - Lifetime CN1086122C (zh) | 1995-04-10 | 1995-04-10 | 一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1086122C (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100389726C (zh) * | 2005-02-06 | 2008-05-28 | 微星科技股份有限公司 | 颅部超声波成像方法及系统 |
CN101160098B (zh) * | 2005-05-09 | 2011-01-05 | 株式会社日立医药 | 超声波诊断装置 |
CN101313855B (zh) * | 2007-06-01 | 2010-06-16 | 深圳市德力凯电子有限公司 | 一种用于自动检测脑血流的方法 |
CN100566663C (zh) * | 2007-06-01 | 2009-12-09 | 深圳市德力凯电子有限公司 | 一种脑血流检测探头架 |
CN100577111C (zh) * | 2007-12-29 | 2010-01-06 | 中国人民解放军第四军医大学 | 利用超声多普勒信息显示血流图像的方法及装置 |
JP5398514B2 (ja) * | 2009-12-21 | 2014-01-29 | 株式会社東芝 | カラードプラ超音波診断装置 |
CN103549977A (zh) * | 2013-11-05 | 2014-02-05 | 深圳大学 | 一种经颅多普勒平面环形相控阵探头 |
CN104622505B (zh) * | 2013-11-13 | 2017-04-19 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声颅内血流检测系统和方法 |
CN104665876A (zh) * | 2015-02-10 | 2015-06-03 | 深圳大学 | 经颅多普勒血流检测方法 |
CN105030278B (zh) * | 2015-05-21 | 2017-11-21 | 深圳市德力凯医疗设备股份有限公司 | 一种用于自动扫描颅内脑血管的方法及系统 |
US10206651B2 (en) * | 2015-09-30 | 2019-02-19 | General Electric Company | Methods and systems for measuring cardiac output |
CN109770946A (zh) * | 2019-03-22 | 2019-05-21 | 苏州大学附属第一医院 | 一种微型探头应用于血管吻合术中评估血流动力学的方法 |
-
1995
- 1995-04-10 CN CN 95114858 patent/CN1086122C/zh not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
CN1133166A (zh) | 1996-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1086122C (zh) | 一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法 | |
US6050948A (en) | Ultrasound Doppler diagnostic apparatus | |
US4754760A (en) | Ultrasonic pulse temperature determination method and apparatus | |
CN101176675B (zh) | 超声波诊断装置及超声波图像处理方法 | |
US4622978A (en) | Ultrasonic diagnosing apparatus | |
US4011750A (en) | Method and apparatus for ultrasonic examination of objects | |
US7666142B2 (en) | Ultrasound doppler diagnostic apparatus and image data generating method | |
CN1296012C (zh) | 全数字超声频谱多普勒成像方法及装置 | |
JP3410821B2 (ja) | 超音波診断装置 | |
US6758815B2 (en) | Apparatus and method for indicating mechanical stiffness properties of body tissue | |
JPH06327672A (ja) | 超音波ドプラ診断装置 | |
CN108852416B (zh) | 一种剪切波传播速度的确定方法及装置 | |
JP2640657B2 (ja) | 超音波ドプラ計 | |
JPH021273B2 (zh) | ||
CN1257695A (zh) | 双超声束多普勒血流速度测量方法 | |
JP2916219B2 (ja) | 超音波診断装置 | |
JP2000041983A (ja) | 超音波診断装置 | |
CN1095356C (zh) | 超声血流测量中的“互相关-频谱图”方法 | |
JP2719710B2 (ja) | 超音波パルスドプラ血流診断装置 | |
JP3352211B2 (ja) | 超音波ドプラ診断装置 | |
EP0605720A1 (en) | Inspection apparatus | |
JPH10314171A (ja) | 超音波診断装置 | |
CN116196030A (zh) | 一种档位选择方法及组件、一种超声扫查方法及组件 | |
JP2526623B2 (ja) | 超音波診断装置 | |
EP4362811A1 (en) | System and method for non-invasive determination of pressure in a biological compartment |
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 | ||
PP01 | Preservation of patent right |
Effective date of registration: 20070927 Pledge (preservation): Preservation |
|
PD01 | Discharge of preservation of patent |
Date of registration: 20071019 Pledge (preservation): Preservation |
|
C17 | Cessation of patent right | ||
CX01 | Expiry of patent term |
Expiration termination date: 20150410 Granted publication date: 20020612 |
|
DD01 | Delivery of document by public notice |
Addressee: Ke Xiahe Document name: Notification of Expiration of Patent Right Duration |