一种骨骼肌体积变化的自动测量方法和装置制造方法及图纸

本发明专利技术提供一种骨骼肌体积变化的自动测量方法和装置。包括以下几个模块：数据采集模块：数据采集模块结合超声影像设备，采集骨骼肌横截面超声图像；图像预处理模块：用来降低超声图像有较大的散斑噪声，增强图像的兴趣区域轮廓；测量模块：提取预处理后的超声图像的兴趣区域轮廓；量化评估模块：量化评估模块是对自动提取的区域轮廓进行量化分析。本发明专利技术能够实现骨骼肌体积变化的自动测量，具有重要的临床应用价值，如可为临床上骨骼肌减少症患者的诊断和康复评估提供量化依据。

Automatic measuring method and device for volume change of skeletal muscle

The invention provides an automatic measuring method and device for volume change of skeletal muscle. Includes the following modules: data acquisition module, data acquisition module and ultrasonic imaging equipment, acquisition of skeletal muscle cross-sectional ultrasound image; image preprocessing module is used to reduce the speckle noise in ultrasound images is greatly enhanced, the region of interest image contours; contour measurement module: extracting region of interest ultrasound after image preprocessing; quantization evaluation module: quantitative evaluation module is the outline of the regional automatic extraction of quantitative analysis. The invention can automatically measure the volume change of skeletal muscles and has important clinical application value, such as providing quantitative basis for the diagnosis and rehabilitation evaluation of skeletal muscle wasting patients.

【技术实现步骤摘要】
一种骨骼肌体积变化的自动测量方法和装置
本专利技术涉及骨骼肌体积变化的测量方法和装置。
技术介绍
临床上医生通过使用卷尺测量骨骼肌减少症患者患肢围长作为骨骼肌变化的定性评估方法。目前骨骼肌减少症主要的定量诊断方法有双能x线吸收法(DXA)、CT、MRI，测量肌力等。用卷尺测量患肢围长的定性评估方法主观性强，不精确。双能x线吸收法测量系统复杂、有辐射，不能精确测量肌肉的横截面积和脂肪成分；CT、MRI是活体断层面肌量测量、肌肉密度和脂肪组织的测量、肌肉面积的评估，是目前最准确的测量方法，可作为诊断的金标准，但检查费用昂贵、测量分析复杂，且CT放射剂量比双能x线高，不适合频繁使用，限制了其临床应用。测量肌力的方法不精确，不能显示肌肉的形态学信息。在开展临床和科研研究时，应该综合考虑所查时间、费用、辐射剂量、测量重复性、准确性等多种因素，合理选择测量方法。
技术实现思路
为了解决以上技术问题，本专利技术提供一种骨骼肌体积变化的自动测量装置，包括以下几个模块：数据采集模块：数据采集模块结合超声影像设备，采集骨骼肌横截面超声图像；图像预处理模块：用来降低超声图像有较大的散斑噪声，增强图像的兴趣区域轮廓；测量模块：提取预处理后的超声图像的兴趣区域轮廓；量化评估模块：量化评估模块是对自动提取的区域轮廓进行量化分析。骨骼肌是人体中分布最广、数量最多的肌肉，约占人体体重的40％，在人体运动中扮演重要角色，负责人体的基本日常活动。骨骼肌减少症是指因骨骼肌体积缩小，质量和功能下降的一种常见症状，常发生于人体衰老、营养不良、肌肉不活动和各种疾病(包括神经肌肉疾病、癌症、细菌和病毒感染、慢性肺和肾疾病、糖尿病和药物副作用等)中。骨骼肌减少症具有患病率高、致残率高等特点。研究表明在临床门诊中，大约20％的门诊患者受肌肉减少症的影响。65-70岁老年人中，骨骼肌减少症的患病率为13％-24％，80岁以上老年人中患病率大于50％。骨骼肌是人体运动的动力，而骨骼肌质量和力量的降低，肌体活动功能下降，会导致跌倒、残疾甚至死亡等不良事件。而对骨骼肌体积变化的精确评估，是骨骼肌体积变化预防和康复过程中非常重要的一部分。肌肉体积(musclevolume,MV)能直接反映肌肉产生肌力的能力，肌肉体积减少会引起肌肉功能和物理性能的降低。肌肉横截面积是一个重要且可靠的肌肉体积计算方法，也是一个能直接反映肌肉产生肌力的能力的重要指标。如临床实践及相关科学研究通过卷尺测量患肢围长作为肌肉体积变化的一个粗略的诊断方法。本专利技术采用以上技术方案，其优点在于，该装置可以对骨骼肌的体积变化进行量化评估，其临床应用范围广泛，如为骨骼肌减少症的诊断和康复评定提供依据。基于此，本专利技术还提供一种骨骼肌体积变化的自动测量方法，包括以下步骤：步骤A:数据预处理：对采集的肌肉形态结构信息进行初步增强处理，强化兴趣区域的轮廓；步骤B:标记肌肉边界；步骤C：以标记结果为基准，创建图像训练集、训练标记集和测试集，提取训练集图像特征，将训练集特征和标记集共同训练成一个可以反映骨骼肌边缘轮廓与图像特征之间关系的机器学习模型，使之可自动分割骨骼肌边缘，计算肌肉横截面积；采用与训练集特征提取相同的类型和方法提取测试集特征，输入到训练集训练得到的机器学习模型中，提取骨骼肌边缘轮廓；步骤D：受试者骨骼肌采用上述机器学习方法提取出的横截面积，并从像素转化为物理尺寸，从而实现肌肉横截面积和体积的量化。优选的，所述步骤C中，包括步骤C1:创建数据样本，在预处理后的骨骼肌图像中，构建图像训练集和测试集；步骤C2:选择特征，在构建的训练集图像中提取设定区域作为感兴趣区域，对感兴趣区域进行轮廓标记构建训练图像轮廓标记集；步骤C3:训练分类器，从感兴趣区域中提取图像特征，用提取的图像特征与训练图像的标记集共同训练成一个可以反映骨骼肌的边缘轮廓与图像特征之间关系的机器学习模型；；步骤C4:测试，对测试集提取与训练集一致的设定区域作为感兴趣区域，提取测试集的图像特征，将提取到的测试图像特征输入到训练集训练得到的机器学习模型中，得到骨骼肌边缘轮廓。本专利技术采用以上技术方案，其优点在于，可以采用手动标记肌肉边界，手动标记的方法主观性较强、耗时，不适用于样本量较多的超声图像中。由于超声图像存在散斑噪声，传统的不基于学习的自动提取过程存在一定的难度，导致提取的精确度降低。而机器学习是一种让计算机更聪明、更个性的人工智能算法，已经成功应用于语音识别、计算机视觉、生物监测等领域，并且受到越来越多的研究学者的重视。通过对机器学习的研究，计算机能识别现有知识获取新知识并不断改善性能实现自我完善，机器学习算法具有较强的自适应性。本专利技术以标记结果为基准点，采用机器学习方法训练分类器，完善调整参数，使之能自动计算肌肉的横截面积。机器学习算法与传统的自动方法，它可以提供更高与手动结果的相似性、精确度，和手工标记方法相比，也改进了大样本数据处理的时间，提高了鲁棒性。优选的，所述步骤A中，采用自适应双边滤波对降低图像的散斑噪声，双边滤波的定义如下：I为原始图像，为平滑滤波后的输出图像，wD为空间域权系数，wR为灰度域权系数，N(m)表示m的邻域范围，n表示邻域的位置，其中，归一化函数Z为σd和σr为空间方差和灰度方差，是决定双边滤波权系数的参数。通过自适应的方式选择滤波参数，将双边滤波转换为自适应双边滤波，由于调节空间参数对噪声不敏感，通过自适应的方式选择灰度方差σr，因此，定义σr为而表示输入图像的估计噪声方差。优选的，所述步骤A中，数据预处理中，采用多尺度增强滤波，图像I(x，y)与高斯滤波器的二阶偏导数的卷积高斯函数G(x，y)为：图像上的每一个像素点f(x，y）的二阶偏导来构造Hessian矩阵：其中fxx、fxy、fyx、fyy分别表示二维灰度图像上像素点f(x，y)的四个二阶偏导数；根据偏导数性质：fxy＝fyx，那么Hessian的特征值有λ1、λ2(λ1＜λ2)，在尺度σ下图像的点p，基于Hessian矩阵的多尺度增强滤波函数定义为：其中参数β用来区别线状和块状物体，参数c和γ为平滑参数。超声图像预处理的方法有很多，例如中值滤波、自适应双边滤波和多尺度增强滤波在降低超声图像的散斑噪声时具有较好地效果。本专利技术进一步采用以上技术特征，其优点在于，由于超声图像具有较大的散斑噪声，需要对采集到的超声图像进行预处理以增强图像兴趣区域轮廓，降低超声图像的散斑噪声。在临床上，临床医生通过卷尺测量骨骼肌围长作为一个量化评估方法，目前并无自动测量方法的报道。本专利技术中的测量装置具有实时、便携、无辐射、成本低等优势。本专利技术能够实现骨骼肌体积变化的自动测量，具有重要的临床应用价值，如可为临床上骨骼肌减少症患者的诊断和康复评估提供量化依据。附图说明图1是本专利技术骨骼肌体积变化自动提取装置示意图。图2是本专利技术股四头肌横截面超声图像。图3是本专利技术手动标记股四头肌横截面边缘图像。图4是本专利技术数据预处理中的灰度调节后图像。图5是本专利技术数据预处理中的双边滤波图像。图6双边滤波后MVEF的图像。图7骨骼肌体积变化量化评估算法流程图。具体实施方式下面结合附图，对本专利技术的较优的实施例作进一步的详细说明：实施例1以股四头肌萎缩为例：股四头肌是人体中最大的肌肉，在人体本文档来自技高网...

【技术保护点】
一种骨骼肌体积变化的自动测量装置，其特征在于，包括以下几个模块：数据采集模块：数据采集模块结合超声影像设备，采集骨骼肌横截面超声图像；图像预处理模块：用来降低超声图像有较大的散斑噪声，增强图像的兴趣区域轮廓；机器学习模块：该模块通过及其学习方法，提取预处理后的超声图像的兴趣区域轮廓；量化评估模块：量化评估模块是对自动提取的区域轮廓进行量化分析。

【技术特征摘要】
1.一种骨骼肌体积变化的自动测量装置，其特征在于，包括以下几个模块：数据采集模块：数据采集模块结合超声影像设备，采集骨骼肌横截面超声图像；图像预处理模块：用来降低超声图像有较大的散斑噪声，增强图像的兴趣区域轮廓；机器学习模块：该模块通过及其学习方法，提取预处理后的超声图像的兴趣区域轮廓；量化评估模块：量化评估模块是对自动提取的区域轮廓进行量化分析。2.一种骨骼肌体积变化的自动测量方法，其特征在于，包括以下步骤：步骤A:数据预处理：对采集的肌肉形态结构信息进行初步增强处理，强化兴趣区域的轮廓；步骤B:标记肌肉边界；步骤C：以标记结果为基准，采用机器学习方法训练分类器，使之可自动计算肌肉横截面积；步骤D：对受试者肌肉采用上述机器学习得到的模型，提取横截面积，并从像素转化为物理尺寸，从而实现肌肉横截面积和体积的量化。3.如权利要求2所述的方法，其特征在于，所述步骤A中，采用自适应双边滤波对降低图像的散斑噪声，双边滤波的定义如下：I为原始图像，为平滑滤波后的输出图像，wD为空间域权系数，wR为灰度域权系数，N(m)表示m的邻域范围，n表示邻域的位置，其中，归一化函数Z为σd和σr为空间方差和灰度方差，是决定双边滤波权系数的参数。通过自适应的方式选择滤波参数，将双边滤波转换为自适应双边滤波，由于调节空间参数对噪声不敏感，通过自适应的方式选择灰度方差σr，因此，定义σr为而表示输入图像的估计噪声方差。4.如权利...

【专利技术属性】
技术研发人员：周永进，石文秀，杨晓娟，张树，徐井旭，
申请(专利权)人：深圳大学，
类型：发明
国别省市：广东,44