CN103845074B - 一种超声弹性成像系统和方法 - Google Patents
一种超声弹性成像系统和方法 Download PDFInfo
- Publication number
- CN103845074B CN103845074B CN201210495070.XA CN201210495070A CN103845074B CN 103845074 B CN103845074 B CN 103845074B CN 201210495070 A CN201210495070 A CN 201210495070A CN 103845074 B CN103845074 B CN 103845074B
- Authority
- CN
- China
- Prior art keywords
- elastic image
- compression parameters
- compression
- parameter
- quality
- 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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
- A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
-
- 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
- Health & Medical Sciences (AREA)
- Life Sciences & Earth 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)
- Vascular Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本申请公开了一种超声弹性成像系统和方法,其中系统包括:弹性检测模块,用于对波束合成输出的信号进行处理,计算反映待检测目标的弹性的物理量,根据该物理量生成相应的弹性图像;参数计算模块,用于计算与所述弹性图像对应的压缩参数,所述压缩参数用于反映探头对待检测目标的压缩程度;显示模块,用于显示弹性图像及其对应的压缩参数。本申请通过在显示弹性图像的同时显示其对应的压缩参数,用户可以直观地根据显示结果来判断探头对待检测目标的压缩是否需要调整,从而给用户尤其是操作经验不足的用户提供一定的指导,使得操作更简便。
Description
技术领域
本申请涉及超声成像技术,尤其涉及一种超声弹性成像系统和方法。
背景技术
超声弹性成像是一种常用的超声成像技术,其基本原理为:将探头轻微压缩目标组织或者借助人体自身的呼吸、血管搏动等过程对组织形成一定的压力,获取压缩前、后两帧超声回波信号,组织被压缩时,组织内将产生一个沿压缩方向的应变,如果组织内部杨氏模量分布不均匀,组织内的应变分布也将有所差异;然后通过一些方法检测出组织的应变信息,输出至界面并以弹性图像形式直观显示出来。输出所得弹性图像的质量好坏受压缩过程的影响。
在采集弹性图像的过程中,为获得更好的图像,用户需要调整探头压缩组织的力度、频率等,其需要依靠大量的探头使用经验或者培训;由于目前超声弹性成像系统的限制,用户不能简便地判断探头对组织的压缩力度、频率等情况,从而用户对超声弹性成像系统的使用存在一定的困难。
发明内容
根据本申请的第一方面,提供一种超声弹性成像系统,包括:弹性检测模块,用于对波束合成输出的信号进行处理,计算反映待检测目标的弹性的物理量,根据该物理量生成相应的弹性图像;参数计算模块,用于计算与所述弹性图像对应的压缩参数,所述压缩参数用于反映探头对待检测目标的压缩程度;显示模块,用于显示弹性图像及其对应的压缩参数。
根据本申请的第二方面,提供一种超声弹性成像方法,包括:弹性检测步骤,对波束合成输出的信号进行处理,计算反映待检测目标的弹性的物理量,根据该物理量生成相应的弹性图像;参数计算步骤,计算与所述弹性图像对应的压缩参数,所述压缩参数用于反映探头对待检测目标的压缩程度;显示步骤,显示弹性图像及其对应的压缩参数。
本申请的有益效果是:通过在显示弹性图像的同时显示其对应的压缩参数,用户可以直观地根据显示结果来判断探头对待检测目标的压缩是否需要调整,从而给用户尤其是操作经验不足的用户提供一定的指导,使得操作更简便。
附图说明
图1为本申请一种实施例的超声弹性成像系统的原理性结构图;
图2为本申请一种实施例的超声弹性成像系统中对压缩参数的一种显示的示意图;
图3为本申请一种实施例的超声弹性成像系统中对压缩参数的另一种显示的示意图;
图4为本申请一种实施例的超声弹性成像系统中对压缩参数的又一种显示的示意图;
图5为本申请另一种实施例的超声弹性成像系统的原理性结构图;
图6为本申请一种实施例的超声弹性成像系统中对压缩参数和质量评估参数的一种显示的示意图;
图7为本申请一种实施例的超声弹性成像系统中对压缩参数和质量评估参数的另一种显示的示意图;
图8为本申请一种实施例的超声弹性成像系统中对压缩参数和质量评估参数的又一种显示的示意图。
具体实施方式
下面通过具体实施方式结合附图对本申请作进一步详细说明。
实施例1
本实施例的超声弹性成像系统10的原理性结构如图1所示,包括:超声探头及其发射接收模块、信号预处理模块101、B信号处理模块102、弹性处理模块103和显示模块104。探头以系统预先设定好的扫描规则进行超声发射并接收超声回波信号;接受的回波信号通过信号预处理模块101进行信号预处理,信号预处理包括波束合成处理,还可以有如信号放大、模数转换、正交分解等;波束合成后的射频信号被送往多个并行的模块进行处理,包括B信号处理模块102和弹性处理模块103,其它实施例还可以有其它并行处理模块如血流信号处理等;经B信号处理装置102和弹性处理模块103等并行处理后的图像信号被送往显示模块104进行显示,显示模块104可以根据用户的选择显示相应的内容,例如仅显示经B信号处理模块处理后的人体组织的灰阶图像,或是仅显示经弹性处理模块后得到的反映弹性信息的弹性图像,或是同时显示灰阶图像和弹性图像等等。本实施例中,探头的发射接收模块、信号预处理模块和B信号处理模块以及显示模块对灰阶图像和弹性图像的显示可以采用常用的相关超声技术实现,此外也可以增加本领域技术人员所熟知的其它处理模块,在此不做详述。当然本实施例的超声弹性成像系统也可以不包括B信号处理模块等其他与弹性处理模块并行的处理。
对于本实施例的弹性处理模块103,其包括弹性检测模块111和参数计算模块112,显示模块101可以在弹性图像121的同时显示该弹性图像对应的压缩参数122。
弹性检测模块111用于对波束合成输出的信号进行处理,计算反映待检测目标的弹性的物理量,根据该物理量生成相应的弹性图像,也就是说,通过弹性检测模块111提取出可反映待检测目标的弹性信息。可采用多种常用的弹性信息提取方法予以实现;例如,一种常用的弹性信息提取方法基于射频信号互相关性,使用绝对差值和(SAD,Sum ofAbsolute Difference)算法快速检测出相邻两帧射频信号间的位移(displacement)信息,进而对位移场求纵向(即沿超声波传播方向)梯度,获得应变(strain)信息;其他例子中可采用别的算法检测位移信息,如平方误差和(SSD,Sum of Squared Difference)等。对弹性信息检测模块得到的弹性信息即应变信息进行输出成像为弹性图像,送到显示模块104显示,从而弹性特征不同的组织之间即可直观分辨。
参数计算模块112用于计算与弹性图像对应的压缩参数,压缩参数反映探头对待检测目标的压缩程度。参数计算模块112可以是在弹性检测模块提取弹性信息的同时进行参数计算,也可以是在弹性检测模块处理结束后进行,本实施例中参数计算模块和弹性检测模块是并行处理。压缩参数包括累积应变量,累积应变量为从第一帧弹性图像的应变数据至当前帧弹性图像的应变数据的积分,具体而言,可分别基于位移信息和应变信息来计算。
(1)基于应变信息计算压缩参数
对于用户选定的ROI(感兴趣区域),获得ROI内各采样位置当前帧(假设其为第k帧)的应变数据,求平均值得到平均应变量Str_mean(k),由于探头压缩组织的上下过程,该平均应变量可能为正数,也可能为负数。
从计算的第一帧的应变数据开始,求积分,得到总的累积应变量Str_integral,如下公式所示:
该累积应变量Str_integral即可反映组织的压缩程度或形变程度,该值的绝对值越大,说明组织相对于原始位置形变越大(压缩或扩张),该值的正负符号表示组织性变相对原始位置的方向。
以上计算虽然是以用户选定的ROI内的Str_integral作为反映压缩程度的参数,但也可以使用整个扫描平面区域内的累积应变量Str_integral。
(2)基于位移信息计算压缩参数
对于用户选定的感兴趣区域ROI,获得ROI内各采样位置当前帧(假设其为第k帧)的位移数据,选出ROI中某一深度Depth的所有位移数据,比如可选用最深的(即离探头最远的)采样位置处的位移数据,求平均值得到平均位移量Dis_mean(k),由于探头压缩组织的上下过程,该平均位移量可能为正数,也可能为负数。接着,结合该深度Depth的大小,算出参考应变量Str_ref(k),如下公式所示:
Str_ref(k)=Dis_mean/Depth
其中Str_ref(k)为当前帧(即第k帧)的计算出的参考应变量。
公式中Depth不一定是物理单位,可以只是与实际物理深度对应的某个数,比如可以为与物理深度大小对应的射频采样点数。
从计算的第一帧的应变数据开始,求积分,得到总的累积应变量Str_integral,如下公式所示:
该累积应变量Str_integral即可反映组织的压缩程度或形变程度,该值的绝对值越大,说明组织相对于原始位置形变越大(压缩或扩张),该值的正负符号表示组织性变相对原始位置的方向。
同样地,以上计算虽然是以用户选定的ROI内的Str_integral作为反映压缩程度的参数,但也可以使用整个扫描平面区域内的累积应变量Str_integral。
上述压缩参数计算完毕后,将压缩参数连同与之相应的弹性信息(弹性图像)送入显示模块。显示模块在同一画面中显示弹性图像及其相应的压缩参数,其中,可采用常用的相关技术实现弹性图像的显示,而对于压缩参数的显示,可以通过多种方式显示。例如,如图2所示,在同一画面200中除了显示弹性图像201外,压缩参数可显示为与该压缩参数的数值对应的长度的条形图202(见图2中的A图)或线条203(见图2中的B图),图2中可通过增加坐标轴或是在坐标轴上设定合适的压缩参数范围以便用户能直观感知压缩程度;又例如,如图3所示,将压缩参数显示为用针的旋转角表示压缩参数的数值的模拟仪表;又例如,如图4所示,在同一画面400中显示弹性图像401,同时通过二维坐标图进行显示压缩参数,二维坐标图的横坐标为时间t,纵坐标为与时间对应的压缩参数的数值,压缩参数显示为与该压缩参数的数值对应的长度的条形图402(见图4中的A图)或线条403(见图4中的B图)。
本实施例中,显示模块中弹性图像及其压缩参数的显示随着图像的实时采集同步更新,这样,用户可以实时看到当前最新一帧弹性图像对应的压缩情况,从而用户可以判断是否需要调节手法重新采集图像,尽量将压缩程度控制在最优范围内;最终,给用户尤其是操作经验不足的用户提供一定的指导,使得操作更简便。
本实施例的超声弹性成像系统还可以作进一步改进,可还包括告警模块,用于判断连续预设时间内各帧弹性图像的压缩参数是否满足预设合适压缩强度,如果不满足则通过图像和/或声音和/或文字提示用户应当调整探头对待检测目标的压缩。
本申请中超声弹性成像方法的一种实施例与上述超声弹性成像系统的实施例相应,包括:
发射接收步骤11、在弹性成像模式下,探头以系统预先设定好的扫描规则进行超声发射并接收超声回波信号;
信号预处理步骤12、对接收到的超声回波信号进行信号预处理,预处理包括波束合成等;
弹性检测步骤13、对波束合成输出的信号进行处理,计算反映待检测目标的弹性的物理量,根据该物理量生成相应的弹性图像;
参数计算步骤14、计算与弹性图像对应的压缩参数,压缩参数用于反映探头对待检测目标的压缩程度;
显示步骤14、显示弹性图像及其对应的压缩参数。
以上各步骤的具体实现可参考上述超声弹性成像系统的实施例中各模块对应的实现过程,在此不再重述。此外上述方法实施例还可以包括对B信号进行处理以形成待检测目标的灰阶图像的步骤,具体可以参考常用B信号处理技术,在此不做详述。
实施例2
本实施例的超声弹性成像系统50的原理性结构如图5所示,包括:超声探头及其发射接收模块、信号预处理模块501、B信号处理模块502、弹性处理模块503和显示模块504,其中弹性处理模块503包括弹性检测模块511和参数计算模块512。其中,探头及其发射接收模块、信号预处理模块501和B信号处理模块502、弹性检测模块511分别与实施例1的探头及其发射接收模块、信号预处理模块101、B信号处理模块102、弹性检测模块111类似,不再赘述。
本实施例中参数计算模块512不仅用于计算与弹性图像对应的、用于反映探头对待检测目标的压缩程度的压缩参数,还用于计算反映弹性图像质量的质量评估参数,其中压缩参数的计算与实施例1中压缩参数的计算类似,不再赘述。对于质量参数,本实施例的质量参数包括形变程度参数和/或互相关检测质量参数,形变程度参数为弹性图像中感兴趣区域内各采样位置的应变量的平均值,互相关检测质量参数为根据采用的位移检测算法所对应的评分标准得到的弹性图像的评分。
(1)形变程度参数
对于弹性检测模块而言,组织的形变量太小,则位移信息太小,影响图像信噪比;若组织的形变量太大,则可能引起压缩前后信号相关性减弱,导致弹性信息的检测不准确度增加。此外,探头采集图像信号时,其施加给组织的压缩操作是一个连续过程,如果在一次连续压缩中,组织的形变量相差太大,则弹性应变信息差别也太大,将导致相邻多帧弹性图像之间差异较大,图像不稳定。因此,本实施例将形变程度视为考核各帧弹性图像的参数之一。
形变程度参数为实时计算出的当前帧弹性图像对应的平均应变值,即取出当前帧的感兴趣区域ROI内或者整个扫描平面区域内的各采样位置的应变数据,求得平均值,则可得到平均应变值Str_mean。平均应变值Str_mean的大小(即Str_mean的绝对值)若落在系统规定的范围内(如Str_mean小于根据经验设定的系统预设阈值),则表示形变程度合适。
(2)互相关检测质量参数
由于弹性检测模块可基于相邻两帧超声回波信号间的互相关性来检测位移信息,然后对位移信息求纵向梯度得到应变信息;因此,位移信息的准确性影响应变信息的准确性,从而会影响最终弹性图像的信噪比、对比度等。如果两帧信号互相关性较大,则检测信噪比更高,检测结果更为准确;如果两帧信号几乎不相关,则检测结果不准确。基于此,本实施例将互相关检测质量视为考核各帧弹性图像的参数之一,互相关检测质量参数为根据弹性信息检测模块采用的位移检测算法所选择相应的评分标准得到的当前帧图像的评分。
位移检测时,对于其中一帧超声回波信号中某个采样位置的信号,需要在另一帧超声回波信号中某个搜索区域内寻找与之最为相关的位置,以采用SAD作为互相关判断为例,该最为相关的位置即为SAD值最小的位置,该位置相对于原采样位置的差异即为该采样位置的位移值,类似图像匹配算法。可以理解,对于采用如SSD作为互相关判断时,该最为相关的位置即为SSD值最小的位置,而对于以互相关系数(CC,Correlation Coefficient)作为互相关判断时,该最为相关的位置即为CC值最大的位置。
这里以SAD作为互相关判断为例进行互相关检测质量参数的说明。对于每帧信号的每个采样位置,记录搜索区域内各处对应的SAD值的最大值SAD_max和最小值SAD_min,计算搜索区域的质量评分,计算步骤(即评分标准)包括:
I.系统预先设置SAD分布上下限,即为[SAD_Low,SAD_High],SAD_Low<SAD_High;
II.计算第一项得分score1,其为[0,1]之间的值,用来评估当前搜索区域中某处的SAD最大值与上限的距离,距离越近,得分越高。比如,令:score1=(SAD_max-SAD_min)/(SAD_High-SAD_min);
III.计算第二项得分score2,其为[0,1]之间的值,用来评估当前搜索区域中某处SAD最小值与下限的距离,距离越近,得分越高。比如,令score2=(SAD_max-SAD_min)/(SAD_max-SAD_Low);
IV.取score1和score2之间的加权结果作为本次搜索的质量得分score_SAD,比如令:score_SAD=score1*p+score2*(1-p),其中p为系统预先设定的参数,p的取值在0~1之间。该加权结果为[0,1]之间的值,再将score_SAD乘以100,拉伸到[0,100]之间。当然,质量得分也可不进行拉伸或拉伸至其他区间范围内,视使用习惯而定。
V.对当前帧信号的所有采样位置的质量得分求平均值,得到该帧的最终质量评分Score_mean,分值越高,意味着搜索质量越好。系统可预设分数阈值,若评分高于该阈值,则认为该帧信号的位移检测满足系统要求。
上述为基于SAD进行位移检测的方法,本领域技术人员从前述的描述可以理解,还可以根据实际选用的其它位移检测算法来选择相应的评分方法以对互相关检测质量进行评分。前述关于评分的具体计算是为了清楚说明本实施例的目的是要得到互相关检测质量的评分,而不是对本申请的限定。此外,前述所提及的如系统设定分数阈值、SAD分布上下限、系统预先设定的参数等可以是超声系统默认自动设定,也可以是由用户根据需要通过用户界面直接设定。
本实施例可以选用形变程度参数和互相关检测质量参数中的任意项来判断当前帧信号的质量是否满足预设质量要求,即,计算出的Str_mean的绝对值落在系统规定的某范围内则认为满足预设质量要求,或者计算出的Score_mean值高于系统规定的某分数阈值则认为满足预设质量要求,或者计算出的Str_mean的绝对值落在系统规定的某范围内的同时Score_mean值高于系统规定的某分数阈值,则认为该帧信号满足预设质量要求。
上述参数计算完毕后,送入显示模块。显示模块在同一画面中显示弹性图像及其相应的压缩参数和/或质量评估参数。如果选择显示弹性图像及其相应的压缩参数,则显示方式类似实施例1相应部分的描述,不再赘述。
如果选择显示弹性图像及其相应的质量评估参数,弹性图像的显示类似实施例1相应部分的描述,不再赘述,而质量评估参数的显示可以是:将满足预设质量要求的质量评估参数以同一种颜色或图案的条形图或同一种颜色或形状的线条显示,将不满足预设质量要求的质量评估参数以另一种颜色或图案的条形图或另一种颜色或形状的线条显示。例如,如图6所示,在同一画面600中除了显示弹性图像601外,质量评估参数可显示为:用同一种颜色或图案(如图6的A图中的斜线图案612)的条形图或同一种颜色或形状(如图6的B图中的虚线形状613)的线条表示满足预设质量要求的质量评估参数,用另一种颜色或图案(如其它不同于斜线图案612的图案)的条形图或同一种颜色或形状(如实线形状)的线条表示不满足预设质量要求的质量评估参数。
如果选择显示弹性图像及其相应的压缩参数和质量评估参数,弹性图像的显示类似实施例1相应部分的描述,不再赘述,而压缩参数和质量评估参数的显示可以通过同一个二维坐标图进行显示,二维坐标图的横坐标为时间,纵坐标为与时间对应的压缩参数的数值,压缩参数显示为与该压缩参数的数值对应的长度的条形图或线条,质量评估参数显示为该条形图的颜色或图案或该线条的形状。如图7的A图中,同一画面700同时显示某帧弹性图像701及其压缩参数和质量评估参数,其中压缩参数采用一条条形图702表示,而质量评估参数因满足预设质量要求而表示为斜线图案712,或者压缩参数采用一条线条703表示,而质量评估参数则表示虚线形状713。
本实施例在同时显示压缩参数和质量评估参数时,还可以作进一步改进,例如,在显示当前最新一帧弹性图像对应的压缩和质量的同时,保存之前一段固定时间内的压缩和质量情况,这样,用户既可以实时看到当前最新一帧弹性图像对应压缩质量,也可以看到之前一段固定时间内的压缩质量情况,例如还可在二维坐标图中设有帧指示标志用于表示系统指定的某帧图像(如当前帧)所处时间的位置;又例如,还可以在二维坐标图中预设一定的高度范围,当表示压缩参数的条形图或线条处于该预设高度范围内,则表示该压缩参数对应为合适的压缩程度。
本实施例的超声弹性成像系统还可以作进一步改进,可还包括告警模块,用于判断连续预设时间内各帧弹性图像的压缩参数是否满足预设合适压缩强度、以及判断连续预设时间内各帧弹性图像的质量评估参数是否满足预设质量要求,,如果均不满足则通过图像和/或声音和/或文字提示用户应当调整探头对待检测目标的压缩力度和频率。
本实施例中,显示模块中弹性图像及其压缩参数和质量评估参数的显示随着图像的实时采集同步更新,这样,用户可以实时看到当前最新一帧弹性图像对应的压缩情况和所采集的图像的质量情况,从而用户可以判断是否需要调节手法重新采集图像,尽量将压缩程度和图像质量控制在最优范围内;最终,给用户尤其是操作经验不足的用户提供一定的指导,使得操作更简便。
本申请中超声弹性成像方法的一种实施例与上述超声弹性成像系统的实施例2相应,包括:
步骤21、在弹性成像模式下,探头以系统预先设定好的扫描规则进行超声发射并接收超声回波信号;
步骤22、对接收到的超声回波信号进行信号预处理,预处理包括波束合成等;
步骤23、对波束合成输出的信号进行处理,计算反映待检测目标的弹性的物理量,根据该物理量生成相应的弹性图像;
步骤24、计算与弹性图像对应的压缩参数,压缩参数用于反映探头对待检测目标的压缩程度;还计算反映弹性图像质量的质量评估参数;
步骤25、显示弹性图像及其对应的压缩参数和质量评估参数。
以上各步骤的具体实现可参考上述超声弹性成像系统的实施例2中各模块对应的实现过程,在此不再重述。此外上述方法实施例还可以包括对B信号进行处理以形成待检测目标的灰阶图像的步骤。
实施例3
本实施例的超声弹性成像系统包括:超声探头及其发射接收模块、信号预处理模块、B信号处理模块、弹性处理模块和显示模块,其中弹性处理模块包括弹性检测模块和参数计算模块。其中,探头及其发射接收模块、信号预处理模块和B信号处理模块、弹性检测模块分别与实施例1的探头及其发射接收模块、信号预处理模块、B信号处理模块、弹性检测模块类似,不再赘述。参数计算模块仅计算反映弹性图像质量的质量评估参数,显示模块显示弹性图像及其对应的质量评估参数。具体的质量评估参数的计算及其显示可参考实施例2中相应部分的描述,不再赘述。
本实施例中,显示模块中弹性图像及其质量评估参数的显示随着图像的实时采集同步更新,这样,用户可以实时看到当前最新一帧弹性图像对应的图像质量情况,从而用户可以判断是否需要重新采集图像,尽量将图像质量控制在最优范围内;最终,给用户尤其是操作经验不足的用户提供一定的指导,使得操作更简便。
下面通过一个实例来说明进一步说明对压缩参数和质量评估参数的显示。
如图8所示,压缩参数和质量评估参数通过一个二维坐标图显示,横坐标对应一定长度的时间,该时间长度由系统指定,纵坐标对应组织相对初始位置的平均偏移量。对于每帧弹性图像,使用一个条形图(bar)来评价其压缩与质量。压缩参数(即前面计算得到的累计应变量Str_integral)对应于与该bar的高度。该Str_intergral有正负区别,表示组织相对于初始位置有上下偏移,指示bar也会分布在坐标轴的上下两侧。当然,也可以直接采用Str_integral的绝对值来与bar的高度对应,此时bar的高度仅表示组织相对初始位置的偏移量,而不指明偏移方向,此时指示bar均分布在坐标轴的同一侧。系统还可指定坐标轴上特定高度范围内为合适的压缩程度,比如图8中的虚线框区域,以帮助用户调节力度,尽量将压缩程度控制在最优范围内。质量评估参数选用形变程度参数和互相关检测质量参数,形变程度Str_mean的绝对值落在系统指定的范围内,同时,互相关检测质量评分Score_mean高于系统指定的分值,则将该帧对应的小bar的颜色或图案指定为某种特殊颜色或图案(比如图8中的斜线图案);否则,则将小bar的颜色或图案指定为其他颜色或图案(比如图8中的无图案bar)。
显示随着图像的实时采集同步更新,这样用户既可以实时看到当前最新一帧弹性图像对应压力质量,也可以看到之前一段固定时间内的压力情况。最终,根据一段时间内小bar的颜色或图案分布情况,用户可初步判断图像质量是否满足系统要求,是否需要调节手法重新采集图像。比如操作后,bar的颜色或图案绝大多数为系统不建议的颜色或图案,则说明所得图像可信度降低,用户手法需要调整;反之,如果一定时间内bar的颜色或图案绝大多数为系统建议的颜色或图案,则说明所得图像可信度较高。根据一段时间内小bar的高度变化,用户可根据需要来适度增大或减小压缩组织的幅度,或者适度加快或减慢压放组织的频率。
成功采集到满足需求的图像后,在用户回放观察已采集图像的同时,二维坐标图中坐标轴随之同步滚动刷新,通过系统设定的帧指示标志,用户可识别出与当前观察帧所对应的bar(比如图8中的小三角标志)。即用户可方便的回放观察已采集的图像及其对应的压力bar。
综上,当用户实时采集弹性数据时,使用一个二维坐标平面上的条形图bar(如彩色bar)来提示压力质量,每获得一帧弹性图像,即在该二维坐标平面上当前时间位置实时显示一个bar,该bar随采集过程实时更新,且与当前用户所见的弹性图像保持同步。该坐标平面对应的时间长度由系统指定。成功采集到满足需求的图像后,在用户回放观察已采集图像的同时,坐标轴随之同步滚动刷新,通过系统设定的帧指示标志,用户可识别出与当前观察帧所对应的bar。彩色bar的纵向长度(或高度)表示压缩程度的大小,通过实时计算出压放过程中的累积应变量Str_integral参数来指示压缩程度。Str_integral有正有负,则代表了组织相对于初始位置的偏移量有上有下。也可采用Str_integral的绝对值来表示压缩程度,则其高度直接表示组织相对于初始位置偏移量的大小,不区分偏移方向。系统还可指定一定的bar高度建议范围,如果一定时间内绝大多数的bar高度均落在建议范围内,说明用户压放大小合适,否则需要适当增加或减小压放幅度。上述彩色bar的颜色表示对图像质量的评价,通过形变程度和互相关检测质量两方面来评估图像质量是否满足系统要求。实时计算出当前帧对应的平均应变值Str_mean值,通过判断其大小(即Strain_mean的绝对值)是否落在系统指定的范围内来表示其形变程度是否合适;实时计算出当前帧对应的位移检测评分Score_mean,通过判断其是否高于系统指定的分数来表示其位移检测是否可信。如果两方面均满足要求,则将该帧对应的bar的颜色设置为特定颜色(比如绿色,由系统指定),否则,则设置为其他颜色。如果一定连续时间内,绝大部分bar的颜色均为系统建议的颜色,则说明图像质量满足系统要求,否则用户需要调整压放的幅度或频率,重新采集图像。
按本申请提供的系统/方法,用户在采集弹性图像过程中,可以根据bar的颜色或图案来判断图像质量是否合适,压缩是否需要调整,根据坐标平面上bar的高度来调整探头压缩组织的力度,根据一定时间长度内bar形成的轮廓变化来调整压缩频率。最终,给用户尤其是首次操作的用户提供一定的指导,使得操作更简便。
本领域技术人员可以理解,上述实施方式中各种方法的全部或部分步骤可以通过程序来指令相关硬件完成,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:只读存储器、随机存储器、磁盘或光盘等。
以上内容是结合具体的实施方式对本申请所作的进一步详细说明,不能认定本申请的具体实施只局限于这些说明。对于本申请所属技术领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干简单推演或替换。
Claims (14)
1.一种超声弹性成像系统,其特征在于,包括:
弹性检测模块,用于对波束合成输出的信号进行处理,计算反映待检测目标的弹性的物理量,根据该物理量生成相应的弹性图像;
参数计算模块,用于计算与所述弹性图像对应的压缩参数,所述压缩参数用于反映探头对待检测目标的压缩程度;
显示模块,用于显示弹性图像及其对应的压缩参数;
其中,所述参数计算模块还用于计算反映弹性图像质量的质量评估参数,所述质量评估参数包括互相关检测质量参数,所述互相关检测质量参数为根据采用的位移检测算法所对应的评分标准得到的弹性图像的评分,所述采用的位移检测算法为绝对值和SAD算法,所述对应的评分标准包括:
计算第一项得分,所述第一项得分用于评估位移检测中当前搜索区域的SAD最大值和预设的SAD上限的距离;
计算第二项得分,所述第二项得分用于评估所述当前搜索区域的SAD最小值与预设的SAD下限的距离;
取第一项得分和第二项得分的加权结果作为当前搜索的质量的分;
对弹性图像遍历得到的质量得分求平均值即得到该弹性图像的评分。
2.如权利要求1所述的超声弹性成像系统,其特征在于,所述显示模块还用于显示与弹性图像对应的质量评估参数。
3.如权利要求2所述的超声弹性成像系统,其特征在于,
所述显示模块在显示弹性图像的同时,在同一画面中显示该弹性图像的压缩参数或质量评估参数,所述压缩参数的显示包括:将压缩参数显示为与该压缩参数的数值对应的长度的条形图或线条,或者将压缩参数显示为用针的旋转角表示所述压缩参数的数值的模拟仪表,或者通过二维坐标图进行显示,所述二维坐标图的横坐标为时间,纵坐标为与时间对应的压缩参数的数值,所述压缩参数显示为与该压缩参数的数值对应的长度的条形图或线条;所述质量评估参数的显示包括:将满足预设质量要求的质量评估参数以同一种颜色或图案的条形图或同一种颜色或形状的线条显示,将不满足预设质量要求的质量评估参数以另一种颜色或图案的条形图或另一种颜色或形状的线条显示;
或者,所述显示模块在显示弹性图像的同时,在同一画面中显示该弹性图像的压缩参数和质量评估参数,所述压缩参数和质量评估参数的显示通过同一个二维坐标图进行显示,所述二维坐标图的横坐标为时间,纵坐标为与时间对应的压缩参数的数值,所述压缩参数显示为与该压缩参数的数值对应的长度的条形图或线条,所述质量评估参数显示为所述条形图的颜色或图案或所述线条的颜色或形状。
4.如权利要求3所述的超声弹性成像系统,其特征在于,所述二维坐标图中设有帧指示标志用于表示系统指定的某帧图像所处时间的位置;和/或,所述二维坐标图中设有预设高度范围,显示所述压缩参数的条形图或线条处于所述预设高度范围内表示所述压缩参数为合适的压缩程度。
5.如权利要求4所述的超声弹性成像系统,其特征在于,还包括:告警模块,用于判断连续预设时间内各帧弹性图像的压缩参数是否满足预设合适压缩强度,如果不满足则通过图像和/或声音和/或文字进行提示,所述告警模块还用于判断连续预设时间内各帧弹性图像的质量评估参数是否满足预设质量要求,如果不满足则通过图像和/或声音和/或文字进行提示。
6.如权利要求1-5任一项所述的超声弹性成像系统,其特征在于,所述压缩参数包括累积应变量,所述累积应变量为从第一帧弹性图像的应变数据至当前帧弹性图像的应变数据的积分;
所述第一帧弹性图像的应变数据为第一帧弹性图像中感兴趣区域内各采样位置的应变量的平均值,所述当前帧弹性图像的应变数据为当前帧弹性图像中感兴趣区域内各采样位置的应变量的平均值;
或者,所述第一帧弹性图像的应变数据为第一帧弹性图像中感兴趣区域内某一深度的所有位移量的平均值与该深度的比值,所述当前帧弹性图像的应变数据为当前帧弹性图像中感兴趣区域内某一深度的所有位移量的平均值与该深度比值。
7.如权利要求6所述的超声弹性成像系统,其特征在于,所述质量评估参数还包括形变程度;所述形变程度参数为弹性图像中感兴趣区域内各采样位置的应变量的平均值。
8.一种超声弹性成像方法,其特征在于,包括:
弹性检测步骤,对波束合成输出的信号进行处理,计算反映待检测目标的弹性的物理量,根据该物理量生成相应的弹性图像;
参数计算步骤,计算与所述弹性图像对应的压缩参数,所述压缩参数用于反映探头对待检测目标的压缩程度;
显示步骤,显示弹性图像及其对应的压缩参数;
其中,所述参数计算步骤还计算反映弹性图像质量的质量评估参数,所述质量评估参数包括互相关检测质量参数,所述互相关检测质量参数为根据采用的位移检测算法所对应的评分标准得到的弹性图像的评分,所述采用的位移检测算法为绝对值和SAD算法,所述对应的评分标准包括:
计算第一项得分,所述第一项得分用于评估位移检测中当前搜索区域的SAD最大值和预设的SAD上限的距离;
计算第二项得分,所述第二项得分用于评估所述当前搜索区域的SAD最小值与预设的SAD下限的距离;
取第一项得分和第二项得分的加权结果作为当前搜索的质量的分;
对弹性图像遍历得到的质量得分求平均值即得到该弹性图像的评分。
9.如权利要求8所述的超声弹性成像方法,其特征在于,所述显示步骤还显示与弹性图像对应的质量评估参数。
10.如权利要求9所述的超声弹性成像方法,其特征在于,
在显示弹性图像的同时,在同一画面中显示该弹性图像的压缩参数或质量评估参数,所述压缩参数的显示包括:将压缩参数显示为与该压缩参数的数值对应的长度的条形图或线条,或者将压缩参数显示为用针的旋转角表示所述压缩参数的数值的模拟仪表,或者通过二维坐标图进行显示,所述二维坐标图的横坐标为时间,纵坐标为与时间对应的压缩参数的数值,所述压缩参数显示为与该压缩参数的数值对应的长度的条形图或线条;所述质量评估参数的显示包括:将满足预设质量要求的质量评估参数以同一种颜色或图案的条形图或同一种颜色或形状的线条显示,将不满足预设质量要求的质量评估参数以另一种颜色或图案的条形图或另一种颜色或形状的线条显示;
或者,在显示弹性图像的同时,在同一画面中显示该弹性图像的压缩参数和质量评估参数,所述压缩参数和质量评估参数的显示通过同一个二维坐标图进行显示,所述二维坐标图的横坐标为时间,纵坐标为与时间对应的压缩参数的数值,所述压缩参数显示为与该压缩参数的数值对应的长度的条形图或线条,所述质量评估参数显示为所述条形图的颜色或图案或所述线条的颜色或形状。
11.如权利要求10所述的超声弹性成像方法,其特征在于,所述二维坐标图中设有帧指示标志用于表示系统指定的某帧图像所处时间的位置;和/或,所述二维坐标图中设有预设高度范围,显示所述压缩参数的条形图或线条处于所述预设高度范围内表示所述压缩参数为合适的压缩程度。
12.如权利要求11所述的超声弹性成像方法,其特征在于,还包括:告警步骤,判断连续预设时间内各帧弹性图像的压缩参数是否满足预设合适压缩强度,如果不满足则通过图像和/或声音和/或文字进行提示;所述告警步骤还判断连续预设时间内各帧弹性图像的质量评估参数是否满足预设质量要求,如果不满足则通过图像和/或声音和/或文字进行提示。
13.如权利要求8-12任一项所述的超声弹性成像方法,其特征在于,所述压缩参数包括累积应变量,所述累积应变量为从第一帧弹性图像的应变数据至当前帧弹性图像的应变数据的积分;
所述第一帧弹性图像的应变数据为第一帧弹性图像中感兴趣区域内各采样位置的应变量的平均值,所述当前帧弹性图像的应变数据为当前帧弹性图像中感兴趣区域内各采样位置的应变量的平均值;
或者,所述第一帧弹性图像的应变数据为第一帧弹性图像中感兴趣区域内某一深度的所有位移量的平均值与该深度的比值,所述当前帧弹性图像的应变数据为当前帧弹性图像中感兴趣区域内某一深度的所有位移量的平均值与该深度比值。
14.如权利要求13所述的超声弹性成像方法,其特征在于,所述质量评估参数还包括形变程度;所述形变程度参数为弹性图像中感兴趣区域内各采样位置的应变量的平均值。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210495070.XA CN103845074B (zh) | 2012-11-28 | 2012-11-28 | 一种超声弹性成像系统和方法 |
PCT/CN2013/083878 WO2014082482A1 (zh) | 2012-11-28 | 2013-09-22 | 一种超声弹性成像系统和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210495070.XA CN103845074B (zh) | 2012-11-28 | 2012-11-28 | 一种超声弹性成像系统和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103845074A CN103845074A (zh) | 2014-06-11 |
CN103845074B true CN103845074B (zh) | 2017-12-12 |
Family
ID=50827153
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210495070.XA Active CN103845074B (zh) | 2012-11-28 | 2012-11-28 | 一种超声弹性成像系统和方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN103845074B (zh) |
WO (1) | WO2014082482A1 (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106572838B (zh) * | 2014-09-03 | 2019-09-06 | 深圳迈瑞生物医疗电子股份有限公司 | 弹性测量检测方法及系统 |
CN118000784A (zh) * | 2016-08-01 | 2024-05-10 | 深圳迈瑞生物医疗电子股份有限公司 | 剪切波弹性成像测量显示方法及系统 |
CN114360727A (zh) * | 2017-04-21 | 2022-04-15 | 深圳迈瑞生物医疗电子股份有限公司 | 超声弹性成像装置及弹性成像结果评价方法 |
JP6810005B2 (ja) * | 2017-09-14 | 2021-01-06 | 株式会社日立製作所 | 超音波診断装置 |
CN107833218B (zh) * | 2017-11-22 | 2021-01-26 | 深圳中科乐普医疗技术有限公司 | 一种生物组织剪切波弹性图像质量评估显示方法 |
CN108175440A (zh) * | 2017-12-21 | 2018-06-19 | 飞依诺科技(苏州)有限公司 | 一种用于超声扫描设备的弹性成像方法和装置 |
CN112534468A (zh) * | 2018-08-24 | 2021-03-19 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性成像装置及对弹性图像进行处理的方法 |
CN109978823B (zh) * | 2019-02-15 | 2021-08-06 | 无锡海斯凯尔医学技术有限公司 | 弹性成像图像分析及组织粘弹性检测方法及装置 |
CN117224161A (zh) * | 2019-12-06 | 2023-12-15 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像方法以及超声成像系统 |
CN113397588A (zh) * | 2020-03-16 | 2021-09-17 | 深圳市理邦精密仪器股份有限公司 | 弹性成像方法、装置及医疗设备 |
CN113476075A (zh) * | 2020-03-16 | 2021-10-08 | 深圳市理邦精密仪器股份有限公司 | 一种超声弹性成像方法、图像数据的筛选方法及装置 |
CN113040816A (zh) * | 2021-04-06 | 2021-06-29 | 无锡海斯凯尔医学技术有限公司 | 超声弹性成像方法、装置、电子设备及存储介质 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6511427B1 (en) * | 2000-03-10 | 2003-01-28 | Acuson Corporation | System and method for assessing body-tissue properties using a medical ultrasound transducer probe with a body-tissue parameter measurement mechanism |
JP4465535B2 (ja) * | 2004-06-09 | 2010-05-19 | 株式会社日立メディコ | 弾性画像表示方法及び超音波診断装置 |
US8353831B2 (en) * | 2004-06-22 | 2013-01-15 | Hitachi Medical Corporation | Diagnostic ultrasound system and method of displaying elasticity image |
US20090177084A1 (en) * | 2006-05-25 | 2009-07-09 | Takeshi Matsumura | Ultrasonic Diagnostic Apparatus |
JP5199690B2 (ja) * | 2008-02-07 | 2013-05-15 | 株式会社日立メディコ | 超音波診断装置 |
JP4999969B2 (ja) * | 2010-07-13 | 2012-08-15 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置及びその制御プログラム |
US20130158900A1 (en) * | 2010-08-31 | 2013-06-20 | Hitachi Medical Corporation | Ultrasonic diagnostic apparatus and evaluation calculation method |
-
2012
- 2012-11-28 CN CN201210495070.XA patent/CN103845074B/zh active Active
-
2013
- 2013-09-22 WO PCT/CN2013/083878 patent/WO2014082482A1/zh active Application Filing
Also Published As
Publication number | Publication date |
---|---|
WO2014082482A1 (zh) | 2014-06-05 |
CN103845074A (zh) | 2014-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103845074B (zh) | 一种超声弹性成像系统和方法 | |
CN103845081B (zh) | 超声弹性成像系统和方法、实时动态帧间处理方法 | |
CN108733857B (zh) | 超声弹性成像装置及弹性成像结果评价方法 | |
KR100893286B1 (ko) | 영상 처리 시스템 및 방법 | |
EP3389498B1 (en) | Fetal monitoring system and method | |
US8144961B2 (en) | Ultrasound diagnostic apparatus and method for measuring a size of a target object | |
US5495852A (en) | Method and apparatus for estimating diameter of an artery using B-mode ultrasonic images | |
EP1889571A1 (en) | Ultrasonograph | |
WO2011093193A1 (ja) | 超音波診断装置およびその計測点追跡方法 | |
US20190159762A1 (en) | System and method for ultrasound elastography and method for dynamically processing frames in real time | |
JP2009207899A (ja) | 超音波映像を処理する超音波システム及び方法 | |
JP2009540915A (ja) | アレイ圧力センサーを利用した脈波分析装置 | |
US10937546B2 (en) | Degree-of-distraction estimation system based on unconstrained bio-information | |
US20150272541A1 (en) | Ultrasonic measuring apparatus and ultrasonic measuring method | |
US20150238169A1 (en) | Ultrasonic measurement apparatus and ultrasonic measurement method | |
CN102930555A (zh) | 一种对超声图像中感兴趣区域进行跟踪的方法及装置 | |
CN108403101A (zh) | 一种模型训练方法及装置、测量血压的方法及装置 | |
CN104739442A (zh) | 压力弹性成像位移检测方法、装置和超声成像设备 | |
US20080139936A1 (en) | Ultrasound diagnostic system and method of displaying multiple growth trend charts | |
US11457823B2 (en) | Wearable blood pressure detecting device and detecting method thereof | |
CN113456033A (zh) | 生理指标特征值数据处理方法、系统及计算机设备 | |
WO2013061911A1 (ja) | 測定装置、測定方法および測定プログラムを格納した記録媒体 | |
CN114652353A (zh) | 超声成像系统以及颈动脉斑块稳定性评估方法 | |
CN112534468A (zh) | 一种超声弹性成像装置及对弹性图像进行处理的方法 | |
JP2010068955A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20140611 Assignee: Shenzhen Mindray Animal Medical Technology Co.,Ltd. Assignor: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd. Contract record no.: X2022440020009 Denomination of invention: Ultrasound elastography system and method Granted publication date: 20171212 License type: Common License Record date: 20220804 |