CN103767737B - 提高彩色血流检测灵敏度的方法和装置 - Google Patents

提高彩色血流检测灵敏度的方法和装置 Download PDF

Info

Publication number
CN103767737B
CN103767737B CN201210407325.2A CN201210407325A CN103767737B CN 103767737 B CN103767737 B CN 103767737B CN 201210407325 A CN201210407325 A CN 201210407325A CN 103767737 B CN103767737 B CN 103767737B
Authority
CN
China
Prior art keywords
variance
energy
speed
index
component signals
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
Application number
CN201210407325.2A
Other languages
English (en)
Other versions
CN103767737A (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.)
Shenzhen Lanying Medical Technology Co ltd
Original Assignee
Shenzhen Landwind Industry 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 Shenzhen Landwind Industry Co Ltd filed Critical Shenzhen Landwind Industry Co Ltd
Priority to CN201210407325.2A priority Critical patent/CN103767737B/zh
Publication of CN103767737A publication Critical patent/CN103767737A/zh
Application granted granted Critical
Publication of CN103767737B publication Critical patent/CN103767737B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开一种提高彩色血流检测灵敏度的方法和装置。其中所述方法包括步骤:将超声诊断设备的超声射频信号处理得到彩色血流初始的速度V、方差T和能量P的分量信号;对速度V、方差T和能量P这三个分量信号的至少其中之一进行约束处理,计算得到相应的索引下标;根据索引下标在预设的索引总表中查询得到索引值;将初始的速度V、方差T和能量P的分量信号与各自对应的索引值相乘,最终的速度V、能量P和方差T的分量信号;利用最终的速度V、能量P和方差T的分量信号显示的彩色血流图像。本发明对血流和非血流的判断更准确,低速血流的检测灵敏度更高。

Description

提高彩色血流检测灵敏度的方法和装置
技术领域
本发明涉及一种涉及医学超声诊断设备的彩色血流成像,尤其是涉及一种提高彩色血流检测灵敏度的方法和装置。
背景技术
超声诊断设备的彩色血流成像是通过检测血流的多普勒频移而得到。在血流检测中,超声设备通过探头向被检测体发射超声波,然后接收被检测体的回波,对接收到的射频回波数据做正交解调、低通滤波、壁滤波、自相关估计等处理,产生被检测部位的速度V、方差T和能量P。
一幅理想的血流图像,能够显示所有流速和能量的血流,同时不显示其它非血流信号(包括组织运动、噪声等)。为了达到这个目的,一般在计算出被检测部位的速度V、方差T和能量P后,会根据血流和非血流信号在速度、能量、方差三个分量各自的特征,直接对三个分量设置不同的阈值,来区分血流和非血流信号。
一般情况下,认为速度小于速度阈值、方差大于方差阈值或者能量小于能量阈值的信号是非血流信号,其它的则为血流信号。这种一刀切的判断方式,容易把低速或者低能量血流误判为非血流信号,从而损失了检测灵敏度。
发明内容
本发明提出一种提高彩色血流检测灵敏度的方法和装置,以解决目前超声诊断设备的彩色血流成像处理过程中低速血流灵敏度不够的技术问题。
本发明采用如下技术方案实现:1、一种提高彩色血流检测灵敏度的方法,包括步骤:
A、将超声诊断设备的超声射频信号处理得到彩色血流初始的速度V、方差T和能量P三个分量信号;
B、对速度V、方差T和能量P这三个分量信号的至少其中之一进行约束处理,计算得到相应的索引下标;
C、根据索引下标在预设的索引总表中查询得到索引值;
D、将初始的速度V、方差T和能量P的分量信号与各自对应的索引值相乘,得到最终的速度V、能量P和方差T的分量信号;
E、用最终的速度V、能量P和方差T来显示彩色血流图像。
其中,步骤B进行约束处理的步骤具体包括:
B1、分别确定速度V、方差T和能量P这三个分量信号各自的感兴趣区域对应的数值分别为[rLowerV,rUpperV]、[rLowerT,rUpperT]和[rLowerP,rUpperP];
B2、将速度V、方差T和能量P这三个分量信号分别约束处理在各自的感兴趣区域。
其中,步骤B2具体包括:
B21、对速度V、能量P和方差T这三个分量信号分别取绝对值;
B22、对速度V、能量P和方差T这三个分量信号的绝对值进行约束处理,约束在感兴趣区域内。
其中,步骤B22中约束处理的步骤具体包括:
对于小于阈值下限rLowerV、rLowerT和rLowerP的速度V、能量P和方差T的分量信号,分别取值为阈值下限rLowerV、rLowerT和rLowerP;
对于大于阈值上限rUpperV、rUpperT和rUpperP的速度V、能量P和方差T的分量信号,则分别取值为阈值上限rUpperV、rUpperT和rUpperP;
对于处于阈值下限与阈值上限之间的速度V、能量P和方差T的分量信号,不作处理。
其中,步骤B中计算得到相应的索引下标的步骤具体包括:按照下列公式进行处理分别得到对应速度V、方差T和能量P的分量信号的索引下标indexV、indexT和indexP:
indexV = [ ( V - rLowerV ) * sV ( rUpperV - rLowerV ) ] ,
indexT = [ ( T - rLowerT ) * sT ( rUpperT - rLowerT ) ] ,
indexP = [ ( P - rLowerP ) * sP ( rUpperP - rLowerP ) ] ,
其中,sV、sT、sP均为大于1的正整数。
其中,索引总表中包括分别对应速度V、方差T和能量P的分量信号的3个相互独立的一维索引表,这3个一维索引表的下标的步长是1,各自下标的取值范围分别为[0,sV-1]、[0,sT-1]和[0,sP-1],且这3个一维索引表中每一点的数值是0或1;其中,sV、sT、sP均为大于1的正整数。
7、一种提高彩色血流检测灵敏度的装置,该装置包括:
分量信号获取模块,用于将超声诊断设备的超声射频信号处理得到彩色血流的速度V、方差T和能量P的分量信号;
索引下标计算模块,用于对速度V、方差T和能量P这三个分量信号的至少其中之一进行约束处理,计算得到相应的索引下标;
索引值查询模块,用于根据索引下标在预设的索引总表中查询得到索引值;
分量信号计算确定模块,用于将初始的速度V、方差T和能量P的分量信号与各自对应的索引值相乘,得到最终的速度V、能量P和方差T三个分量信号;
血流图像显示模块,利用最终的速度V、能量P和方差T三个分量信号来显示彩色血流图像。
其中,索引总表中包括分别对应速度V、方差T和能量P的分量信号的3个相互独立的一维索引表,这3个一维索引表的下标的步长是1,各自下标的取值范围分别为[0,sV-1]、[0,sT-1]和[0,sP-1],且这3个一维索引表中每一点的数值是0或1;其中,sV、sT、sP均为大于1的正整数。
与现有技术相比,本发明具有如下有益效果:
本发明对自相关计算之后的速度V、方差T和能量P,先进行预处理,然后再采用查索引总表的方法确定其是否血流信号。由于索引总表可以任意编辑进行预设和调整,对不同检测部位可设置不同的索引总表,因此,对血流和非血流的判断更准确,低速血流的检测灵敏度更高。
附图说明
图1是本发明一个实施例的流程示意图;
图2是对速度和能量分量的感兴趣区域进行细分的示意图;
图3是本发明一个实施例的结构示意图。
具体实施方式
针对传统方法中低速血流灵敏度不够的问题,本发明对自相关计算之后的速度V、方差T和能量P,先进行预处理,然后再采用查索引总表的方法确定其是否血流信号。由于索引总表可以任意编辑进行预设和调整,对不同检测部位可设置不同的索引总表,因此,对血流和非血流的判断更准确,低速血流的检测灵敏度更高。
如图1所示,本发明的一个实施例具体包括如下实现步骤:
步骤S1、将超声诊断设备的超声射频信号RF处理得到彩色血流的速度V、方差T和能量P这三个分量信号,又称为初始的速度V、方差T和能量P的分量信号。
超声射频信号RF先经过正交解调和低通滤波后得到I(In-phaseComponent)、Q(Quadrate Component)两路信号;然后分别对I、Q信号进行壁滤波、自相关估计后,得到彩色血流的速度V、方差T和能量P这三个分量信号。此步骤为现有技术,不再详细描述具体的处理过程。
步骤S2、从速度V、方差T和能量P这三个分量信号的数值范围中,各自选择感兴趣区域分别对速度V、方差T和能量P这三个分量信号进行约束处理,并计算得到对应速度V、方差T和能量P这三个分量信号各自的索引下标indexV、indexT和indexP。
感兴趣区域一般选择血流信号和非血流信号难以区分的区域。速度V、方差T和能量P这三个分量信号各自的感兴趣区域对应的数值分别为[rLowerV,rUpperV]、[rLowerT,rUpperT]和[rLowerP,rUpperP]。其中,rLowerV、rUpperV、rLowerT、rUpperT、rLowerP和rUpperP这6个参数都是独立可调的,且rLowerV<rUpperV,rLowerT<rUpperT,rLowerP<rUpperP。
由于对速度V、方差T和能量P的分量信号进行约束处理是独立的。根据实际需要,可以只选择这三个分量信号的其中之一、其中之二或这三个分量进行约束处理。比如如图3所示,是以横坐标表示速度V的分量信号、纵坐标表示能量P的分量信号为例,划分出了速度V能量P和的感兴趣区域。
以同时对速度V、方差T和能量P这三个分量信号进行根据划分出的感兴趣区域进行约束处理为例进一步说明,具体的约束处理的步骤包括:首先,对速度V、能量P和方差T这三个分量信号分别取绝对值;然后,对速度V、能量P和方差T这三个分量信号的绝对值按照上述给定的感兴趣区域进行约束处理,即把速度V、能量P和方差T分别限制到[rLowerV,rUpperV]、[rLowerT,rUpperT]、[rLowerP,rUpperP]之内:对于小于阈值下限rLowerV、rLowerT和rLowerP的速度V、能量P和方差T的分量信号,则分别保护为阈值下限rLowerV、rLowerT和rLowerP;对于大于阈值上限rUpperV、rUpperT和rUpperP的速度V、能量P和方差T的分量信号,则分别保护为阈值上限rUpperV、rUpperT和rUpperP;对于处于阈值下限与阈值上限之间的速度V、能量P和方差T的分量信号则保护为该速度V、能量P和方差T的分量信号的实际数值。
以速度V为例,约束处理方法如下:
V = rLowerV , V < rLowerV rUpperV , V > rUpperV V , rLowerV &le; V &le; rUpperV
经过约束处理后速度V、方差T和能量P,分别按照下列公式分别进行处理,得到各自的索引下标indexV、indexT和indexP:
indexV = [ ( V - rLowerV ) * sV ( rUpperV - rLowerV ) ]
indexT = [ ( T - rLowerT ) * sT ( rUpperT - rLowerT ) ]
indexP = [ ( P - rLowerP ) * sP ( rUpperP - rLowerP ) ]
式中,sV、sT、sP均为预定义的值,均为大于1的正整数。
步骤S3、在预定义的索引总表中查找各自的索引下标indexV、indexT和indexP对应的数值,该数值就为各自的索引下标indexV、indexT和indexP分别对应在索引总表中的索引值。
索引总表中包括分别对应速度V、方差T和能量P这三个分量信号的3个一维索引表,这3个一维索引表是分别独立的,且这3个一维索引表的下标取值范围分别为[0,sV-1]、[0,sT-1]和[0,sP-1],且这3个一维索引表中每一点的数值是0或1。
当然,当实际中需要对速度V、方差T和能量P这三个分量信号同时进行约束处理后按步骤S3进行索引值查询时,此时,索引总表中的3个一维索引表可以建立一个以XYZ为坐标轴的三维立方体,XYZ分别为速度V、方差T和能量P三个分量各自对应的一维索引表的下标,因此XYZ的取值范围分别为[0,sV-1]、[0,sT-1]和[0,sP-1],且索引总表中任何一点的数值是0或1,其中0代表非血流,1代表血流。
步骤S4、将索引值与步骤S1得到的速度V、能量P和方差T这三个分量信号相乘,得到最终的速度V、能量P和方差T的分量信号。
步骤S5、将最终的速度V、能量P和方差T这三个分量信号进行常规的时间平滑、空间平滑和优先权编码等处理后,产生用于显示的彩色血流。此步骤为现有技术,不再详述。
结合图3所示,本发明提出提高彩色血流检测灵敏度的装置包括:用于实现步骤S1的分量信号获取模块;用于实现步骤S2的索引下标计算模块;用于实现步骤S3的索引值查询模块;用于实现步骤S4的分量信号计算确定模块;用于实现步骤S5的血流图像显示模块。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种提高彩色血流检测灵敏度的方法,其特征在于,包括步骤:
A、将超声诊断设备的超声射频信号处理得到彩色血流初始的速度V、方差T和能量P三个分量信号;
B、对速度V、方差T和能量P这三个分量信号的至少其中之一进行约束处理,计算得到相应的索引下标;
C、根据索引下标在预设的索引总表中查询得到索引值;
D、将初始的速度V、方差T和能量P的分量信号与各自对应的索引值相乘,得到最终的速度V、能量P和方差T的分量信号;
E、用最终的速度V、能量P和方差T来显示彩色血流图像。
2.根据权利要求1所述提高彩色血流检测灵敏度的方法,其特征在于,步骤B进行约束处理的步骤具体包括:
B1、分别确定速度V、方差T和能量P这三个分量信号各自的感兴趣区域对应的数值分别为[rLowerV,rUpperV]、[rLowerT,rUpperT]和[rLowerP,rUpperP];
B2、将速度V、方差T和能量P这三个分量信号分别约束处理在各自的感兴趣区域。
3.根据权利要求2所述提高彩色血流检测灵敏度的方法,其特征在于,步骤B2具体包括:
B21、对速度V、能量P和方差T这三个分量信号分别取绝对值;
B22、对速度V、能量P和方差T这三个分量信号的绝对值进行约束处理,约束在感兴趣区域内。
4.根据权利要求3所述提高彩色血流检测灵敏度的方法,其特征在于,步骤B22中约束处理的步骤具体包括:
对于小于阈值下限rLowerV、rLowerT和rLowerP的速度V、能量P和方差T的分量信号,分别取值为阈值下限rLowerV、rLowerT和rLowerP;
对于大于阈值上限rUpperV、rUpperT和rUpperP的速度V、能量P和方差T的分量信号,则分别取值为阈值上限rUpperV、rUpperT和rUpperP;
对于处于阈值下限与阈值上限之间的速度V、能量P和方差T的分量信号,不作处理。
5.根据权利要求2所述提高彩色血流检测灵敏度的方法,其特征在于,步骤B中计算得到相应的索引下标的步骤具体包括:按照下列公式进行处理分别得到对应速度V、方差T和能量P的分量信号的索引下标indexV、indexT和indexP:
indexV = [ ( V - rLowerV ) * sV ( rUpperV - rLowerV ) ] ,
indexT = [ ( T - rLowerT ) * sT ( rUpperT - rLowerT ) ] ,
indexP = [ ( P - rLowerP ) * sP ( rUpperP - rLowerP ) ] ,
其中,sV、sT、sP均为大于1的正整数。
6.根据权利要求1所述提高彩色血流检测灵敏度的方法,其特征在于,索引总表中包括分别对应速度V、方差T和能量P的分量信号的3个相互独立的一维索引表,这3个一维索引表的下标的步长是1,各自下标的取值范围分别为[0,sV-1]、[0,sT-1]和[0,sP-1],且这3个一维索引表中每一点的数值是0或1;其中,sV、sT、sP均为大于1的正整数。
7.一种提高彩色血流检测灵敏度的装置,其特征在于,该装置包括:
分量信号获取模块,用于将超声诊断设备的超声射频信号处理得到彩色血流的速度V、方差T和能量P的分量信号;
索引下标计算模块,用于对速度V、方差T和能量P这三个分量信号的至少其中之一进行约束处理,计算得到相应的索引下标;
索引值查询模块,用于根据索引下标在预设的索引总表中查询得到索引值;
分量信号计算确定模块,用于将初始的速度V、方差T和能量P的分量信号与各自对应的索引值相乘,得到最终的速度V、能量P和方差T三个分量信号;
血流图像显示模块,利用最终的速度V、能量P和方差T三个分量信号来显示彩色血流图像。
8.根据权利要求7所述提高彩色血流检测灵敏度的装置,其特征在于,索引总表中包括分别对应速度V、方差T和能量P的分量信号的3个相互独立的一维索引表,这3个一维索引表的下标的步长是1,各自下标的取值范围分别为[0,sV-1]、[0,sT-1]和[0,sP-1],且这3个一维索引表中每一点的数值是0或1;其中,sV、sT、sP均为大于1的正整数。
CN201210407325.2A 2012-10-23 2012-10-23 提高彩色血流检测灵敏度的方法和装置 Active CN103767737B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210407325.2A CN103767737B (zh) 2012-10-23 2012-10-23 提高彩色血流检测灵敏度的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210407325.2A CN103767737B (zh) 2012-10-23 2012-10-23 提高彩色血流检测灵敏度的方法和装置

Publications (2)

Publication Number Publication Date
CN103767737A CN103767737A (zh) 2014-05-07
CN103767737B true CN103767737B (zh) 2015-09-23

Family

ID=50560815

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210407325.2A Active CN103767737B (zh) 2012-10-23 2012-10-23 提高彩色血流检测灵敏度的方法和装置

Country Status (1)

Country Link
CN (1) CN103767737B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101161205A (zh) * 2006-10-13 2008-04-16 深圳迈瑞生物医疗电子股份有限公司 多普勒血流声音信号抗混叠的方法及装置
CN101209210A (zh) * 2006-12-27 2008-07-02 株式会社东芝 超声波图像取得装置
CN101371794A (zh) * 2008-08-05 2009-02-25 深圳市蓝韵实业有限公司 一种彩色多普勒超声诊断仪编码装置和方法
CN101972153A (zh) * 2010-11-10 2011-02-16 中国医学科学院生物医学工程研究所 彩色超声多普勒血流速度估算反正切优化算法
US8260836B2 (en) * 2006-11-28 2012-09-04 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Method and device for generating a filter coefficient in real time

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101138507B (zh) * 2006-09-05 2010-05-12 深圳迈瑞生物医疗电子股份有限公司 多普勒血流声音信号的处理方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101161205A (zh) * 2006-10-13 2008-04-16 深圳迈瑞生物医疗电子股份有限公司 多普勒血流声音信号抗混叠的方法及装置
US8260836B2 (en) * 2006-11-28 2012-09-04 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Method and device for generating a filter coefficient in real time
CN101209210A (zh) * 2006-12-27 2008-07-02 株式会社东芝 超声波图像取得装置
CN101371794A (zh) * 2008-08-05 2009-02-25 深圳市蓝韵实业有限公司 一种彩色多普勒超声诊断仪编码装置和方法
CN101972153A (zh) * 2010-11-10 2011-02-16 中国医学科学院生物医学工程研究所 彩色超声多普勒血流速度估算反正切优化算法

Also Published As

Publication number Publication date
CN103767737A (zh) 2014-05-07

Similar Documents

Publication Publication Date Title
US8523776B2 (en) Ultrasonic doppler imaging apparatus and method with blood velocity waveform processing
US20130006111A1 (en) Ultrasonic diagnostic apparatus
CN105476665A (zh) 一种基于超声的血流速度测量及血流流量测量方法
CN104034733A (zh) 基于双目视觉监测与表面裂纹图像识别的寿命预测方法
EP2790588A1 (en) Automated doppler pulse cycle selection
US20090227869A1 (en) Volume Measurement In An Ultrasound System
CN107525674A (zh) 基于脊线概率分布和局部波动的转频估计方法及检测装置
JP2010075586A (ja) 超音波診断装置
CN105433986A (zh) 超声设备的自动校准方法及系统
CN104873221B (zh) 基于超声波的肝脏脂肪定量方法及系统
CN204595094U (zh) 双通道相位差测量系统
CN111354006A (zh) 超声图像中目标组织的描迹方法及装置
CN105708496A (zh) 一种基于超声的血流信息多维成像系统
CN104127207A (zh) 基于超声成像检测小血流信号灵敏度的方法及系统
EP2248468B1 (en) Ultrasound system and method for adaptively performing clutter filtering
WO2022213938A1 (zh) 超声弹性成像方法、装置、电子设备及存储介质
CN102809687A (zh) 一种交流电频率的数字化测量方法
CN105395214B (zh) 基于rf数据超声成像处理方法及系统
CN103767737B (zh) 提高彩色血流检测灵敏度的方法和装置
CN103505246A (zh) 多普勒参数实时自动标记的方法
US8189427B2 (en) Clutter signal filtering for doppler signal
KR20100119496A (ko) 적응적 클러터 필터링을 수행하는 초음파 시스템 및 방법
CN110801214A (zh) 一种心率实时检测方法及系统
CN115670517A (zh) 胎儿脑室比率的测量方法及超声成像系统
CN104062635A (zh) 一种海杂波图像下基于最小区间思想的海浪主波周期反演方法

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
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20151211

Address after: 518000, Guangdong, Shiyan City, Shenzhen province Baoan District Stone Street neighborhood committee, brick village, Liyuan Industrial Zone, hung long hi tech Industrial Park (blue rhyme Industrial Park), 2, 1 floor

Patentee after: Shenzhen blue ribbon medical imaging Co., Ltd.

Address before: Futian District King Road Shenzhen city Guangdong province 518000 No. 81 country garden building E 601

Patentee before: Landwind Co., Ltd.

CP03 Change of name, title or address

Address after: No.103, baguang District Service Center, No.2 BaiShaWan Road, baguang community, Kuiyong street, Dapeng New District, Shenzhen, Guangdong 518000

Patentee after: Shenzhen Lanying Medical Technology Co.,Ltd.

Address before: 518000 Shenzhen, Baoan District, Guangdong Shiyan street, Lang Xin neighborhood committee, brick factory village, Liyuan Industrial Zone, hung lung hi tech Industrial Park (blue Yun Industrial Park), 1 1 buildings.

Patentee before: SHENZHEN LANYUN MEDICAL IMAGE CO.,LTD.

CP03 Change of name, title or address