CN107811652B - 自动调整参数的超声成像方法及系统 - Google Patents

自动调整参数的超声成像方法及系统 Download PDF

Info

Publication number
CN107811652B
CN107811652B CN201710972718.0A CN201710972718A CN107811652B CN 107811652 B CN107811652 B CN 107811652B CN 201710972718 A CN201710972718 A CN 201710972718A CN 107811652 B CN107811652 B CN 107811652B
Authority
CN
China
Prior art keywords
current
time window
current time
time
scanning
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
CN201710972718.0A
Other languages
English (en)
Other versions
CN107811652A (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.)
Feiyinuo Technology Co ltd
Original Assignee
Vinno Technology Suzhou 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 Vinno Technology Suzhou Co Ltd filed Critical Vinno Technology Suzhou Co Ltd
Priority to CN201710972718.0A priority Critical patent/CN107811652B/zh
Publication of CN107811652A publication Critical patent/CN107811652A/zh
Application granted granted Critical
Publication of CN107811652B publication Critical patent/CN107811652B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42Details of probe positioning or probe attachment to the patient
    • A61B8/4245Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient
    • A61B8/4254Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient using sensors mounted on the probe
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54Control of the diagnostic device

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

自动调整参数的超声成像方法及系统
技术领域
本发明涉及超声领域,尤其涉及一种自动调整参数的超声成像方法及系统。
背景技术
随着电子学、计算机、材料科学等相关领域技术的发展;超声成像因其无创性、实时性、操作方便、价格便宜等诸多优势,使其成为临床上应用最为广泛的辅助诊断的手段之一;
在现有超声成像系统中,针对不同的应用场景通常需要应用不同的成像参数来获取更好的图像效果。
例如:在特定的应用环境下,使用同一探头对变化速度较慢的组织器官:如肌肉,骨骼,神经,静脉,腹腔器官,和对变化速度较快的组织器官:如心脏,动脉,进行成像时,采用的成像策略、参数均不相同。
为了获得精度更高的超声图像,现有技术通常由操作者在操作超声成像系统进行扫查前、手动通过用户界面对场景和参数的进行选择和调整。
发明内容
本发明的目的在于提供一种自动调整参数的超声成像方法及系统。
为实现上述目的之一,本发明一实施方式的自动调整参数的超声成像方法,所述方法包括:
M1、设定时间窗口的间隔时间长度阈值;
M4、通过超声探头扫查成像过程中,实时获取每个时刻下的回波信号;
实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的回波信号,判断当前时间窗口内的回波信号是否发生相对变化;
若是,调整当前成像参数,以使当前的扫查模式调整为动态扫查模式;
若否,调整当前成像参数,以使当前的扫查模式调整为静态扫查模式。
作为本发明一实施方式的进一步改进,步骤M1和步骤M4之间还包括:
M2、通过超声探头扫查成像过程中,实时获取每个时刻下探头的运动状态;
M3、实时以当前时刻作为时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的探头运动状态,判断当前时间窗口内探头的运动状态是否发生相对变化;
若是,保持当前的扫查参数继续扫查;
若否,进入步骤M4。
作为本发明一实施方式的进一步改进,所述超声探头中集成运动传感器,以用于实时获取超声探头的运动状态;所述运动传感器为:加速度传感器、陀螺仪或磁强计中的至少一种。
作为本发明一实施方式的进一步改进,所述运动传感器为三轴加速度传感器,
所述步骤M3具体包括:实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的加速度矢量判断当前时间窗口内探头的运动状态;当三轴加速度传感器的加速度矢量满足以下条件至少其中之一时,确认所述探头在当前时间窗口内的运动状态发生变化;
条件1:|Gt|-|Gs|>静态阈值;
条件2:|Gt-Gt-t0|>动态阈值;
其中,Gt表示当前时刻测得的加速度矢量,t表示当前时刻,Gt-t0表示在t-t0时刻测得的加速度矢量,t0的值小于等于时间窗口对应的间隔时间长度阈值;Gs表示当前地理位置标定的静止状态下的加速度矢量;静态阈值、动态阈值均为常数。
作为本发明一实施方式的进一步改进,所述方法还包括:采用SAD算法、SSD算法、Census变换匹配算法、图像相关性算法、SIFT特征识别算法中的至少一种算法判断当前时间窗口内的回波信号是否发生相对变化;
其中,所述图像相关性算法为:
以当前时刻作为当前时间窗口的结束时刻,获取当前时刻以及其上一时刻的帧图像中每根扫查线所对应的相关系数,并根据其判断当前时间窗口内的回波信号是否发生相对变化;
Figure GDA0002532103450000031
X表示当前时刻帧图像的其中一条扫查线,Y表示上一时刻帧图像中对应X的扫查线r(X,Y)表示对应扫查线的相关系数,Cov(X,Y)表示X与Y的协方差,Var[X]为X的方差,Var[Y]为Y的方差。
为了实现上述发明目的另一,本发明一实施方式提供一种自动调整参数的超声成像系统,所述系统包括:计数模块,用于设定时间窗口的间隔时间长度阈值;
获取模块,用于在通过超声探头扫查成像过程中,实时获取每个时刻下的回波信号;
处理模块,用于实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的回波信号,判断当前时间窗口内的回波信号是否发生相对变化;
若是,调整当前成像参数,以使当前的扫查模式调整为动态扫查模式;
若否,调整当前成像参数,以使当前的扫查模式调整为静态扫查模式。
作为本发明一实施方式的进一步改进,所述处理模块还用于:通过超声探头扫查成像过程中,实时获取每个时刻下探头的运动状态;
实时以当前时刻作为时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的探头运动状态,判断当前时间窗口内探头的运动状态是否发生相对变化;
若是,保持当前的扫查参数继续扫查;
若否,实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的回波信号,判断当前时间窗口内的回波信号是否发生相对变化。
作为本发明一实施方式的进一步改进,所述超声探头中集成运动传感器,以用于实时获取超声探头的运动状态;所述运动传感器为:加速度传感器、陀螺仪或磁强计中的至少一种。
作为本发明一实施方式的进一步改进,所述运动传感器为三轴加速度传感器,
所述处理模块还用于:实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的加速度矢量判断当前时间窗口内探头的运动状态;当三轴加速度传感器的加速度矢量满足以下条件至少其中之一时,确认所述探头在当前时间窗口内的运动状态发生变化;
条件1:|Gt|-|Gs|>静态阈值;
条件2:|Gt-Gt-t0|>动态阈值;
其中,Gt表示当前时刻测得的加速度矢量,t表示当前时刻,Gt-t0表示在t-t0时刻测得的加速度矢量,t0的值小于等于时间窗口对应的间隔时间长度阈值;Gs表示当前地理位置标定的静止状态下的加速度矢量;静态阈值、动态阈值均为常数。
作为本发明一实施方式的进一步改进,所述处理模块采用SAD算法、SSD算法、Census变换匹配算法、图像相关性算法、SIFT特征识别算法中的至少一种算法判断当前时间窗口内的回波信号是否发生相对变化;
其中,所述图像相关性算法为:
以当前时刻作为当前时间窗口的结束时刻,获取当前时刻以及其上一时刻的帧图像中每根扫查线所对应的相关系数,并根据其判断当前时间窗口内的回波信号是否发生相对变化;
Figure GDA0002532103450000051
X表示当前时刻帧图像的其中一条扫查线,Y表示上一时刻帧图像中对应X的扫查线r(X,Y)表示对应扫查线的相关系数,Cov(X,Y)表示X与Y的协方差,Var[X]为X的方差,Var[Y]为Y的方差。
与现有技术相比,本发明的有益效果是:本发明的自动调整参数的超声成像方法及系统,能够适时对成像参数进行调整,简化用户操作,提升在系统在单一工作状态下的适应能力和图像质量。
附图说明
图1是本发明一实施方式中自动调整参数的超声成像方法的流程示意图。
图2是本本发明一实施方式中自动调整参数的超声成像系统的模块示意图。
具体实施方式
以下将结合附图所示的具体实施方式对本发明进行详细描述。但这些实施方式并不限制本发明,本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的变换均包含在本发明的保护范围内。
通过探头扫查以进行超声成像过程中,通常情况下,其扫查模式为固定扫查模式,该固定扫查模式下,扫查过程中,探头的参数固定不变;而本发明应用的扫查模式为自适应扫查模式,该自适应扫查模式下,扫查过程中,探头的参数会随着探头以及回波信号的改变而进行自适应的调整,如此,以使扫查结果更加精准。
如图1所示,本发明一实施方式中,自动调整参数的超声成像方法,所述方法包括:
M1、设定时间窗口的间隔时间长度阈值;
时间窗口在本申请中为一时间监测周期,实际应用中为一时间固定值,其大小可根据需要具体调整,本发明中其值大小以间隔时间长度阈值表示;由于超声的扫查精度及对时限精度要求较高,通常情况下,该时间窗口的间隔时间长度阈值以秒为单位,例如:3s、5s。
M2、通过超声探头扫查成像过程中,实时获取每个时刻下探头的运动状态;
M3、实时以当前时刻作为时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的探头运动状态,判断当前时间窗口内探头的运动状态是否发生相对变化;
若是,保持当前的扫查参数继续扫查;
若否,进入下一步骤。
该步骤中,可通过多种方式获得探头的运动状态;本发明具体实施方式中,在所述超声探头中集成运动传感器,以用于实时获取超声探头的运动状态,并根据超声探头在各个时刻下的运动状态判断其是否发生相对变化。
所述运动传感器可为加速度传感器、陀螺仪或磁强计中的至少一种,相应的,运动传感器会随着探头的运动而改变其姿态和/或位置;所述姿态的变化例如:探头运动过程中导致的运动传感器的角度变化;所述位置的变化例如:探头运动过程中导致的运动传感器位置的移动等,在此不做具体赘述。
本发明具体实施方式中,所述运动传感器为三轴加速度传感器,通过所述三轴加速度传感器判断探头运动状态是否发生变化。
具体的,
实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的加速度矢量判断当前时间窗口内探头的运动状态;当三轴加速度传感器的加速度矢量满足以下条件至少其中之一时,确认所述探头在当前时间窗口内的运动状态发生变化;
条件1:|Gt|-|Gs|>静态阈值;
条件2:|Gt-Gt-t0|>动态阈值;
其中,Gt表示当前时刻测得的加速度矢量,t表示当前时刻,Gt-t0表示在t-t0时刻测得的加速度矢量,t0的值小于等于时间窗口对应的间隔时间长度阈值;Gs表示当前地理位置标定的静止状态下的加速度矢量;静态阈值、动态阈值均为常数。
该实施方式中,三轴加速度传感器固定设置于探头中,当判断三轴加速传感器的运动状态变化时,即表示探头的运动状态发生同步变化。
该实施方式中,当满足条件1时,三轴加速度传感器处于加速运动状态,当满足条件2时,三轴加速度传感器处于加速运动状态或其姿态迅速发生变化。
如上两个条件满足其中至少一个时,即可认为探头的运动状态发生变化,当然,如果以两个条件同时满足作为判断标准,则获得判断结果会更加精准。
该实施方式中,静态阈值、动态阈值均为一个常数,其可以根据需要具体调整,通常情况下,设备出厂时会为其指定数值;当用户需要时,还可以具体调整其大小,在此不做详细赘述。
当然,在本发明的其他实施方式中,还可以在探头上安装光学或磁场辅助定位装置以替代运动传感器,其同样可以通过运动传感器辅助检测探头运动状态的目的,在此不做详细赘述。
在固定的时间窗口内,判断探头的运动状态未发生变化后,可进一步通过回波信号判断是否需要对扫查模式进行调整;而判断探头的运动状态发生变化时,保持当前的扫查参数继续扫查;当然,在本发明其他实施方式中,当探头中未集成运动传感器时,可直接通过回波信号调整扫查模块,即在步骤M1后,直接执行步骤M4。
本发明优选实施方式中,在执行步骤M3后,再根据其结果判断是否需要对回波信号进行分析。
M4、通过超声探头扫查成像过程中,实时获取每个时刻下的回波信号;
超声成像过程中,超声探头发射声波对待测组织进行逐次的扫查,每次的扫查均有一定的方向性,且没逐次扫查过程中,均可获得多个方向的回波信号,进一步的,将逐次扫查获得的多个方向的回波信息集合到一起,就可以得到一帧扫查图像。
进一步的,所示步骤M4后,所示方法还包括:
M5、实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的回波信号,判断当前时间窗口内的回波信号是否发生相对变化;
若是,调整当前成像参数,以使当前的扫查模式调整为动态扫查模式;
若否,调整当前成像参数,以使当前的扫查模式调整为静态扫查模式。
所述动态扫查模式以及静态扫查模式均对应不同的成像参数,如此,通过实时调整成像参数,即可以自适应调整探头的扫查模式,以使获得图像更为精准。
本发明的步骤M5在实现过程中,可以采用多种方式判断当前时间窗口内的回波信号是否发生相对变化;本发明具体实施方式中,可在两个层面实现步骤M4。
扫查线层面,每条扫查线包含一次发射和接收的回波信号的包络,假设图像帧频为10帧/s,则一帧图像中的扫查线最慢可以以1/10=0.1秒的时间间隔重复,通过对比相隔0.1秒先后获得的扫查线包络信息的变化,就可以得出图像信息是否发生改变。常用的基于区域的局部匹配准则主要有图像序列中对应像素差的绝对值SAD(sum of the absolutedifferences)算法、图像序列中对应像素差的平方和SSD(sum of squared differences)算法、Census变换匹配算法、图像相关性算法。
SAD算法:
SAD算法是一种最简单的匹配算法,此种方法就是以左目图像的源匹配点为中心,定义一个窗口D,其大小为(2m+1)(2n+1),统计其窗口的灰度值的和,然后在右目图像中逐步计算其左右窗口的灰度和的差值,最后搜索到的差值最小的区域的中心像素即为匹配点。
SSD匹配算法:SSD算法与SAD算法相似,在此不做详细赘述。
Census变换匹配算法:
Census变换算法是先对单幅图像的窗口区域进行变换,然后再按照Hamming距离来比较它们的相似性进行匹配。
图像相关性算法,是计算两幅图像匹配区域的互相关性,以当前时刻作为当前时间窗口的结束时刻,获取当前时刻以及其上一时刻的帧图像中每根扫查线所对应的相关系数,并根据其判断当前时间窗口内的回波信号是否发生相对变化。
其中,相关系数以公式
Figure GDA0002532103450000101
表示,该式中,X表示当前时刻帧图像的其中一条扫查线,Y表示上一时刻帧图像中对应X的扫查线r(X,Y)表示对应扫查线的相关系数,Cov(X,Y)表示X与Y的协方差,Var[X]为X的方差,Var[Y]为Y的方差。
该相关系数的值通常介于0~1之间,当其值越接近0时,表示两幅图像相关,当其值越接近1时,表示两幅图像不相关。
进一步的,还可优选的结合当前帧图像中多根扫查线和上帧图像中多帧扫查线对比的情况,以判断回波信息是否发生大幅度的改变,使检测结果更加精准,在此不再继续赘述。
图像层面,扫查线层面的方式同样可以以二维的数学形式应用到二维的图像差异判别上,除此之外类似SIFT(scale invariant feature transform)的特征识别算法也可以被用来捕捉图像中的特征点,通过比对特征点的位置来判断图像是否发生改变。
如此之后,针对超声扫查过程中,扫查对象多变,以及用户针对不同扫查反复调整参数比较麻烦的问题;采用上述方法后,可以在超声扫查过程中,恰当的调整扫查参数,以调显著提升超声图像质量;另外,扫查参数的自动调整,能够明显简化超声设备尤其是便携式和手持式超声设备的操作流程,提升产品易用性和用户体验。
需要说明的是,本发明的优选实施方式中,所述方法的执行步骤可根据需要进行调整,其可按照步骤M1、M2、M3、M4、M5的顺序执行,也可以先执行步骤M1,之后同时执行步骤M2、M4,再根据M3的结果判断是否执行步骤M5。
结合图2所示,本发明一实施方式,提供的自动调整参数的超声成像系统包括:计数模块100、获取模块200、处理模块300。
计数模块100用于设定时间窗口的间隔时间长度阈值;
时间窗口在本申请中为一时间监测周期,实际应用中为一时间固定值,其大小可根据需要具体调整,本发明中其值大小以间隔时间长度阈值表示;由于超声的扫查精度及对时限精度要求较高,通常情况下,该时间窗口的间隔时间长度阈值以秒为单位,例如:3s、5s。
该计数模块中可以设置一计数器,以用于记录各个时间节点,在此不做详细赘述。
获取模块200用于通过超声探头扫查成像过程中,实时获取每个时刻下探头的运动状态;
处理模块300实时以当前时刻作为时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的探头运动状态,判断当前时间窗口内探头的运动状态是否发生相对变化;
若是,保持当前的扫查参数继续扫查;
若否,调整当前成像参数,以使当前的扫查模式调整为静态扫查模式。
探头的运动状态可通过多种方式获得;本发明具体实施方式中,在所述超声探头中集成运动传感器,以用于实时获取超声探头的运动状态,并根据超声探头在各个时刻下的运动状态判断其是否发生相对变化。
所述运动传感器可为加速度传感器、陀螺仪或磁强计中的至少一种,相应的,运动传感器会随着探头的运动而改变其姿态和/或位置;所述姿态的变化例如:探头运动过程中导致的运动传感器的角度变化;所述位置的变化例如:探头运动过程中导致的运动传感器位置的移动等,在此不做具体赘述。
本发明具体实施方式中,所述运动传感器为三轴加速度传感器,通过所述三轴加速度传感器判断探头运动状态是否发生变化。
具体的,
处理模块300还用于实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的加速度矢量判断当前时间窗口内探头的运动状态;当三轴加速度传感器的加速度矢量满足以下条件至少其中之一时,确认所述探头在当前时间窗口内的运动状态发生变化;
条件1:|Gt|-|Gs|>静态阈值;
条件2:|Gt-Gt-t0|>动态阈值;
其中,Gt表示当前时刻测得的加速度矢量,t表示当前时刻,Gt-t0表示在t-t0时刻测得的加速度矢量,t0的值小于等于时间窗口对应的间隔时间长度阈值;Gs表示当前地理位置标定的静止状态下的加速度矢量;静态阈值、动态阈值均为常数。
该实施方式中,三轴加速度传感器固定设置于探头中,当判断三轴加速传感器的运动状态变化时,即表示探头的运动状态发生同步变化。
该实施方式中,当满足条件1时,三轴加速度传感器处于加速运动状态,当满足条件2时,三轴加速度传感器处于加速运动状态或其姿态迅速发生变化。
如上两个条件满足其中至少一个时,即可认为探头的运动状态发生变化,当然,如果以两个条件同时满足作为判断标准,则获得判断结果会更加精准。
该实施方式中,静态阈值、动态阈值均为一个常数,其可以根据需要具体调整,通常情况下,设备出厂时会为其指定数值;当用户需要时,还可以具体调整其大小,在此不做详细赘述。
当然,在本发明的其他实施方式中,还可以在探头上安装光学或磁场辅助定位装置以替代运动传感器,其同样可以通过运动传感器辅助检测探头运动状态的目的,在此不做详细赘述。
在固定的时间窗口内,处理模块300判断探头的运动状态未发生变化后,可进一步通过回波信号判断是否需要对扫查模式进行调整;而判断探头的运动状态发生变化时,保持当前的扫查参数继续扫查;当然,在本发明其他实施方式中,当探头中未集成运动传感器等零部件时,也可以直接通过回波信号的变化判断探头的工作状态,进而判断是否更改探头的工作模式,以下内容中将会详细描述。
进一步的,获取模块200还用于通过超声探头扫查成像过程中,实时获取每个时刻下的回波信号;
超声成像过程中,超声探头发射声波对待测组织进行逐次的扫查,每次的扫查均有一定的方向性,且没逐次扫查过程中,均可获得多个方向的回波信号,进一步的,将逐次扫查获得的多个方向的回波信息集合到一起,就可以得到一帧扫查图像。
进一步的,处理模块300还用于实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的回波信号,判断当前时间窗口内的回波信号是否发生相对变化;
若是,调整当前成像参数,以使当前的扫查模式调整为动态扫查模式;
若否,调整当前成像参数,以使当前的扫查模式调整为静态扫查模式。
所述动态扫查模式以及静态扫查模式均对应不同的成像参数,如此,通过实时调整成像参数,即可以自适应调整探头的扫查模式,以使获得图像更为精准。
本发明在实现过程中,处理模块300可以采用多种方式判断当前时间窗口内的回波信号是否发生相对变化;本发明具体实施方式中,可在两个层面实现。
扫查线层面,每条扫查线包含一次发射和接收的回波信号的包络,假设图像帧频为10帧/s,则一帧图像中的扫查线最慢可以以1/10=0.1秒的时间间隔重复,通过对比相隔0.1秒先后获得的扫查线包络信息的变化,就可以得出图像信息是否发生改变。处理模块300常用的基于区域的局部匹配准则主要有图像序列中对应像素差的绝对值SAD(sum ofthe absolute differences)算法、图像序列中对应像素差的平方和SSD(sum of squareddifferences)算法、Census变换匹配算法、图像相关性算法。
SAD算法:
SAD算法是一种最简单的匹配算法,此种方法就是以左目图像的源匹配点为中心,定义一个窗口D,其大小为(2m+1)(2n+1),统计其窗口的灰度值的和,然后在右目图像中逐步计算其左右窗口的灰度和的差值,最后搜索到的差值最小的区域的中心像素即为匹配点。
SSD匹配算法:SSD算法与SAD算法相似,在此不做详细赘述。
Census变换匹配算法:
Census变换算法是先对单幅图像的窗口区域进行变换,然后再按照Hamming距离来比较它们的相似性进行匹配。
图像相关性算法,是计算两幅图像匹配区域的互相关性,以当前时刻作为当前时间窗口的结束时刻,获取当前时刻以及其上一时刻的帧图像中每根扫查线所对应的相关系数,并根据其判断当前时间窗口内的回波信号是否发生相对变化。
其中,相关系数以公式
Figure GDA0002532103450000151
表示,该式中,X表示当前时刻帧图像的其中一条扫查线,Y表示上一时刻帧图像中对应X的扫查线r(X,Y)表示对应扫查线的相关系数,Cov(X,Y)表示X与Y的协方差,Var[X]为X的方差,Var[Y]为Y的方差。
该相关系数的值通常介于0~1之间,当其值越接近0时,表示两幅图像相关,当其值越接近1时,表示两幅图像不相关。
进一步的,还可优选的结合当前帧图像中多根扫查线和上帧图像中多帧扫查线对比的情况,以判断回波信息是否发生大幅度的改变,使检测结果更加精准,在此不再继续赘述。
图像层面,扫查线层面的方式同样可以以二维的数学形式应用到二维的图像差异判别上,除此之外类似SIFT(scale invariant feature transform)的特征识别算法也可以被用来捕捉图像中的特征点,通过比对特征点的位置来判断图像是否发生改变。
综上所述,本发明的自动调整参数的超声成像方法及系统,能够适时对成像参数进行调整,简化用户操作,提升在系统在单一工作状态下的适应能力和图像质量。
为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本发明时可以把各模块的功能在同一个或多个软件和/或硬件中实现。
以上所描述的装置实施方式仅仅是示意性的,其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施方式方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施方式中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。

Claims (8)

1.一种自动调整参数的超声成像方法,其特征在于,所述方法包括:
M1、设定时间窗口的间隔时间长度阈值;
M2、通过超声探头扫查成像过程中,实时获取每个时刻下探头的运动状态;
M3、实时以当前时刻作为时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的探头运动状态,判断当前时间窗口内探头的运动状态是否发生相对变化;
若是,保持当前的扫查参数继续扫查;
若否,进入步骤M4;
M4、通过超声探头扫查成像过程中,实时获取每个时刻下的回波信号;
实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的回波信号,判断当前时间窗口内的回波信号是否发生相对变化;
若是,调整当前成像参数,以使当前的扫查模式调整为动态扫查模式;
若否,调整当前成像参数,以使当前的扫查模式调整为静态扫查模式。
2.根据权利要求1所述的自动调整参数的超声成像方法,其特征在于,
所述超声探头中集成运动传感器,以用于实时获取超声探头的运动状态;所述运动传感器为:加速度传感器、陀螺仪或磁强计中的至少一种。
3.根据权利要求2所述的自动调整参数的超声成像方法,其特征在于,所述运动传感器为三轴加速度传感器,
所述步骤M3具体包括:实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的加速度矢量判断当前时间窗口内探头的运动状态;当三轴加速度传感器的加速度矢量满足以下条件至少其中之一时,确认所述探头在当前时间窗口内的运动状态发生变化;
条件1:|Gt|-|Gs|>静态阈值;
条件2:|Gt-Gt-t0|>动态阈值;
其中,Gt表示当前时刻测得的加速度矢量,t表示当前时刻,Gt-t0表示在t-t0时刻测得的加速度矢量,t0的值小于等于时间窗口对应的间隔时间长度阈值;Gs表示当前地理位置标定的静止状态下的加速度矢量;静态阈值、动态阈值均为常数。
4.根据权利要求1所述的自动调整参数的超声成像方法,其特征在于,所述方法还包括:
采用SAD算法、SSD算法、Census变换匹配算法、图像相关性算法、SIFT特征识别算法中的至少一种算法判断当前时间窗口内的回波信号是否发生相对变化;
其中,所述图像相关性算法为:
以当前时刻作为当前时间窗口的结束时刻,获取当前时刻以及其上一时刻的帧图像中每根扫查线所对应的相关系数,并根据其判断当前时间窗口内的回波信号是否发生相对变化;
Figure FDA0002729867860000021
X表示当前时刻帧图像的其中一条扫查线,Y表示上一时刻帧图像中对应X的扫查线, r(X,Y)表示对应扫查线的相关系数,Cov(X,Y)表示X与Y的协方差,Var[X]为X的方差,Var[Y]为Y的方差。
5.一种自动调整参数的超声成像系统,其特征在于,所述系统包括:
计数模块,用于设定时间窗口的间隔时间长度阈值;
获取模块,用于在通过超声探头扫查成像过程中,实时获取每个时刻下的回波信号;
处理模块,用于通过超声探头扫查成像过程中,实时获取每个时刻下探头的运动状态;
实时以当前时刻作为时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的探头运动状态,判断当前时间窗口内探头的运动状态是否发生相对变化;
若发生相对变化,则保持当前的扫查参数继续扫查;
若未发生相对变化,则实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的回波信号,判断当前时间窗口内的回波信号是否发生相对变化;若是,调整当前成像参数,以使当前的扫查模式调整为动态扫查模式;
若否,调整当前成像参数,以使当前的扫查模式调整为静态扫查模式。
6.根据权利要求5所述的自动调整参数的超声成像系统,其特征在于,
所述超声探头中集成运动传感器,以用于实时获取超声探头的运动状态;所述运动传感器为:加速度传感器、陀螺仪或磁强计中的至少一种。
7.根据权利要求6所述的自动调整参数的超声成像系统,其特征在于,
所述运动传感器为三轴加速度传感器,
所述处理模块还用于:实时以当前时刻作为当前时间窗口的结束时刻,根据当前时间窗口内各个时刻分别对应的加速度矢量判断当前时间窗口内探头的运动状态;当三轴加速度传感器的加速度矢量满足以下条件至少其中之一时,确认所述探头在当前时间窗口内的运动状态发生变化;
条件1:|Gt|-|Gs|>静态阈值;
条件2:|Gt-Gt-t0|>动态阈值;
其中,Gt表示当前时刻测得的加速度矢量,t表示当前时刻,Gt-t0表示在t-t0时刻测得的加速度矢量,t0的值小于等于时间窗口对应的间隔时间长度阈值;Gs表示当前地理位置标定的静止状态下的加速度矢量;静态阈值、动态阈值均为常数。
8.根据权利要求5所述的自动调整参数的超声成像系统,其特征在于,所述处理模块采用SAD算法、SSD算法、Census变换匹配算法、图像相关性算法、SIFT特征识别算法中的至少一种算法判断当前时间窗口内的回波信号是否发生相对变化;
其中,所述图像相关性算法为:
以当前时刻作为当前时间窗口的结束时刻,获取当前时刻以及其上一时刻的帧图像中每根扫查线所对应的相关系数,并根据其判断当前时间窗口内的回波信号是否发生相对变化;
Figure FDA0002729867860000041
X表示当前时刻帧图像的其中一条扫查线,Y表示上一时刻帧图像中对应X的扫查线, r(X,Y)表示对应扫查线的相关系数,Cov(X,Y)表示X与Y的协方差,Var[X]为X的方差,Var[Y]为Y的方差。
CN201710972718.0A 2017-10-18 2017-10-18 自动调整参数的超声成像方法及系统 Active CN107811652B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710972718.0A CN107811652B (zh) 2017-10-18 2017-10-18 自动调整参数的超声成像方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710972718.0A CN107811652B (zh) 2017-10-18 2017-10-18 自动调整参数的超声成像方法及系统

Publications (2)

Publication Number Publication Date
CN107811652A CN107811652A (zh) 2018-03-20
CN107811652B true CN107811652B (zh) 2021-01-26

Family

ID=61608103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710972718.0A Active CN107811652B (zh) 2017-10-18 2017-10-18 自动调整参数的超声成像方法及系统

Country Status (1)

Country Link
CN (1) CN107811652B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111685793A (zh) * 2019-03-15 2020-09-22 通用电气公司 用于对成像系统参数进行基于图像的控制的装置和方法
CN110613479A (zh) * 2019-10-22 2019-12-27 无锡祥生医疗科技股份有限公司 超声设备降功耗的方法和装置
CN114073547B (zh) * 2020-08-14 2024-06-07 深圳迈瑞生物医疗电子股份有限公司 用于超声成像的方法、装置和存储介质
CN113536964B (zh) * 2021-06-25 2023-09-26 合肥合滨智能机器人有限公司 一种超声视频的分类提取方法
CN113576524A (zh) * 2021-08-02 2021-11-02 深圳技术大学 一种成像模式自动切换方法、装置及超声成像系统

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9078592B2 (en) * 2008-10-27 2015-07-14 Wisconsin Alumni Research Foundation Ultrasonic strain imaging device with selectable cost-function
CN103181778B (zh) * 2011-12-31 2016-05-25 深圳迈瑞生物医疗电子股份有限公司 超声成像方法及系统
JP6103857B2 (ja) * 2012-08-28 2017-03-29 キヤノン株式会社 被検体情報取得装置、表示方法、及びプログラム
CN104739452B (zh) * 2013-12-30 2019-02-12 深圳迈瑞生物医疗电子股份有限公司 一种超声成像装置及方法
US9918701B2 (en) * 2014-09-03 2018-03-20 Contextvision Ab Methods and systems for automatic control of subjective image quality in imaging of objects
CN104811619B (zh) * 2015-04-23 2018-12-25 小米科技有限责任公司 成像校正方法及装置、电子设备
US20160310110A1 (en) * 2015-04-23 2016-10-27 Siemens Medical Solutions Usa, Inc. Acquisition control for mixed mode ultrasound imaging
CN106361375B (zh) * 2016-09-14 2019-03-19 飞依诺科技(苏州)有限公司 用于超声脉冲多普勒成像的孔径自动调节方法及系统

Also Published As

Publication number Publication date
CN107811652A (zh) 2018-03-20

Similar Documents

Publication Publication Date Title
CN107811652B (zh) 自动调整参数的超声成像方法及系统
CN111432733B (zh) 用于确定超声探头的运动的设备和方法
CN109044398B (zh) 超声系统成像方法、装置及计算机可读存储介质
US5795296A (en) Pipeline process for automatically measuring object boundary from ultrasound image samples
US20150023561A1 (en) Dynamic ultrasound processing using object motion calculation
CN110595602B (zh) 振动检测方法及相关产品
US8867813B2 (en) Ultrasonic imaging device, ultrasonic imaging method and program for ultrasonic imaging
CN111602139A (zh) 图像处理方法、装置、控制终端及可移动设备
KR20150120774A (ko) 관심영역 검출 시스템 및 방법
US20100138191A1 (en) Method and system for acquiring and transforming ultrasound data
WO2016015994A1 (en) Ultrasound imaging apparatus
EP1513104A2 (en) Improvements in or relating to similarity measures
US20160213353A1 (en) Ultrasound imaging apparatus, ultrasound imaging method and ultrasound imaging program
CN107405134B (zh) 超声成像装置
CN109452954B (zh) 超声成像方法及装置
CN113168696A (zh) 参数选定装置、参数选定方法以及参数选定程序
US20210219961A1 (en) Compounding and non-rigid image registration for ultrasound speckle reduction
US20200301004A1 (en) Compensating radio tracking with comparison to image based tracking
CN115204221A (zh) 生理参数的检测方法、设备及存储介质
US11430120B2 (en) Ultrasound image evaluation apparatus, ultrasound image evaluation method, and computer-readable non-transitory recording medium storing ultrasound image evaluation program
US11118915B2 (en) Position estimation device, moving-object control system, position estimation method, and computer program product
JPH0926473A (ja) 時空間微分法を用いた計測装置
CN112237439A (zh) 调整探头位置的方法、探头、超声设备和存储介质
CN114767074B (zh) 一种生命体征测量方法、设备及存储介质
US20210038184A1 (en) Ultrasound diagnostic device and ultrasound image processing method

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
CP01 Change in the name or title of a patent holder

Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee after: Feiyinuo Technology Co.,Ltd.

Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee before: Feiyinuo Technology (Suzhou) Co.,Ltd.

Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee after: Feiyinuo Technology (Suzhou) Co.,Ltd.

Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee before: VINNO TECHNOLOGY (SUZHOU) Co.,Ltd.

CP01 Change in the name or title of a patent holder