CN110368031A - 弹性测量检测方法及系统 - Google Patents
弹性测量检测方法及系统 Download PDFInfo
- Publication number
- CN110368031A CN110368031A CN201910741685.8A CN201910741685A CN110368031A CN 110368031 A CN110368031 A CN 110368031A CN 201910741685 A CN201910741685 A CN 201910741685A CN 110368031 A CN110368031 A CN 110368031A
- Authority
- CN
- China
- Prior art keywords
- target area
- shearing wave
- time
- data
- propagation distance
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52019—Details of transmitters
- G01S7/5202—Details of transmitters for pulse systems
- G01S7/52022—Details of transmitters for pulse systems using a sequence of pulses, at least one pulse manipulating the transmissivity or reflexivity of the medium
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明提供了一种弹性测量检测方法,包括以下步骤,于目标区域产生剪切波;向目标区域发射追踪脉冲并接收所述追踪脉冲的回波数据;根据所述回波数据进行数据计算以获取测量计算结果;显示所述测量计算结果。本发明的弹性测量检测方法不需要操作者控制按压探头的力度即可对反映组织弹性的相关参数进行计算获取,使得测量的重复性和稳定性得到保证。本发明还公开一种弹性测量检测系统。
Description
技术领域
本发明涉及医疗检测领域,尤其涉及一种弹性测量检测方法及系统。
背景技术
超声弹性成像用于反映组织(Tissue)的弹性或软硬程度,在组织癌症病变的辅助检测、良恶性判别、预后恢复评价等方面得到越来越多应用。常规超声弹性成像方法主要通过探头按压组织产生一定的形变,再将应变量、应变率等与组织弹性相关的参数计算出来并成像,间接反映感兴趣区域中不同组织之间的弹性差异。但是,由于应变参数与弹性参数之间的对应关系受到压力大小的影响,即使是相同的组织,或者具有相同的弹性参数的不同组织,如果施加的压力不同,产生的应变也不同,可能造成对弹性差异的误判。而操作者使用时需要人为控制每次操作的压力都保持一致,且维持均匀稳定,这使得图像的重复性和稳定性得到挑战。
发明内容
提供一种不受操作要求限制,具有稳定测量结果的弹性测量检测方法及系统。
一种弹性测量检测方法,包括以下步骤,
于目标区域产生剪切波;
向所述目标区域发射追踪脉冲并接收所述追踪脉冲的回波数据;
根据所述回波数据进行数据计算以获取目标区域内的剪切波传播距离;以及
显示所述剪切波传播距离。
进一步的,根据所述回波数据进行数据计算以获取目标区域内的剪切波传播距离,包括:获取所述目标区域内各位置处一段时间内的剪切波传播距离;
显示所述剪切波传播距离包括:
分别显示所述目标区域内各位置处一段时间内的剪切波传播距离,形成剪切波传播距离分布图;或者,根据所述回波数据获得目标区域的超声图像,在所述超声图像上叠加显示所述目标区域内各位置处一段时间内的剪切波传播距离,形成剪切波传播距离分布图;或者,根据所述目标区域内各位置处一段时间内的剪切波传播距离,显示所述目标区域内的各位置处的平均传播距离。
进一步的,获取所述目标区域内各位置处一段时间内的剪切波传播距离,包括:
确定剪切波经过所述目标区域中各位置的任一位置的时刻;
计算剪切波经过该位置的时刻附近的一段时间内的剪切波传播距离。
进一步的,确定剪切波经过所述目标区域中各位置的任一位置的时刻,包括:
根据所述回波数据进行数据计算获得目标区域的位移数据或质点速度数据;
取目标区域中任一位置对应的各个时刻的位移数据或质点速度数据,进而得到剪切波经过该位置的时刻。
进一步的,采用灰阶或彩色编码模式,显示各位置处一段时间内的所述剪切波传播距离和平均传播距离。
进一步的,根据所述回波数据进行数据计算以获取目标区域内的剪切波传播距离,包括:
取多个时刻;
分别获取剪切波于所述目标区域中在各个时刻附近的一段时间内的剪切波传播距离。
进一步的,显示所述剪切波传播距离包括:以曲线形式显示各个时刻附近的一段时间内的剪切波传播距离。
进一步的,显示所述剪切波传播距离包括:根据剪切波于所述目标区域内在各个时刻附近的一段时间内的剪切波传播距离,显示剪切波传播过程中一段时间内的平均传播距离。
进一步的,还包括:累积多个时刻附近的一段时间的剪切波传播距离,得到各个时刻目标区域内的剪切波的总传播距离;
显示所述剪切波传播距离包括:以曲线形式显示各个时刻目标区域内的剪切波的总传播距离。
进一步的,所述经过任一位置的时刻附近的一段时间包括经过任一位置的时刻的向前一段时间、向后一段时间或前后一段时间;所述各个时刻附近的一段时间包括所述各个时刻向前的一段时间、各个时刻向后的一段时间、或各个时刻前后的一段时间。
进一步的,所述弹性测量检测方法于目标区域产生剪切波前还包括以下步骤,
向目标区域发射参考脉冲,接收所述参考脉冲于所述目标区域内的回波信息,获取所述参考脉冲的回波数据。
进一步的,当根据所述回波数据进行数据计算以获取测量计算结果时,通过数据计算获取位移数据,当获取位移数据时包括以下步骤,
取目标区域内任一位置的任一时刻的回波数据作为基准回波数据;
将基准回波数据分为多段基准回波核数据,设置所述基准回波核数据的基准回波核数据中心;
获取各个时刻的回波数据上与所述基准回波核数据中心之间互相关性最大的位置,取该互相关性最大的位置与所述基准回波核数据的中心位置之差为各个时刻的位移数据。
进一步的,当取目标区域内任一位置的任一时刻的回波数据作为基准回波数据时,所述基准回波数据采用任一时刻的追踪脉冲的回波数据或任一时刻的参考脉冲的回波数据。
进一步的,当根据所述回波数据进行数据计算以获取测量计算结果时,通过数据计算获取不同时刻的质点速度,当获取质点速度时包括以下步骤,
获取目标区域的组织的位移数据后,取同一个位置不同时刻的位移数据;
沿时间方向求梯度以获取该位置处不同时刻的质点速度曲线。
进一步的,当根据所述回波数据进行数据计算以获取测量计算结果时,通过数据计算获取一段时间段内的剪切波的传播距离,当获取所述传播距离进行计算时包括以下步骤,
取目标区域中任一位置对应的各个时刻的位移数据,形成位移-时间曲线,找出该曲线上峰值对应的时刻;
获取目标区域内的所述峰值对应的时刻前后一段时间段内的位移-横向位置曲线或者质点速度-横向位置曲线,对所述时间段内各个时刻的位移-横向位置曲线或者质点速度-横向位置曲线进行互相关判断,获取所述时间段内的剪切波的传播距离。
进一步的,当根据所述回波数据进行数据计算以获取测量计算结果时,获取不同目标区域之间的杨氏模量比值,当获取所述杨氏模量比值时,满足以下公式,
E1/E2≈(c1*Tlong)2/(c2*Tlong)2≈d1 2/d2 2
式中,d1、d1表示不同的目标区域在一段时间Tlong内的所述剪切波的传播距离,c1和c2表示不同的目标区域的剪切波的传播速度。
进一步的,当显示所述测量计算结果时,采用波前电影图、传播距离分布图、传播时间分布图、传播距离曲线图、传播距离示意图、弹性比值示意图、传播速度平方比示意图、传播距离平方比示意图中的至少一种。
进一步的,当显示所述测量计算结果时,采用灰阶或彩色编码模式。
一种弹性测量检测系统,包括超声探头、控制模块、信号处理模块、计算模块及显示装置,所述超声探头设有收发模块,所述超声探头的收发模块、信号处理模块、计算模块、显示装置依次连接,所述控制模块连接于所述收发模块,
所述收发模块用于发射推动脉冲、追踪脉冲及参考脉冲,并接收所述追踪脉冲与所述参考脉冲的回波数据;
所述控制模块用于控制所述收发模块发射推动脉冲、追踪脉冲及参考脉冲;
所述信号处理模块用于对回波数据进行信号预处理;
所述计算模块用于对波束合成输出的信号进行处理计算;
所述显示装置用于将所述计算模块生成的弹性图像进行显示。
本发明的方法不需要操作者控制按压探头的力度即可对反映组织弹性的相关参数进行计算获取,使得测量的重复性和稳定性得到保证。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是第一较佳实施方式提供的弹性测量检测方法流程示意图;
图2至图4是第一较佳实施方式提供的弹性测量检测方法中产生剪切波传播的推动波发射示意图;
图5是第一较佳实施方式提供的弹性测量检测方法中的位移计算示意图;
图6是第一较佳实施方式提供的弹性测量检测方法中不同目标区域内的传播距离-时间曲线示意图;
图7是第二较佳实施方式提供的弹性测量检测方法流程示意图;
图8是本发明提供的弹性测量检测系统的构成示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,本发明第一较佳实施方式提供一种弹性测量检测方法,包括以下步骤:
步骤S101,于目标区域产生剪切波。本实施例中,目标区域可根据弹性测量的需求确定,目标区域后向目标区域发射推动脉冲(push脉冲)以产生剪切波。所述目标区域可通过任意适用方式选择确定,如可采用常规二维B模式成像、常规弹性成像E模式等各类适用成像检测方式初步检测后确定,也可根据检测需求选定。
所述目标区域的数量可以为一个或者多个。当目标区域的数量为多个时,多个目标区域的各自纵向深度或横向位置可以不同。可以理解的是,当同时对多个目标区域进行测量时,可通过后续步骤获取多个目标区域之间的平均距离比值,从而反映所述多个目标区域之间的弹性差异。
确定目标区域后,即可通过探头向目标区域附近发射特定脉冲,即推动脉冲。所述推动脉冲发射时长比常规超声发射脉冲长,约为几十到几百微秒量级。所述推动脉冲可以在发射位置产生剪切波源,自剪切波源发射剪切波并沿着与所述推动脉冲发射方向不同的方向传播。
如图2至图4所示,在本实施例中,发射所述推动脉冲可以采用聚焦发射或无聚焦发射,所述推动脉冲也可采用发射一次或者连续发射多次,还可以调整连续多次发射的聚焦深度或/和横向位置,以形成特定的剪切波传播方向和拓宽剪切波的传播范围。
当确认并设置多个目标区域时,各个目标区域可以发射具有相同或不同的发射聚焦方式、连续发射次数、聚焦位置的多组推动脉冲,每组推动脉冲针对各自目标区域发射并在该目标区域附近产生剪切波,以便对相应剪切波的回波数据进行分析计算。由于推动脉冲发射需要较大的声场能量以提高剪切波强度,多组推动脉冲发射之间可能需要间隔一定时间以保证能量安全性,每组推动脉冲的发射的强度和各组推动脉冲的发射的时间间隔由系统预先确定。
多个目标区域也可以共用同一组推动脉冲发射,即一组推动脉冲发射后,产生剪切波传播,剪切波传播会经过多个目标区域,然后获取上述经过的多个目标区域的回波数据进行分析计算。系统可以根据目标区域深度和横向位置的远近来决定是否需要共用推动脉冲发射。
当向目标区域内发射多组脉冲时,每组发射可包括多次连续发射,每次发射的脉冲可具有不同的聚焦位置、纵向深度、横向位置。可通过调节各次发射的位置和/或发射间隔时间,来控制其产生的剪切波的传播方向、传播宽度。多组不同方向或者不同位置的剪切波传播后的计算结果,可以进行复合、加权等处理,以增加结果准确性。
步骤S102,向目标区域发射追踪脉冲(tracking波)并接收其回波信息,并获取回波数据。
在本步骤中,所述回波数据指追踪脉冲的追踪脉冲的回波数据。向目标区域发射推动脉冲并产生剪切波后,发射追踪脉冲并接收其回波信息,从而得到目标区域内一段时间内一段传播范围内的追踪脉冲的回波数据。追踪脉冲的发射间隔时间可预先确定。所述追踪脉冲的回波数据记录了剪切波传播过程中,上述一段传播范围内各位置上的组织信息。
步骤S103,根据所述回波数据进行数据计算,获取测量计算结果。在本实施例中,可通过追踪脉冲的回波数据进行数据计算并获取位移、质点速度、一段时间内各时刻剪切波的传播距离、一段距离内剪切波的传播时间等计算结果,并加以显示。
如图5所示,剪切波经过目标区域时,目标区域的组织的相应位置会产生纵向位移,即位置变化。所述位移计算即指对目标区域的组织位移进行计算并获取位移数据。位移计算可进一步包括以下步骤:
步骤S1031,取目标区域内任一位置的任一时刻的回波数据作为基准回波数据。
步骤S1032,将基准回波数据分段,每段设为基准回波核数据(kernel),设置所述基准回波核数据的基准回波核数据中心。可以理解的是,每段所述基准回波核数据的长度与所述基准回波核数据中心之间的间隔可由系统预先确定。
步骤S1033,获取各个时刻的回波数据上与所述基准回波核数据中心之间互相关性最大的位置,取该互相关性最大的位置与所述基准回波核数据的中心位置之差为各个时刻的位移数据。
本步骤中,所述各个时刻的回波数据上与所述基准回波核数据中心互相关性最大的位置与所述基准回波核数据的中心位置之差即为该基准回波核数据位置对应的位移。本步骤中,可取目标区域内某一横向位置上不同时刻的追踪脉冲的回波数据,分别与该基准回波数据进行互相关比较,可以得到各个深度位置上的组织在不同时刻的位移数据。该位移反映了剪切波传播导致的组织内质点相对于上述传播时刻的相对纵向位移。
在步骤S1033中,基准回波数据与不同时刻的推动脉冲的回波数据可采用块匹配(block-matching)的方式进行互相关比较,以获取位置变化的位移值。可以理解的是,块匹配的算法选择及块匹配位移的匹配标准也可自行设定。并且由于剪切波传播产生的位移非常小,也可以将回波数据上相同位置的分段数据进行相移(phase shift)计算,直接通过相移计算出位移值,该位置值即为位移数据。
对目标区域内的一个或多个位置,都可以进行上述位移计算过程。对所得位移数据,可以在空间或时间维度进行一定的滤波处理以降低噪声。
步骤S103中根据所述回波数据进行数据计算时,可包括质点速度计算。所述质点速度计包括以下步骤:
步骤S1034,获取目标区域的组织的位移数据后,取同一个位置不同时刻的位移数据,沿时间方向求梯度,即可得到该位置处不同时刻的质点速度曲线。
该质点速度曲线反映了剪切波传播过程中各个位置上质点在不同时刻的运动瞬时速度。对目标区域内的一个或多个位置,都可以进行上述质点速度计算过程。对所得质点速度数据,可以在空间或时间维度进行一定的滤波处理以降低噪声。
步骤S103中根据所述回波数据进行数据计算时,可包括传播距离计算。所述传播距离计算包括以下步骤:
步骤S1035,取目标区域中任一位置对应的各个时刻的位移数据,形成位移-时间曲线,找出该曲线上峰值对应的时刻。该峰值对应的时刻反映了剪切波经过该位置的时刻。可以理解的是,剪切波经过该位置的时刻也为质点速度-时间曲线的峰值对应的时刻,也可以利用该位置的质点速度-时间曲线来寻找对应的剪切波经过该位置的时刻。本步骤中,可以直接取出多个特定时刻附近一段时间ΔT内的位移或速度数据,进行上述传播距离计算。区别仅在于计算出来的距离难以准确判断位于目标区域内的哪个位置,但能判断位移剪切波传播多长时间之后。
步骤S1036,获取目标区域内的所述峰值对应的时刻前后一段时间段ΔT内的位移-横向位置曲线或者质点速度-横向位置曲线,对所述时间段ΔT内各个时刻的位移-横向位置曲线或者质点速度-横向位置曲线进行互相关判断,获取所述时间段内的剪切波横向位置的差异,所述差异反映了剪切波在上述一段时间ΔT内的传播距离。对目标区域内的所有位置均进行上述计算处理,可以得到剪切波经过目标区域内各位置附近时在一段时间ΔT内的传播距离。在本步骤中,所述时间段的长度ΔT可自行确定,该时间段可以为上述时刻之前至之后、上述时刻之前或上述时刻之后。
可以理解的是,在步骤S1035中,可直接取出多个特定时刻附近一段时间ΔT内的位移或速度数据,进行上述传播距离计算。对上述传播距离求均值,可以得到目标区域内的一段时间ΔT内的平均传播距离,或者剪切波传播过程中一段时间ΔT内的平均传播距离。
将多个连续时刻对应的ΔT内的传播距离取出,如图6所示,可以得到传播距离-时间曲线,也可以累积起来得到较长一段时间Tlong内的总传播距离。还可以取出两个不同位置对应的各个时刻的位移数据,进行互相关比较,找出两者的时间差异,则可以得到剪切波在这两个位置之间的传播时间。
步骤S103中根据所述回波数据进行数据计算时,可包括弹性比值计算,所述弹性比值主要指系统有多个目标区域时,可以计算出各个目标区域的杨氏模量比值,以反映各目标区域间的硬度差异程度。
在一定条件下,剪切波的传播速度与组织硬度有近似固定的关系:
E=3ρc2
上式中,ρ表示组织密度,E表示组织的杨氏模量值,c表示剪切波的传播速度。一定条件下,杨氏模量越大,意味着组织硬度越大。假设两个组织间的杨氏模量不同,分别为E1和E2,则其比值满足:
E1/E2=3ρ1c1 2/(3ρ2c2 2)
假设组织密度相近,即ρ1≈ρ2,则:
E1/E2≈c1 2/c2 2
可以使用一段时间Tlong内的剪切波传播距离的平方比值来计算,即:
E1/E2≈(c1*Tlong)2/(c2*Tlong)2≈d1 2/d2 2
式中,d1、d1代表不同的目标区域在一段时间Tlong内的剪切波传播距离,c1和c2表示不同的目标区域的剪切波的传播速度。
当然,也可以直接计算出各个目标区域内的传播速度,然后计算速度平方之比来得到弹性比值。
步骤S103中根据所述回波数据进行数据计算时,可包括传播速度计算。
对各个目标区域中的任意深度,取出其各横向位置的质点速度-时间曲线或位移-时间曲线,利用互相关比较找到各曲线之间的时移,该时移代表剪切波经过上述参与互相关比较的两个横向位置之间的时间差。横向位置差与时移的比值即代表相应位置处的传播速度。还比如可以利用波的传播方程直接推导出近似计算公式如下:
式中,c表示传播速度,uz可以视为纵向位移数据或纵向速度数据,x代表横向坐标,z代表纵向坐标。
步骤S104,显示测量计算结果。经过上述弹性测量计算后,有多种方式将结果显示。
如采用波前电影图,即将某时刻的位移数据或质点速度数据取出,形成当前时刻在目标区域各位置的分布图,该分布图可以反映当前时刻下剪切波的传播位置。当多个时刻的位移数据或质点速度数据分布图连续播放时,即可形成剪切波波前传播电影图,直观表现剪切波的传播过程。
如采用传播距离分布图,即将目标区域内各个位置处的ΔT时间内的传播距离进行成像显示。
如采用传播时间分布图,即将目标区域内各位置附近某固定距离之间的传播时间进行成像显示。当然也可以曲线形式显示传播时间。
如采用传播距离曲线图,即将目标区域内从某个时刻开始往后一段时间内的传播距离-时间曲线绘制出来,多个目标区域时可显示多条曲线。
如直接显示目标区域内的各位置处平均传播距离示意图,或者直接显示剪切波在各个不同时刻时某小段时间内的平均传播距离示意图。
如当有多个目标区域时,显示多个目标区域内的弹性比值示意图、或传播速度平方比示意图、或传播距离平方比示意图。
上述显示过程中,可以使用灰阶或彩色编码,或与其他模式的图像叠加、融合等方式加强显示效果。
如图7所示,本发明的第二较佳实施例提供一种弹性测量检测方法,本实施例的弹性测量检测方法与第一较佳实施例大致相同,包括以下步骤:
步骤S201,向目标区域发射参考脉冲(reference波),并接收所述参考脉冲于所述目标区域内的回波信息,获取参考脉冲的回波数据,取参考脉冲的回波数据为基准回波数据。在本实施例中,在每组推动脉冲发射前发射参考脉冲并接受其在目标区域内的回波信息,获取参考脉冲的回波数据并设其为基准回波数据。所述参考脉冲的回波信息记录了在产生剪切波并开始传播前的目标区域的组织信息。
步骤S202,于目标区域产生剪切波。
步骤S203,向目标区域发射追踪脉冲并接收其回波信息,并获取回波数据。在本实施例的本步骤中,所述回波数据指追踪脉冲的追踪脉冲的回波数据。
步骤S204,根据所述回波数据进行数据计算,获取测量计算结果。在本实施例中,将参考脉冲的回波数据取为基准回波数据,通过追踪脉冲的回波数据与作为基准回波数据的参考脉冲的回波数据进行数据计算并获取位移、质点速度、一段时间内各时刻剪切波的传播距离、一段距离内剪切波的传播时间等计算结果。
取目标区域内某一横向位置上不同时刻的追踪脉冲的回波数据,分别与该横向位置上的基准回波数据进行互相关比较,可以得到各个深度位置上的组织在不同时刻的位移数据。该位移主要反映了剪切波传播导致的组织内质点相对于剪切波传播前的组织内质点的纵向位移。质点速度、一段时间内各时刻剪切波的传播距离、一段距离内剪切波的传播时间等计算结果等计算过程于第一较佳实施大致相同,在此不再赘述。
步骤S205,显示测量计算结果。
本实施例中,在于目标区域产生剪切波之前先行向目标区域发射参考脉冲,并获取该参考脉冲的回波数据作为基准回波数据,从而获取剪切波传播导致的组织内质点相对于剪切波传播前的组织内质点的纵向位移、相应质点速度、一段时间内各时刻剪切波的传播距离、一段距离内剪切波的传播时间等计算结果等计算结果。本实施例中的其他各个步骤与第一较佳实施例大致相同,在此不再赘述。
如图8所示,本发明还提供一种超声弹性测量系统,包括超声探头11、控制模块12、信号处理模块13、计算模块15及显示装置17,所述超声探头11设有收发模块110,所述超声探头11的收发模块110、信号处理模块13、计算模块15、显示装置17依次连接,所述控制模块12连接于所述收发模块110。其中:
所述收发模块110用于发射推动脉冲、追踪脉冲及参考脉冲,并接收所述追踪脉冲与所述参考脉冲的回波数据。
所述控制模块12用于控制所述收发模块110发射推动脉冲、追踪脉冲及参考脉冲。实际使用中,所述控制模块12按照预设时间间隔发射由推动脉冲、参考脉冲、追踪脉冲构成的特定的超声序列,从而便于在目标区域产生剪切波并供所述超声探头11的收发模块110接收相应的回波数据。
所述信号处理模块13用于对回波数据进行信号预处理,从而便于计算模块15进行后续计算,所述信号预处理可包括波束合成处理,还可包括如信号放大、模数转换、正交分解等。
计算模块15用于对波束合成输出的信号进行处理计算。所述计算模块15可用于计算剪切波的各类传播参数,具体的,所述计算模块15包括:
质点速度计算单元151,用于计算目标组织的质点速度;
传播距离计算单元152,用于计算剪切波的传播距离;
杨氏模量比值计算单元153,用于计算目标组织的杨氏模量比值。
计算模块15通过计算某个位置开始各个时刻的剪切波传播距离,或某段时间内的平均距离,或不同位置开始某段时间内的平均距离比值,或几个目标区域内的平均传播速度比值等,从而反映目标区域的物理特性,并生成相应的弹性图像。
显示装置17用于将所述计算模块15生成的弹性图像进行显示。
本发明的弹性测量检测方法及系统,通过向待检测组织内部的目标区域内发射一系列特定脉冲产生剪切波的传播,然后在剪切波传播路径范围内不断发射接收回波信号,从该回波信号中提取出一段时间内剪切波的传播距离曲线,且计算出平均距离,最终显示。如果相同发射条件下,在同一时间内,剪切波传播的越远,则说明组织的硬度越高,否则说明组织越软。
本发明的方法不需要操作者控制按压探头的力度即可对反映组织弹性的相关参数进行计算获取,使得测量的重复性和稳定性得到保证。另一方面,本发明可以得到定量的参数值,特别适用于对均匀组织的弹性测量。同时,本发明还可提供不同组织间的弹性差异参数,用以反映不同组织之间的弹性差异程度。
以上所揭露的仅为本发明一种较佳实施例而已,当然不能以此来限定本发明之权利范围,本领域普通技术人员可以理解实现上述实施例的全部或部分流程,并依本发明权利要求所作的等同变化,仍属于发明所涵盖的范围。
Claims (17)
1.一种弹性测量检测方法,其特征在于:包括以下步骤,
于目标区域产生剪切波;
向所述目标区域发射追踪脉冲并接收所述追踪脉冲的回波数据;
根据所述回波数据进行数据计算以获取目标区域内的剪切波传播距离;以及
显示所述剪切波传播距离。
2.根据权利要求1所述的弹性测量检测方法,其特征在于,根据所述回波数据进行数据计算以获取目标区域内的剪切波传播距离,包括:获取所述目标区域内各位置处一段时间内的剪切波传播距离;
显示所述剪切波传播距离包括:
分别显示所述目标区域内各位置处一段时间内的剪切波传播距离,形成剪切波传播距离分布图;或者,根据所述回波数据获得目标区域的超声图像,在所述超声图像上叠加显示所述目标区域内各位置处一段时间内的剪切波传播距离,形成剪切波传播距离分布图;或者,根据所述目标区域内各位置处一段时间内的剪切波传播距离,显示所述目标区域内的各位置处的平均传播距离。
3.根据权利要求2所述的弹性测量检测方法,其特征在于,获取所述目标区域内各位置处一段时间内的剪切波传播距离,包括:
确定剪切波经过所述目标区域中各位置的任一位置的时刻;
计算剪切波经过该位置的时刻附近的一段时间内的剪切波传播距离。
4.根据权利要求3所述的弹性测量检测方法,其特征在于,确定剪切波经过所述目标区域中各位置的任一位置的时刻,包括:
根据所述回波数据进行数据计算获得目标区域的位移数据或质点速度数据;
取目标区域中任一位置对应的各个时刻的位移数据或质点速度数据,进而得到剪切波经过该位置的时刻。
5.根据权利要求2所述的弹性测量检测方法,其特征在于:采用灰阶或彩色编码模式,显示各位置处一段时间内的所述剪切波传播距离和平均传播距离。
6.根据权利要求1所述的弹性测量检测方法,其特征在于,根据所述回波数据进行数据计算以获取目标区域内的剪切波传播距离,包括:
取多个时刻;
分别获取剪切波于所述目标区域中在各个时刻附近的一段时间内的剪切波传播距离。
7.根据权利要求6所述的弹性测量检测方法,其特征在于,显示所述剪切波传播距离包括:以曲线形式显示各个时刻附近的一段时间内的剪切波传播距离。
8.根据权利要求6所述的弹性测量检测方法,其特征在于,显示所述剪切波传播距离包括:根据剪切波于所述目标区域内在各个时刻附近的一段时间内的剪切波传播距离,显示剪切波传播过程中一段时间内的平均传播距离。
9.根据权利要求6所述的弹性测量检测方法,其特征在于,还包括:累积多个时刻附近的一段时间的剪切波传播距离,得到各个时刻目标区域内的剪切波的总传播距离;
显示所述剪切波传播距离包括:以曲线形式显示各个时刻目标区域内的剪切波的总传播距离。
10.根据权利要求2或6所述的弹性测量检测方法,其特征在于,所述经过任一位置的时刻附近的一段时间包括经过任一位置的时刻的向前一段时间、向后一段时间或前后一段时间;所述各个时刻附近的一段时间包括所述各个时刻向前的一段时间、各个时刻向后的一段时间、或各个时刻前后的一段时间。
11.根据权利要求1所述的弹性测量检测方法,其特征在于:
所述弹性测量检测方法于目标区域产生剪切波前还包括以下步骤,
向目标区域发射参考脉冲,接收所述参考脉冲于所述目标区域内的回波信息,获取所述参考脉冲的回波数据。
12.根据权利要求11所述的弹性测量检测方法,其特征在于:所述根据所述回波数据进行数据计算获得目标区域的位移数据包括以下步骤:
取目标区域内任一位置的任一时刻的回波数据作为基准回波数据;
将基准回波数据分为多段基准回波核数据,设置所述基准回波核数据的基准回波核数据中心;
获取各个时刻的回波数据上与所述基准回波核数据中心之间互相关性最大的位置,取该互相关性最大的位置与所述基准回波核数据的中心位置之差为各个时刻的位移数据。
13.根据权利要求12所述的弹性测量检测方法,其特征在于:当取目标区域内任一位置的任一时刻的回波数据作为基准回波数据时,所述基准回波数据采用任一时刻的追踪脉冲的回波数据或任一时刻的参考脉冲的回波数据。
14.根据权利要求4所述的弹性测量检测方法,其特征在于:还能采用以下方法获得剪切波经过所述位置的时刻,包括以下步骤:
取目标区域中任一位置对应的各个时刻的位移数据,沿时间方向求梯度以获取该位置处各个时刻的质点速度数据,进而得到剪切波经过该位置的时刻。
15.根据权利要求1-4中任一项所述的弹性测量检测方法,其特征在于:获取目标区域内的剪切波传播距离包括以下步骤,
取目标区域中任一位置对应的各个时刻的位移数据,形成位移-时间曲线,找出该曲线上峰值对应的时刻,即为剪切波经过该位置的时刻;
获取目标区域内的所述峰值对应的时刻前后一段时间段内的位移-横向位置曲线或者质点速度-横向位置曲线,对所述时间段内各个时刻的位移-横向位置曲线或者质点速度-横向位置曲线进行互相关判断,获取所述时间段内的剪切波的传播距离。
16.根据权利要求15所述的弹性测量检测方法,其特征在于:当根据所述回波数据进行数据计算以获取测量计算结果时,获取不同目标区域之间的杨氏模量比值,当获取所述杨氏模量比值时,满足以下公式,
E1/E2≈(c1*Tlong)2/(c2*Tlong)2≈d1 2/d2 2
式中,d1、d1表示不同的目标区域在一段时间Tlong内的所述剪切波的传播距离,c1和c2表示不同的目标区域的剪切波的传播速度。
17.一种弹性测量检测系统,其特征在于:包括超声探头、控制模块、信号处理模块、计算模块及显示装置,所述超声探头设有收发模块,所述超声探头的收发模块、信号处理模块、计算模块、显示装置依次连接,所述控制模块连接于所述收发模块,
所述收发模块用于发射推动脉冲、追踪脉冲及参考脉冲,并接收所述追踪脉冲与所述参考脉冲的回波数据;
所述控制模块用于控制所述收发模块发射推动脉冲、追踪脉冲及参考脉冲;
所述信号处理模块用于对回波数据进行信号预处理;
所述计算模块用于对波束合成输出的信号进行处理计算,获得目标区域内的剪切波传播距离;
所述显示装置用于将所述计算模块生成的弹性图像进行显示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910741685.8A CN110368031B (zh) | 2014-09-03 | 2014-09-03 | 弹性测量检测方法及系统 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910741685.8A CN110368031B (zh) | 2014-09-03 | 2014-09-03 | 弹性测量检测方法及系统 |
PCT/CN2014/085821 WO2016033752A1 (zh) | 2014-09-03 | 2014-09-03 | 弹性测量检测方法及系统 |
CN201480081082.5A CN106572838B (zh) | 2014-09-03 | 2014-09-03 | 弹性测量检测方法及系统 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201480081082.5A Division CN106572838B (zh) | 2014-09-03 | 2014-09-03 | 弹性测量检测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110368031A true CN110368031A (zh) | 2019-10-25 |
CN110368031B CN110368031B (zh) | 2022-10-04 |
Family
ID=55439001
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910741683.9A Active CN110432926B (zh) | 2014-09-03 | 2014-09-03 | 弹性测量检测方法及系统 |
CN201910741685.8A Active CN110368031B (zh) | 2014-09-03 | 2014-09-03 | 弹性测量检测方法及系统 |
CN201480081082.5A Active CN106572838B (zh) | 2014-09-03 | 2014-09-03 | 弹性测量检测方法及系统 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910741683.9A Active CN110432926B (zh) | 2014-09-03 | 2014-09-03 | 弹性测量检测方法及系统 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201480081082.5A Active CN106572838B (zh) | 2014-09-03 | 2014-09-03 | 弹性测量检测方法及系统 |
Country Status (2)
Country | Link |
---|---|
CN (3) | CN110432926B (zh) |
WO (1) | WO2016033752A1 (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6996035B2 (ja) * | 2017-11-02 | 2022-01-17 | 富士フイルムヘルスケア株式会社 | 超音波診断装置、および、生体組織の物性評価方法 |
CN108158610B (zh) * | 2018-01-16 | 2024-02-06 | 苏州国科昂卓医疗科技有限公司 | 一种弹性成像方法、装置、设备及超声成像探头 |
WO2019196033A1 (zh) * | 2018-04-11 | 2019-10-17 | 深圳迈瑞生物医疗电子股份有限公司 | 超声弹性成像方法和系统 |
CN111401108B (zh) * | 2019-01-02 | 2021-04-06 | 无锡海斯凯尔医学技术有限公司 | 斜率参数提取方法、装置及计算机可读存储介质 |
CN111399038B (zh) * | 2019-01-02 | 2023-04-07 | 无锡海斯凯尔医学技术有限公司 | 斜率参数提取方法、装置及计算机可读存储介质 |
CN111399039B (zh) * | 2019-01-02 | 2023-04-07 | 无锡海斯凯尔医学技术有限公司 | 斜率参数提取方法、装置及存储介质 |
CN109893172B (zh) * | 2019-02-22 | 2020-06-19 | 清华大学 | 基于弹性成像的力学参数的确定方法及装置、计算机设备 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101065067A (zh) * | 2004-08-05 | 2007-10-31 | 株式会社日立医药 | 弹性图像显示方法以及超声波诊断装置 |
US20080249408A1 (en) * | 2007-02-09 | 2008-10-09 | Palmeri Mark L | Methods, Systems and Computer Program Products for Ultrasound Shear Wave Velocity Estimation and Shear Modulus Reconstruction |
CN101869485A (zh) * | 2010-06-23 | 2010-10-27 | 深圳大学 | 超声成像方法及装置 |
US20110319756A1 (en) * | 2009-03-09 | 2011-12-29 | Mayo Foundation For Medical Education And Research | Method for Ultrasound Vibrometry Using Orthogonal Basis Functions |
CN102469989A (zh) * | 2009-07-07 | 2012-05-23 | 株式会社日立医疗器械 | 超声波诊断装置和超声波测量方法 |
CN102667522A (zh) * | 2009-11-25 | 2012-09-12 | 皇家飞利浦电子股份有限公司 | 采用聚焦扫描线波束形成的超声剪切波成像 |
CN103269639A (zh) * | 2010-12-22 | 2013-08-28 | 皇家飞利浦电子股份有限公司 | 利用质心估计剪切波速度 |
CN103458800A (zh) * | 2011-03-30 | 2013-12-18 | 日立阿洛卡医疗株式会社 | 用于超声波成像的方法和装置 |
CN103720489A (zh) * | 2013-12-30 | 2014-04-16 | 中国科学院深圳先进技术研究院 | 病变组织生长监测方法和系统 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7374538B2 (en) * | 2000-04-05 | 2008-05-20 | Duke University | Methods, systems, and computer program products for ultrasound measurements using receive mode parallel processing |
US6702744B2 (en) * | 2001-06-20 | 2004-03-09 | Advanced Cardiovascular Systems, Inc. | Agents that stimulate therapeutic angiogenesis and techniques and devices that enable their delivery |
CN102481137B (zh) * | 2009-06-30 | 2015-03-18 | 皇家飞利浦电子股份有限公司 | 用于剪切波分散振动测定的超声设备和方法 |
RU2580419C2 (ru) * | 2009-07-17 | 2016-04-10 | Конинклейке Филипс Электроникс, Н.В. | Дискретизация ультразвуковой виброметрии рассеяния поперечных волн с высоким пространственным разрешением |
US8500639B2 (en) * | 2009-09-11 | 2013-08-06 | Mr Holdings (Hk) Limited | Systems and methods for shear wave field formation |
CN101699280B (zh) * | 2009-10-15 | 2011-08-17 | 北京索瑞特医学技术有限公司 | 超声无损检测粘弹性介质弹性的方法及其装置 |
CN103250072B (zh) * | 2010-04-27 | 2017-02-22 | 哈里伯顿能源服务公司 | 通过使用钻头地震对断裂的干涉测量钻头成像、时间反演成像来进行断裂表征、和经由时间反演的声学和电震对断裂产生的监测 |
EP2638388A1 (en) * | 2010-11-14 | 2013-09-18 | University Of Utah Research Foundation | Ultrasonic method and system for determining tissue pathology |
BR112013021791B1 (pt) * | 2011-02-25 | 2020-11-17 | Mayo Foundation For Medical Education Ano Research | método para medir uma propriedade mecânica de um paciente com um sistema de ultrassom |
CN102824194B (zh) * | 2011-06-14 | 2016-08-03 | 深圳迈瑞生物医疗电子股份有限公司 | 一种弹性成像中的位移检测方法及装置 |
CN102525568A (zh) * | 2012-01-17 | 2012-07-04 | 北京索瑞特医学技术有限公司 | 一种弹性减影成像方法 |
US10357226B2 (en) * | 2012-03-12 | 2019-07-23 | Mayo Foundation For Medical Education And Research | System and method for model-independent quantification of tissue viscoelastic properties using ultrasound |
CN103300890B (zh) * | 2012-03-16 | 2016-06-08 | 通用电气公司 | 用于测量组织机械特性的系统及方法 |
US9220479B2 (en) * | 2012-03-30 | 2015-12-29 | Hitachi Aloka Medical, Ltd. | Methods and apparatus for ultrasound imaging |
JP5733835B2 (ja) * | 2012-08-27 | 2015-06-10 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 計測装置及び超音波診断装置 |
CN103845074B (zh) * | 2012-11-28 | 2017-12-12 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性成像系统和方法 |
CN103908289B (zh) * | 2012-12-31 | 2019-11-12 | Ge医疗系统环球技术有限公司 | 消除剪切波中的背景噪声的方法和相应的超声成像系统 |
CN103462643B (zh) * | 2013-09-29 | 2015-02-11 | 深圳市开立科技有限公司 | 一种剪切波速度测量方法、装置和系统 |
-
2014
- 2014-09-03 CN CN201910741683.9A patent/CN110432926B/zh active Active
- 2014-09-03 CN CN201910741685.8A patent/CN110368031B/zh active Active
- 2014-09-03 CN CN201480081082.5A patent/CN106572838B/zh active Active
- 2014-09-03 WO PCT/CN2014/085821 patent/WO2016033752A1/zh active Application Filing
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101065067A (zh) * | 2004-08-05 | 2007-10-31 | 株式会社日立医药 | 弹性图像显示方法以及超声波诊断装置 |
US20080249408A1 (en) * | 2007-02-09 | 2008-10-09 | Palmeri Mark L | Methods, Systems and Computer Program Products for Ultrasound Shear Wave Velocity Estimation and Shear Modulus Reconstruction |
US20110319756A1 (en) * | 2009-03-09 | 2011-12-29 | Mayo Foundation For Medical Education And Research | Method for Ultrasound Vibrometry Using Orthogonal Basis Functions |
CN102469989A (zh) * | 2009-07-07 | 2012-05-23 | 株式会社日立医疗器械 | 超声波诊断装置和超声波测量方法 |
CN102667522A (zh) * | 2009-11-25 | 2012-09-12 | 皇家飞利浦电子股份有限公司 | 采用聚焦扫描线波束形成的超声剪切波成像 |
CN101869485A (zh) * | 2010-06-23 | 2010-10-27 | 深圳大学 | 超声成像方法及装置 |
CN103269639A (zh) * | 2010-12-22 | 2013-08-28 | 皇家飞利浦电子股份有限公司 | 利用质心估计剪切波速度 |
CN103458800A (zh) * | 2011-03-30 | 2013-12-18 | 日立阿洛卡医疗株式会社 | 用于超声波成像的方法和装置 |
CN103720489A (zh) * | 2013-12-30 | 2014-04-16 | 中国科学院深圳先进技术研究院 | 病变组织生长监测方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110432926B (zh) | 2022-06-07 |
CN110368031B (zh) | 2022-10-04 |
CN106572838B (zh) | 2019-09-06 |
CN106572838A (zh) | 2017-04-19 |
CN110432926A (zh) | 2019-11-12 |
WO2016033752A1 (zh) | 2016-03-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106572838B (zh) | 弹性测量检测方法及系统 | |
CN104203112B (zh) | 超声波诊断装置以及弹性评价方法 | |
CN103988092B (zh) | 位置测量装置 | |
CN105816205B (zh) | 声辐射力脉冲成像中的稀疏跟踪 | |
CN103462643B (zh) | 一种剪切波速度测量方法、装置和系统 | |
EP3153884B1 (en) | Detection apparatus, fish finder, and radar | |
CN106419961A (zh) | 声学辐射力成像中的自适应运动估计 | |
CN107072639B (zh) | 使用切变波进行测量的超声系统和方法 | |
US20140276003A1 (en) | Systems and Methods to Detect and Present Interventional Devices via Ultrasound Imaging | |
CN103239258A (zh) | 采用超声波的同轴切变波表征 | |
FR3072870A1 (fr) | Estimation viscoelastique d'un tissu a partir d'une vitesse de cisaillement dans une imagerie medicale a ultrasons | |
CN106659471A (zh) | 超声剪切波弹性摄影中的自相关引导的互相关 | |
CN105877783B (zh) | 二维剪切波弹性成像方法和装置 | |
CN101828929B (zh) | 利用视在位移的多普勒血流速度矢量测量方法 | |
US6804167B2 (en) | Bi-directional temporal correlation SONAR | |
CN106955125A (zh) | 声学辐射力脉冲成像中的运动无关性 | |
KR101055500B1 (ko) | Bc-모드 영상을 형성하는 초음파 시스템 및 방법 | |
CN109259801A (zh) | 一种剪切波弹性成像方法及装置 | |
EP2901939A1 (en) | Ultrasonic measurement apparatus and ultrasonic measurement method | |
CN107296629A (zh) | 超声波诊断装置和超声波诊断方法 | |
US11452503B2 (en) | Shear wave imaging based on ultrasound with increased pulse repetition frequency | |
RU2338220C1 (ru) | Способ измерения скорости снаряда | |
CN111050661B (zh) | 一种剪切波弹性测量方法及剪切波弹性成像系统 | |
RU167164U1 (ru) | Гидроакустическая многолучевая система для определения численности и параметров рыб | |
RU2798390C1 (ru) | Способ пассивного определения координат шумящего объекта |
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 |