CN112168211B - 一种腹部超声图像的脂肪厚度和肌肉厚度测量方法及系统 - Google Patents
一种腹部超声图像的脂肪厚度和肌肉厚度测量方法及系统 Download PDFInfo
- Publication number
- CN112168211B CN112168211B CN202010223512.XA CN202010223512A CN112168211B CN 112168211 B CN112168211 B CN 112168211B CN 202010223512 A CN202010223512 A CN 202010223512A CN 112168211 B CN112168211 B CN 112168211B
- Authority
- CN
- China
- Prior art keywords
- image
- layer
- fat
- muscle
- layer image
- 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
- 210000003205 muscle Anatomy 0.000 title claims abstract description 97
- 238000000034 method Methods 0.000 title claims abstract description 40
- 210000001015 abdomen Anatomy 0.000 title claims abstract description 19
- 210000001519 tissue Anatomy 0.000 claims abstract description 69
- 230000000007 visual effect Effects 0.000 claims abstract description 21
- 238000010586 diagram Methods 0.000 claims abstract description 20
- 230000003187 abdominal effect Effects 0.000 claims abstract description 19
- 210000003195 fascia Anatomy 0.000 claims description 42
- 230000011218 segmentation Effects 0.000 claims description 22
- 238000012545 processing Methods 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 10
- 238000013135 deep learning Methods 0.000 claims description 10
- 238000013507 mapping Methods 0.000 claims description 10
- 238000012549 training Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 238000009499 grossing Methods 0.000 claims description 9
- 238000010606 normalization Methods 0.000 claims description 9
- 238000002604 ultrasonography Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 5
- 239000003086 colorant Substances 0.000 description 6
- 210000004003 subcutaneous fat Anatomy 0.000 description 6
- 210000000577 adipose tissue Anatomy 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004040 coloring Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 208000008589 Obesity Diseases 0.000 description 1
- 210000003489 abdominal muscle Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 235000020824 obesity Nutrition 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
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/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- 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)
- Image Processing (AREA)
Abstract
本发明公开了一种腹部超声图像的脂肪厚度和肌肉厚度测量方法，包括采集腹部超声图像，选取部分作为待处理图像，进行缩放得到待分割图像，分割待分割图像的组织结构，得到第一组织结构图像，将第一组织结构图像进行缩放处理，得到与待处理图像大小一致的第二组织组织图像，利用第二组织结构图像中的第二脂肪层图像和第二肌肉层图像，计算脂肪层厚度和肌肉层厚度；还包括对待处理图像和第二组织结构图像进行伪彩化处理，得到可视化结构图。本发明能够简单高效的计算脂肪层厚度和肌肉层厚度，还能得到清晰直观的展示腹部组织结构的可视化结构图。
Description
技术领域
本发明涉及图像处理技术领域，特别是一种腹部超声图像的脂肪厚度和肌肉厚度测量方法及系统。
背景技术
肥胖及其导致的疾病越来越为人们所重视，对自己身体的脂肪含量也越来越关注。同时人们对健身越来越热爱，在健身中每日关注自己的皮下脂肪和腹肌厚度可以有提高健身的积极性，也有利于有效安排训练计划。皮下脂肪占身体总脂肪的40–60％，因而检测皮下脂肪厚度可以反映身体各部分的脂肪分布。超声成像技术可以较准确测量脂肪组织厚度，但是超声图像专业性很强，普通没有经过训练的人很难分清各层组织结构。
专利号为CN109009224A的发明专利，公开了一种方便快捷的皮下脂肪测量仪，通过设计硬件电路，利用超声波在皮下脂肪层和肌肉交界处的反射，利用超声波发射和回波的时间间隔和在脂肪中的传播速度计算皮下脂肪厚度。但是只能测量脂肪的厚度，不能测量其他组织结构的厚度。而且因为是在信号上做的，相应的信号传输还依赖于硬件电路的功能。
专利号为CN103027713A的发明专利，公开了一种基于超声图像的肌肉厚度测量方法和系统，利用通过对选取的多个跟踪窗口进行跟踪，通过跟踪算法确定在后续每帧图像中的多个跟踪窗口的位置，计算每帧图像中的经过中心点处理的跟踪窗口的位置与边缘检测法处理后的跟踪窗口的位置之间的最大垂直距离，作为肌肉厚度值。测量复杂，且数据单一，不能全面的反映被测量位置的组织结构情况。
发明内容
本发明的目的在于：针对现有技术存在的无法分清结构和测量数据单一的问题，提供一种腹部超声图像的脂肪厚度和肌肉厚度测量方法，利用腹部超声图像的组织结构分割图像计算脂肪层厚度和肌肉层厚度，还能够将腹部超声图像处理成用不同色彩区分不同组织结构的可视化结构图。
为了实现上述目的，本发明采用的技术方案为：
一种腹部超声图像的脂肪厚度和肌肉厚度测量方法，包括以下步骤：
A1、采集腹部超声图像，选取部分图像作为待处理图像，对所述待处理图进行缩放，得到大小一致的待分割图像；
A2、利用基于深度学习训练得到的超声腹部组织分割模型对所述待分割图像中的组织结构进行分割，分别得到第一皮肤层图像、第一脂肪层图像、第一筋膜层图像、第一肌肉层图像和第一背景层图像；
A3、将所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像进行缩放处理，得到与所述待处理图像大小一致的第二皮肤层图像、第二脂肪层图像、第二筋膜层图像、第二肌肉层图像和第二背景层图像；
A4、分别利用所述所述第二脂肪层图像和所述第二肌肉层图像，计算脂肪层厚度和肌肉层厚度。
作为本发明的优选方案，所述步骤A2分割得到的所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像均为大小一致的二值图像，并且图像中包含属于该层组织结构的组织像素信息。
作为本发明的优选方案，所述步骤A3包括：利用所述第二脂肪层图像中的脂肪像素信息和图像的宽度信息计算所述脂肪厚度；以及利用所述第二肌肉层图像中的肌肉像素信息和宽度信息计算所述脂肪厚度。
作为本发明的优选方案，还包括对所述待处理图像进行伪彩色处理，得到可视化结构图；所述伪彩色处理：
对所述待处理图像归一化处理：首先计算所述待处理图像Img的均值mean和标准差std，根据下列公式进行归一化，得到归一化图像Img_nor：
Img_nor(i,j)＝(Img(i,j)-mean)/std
其中，Img_nor(i,j)为所述归一化图像的像素值；Img(i,j)为所述待处理图像的像素值，mean为所述待处理图像的均值；std为所述待处理图像的标准差；
对所述归一化图像Img_nor进行Gamma校正，得到校正图像Img_gamma：
Img_gamma＝(Img_nor)γ
采用低通滤波，平滑所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像的边缘结构，得到各层组织结构对应的平滑图像；
将所述校正图像以各层对应的平滑图像为权重，为各层组织结构应用不同伪彩色映射；得到所述待处理图像对应的所述可视化结构图。
本发明还提供一种腹部超声图像的脂肪厚度和肌肉厚度测量系统，包括：
图像采集模块：用于采集腹部超声图像，选取部分图像作为待处理图像，对所述待处理图进行缩放，得到大小一致的待分割图像；
深度学习组织分割模块：用于分割所述待分割图像的各层组织结构，分别得到第一皮肤层图像、第一脂肪层图像、第一筋膜层图像、第一肌肉层图像和第一背景层图像；
缩放模块，所述缩放模块用于将所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像进行缩放处理，得到与所述待处理图像大小一致的第二皮肤层图像、第二脂肪层图像、第二筋膜层图像、第二肌肉层图像和第二背景层图像；
计算模块：用于根据所述第二脂肪层图像和所述第二肌肉层图像，分别计算脂肪层厚度肌肉层厚度。
作为本发明的优选方案，所述系统还包括伪彩化模块，所述伪彩化模块包括：
归一化模块：用于对所述待处理图像进行归一化处理，得到归一化图像；
Gamma校正模块：用于对所述归一化图像进行校正，得到校正图像；
低通滤波模块：用于平滑所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像的边缘结构；
伪彩映射模块：用于将所述校正图像以所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像平滑边缘后的模板为权重，为各个组织层结构应用不同的伪彩色映射，得到可视化结构图。
作为本发明的优选方案，所述系统还包括显示模块：用于显示所述可视化结构图、所述脂肪层厚度和所述肌肉层厚度。
本发明还提供一种电子设备，包括至少一个处理器，以及与所述至少一个处理器通信连接的存储器；所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够实现上述方法。
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述方法。
综上所述，由于采用了上述技术方案，本发明的有益效果是：
1、本发明利用腹部超声图像分割出的组织结构图像，计算脂肪层厚度和肌肉层厚度，计算简单、高效。
2、本发明计算脂肪层厚度、肌肉层厚度的同时，还可通过伪彩化处理得到组织结构的可视化结构图，能够更直观的表明脂肪和肌肉的分布情况。
附图说明
图1是本发明的方法流程图。
图2是本发明的待处理图像。
图3是本发明的皮肤层图像
图4是本发明的脂肪层图像。
图5是本发明的筋膜层图像
图6是本发明的肌肉层图像。
图7是本发明的可视化结构图。
图8是本发明的系统结构示意图。
图9是本发明提供的电子设备终端示意图。
图标：S1-皮肤层；S2-脂肪层；S3-筋膜层；S4-肌肉层；S5-背景层。
具体实施方式
下面结合附图，对本发明作详细的说明。
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
实施例1
参照图1，一种腹部超声图像的脂肪厚度和肌肉厚度测量方法，包括以下步骤：
采集腹部超声图像，提取待处理图像；
A1、采集腹部超声图像，选取部分图像作为待处理图像，对所述待处理图进行缩放，得到大小一致的待分割图像；
在采集到的腹部超声图像中选取部分图像，缩放成与所述超声腹部组织分割模型的输入所需的大小一致的单通道图像，参照图2，为本实施例选取的所述待处理图像，因为超声图像为灰度图像，图中各层组织结构难以辨别。本实施例的待处理图像为1024*512，经缩放处理后，得到512*512的所述待分割图像。
分割待分割图像的组织结构；
A2、利用基于深度学习得到的超声腹部组织分割模型对所述待分割图像中的组织结构进行分割，分别得到第一皮肤层图像、第一脂肪层图像、第一筋膜层图像、第一肌肉层图像和第一背景层图像；
分割得到的所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像均为与所述待分割图像大小一致的二值图像，并且图像中包含属于该层组织结构的组织像素信息。
本实施例所述的基于深度学习得到的超声图像腹部组织分割模型，首先搭建基于U-Net的深度学习网络，利用腹部超声图像和对应的各层分割图像的二值图像进行训练，得到训练好的超声腹部组织分割模型。网络结构采用改进的U-Net网络架构，主体结构与U-Net网络类似。与U-net网络不同的是：
1、网络输入为一个512*512*1单通道的图像，输出为一个512*512*5的张量。
2、网络包含4个下采样层与4个上采样层，其中每个下采样层的卷积层均采用与MobileNetV2一致的深度可分离卷积的结构。
3、上采样层与下采样层的每个卷积操作后都接一个BatchNormalization层。
训练网络采用的损失函数为：
Total_loss＝Dice_loss+0.5*Focal_loss
其中：Dice_loss＝Dice_skin+Dice_fat+Dice_fascia+Dice_muscle+0.5*Dice_background；Focal_loss＝Focal_skin+Focal_fat+Focal_fascia+Focal_muscle+0.5*Focal_background。
其中：Dice_skin，Dice_fat，Dice_fascia，Dice_muscle，Dice_background均为标准的Dice coefficient。
Dice coefficient定义如下：
Dice系数，根据Lee Raymond Dice命名，是一种集合相似度度量函数，通常用于计算两个样本的相似度(值范围为[0,1])：
|X∩Y|是X和Y之间的交集；|X|和|Y|分别表示X和Y的元素个数，其中，分子中的系数2，是因为分母存在重复计算X和Y之间的共同元素的原因，
在本实施例中，X是标注的各层组织结构的区域，Y是网络输出的预测结果。
Focal_skin，Focal_fat，Focal_fascia，Focal_muscle，Focal_background均为标准的Focal Loss,其中公式中的alpha设置为0.25,gamma设置为2.0。
基于深度学习训练而获得的超声腹部组织分割模型。此模型的输入为所述待分割图像，输出为所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像。
本实施例每一输出层都为与输入层大小一致的二值图像，图像中为1的像素表征该像素属于对应的结构层，为0的像素表示不属于对应的结构层；
将分割图像进行缩放处理，得到与待处理图像大小一致的组织图像
A3、将所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像进行缩放处理，得到与所述待处理图像大小一致的第二皮肤层图像、第二脂肪层图像、第二筋膜层图像、第二肌肉层图像和第二背景层图像。
该步骤为所述步骤A1中缩放处理的逆缩放，将分割得到的512*512的所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像，再次缩放为与所述待处理图像大小一致的1024*512的所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像。
参照图3、图4图5和图6，分别为本实施例分割得到的所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像和所述第二肌肉层图像。
A4、利用所述第二脂肪层图像和所述第二肌肉层图像分别计算脂肪厚度和肌肉厚度；
按照以下公式进行计算：
Count/W*mm_pixel
其中，Count为要计算的层对应的图像中属于该层组织结构的像素个数；W为该层组织结构图像的宽度，即为该层组织结构的宽度；mm_pixel为一个像素对应的毫米。
本实施例采用像素值1表示分别统计所述第二脂肪层图像和所述第二肌肉层图像中像素为1的像素个数，分别得到脂肪层的像素个数和肌肉层的像素个数。
则脂肪层的厚度为：
脂肪层的厚度＝Fat_Count/W*mm_pixel
其中，Fat_Count为第二脂肪层图像中像素值为1的像素个数；W为第二脂肪层图像的宽度；mm_pixel为一个像素对应多少毫米。
肌肉层的厚度为：
肌肉层的厚度＝Muscle_Count/W*mm_pixel
其中，Muscle_Count第二肌肉层图像中像素值为1的像素个数；W为第二肌肉层图像的宽度；mm_pixel为一个像素对应的毫米。
本方法分割得到的各层组织结构图像进行缩放处理后的大小均一致，所以各层组织结构图像的宽度信息相同，即各层组织结构的宽度相同。
本方法还包括利用所述待处理图像Img、所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像进行伪彩色处理，得到可视化结构图；
将所述待处理图像做归一化处理：首先计算所述待处理图像Img的均值mean和标准差std，根据下列公式进行归一化，得到归一化图像Img_nor：
Img_nor(i,j)＝(Img(i,j)-mean)/std
其中，Img_nor(i,j)所述归一化图像的像素值；Img(i,j)为所述待处理图像的像素值，mean为所述待处理图像的均值；std为所述待处理图像的标准差；
为了增加各个组织层次的对比度，对所述归一化图像Img_nor进行Gamma校正，得到校正图像Img_gamma：
Img_gamma＝(Img_nor)γ
本实施例中各个层次的γ取值为：皮肤层1.2，脂肪层1.1，筋膜层1.0，肌肉层1.3，背景层1.07。
采用低通滤波，平滑所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像的边缘结构。
将所述校正图像Img_gamma以平滑边缘结构后的各层组织结构图像为权重，为各层组织结构应用不同伪彩色映射；不同组织层结构的伪彩选择以接近现实的颜色为原则，得到所述可视化结构图。
上述实施例所述的方法能够得到得到脂肪层厚度和肌肉层厚度，参照图7，为本实施例的可视化结构图，用不同颜色表示不同组织结构的可视化结构示意图，可清楚的分辨图中的皮肤层(S1)、脂肪层(S2)、筋膜层(S3)、肌肉层(S4)和背景层(S5)。
参照图8，一种腹部超声图像的脂肪厚度和肌肉厚度测量系统，包括：
图像采集模块：用于采集腹部超声图像，选取部分图像作为待处理图像，对所述待处理图进行缩放，得到大小一致的待分割图像
本实施例利用现有的超声设备，进行图像采集，将超声探头置于肚脐上1cm，垂直于涂有超声耦合剂的皮肤表面(不要用力按压)，扫描获得一系列超声腹部图像后冻结超声系统，从超声电影回放中选取一帧组织结构清晰的超声图像，将选取的图像作为所述待处理图像，将所述待处理图像缩放成大小一致的单通道图像，作为组织结构分割模型的输入图像，即为所述待分割图像。
深度学习组织分割模块：用于分割所述待处理图像的各层组织结构，分别得到第一皮肤层图像、第一脂肪层图像、第一筋膜层图像、第一肌肉层图像和第一背景层图像；
缩放模块，所述缩放模块用于将所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像进行缩放处理，得到与所述待处理图像大小一致的第二皮肤层图像、第二脂肪层图像、第二筋膜层图像、第二肌肉层图像和第二背景层图像；
伪彩化模块：对所述待处理图像和所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像进行伪彩色处理，得到用不同伪彩色区分皮肤层、脂肪层、筋膜层、肌肉层和背景层的可视化结构图；
所述伪彩化模块包括归一化模块、Gamma校正模块、低通滤波模块和伪彩映射模块；
所述归一化模块，利用所述待处理图像的均值和标准差进行归一化处理，得到归一化图像。该模块减小医学图片由于光线不均匀造成的干扰。
所述Gamma模块对所述归一化图像进行校正，得到校正图像，增加了各个组织层次的对比度。
所述低通滤波模块，用于处理所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像，平滑各层的边缘结构；有利于更美观的可视化。
所述伪彩映射模块，用于将所述校正图像以所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像平滑边缘后的模板为权重，为各个组织层结构应用不同的伪彩色映射，得到可视化结构图。
计算模块，用于计算所述脂肪厚度和所述肌肉厚度。直接在所述第二脂肪层图像和所述第二肌肉层图像中，进行计算。
显示模块：用于显示所述伪彩化模块的所述可视化结构图、所述计算模块计算得到的脂肪增厚度和肌肉层厚度。
本实施例所述的显示模块可以是显示屏幕、平板电脑或者手机设备等。
综上所述，本实施例所述的方法及系统可以将腹部的超声图像，进行组织结构分割，利用分割得到的组织结构图计算可得到脂肪厚度和肌肉厚度，还可对超声图像进行伪彩化处理，得到可视化结构图，图中不同的色彩代表不同的层次结构，便于辨别组织层次，能够更加直观的了解脂肪和肌肉的分布状况。
如图9所示，本发明还提供了一种电子设备，包括至少一个处理器，以及与所述至少一个处理器通信连接的存储器；所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行前述任一实施例所述的方法。所述输入输出接口可以包括显示器、键盘、鼠标、以及USB接口，用于输入输出数据；电源用于为电子设备提供电能。
本发明还提供一种计算机可读存储介质，本申请实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random AccessMemory)等。
需要说明的是，对于装置及电子设备实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
最后需要注意的是，公布实施例的目的在于帮助进一步理解本发明，但是本领域的技术人员可以理解：在不脱离本发明及所附的权利要求的精神和范围内，各种替换和修改都是可能的。因此，本发明不应局限于实施例所公开的内容，本发明要求保护的范围以权利要求书界定的范围为准。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种腹部超声图像的脂肪厚度和肌肉厚度测量方法，其特征在于，包括以下步骤：
A1、采集腹部超声图像，选取部分图像作为待处理图像，对所述待处理图像进行缩放，得到大小一致的待分割图像；
A2、利用基于深度学习训练得到的超声腹部组织分割模型对所述待分割图像中的组织结构进行分割，分别得到第一皮肤层图像、第一脂肪层图像、第一筋膜层图像、第一肌肉层图像和第一背景层图像；
A3、将所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像进行缩放处理，得到与所述待处理图像大小一致的第二皮肤层图像、第二脂肪层图像、第二筋膜层图像、第二肌肉层图像和第二背景层图像；
A4、分别利用所述第二脂肪层图像和所述第二肌肉层图像，计算脂肪层厚度和肌肉层厚度；
所述超声腹部组织分割模型的网络结构采用改进的U-net网络结构，其中，所述超声腹部组织分割模型的网络输入为一个512*512*1单通道的图像，输出为一个512*512*5的张量；包含4个下采样层与4个上采样层，其中每个下采样层的卷积层均采用与MobileNet V2一致的深度可分离卷积的结构；上采样层与下采样层的每个卷积操作后都接一个BatchNormalization 层；
训练网络采用的损失函数为：
Total_loss=Dice_loss+0.5*Focal_loss
其中：Dice_loss=Dice_skin+Dice_fat+Dice_fascia+Dice_muscle +0.5*Dice_background；Focal_loss=Focal_skin+Focal_fat+ Focal_fascia+ Focal_muscle+0.5*Focal_background，
其中：Dice_skin，Dice_fat，Dice_fascia，Dice_muscle，Dice_background均为标准的Dice coefficient；
Focal_skin，Focal_fat，Focal_fascia，Focal_muscle，Focal_background，均为标准的Focal Loss；
步骤A4中，脂肪层厚度计算公式为：
脂肪层的厚度=（Fat_Count/W1）* mm_pixel
其中，Fat_Count为第二脂肪层图像中像素值为1 的像素个数；W1 为第二脂肪层图像的宽度；mm_pixel 为一个像素对应多少毫米；
肌肉层的厚度为：
肌肉层的厚度=（Muscle_Count/W2）* mm_pixel
其中，Muscle_Count 第二肌肉层图像中像素值为1的像素个数；W2为第二肌肉层图像的宽度；mm_pixel为一个像素对应多少毫米；
还包括对所述待处理图像进行伪彩色处理，得到可视化结构图；所述伪彩色处理：
对所述待处理图像归一化处理：首先计算所述待处理图像的均值mean 和标准差std，根据下列公式进行归一化，得到归一化图像Img_nor ：
Img_nor（i,j）=（Img（i,j）-mean）/std
其中，Img_nor（i,j）为所述归一化图像的像素值；Img（i,j）为所述待处理图像的像素值，mean为所述待处理图像的均值；std 为所述待处理图像的标准差；
对所述归一化图像Img_nor 进行Gamma 校正，得到校正图像Img_gamma ：
；
的取值为皮肤层1.2，脂肪层1.1，筋膜层1.0，肌肉层1.3，背景层1.07；采用低通滤波，平滑所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像的边缘结构，得到各层组织结构对应的平滑图像；
将所述校正图像以各层对应的平滑图像为权重，为各层组织结构应用不同伪彩色映射；得到所述待处理图像对应的所述可视化结构图。
2.根据权利要求1所述的方法，其特征在于，所述步骤A2分割得到的所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像均为大小一致的二值图像，并且各个层图像中包含属于该层组织结构的组织像素信息。
3.一种腹部超声图像的脂肪厚度和肌肉厚度测量系统，其特征在于，包括：
图像采集模块，用于采集腹部超声图像，选取部分图像作为待处理图像，对所述待处理图像进行缩放，得到大小一致的待分割图像；
深度学习组织分割模块，用于分割所述待分割图像的各层组织结构，分别得到第一皮肤层图像、第一脂肪层图像、第一筋膜层图像、第一肌肉层图像和第一背景层图像；
缩放模块，所述缩放模块用于将所述第一皮肤层图像、所述第一脂肪层图像、所述第一筋膜层图像、所述第一肌肉层图像和所述第一背景层图像进行缩放处理，得到与所述待处理图像大小一致的第二皮肤层图像、第二脂肪层图像、第二筋膜层图像、第二肌肉层图像和第二背景层图像；
计算模块，用于根据所述第二脂肪层图像和所述第二肌肉层图像，分别计算脂肪层厚度肌肉层厚度；
所述深度学习组织分割模块利用基于深度学习训练得到的超声腹部组织分割模型对所述待分割图像中的组织结构进行分割，
所述超声腹部组织分割模型的网络结构采用改进的U-net网络结构，其中，所述超声腹部组织分割模型的网络输入为一个512*512*1单通道的图像，输出为一个512*512*5的张量；包含4个下采样层与4个上采样层，其中每个下采样层的卷积层均采用与MobileNet V2一致的深度可分离卷积的结构；上采样层与下采样层的每个卷积操作后都接一个BatchNormalization 层；
训练网络采用的损失函数为：
Total_loss=Dice_loss+0.5*Focal_loss
其中：Dice_loss=Dice_skin+Dice_fat+Dice_fascia+ Dice_muscle+0.5 *Dice_background；Focal_loss=Focal_skin+Focal_fat+Focal_fascia+ Focal_muscle+0.5*Focal_background；
其中：Dice_skin ，Dice_fat ，Dice_fascia，Dice_muscle ，Dice_background均为标准的Dice coefficient；
Focal_skin，Focal_fat，Focal_fascia，Focal_muscle，Focal_background，均为标准的Focal Loss；
步骤A4中，脂肪层厚度计算公式为：
脂肪层的厚度= （Fat_Count/W1）*mm_pixel
其中，Fat_Count为第二脂肪层图像中像素值为1的像素个数；W1为第二脂肪层图像的宽度；mm_pixel为一个像素对应多少毫米；
肌肉层的厚度为：
肌肉层的厚度=（Muscle_Count/W2）*mm_pixel
其中，Muscle_Count第二肌肉层图像中像素值为1的像素个数；W为第二肌肉层图像的宽度；mm_pixel为一个像素对应多少毫米；
系统还包括伪彩化模块，所述伪彩化模块包括：
归一化模块：用于对所述待处理图像进行归一化处理，得到归一化图像；
对所述待处理图像归一化处理：首先计算所述待处理图像Img 的均值mean 和标准差std，根据下列公式进行归一化，得到归一化图像Img_nor ：
Img_nor（i,j）=（Img（i,j）-mean）/std
其中，Img_nor（i,j）为所述归一化图像的像素值；Img（i,j）为所述待处理图像的像素值，mean为所述待处理图像的均值；std为所述待处理图像的标准差；
Gamma校正模块：用于对所述归一化图像进行校正，得到校正图像，具体步骤包括：对所述归一化图像Img_nor 进行Gamma 校正，得到校正图Img_gamma ：
；
的取值为皮肤层1.2，脂肪层1.1，筋膜层1.0，肌肉层1.3，背景层1.07；
低通滤波模块：用于平滑所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像的边缘结构，得到各层组织结构对应的平滑图像；
伪彩映射模块：用于将所述校正图像以所述第二皮肤层图像、所述第二脂肪层图像、所述第二筋膜层图像、所述第二肌肉层图像和所述第二背景层图像平滑边缘后的模板为权重，为各个组织层结构应用不同的伪彩色映射，得到可视化结构图。
4.根据权利要求3所述的系统，其特征在于，所述系统还包括显示模块：用于显示所述可视化结构图、所述脂肪层厚度和所述肌肉层厚度。
5.一种电子设备，其特征在于，包括至少一个处理器，以及与所述至少一个处理器通信连接的存储器；所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1 至2中任一项所述的方法。
6.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1 至2任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010223512.XA CN112168211B (zh) | 2020-03-26 | 2020-03-26 | 一种腹部超声图像的脂肪厚度和肌肉厚度测量方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010223512.XA CN112168211B (zh) | 2020-03-26 | 2020-03-26 | 一种腹部超声图像的脂肪厚度和肌肉厚度测量方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112168211A CN112168211A (zh) | 2021-01-05 |
| CN112168211B true CN112168211B (zh) | 2024-08-09 |
Family
ID=73919091
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010223512.XA Active CN112168211B (zh) | 2020-03-26 | 2020-03-26 | 一种腹部超声图像的脂肪厚度和肌肉厚度测量方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112168211B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113643319B (zh) * | 2021-10-18 | 2021-12-10 | 成都汇声科技有限公司 | 在图像上标定脂肪厚度起算线及脂肪厚度测量方法 |
| CN116309385B (zh) * | 2023-02-27 | 2023-10-10 | 之江实验室 | 基于弱监督学习的腹部脂肪与肌肉组织测量方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102580260A (zh) * | 2010-12-17 | 2012-07-18 | 通用电气公司 | 超声图像中自动识别和分割不同组织类型的系统和方法 |
| CN106097283A (zh) * | 2016-07-20 | 2016-11-09 | 南京邮电大学 | 一种基于人类视觉系统特性的多尺度x线图像增强方法 |
| CN110211140A (zh) * | 2019-06-14 | 2019-09-06 | 重庆大学 | 基于3D残差U-Net和加权损失函数的腹部血管分割方法 |
Family Cites Families (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100036246A1 (en) * | 2008-08-07 | 2010-02-11 | Leonid Kushculey | Automatic fat thickness measurements |
| CN103827874B (zh) * | 2011-09-26 | 2017-02-22 | 皇家飞利浦有限公司 | 医学图像系统和方法 |
| CN102920436A (zh) * | 2012-10-25 | 2013-02-13 | 深圳先进技术研究院 | 利用霍夫变换的肌肉图像处理方法 |
| US9642593B2 (en) * | 2014-09-19 | 2017-05-09 | MuscleSound, LLC | System and method for non-invasive determination of human body fat |
| CN106355574B (zh) * | 2016-08-31 | 2018-11-16 | 上海交通大学 | 一种基于深度学习的腹内脂肪分割方法 |
| CN108665461B (zh) * | 2018-05-09 | 2019-03-12 | 电子科技大学 | 一种基于fcn和迭代声影校正的乳腺超声图像分割方法 |
| CN109325954B (zh) * | 2018-09-18 | 2021-08-10 | 北京旷视科技有限公司 | 图像分割方法、装置及电子设备 |
| CN109448006B (zh) * | 2018-11-01 | 2022-01-28 | 江西理工大学 | 一种注意力机制u型密集连接视网膜血管分割方法 |
| CN109523561B (zh) * | 2018-12-19 | 2021-04-02 | 睿佳(武汉)软件科技有限公司 | 自动腹内肌肉和脂肪图像分割方法 |
| CN109993735A (zh) * | 2019-03-29 | 2019-07-09 | 成都信息工程大学 | 基于级联卷积的图像分割方法 |
| CN110223300A (zh) * | 2019-06-13 | 2019-09-10 | 北京理工大学 | Ct图像腹部多器官分割方法及装置 |
| CN110310292B (zh) * | 2019-06-28 | 2021-02-02 | 浙江工业大学 | 一种手腕部参照骨分割方法 |
| CN110610497B (zh) * | 2019-08-05 | 2023-06-16 | 史记生物技术(南京)有限公司 | 一种基于ct图像处理的活体猪胴体组织含量测定方法 |
| CN110415246B (zh) * | 2019-08-06 | 2023-05-16 | 东北大学 | 一种腹部脂肪成分的分析方法 |
| CN110517235B (zh) * | 2019-08-19 | 2021-10-19 | 苏州大学 | 一种基于GCS-Net进行OCT图像脉络膜自动分割方法 |
| CN110853049A (zh) * | 2019-10-17 | 2020-02-28 | 上海工程技术大学 | 一种腹部超声图像分割方法 |
-
2020
- 2020-03-26 CN CN202010223512.XA patent/CN112168211B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102580260A (zh) * | 2010-12-17 | 2012-07-18 | 通用电气公司 | 超声图像中自动识别和分割不同组织类型的系统和方法 |
| CN106097283A (zh) * | 2016-07-20 | 2016-11-09 | 南京邮电大学 | 一种基于人类视觉系统特性的多尺度x线图像增强方法 |
| CN110211140A (zh) * | 2019-06-14 | 2019-09-06 | 重庆大学 | 基于3D残差U-Net和加权损失函数的腹部血管分割方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112168211A (zh) | 2021-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11373303B2 (en) | Systems and methods for ultrasound imaging | |
| US11229419B2 (en) | Method for processing 3D image data and 3D ultrasonic imaging method and system | |
| US9679375B2 (en) | Ovarian follicle segmentation in ultrasound images | |
| CN109949254B (zh) | 穿刺针超声图像增强方法及装置 | |
| US9585636B2 (en) | Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method | |
| CN108805946B (zh) | 用于为二维超声图像绘阴影的方法和系统 | |
| US20170124700A1 (en) | Method and system for measuring a volume from an ultrasound image | |
| CN111820948B (zh) | 胎儿生长参数测量方法、系统及超声设备 | |
| CN112168211B (zh) | 一种腹部超声图像的脂肪厚度和肌肉厚度测量方法及系统 | |
| US10548564B2 (en) | System and method for ultrasound imaging of regions containing bone structure | |
| CN105266849B (zh) | 实时超声弹性成像方法和系统 | |
| CN104143047A (zh) | 血管内超声灰阶图像的自动组织标定方法 | |
| CN109919953A (zh) | 颈动脉内中膜厚度测量的方法、系统和设备 | |
| CN107970042A (zh) | 一种超声非线性定量实时成像方法与系统 | |
| US11484286B2 (en) | Ultrasound evaluation of anatomical features | |
| Zhang et al. | Boundary delineation in transrectal ultrasound image for prostate cancer | |
| CN110930394B (zh) | 测量肌肉肌纤维束线斜率和羽状角的方法及终端设备 | |
| CN114708283A (zh) | 图像目标的分割方法、装置、电子设备及存储介质 | |
| CN113229850A (zh) | 超声盆底成像方法和超声成像系统 | |
| US20240119705A1 (en) | Systems, methods, and apparatuses for identifying inhomogeneous liver fat | |
| Mohamed et al. | Advancing Cardiac Image Processing: An Innovative Model Utilizing Canny Edge Detection For Enhanced Diagnostics | |
| CN116777879A (zh) | 一种基于YOLOv5的前列腺超声图像辅助诊断方法及系统 | |
| CN117408912A (zh) | 一种超声图像增强和去噪的方法及设备 | |
| CN116650006A (zh) | 用于自动超声检查的系统和方法 | |
| CN115708694A (zh) | 超声图像处理方法及设备 |
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 |