CN105708494B - 一种基于超声的血压测量方法 - Google Patents
一种基于超声的血压测量方法 Download PDFInfo
- Publication number
- CN105708494B CN105708494B CN201610057604.9A CN201610057604A CN105708494B CN 105708494 B CN105708494 B CN 105708494B CN 201610057604 A CN201610057604 A CN 201610057604A CN 105708494 B CN105708494 B CN 105708494B
- Authority
- CN
- China
- Prior art keywords
- signal
- sub
- blood
- sampling
- blood flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/04—Measuring blood pressure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Hematology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明涉及超声成像领域,特别涉及一种基于超声的血压测量方法。传统的D模式成像方法(脉冲波多普勒频谱成像方法)相比,本发明提供的血压测量方法通过N个包含血管全深度范围的子采样窗口对采集的IQ信号进行子采样划分,并进一步的通过将每个子采样窗口内的IQ信号进行时间域复合、频域转换、频域复合等操作得血管全深度范围内各处的血流速度信息,进而计算被测患者的血压信息及血压抖动指数;由于本方法同时对血管全深度范围内进行子窗口采样,获取血管全深度范围内各个位置的血流信息,因此,用户可以选择获得整个血管的血流信息分布,从而可以方便精确的计算出被测患者的血压信息及血压抖动指数。
Description
技术领域
本发明涉及超声成像领域,特别涉及一种基于超声的血压测量方法。
背景技术
血压与很多血流疾病直接相关,例如动脉硬化,血管拥塞,脑供血不足等,而传统的超声血液诊断中,医生多采用超声B模式成像结果或D模式成像(又称脉冲波多普勒频谱成像)结果对血流疾病进行诊断,由于技术限制,现有的D模式成像大多仅包含1个子采样窗口,同时该子采样窗口包含的范围非常小(D模式下,只能观测到子采样窗口内的血流信息),因此其仅能获得子窗口局部血流流量,因此无法对被测患者的血压值进行精确计算;同时一些传统的血压测量方法(如汞柱测量法)均为间接测量方法,其测量精度有限。
发明内容
本发明的目的在于克服现有的D模式成像(脉冲波多普勒频谱成像)技术中,由于技术限制,仅有一个子采样窗口而导致仅能获得子采样窗口范围内的血流信息,进而不能够对对被测患者的血压值进行精确计算问题,提供一种包含血管全深度信息的多子采样窗口的血压测量方法。
本文中,血管深度指被测血管横切面直径上任一点到该直径的一个选定端头的距离值。
为了实现上述发明目的,本发明提供了以下技术方案:
一种基于超声的血压测量方法,
包含接收采集自被测血管内的反馈射频信号,并将该射频信号经希尔伯特变换后获得离散IQ信号的步骤;
还包含如下步骤:
将离散IQ信号通过N个子采样门进行子采样划分,N个子采样门包含血管的全部深度范围;N个子采样门可以互不交叉或者略有交叉;N为2以上自然数;
将每个子采样门内的IQ信号顺序进行时间域信号复合、壁滤波器滤波、频域转换,生成频域信号;
从获得的频域信号中计算幅度频谱;
根据幅度频谱计算tj时刻的血流流速分布信息V(k,tj);
将血流流速分布信息V(k,tj)积分得出血流流量;
根据血流流量并结合被测血管的实时截面面积计算血压。
进一步的,子采样门的个数N由血管深度及所述脉冲波波长确定,即由于声学物理条件约束,在一个脉冲波波长内的速度估计是不可再分的,所以目前所能实现的最准确的速度估计即为一个脉冲波波长,所以一个子采样门,为至少一个脉冲波波长的长度或多个脉冲波波长的长度。
优选的,子采样门的个数取符合条件的最大偶数,由上可知,子采样门的最佳长度应为一个脉冲波波长的长度,但为了便于后续的图像显示及计算,优选的将血管深度对称划分,即从血管中心向两侧各有N/2个子采样门。
进一步的,子采样门的大小(即一个脉冲波长内的采样离散信号个数)由公式获得,其中,m是用于检测的脉冲波内包含的脉冲信号个数,其为1以上自然数,Fs是采样频率,Fc是脉冲信号的中心频率,对于每个脉冲信号来说,每个波长内的采样点数=采样频率/(2*信号中心频率);因此对于包含m个脉冲信号的脉冲波来说,当子采样门的长度为信号波长时,每个子采样门内包含的离散采样点数为个,一些实施例下,可以通过公式来确定子采样门N的个数,其中,Nnum为采集到的射频信号的长度,即自血管中采集到的离散信号的个数;在Fs、Fc及被测血管直径相同的前提下,
进一步的,所述时间域信号复合采用公式sk(i)=∑s(i,k)进行,其中i表示IQ信号在时间上的离散序号,k是子采样门的序号;s(i,k)表示第k个子采样门中i时刻的IQ原始信号,sk(i)表示经过时间域符合操作后的信号。
进一步的,将IQ信号转换为频域信号采用公式进行;其中,W表示短时窗口傅里叶变换的窗口大小,表示转换后的频域信号。
进一步的,获取幅度频谱成像的步骤中,所述幅度频谱通过公式获得。
进一步的,根据幅度频谱计算任意tj时刻的血流流速分布信息V(k,tj)采用如下公式:其中,Vi表示tj时刻第k个子采样门幅度频谱Pk(T,ω)对应的血速;
进一步的,将血流流速分布信息V(k,T)积分得出血流流量的步骤中,血流流量:离散表示为其中,r为子采样
门宽度,R为被测血管半径。
进一步的,还包括对血流流量进行角度修正的步骤;修正后的血流流量s表示单心脏跳动周期中的心脏收缩时刻,d表示心脏跳动周期中的心脏舒张时刻;α表示血管壁与子采样门之间的夹角。
进一步的,根据血流流量并结合被测血管实时截面面积计算血压的步骤中,被测血管实时截面面积获得,其中d(t)为被测血管直径随测量时间变化函数;
被测血管在测量时间任意时刻的血压式中,ρ为血流密度,A0为初始血管截面面积,p0为初始血压。
与现有技术相比,本发明的有益效果:与传统的D模式成像方法(脉冲波多普勒频谱成像方法)相比,本发明提供的血压测量方法通过N个包含血管全深度范围的子采样窗口对采集的IQ信号进行子采样划分,并进一步的通过将每个子采样窗口内的IQ信号进行时间域复合、频域转换、频域复合等操作得血管全深度范围内各处的血流速度信息,进而计算被测患者的血压信息及血压抖动指数;由于本方法同时对血管全深度范围内进行子窗口采样,获取血管全深度范围内各个位置的血流信息,因此,用户可以选择获得整个血管的血流信息分布,从而可以方便精确的计算出被测患者的血压信息及血压抖动指数;
与传统的汞柱式血压测量方法等间接测压方法相比,本发明提供的血压测量方法是对血压的直接测量,因此误差更小,抗干扰能力更强。
附图说明:
图1为本发明具体实施例的流程图。
图2即为计算得出的血压与测量时间对应图示例。
具体实施方式
下面结合附图及具体实施例对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
如图1所示,本实施例提供一种基于超声的血压测量方法;包含如下步骤:
S100:接收采集自被测血管内的反馈射频信号,并将该射频信号经希尔伯特变换后获得离散IQ信号的步骤;该反馈射频信号,是按照预设脉冲重复频率发射的超声脉冲波对血管进行检测后的反馈射频信号。
本步骤中,通过希尔伯特变换获得离散IQ信号的过程如下公式表达:
其中,Fc是脉冲波中各个脉冲信号的中心频率,Fs是采样频率,Nnum为采集的射频信号的长度,即自血管中采集到的离散信号的个数;如,当采样频率为20MHz(每秒钟采样20M个采样点)时,声波脉冲波的速率是1540m/s(声音在人体中的传播速度),那么在该采样频率下,每厘米会有260个采样点,假设颈动脉宽度为0.8cm,那么此时,采集的射频信号的长度Nnum=0.8*260=208个采样点。即IQ信号的长度为208。
S101:将所述离散IQ信号通过低通滤波器进行低通滤波,用于过滤掉信号中的直流分量和无效频率。可选的,可采用低通IIR滤波器(Low pass IIR filter)完成本步骤。
S200:将离散IQ信号通过N个子采样门进行子采样划分,N个子采样门包含血管的全部深度范围;N为2以上自然数;优选的,本实施例中,子采样门的个数为满足最大的偶数,且N个子采样门互不交叉;由于血管深度远远大于超声脉冲波波长,因此可以认为本实施例中每个子采样门的门宽即为脉冲波波长;子采样门个数为偶数的优点是,方便后续过程中的计算及图像的显示,优选的将血管深度对称划分,即从血管中心向两侧各有N/2个子采样门。
经过本步骤子采样划分后,每个子采样门中包含的采样点公式获得,其中,m是用于检测的脉冲波内包含的脉冲信号个数,其为1以上自然数,Fs是采样频率,Fc是脉冲信号的中心频率,对于每个脉冲信号来说,每个波长内的采样点数=采样频率/(2*信号中心频率);因此对于包含m个脉冲信号的脉冲波来说,当子采样门的长度为信号波长时,每个子采样门内包含的离散采样点数为
S300:将每个子采样门内的IQ信号顺序进行时间域信号复合、壁滤波器滤波、频域转换,生成频域信号;采用壁滤波器进行滤波的目的为提出信号中的低速血流信息及低速组织运动信息,从而提高后续频谱信号的信噪比。
所述时间域信号复合采用公式sk(i)=∑s(i,k)进行,其中i表示IQ信号在时间上的离散序号,k是子采样门的序号;s(i,k)表示第k个子采样门中i时刻的IQ原始信号,sk(i)表示经过时间域符合操作后的信号。
将IQ信号转换为频域信号为采用快速傅里叶变换实现,具体公式为:其中,W表示短时窗口傅里叶变换的窗口大小,表示转换后的频域信号。
S400:从获得的频域信号中计算幅度频谱后成像;所述幅度频谱通过公式获得。
S500:根据幅度频谱计算tj时刻的血流流速分布信息其中,Vi表示tj时刻第k个子采样门幅度频谱Pk(T,ω)对应的血速;;
S600:将血流流速分布信息V(k,tj)积分得出血流流量离散表示为其中,r为子采样门宽度,R为被测血管半径,优
S601:对血流流量进行角度修正;修正后的血流流量s表示单心脏收缩周期的收缩时刻,d表示但心脏收缩周期的舒张时刻,α表示血管壁与子采样门之间的夹角。
S602:计算血流流量抖动指数,血流流量抖动指数其中Nh为心脏周期数,xi-bf表示单个心脏周期数时的血流流量,ubf为整个测量期间的血流流量平均值,当Indexbf的取值可以表明血流流量是否稳定,如,当Indexbf较大时(例如大于90、100或者120;视具体情况决定),可以认为血流流量不稳定。
S700:根据血流流量并结合被测血管的实时截面面积计算血压;本步骤中,被测血管实时截面面积获得,其中d(t)为被测血管直径随测量时间变化函数;
被测血管在测量时间任意时刻的血压式中,ρ为血流密度,A0为初始血管截面面积,p0为初始血压。
其中,d(t)采用如下方法获得:定义Cm(t)为第m个需要估计位置的血管壁边界,其寻找公式为其中λ为一常量,用于归一化数据,通过对上式的最大值求解即可得到d(t)函数分布情况。图2即为计算得出的血压与测量时间对应图示例。
S701:还包括计算血流压力抖动指数的步骤,从而定量表达被测患者的心脏血压是否稳定,其公式如下其中,xi-bp表示某时刻的血压值,ubp表面检测期间的血压平均值,Nh-bp表示观测次数;Indexbp越大,表明血压越不稳定。
Claims (5)
1.一种基于超声的血流流速分布信息计算方法,
包含接收采集自被测血管内的反馈射频信号,并将该射频信号经希尔伯特变换后获得离散IQ信号的步骤;
其特征在于,还包含如下步骤:
将离散IQ信号通过N个子采样门进行子采样划分,N个子采样门包含被测血管的全部深度范围;N为2以上自然数;
将每个子采样门内的IQ信号顺序进行时间域信号复合、壁滤波器滤波、频域转换,生成频域信号;
从获得的频域信号中计算幅度频谱;
根据幅度频谱计算tj时刻的血流流速分布信息V(k,tj);
其中,子采样门的个数N由血管深度及脉冲波波长确定,即
子采样门的个数取符合条件的最大偶数;
子采样门的大小由公式获得,其中,m是用于检测的脉冲波内包含的脉冲信号个数,其为1以上自然数,Fs是采样频率,Fc是脉冲信号的中心频率。
2.如权利要求1所述的血流流速分布信息计算方法,其特征在于,所述时间域信号复合采用公式sk(i)=∑s(i,k)进行,其中i表示IQ信号在时间上的离散序号,k是子采样门的序号;s(i,k)表示第k个子采样门中i时刻的IQ原始信号,sk(i)表示经过时间域复合操作后的信号。
3.如权利要求2所述的血流流速分布信息计算方法,其特征在于,所述频域转换的步骤为采用短时窗口快速傅里叶变换实现,具体公式为:其中,w表示短时窗口傅里叶变换的窗口大小,表示转换后的频域信号。
4.如权利要求3所述的血流流速分布信息计算方法,其特征在于,获取幅度频谱成像的步骤中,所述幅度频谱通过公式获得。
5.如权利要求4所述的血流流速分布信息计算方法,其特征在于,根据幅度频谱计算任意tj时刻的血流流速分布信息V(k,tj)采用如下公式:其中,Vi表示tj时刻第k个子采样门幅度频谱Pk(t,ω)对应的血速。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610057604.9A CN105708494B (zh) | 2016-01-27 | 2016-01-27 | 一种基于超声的血压测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610057604.9A CN105708494B (zh) | 2016-01-27 | 2016-01-27 | 一种基于超声的血压测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105708494A CN105708494A (zh) | 2016-06-29 |
CN105708494B true CN105708494B (zh) | 2019-04-26 |
Family
ID=56155245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610057604.9A Active CN105708494B (zh) | 2016-01-27 | 2016-01-27 | 一种基于超声的血压测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105708494B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180092622A1 (en) * | 2016-09-30 | 2018-04-05 | Robert Bosch Gmbh | Phased Array for Detecting Artery Location to Measure Blood Velocity |
CN108078590B (zh) * | 2018-01-03 | 2021-02-09 | 声泰特(成都)科技有限公司 | 基于超声频谱多普勒的血流动力学可视化方法与*** |
CN112515704B (zh) * | 2020-12-01 | 2022-07-19 | 声泰特(成都)科技有限公司 | 一种基于超声的血管硬度测量方法 |
CN114287967B (zh) * | 2021-11-24 | 2023-12-29 | 中国科学院深圳先进技术研究院 | 一种基于单阵元超声脉冲回波的血压测量***及方法 |
WO2023092329A1 (zh) * | 2021-11-24 | 2023-06-01 | 中国科学院深圳先进技术研究院 | 一种基于单阵元超声脉冲回波的血压测量***及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011240006A (ja) * | 2010-05-20 | 2011-12-01 | Hitachi Aloka Medical Ltd | 超音波診断装置 |
CN102613990A (zh) * | 2012-02-03 | 2012-08-01 | 声泰特(成都)科技有限公司 | 三维超声频谱多普勒的血流速度及其空间分布显示方法 |
CN104095656A (zh) * | 2014-07-25 | 2014-10-15 | 声泰特(成都)科技有限公司 | 一种基于超声多普勒频谱的彩色血流成像及其显示方法 |
CN104146731A (zh) * | 2014-07-21 | 2014-11-19 | 深圳先进技术研究院 | 一种估算流体速度的***和方法 |
CN105476665A (zh) * | 2016-01-27 | 2016-04-13 | 成都欣声科技有限公司 | 一种基于超声的血流速度测量及血流流量测量方法 |
-
2016
- 2016-01-27 CN CN201610057604.9A patent/CN105708494B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011240006A (ja) * | 2010-05-20 | 2011-12-01 | Hitachi Aloka Medical Ltd | 超音波診断装置 |
CN102613990A (zh) * | 2012-02-03 | 2012-08-01 | 声泰特(成都)科技有限公司 | 三维超声频谱多普勒的血流速度及其空间分布显示方法 |
CN104146731A (zh) * | 2014-07-21 | 2014-11-19 | 深圳先进技术研究院 | 一种估算流体速度的***和方法 |
CN104095656A (zh) * | 2014-07-25 | 2014-10-15 | 声泰特(成都)科技有限公司 | 一种基于超声多普勒频谱的彩色血流成像及其显示方法 |
CN105476665A (zh) * | 2016-01-27 | 2016-04-13 | 成都欣声科技有限公司 | 一种基于超声的血流速度测量及血流流量测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105708494A (zh) | 2016-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105476665B (zh) | 一种基于超声的血流速度测量及血流流量测量方法 | |
CN105708494B (zh) | 一种基于超声的血压测量方法 | |
US7798968B2 (en) | Automatic detection system and method of spectral Doppler blood flow velocity | |
US7611467B2 (en) | Method and apparatus for extracting an envelope curve of a spectrogram | |
KR102025328B1 (ko) | 평면파 합성을 이용한 초음파 벡터 도플러 영상의 생성 장치 및 방법 | |
CN105708496B (zh) | 一种基于超声的血流信息多维成像*** | |
KR101433032B1 (ko) | 평면파를 이용한 기능성 혈류 영상 생성 방법 및 장치 | |
US20090062654A1 (en) | Method and apparatus for automatic optimization of doppler imaging parameters | |
JPS63500148A (ja) | 自己回帰分析を使用して平均反射係数から血液の流れを非侵入的に測定するシステム及び方法 | |
WO2012021765A2 (en) | Methods and apparatus for determining arterial pulse wave velocity | |
JPH0732774B2 (ja) | 心臓サイクルの同期を利用する非侵入血液流測定 | |
US20110270059A1 (en) | Signal processing for pulse oximetry | |
JP2005506120A (ja) | 流路の寸法に関する情報を取得するための方法及びシステム | |
CN108186050A (zh) | 一种基于超声通道数据的多普勒血流速度成像方法和*** | |
JP3745672B2 (ja) | 生体信号計測装置及び超音波診断装置 | |
Thorne et al. | In vivo Doppler ultrasound quantification of turbulence intensity using a high-pass frequency filter method | |
CN106955098B (zh) | 一种血管流速计算方法及装置 | |
JP2013208311A (ja) | 拍動検出装置、電子機器及びプログラム | |
Missaridis et al. | The effect of hemodynamics, vessel wall compliance and hematocrit on ultrasonic Doppler power: an in vitro study | |
Lee et al. | Digital envelope detector for blood pressure measurement using an oscillometric method | |
CN113660901B (zh) | 颅内血压估计方法和设备 | |
JP2003235820A (ja) | 循環動態測定装置 | |
JP7175304B2 (ja) | シャント音解析装置及び方法、コンピュータプログラム並びに記憶媒体 | |
JP7496320B2 (ja) | 心臓の血流と組織の4d超高速ドップラー撮像及び定量化パラメータの取得を同時に行う方法並びに装置。 | |
JP2004222754A (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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20170713 Address after: 610041, No. 4, building 1, No. 1908, 1700, Tianfu Avenue, north section, Tianfu Road, Chengdu, Sichuan. 19 Applicant after: Chengdu Polytechnic Medical Technology Co Ltd Address before: The middle high tech Zone Yizhou road in Chengdu city of Sichuan province 610041 722 No. 2 Building 1 unit 6 floor No. 612 Applicant before: CHENGDU XINSHENG TECHNOLOGY CO., LTD. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |