CN107961038A - 一种根据超声弹性肌动图获取生物力学参数的方法及装置 - Google Patents

一种根据超声弹性肌动图获取生物力学参数的方法及装置 Download PDF

Info

Publication number
CN107961038A
CN107961038A CN201711322506.4A CN201711322506A CN107961038A CN 107961038 A CN107961038 A CN 107961038A CN 201711322506 A CN201711322506 A CN 201711322506A CN 107961038 A CN107961038 A CN 107961038A
Authority
CN
China
Prior art keywords
muscle
skeletal muscle
muscular
parameter
dynamic
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
Application number
CN201711322506.4A
Other languages
English (en)
Other versions
CN107961038B (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.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
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 Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201711322506.4A priority Critical patent/CN107961038B/zh
Priority to US16/482,802 priority patent/US11452502B2/en
Priority to EP17934679.6A priority patent/EP3566655B1/en
Priority to PCT/CN2017/118725 priority patent/WO2019114034A1/zh
Publication of CN107961038A publication Critical patent/CN107961038A/zh
Application granted granted Critical
Publication of CN107961038B publication Critical patent/CN107961038B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本申请实施例提供了一种根据超声弹性肌动图获取生物力学参数的方法及装置，其中，该方法包括：同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列；分别获取所述动态肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数，以及，所述动态弹性图像序列中各弹性图像对应的弹性模量值；以所述肌肉动力学参数为横坐标，以同步的所述弹性模量值为纵坐标生成超声弹性肌动图曲线，根据所述超声弹性肌动图曲线估算肌肉生物力学参数。本申请实施例能够获得动态变化的生物力学参数，且获得的肌肉生物力学参数精确度较高。

Description

一种根据超声弹性肌动图获取生物力学参数的方法及装置
技术领域
本申请涉及测试技术领域，具体而言，涉及一种根据超声弹性肌动图获取生物力学参数的方法及装置。
背景技术
肌肉是构成人体的重要组织，占体重的40％～45％，分为骨骼肌、平滑肌和心肌三种类型，分布在各器官和骨骼的周围，具有不同的构造、力学性质和功能。骨肌系统依靠肌肉的约束力保持稳定，依靠肌肉的驱动力形成行为运动，在人体运动中扮演者至关重要的作用。
目前，对于肌肉生物力学参数的测试主要是基于超声剪切波弹性成像技术来完成的，超声剪切波弹性成像技术实现了对肌肉组织弹性真正的定量测量，其基本原理是：由超声换能器产生的辐射力在聚焦位置附近产生剪切波，然后测量剪切波在组织中的传播速度，计算出组织的弹性模量。
但是，现有技术中，主要是基于超声剪切波弹性成像技术按照以下方式测量肌肉生物力学参数的：(1)测量松弛状态下骨骼肌的弹性模量值，通过单一的弹性模量值差异评估，肌肉是人体完成各种动作的动力来源，运动过程中肌肉力会引起弹性模量变化，肌肉松弛状态下的弹性模量值智能反映肌力为零时的硬度，无法充分表征其生物力学特性。(2)分别测量松弛状态和紧张状态(主动收缩)下骨骼肌的弹性模量值，用以评估肌肉的功能、张力情况，但是主动收缩产生肌力难以做出定量的评估。(3)测量不同关节角度下骨骼肌的弹性模量值，研究拉伸状态下肌肉租住的生物力学特性，但是不同的关节角度仍然是定性的指标参数，因人而异，不能对肌力做出定量的评估。
综上，现有的基于超声剪切波弹性成像技术测试肌肉生物力学参数的方法无法获取动态信息，且无法定量获取，所以获得的肌肉生物力学参数精确度不高。
发明内容
有鉴于此，本申请的目的在于提供一种肌肉生物力学参数的测试方法及装置，获取动态变化的生物力学参数及提高获取肌肉生物力学参数的精确度。
第一方面，本申请实施提供了一种根据超声弹性肌动图获取生物力学参数的方法，包括：
同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列；
分别获取所述动态肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数，以及，所述动态弹性图像序列中各弹性图像对应的弹性模量值；
以所述肌肉动力学参数为横坐标，以同步的所述弹性模量值为纵坐标生成超声弹性肌动图曲线，根据所述超声弹性肌动图曲线估算肌肉生物力学参数。
结合第一方面，本申请实施例提供了第一方面的第一种可能的实施方式，所述肌肉动力学图像包括位移图像，所述肌肉动力学参数包括肌肉拉伸长度，在所述同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列之前，所述方法还包括：
确定初始状态下所述单块骨骼肌的肌肉组织伸长端的特征点；
计算所述特征点至所述单块骨骼肌的肌肉组织伸长端的对端的距离，得到所述骨骼肌的初始长度。
结合第一方面的第一种可能的实施方式，本申请实施提供了第一方面的第二种可能的实施方式，确定所述单块骨骼肌的肌肉组织伸长端的特征点包括：
控制超声诊断仪通过成像模式扫描所述单块骨骼肌沿肌束方向的纵切面，将检测到的所述单块骨骼肌与肌腱的结合点作为所述特征点。
结合第一方面的第一种可能的实施方式或第一方面的第二种可能的实施方式，本申请实施提供了第一方面的第三种可能的实施方式，所述获取所述动态肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数包括：
获取所述特征点在所述位移图像中相对初始状态下的特征点位置的位移变化的有效值；
将所述初始长度和所述位移变化的有效值之和记为所述位移图像对应的肌肉拉伸长度。
结合第一方面的第三种可能的实施方式，本申请实施提供了第一方面的第四种可能的实施方式，所述获取所述特征点在所述位移图像中相对初始状态下的特征点位置的位移变化的有效值包括：
控制肌电仪监控所述单块骨骼肌及周围肌群的肌电信号；
当所述肌电信号的幅值小于设定值时，将所述位移图像中相对初始状态下的特征点位置的位移变化作为所述有效值。
结合第一方面，本申请实施提供了第一方面的第五种可能的实施方式，在所述同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列之前，所述方法还包括：
确定所述单块骨骼肌的肌腹位置，以在所述肌腹位置处采集所述动态弹性图像序列。
结合第一方面的第五种可能的实施方式，本申请实施提供了第一方面的第六种可能的实施方式，所述确定所述单块骨骼肌的肌腹位置，包括：
控制超声诊断仪通过成像模式扫描所述骨骼肌中垂直于肌束方向的横切面，将横切面最大的位置作为所述单块骨骼肌的肌腹位置。
结合第一方面的第五种可能的实施方式或第一方面的第六种可能的实施方式，本申请实施提供了第一方面的第七种可能的实施方式，所述获取动态弹性图像序列中各动态弹性图像对应的弹性模量值，包括：
在所述单块骨骼肌的肌腹位置沿肌束方向输出所述骨骼肌在连续拉伸时的所述动态弹性图像序列；
控制所述超声诊断仪通过分析模式对所述动态弹性图像序列中的每一幅弹性图像中的设定区域进行定量分析后，得到各弹性图像对应的弹性模量值。
结合第一方方面至第一方面的第二种可能的实施方式、第一方面的第五种可能的实施方式或第一方面的第六种可能的实施方式中的任一项，本申请实施例提供了第一方面的第八种可能的实施方式，在所述以所述肌肉动力学参数为横坐标，以同步的所述弹性模量值为纵坐标生成超声弹性肌动图曲线之后，根据所述超声弹性肌动图曲线估算肌肉生物力学参数之前，所述方法还包括：
确定所述超声弹性肌动图曲线与设定的力学模型曲线是否吻合，若吻合，执行所述根据所述超声弹性肌动图曲线估算肌肉生物力学参数的步骤。
结合第一方面的八种可能的实施方式，本申请实施例提供了第一方面的第九种可能的实施方式，所述肌肉动力学图像包括位移图像，所述肌肉动力学参数包括肌肉拉伸长度，所述根据所述超声弹性肌动图曲线估算肌肉生物力学参数包括：
获取所述超声弹性肌动图曲线对应的模型公式，以及所述模型公式中的各个参数；
根据所述参数确定所述单块骨骼肌在松弛状态下的长度、与所述单块骨骼肌在松弛状态下的长度对应的弹性模量值以及所述单块骨骼肌的延展度，所述延展度为所述单块骨骼肌的最大拉伸长度与所述单块骨骼肌在松弛状态下的长度之差；
在所述骨骼肌长度大于所述单块骨骼肌在松弛状态下的长度时，对所述模型公式按照高斯牛顿迭代法，确定所述单块骨骼肌的硬度指标。
第二方面，本申请实施例提供了一种根据超声弹性肌动图获取肌肉生物力学参数的装置，包括：
采集模块，用于同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列；
获取模块，用于分别获取所述动态肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数，以及，所述动态弹性图像序列中各弹性图像对应的弹性模量值；
确定模块，用于以所述肌肉动力学参数为横坐标，以同步的所述弹性模量值为纵坐标生成超声弹性肌动图曲线，根据所述超声弹性肌动图曲线估算肌肉生物力学参数。
与现有技术相比，本申请实施例是通过同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列，分别获取肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数，以及，动态弹性图像序列中各动态弹性图像对应的弹性模量值；以所述肌肉动力学参数为横坐标，以同步的所述弹性模量值为纵坐标生成超声弹性肌动图曲线，根据所述超声弹性肌动图曲线估算肌肉生物力学参数，可见，本申请获得的肌肉动力学参数和弹性模量值是在单块骨骼肌连续拉伸时测得的变化参数，且最终生超声弹性肌动图曲线，可以由该曲线估算单块骨骼肌在连续拉伸过程中的生物力学参数，这样得到的生物力学参数是动态的且是定量的。
为使本申请的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本申请的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1示出了本申请实施例所提供的一种根据超声弹性肌动图获取肌肉生物力学参数的方法流程图；
图2示出了本申请实施例所提供的一种获取骨骼肌初始长度的方法流程图；
图3示出了本申请实施例所提供的一种位移图像序列图；
图4示出了本申请实施例所提供的一种弹性模量序列图；
图5示出了本申请实施例所提供的一种获取骨骼肌肌肉拉伸长度方法流程图；
图6示出了本申请实施例所提供的得到各动态弹性图像对应的弹性模量值流程图；
图7示出了本申请实施例所提供的估算肌肉生物力学参数的方法流程图；
图8示出了本申请实施例所提供的弹性肌动图曲线示意图；
图9示出了本申请实施例所提供一种根据超声弹性肌动图获取肌肉生物力学参数的装置结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围，而是仅仅表示本申请的选定实施例。基于本申请的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。
实施例1
本申请实施例1提供了一种根据超声弹性肌动图获取肌肉生物力学参数的方法，其方法流程图如图1所述，具体步骤如下：
S100，同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列。
这里的同步采集，是指每一幅肌肉动力学图像对对应一幅弹性图像，即使得之后获取的弹性模量值随着肌肉动力学参数的变化而变化。
S110，分别获取动态肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数，以及，动态弹性图像序列中各弹性图像对应的弹性模量值。
S120，以肌肉动力学参数为横坐标，以同步的弹性模量值为纵坐标生成超声弹性肌动图曲线，根据超声弹性肌动图曲线估算肌肉生物力学参数。
其中，肌肉动力学参数可以包括肌肉拉伸长度、肌力、应变和关节角度等，相应的肌肉动力学图像序列可以包括位移图像序列、肌力图像序列、应变图像序列和关节角度图像序列，本申请实施例中以位移图像序列为例进行详细说明。
一种较佳的实施方式，在本申请实施例1提出的技术方案中，肌肉动力学图像包括位移图像，肌肉动力学参数包括肌肉拉伸长度，在S100步骤中，即在同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列之前，该方法还包括如下步骤，流程图如图2所示：
S200，确定初始状态下单块骨骼肌的肌肉组织伸长端的特征点。
要求被测者在测试前2小时内禁止进行剧烈活动，确保能够采集到被测者在初始状态下的单块骨骼肌的肌肉组织伸长端的特征点。
可选地，确定单块骨骼肌的肌肉组织伸长端的特征点包括：
控制超声诊断仪通过成像模式扫描单块骨骼肌沿肌束方向的纵切面，将检测到的单块骨骼肌与肌腱的结合点作为特征点。
S210，计算特征点至单块骨骼肌的肌肉组织伸长端的对端的距离，得到骨骼肌的初始长度。
这里的初始长度即单块骨骼肌在拉伸前的长度，即特征点位于该单块骨骼肌两端的一个端点，初始长度为该单块骨骼肌的一个端点到对端的长度。
如图3所示，即为实验获取的位移图像序列，以及，图4为实验同步获取的弹性图像序列。
一种较佳的实施方式，在本申请实施例1提出的技术方案中，在步骤S110中，即获取动态肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数包括以下步骤，流程图如图5所示：
S300，获取特征点在位移图像中相对初始状态下的特征点位置的位移变化的有效值。
其中，获取特征点在位移图像中相对初始状态下的特征点位置的位移变化的有效值包括：
控制肌电仪监控单块骨骼肌及周围肌群的肌电信号。
当肌电信号的幅值小于设定值时，将位移图像中相对初始状态下的特征点位置的位移变化作为有效值。
为了避免被测者肌肉无意识主动收缩产生伪像，测试时采集被测单块骨骼肌及周围肌群的肌电信号进行监控，当肌电信号幅度小于最大收缩时幅度的1％时，认为采集的位移变化作为有效值进行记录。
S310，将初始长度和位移变化的有效值之和记为位移图像对应的肌肉拉伸长度。
若初始长度记为linit,位移变化的有效值记为Δl，则肌肉拉伸长度l按照以下公式(1)计算：
l＝linit+Δl (1)
一种较佳的实施方式，在本申请实施例1提出的技术方案中，在步骤S100中，即在同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列之前，该方法还包括：
确定单块骨骼肌的肌腹位置，以在肌腹位置处采集所述动态弹性图像序列。
按照以下方式确定所述单块骨骼肌的肌腹位置：
控制超声诊断仪通过成像模式扫描所述骨骼肌中垂直于肌束方向的横切面，将横切面最大的位置作为所述单块骨骼肌的肌腹位置。
本申请实施例中控制超声诊断仪的成像模式扫查垂直与待测骨骼肌肌束的横截面，在最大横截面对应的位置旋转超声诊断仪的探头，沿肌束方向纵切面检查骨骼肌长轴，启动弹性成像模式来获取动态弹性图像序列中各弹性图像对应的弹性模量值。
其中，获取弹性图像序列中各弹性图像对应的弹性模量值，包括以下步骤，如图6所示为其流程示意图：
S400，在单块骨骼肌的肌腹位置沿肌束方向输出骨骼肌在连续拉伸时的动态弹性图像序列。
S410，控制超声诊断仪通过分析模式对动态弹性图像序列中的每一幅弹性图像中的设定区域进行定量分析后，得到各动态弹性图像对应的弹性模量值。
控制超生诊断仪在找到骨骼肌的肌腹位置后，启动弹性成像模式，固定探头位置，按照设定频率输出被测骨骼肌在被连续拉伸时的动态弹性图像序列，如图4所示，在确立每一幅弹性图像的感兴趣区即测试区，比如10mm×10mm的区域后,启动超声诊断仪的分析模式对该区域进行定量分析后，输出与各弹性图像对应的弹性模量值。
较佳地，在测量待测单块骨骼肌的肌肉拉伸长度以及对该肌肉拉伸长度对应的弹性模量值时，为了降低测量误差，可以采取测量被测者多次，取多次的平均值作为最终的测量结果，比如测量3次。
一种较佳的实施方式，在本申请实施例1提出的技术方案中，在步骤S120中，即在以肌肉动力学参数为横坐标，以同步的弹性模量值为纵坐标生成超声弹性肌动图曲线之后，根据超声弹性肌动图曲线估算肌肉生物力学参数之前，该方法还包括：
确定超声弹性肌动图曲线与设定的力学模型曲线是否吻合，若吻合，执行根据超声弹性肌动图曲线估算肌肉生物力学参数的步骤。
在以肌肉拉伸长度为横坐标，以弹性模量值以纵坐标生成超声弹性肌动图曲线后，根据该曲线的形状，判断是否与设定的力学模型曲线吻合，在吻合时，执行后续的根据超声弹性肌动图曲线估算肌肉生物力学参数的步骤。
一种较佳的实施方式，在本申请实施例1提出的技术方案中，肌肉动力学图像包括位移图像，肌肉动力学图像包括位移图像，在步骤S120中，根据超声弹性肌动图曲线估算肌肉生物力学参数包括以下步骤，流程图如图7所示：
S500，获取超声弹性肌动图曲线对应的模型公式，以及模型公式中的各个参数。
S510，根据参数确定单块骨骼肌在松弛状态下的长度、与单块骨骼肌在松弛状态下的长度对应的弹性模量值以及单块骨骼肌的延展度，延展度为单块骨骼肌的最大拉伸长度与所述单块骨骼肌在松弛状态下的长度之差。
S520，在骨骼肌长度大于单块骨骼肌在松弛状态下的长度时，对模型公式按照高斯牛顿迭代法，确定单块骨骼肌的硬度指标。
单块骨骼肌的松弛状态是指该单块骨骼肌从初始状态开始被拉伸过程中，肌力为零的一个状态，在该松弛状态下，该单块骨骼肌的弹性模量几乎不变。
在以肌肉拉伸长度为横坐标，以弹性模量值以纵坐标生成超声弹性肌动图曲线如图8所示，得到的该超声弹性肌动图曲线对应的模型公式如下式(2)所示：
上式中，E0为肌肉组织松弛弹性模量值；l0为肌肉松弛长度；α为硬度指标；k为比例系数；弹性肌动图曲线中纵坐标对应的弹性模量值在肌肉处于松弛状态时一直是平坦不变的，其显著增大的变化结点，即开始产生肌力对应的横坐标长度作为l0的估计值，如图8所示；弹性肌动图中横坐标小于l0的区间对应的所有测量弹性模量的均值Emean作为E0的估计值，如图8所示；延展度le定义为测量的最大拉伸长度lmax-l0；α需要从横坐标大于等于l0序列区间(l≥l0)中估计出来，具体步骤如下：
令
其中需要估计的未知参数为a＝[α,k]，使用最小二乘判别原则，采用基于泰勒级数的迭代法，记初始的试验系数为a0，即：
y＝f(x,a0)； (4)
初始参数通常不满足误差的平方和最小，通过调整系数a0使误差的平方和最小，记改进的系数为a1，则有：
a1＝a0+△a； (5)
将公式(5)代入公式(4)得到：
y＝f(x,a1)＝f(x,a0+△a)； (6)
将函数按泰勒级数展开，保留一阶导数项，得到：
令与yi之间的误差为：
误差的平方和为：
需要求误差平方和的极小值，则：
即：
可用矩阵记号表示为：K(△a)＝b(12)
其中△a的元素为Δap(p＝0,1)，且:
上述方程组采用高斯消元法求解△a，继而算出a1的值，最后，将算得的a1作为初值，重复上述步骤，反复迭代和修正△a的值，直至△a的范数小于一个指定的阈值。
由图8所示，该图为正常人腓肠肌的弹性肌动图，以单块骨骼肌的伸展长度为横坐标，以对应的弹性模量值为纵坐标，从图中估计出肌肉松弛弹性模量E0＝8.5kPa；松弛长度l0＝0.427m；延展度le＝0.026m；硬度指标α＝158.8。
实施例2
本申请实施例2提供了一种根据超声弹性肌动图获取肌肉生物力学参数的装置，如图9所示，包括：采集模块901、获取模块902和确定模块903。
采集模块901，用于同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列。
获取模块902，用于分别获取动态肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数，以及，动态弹性图像序列中各弹性图像对应的弹性模量值。
确定模块903，用于以所述肌肉动力学参数为横坐标，以同步的弹性模量值为纵坐标生成超声弹性肌动图曲线，根据超声弹性肌动图曲线估算肌肉生物力学参数。
较佳地，肌肉动力学图像包括位移图像，肌肉动力学参数包括肌肉拉伸长度，采集模块901在同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列之前，还用于：
确定初始状态下单块骨骼肌的肌肉组织伸长端的特征点。
计算特征点至单块骨骼肌的肌肉组织伸长端的对端的距离，得到骨骼肌的初始长度。
较佳地，采集模块901具体用于：
控制超声诊断仪通过成像模式扫描所述单块骨骼肌沿肌束方向的纵切面，将检测到的单块骨骼肌与肌腱的结合点作为特征点。
较佳地，获取模块902具体用于：
获取特征点在位移图像中相对初始状态下的特征点位置的位移变化的有效值。
将初始长度和位移变化的有效值之和记为位移图像对应的肌肉拉伸长度。
较佳地，获取模块902按照以下方式确定有效值：
控制肌电仪监控单块骨骼肌及周围肌群的肌电信号。
当肌电信号的幅值小于设定值时，将位移图像中相对初始状态下的特征点位置的位移变化作为有效值。
较佳地，采集模块901在同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列之前，还用于：
确定单块骨骼肌的肌腹位置，以在肌腹位置处采集所述动态弹性图像序列。
较佳地，采集模块901具体根据以下方式确定单块骨骼肌的肌腹位置：
控制超声诊断仪通过成像模式扫描所述骨骼肌中垂直于肌束方向的横切面，将横切面最大的位置作为单块骨骼肌的肌腹位置。
较佳地，获取模块902具体按照以下方式获取弹性模量值：
在单块骨骼肌的肌腹位置沿肌束方向输出所述骨骼肌在连续拉伸时的动态弹性图像序列。
控制超声诊断仪通过分析模式对动态弹性图像序列中的每一幅弹性图像中的设定区域进行定量分析后，得到各动态弹性图像对应的弹性模量值。
较佳地，确定模块903在以所述肌肉动力学参数为横坐标，以同步的弹性模量值为纵坐标生成超声弹性肌动图曲线之后，根据超声弹性肌动图曲线估算肌肉生物力学参数之前，还用于：
确定超声弹性肌动图曲线与设定的力学模型曲线是否吻合，若吻合，执行根据所述超声弹性肌动图曲线估算肌肉生物力学参数的步骤。
较佳地，肌肉动力学图像包括位移图像，所述肌肉动力学参数包括肌肉拉伸长度，确定模块903按照以下方式估算肌肉生物力学参数：
获取超声弹性肌动图曲线对应的模型公式，以及模型公式中的各个参数。
根据参数确定所述单块骨骼肌在松弛状态下的长度、与单块骨骼肌在松弛状态下的长度对应的弹性模量值以及单块骨骼肌的延展度，延展度为单块骨骼肌的最大拉伸长度与单块骨骼肌在松弛状态下的长度之差。
在骨骼肌长度大于单块骨骼肌在松弛状态下的长度时，对模型公式按照高斯牛顿迭代法，确定单块骨骼肌的硬度指标。
与现有技术相比，本申请实施例是通过同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列，分别获取肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数，以及，动态弹性图像序列中各动态弹性图像对应的弹性模量值；以所述肌肉动力学参数为横坐标，以同步的所述弹性模量值为纵坐标生成超声弹性肌动图曲线，根据所述超声弹性肌动图曲线估算肌肉生物力学参数，可见，本申请实施例获得的肌肉动力学参数和弹性模量值是在单块骨骼肌连续拉伸时测得的变化参数，且最终生超声弹性肌动图曲线，可以由该曲线估算单块骨骼肌在连续拉伸过程中的生物力学参数，这样得到的生物力学参数是动态的且是定量的。
本申请实施例所提供的进行的根据超声弹性肌动图获取肌肉生物力学参数的方法计算机程序产品，包括存储了程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。
本申请实施例所提供的根据超声弹性肌动图获取肌肉生物力学参数的装置可以为设备上的特定硬件或者安装于设备上的软件或固件等。本申请实施例所提供的装置，其实现原理及产生的技术效果和前述方法实施例相同，为简要描述，装置实施例部分未提及之处，可参考前述方法实施例中相应内容。所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，前述描述的系统、装置和单元的具体工作过程，均可以参考上述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的实施例中，应该理解到，所揭露装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请提供的实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释，此外，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
最后应说明的是：以上所述实施例，仅为本申请的具体实施方式，用以说明本申请的技术方案，而非对其限制，本申请的保护范围并不局限于此，尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本申请实施例技术方案的精神和范围。都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应所述以权利要求的保护范围为准。

Claims (11)

1.一种根据超声弹性肌动图获取肌肉生物力学参数的方法，其特征在于，包括：
同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列；
分别获取所述动态肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数，以及，所述动态弹性图像序列中各弹性图像对应的弹性模量值；
以所述肌肉动力学参数为横坐标，以同步的所述弹性模量值为纵坐标生成超声弹性肌动图曲线，根据所述超声弹性肌动图曲线估算肌肉生物力学参数。
2.根据权利要求1所述的方法，其特征在于，所述肌肉动力学图像包括位移图像，所述肌肉动力学参数包括肌肉拉伸长度，在所述同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列之前，所述方法还包括：
确定初始状态下所述单块骨骼肌的肌肉组织伸长端的特征点；
计算所述特征点至所述单块骨骼肌的肌肉组织伸长端的对端的距离，得到所述骨骼肌的初始长度。
3.根据权利要求2所述的方法，其特征在于，确定所述单块骨骼肌的肌肉组织伸长端的特征点包括：
控制超声诊断仪通过成像模式扫描所述单块骨骼肌沿肌束方向的纵切面，将检测到的所述单块骨骼肌与肌腱的结合点作为所述特征点。
4.根据权利要求2或3所述的方法，其特征在于，所述获取所述动态肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数包括：
获取所述特征点在所述位移图像中相对初始状态下的特征点位置的位移变化的有效值；
将所述初始长度和所述位移变化的有效值之和记为所述位移图像对应的肌肉拉伸长度。
5.根据权利要求4所述的方法，其特征在于，所述获取所述特征点在所述位移图像中相对初始状态下的特征点位置的位移变化的有效值包括：
控制肌电仪监控所述单块骨骼肌及周围肌群的肌电信号；
当所述肌电信号的幅值小于设定值时，将所述位移图像中相对初始状态下的特征点位置的位移变化作为所述有效值。
6.根据权利要求1所述的方法，其特征在于，在所述同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列之前，所述方法还包括：
确定所述单块骨骼肌的肌腹位置，以在所述肌腹位置处采集所述动态弹性图像序列。
7.根据权利要求6所述的方法，其特征在于，所述确定所述单块骨骼肌的肌腹位置，包括：
控制超声诊断仪通过成像模式扫描所述骨骼肌中垂直于肌束方向的横切面，将横切面最大的位置作为所述单块骨骼肌的肌腹位置。
8.根据权利要求6或7所述的方法，其特征在于，所述获取动态弹性图像序列中各弹性图像对应的弹性模量值，包括：
在所述单块骨骼肌的肌腹位置沿肌束方向输出所述骨骼肌在连续拉伸时的所述动态弹性图像序列；
控制所述超声诊断仪通过分析模式对所述动态弹性图像序列中的每一幅弹性图像中的设定区域进行定量分析后，得到各弹性图像对应的弹性模量值。
9.根据权利要求1至3、6或7任一项所述的方法，其特征在于，在所述以所述肌肉动力学参数为横坐标，以同步的所述弹性模量值为纵坐标生成超声弹性肌动图曲线之后，根据所述超声弹性肌动图曲线估算肌肉生物力学参数之前，所述方法还包括：
确定所述超声弹性肌动图曲线与设定的力学模型曲线是否吻合，若吻合，执行所述根据所述超声弹性肌动图曲线估算肌肉生物力学参数的步骤。
10.根据权利要求9所述的方法，其特征在于，所述肌肉动力学图像包括位移图像，所述肌肉动力学参数包括肌肉拉伸长度，所述根据所述超声弹性肌动图曲线估算肌肉生物力学参数包括：
获取所述超声弹性肌动图曲线对应的模型公式，以及所述模型公式中的各个参数；
根据所述参数确定所述单块骨骼肌在松弛状态下的长度、与所述单块骨骼肌在松弛状态下的长度对应的弹性模量值以及所述单块骨骼肌的延展度，所述延展度为所述单块骨骼肌的最大拉伸长度与所述单块骨骼肌在松弛状态下的长度之差；
在所述骨骼肌长度大于所述单块骨骼肌在松弛状态下的长度时，对所述模型公式按照高斯牛顿迭代法，确定所述单块骨骼肌的硬度指标。
11.一种根据超声弹性肌动图获取肌肉生物力学参数的装置，其特征在于，包括：
采集模块，用于同步采集单块骨骼肌连续拉伸时的动态肌肉动力学图像序列和动态弹性图像序列；
获取模块，用于分别获取所述动态肌肉动力学图像序列中各肌肉动力学图像对应的肌肉动力学参数，以及，所述动态弹性图像序列中各弹性图像对应的弹性模量值；
确定模块，用于以所述肌肉动力学参数为横坐标，以同步的所述弹性模量值为纵坐标生成超声弹性肌动图曲线，根据所述超声弹性肌动图曲线估算肌肉生物力学参数。
CN201711322506.4A 2017-12-12 2017-12-12 一种根据超声弹性肌动图获取生物力学参数的方法及装置 Active CN107961038B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201711322506.4A CN107961038B (zh) 2017-12-12 2017-12-12 一种根据超声弹性肌动图获取生物力学参数的方法及装置
US16/482,802 US11452502B2 (en) 2017-12-12 2017-12-26 Method and device for acquiring biomechanical parameters based on ultrasonic elastomyogram
EP17934679.6A EP3566655B1 (en) 2017-12-12 2017-12-26 Method and apparatus for acquiring biomechanical parameter according to ultrasonic elasticity myogram
PCT/CN2017/118725 WO2019114034A1 (zh) 2017-12-12 2017-12-26 一种根据超声弹性肌动图获取生物力学参数的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711322506.4A CN107961038B (zh) 2017-12-12 2017-12-12 一种根据超声弹性肌动图获取生物力学参数的方法及装置

Publications (2)

Publication Number Publication Date
CN107961038A true CN107961038A (zh) 2018-04-27
CN107961038B CN107961038B (zh) 2020-04-14

Family

ID=61994275

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711322506.4A Active CN107961038B (zh) 2017-12-12 2017-12-12 一种根据超声弹性肌动图获取生物力学参数的方法及装置

Country Status (4)

Country Link
US (1) US11452502B2 (zh)
EP (1) EP3566655B1 (zh)
CN (1) CN107961038B (zh)
WO (1) WO2019114034A1 (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109620270A (zh) * 2018-12-14 2019-04-16 深圳先进技术研究院 一种肌张力评估方法及装置
CN109875609A (zh) * 2019-03-01 2019-06-14 清华大学 肌肉力学参数的测量装置及方法、弹性成像方法
CN110755095A (zh) * 2019-10-31 2020-02-07 清华大学 一种骨骼肌运动机能测试系统及其测试方法
WO2020118710A1 (zh) * 2018-12-14 2020-06-18 深圳先进技术研究院 一种肌张力评估方法及装置
CN111449641A (zh) * 2020-04-20 2020-07-28 浙江大学 一种基于光电信号检测的肌肉功能状态的评估装置和评估方法
CN112674791A (zh) * 2020-11-30 2021-04-20 深圳大学 肌肉超声弹性成像的优化方法及系统
CN113084814A (zh) * 2021-04-13 2021-07-09 中国科学院自动化研究所 基于分布位置优化实现肌肉骨骼机器人运动控制的方法
CN113693633A (zh) * 2021-10-28 2021-11-26 深圳高性能医疗器械国家研究院有限公司 骨骼肌被动弹性系数获取装置、设备及存储介质
CN113712597A (zh) * 2021-11-02 2021-11-30 深圳迈瑞动物医疗科技有限公司 一种马跟腱的检查方法和超声成像设备

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115131276B (zh) * 2022-03-03 2023-07-25 中国人民解放军总医院第四医学中心 骨骼受力分布的获取方法、装置、设备及存储介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090177084A1 (en) * 2006-05-25 2009-07-09 Takeshi Matsumura Ultrasonic Diagnostic Apparatus
JP2009219834A (ja) * 2008-03-13 2009-10-01 Seikosha:Kk Ａモード超音波画像二画像の差分評価による筋力評価装置
CN102920485A (zh) * 2012-10-30 2013-02-13 浙江大学 一种超声弹性成像中生物组织二维位移场的估计方法
CN103584884A (zh) * 2013-11-07 2014-02-19 中国科学院深圳先进技术研究院 肌力评估方法及装置、肌肉康复训练跟踪评估方法及系统
CN103690192A (zh) * 2013-12-11 2014-04-02 中国科学院深圳先进技术研究院 基于超声影像的跟腱形态学分析的系统和方法
CN105997079A (zh) * 2016-06-29 2016-10-12 中国人民解放军第二军医大学 一种通过超声融合mri成像的骨骼肌力在体测量方法
CN106202739A (zh) * 2016-07-14 2016-12-07 哈尔滨理工大学 一种骨骼肌力学行为多尺度建模方法
CN106264573A (zh) * 2016-07-26 2017-01-04 吉林大学 便携式肌肉力学参数及肌肉力在体超声检测装置及方法
CN106821308A (zh) * 2016-11-15 2017-06-13 深圳大学 一种骨骼肌刺激下收缩活动评估的方法和装置
CN106890009A (zh) * 2017-03-17 2017-06-27 深圳大学 一种骨骼肌体积变化的自动测量方法和装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FI20045484A0 (fi) * 2004-12-15 2004-12-15 Nomir Oy Menetelmä MRI-kuvantamisessa ja MRI-laite, jossa on laukaisin
US8147410B2 (en) 2009-03-23 2012-04-03 The Hong Kong Polytechnic University Method and apparatus for ultrasound imaging and elasticity measurement
WO2013148990A1 (en) * 2012-03-28 2013-10-03 University Of Washington Through Its Center For Commercialization Methods and systems for determining mechanical properties of a tissue
WO2015104582A1 (en) * 2014-01-08 2015-07-16 Amid S.R.L. Method and device for estimation of the elastic properties of tissues, particularly muscle tissues
CN104622511B (zh) 2015-01-26 2017-06-20 首都医科大学附属北京天坛医院 一种利用超声波检测骨骼肌机械力学参数的装置及超声波检测方法
CN106175831B (zh) 2016-09-28 2018-10-26 重庆理工大学 一种基于剪切波幅度和相位检测的组织粘弹性测量方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090177084A1 (en) * 2006-05-25 2009-07-09 Takeshi Matsumura Ultrasonic Diagnostic Apparatus
JP2009219834A (ja) * 2008-03-13 2009-10-01 Seikosha:Kk Ａモード超音波画像二画像の差分評価による筋力評価装置
CN102920485A (zh) * 2012-10-30 2013-02-13 浙江大学 一种超声弹性成像中生物组织二维位移场的估计方法
CN103584884A (zh) * 2013-11-07 2014-02-19 中国科学院深圳先进技术研究院 肌力评估方法及装置、肌肉康复训练跟踪评估方法及系统
CN103690192A (zh) * 2013-12-11 2014-04-02 中国科学院深圳先进技术研究院 基于超声影像的跟腱形态学分析的系统和方法
CN105997079A (zh) * 2016-06-29 2016-10-12 中国人民解放军第二军医大学 一种通过超声融合mri成像的骨骼肌力在体测量方法
CN106202739A (zh) * 2016-07-14 2016-12-07 哈尔滨理工大学 一种骨骼肌力学行为多尺度建模方法
CN106264573A (zh) * 2016-07-26 2017-01-04 吉林大学 便携式肌肉力学参数及肌肉力在体超声检测装置及方法
CN106821308A (zh) * 2016-11-15 2017-06-13 深圳大学 一种骨骼肌刺激下收缩活动评估的方法和装置
CN106890009A (zh) * 2017-03-17 2017-06-27 深圳大学 一种骨骼肌体积变化的自动测量方法和装置

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109620270A (zh) * 2018-12-14 2019-04-16 深圳先进技术研究院 一种肌张力评估方法及装置
WO2020118710A1 (zh) * 2018-12-14 2020-06-18 深圳先进技术研究院 一种肌张力评估方法及装置
CN109875609A (zh) * 2019-03-01 2019-06-14 清华大学 肌肉力学参数的测量装置及方法、弹性成像方法
CN109875609B (zh) * 2019-03-01 2020-12-22 清华大学 肌肉力学参数的测量装置及方法、弹性成像方法
CN110755095A (zh) * 2019-10-31 2020-02-07 清华大学 一种骨骼肌运动机能测试系统及其测试方法
CN111449641A (zh) * 2020-04-20 2020-07-28 浙江大学 一种基于光电信号检测的肌肉功能状态的评估装置和评估方法
CN112674791A (zh) * 2020-11-30 2021-04-20 深圳大学 肌肉超声弹性成像的优化方法及系统
CN112674791B (zh) * 2020-11-30 2023-08-29 深圳大学 肌肉超声弹性成像的优化方法及系统
CN113084814A (zh) * 2021-04-13 2021-07-09 中国科学院自动化研究所 基于分布位置优化实现肌肉骨骼机器人运动控制的方法
CN113084814B (zh) * 2021-04-13 2022-05-10 中国科学院自动化研究所 基于分布位置优化实现肌肉骨骼机器人运动控制的方法
CN113693633A (zh) * 2021-10-28 2021-11-26 深圳高性能医疗器械国家研究院有限公司 骨骼肌被动弹性系数获取装置、设备及存储介质
CN113693633B (zh) * 2021-10-28 2022-02-08 深圳高性能医疗器械国家研究院有限公司 骨骼肌被动弹性系数获取装置、设备及存储介质
CN113712597A (zh) * 2021-11-02 2021-11-30 深圳迈瑞动物医疗科技有限公司 一种马跟腱的检查方法和超声成像设备

Also Published As

Publication number Publication date
US20200229794A1 (en) 2020-07-23
US11452502B2 (en) 2022-09-27
EP3566655A1 (en) 2019-11-13
WO2019114034A1 (zh) 2019-06-20
EP3566655B1 (en) 2023-06-28
CN107961038B (zh) 2020-04-14
EP3566655A4 (en) 2020-06-03

Similar Documents

Publication Publication Date Title
CN107961038A (zh) 一种根据超声弹性肌动图获取生物力学参数的方法及装置
Ahn et al. Measurement and characterization of soft tissue behavior with surface deformation and force response under large deformations
Fung Elasticity of soft tissues in simple elongation
Stone et al. Modeling the motion of the internal tongue from tagged cine-MRI images
Sikdar et al. Dynamic ultrasound imaging applications to quantify musculoskeletal function
WO2007066300A2 (en) Method for in-vivo measurement of biomechanical properties of internal tissues
CN101835421B (zh) 特征量候选生成装置和特征量候选生成方法
Gerus et al. Ultrasound-based subject-specific parameters improve fascicle behaviour estimation in Hill-type muscle model
CN108113709A (zh) 一种肌肉萎缩程度的检测方法、装置、设备及存储介质
Shan et al. Inhomogeneous and anisotropic mechanical properties of the triceps surae muscles and aponeuroses in vivo during submaximal muscle contraction
Egorov et al. Biomechanical characterization of the pelvic floor using tactile imaging
CN106333697A (zh) 一种基于肌骨生物动力学的人体腹内压测试方法
JP2012213624A (ja) 身体能力判定装置及び身体能力判定方法
Ahmed et al. Corneal biomechanics assessment using high frequency ultrasound b-mode imaging
CN103690192A (zh) 基于超声影像的跟腱形态学分析的系统和方法
d’Aulignac et al. Towards a realistic echographic simulator
Stone et al. Modelling the internal tongue using principal strains
Mountjoy et al. Use of the fast orthogonal search method to estimate optimal joint angle for upper limb Hill-muscle models
JP4054879B2 (ja) 運動学習支援装置及び方法、運動学習支援プログラム及び該プログラムを記録した記録媒体
Gladilin et al. Towards a Realistic Simulation of Individual Facial Mimics.
Kaczmarek et al. Experimentally verified model of mechanomyograms recorded during single motor unit contractions
Medved Towards a virtual reality-assisted movement diagnostics-an outline
Otsuka et al. Validity of block start performance without arm forces or by kinematics-only methods
Vaz Muscle Synergy Analysis of Lower-Limb Movements
US20220230732A1 (en) Characterizing soft tissue stress for ameliorating injury in performing a process

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