本发明涉及肌肉测量技术领域，特别是涉及一种实时检测肌肉形态参数的超声成像设备及方法。
背景技术
肌肉是构成人体的重要组织，它密布在各组织器官和骨骼的周围，人体总共有600多块肌肉，按结构和功能肌肉可以分为骨骼肌、心肌和平滑肌三种。肌肉通过收缩和舒张可以帮助我们完成食物在消化道的顺利运输、心脏的正常搏动以及和骨骼配合让我们完成各种复杂的运动。肌肉在人体的生命活动中扮演着重要的作用，然而肌肉的构成又非常的复杂，定量分析和评估肌肉的各项功能是当前临床医学研究的难点和热点。
测量肌肉运动的最主要方法是表面肌电信号(surfaceelectromyography，semg)，它在康复医学、人体动力学、生物肌肉力学等方面都有广泛的应用。肌电信号的时频特性直接反映肌肉力的大小、肌肉的功能和状态、肌群的相互协作协调等特性。然而，肌电信号容易受到各种潜在因素的影响，如电极位置、肌肉类型、邻近肌肉干扰等，这些都制约了肌电信号在肌肉评估中的应用。
超声成像(ultrasound,us)相对于现在比较流行的x线成像、计算机断层成像(computertomography,ct)以及磁共振成像(magneticresonanceimaging,mri)有着实时、快速、无辐射、价格低廉等优势；同时超声成像可获得肌肉收缩和舒张时的实时动态图像，已被广泛应用于研究肌肉的形态参数，这些形态参数和肌肉组织的力学特性和生理学功能直接相关，从而超声成像为临床观察和诊断肌肉组织特性提供了便捷途径。
现有超声成像技术虽然可以实时检测肌肉在收缩过程中的形态结构变化，但这种检测是定性的，只能凭医生的经验来观察。现有的超声设备，都没有动态测量肌肉各项形态参数的功能。现有的对超声肌肉数据的研究全部对离线数据的分析。但这种方法中间过程很复杂，需要配合其他一些设备，例如图像采集卡，无法做到实时测量，导致测量结果滞后。
技术实现要素：
本发明的目的是提供一种实时检测肌肉形态参数的超声成像设备及方法，可实现实时测量肌肉动态变化。
为实现上述目的，本发明提供了如下方案：
一种实时检测肌肉形态参数的超声成像设备，所述超声成像设备包括：
超声成像模块，用于获取肌肉运动过程中的超声图像序列；
图像序列处理模块，用于根据所述超声图像序列中的每一帧图像，获得肌肉形态参数随时间变化的动态曲线。
可选的，所述超声成像设备还包括：
显示模块，分别与超声成像模块和/或所述图像序列处理模块连接，用于显示肌肉超声图像，和/或肌肉形态参数随时间变化的动态曲线。
可选的，所述超声成像模块包括：
探头，用于向待检测区域发射超声波信号，并接收来自所述待检测区域的回波信号；
发射电路，用于激励所述探头向扫描目标区域发射超声波；
收发开关，用于将所述发射电路产生的激励信号发送至所述超声探头，并将所述超声回波信号发送至所述接收电路；
接收电路，用于接收由目标区域反射的超声回波信号，并形成数字信号；
波束合成电路，用于对所述数字信号进行动态聚焦、动态加权和求和处理得到超声扫描线；
信号处理单元，用于对超声扫描线进行动态滤波处理，得到滤波后的超声扫描线；
图像处理单元，用于对所述经过信号处理单元的超声扫描线进行包络检测、对数压缩和数字扫描变换，得到所述目标区域的超声图像；
图像存储单元，用于与所述图像处理单元和图像预处理单元连接，用于存储所述超声图像。
可选的，所述图像序列处理模块还包括：
图像预处理单元，与所述图像存储单元和参数提取单元连接，用于对所述图像序列进行数据降噪、滤波处理，得到预处理图像信息；
参数提取单元，与所述图像预处理单元和显示模块连接，用于提取当前超声图像中的肌肉形态参数，得到肌肉参数随时间变化的动态曲线，并发送至所述显示模块进行显示。
可选的，所述图像处理方法还包括：
利用hough变换和radon变换算法，得到肌肉腱膜位置和肌束方向的相关信息，以及利用三角函数公式得到肌肉各参数的值。
可选的，所述肌肉形态参数包括肌束的羽状角、肌束长度、肌肉厚度及肌肉横截面积中至少一者。
为实现上述目的，本发明提供了如下方案：
一种实时检测肌肉形态参数的超声成像方法，所述超声成像方法包括：
获取肌肉运动过程中的超声图像序列；
根据所述超声图像序列中的每一帧图像，获得肌肉形态参数随时间变化的动态曲线。
可选的，所述超声成像方法还包括：
显示肌肉超声图像，和/或肌肉形态参数随时间变化的动态曲线。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明通过超声成像模块和图像序列处理模块配合，实时获取肌肉运动过程中的超声图像序列，并根据所述超声图像序列中的每一帧图像，获得肌肉形态参数随时间变化的动态曲线，从而实现对肌肉动态变化的实时检测。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例实时检测肌肉形态参数的超声成像设备的模块结构示意图；
图2为本发明实施例实时检测肌肉形态参数的超声成像设备的实施例结构示意图；
图3为肌肉形态结构示意图；
图4为本发明实施例实时检测肌肉形态参数的超声成像方法的流程图；
图5为肌肉形态参数提取的流程图。
符号说明：
超声成像模块—1，发射电路—11，收发开关—12，超声波探头—13，接收电路—14，波束合成单元—15，信号处理单元—16，图像处理单元—17，图像存储单元—18，图像序列处理模块—2，图像预处理单元—21，参数提取单元—22，显示模块—3。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种实时检测肌肉形态参数的超声成像设备，通过超声成像模块和图像序列处理模块配合，实时获取肌肉运动过程中的超声图像序列，并根据所述超声图像序列中的每一帧图像，获得肌肉形态参数随时间变化的动态曲线，从而实现对肌肉动态变化的实时检测。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示，本发明实时检测肌肉形态参数的超声成像设备包括超声成像模块1、图像序列处理模块2及显示模块3。
其中，所述超声成像模块1用于获取肌肉运动过程中的超声图像序列；所述图像序列处理模块2用于根据所述超声图像序列中的每一帧图像，获得肌肉形态参数随时间变化的动态曲线。
进一步地，本发明超声成像设备还包括显示模块3，所述显示模块3分别与超声成像模块1和/或所述图像序列处理模块2连接，所述显示模块3用于显示肌肉超声图像，和/或肌肉形态参数随时间变化的动态曲线。
如图2所示，所述超声成像模块1包括发射电路11、收发开关12、超声波探头13、接收电路14、波束合成单元15、信号处理单元16、图像处理单元17及图像存储单元18。
其中，所述发射电路11用于触发有序脉冲，产生激励信号。所述收发开关12分别与超声波探头13、发射电路11及接收电路14连接，用于将所述激励信号发送至所述超声波探头13，以控制所述超声波探头13向待检测区域发射超声波信号，并将超声回波信号发送至所述接收电路14。所述接收电路14与波束合成单元15连接，用于接收超声波探头13反射回来的超声回波信号，并形成数字信号。所述波束合成单元15与信号处理单元16连接，用于对所述数字信号进行动态聚焦、动态加权和求和处理得到超声扫描线。所述信号处理单元16与图像处理单元17连接，用于对所述超声扫描线进行进一步处理。所述图像处理单元17与信号处理单元16和图像存储单元18连接，用于对所述经过信号处理单元的超声扫描线进行包络检测、对数压缩和扫描变换，得到所述目标区域的超声图像。所述图像存储单元18与图像处理单元17、图像预处理单元21和显示模块连接，用于存储所述超声图像。
其中，所述超声波探头13是超声信号发射和接收的主要部件，所述超声波探头13的性能直接影响超声波的特性从而影响成像效果。在本实施例中，成像部位是肌肉在人体的浅表位置，为了后续图像处理和观察方便，所述超声波探头为频率为5.0-12mhz的线阵探头，这样既能满足对测量深度的要求也能满足对分辨率的要求。
具体地，所述超声波探头13包括电极、压电晶体、吸声块、声学绝缘层、匹配保护层和信号线。
此外，所述超声波探头13还包括壳体，所述壳体内容易有所述电极、压电晶体、吸声块、声学绝缘层、匹配保护层和信号线。
具体地，所述波束合成单元包括：动态聚焦模块、动态加权模块及求和模块。
所述动态聚焦模块与接收电路14连接，所述动态聚焦模块用于根据数字信号的深度，由浅到深地改变焦距，得到聚焦后的超声回波信号；
所述动态加权模块与所述动态聚焦模块连接，所述动态加权模块用于将延时之后的数据乘上不同的加权系数，然后累加，得到加权后的信号；
所述求和模块与所述动态加权模块连接，所述求和模块用于将各个通道加权后的数据进行求和处理，得到超声扫描线。
具体地，所述信号处理单元16还包括动态滤波模块。
所述动态滤波模块与所述波束合成单元15连接，所述动态滤波模块用于动态捕捉携带有用信息的回波信号，滤除皮肤表面回波、近场低频回波和远场高频干扰，得到高信噪比的超声回波信号。
具体地，所述图像处理单元17还包括：包络检测模块、对数压缩模块及数字扫描变换模块。
所述包络检测模块与所述动态滤波模块连接，所述包络检测模块用于把有用的人体组织结构信息从载波信号里解调出来，得到人体组织结构的包络信息。
所述对数压缩模块与所述包络检测模块连接，所述对数压缩模块用于压缩超声回波信号的动态范围，得到动态范围20db左右的超声回波信号。
所述扫数字描变换模块与所述对数压缩模块连接，所述扫数字描变换模块用于将超声扫描线数据变成可以在屏幕上显示的图像数据。
如图2所示，所述图像序列处理模块2包括图像预处理单元21、及参数提取单元22。
所述图像预处理单元21，与所述图像存储单元18及参数提取单元22连接，用于对所述图像序列进行数据降噪、滤波处理，得到预处理图像信息。
所述参数提取单元22，与所述图像预处理单元21及显示模块3连接，用于提取当前超声图像中的肌肉形态参数，得到肌肉参数随时间变化的动态曲线，并发送至所述显示模块3进行显示。
其中，如图3所示，所述肌肉形态参数包括肌束的羽状角(α)，肌束长度(fasciclelength,fl)，肌肉厚度(musclethickness,mt)，肌肉横截面积(psca)中至少一者。图5所示为所述参数提取单元22提取当前超声图像中的肌肉形态参数的具体步骤，其中，迭代hough变换步骤主要用于找出肌肉超声图像中的深浅腱膜，确定深浅腱膜位置后，计算出深浅腱膜距离的平均值，即为肌肉厚度t。肌肉横截面积p近似于以t为直径的圆的面积，即
频域radon变换步骤主要是先对肌肉超声图像进行二维傅立叶变换，然后进行radon变换，从而确定肌束的方向。深层腱膜和肌束方向形成的夹角即为羽状角θ。根据图3可以看出肌肉厚度、肌束长度和羽状角分别为直角三角形的一条直角边、斜边和一个锐角。肌肉长度即为
此外，本发明还提供一种实时检测肌肉形态参数的超声成像方法。如图4所示，本发明实时检测肌肉形态参数的超声成像方法包括：
步骤100：获取肌肉运动过程中的超声图像序列；
步骤200：根据所述超声图像序列中的每一帧图像，获得肌肉形态参数随时间变化的动态曲线；
步骤300：显示肌肉超声图像，和/或肌肉形态参数随时间变化的动态曲线。
相对于现有技术，本发明实时检测肌肉形态参数的超声成像方法与上述实时检测肌肉形态参数的超声成像设备的方法相同，在此不再赘述。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。