CN109259801A - 一种剪切波弹性成像方法及装置 - Google Patents
一种剪切波弹性成像方法及装置 Download PDFInfo
- Publication number
- CN109259801A CN109259801A CN201811063847.9A CN201811063847A CN109259801A CN 109259801 A CN109259801 A CN 109259801A CN 201811063847 A CN201811063847 A CN 201811063847A CN 109259801 A CN109259801 A CN 109259801A
- Authority
- CN
- China
- Prior art keywords
- interest
- area
- target component
- value
- shearing wave
- 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
Links
- 238000010008 shearing Methods 0.000 title claims abstract description 152
- 238000003384 imaging method Methods 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000001514 detection method Methods 0.000 claims description 105
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 210000001519 tissue Anatomy 0.000 description 47
- 230000005284 excitation Effects 0.000 description 21
- 230000003902 lesion Effects 0.000 description 15
- 230000005855 radiation Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 230000000644 propagated effect Effects 0.000 description 8
- 239000000523 sample Substances 0.000 description 7
- 210000004872 soft tissue Anatomy 0.000 description 7
- 238000003745 diagnosis Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000003786 synthesis reaction Methods 0.000 description 4
- 238000002604 ultrasonography Methods 0.000 description 4
- 208000009453 Thyroid Nodule Diseases 0.000 description 3
- 208000024770 Thyroid neoplasm Diseases 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 235000017284 Pometia pinnata Nutrition 0.000 description 1
- 240000007653 Pometia tomentosa Species 0.000 description 1
- 241000219000 Populus Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000002091 elastography Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000009940 knitting Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Clinical applications
-
- 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/485—Diagnostic techniques involving measuring strain or elastic properties
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本申请公开了一种剪切波弹性成像方法及装置,其中,方法包括:获取目标参数的目标取值;所述目标参数为用于对感兴趣区域进行剪切波成像的参数;所述目标取值为与所述感兴趣区域的组织特性相匹配的取值;将所述目标参数的当前取值变更为所述目标取值;基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像。通过本申请实施例,使得依据目标取值对感兴趣区域进行剪切波成像得到的成像结果可靠性更高。
Description
技术领域
本申请涉及医学成像领域,特别是涉及一种剪切波弹性成像方法及装置。
背景技术
基于声辐射力的剪切波弹性成像是一种评估组织弹性的超声弹性成像技术,其被广泛应用于人体软组织病变的分析和诊断中。基于声辐射力的剪切波弹性成像技术的原理为:由探头向人体的软组织发射高能量的超声波后,在声辐射力和组织的剪切应力的作用下,特定区域内的软组织会产生向四周传播的振动,从而产生剪切波。而由于人体软组织的弹性与剪切波的速度之间存在着密切的关系,因此可以通过检测剪切波的传播速度来评估软组织的弹性,进而确定是否存在软组织的相关病变。
目前的剪切波弹性成像设备在用户设定好感兴趣区域后,系统再执行既定的超声声辐射力激励和剪切波检测发射序列,从而获得包含组织运动信息的回波数据。随后利用组织运动估计算法得到包含剪切波运动信息的数据,然后再利用滤波器滤除运动噪声,最后利用时间飞跃法和/或互相关等算法求得剪切波经过感兴趣区域中各位置点的传播速度,并依据关系式,将剪切波速度转换为剪切模量和杨氏模量等参数,并依据剪切模量和杨氏模量评估软组织的弹性。
但是,现有技术中对感兴趣区域进行剪切波弹性成像的结果可靠性较低。
发明内容
基于此,本申请提出了一种剪切波弹性成像方法,用以提高对感兴趣区域的剪切波弹性成像结果的可靠性。
本申请还提供了一种剪切波弹性成像装置,用以保证上述方法在实际中的实现及应用。
本申请提供的技术方案为:
本申请公开了一种剪切波弹性成像方法,包括:
获取目标参数的目标取值;所述目标参数为用于对感兴趣区域进行剪切波成像的参数;所述目标取值为与所述感兴趣区域的组织特性相匹配的取值;
将所述目标参数的当前取值变更为所述目标取值;
基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像。
其中,所述目标参数包括:用于将所述感兴趣区域划分为检测区域的行数和列数;
所述获取目标参数的目标取值,包括:
确定所述感兴趣区域的长度和宽度;所述长度为所述感兴趣区域在剪切波的传播方向上的尺度;所述宽度为所述感兴趣区域在轴向方向上的尺度;所述轴向方向为平行于所述感兴趣区域且垂直于所述传播方向的方向;
依据剪切波在传播方向上的预设传播距离,以及在轴向方向上的预设聚焦范围,确定将所述感兴趣区域划分为满足预设条件的检测区域的行数和列数;所述预设条件包括:在所述传播方向上的长度与所述预设传播距离间的差值小于预设第一阈值,且在所述轴向方向上的宽度与所述预设聚焦范围间的差值小于预设第二阈值。
其中,所述获取目标参数的目标取值,包括:
接收用户对所述目标参数的第一调整数值;
确定所述第一调整数值为所述目标参数的目标取值。
其中,在所述基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像之后,还包括:
在接收到调整指令的情况下,接收用户对所述目标参数的第二调整数值;
将所述目标参数的当前取值变更为所述第二调整数值;
基于所述第二调整数值,对所述感兴趣区域进行剪切波弹性成像。
其中,在所述接收用户对所述目标参数的第二调整数值之前,还包括:
显示所述目标参数的当前取值。
其中,在所述获取目标参数的目标取值之前,还包括:
对待检测组织进行多模式成像,得到所述待检测组织的多模式图像;
显示所述多模式图像,使得用户依据所述多模式图像确定所述感兴趣区域。
本申请还公开了一种剪切波弹性成像装置,包括:
获取单元,用于获取目标参数的目标取值;所述目标参数为用于对感兴趣区域进行剪切波成像的参数;所述目标取值为与所述感兴趣区域的组织特性相匹配的取值;
第一变更单元,用于将所述目标参数的当前取值变更为所述目标取值;
第一成像单元,用于基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像。
其中,所述获取单元所获取的目标参数包括:用于将所述感兴趣区域划分为检测区域的行数和列数;
所述获取单元,包括:
第一确定子单元,用于确定所述感兴趣区域的长度和宽度;所述长度为所述感兴趣区域在剪切波的传播方向上的尺度;所述宽度为所述感兴趣区域在轴向方向上的尺度;所述轴向方向为平行于所述感兴趣区域且垂直于所述传播方向的方向;
第二确定子单元,用于依据剪切波在传播方向上的预设传播距离,以及在轴向方向上的预设聚焦范围,确定将所述感兴趣区域划分为满足预设条件的检测区域的行数和列数;所述预设条件包括:在所述传播方向上的长度与所述预设传播距离间的差值小于预设第一阈值,且在所述轴向方向上的宽度与所述预设聚焦范围间的差值小于预设第二阈值。
其中,所述获取单元还包括:
接收子单元,用于接收用户对所述目标参数的第一调整数值;
第三确定子单元,用于确定所述第一调整数值为所述目标参数的目标取值。
其中,还包括:
接收单元,用于在所述第一成像单元基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像之后,在接收到调整指令的情况下,接收用户对所述目标参数的第二调整数值;
第二变更单元,用于将所述目标参数的当前取值变更为所述第二调整数值;
第二成像单元,用于基于所述第二调整数值,对所述感兴趣区域进行剪切波弹性成像。
其中,还包括:
第一显示单元,用于在所述接收单元接收用户对所述目标参数的第二调整数值之前,显示所述目标参数的当前取值。
其中,还包括:
多模式成像单元,用于对待检测组织进行多模式成像,得到所述待检测组织的多模式图像;
第二显示单元,用于显示所述多模式图像,使得用户依据所述多模式图像确定所述感兴趣区域。
本申请的有益效果为:
剪切波弹性成像结果的可靠性的原因为:现有的剪切波弹性成像设备采用固定的激励与检测参数(成像参数),固定的激励与检测参数不能适用于所有病人的组织特性,进而对于某些病人的感兴趣区域不能呈现出用户预估的成像效果(用户依据多模式图像对感兴趣组织区域病变程度,预估剪切波弹性成像的结果),进而导致成像结果的可靠性较低。
本申请实施例中,将目标参数的当前取值变更为目标取值,表明本申请实施例中的剪切波弹性成像设备中的目标参数的取值可调,因此,针对任意的感兴趣区域,都可以获取目标参数的目标取值,并将目标参数的当前取值变更为与当前感兴趣区域的组织特性相匹配的目标取值,使得依据目标取值对感兴趣区域进行剪切波成像得到的成像结果可靠性更高。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本申请中一种剪切波弹性成像方法实施例的流程图;
图2为本申请中又一种剪切波弹性成像方法实施例的流程图;
图3为本申请中一种超声声辐射力的产生示意图;
图4为本申请中传统单焦点记录的发射声场示意图;
图5为本申请中又一种剪切波弹性成像方法实施例的流程图;
图6为本申请中一种将感兴趣区域划分为9个检测区域的结果示意图;
图7为本申请中一种将感兴趣区域划分为4个检测区域的结果示意图;
图8为本申请中一种剪切波弹性成像装置实施例的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
本申请实施例可以应用于剪切波弹性成像设备,也可以独立设置。
图1为本申请中一种剪切波弹性成像方法,可以包括:
步骤101:获取目标参数的目标取值。
在进入剪切波弹性成像后,执行本实施例的动作。具体的,在本步骤中,目标参数为用于对感兴趣区域进行剪切波弹性成像的参数,可以包括:将感兴趣区域划分为检测区域的行数和列数、对检测区域进行激励的激励参数,以及对检测区域进行检测的检测参数。其中,激励参数可以包括:激励脉冲周期数、激励发射频率和激励发射F值等;检测参数可以包括检测脉冲周期数、检测发射频率、检测发射F值以及目标检测区域的数量等。
在本步骤中,目标参数的目标取值为与感兴趣区域的组织特性相匹配的取值。对于目标参数的目标取值的获取方式,可以依据感兴趣区域确定目标参数,也可以接收用户设定的目标取值。对于依据感兴趣区域确定目标参数的目标取值的过程在图2对应的实施例中进行详细介绍;对于接收用户确定的目标取值的过程在图5对应的实施例中进行详细介绍。
步骤102:将目标参数的当前取值变更为目标取值。
在本实施例中,剪切波弹性成像设备中对每个目标参数都设置有当前取值,在本步骤中,将目标参数的当前取值变更为目标取值。具体的,剪切波弹性成像设备将目标参数的当前取值变更为目标取值的过程可以包括:将目标参数的当前取值与目标参数的目标取值进行对比,将目标参数中不同于目标取值的当前取值变更为目标取值;还可以包括:直接将目标参数的当前取值全部变更为目标取值。本实施例不对具体的实现方式作限定。
步骤103:基于目标参数的目标取值,对感兴趣区域进行剪切波成像。
在本步骤中,将感兴趣区域的长度按照列数进行均分,并将感兴趣区域的宽度按照行数进行均分。例如,感兴趣区域的长度为10mm,宽度为20mm,行数为4,列数为2,此时,剪切波弹性成像设备将长度均分为2份,将宽度均分为4份。
接着,按照对检测区域进行激励与检测的预设顺序,依据激励参数的目标取值以及检测参数的目标取值,对每个检测区域进行激励,使得产生向检测区域传播的剪切波;以及按照检测参数的目标取值对检测区域进行检测,得到检测区域包含组织运动信息的超声回波数据。在得到超声脉冲检测波束反射回的剪切波经过感兴趣区域内各位置的超声回波数据后,依据超声脉冲检测波束反射的超声回波数据,确定剪切波经过感兴趣区域内的各位置时的传播速度。
具体的,依据感兴趣区域的超声回波数据,确定剪切波经过感兴趣区域中各位置点的传播速度的过程可以包括:对感兴趣区域的超声回波数据进行硬件或软件波束合成,得到波束合成后的RF或者IQ数据。利用基于RF数据的互相关或基于IQ数据的相移等算法得到轴向位移或轴向速度数据。在实际应用中,由于IQ数据相比于RF来说,具有在不损失原始信息的情况下,不仅数据量大大缩小,而且还含有相位信息的优点,因此,现有的大多数超声诊断设备利用IQ数据进行下一步处理。
若波束合成得到的数据为RF数据时,采用基于RF数据的互相关算法来确定轴向速度或轴向位移;若波束合成得到的数据为IQ数据时,采用基于IQ数据的组织过运动估算算法,常用的是Kasai的1D自相关算法和Loupas的2D自相关算法。
在人体中,由于其复杂的组织结构,因此在剪切波的传播过程中,当遇到声阻抗性质不一致的交界处时会产生发射等现象。由于反射等现象的存在,其会导致信号之间不相关,影响组织弹性的准确评估。因此,需要利用方向滤波器,将沿着同一个方向传播的剪切波进行滤出,同时将沿着其他方向传播的剪切波进行滤除,即将入射波和反射波分开。
然后再依据剪切波的传播特性,互相关和/或时间飞跃法等算法处理后得到剪切波经过感兴趣区域的各位置时的传播速度。
在得到剪切波经过感兴趣区域内各位置时的传播速度,即得到感兴趣区域对应的一个剪切波传播速度矩阵,一般情况下,剪切波传播速度矩阵中的各数值类型为浮点型。在本实施例中可以利用线性映射的方法,将每个浮点型数值转换成0-255范围内8位整型数值,得到整型剪切波传播速度矩阵;再利用灰度映射、伪彩或彩色算法将整型剪切波传播数据矩阵进行显示。
在实际应用中,通过剪切波经过感兴趣区域内各位置点时的传播速度,确定感兴趣组织的弹性的过程一般包括:将剪切波经过感兴趣区域内的各位置点的传播速度转换为剪切模量和杨式模量,再通过剪切模量与杨式模量分别与组织弹性间的关系,确定感兴趣组织的弹性。
其中,在理想条件下,可依据以下公式(1)公式(2)得到感兴趣区域内各位置处的剪切模量和杨氏模量:
G=ρ*vs 2 (1)
E=3*ρ*vs 2 (2)
式中,ρ表示组织密度,vs表示剪切波传播速度,G表示剪切模量,E表示杨氏模量。
本实施例中,将目标参数的当前取值变更为目标取值,体现了本申请实施例中的剪切波弹性成像设备中的目标参数的取值可调,因此,针对任意的感兴趣区域,都可以将目标参数的当前取值调整为与当前感兴趣区域的组织特性相匹配的目标取值,使得依据目标取值对感兴趣区域进行剪切波成像得到的成像结果可靠性更高。
图2为本申请中又一种剪切波弹性成像方法,可以包括:
步骤201:获取在多模式成像下的待检测组织中设定的感兴趣区域。
病变区域的组织弹性与病变区域的组织病理状态密切相关,因此病变区域组织弹性可以为疾病的诊断提供重要依据。如乳腺癌、甲状腺结节等疾病发生时,病变区域的组织弹性通常会发生改变。但是由于诊断对象的个体差异性和诊断项的不同,特别是病变区域复杂的组织成分及边界不清晰等特性,如果不能准确地得到病变区域的轮廓、大小、位置等信息,那么将会直接影响诊断的结果,因此,病变区域的轮廓、大小和位置等信息的确定准确性,直接影响诊断结果的准确性。
因此,在本步骤中,对待检测区域进行多模式成像。其中,待检测区域为待检测的组织区域。例如,需要确定甲状腺结节病变时,甲状腺结节就是待检测组织。多模式成像为包括B模式以及彩色血流模式的成像,在本步骤中,对待检测组织进行多模式成像,具体可以为对待检测组织进行B模式成像和/或彩色血流模式成像,得到B模式图像和/或彩色血流图像。
在本实施例中,用户依据多模式图像确定病变区域,以及感兴趣区域。由于B模式成像利用组织之间声阻抗的差异,在不同的组织界面上反射回不同强度的回波信号,从而可以得到描述不同组织结构的图像,因此,B模式成像反映的是待检测组织的结构信息。彩色血流模式成像利用多普勒效应,反映待检测组织的血流信息。因此,基于B模式图像和/或彩色血流模式图像,用户可以较准确地确定出待检测组织中病变区域的大小、位置以及轮廓等。
用户手动操作进入剪切波弹性成像模式,然后根据病变区域的范围设置感兴趣区域。在实际中,需要将感兴趣区域设置的比病变区域大,一方面是为了综合地分析病变区域内的弹性信息,另一方面则是为了包含一部分正常的组织,以此更好地进行正常与病变组织之间的对比分析。
步骤202:确定用于将感兴趣区域划分为检测区域的行数和列数。
用户在选定好感兴趣区域之后,剪切波弹性成像设备就依据目标参数的当前取值进行剪切波弹性成像。在本实施例中,对于目标参数中的行数和列数需要依据感兴趣区域的特点进行确定,原因如下:
剪切波弹性成像设备通过发射高能量的超声脉冲于组织中,使感兴趣区域内产生剪切波。产生声辐射力的方法是发射高能量的超声脉冲,该超声脉冲可达几百至几千个周期,即可达几十至几百微秒的时长,进而产生声辐射力。如图3所示,为本申请中一种超声声辐射力的产生示意图,由图3中的5超声探头发射聚焦长脉冲1,聚焦于某一深度2,该深度一般位于感兴趣区域的中心位置。通过动量原理,即超声声辐射力与组织相互作用,该超声声辐射力被组织大量的吸收,使组织内部局部发生振动而产生剪切形变,从而在感兴趣区域内产生剪切波3。其中,发射探头可以是采用和剪切波检测同一个探头,也可以是采用和剪切波检测不同的探头。然后发射剪切波跟踪脉冲4跟踪剪切波的传播。
对于传统的单焦点激励方法来说,超声声辐射力的计算公式如下公式(3)所示:
其中,F是超声声辐射力,Wabsorbed是给定空间位置处所吸收的能量,c为常数。也即,当用于激励的脉冲时长变大时,其产生的声辐射力则越大,那么所产生的剪切波的信号强度就越大。但由于组织的衰减的影响,因此剪切波在组织中的有效传播距离是有限的。
然后,在剪切波的传播路径上,重复多次发射覆盖范围较大的剪切波检测超声脉冲波束。这一类超声波束可以是非聚集单角度或多角度平面超声波束,也可以是聚焦超声波束,也可以是宽超声波束,从而采集感兴趣区域内针对多次发射的超声脉冲检测波束反射回来的包含剪切波传播信息的回波数据。
如图4所示是传统单焦点激励的发射声场示意图。其中,图4中的6表示聚焦深度,图4中的7表示用于发射的孔径大小,图4中的8表示波束宽度,图4中的9表示聚焦范围。由图4中的6和7可以定义发射的F值,如下公式(4)式所示:
其中,Z就表示为聚焦深度,D表示为孔径的大小。可以得出,在一定的聚焦深度下,F值和孔径大小呈反比,即F值越小,用于超声声辐射力激励的孔径就越大,则超声声辐射力就越强;反之则超声声辐射力越弱。而图4中的8则可由如下公式(5)决定:
BeamWidth=1.4*F*wavelength (5)
其中,BeamWidth是波束宽度,wavelength是发射超声波的波长。由于在此波束宽度范围内无剪切波的传播,考虑到波束宽度是以焦点左右对称的,因此检测剪切波传播波束的位置应至少远离一半的波束宽度。而对于图4中的聚焦范围,则可由如下公式(6)决定:
DepthOFiel d=7.1*F*F*wavelength s (6)
其中,DepthOfField是指聚焦范围。由于传统单点声辐射力激励所产生的聚焦范围是有限的,也就是说,在该聚焦范围内,所采集的包含剪切波运动的数据是信噪比较高的。
综上可得,由于剪切波在感兴趣区域内传播的过程中,剪切波的信号强度会产生衰减,在现有技术中,通过单点激励在感兴趣区域内所产生的剪切波信号强度随着传播距离的增大迅速衰减,并且单点激励所产生的聚焦范围有限,使得聚焦范围外的剪切波信号强度发生衰减。当剪切波信号强度衰减到可检测强度范围外时,超声脉冲检测波束所反射的超声回波数据的信噪比较低,使得依据超声回波数据所确定出的剪切波传播速度的准确性降低。因此,本实施例中,剪切波弹性成像设备需要将感兴趣区域划分为多个子区域(检测区域)。具体的,将感兴趣区域划分为多个检测区域,需要确定将感兴趣区域划分为几行和几列,行列所组成的每个区域就是检测区域。
具体的,本实施例中,确定用于将感兴趣区域划分为检测区域的行数和列数的过程可以包括步骤A1~步骤A4:
步骤A1:获取感兴趣区域在剪切波传播方向上的长度,以及感兴趣区域在轴向方向上的宽度。
在实际应用中,当获取到感兴趣区域后,就可以按照传统产生剪切波的方式,确定出后续在感兴趣区域内传播的剪切波的传播方向。在本步骤中,轴向方向为平行于感兴趣区域的表面且垂直于传播方向的方向。
具体的,用户确定出待检测组织中的感兴趣区域后,一般通过一个矩形标注出感兴趣区域,因此,在本步骤中,基于用户标注的矩形,确定矩形在传播方向上的长度以及在轴向方向上的宽度。
需要说明的是,本步骤只是给出了一种用户对感兴趣区域的具体标注形式,当然,在实际应用中,用户还可以通过其他形状对感兴趣区域进行标注,本实施例不对具体的标注方式作限定。只要可以识别出用户所标注出的感兴趣区域,以及确定出感兴趣区域在传播方向上的长度,以及在轴向方向上的宽度即可。
步骤A2:依据剪切波在传播方向上的预设传播距离与感兴趣区域在传播方向上的长度,确定将感兴趣区域在传播方向上划分成的列数,使得所划分得到的每列长度与预设传播距离间的差值小于预设第一阈值。
在本实施例中,预设传播距离为单点激励中剪切波在组织中的有效传播距离。在本步骤中,依据预设传播距离与所确定出的感兴趣区域在传播方向上的长度,将感兴趣区域在传播方向上划分为多列,其中,所划分得到的每列的长度与预设传播距离间的差值不大于预设第一阈值。在本实施例中,为了描述方便,将感兴趣区域在传播方向上的长度所划分成的线段数量,统称为列数。
例如,假设单点激励所产生的剪切波在组织中的有效传播距离为7mm,即剪切波的预设传播距离为7mm;并且,感兴趣区域在传播方向上的长度为14mm,此时,在本步骤中,将感兴趣区域在传播方向上的长度划分为2份,即列数为2。
步骤A3:依据剪切波在轴向方向上的预设聚焦范围与确定出感兴趣区域在轴向方向上的宽度,确定将感兴趣区域在轴向方向上划分成的行数,使得划分得到的每行宽度与预设聚焦范围间的差值小于预设第二阈值。
在本步骤中,依据预设聚焦范围与所确定出的感兴趣区域在轴向方向上的宽度,将感兴趣区域在轴向方向上划分为多段,其中,所划分得到的每段的宽度与预设聚焦范围间的差值不大于预设第二阈值。在本实施例中,为了描述方便,将感兴趣区域在轴向方向上的宽度所划分成的线段数量,统称为行数。
例如,假设单点激励所产生的剪切波的聚焦范围为10mm,感兴趣区域在轴向方向上的宽度为20mm,此时,在本步骤中,将感兴趣区域在轴向方向上的宽度划分为两,即行数为2。
步骤A4:依据行数与列数,将感兴趣区域划分为至少两个初始子区域。
还以单点激励所产生的剪切波聚焦范围为10mm,预设传播距离为7mm,感兴趣区域的长度为14mm,宽度为20mm为例,得到将感兴趣区域在传播方向上划分为两段,在轴向方向上也划分为两段;在本步骤中,将行数与列数的乘积,确定为将感兴趣区域划分成子区域的数量,即将感兴趣区域划分为4个子区域。在本实施例中,为了描述方便,将感兴趣区域所划分后的每个子区域,统称为检测区域。
需要说明的是,根据上述步骤A1~步骤A4所划分得到的每个检测区域在传播方向上的长度与预设距离间的差值为0,该差值是小于预设第一阈值,并且,所划分得到的每个检测区域在轴向方向上的宽度与预设聚焦范围间的差值为0,该差值同样小于预设第二阈值。本实施例只是给出了一种划分方式,在实际应用中,还可以采用其他划分方式对感兴趣区域进行划分,本实施例并不对具体的划分方式作限定,只要所划分得到的每个检测区域的长度与预设传播距离间的差值小于预设第一阈值,每个检测区域的宽度与预设聚焦范围间的差值小于预设第二阈值即可。
上述步骤A1~步骤A4确定出用于将感兴趣区域划分为检测区域的行数和列数;对于目标参数中的除了行数和列数之外的其他参数的目标取值可以为默认值。
步骤203:将目标参数的当前取值变更为目标取值。
具体过程可以参考步骤102,这里不再赘述。
步骤204:基于目标参数的目标取值,对感兴趣区域进行剪切波弹性成像。
具体过程,可以参考步骤103,这里不再赘述。
在本申请实施例中,将感兴趣区域划分成的每个检测区域在水平方向上的长度与预设传播距离间的差值小于预设第一阈值,并且,在轴向方向上的宽度与预设聚焦范围间的差值小于预设第二阈值;因此,本申请实施例中分别对每个检测区域对应的激励位置进行激励与检测,此时,超声脉冲检测波束所反射回的超声回波数据的信噪比得到提高,此时,依据该位置点对应的超声回波数据所确定出的剪切波传播数据的准确性较高;因此,本申请实施例将感兴趣区域划分为多个检测区域,并将每个检测区域进行激励与检测的方式,可以提高所确定出感兴趣区域中各位置处的剪切波传播速度的可靠性,进而,提高所确定出的组织弹性。
参考图5示出了本申请中又一种剪切波弹性成像方法,可以包括:
步骤501:获取在多模式成像下的待检测组织中设定的感兴趣区域。
具体过程可以参考步骤201,这里不再赘述。
步骤502:接收用户对目标参数的第一调整数值。
在用户选定好感兴趣区域后,剪切波弹性成像设备进行剪切波弹性成像前,在本步骤中,接收用户对目标参数的调整数值,为了描述方便,将接收用户对目标参数的调整数值称为第一调整数值。
具体的,在实际应用中,用户对目标参数调整的一种场景可以为:用户输入将感兴趣区域划分为目标区域的行数和列数(如4*16,其中,4表示行数,16表示列数);在实际应用中,用户依据增加检测区域的个数可以保证每个检测区域的剪切波信号的信噪比高的原则,对行数与列数进行调整,从而可以针对病灶区域的结果进行增强处理,得到更加可靠的成像结果。
用户还可以输入对检测区域进行激励所需的激励参数,对检测区域进行检测的检测参数;具体的,用户可以依据以下原则,对激励参数与检测参数进行调整:在需要提高激励能量的情况下,增加激励脉冲周期数和/或降低激励发射F值;激励发射频率的取值范围可以为以探头中心频率的-6dB带宽的下边界为中心的第一预设范围。在需要增加超声波的穿透深度的情况下,增加检测脉冲周期数;检测发射频率的取值范围可以为以探头中心频率为中心的第二预设范围;
在实际应用中,用户还可以按照划分后的检测区域在感兴趣区域中的分布,输入对检测区域进行激励与检测的先后顺序。其中,顺序的具体形式,可以根据实际的场景进行设定,本实施例不对顺序的具体形式作限定。
步骤503:将目标参数的当前取值变更为目标取值。
具体过程可以参考步骤103,这里不再赘述。
步骤504:基于目标参数的目标取值,对感兴趣区域进行剪切波弹性成像。
本实施例给出了一个剪切波弹性成像设备依据目标参数的目标取值,将感兴趣区域划分成目标检测区域的示意图。具体的,如图6所示,为本申请中一种将感兴趣区域划分为9个检测区域的结果示意图。
目标参数中将感兴趣区域划分为目标检测区域的行数和列数都为3,即将感兴趣区域划分为9个目标检测区域。在该图6中,将感兴趣区域划分为18、20、22、24、26、28、30、32和34这9个检测区域,并且,19表示对应检测区域18的激励位置,21表示对应检测区域20的激励位置,23表示对应检测区域22的激励位置,25表示对应检测区域24的激励位置,27表示对应检测区域26的激励位置,29表示对应检测区域28的激励位置,31表示对应检测区域30的激励位置,33表示对应检测区域32的激励位置,35表示对应检测区域34的激励位置。
具体的,以上述将感兴趣区域划分为9个检测区域为例,对每个检测区域进行激励与检测的顺序依次为:18、20、22、24、26、28、30、32和34。对于任意一个检测区域,以标号为18的检测区域为例,在标号18的检测区域对应的激励位置19上,按照标号18的检测区域对应的激励参数进行激励,产生向标号18的检测区域传播的剪切波,接着,按照用户对检测区域设置的检测参数,对标号18的检测区域进行剪切波检测,得到超声脉冲检测波束发射的标号18的检测区域的超声回波数据。接着,按照相同的激励与检测逻辑,对标号20的检测区域进行激励与检测,依次循环,使得对每个检测区域都产生剪切波以及得到检测波束反射的包含剪切波传播信息的超声回波数据。
需要说明的是,本实施例只是给出了一种对多个检测区域进行激励与检测的顺序,在实际应用中,用户在设定顺序时不一定要按照18、20、22、24、26、28、30、32和34的顺序设定。
在本实施例中,剪切波弹性成像设备为用户设置目标参数的目标取值提供了可能,使得剪切波弹性成像设备可以按照用户确定的目标参数的目标取值,对感兴趣区域进行剪切波弹性成像。
在图2和图5所对应的实施例中,在依据目标参数的目标取值,对感兴趣区域进行剪切波弹性成像之后,还可以包括:
在依据目标参数的目标取值,对感兴趣区域进行剪切波弹性成像之后,用户可能觉得当前的成像结果与预估结果(依据多模式成像对剪切波弹性成像的预估结果)有差距。在这种情况下,用户可以对目标参数进行调整,具体的,用户向剪切波弹性成像设备发送调整指令,此时,剪切波弹性成像设备显示用于用户调整目标参数的界面,用户可以通过该界面对目标参数进行调整,此时,剪切波弹性成像设备接收用户调整后的目标参数的目标取值,为了描述方便,将本次用户对目标参数的调整值称为第二调整数值。具体调整过程可以参考步骤502,这里不再赘述。
在实际应用中,为了方便用户对目标参数的调整,剪切波弹性成像设备在接收到用户的调整指令后,在显示界面上显示目标参数的当前取值,以便用户依据目标参数的当前取值确定目标参数的目标取值。
具体的,显示将感兴趣区域所划分成的子区域分布图,以及每个检测区域对应的激励参数与检测参数。如图7所示,为本申请中一种将感兴趣区域划分为4个检测区域的结果示意图。在该图7中,示出了将感兴趣区域划分为4个大小相同的检测区域的分布情况,以及,针对每个检测区域的激励位置。其中,10、12、14和16分别对应的区域为划分得到的检测区域,11表示对应检测区域10的激励位置,13表示对应检测区域12的激励位置,15表示对应检测区域14的激励位置,17则表示对应检测区域16的激励位置。
其中,对于任一激励位置的宽度表示:用于在激励位置激励剪切波的激励波束的宽度,并且,激励位置与检测区域间的距离不小于激励波束的宽度的一半。
在本实施例中,显示出将感兴趣区域划分为的检测区域,以及每个检测区域的激励位置、激励参数与检测参数后,用户可以根据感兴趣区域内的组织的具体情况确定目标参数的目标取值。
图8为本申请中一种剪切波弹性成像装置实施例,可以包括:
获取单元801,用于获取目标参数的目标取值;所述目标参数为用于对感兴趣区域进行剪切波成像的参数;所述目标取值为与所述感兴趣区域的组织特性相匹配的取值;
第一变更单元802,用于将所述目标参数的当前取值变更为所述目标取值;
第一成像单元803,用于基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像。
其中,所述获取单元所获取的目标参数包括:用于将所述感兴趣区域划分为检测区域的行数和列数;
所述获取单元801,可以包括:
第一确定子单元,用于确定所述感兴趣区域的长度和宽度;所述长度为所述感兴趣区域在剪切波的传播方向上的尺度;所述宽度为所述感兴趣区域在轴向方向上的尺度;所述轴向方向为平行于所述感兴趣区域且垂直于所述传播方向的方向;
第二确定子单元,用于依据剪切波在传播方向上的预设传播距离,以及在轴向方向上的预设聚焦范围,确定将所述感兴趣区域划分为满足预设条件的检测区域的行数和列数;所述预设条件包括:在所述传播方向上的长度与所述预设传播距离间的差值小于预设第一阈值,且在所述轴向方向上的宽度与所述预设聚焦范围间的差值小于预设第二阈值。
其中,所述获取单元801还可以包括:
接收子单元,用于接收用户对所述目标参数的第一调整数值;
第三确定子单元,用于确定所述第一调整数值为所述目标参数的目标取值。
其中,该装置实施例还包括:
接收单元,用于在所述第一成像单元基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像之后,在接收到调整指令的情况下,接收用户对所述目标参数的第二调整数值;
第二变更单元,用于将所述目标参数的当前取值变更为所述第二调整数值;
第二成像单元,用于基于所述第二调整数值,对所述感兴趣区域进行剪切波弹性成像。
其中,该装置实施例还可以包括:
第一显示单元,用于在所述接收单元接收用户对所述目标参数的第二调整数值之前,显示所述目标参数的当前取值。
其中,该装置实施例还可以包括:
多模式成像单元,用于对待检测组织进行多模式成像,得到所述待检测组织的多模式图像;
第二显示单元,用于显示所述多模式图像,使得用户依据所述多模式图像确定所述感兴趣区域。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其它实施例的不同之处,各个实施例之间相同相似部分互相参见即可。在本文中,诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。在文中的“包括”、“包含”等词语解释为包含的含义而不是排他或穷举的含义;也就是说,是“包含但不限于”的含义。在不脱离本发明构思的前提下,还可以做出变形、同等替换、改进等,这些都属于本发明的保护范围。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下,在其它实施例中实现。因此,本申请将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (12)
1.一种剪切波弹性成像方法,其特征在于,包括:
获取目标参数的目标取值;所述目标参数为用于对感兴趣区域进行剪切波成像的参数;所述目标取值为与所述感兴趣区域的组织特性相匹配的取值;
将所述目标参数的当前取值变更为所述目标取值;
基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像。
2.根据权利要求1所述的方法,其特征在于,所述目标参数包括:用于将所述感兴趣区域划分为检测区域的行数和列数;
所述获取目标参数的目标取值,包括:
确定所述感兴趣区域的长度和宽度;所述长度为所述感兴趣区域在剪切波的传播方向上的尺度;所述宽度为所述感兴趣区域在轴向方向上的尺度;所述轴向方向为平行于所述感兴趣区域且垂直于所述传播方向的方向;
依据剪切波在传播方向上的预设传播距离,以及在轴向方向上的预设聚焦范围,确定将所述感兴趣区域划分为满足预设条件的检测区域的行数和列数;所述预设条件包括:在所述传播方向上的长度与所述预设传播距离间的差值小于预设第一阈值,且在所述轴向方向上的宽度与所述预设聚焦范围间的差值小于预设第二阈值。
3.根据权利要求1所述的方法,其特征在于,所述获取目标参数的目标取值,包括:
接收用户对所述目标参数的第一调整数值;
确定所述第一调整数值为所述目标参数的目标取值。
4.根据权利要求1~3任意一项所述的方法,其特征在于,在所述基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像之后,还包括:
在接收到调整指令的情况下,接收用户对所述目标参数的第二调整数值;
将所述目标参数的当前取值变更为所述第二调整数值;
基于所述第二调整数值,对所述感兴趣区域进行剪切波弹性成像。
5.根据权利要求4所述的方法,其特征在于,在所述接收用户对所述目标参数的第二调整数值之前,还包括:
显示所述目标参数的当前取值。
6.根据权利要求1所述的方法,其特征在于,在所述获取目标参数的目标取值之前,还包括:
对待检测组织进行多模式成像,得到所述待检测组织的多模式图像;
显示所述多模式图像,使得用户依据所述多模式图像确定所述感兴趣区域。
7.一种剪切波弹性成像装置,其特征在于,包括:
获取单元,用于获取目标参数的目标取值;所述目标参数为用于对感兴趣区域进行剪切波成像的参数;所述目标取值为与所述感兴趣区域的组织特性相匹配的取值;
第一变更单元,用于将所述目标参数的当前取值变更为所述目标取值;
第一成像单元,用于基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像。
8.根据权利要求7所述的装置,其特征在于,所述获取单元所获取的目标参数包括:用于将所述感兴趣区域划分为检测区域的行数和列数;
所述获取单元,包括:
第一确定子单元,用于确定所述感兴趣区域的长度和宽度;所述长度为所述感兴趣区域在剪切波的传播方向上的尺度;所述宽度为所述感兴趣区域在轴向方向上的尺度;所述轴向方向为平行于所述感兴趣区域且垂直于所述传播方向的方向;
第二确定子单元,用于依据剪切波在传播方向上的预设传播距离,以及在轴向方向上的预设聚焦范围,确定将所述感兴趣区域划分为满足预设条件的检测区域的行数和列数;所述预设条件包括:在所述传播方向上的长度与所述预设传播距离间的差值小于预设第一阈值,且在所述轴向方向上的宽度与所述预设聚焦范围间的差值小于预设第二阈值。
9.根据权利要求7所述的装置,其特征在于,所述获取单元还包括:
接收子单元,用于接收用户对所述目标参数的第一调整数值;
第三确定子单元,用于确定所述第一调整数值为所述目标参数的目标取值。
10.根据权利要求7~9任意一项所述的装置,其特征在于,还包括:
接收单元,用于在所述第一成像单元基于所述目标参数的目标取值,对所述感兴趣区域进行剪切波弹性成像之后,在接收到调整指令的情况下,接收用户对所述目标参数的第二调整数值;
第二变更单元,用于将所述目标参数的当前取值变更为所述第二调整数值;
第二成像单元,用于基于所述第二调整数值,对所述感兴趣区域进行剪切波弹性成像。
11.根据权利要求10所述的装置,其特征在于,还包括:
第一显示单元,用于在所述接收单元接收用户对所述目标参数的第二调整数值之前,显示所述目标参数的当前取值。
12.根据权利要求7所述的装置,其特征在于,还包括:
多模式成像单元,用于对待检测组织进行多模式成像,得到所述待检测组织的多模式图像;
第二显示单元,用于显示所述多模式图像,使得用户依据所述多模式图像确定所述感兴趣区域。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811063847.9A CN109259801B (zh) | 2018-09-12 | 2018-09-12 | 一种剪切波弹性成像方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811063847.9A CN109259801B (zh) | 2018-09-12 | 2018-09-12 | 一种剪切波弹性成像方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109259801A true CN109259801A (zh) | 2019-01-25 |
| CN109259801B CN109259801B (zh) | 2021-05-14 |
Family
ID=65188801
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811063847.9A Active CN109259801B (zh) | 2018-09-12 | 2018-09-12 | 一种剪切波弹性成像方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109259801B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111671464A (zh) * | 2020-06-19 | 2020-09-18 | 深圳开立生物医疗科技股份有限公司 | 一种图像标记方法、装置、设备及可读存储介质 |
| CN113261991A (zh) * | 2019-09-27 | 2021-08-17 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性成像方法、系统及计算机可读存储介质 |
| WO2022147690A1 (zh) * | 2021-01-06 | 2022-07-14 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性成像方法和超声成像系统 |
| WO2022147689A1 (zh) * | 2021-01-06 | 2022-07-14 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性测量方法、装置和存储介质 |
| CN119818085A (zh) * | 2023-10-12 | 2025-04-15 | 青岛海信医疗设备股份有限公司 | 超声图像的质量评估方法及超声设备 |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6508768B1 (en) * | 2000-11-22 | 2003-01-21 | University Of Kansas Medical Center | Ultrasonic elasticity imaging |
| US20040034304A1 (en) * | 2001-12-21 | 2004-02-19 | Chikayoshi Sumi | Displacement measurement method and apparatus, strain measurement method and apparatus elasticity and visco-elasticity constants measurement apparatus, and the elasticity and visco-elasticity constants measurement apparatus-based treatment apparatus |
| CN1875881A (zh) * | 2005-06-07 | 2006-12-13 | 西门子公司 | 拍摄、分析和显示医学图像数据组的方法 |
| CN102081697A (zh) * | 2009-11-27 | 2011-06-01 | 深圳迈瑞生物医疗电子股份有限公司 | 一种在超声成像空间中定义感兴趣容积的方法及其装置 |
| CN102283679A (zh) * | 2011-08-04 | 2011-12-21 | 中国科学院深圳先进技术研究院 | 弹性测量的超声成像系统及测量生物组织弹性的方法 |
| US8175350B2 (en) * | 2007-01-15 | 2012-05-08 | Eigen, Inc. | Method for tissue culture extraction |
| CN102940510A (zh) * | 2012-08-31 | 2013-02-27 | 华南理工大学 | 一种超声弹性成像的自动对焦方法 |
| CN103156636A (zh) * | 2011-12-15 | 2013-06-19 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像装置和方法 |
| US20130329030A1 (en) * | 2007-03-08 | 2013-12-12 | Sync-Rx, Ltd. | Apparatus and methods for masking a portion of a moving image stream |
| CN104688266A (zh) * | 2013-12-04 | 2015-06-10 | 美国西门子医疗解决公司 | 三维弹性超声成像中的运动修正 |
| CN105844645A (zh) * | 2016-03-31 | 2016-08-10 | 深圳大学 | 一种瞬时弹性成像编码检测方法及系统 |
| CN107049360A (zh) * | 2017-01-26 | 2017-08-18 | 清华大学 | 剪切波弹性成像方法和装置 |
| CN107970043A (zh) * | 2017-12-28 | 2018-05-01 | 深圳开立生物医疗科技股份有限公司 | 一种剪切波的检测方法及装置 |
-
2018
- 2018-09-12 CN CN201811063847.9A patent/CN109259801B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6508768B1 (en) * | 2000-11-22 | 2003-01-21 | University Of Kansas Medical Center | Ultrasonic elasticity imaging |
| US20040034304A1 (en) * | 2001-12-21 | 2004-02-19 | Chikayoshi Sumi | Displacement measurement method and apparatus, strain measurement method and apparatus elasticity and visco-elasticity constants measurement apparatus, and the elasticity and visco-elasticity constants measurement apparatus-based treatment apparatus |
| CN1875881A (zh) * | 2005-06-07 | 2006-12-13 | 西门子公司 | 拍摄、分析和显示医学图像数据组的方法 |
| US8175350B2 (en) * | 2007-01-15 | 2012-05-08 | Eigen, Inc. | Method for tissue culture extraction |
| US20130329030A1 (en) * | 2007-03-08 | 2013-12-12 | Sync-Rx, Ltd. | Apparatus and methods for masking a portion of a moving image stream |
| CN102081697A (zh) * | 2009-11-27 | 2011-06-01 | 深圳迈瑞生物医疗电子股份有限公司 | 一种在超声成像空间中定义感兴趣容积的方法及其装置 |
| CN102283679A (zh) * | 2011-08-04 | 2011-12-21 | 中国科学院深圳先进技术研究院 | 弹性测量的超声成像系统及测量生物组织弹性的方法 |
| CN103156636A (zh) * | 2011-12-15 | 2013-06-19 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像装置和方法 |
| CN102940510A (zh) * | 2012-08-31 | 2013-02-27 | 华南理工大学 | 一种超声弹性成像的自动对焦方法 |
| CN104688266A (zh) * | 2013-12-04 | 2015-06-10 | 美国西门子医疗解决公司 | 三维弹性超声成像中的运动修正 |
| CN105844645A (zh) * | 2016-03-31 | 2016-08-10 | 深圳大学 | 一种瞬时弹性成像编码检测方法及系统 |
| CN107049360A (zh) * | 2017-01-26 | 2017-08-18 | 清华大学 | 剪切波弹性成像方法和装置 |
| CN107970043A (zh) * | 2017-12-28 | 2018-05-01 | 深圳开立生物医疗科技股份有限公司 | 一种剪切波的检测方法及装置 |
Non-Patent Citations (4)
| Title |
|---|
| ANDREW EVANS, PATSY WHELEHAN, KIM THOMSON,ET AL: "Quantitative shear wave ultrasound elastography:initial experience in solid breast masses", 《BREAST CANCER RESEARCH》 * |
| GEORGIANA NAGY, ET AL: "Computerized ultrasound image analysis for noninvasive evaluation of hepatic steatosis", 《MED ULTRASON》 * |
| 王丽婷 姜翔飞 张嵩 王丛知 郑海荣: "超声声辐射力弹性成像中基于拉东变换的剪切波速度估计方法的改进研究", 《集成技术》 * |
| 黄燕平,郑永平: "基于超声的组织弹性测量之剪切波传播法", 《专论》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113261991A (zh) * | 2019-09-27 | 2021-08-17 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性成像方法、系统及计算机可读存储介质 |
| CN111671464A (zh) * | 2020-06-19 | 2020-09-18 | 深圳开立生物医疗科技股份有限公司 | 一种图像标记方法、装置、设备及可读存储介质 |
| WO2022147690A1 (zh) * | 2021-01-06 | 2022-07-14 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性成像方法和超声成像系统 |
| WO2022147689A1 (zh) * | 2021-01-06 | 2022-07-14 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性测量方法、装置和存储介质 |
| CN116133597A (zh) * | 2021-01-06 | 2023-05-16 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性测量方法、装置和存储介质 |
| CN116528771A (zh) * | 2021-01-06 | 2023-08-01 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性成像方法和超声成像系统 |
| CN119818085A (zh) * | 2023-10-12 | 2025-04-15 | 青岛海信医疗设备股份有限公司 | 超声图像的质量评估方法及超声设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109259801B (zh) | 2021-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2473113B1 (en) | Method and apparatus for ultrasonic detection and imaging of hemodynamic information, particularly venous blood flow information | |
| KR101586998B1 (ko) | 의료용 초음파 이미징에서의 전단파 정보의 해석 | |
| US10338203B2 (en) | Classification preprocessing in medical ultrasound shear wave imaging | |
| CN109259801A (zh) | 一种剪切波弹性成像方法及装置 | |
| EP3040033A1 (en) | Diagnostic ultrasound apparatus and elasticity evaluation method | |
| JP2010526626A (ja) | 関心領域の粘弾性の平均値を測定するための方法および装置 | |
| US20130211253A1 (en) | On-axis Shear Wave Characterization with Ultrasound | |
| US20150094579A1 (en) | Shear Wave Detection in Medical Ultrasound Imaging | |
| US10548572B2 (en) | Ultrasound processing device | |
| CN109717899A (zh) | 超声医学成像中根据剪切速率的组织粘弹性估计 | |
| KR20180013956A (ko) | 단일 추적 위치 전단파 탄성 이미징을 위한 방법, 시스템 및 컴퓨터 프로그램 제품 | |
| CN108784744A (zh) | 用于剪切波成像的可变焦点 | |
| JP2006115937A (ja) | 超音波診断装置 | |
| CN109512465A (zh) | 一种声辐射力双向剪切波复合成像方法及装置 | |
| US11963824B2 (en) | Shear wave imaging based on ultrasound with increased pulse repetition interval | |
| US11446007B2 (en) | Ultrasonic shearwave imaging with patient-adaptive shearwave generation | |
| CN108852416A (zh) | 一种剪切波传播速度的确定方法及装置 | |
| Urban et al. | Harmonic pulsed excitation and motion detection of a vibrating reflective target | |
| US20240350115A1 (en) | Method and apparatus for backscattering imaging vascular activity at a microscopic scale | |
| EP2853918B1 (en) | Shear wave detection in medical ultrasound imaging | |
| JP5491671B2 (ja) | 関心領域の粘弾性の平均値を測定するための方法 | |
| Chen | BASICS OF ULTRASOUND AND PHOTOACOUSTIC IMAGING | |
| JP4060412B2 (ja) | 超音波診断装置 | |
| EP1898794A2 (en) | Method and apparatus for use with doppler measurements in medical applications |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |