KR102436035B1 - Apparatus and method for diagnosing skeletal muscle condition from ultrasound images using deep neural network - Google Patents

Apparatus and method for diagnosing skeletal muscle condition from ultrasound images using deep neural network Download PDF

Info

Publication number
KR102436035B1
KR102436035B1 KR1020200079181A KR20200079181A KR102436035B1 KR 102436035 B1 KR102436035 B1 KR 102436035B1 KR 1020200079181 A KR1020200079181 A KR 1020200079181A KR 20200079181 A KR20200079181 A KR 20200079181A KR 102436035 B1 KR102436035 B1 KR 102436035B1
Authority
KR
South Korea
Prior art keywords
layer
skeletal muscle
output
unit
ultrasound image
Prior art date
Application number
KR1020200079181A
Other languages
Korean (ko)
Other versions
KR20220001177A (en
Inventor
이영한
신이랑
양재문
이지숙
Original Assignee
연세대학교 산학협력단
인제대학교 산학협력단
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 연세대학교 산학협력단, 인제대학교 산학협력단 filed Critical 연세대학교 산학협력단
Priority to KR1020200079181A priority Critical patent/KR102436035B1/en
Publication of KR20220001177A publication Critical patent/KR20220001177A/en
Application granted granted Critical
Publication of KR102436035B1 publication Critical patent/KR102436035B1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0858Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/042Knowledge-based neural networks; Logical representations of neural networks
    • G06N3/0427
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/40ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Theoretical Computer Science (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Primary Health Care (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Epidemiology (AREA)
  • Computational Linguistics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Quality & Reliability (AREA)
  • Physiology (AREA)
  • Databases & Information Systems (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

본 발명의 골격근의 상태를 진단하기 위한 장치는 골격근에 대한 초음파 영상이 입력되면, 소정의 형식에 따라 상기 초음파 영상을 가공하는 전처리부와, 입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하는 예측망을 통해 상기 가공된 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 예측하는 예측부를 포함한다. The apparatus for diagnosing the state of skeletal muscle according to the present invention includes a preprocessing unit that processes the ultrasound image according to a predetermined format when an ultrasound image of the skeletal muscle is input, an input layer, at least one pair of alternately repeated convolutional layers; By performing pixel unit operation on the processed ultrasound image through a prediction network including a pooling layer, at least one fully connected layer, and an output layer, it is calculated as a probability whether Sarcopenia is in progress in the skeletal muscle, It includes a predictor for predicting whether sarcopenia progresses according to the calculated probability.

Description

심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치 및 이를 위한 방법{Apparatus and method for diagnosing skeletal muscle condition from ultrasound images using deep neural network}Apparatus and method for diagnosing skeletal muscle condition from ultrasound images using deep neural network}

본 발명은 골격근 상태 진단 기술에 관한 것으로, 보다 상세하게는, 심층신경망(DNN: Deep Neural Network)을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치 및 이를 위한 방법에 관한 것이다. The present invention relates to a technique for diagnosing a skeletal muscle state, and more particularly, to an apparatus for diagnosing a skeletal muscle state from an ultrasound image using a deep neural network (DNN) and a method therefor.

근육감소증(Sarcopenia)은 Rosenberg에 의해 1988년에 도입되었으며, 원래의 정의는 노인들의 맹장 근육 질량에 대한 근육 손실이었다. 2010년에 근육감소증의 정의는 근육량이 낮을 뿐만 아니라 근육 기능 저하되는 것으로 변경되었다. 그리고 근육감소증은 개인의 근육 강도의 차이가 있을 수 있으며 나이가 들어감에 따라 감소한다. 근육감소증을 평가하기 위해 근육량을 측정해야 하며, 이러한 근육량을 측정하기 위한 일반적인 의료 영상 양식은 이중 에너지 X-선 흡수 측정법(DXA: dual-energy X-ray absorptiometry), 컴퓨터 단층 촬영(CT: computed tomography), 자기 공명 영상(MRI: magnetic resonance imaging) 및 초음파(USG: ultrasonography)를 포함한다. Sarcopenia was introduced in 1988 by Rosenberg, and its original definition was muscle loss to cecal muscle mass in the elderly. In 2010, the definition of sarcopenia was changed to be a decrease in muscle function as well as low muscle mass. And sarcopenia may be due to differences in individual muscle strength and decreases with age. To evaluate sarcopenia, muscle mass must be measured, and common medical imaging modalities for measuring muscle mass include dual-energy X-ray absorptiometry (DXA) and computed tomography (CT). ), magnetic resonance imaging (MRI), and ultrasonography (USG).

한국공개특허 제2019-0113089호 2019년 10월 08일 공개 (명칭: 근감소증 분석지원을 위한 인공 신경망 기반의 인체 형태 분석법을 채용하는 영상 처리 장치 및 이를 이용한 영상 처리 방법)Korean Patent Publication No. 2019-0113089 published on October 08, 2019 (Title: Image processing apparatus employing artificial neural network-based human body shape analysis method for sarcopenia analysis support and image processing method using the same)

본 발명의 목적은 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치 및 이를 위한 방법을 제공함에 있다. An object of the present invention is to provide an apparatus for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network and a method therefor.

상술한 바와 같은 목적을 달성하기 위한 본 발명의 바람직한 실시에에 따른 골격근의 상태를 진단하기 위한 장치는 골격근에 대한 초음파 영상이 입력되면, 소정의 형식에 따라 상기 초음파 영상을 가공하는 전처리부와, 입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하는 예측망을 통해 상기 가공된 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 예측하는 예측부를 포함한다. An apparatus for diagnosing the state of skeletal muscle according to a preferred embodiment of the present invention for achieving the above object includes a preprocessing unit that processes the ultrasound image according to a predetermined format when an ultrasound image of the skeletal muscle is input; By performing a pixel unit operation on the processed ultrasound image through a prediction network including an input layer, at least a pair of alternately repeated convolutional and pooling layers, at least one fully connected layer, and an output layer, sarcopenia in the skeletal muscle (Sarcopenia) includes a predictor for calculating whether or not progressing with a probability, and predicting whether sarcopenia progresses according to the calculated probability.

상기 장치는 상기 예측부가 근육감소증을 예측하면, 상기 예측망의 상기 컨벌루션층의 출력값 및 경사도(gradient)를 비교하여 상기 초음파 영상에서 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시하는 시각화 영상을 생성하는 시각화부를 더 포함한다. When the prediction unit predicts sarcopenia, the apparatus compares the output value and the gradient of the convolutional layer of the prediction network to distinguish and display the region causing sarcopenia in the ultrasound image from other regions. It further includes a visualization unit for generating.

상기 시각화부는 상기 예측망의 연산의 결과를 수집한 후, 수학식

Figure 112020066801883-pat00001
에 따라 뉴런 중요도 가중치를 산출하고, 수학식
Figure 112020066801883-pat00002
에 따라 상기 시각화 영상을 생성하며, 상기
Figure 112020066801883-pat00003
는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고, 상기
Figure 112020066801883-pat00004
는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고, 상기
Figure 112020066801883-pat00005
Figure 112020066801883-pat00006
Figure 112020066801883-pat00007
에 대해 가지는 경사도(gradient)이고, 상기
Figure 112020066801883-pat00008
는 상기 뉴런 중요도 가중치이고, 상기
Figure 112020066801883-pat00009
는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 한다. After the visualization unit collects the results of the calculation of the prediction network,
Figure 112020066801883-pat00001
Calculate the neuron importance weight according to the formula
Figure 112020066801883-pat00002
to generate the visualization image according to the
Figure 112020066801883-pat00003
is the c-th value of the output before performing the convolution operation in the convolutional layer and performing the operation by the activation function,
Figure 112020066801883-pat00004
is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer,
Figure 112020066801883-pat00005
Is
Figure 112020066801883-pat00006
go
Figure 112020066801883-pat00007
is the gradient with respect to , and
Figure 112020066801883-pat00008
is the neuron importance weight, and
Figure 112020066801883-pat00009
is the value of the coordinates (i, j) of the visualization image.

상기 장치는 상기 전처리부가 상기 골격근에 대한 초음파 영상으로부터 복수의 방사성 특징을 추출하면, 입력층, 은닉층 및 출력층을 포함하는 식별망을 통해 상기 복수의 방사성 특징에 대한 연산을 수행하여 상기 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 식별하는 식별부를 더 포함한다. When the preprocessor extracts a plurality of radioactive features from the ultrasound image of the skeletal muscle, the apparatus performs an operation on the plurality of radioactive features through an identification network including an input layer, a hidden layer, and an output layer to cause sarcopenia in the skeletal muscle (Sarcopenia) further includes an identification unit for calculating whether or not progressing with a probability, and identifying whether sarcopenia is progressing according to the calculated probability.

상기 장치는 상기 전처리부가 상기 골격근에 대한 초음파 영상을 상기 소정의 형식과 다른 형식으로 상기 초음파 영상을 가공하면, 입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하는 분류망을 통해 상기 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근이 근육에코발생도(Muscle echogenicity)의 각 등급에 속할 확률을 산출하고, 산출된 확률에 따라 상기 골격근의 근육에코발생도의 등급을 분류하는 분류부를 더 포함한다. In the device, when the preprocessor processes the ultrasound image of the skeletal muscle in a format different from the predetermined format, an input layer, at least a pair of convolutional layers and a pooling layer that are repeated alternately, and at least one complete connection A pixel unit operation is performed on the ultrasound image through a classification network including a layer and an output layer to calculate the probability that the skeletal muscle belongs to each class of muscle echogenicity, and according to the calculated probability, It further includes a classification unit for classifying the grade of the muscle echogenicity.

상기 예측망의 상기 출력층이 소프트맥스 연산을 수행하는 출력노드로 이루어진 상기 예측망의 초기 모델을 마련하고, 상기 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행하며, 상기 초기 학습이 완료되면, 상기 예측망의 마지막 완전연결층을 초기화된 완전연결층으로 교체하고, 상기 출력층의 상기 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 상기 예측망의 정규 모델을 마련하고, 상기 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행하는 학습부를 더 포함한다. The output layer of the prediction network prepares an initial model of the prediction network including an output node that performs a softmax operation, and performs initial learning using an image stored in a public database for the initial model as training data, When the initial learning is completed, the last fully connected layer of the prediction network is replaced with an initialized fully connected layer, and the output node performing the softmax operation of the output layer is replaced with an output node performing the sigmoid operation. It further includes a learning unit that prepares a regular model of the prediction network, and performs regular learning on the regular model by using an ultrasound image as training data.

상기 학습부는 상기 예측망에 대해 상기 정규 학습을 수행할 때, 손실함수

Figure 112022026989181-pat00010
를 이용하여 상기
Figure 112022026989181-pat00011
을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하고, 상기
Figure 112022026989181-pat00012
을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행하고, 상기 E(y, f)는 손실함수이고, 상기 yi는 출력층의 출력값이고, 상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고, 상기 gij는 은닉층의 마지막층의 출력값이고, 상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고, 상기 i는 출력층의 노드에 대응하는 인덱스이고, 상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고, 상기
Figure 112022026989181-pat00013
는 하이퍼파라미터인 것을 특징으로 한다. When the learning unit performs the regular learning on the prediction network, a loss function
Figure 112022026989181-pat00010
above using
Figure 112022026989181-pat00011
is set to 0, and classification loss optimization is performed to correct the parameters of the prediction network so that the classification loss, which is the difference between the output value of the prediction network and the classification label, is minimized, and
Figure 112022026989181-pat00012
is set to 0.5, and the composite loss including the classification loss, which is the difference between the output value of the prediction network and the classification label, and the concealment loss, which is the difference between the output value and the hidden label of the last layer of the hidden layer of the prediction network, is predicted to be minimal. Perform complex loss optimization by modifying network parameters, where E(y, f) is a loss function, yi is an output value of an output layer, ei is a classification label corresponding to an output value of an output layer, and gij is a hidden layer is the output value of the last layer of , cij is the hidden label corresponding to the output value of the last layer of the hidden layer, i is the index corresponding to the node of the output layer, and j is the index corresponding to the node of the last layer of the hidden layer, remind
Figure 112022026989181-pat00013
is a hyperparameter.

상술한 바와 같은 목적을 달성하기 위한 본 발명의 바람직한 실시예에 따른 골격근의 상태를 진단하기 위한 장치는 골격근에 대한 초음파 영상이 입력되면, 서로 다른 2개의 형식에 따라 상기 초음파 영상을 가공하고, 상기 초음파 영상으로부터 복수의 방사성 특징을 추출하는 전처리부와, 복수의 계층을 포함하는 분류망을 통해 상기 2개의 형식 중 어느 하나의 형식으로 가공된 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근의 근육에코발생도의 등급을 확률로 산출하고, 산출된 확률에 따라 상기 골격근의 근육에코발생도의 등급을 분류하는 분류부와, 복수의 계층을 포함하는 예측망을 통해 상기 2개의 형식 중 다른 하나의 형식으로 가공된 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 예측하는 예측부와, 복수의 계층을 포함하는 식별망을 통해 상기 복수의 방사성 특징에 대한 연산을 수행하여 상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 식별하는 식별부를 포함한다. In an apparatus for diagnosing a state of skeletal muscle according to a preferred embodiment of the present invention for achieving the above object, when an ultrasound image of skeletal muscle is input, the ultrasound image is processed according to two different formats, and the A preprocessor for extracting a plurality of radioactive features from an ultrasound image, and a classification network including a plurality of layers perform a pixel unit operation on an ultrasound image processed in one of the two formats through a classification network including a plurality of layers A classification unit that calculates a grade of echogenicity as a probability, classifies a grade of a muscle echogenicity of the skeletal muscle according to the calculated probability, and a prediction network including a plurality of layers. A prediction unit that calculates whether sarcopenia is progressing in the skeletal muscle by performing pixel unit operation on the ultrasound image processed in the form with a probability, and predicts whether sarcopenia is progressing according to the calculated probability; and a plurality of layers. It includes an identification unit that calculates whether sarcopenia is in progress in the skeletal muscle by performing an operation on the plurality of radioactive features through an identification network including a probability, and identifies whether sarcopenia is in progress according to the calculated probability.

상기 예측망은 입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하며, 상기 예측부가 근육감소증을 예측하면, 상기 예측망의 상기 컨벌루션층의 출력값 및 경사도(gradient)를 비교하여 상기 초음파 영상에서 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시하는 시각화 영상을 생성하는 시각화부를 더 포함한다. The prediction network includes an input layer, at least a pair of alternately repeated convolutional and pooling layers, at least one fully connected layer and an output layer, and when the prediction unit predicts sarcopenia, the The method further includes a visualization unit that compares an output value and a gradient to generate a visualization image that distinguishes and displays a region causing sarcopenia in the ultrasound image from other regions.

상기 시각화부는 상기 예측망의 연산의 결과를 수집한 후, 수학식

Figure 112020066801883-pat00014
에 따라 뉴런 중요도 가중치를 산출하고, 수학식
Figure 112020066801883-pat00015
에 따라 상기 시각화 영상을 생성하며, 상기
Figure 112020066801883-pat00016
는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고, 상기
Figure 112020066801883-pat00017
는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고, 상기
Figure 112020066801883-pat00018
Figure 112020066801883-pat00019
Figure 112020066801883-pat00020
에 대해 가지는 경사도(gradient)이고, 상기
Figure 112020066801883-pat00021
는 상기 뉴런 중요도 가중치이고, 상기
Figure 112020066801883-pat00022
는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 한다. After the visualization unit collects the results of the calculation of the prediction network,
Figure 112020066801883-pat00014
Calculate the neuron importance weight according to the formula
Figure 112020066801883-pat00015
to generate the visualization image according to the
Figure 112020066801883-pat00016
is the c-th value of the output before performing the convolution operation in the convolutional layer and performing the operation by the activation function,
Figure 112020066801883-pat00017
is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer,
Figure 112020066801883-pat00018
Is
Figure 112020066801883-pat00019
go
Figure 112020066801883-pat00020
is the gradient with respect to , and
Figure 112020066801883-pat00021
is the neuron importance weight, and
Figure 112020066801883-pat00022
is the value of the coordinates (i, j) of the visualization image.

상기 장치는 상기 예측망의 출력층의 출력노드가 소프트맥스 연산을 수행하는 상기 예측망의 초기 모델을 마련하고, 상기 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행하며, 상기 초기 학습이 완료되면, 상기 분류망 및 상기 예측망의 은닉층의 마지막층을 초기화된 완전연결층으로 교체하고, 상기 출력층의 상기 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 상기 예측망의 정규 모델을 마련하고, 상기 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행하는 학습부를 더 포함한다. The apparatus prepares an initial model of the prediction network in which the output node of the output layer of the prediction network performs a softmax operation, and performs initial learning using an image stored in a public database for the initial model as training data, When the initial learning is completed, the last layer of the hidden layer of the classification network and the prediction network is replaced with an initialized fully connected layer, and an output node that performs the softmax operation of the output layer performs a sigmoid operation. It further includes a learning unit that prepares a regular model of the prediction network by replacing it with a node, and performs regular learning on the regular model by using an ultrasound image as training data.

상기 학습부는 상기 예측망에 대해 상기 정규 학습을 수행할 때, 손실함수

Figure 112022026989181-pat00023
를 이용하여 상기
Figure 112022026989181-pat00024
을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하고, 상기
Figure 112022026989181-pat00025
을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행한다. When the learning unit performs the regular learning on the prediction network, a loss function
Figure 112022026989181-pat00023
above using
Figure 112022026989181-pat00024
is set to 0, and classification loss optimization is performed to correct the parameters of the prediction network so that the classification loss, which is the difference between the output value of the prediction network and the classification label, is minimized, and
Figure 112022026989181-pat00025
is set to 0.5, and the composite loss including the classification loss, which is the difference between the output value of the prediction network and the classification label, and the concealment loss, which is the difference between the output value and the hidden label of the last layer of the hidden layer of the prediction network, is predicted to be minimal. Perform complex loss optimization that modifies network parameters.

여기서, 상기 E(y, f)는 손실함수이고, 상기 yi는 출력층의 출력값이고, 상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고, 상기 gij는 은닉층의 마지막층의 출력값이고, 상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고, 상기 i는 출력층의 노드에 대응하는 인덱스이고, 상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고, 상기

Figure 112020066801883-pat00026
는 하이퍼파라미터인 것을 특징으로 한다. where E(y, f) is a loss function, yi is an output value of the output layer, ei is a classification label corresponding to an output value of the output layer, gij is an output value of the last layer of the hidden layer, and cij is the output value of the hidden layer is a hidden label corresponding to the output value of the last layer of , wherein i is an index corresponding to a node of the output layer, and j is an index corresponding to a node of the last layer of the hidden layer, wherein
Figure 112020066801883-pat00026
is a hyperparameter.

상술한 바와 같은 목적을 달성하기 위한 본 발명의 바람직한 실시예에 따른 골격근의 상태를 진단하기 위한 방법은 예측부가 골격근에 대한 초음파 영상이 입력되면, 입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하는 예측망을 통해 상기 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하는 단계와, 상기 예측부가 산출된 확률에 따라 근육감소증 진행 여부를 예측하는 단계를 포함한다. In the method for diagnosing the state of skeletal muscle according to a preferred embodiment of the present invention for achieving the above object, when the prediction unit receives an ultrasound image of skeletal muscle, an input layer and at least one pair of convolutional layers alternately repeated and performing a pixel unit operation on the ultrasound image through a prediction network including a pooling layer, at least one fully connected layer, and an output layer to calculate with a probability whether sarcopenia is in progress in the skeletal muscle; and predicting whether sarcopenia progresses according to the additionally calculated probability.

상기 방법은 상기 예측부가 근육감소증을 예측하면, 시각화부가 상기 예측망의 상기 컨벌루션층의 출력값 및 경사도(gradient)를 비교하여 상기 초음파 영상에서 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시하는 시각화 영상을 생성하는 단계를 더 포함한다. In the method, when the prediction unit predicts sarcopenia, the visualization unit compares the output value and the gradient of the convolutional layer of the prediction network to distinguish and display the region causing sarcopenia in the ultrasound image from other regions The method further includes generating a visualization image.

상기 시각화 영상을 생성하는 단계는 상기 시각화부가 상기 예측망의 연산의 결과를 수집하는 단계와, 상기 시각화부가 수학식

Figure 112020066801883-pat00027
에 따라 뉴런 중요도 가중치를 산출하는 단계와, 상기 시각화부가 수학식
Figure 112020066801883-pat00028
에 따라 상기 시각화 영상을 생성하는 단계를 포함한다. 여기서, 상기
Figure 112020066801883-pat00029
는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고, 상기
Figure 112020066801883-pat00030
는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고, 상기
Figure 112020066801883-pat00031
는 상기
Figure 112020066801883-pat00032
가 상기
Figure 112020066801883-pat00033
에 대해 가지는 경사도(gradient)이고, 상기
Figure 112020066801883-pat00034
는 상기 뉴런 중요도 가중치이고, 상기
Figure 112020066801883-pat00035
는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 한다. The generating of the visualization image includes: the visualization unit collecting the results of the calculation of the prediction network;
Figure 112020066801883-pat00027
calculating the neuron importance weight according to
Figure 112020066801883-pat00028
and generating the visualization image according to the Here, the
Figure 112020066801883-pat00029
is the c-th value of the output before performing the convolution operation in the convolutional layer and performing the operation by the activation function,
Figure 112020066801883-pat00030
is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer,
Figure 112020066801883-pat00031
is said
Figure 112020066801883-pat00032
is reminded
Figure 112020066801883-pat00033
is the gradient with respect to , and
Figure 112020066801883-pat00034
is the neuron importance weight, and
Figure 112020066801883-pat00035
is the value of the coordinates (i, j) of the visualization image.

상기 방법은 식별부가 상기 골격근에 대한 초음파 영상으로부터 추출된 복수의 방사성 특징이 입력되면, 입력층, 은닉층 및 출력층을 포함하는 식별망을 통해 상기 복수의 방사성 특징에 대한 연산을 수행하여 상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하는 단계와, 상기 식별부가 산출된 확률에 따라 근육감소증 진행 여부를 식별하는 단계를 더 포함한다. In the method, when a plurality of radioactive features extracted from the ultrasound image of the skeletal muscle are input by the identification unit, an operation is performed on the plurality of radioactive features through an identification network including an input layer, a hidden layer, and an output layer. The method further includes: calculating whether or not sarcopenia is progressing with a probability; and identifying whether sarcopenia is progressing according to the probability calculated by the identification unit.

상기 방법은 분류부가 상기 골격근에 대한 초음파 영상이 입력되면, 입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하는 분류망을 통해 상기 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근이 근육에코발생도의 각 등급에 속할 확률을 산출하는 단계와, 상기 분류부가 상기 산출된 확률에 따라 상기 골격근의 근육에코발생도의 등급을 분류하는 단계를 더 포함한다. In the method, when the classification unit receives the ultrasound image of the skeletal muscle, the ultrasound image is passed through a classification network including an input layer, at least a pair of alternately repeated convolutional and pooling layers, at least one fully connected layer, and an output layer. calculating the probability that the skeletal muscle belongs to each grade of the muscle echogenicity by performing pixel unit operation on include more

상기 방법은 상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하는 단계 전, 학습부가 상기 예측망의 상기 출력층이 소프트맥스 연산을 수행하는 출력노드로 이루어진 상기 예측망의 초기 모델을 마련하는 단계와, 상기 학습부가 상기 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행하는 단계와, 상기 학습부가 상기 예측망의 마지막 완전연결층을 초기화된 완전연결층으로 교체하고, 상기 출력층의 상기 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 상기 예측망의 정규 모델을 마련하는 단계와, 상기 학습부가 상기 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행하는 단계를 더 포함한다. In the method, before calculating with probability whether sarcopenia is in progress in the skeletal muscle, a learning unit provides an initial model of the prediction network including an output node in which the output layer of the prediction network performs a softmax operation and performing initial learning by the learning unit using an image stored in a public database for the initial model as training data, and the learning unit replacing the last fully connected layer of the prediction network with an initialized fully connected layer, preparing a regular model of the prediction network by replacing an output node that performs the softmax operation of the output layer with an output node that performs a sigmoid operation; It further includes the step of performing regular learning using.

상기 정규 학습을 수행하는 단계는 상기 학습부가 손실함수

Figure 112022026989181-pat00036
를 이용하여 상기
Figure 112022026989181-pat00037
을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하는 단계와, 상기 학습부가 상기
Figure 112022026989181-pat00038
을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행하는 단계를 포함한다. 여기서, 상기 E(y, f)는 손실함수이고, 상기 yi는 출력층의 출력값이고, 상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고, 상기 gij는 은닉층의 마지막층의 출력값이고, 상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고, 상기 i는 출력층의 노드에 대응하는 인덱스이고, 상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고, 상기
Figure 112022026989181-pat00039
는 하이퍼파라미터인 것을 특징으로 한다. In the step of performing the regular learning, the learning unit is a loss function
Figure 112022026989181-pat00036
above using
Figure 112022026989181-pat00037
is set to 0, and then performing classification loss optimization by modifying the parameters of the prediction network so that the classification loss, which is the difference between the output value of the prediction network and the classification label, is minimized;
Figure 112022026989181-pat00038
is set to 0.5, and the composite loss including the classification loss, which is the difference between the output value of the prediction network and the classification label, and the concealment loss, which is the difference between the output value and the hidden label of the last layer of the hidden layer of the prediction network, is predicted to be minimal. and performing complex loss optimization to modify parameters of the network. where E(y, f) is a loss function, yi is an output value of the output layer, ei is a classification label corresponding to an output value of the output layer, gij is an output value of the last layer of the hidden layer, and cij is the output value of the hidden layer is a hidden label corresponding to the output value of the last layer of , wherein i is an index corresponding to a node of the output layer, and j is an index corresponding to a node of the last layer of the hidden layer, wherein
Figure 112022026989181-pat00039
is a hyperparameter.

본 발명에 따르면 복수의 심층신경망을 이용하여 초음파 영상으로 촬상된 근육의 상태에 대해 다양한 형태의 정보를 도출할 수 있다. 이는 2 종의 서로 다른 형식의 초음파 영상(B-mode USG 및 SWE(shear-wave elastography) USG) 및 초음파 영상(B-mode USG)로부터 추출되는 방사성 특징(B-mode USG radiomics features)을 분석하여 얻어지는 근육감소증(Sarcopenia)에 대한 확률과 같이 진단에 직관적인 정보뿐만 아니라 근육에코발생도(Muscle echogenicity) 및 시각화 영상과 같은 간접적인 정보를 포함한다. 이에 따라, 다각적이고 종합적으로 골격근의 상태를 진단 할 수 있다. According to the present invention, it is possible to derive various types of information about the state of a muscle captured by an ultrasound image using a plurality of deep neural networks. This is achieved by analyzing the radiomics features extracted from two different types of ultrasound images (B-mode USG and shear-wave elastography (SWE) USG) and ultrasound images (B-mode USG). It includes indirect information such as muscle echogenicity and visualization images as well as intuitive information for diagnosis, such as the probability of acquired sarcopenia. Accordingly, it is possible to diagnose the condition of skeletal muscle in a multifaceted and comprehensive manner.

도 1은 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 시스템의 구성을 설명하기 위한 도면이다.
도 2는 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 세부적인 구성을 설명하기 위한 도면이다.
도 3은 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 전처리부를 설명하기 위한 도면이다.
도 4 내지 도 6은 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 분류부를 설명하기 위한 도면이다.
도 7 및 도 8은 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 예측부를 설명하기 위한 도면이다.
도 9는 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 식별부를 설명하기 위한 도면이다.
도 10은 본 발명의 실시예에 따른 골격근의 상태를 진단하도록 식별망을 학습시키는 방법을 설명하기 위한 흐름도이다.
도 11은 본 발명의 실시예에 따른 골격근의 상태를 진단하도록 분류망 및 예측망을 학습시키는 방법을 설명하기 위한 흐름도이다.
도 12는 본 발명의 실시예에 따른 골격근의 상태를 진단하도록 분류망 및 예측망에 대한 정규 학습을 수행하는 방법을 설명하기 위한 흐름도이다.
도 13은 본 발명의 실시예에 따른 골격근의 상태를 진단하도록 분류망 및 예측망에 대한 정규 학습을 수행하는 방법을 설명하기 위한 벡터공간의 예이다.
도 14는 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 방법을 설명하기 위한 흐름도이다.
1 is a diagram for explaining the configuration of a system for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention.
2 is a view for explaining the detailed configuration of an apparatus for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention.
3 is a view for explaining a preprocessor of an apparatus for diagnosing a state of a skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention.
4 to 6 are diagrams for explaining a classification unit of an apparatus for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention.
7 and 8 are diagrams for explaining a prediction unit of an apparatus for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention.
9 is a view for explaining an identification unit of an apparatus for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention.
10 is a flowchart illustrating a method of training an identification network to diagnose a state of a skeletal muscle according to an embodiment of the present invention.
11 is a flowchart illustrating a method for training a classification network and a prediction network to diagnose a state of a skeletal muscle according to an embodiment of the present invention.
12 is a flowchart illustrating a method of performing regular learning on a classification network and a prediction network to diagnose a state of a skeletal muscle according to an embodiment of the present invention.
13 is an example of a vector space for explaining a method of performing regular learning for a classification network and a prediction network to diagnose a state of a skeletal muscle according to an embodiment of the present invention.
14 is a flowchart illustrating a method for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention.

본 발명의 상세한 설명에 앞서, 이하에서 설명되는 본 명세서 및 청구범위에 사용된 용어나 단어는 통상적이거나 사전적인 의미로 한정해서 해석되어서는 아니 되며, 발명자는 그 자신의 발명을 가장 최선의 방법으로 설명하기 위해 용어의 개념으로 적절하게 정의할 수 있다는 원칙에 입각하여 본 발명의 기술적 사상에 부합하는 의미와 개념으로 해석되어야만 한다. 따라서 본 명세서에 기재된 실시예와 도면에 도시된 구성은 본 발명의 가장 바람직한 실시예에 불과할 뿐, 본 발명의 기술적 사상을 모두 대변하는 것은 아니므로, 본 출원시점에 있어서 이들을 대체할 수 있는 다양한 균등물과 변형 예들이 있을 수 있음을 이해하여야 한다. Prior to the detailed description of the present invention, the terms or words used in the present specification and claims described below should not be construed as being limited to their ordinary or dictionary meanings, and the inventors should develop their own inventions in the best way. It should be interpreted as meaning and concept consistent with the technical idea of the present invention based on the principle that it can be appropriately defined as a concept of a term for explanation. Accordingly, the embodiments described in this specification and the configurations shown in the drawings are only the most preferred embodiments of the present invention, and do not represent all the technical ideas of the present invention, so various equivalents that can replace them at the time of the present application It should be understood that there may be water and variations.

이하, 첨부된 도면을 참조하여 본 발명의 바람직한 실시예들을 상세히 설명한다. 이때, 첨부된 도면에서 동일한 구성 요소는 가능한 동일한 부호로 나타내고 있음을 유의해야 한다. 또한, 본 발명의 요지를 흐리게 할 수 있는 공지 기능 및 구성에 대한 상세한 설명은 생략할 것이다. 마찬가지의 이유로 첨부 도면에 있어서 일부 구성요소는 과장되거나 생략되거나 또는 개략적으로 도시되었으며, 각 구성요소의 크기는 실제 크기를 전적으로 반영하는 것이 아니다. Hereinafter, preferred embodiments of the present invention will be described in detail with reference to the accompanying drawings. In this case, it should be noted that the same components in the accompanying drawings are denoted by the same reference numerals as much as possible. In addition, detailed descriptions of well-known functions and configurations that may obscure the gist of the present invention will be omitted. For the same reason, some components are exaggerated, omitted, or schematically illustrated in the accompanying drawings, and the size of each component does not fully reflect the actual size.

먼저, 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 시스템에 대해서 설명하기로 한다. 도 1은 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 시스템의 구성을 설명하기 위한 도면이다. 도 1을 참조하면, 본 발명의 실시예에 따른 진단 시스템은 진단장치(10) 및 진단영상장치(20)를 포함한다. First, a system for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention will be described. 1 is a diagram for explaining the configuration of a system for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention. Referring to FIG. 1 , a diagnostic system according to an embodiment of the present invention includes a diagnostic apparatus 10 and a diagnostic imaging apparatus 20 .

진단영상장치(20)는 기본적으로, 초음파(USG: ultrasonography) 영상을 촬영하기 위한 것이다. 본 발명의 실시예에서 진단영상장치(20)는 B 모드 초음파 영상 및 SWE 초음파 영상 양자 모두를 촬영할 수 있는 하나의 장치로 설명될 것이지만, 진단영상장치(20)는 B 모드 초음파 영상을 촬영하는 장치와 SWE 초음파 영상을 촬영하는 장치로 구분되어 2개의 장치로 구현될 수 있다. The diagnostic imaging apparatus 20 is basically for capturing an ultrasonography (USG) image. In the embodiment of the present invention, the diagnostic imaging apparatus 20 will be described as one device capable of photographing both a B-mode ultrasound image and an SWE ultrasound image, but the diagnostic imaging apparatus 20 is a device for photographing a B-mode ultrasound image. It is divided into a device for capturing an SWE ultrasound image and can be implemented as two devices.

진단장치(10)는 진단영상장치(20)로부터 초음파 영상을 입력 받고, 심층신경망, 예컨대, 심층컨벌루션신경망(DCNN: deep convolutional neural network), 다층퍼셉트론(MLP: Multi-Layer Perceptron)을 이용하여 2 종의 서로 다른 형식의 초음파 영상(B-mode USG 및 SWE USG) 및 초음파 영상(B-mode USG)로부터 추출되는 방사성 특징(B-mode USG radiomics features)을 분석하여 근육의 상태에 대해 다양한 형태의 정보를 도출하여 이를 제공하기 위한 것이다. 이러한 진단장치(10)는 인터페이스부(11), 입력부(12), 표시부(13), 저장부(14) 및 제어부(15)를 포함한다. The diagnostic apparatus 10 receives an ultrasound image from the diagnostic imaging apparatus 20, and uses a deep neural network, for example, a deep convolutional neural network (DCNN), a multi-layer perceptron (MLP). By analyzing the different types of ultrasound images (B-mode USG and SWE USG) and from the ultrasound images (B-mode USG) of the species (B-mode USG radiomics features), various types of It is intended to derive information and provide it. The diagnostic apparatus 10 includes an interface unit 11 , an input unit 12 , a display unit 13 , a storage unit 14 , and a control unit 15 .

인터페이스부(11)는 진단영상장치(20)로부터 초음파 영상을 수신하기 위한 것이다. 인터페이스부(11)는 필요에 따라 무선 통신, 유선 통신 등을 통해 진단영상장치(20)와 연결되어 초음파 영상을 수신할 수 있다. The interface unit 11 is for receiving an ultrasound image from the diagnostic imaging apparatus 20 . The interface unit 11 may be connected to the diagnostic imaging apparatus 20 through wireless communication, wired communication, etc. as necessary to receive an ultrasound image.

입력부(12)는 진단장치(10)의 각 종 기능, 동작 등을 제어하기 위한 사용자의 키 조작을 입력받고 입력 신호를 생성하여 제어부(15)에 전달한다. 입력부(12)는 키보드, 마우스 등을 예시할 수 있다. 입력부(12)는 전원 on/off를 위한 전원키, 문자키, 숫자키, 방향키 등을 포함할 수 있다. 입력부(12)의 기능은 표시부(13)가 터치스크린으로 구현된 경우, 표시부(13)에서 이루어질 수 있으며, 표시부(13)만으로 모든 기능을 수행할 수 있는 경우, 입력부(12)는 생략될 수도 있다. The input unit 12 receives a user's key manipulation for controlling various functions and operations of the diagnostic apparatus 10 , generates an input signal, and transmits the generated input signal to the control unit 15 . The input unit 12 may be a keyboard, a mouse, or the like. The input unit 12 may include a power key for power on/off, a character key, a number key, a direction key, and the like. The function of the input unit 12 may be performed on the display unit 13 when the display unit 13 is implemented as a touch screen, and when all functions can be performed only with the display unit 13, the input unit 12 may be omitted. have.

표시부(13)는 진단장치(10)의 메뉴, 입력된 데이터, 기능 설정 정보 및 기타 다양한 정보를 사용자에게 시각적으로 제공한다. 특히, 본 발명의 실시예에 따르면, 제어부(15)는 입력된 초음파 영상에 대한 분석을 수행하여 그 분석 결과를 제공한다. 이러한 경우, 표시부(13)는 제어부(15)의 제어에 따라 분석 결과를 표시부(13)를 통해 표시될 수 있다. 표시부(13)는 진단장치(10)의 부팅 화면, 대기 화면, 메뉴 화면, 등의 화면을 출력하는 기능을 수행한다. 이러한 표시부(13)는 액정표시장치(LCD, Liquid Crystal Display), 유기 발광 다이오드(OLED, Organic Light Emitting Diodes), 능동형 유기 발광 다이오드(AMOLED, Active Matrix Organic Light Emitting Diodes) 등으로 형성될 수 있다. The display unit 13 visually provides a menu of the diagnostic apparatus 10 , input data, function setting information, and other various information to the user. In particular, according to an embodiment of the present invention, the controller 15 analyzes the input ultrasound image and provides the analysis result. In this case, the display unit 13 may display the analysis result through the display unit 13 under the control of the control unit 15 . The display unit 13 performs a function of outputting a boot screen, a standby screen, a menu screen, and the like of the diagnostic apparatus 10 . The display unit 13 may be formed of a liquid crystal display (LCD), an organic light emitting diode (OLED), an active matrix organic light emitting diode (AMOLED), or the like.

저장부(14)는 진단장치(10)의 동작에 필요한 각 종 데이터, 애플리케이션, 진단장치(10)의 동작에 따라 구성되는 각 종 데이터를 저장하는 역할을 수행한다. 이러한 저장부(14)는 크게 프로그램 영역과 데이터 영역을 포함할 수 있다. 프로그램 영역은 진단장치(10)의 부팅(booting) 및 운영(operation)을 위한 운영체제(OS, Operating System), 본 발명의 실시예에 따른 애플리케이션 등을 저장할 수 있다. 데이터 영역은 본 발명의 실시예에 따른 진단을 위한 각 종 데이터를 저장할 수 있다. 저장부(14)에 저장되는 각 종 데이터는 사용자의 조작에 따라, 삭제, 변경, 추가될 수 있다. The storage unit 14 serves to store various types of data necessary for the operation of the diagnostic apparatus 10 , applications, and various types of data configured according to the operation of the diagnostic apparatus 10 . The storage unit 14 may largely include a program area and a data area. The program area may store an operating system (OS) for booting and operation of the diagnostic apparatus 10 , an application according to an embodiment of the present invention, and the like. The data area may store various types of data for diagnosis according to an embodiment of the present invention. Various types of data stored in the storage unit 14 may be deleted, changed, or added according to a user's operation.

제어부(15)는 진단장치(10)의 전반적인 동작 및 진단장치(10)의 내부 블록들 간 신호 흐름을 제어하고, 데이터를 처리하는 데이터 처리 기능을 수행할 수 있다. 이러한 제어부(15)는 중앙처리장치(CPU: Central Processing Unit), 디지털신호처리기(DSP: Digital Signal Processor) 등이 될 수 있다. 또한, 제어부(15)는 추가로 이미지 프로세서(Image processor) 혹은 GPU(Graphic Processing Unit)를 더 구비할 수 있다. 이러한 제어부(15)의 구체적인 동작에 대해서는 아래에서 더 상세하게 설명될 것이다. The controller 15 may control the overall operation of the diagnostic apparatus 10 and signal flow between internal blocks of the diagnostic apparatus 10 , and perform a data processing function of processing data. The control unit 15 may be a central processing unit (CPU), a digital signal processor (DSP), or the like. Also, the controller 15 may further include an image processor or a graphic processing unit (GPU). The specific operation of the control unit 15 will be described in more detail below.

다음으로, 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 구성에 대해 보다 상세하게 설명하기로 한다. 도 2는 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 세부적인 구성을 설명하기 위한 도면이다. 도 3은 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 전처리부를 설명하기 위한 도면이다. 도 4 내지 도 6은 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 분류부를 설명하기 위한 도면이다. 도 7 및 도 8은 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 예측부를 설명하기 위한 도면이다. 그리고 도 9는 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 식별부를 설명하기 위한 도면이다. Next, a configuration of an apparatus for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention will be described in more detail. 2 is a view for explaining the detailed configuration of an apparatus for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention. 3 is a view for explaining a preprocessor of an apparatus for diagnosing a state of a skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention. 4 to 6 are diagrams for explaining a classification unit of an apparatus for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention. 7 and 8 are diagrams for explaining a prediction unit of an apparatus for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention. And FIG. 9 is a view for explaining an identification unit of an apparatus for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention.

먼저, 도 2를 참조하면, 제어부(15)는 학습부(100), 전처리부(200), 분류부(300), 예측부(400), 시각화부(500), 식별부(600) 및 리포트부(700)를 포함한다. First, referring to FIG. 2 , the control unit 15 includes a learning unit 100 , a preprocessing unit 200 , a classification unit 300 , a prediction unit 400 , a visualization unit 500 , an identification unit 600 , and a report. part 700 .

학습부(100)는 본 발명의 실시예에 따른 심층신경망(DNN: Deep Neural Network)을 학습시키기 위한 일련의 프로세스를 수행한다. 한편, 분류부(300), 예측부(400) 및 식별부(600)는 각각 심층신경망으로 분류망(310), 예측망(410) 및 식별망(610)을 포함한다. 학습부(100)는 이러한 분류망(310), 예측망(410) 및 식별망(610)을 학습시킨다. 이러한 학습 프로세스에 대해서는 아래에서 더 상세하게 설명하기로 한다. The learning unit 100 performs a series of processes for learning a deep neural network (DNN) according to an embodiment of the present invention. Meanwhile, the classification unit 300 , the prediction unit 400 , and the identification unit 600 are deep neural networks and include a classification network 310 , a prediction network 410 , and an identification network 610 , respectively. The learning unit 100 trains the classification network 310 , the prediction network 410 , and the identification network 610 . This learning process will be described in more detail below.

도 3을 참조하면, 전처리부(200)는 진단영상장치(20)로부터 인터페이스부(11)를 통해 초음파 영상을 수신하면, 분류부(300)의 분류망(310), 예측부(400)의 예측망(410) 및 식별부(600)의 식별망(610) 중 적어도 하나에 부합하는 형식에 맞춰 수신된 초음파 영상을 가공하거나, 정보를 추출하여 가공된 초음파 영상 혹은 정보를 제공한다. 이를 위하여 전처리부(200)는 증강부(210), 영상처리부(220) 및 특징처리부(230)를 포함한다. Referring to FIG. 3 , when the preprocessor 200 receives an ultrasound image from the diagnostic imaging apparatus 20 through the interface unit 11 , the classification network 310 of the classification unit 300 and the prediction unit 400 are The received ultrasound image is processed according to a format corresponding to at least one of the prediction network 410 and the identification network 610 of the identification unit 600 , or information is extracted to provide the processed ultrasound image or information. To this end, the preprocessor 200 includes an augmentation unit 210 , an image processing unit 220 , and a feature processing unit 230 .

증강부(210)는 본 발명의 실시예에 따른 학습 데이터의 수를 증강시키기 위한 것이다. 일 실시예에 따르면, 증강부(210)는 학습용으로 마련된 초음파 영상이 입력되면, 플립핑(flipping), 시프트(shift), 쉬어(shear), 확대, 축소 및 회전 등의 변형을 통해 학습 데이터의 수를 증강시킬 수 있다. The augmentation unit 210 is for augmenting the number of learning data according to an embodiment of the present invention. According to an embodiment, when an ultrasound image prepared for learning is input, the augmentation unit 210 performs transformations such as flipping, shifting, shearing, enlargement, reduction, and rotation of the learning data. number can be increased.

영상처리부(220)는 진단영상장치(20)로부터 인터페이스부(11)를 통해 초음파 영상을 수신하면, 분류부(300)의 분류망(310) 및 예측부(400)의 예측망(410) 중 적어도 하나에 부합하는 형식에 맞춰 수신된 초음파 영상을 가공하고, 가공된 초음파 영상을 제공한다. 영상처리부(220)는 B 모드 및 SWE 중 어느 하나의 포맷의 초음파 영상이 입력되면, 초음파 영상의 크기를 분류망(310)의 입력층(GIL)의 크기에 맞춰 예컨대, 112x112 크기로 조정하고, DICOM(Digital Imaging and Communications in Medicine) 포맷을 이미지 파일, 예컨대, PNG, JPEG, GIF 등의 포맷으로 변환할 수 있다. 이와 같이, 영상처리부(220)는 초음파 영상에 대해 크기 조정 및 포맷 변환과 같은 가공을 수행한 후, 가공된 초음파 영상을 분류부(300)의 분류망(310)에 제공한다. 또한, 영상처리부(220)는 초음파 영상이 입력되면, B 모드 및 SWE 중 어느 하나의 포맷의 초음파 영상이 입력되면, 입력된 초음파 영상으로부터 소정의 크기, 예컨대, 84x84의 크기로 복수의 영상을 잘라내어 추출하고, 추출된 영상을 DICOM 포맷에서 이미지 파일의 포맷으로 변환할 수 있다. 초음파 영상은 높이와 너비가 서로 다르기 때문에 입력을 만들기 위해 크기를 조정하는 경우, 근육감소증을 분류하는 성능이 상당히 저하된다. 따라서 초음파 영상으로부터 소정의 크기로 잘라낸 복수의 초음파 영상을 예측망(410)의 입력층(EIL)의 크기에 맞춰 재구성한다. 그런 다음, 영상처리부(220)는 전술한 바와 같이 가공된 복수의 초음파 영상을 예측부(400)의 분류망(410)에 제공한다. When the image processing unit 220 receives an ultrasound image from the diagnostic imaging apparatus 20 through the interface unit 11 , one of the classification network 310 of the classification unit 300 and the prediction network 410 of the prediction unit 400 . The received ultrasound image is processed according to a format corresponding to at least one, and the processed ultrasound image is provided. When an ultrasound image in any one format of the B mode and SWE is input, the image processing unit 220 adjusts the size of the ultrasound image to a size of, for example, 112x112 according to the size of the input layer (GIL) of the classification network 310, The DICOM (Digital Imaging and Communications in Medicine) format may be converted into an image file, for example, a format such as PNG, JPEG, or GIF. As described above, the image processing unit 220 provides the processed ultrasound image to the classification network 310 of the classification unit 300 after processing such as size adjustment and format conversion on the ultrasound image. Also, when an ultrasound image is input, when an ultrasound image in any one of the B mode and SWE is input, the image processing unit 220 cuts out a plurality of images to a predetermined size, for example, a size of 84x84 from the input ultrasound image. It is possible to extract and convert the extracted image from the DICOM format to the format of an image file. Because ultrasound images have different heights and widths, the performance of classifying sarcopenia is significantly degraded when resizing to make an input. Accordingly, a plurality of ultrasound images cut out from the ultrasound image to a predetermined size are reconstructed according to the size of the input layer EIL of the prediction network 410 . Then, the image processing unit 220 provides the plurality of ultrasound images processed as described above to the classification network 410 of the prediction unit 400 .

특징처리부(230)는 초음파 영상에서 방사성 특징을 추출하기 위한 것이다. 먼저, 특징처리부(230)는 B-모드 초음파 영상에서 모든 복수의 후보 방사성 특징(Radiomics features)을 추출한다. 그런 다음, 특징처리부(230)는 추출된 복수의 후보 방사성 특징으로부터 LASSO(least absolute shrinkage and selection operator) 알고리즘을 통해 루소 계수(LASSO regression coefficients)가 높은 순으로 식별망의 입력층(PIL)의 크기에 맞춰, 소정 개수, 예컨대, 30개의 방사성 특징을 도출하였다. 이러한 방사성 특징은 다음의 표 1과 같다. The feature processing unit 230 is for extracting radioactive features from the ultrasound image. First, the feature processing unit 230 extracts all of the plurality of candidate radiomics features from the B-mode ultrasound image. Then, the feature processing unit 230 from the plurality of extracted candidate radioactive features through the LASSO (least absolute shrinkage and selection operator) algorithm, the Rousseau coefficient (LASSO regression coefficients) in the order of the highest size of the input layer (PIL) of the identification network In accordance with , a predetermined number, eg, 30, radioactive features were derived. These radioactive characteristics are shown in Table 1 below.

Radiomic FeatureRadiomic Feature CoefficientsCoefficients Radiomic FeatureRadiomic Feature CoefficientsCoefficients original_FirstOrder_Skewnessoriginal_FirstOrder_Skewness 5.991432e-015.991432e-01 wavelet-LH_glszm_SmallAreaHighGrayLevelEmphasiswavelet-LH_glszm_SmallAreaHighGrayLevelEmphasis 1.496917e-011.496917e-01 wavelet-LL_FirstOrder_Skewnesswavelet-LL_FirstOrder_Skewness 5.418939e-015.418939e-01 original_gldm_LargeDependenceLowGrayLevelEmphasisoriginal_gldm_LargeDependenceLowGrayLevelEmphasis 1.229492e-011.229492e-01 wavelet-LL_glszm_GrayLevelVariancewavelet-LL_glszm_GrayLevelVariance 5.129526e-015.129526e-01 original_glcm_Autocorrelation original_glcm_Autocorrelation 1.110921e-011.110921e-01 original_glszm_GrayLevelVarianceoriginal_glszm_GrayLevelVariance 4.527734e-014.527734e-01 wavelet-LH_glszm_GrayLevelVariance wavelet-LH_glszm_GrayLevelVariance 1.080109e-011.080109e-01 original_ngtdm_Complexityoriginal_ngtdm_Complexity 4.306312e-014.306312e-01 wavelet-LH_glrlm_LongRunEmphasis wavelet-LH_glrlm_LongRunEmphasis 1.019890e-011.019890e-01 original_glszm_SmallAreaHighGrayLevelEmphasisoriginal_glszm_SmallAreaHighGrayLevelEmphasis 3.301892e-013.301892e-01 wavelet-HL_glcm_Autocorrelationwavelet-HL_glcm_Autocorrelation 9.387824e-029.387824e-02 wavelet-LH_FirstOrder_Skewnesswavelet-LH_FirstOrder_Skewness 2.804855e-012.804855e-01 wavelet-HH_gldm_LargeDependenceEmphasiswavelet-HH_gldm_LargeDependenceEmphasis 9.138463e-029.138463e-02 original_glrlm_RunVarianceoriginal_glrlm_RunVariance 2.628373e-012.628373e-01 wavelet-HL_glrlm_HighGrayLevelRunEmphasiswavelet-HL_glrlm_HighGrayLevelRunEmphasis 7.951094e-027.951094e-02 wavelet-LL_gldm_SmallDependenceHighGrayLevelEmphasis wavelet-LL_gldm_SmallDependenceHighGrayLevelEmphasis 2.403814e-012.403814e-01 wavelet-LL_glrlm_GrayLevelVariancewavelet-LL_glrlm_GrayLevelVariance 7.916654e-027.916654e-02 wavelet-LL_FirstOrder_RobustMeanAbsoluteDeviationwavelet-LL_FirstOrder_RobustMeanAbsoluteDeviation 2.086700e-012.086700e-01 wavelet-LL_FirstOrder_InterquartileRangewavelet-LL_FirstOrder_InterquartileRange 7.887526e-027.887526e-02 wavelet-HH_glszm_SmallAreaHighGrayLevelEmphasiswavelet-HH_glszm_SmallAreaHighGrayLevelEmphasis 1.927658e-011.927658e-01 wavelet-LL_gldm_LargeDependenceEmphasiswavelet-LL_gldm_LargeDependenceEmphasis 7.362614e-027.362614e-02 wavelet-LH_glrlm_RunVariancewavelet-LH_glrlm_RunVariance 1.841436e-011.841436e-01 wavelet-HH_gldm_LargeDependenceLowGrayLevelEmphasiswavelet-HH_gldm_LargeDependenceLowGrayLevelEmphasis 7.329982e-027.329982e-02 wavelet-LH_gldm_LargeDependenceEmphasiswavelet-LH_gldm_LargeDependenceEmphasis 1.569288e-011.569288e-01 original_glrlm_LongRunEmphasis original_glrlm_LongRunEmphasis 6.850687e-026.850687e-02 wavelet-LH_glrlm_LongRunLowGrayLevelEmphasiswavelet-LH_glrlm_LongRunLowGrayLevelEmphasis 1.565633e-011.565633e-01 original_glrlm_ShortRunHighGrayLevelEmphasisoriginal_glrlm_ShortRunHighGrayLevelEmphasis 6.227870e-026.227870e-02 wavelet-HL_gldm_LargeDependenceLowGrayLevelEmphasiswavelet-HL_gldm_LargeDependenceLowGrayLevelEmphasis 1.552139e-011.552139e-01 original_FirstOrder_InterquartileRangeoriginal_FirstOrder_InterquartileRange 6.074799e-026.074799e-02

다음으로, 도 4 내지 도 6을 참조하여 분류부(300)에 대해 설명한다. 전술한 바와 같이, 분류부(300)는 분류망(310)을 포함하며, 분류망(310)은 심층신경망 중 심층컨벌루션신경망(DCNN: deep convolutional neural network)이 될 수 있다. 본 발명의 분류망(310)은 복수의 계층을 포함한다. 즉, 분류망(310)은 입력층(GIL), 교번으로 반복되는 적어도 한 쌍의 컨벌루션층(GCL)과 풀링층(GPL), 적어도 하나의 완전연결층(GFL1, GFL2) 및 출력층(GOL)을 포함한다. 분류부(300)는 전처리부(200)의 영상처리부(220)로부터 골격근에 대한 초음파 영상을 수신하면, 분류망(310)을 통해 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근이 근육에코발생도(Muscle echogenicity)의 각 등급에 속할 확률을 산출하고, 산출된 확률에 따라 상기 골격근의 근육에코발생도의 등급을 분류한다. 일 실시예에 따르면, 분류망(310)은 도 4에 도시된 바와 같이, 순차로 입력층(GIL), 컨볼루션층(GCL), 풀링층(GPL), 제1 완전연결계층(GFL1), 제2 완전연결계층(GFL2) 및 출력층(GOL)을 포함한다. Next, the classification unit 300 will be described with reference to FIGS. 4 to 6 . As described above, the classification unit 300 includes a classification network 310 , and the classification network 310 may be a deep convolutional neural network (DCNN) among deep neural networks. The classification network 310 of the present invention includes a plurality of layers. That is, the classification network 310 includes an input layer (GIL), at least a pair of alternately repeated convolutional layers (GCL) and a pooling layer (GPL), at least one fully connected layer (GFL1, GFL2), and an output layer (GOL). includes When the classification unit 300 receives the ultrasound image of the skeletal muscle from the image processing unit 220 of the pre-processing unit 200, the skeletal muscle generates a muscle echo by performing a pixel unit operation on the ultrasound image through the classification network 310. The probability of belonging to each grade of the degree (Muscle echogenicity) is calculated, and the grade of the muscle echogenicity of the skeletal muscle is classified according to the calculated probability. According to an embodiment, as shown in FIG. 4 , the classification network 310 sequentially includes an input layer (GIL), a convolution layer (GCL), a pooling layer (GPL), a first fully connected layer (GFL1), It includes a second fully connected layer (GFL2) and an output layer (GOL).

컨볼루션층(GCL) 및 풀링층(GPL)은 적어도 하나의 특징영상(FM: Feature Map)으로 구성된다. 특징영상(FM)은 이전 계층의 값에 대해 가중치가 적용된 연산 수행 결과로 생성된다. 이러한 가중치는 필터 혹은 커널(W)을 통해 적용된다. 본 발명의 실시예에서 컨볼루션층(GCL)의 컨벌루션 연산은 제1 필터(W1)가 사용되며, 풀링층(PL)의 풀링 연산은 제2 필터(W2)가 사용된다. The convolutional layer (GCL) and the pooling layer (GPL) are composed of at least one feature map (FM). The feature image (FM) is generated as a result of performing an operation in which weights are applied to the values of the previous layer. These weights are applied through a filter or kernel (W). In the embodiment of the present invention, the first filter W1 is used for the convolution operation of the convolutional layer GCL, and the second filter W2 is used for the pooling operation of the pooling layer PL.

입력층(GIL)에 전처리부(100)로부터 초음파 영상이 입력되면, 컨볼루션층(GCL)은 입력층(GIL)의 초음파 영상에 대해 제1 필터(W1)를 이용한 컨벌루션(convolution) 연산 및 활성화함수에 의한 연산을 수행하여 적어도 하나의 제1 특징 영상(FM1)을 도출한다. When an ultrasound image is input from the preprocessor 100 to the input layer GIL, the convolutional layer GCL performs a convolution operation and activation on the ultrasound image of the input layer GIL using the first filter W1. At least one first feature image FM1 is derived by performing an operation by a function.

이어서, 풀링층(GPL)은 컨볼루션층(GCL)의 적어도 하나의 제1 특징 영상(FM1)에 대해 제2 필터(W2)를 이용한 풀링(pooling 또는 sub-sampling) 연산을 수행하여 적어도 하나의 제2 특징 영상(FM2)을 도출한다. Subsequently, the pooling layer GPL performs a pooling or sub-sampling operation using the second filter W2 on at least one first feature image FM1 of the convolutional layer GCL to obtain at least one A second feature image FM2 is derived.

제1 완결연결층(GFL1)은 복수의 제1 노드(f1 내지 fm)로 이루어진다. 제1 완결연결층(GFL1)의 복수의 제1 노드(f1 내지 fm)는 풀링층(GPL)의 적어도 하나의 제2 특징영상(FM2)에 대해 활성화함수에 의한 연산을 통해 복수의 제1 노드값을 산출한다. The first complete connection layer GFL1 includes a plurality of first nodes f1 to fm. The plurality of first nodes f1 to fm of the first finalized connection layer GFL1 is formed by calculating the activation function with respect to at least one second feature image FM2 of the pooling layer GPL. Calculate the value.

제2 완결연결층(GFL2)은 복수의 제2 노드(g1 내지 gn)로 이루어진다. 제2 완결연결층(GFL2)의 복수의 제2 노드(f1 내지 fm)는 제1 완전연결층(GFL1)의 복수의 제1 노드값에 대해 활성화함수에 의한 연산을 통해 복수의 제2 노드값을 산출한다. The second complete connection layer GFL2 includes a plurality of second nodes g1 to gn. The plurality of second nodes f1 to fm of the second fully connected layer GFL2 has a plurality of second node values through an operation using an activation function on the plurality of first node values of the first fully connected layer GFL1. to calculate

전술한 컨벌루션층(GCL)과 제1 및 제2 완결연결층(GFL1, GFL2)에서 사용되는 활성화함수는 시그모이드(Sigmoid), 하이퍼볼릭탄젠트(tanh: Hyperbolic tangent), ELU(Exponential Linear Unit), ReLU(Rectified Linear Unit), Leakly ReLU, Maxout, Minout, Softmax 등을 예시할 수 있다. 이러한 활성화함수 중 어느 하나를 선택하여 사용할 수 있다. The activation functions used in the above-described convolutional layer (GCL) and the first and second complete connection layers (GFL1, GFL2) are sigmoid, hyperbolic tangent (tanh), and exponential linear unit (ELU). , ReLU (Rectified Linear Unit), Leakly ReLU, Maxout, Minout, Softmax, etc. may be exemplified. Any one of these activation functions can be selected and used.

도 5 및 도 6을 참조하면, 출력층(GOL)은 2개의 출력노드(H, L)를 포함한다. 도 5를 참조하면, 제2 완결연결층(GFL2)의 복수의 제2 노드(g1 내지 gn) 각각은 가중치(W: weight)를 가지는 채널(점선으로 표시)로 출력층(OL)의 출력노드(H, L)와 연결된다. 다른 말로, 복수의 제2 노드(g1 내지 gn)의 복수의 제2 노드값은 가중치가 적용되어 출력노드(H, L)에 각각 입력된다. 이에 따라, 출력층(GOL)의 출력노드(H, L)는 제2 완결연결층(GFL2)의 가중치가 적용되는 복수의 제2 노드값에 대해 활성화함수에 의한 연산을 통해 입력층(GIL)으로 입력된 초음파 영상의 골격근이 근육에코발생도(Muscle echogenicity)의 각 등급(Grade)에 속할 확률을 산출한다. 5 and 6 , the output layer GOL includes two output nodes H and L. Referring to FIG. 5 , each of the plurality of second nodes g1 to gn of the second complete connection layer GFL2 is a channel (indicated by a dotted line) having a weight W, and an output node ( H, L) are connected. In other words, a weight is applied to the plurality of second node values of the plurality of second nodes g1 to gn and respectively input to the output nodes H and L. Accordingly, the output nodes H and L of the output layer GOL are converted to the input layer GIL through an operation by an activation function for a plurality of second node values to which the weight of the second complete connection layer GFL2 is applied. The probability that the skeletal muscle of the input ultrasound image belongs to each grade of muscle echogenicity is calculated.

한편, 도 6에 본 발명의 실시예에 따른 노드(N)가 도시되었다. 이러한 노드(N)는 분류망(310)의 출력노드(H, L) 중 어느 하나로 가정하여 설명할 것이지만, 이 노드(N)에 대한 설명은 분류망(310)에 포함된 제1 노드(f1 내지 fm), 제2 노드(f1 내지 fm) 및 출력노드(H, L)와 그 밖에 본 발명의 다른 심층신경망(410, 610)에 포함되는 모든 노드에 공통으로 적용될 수 있다. Meanwhile, a node N according to an embodiment of the present invention is illustrated in FIG. 6 . This node N will be described assuming any one of the output nodes H and L of the classification network 310 , but the description of this node N will be described with the first node f1 included in the classification network 310 . to fm), the second nodes f1 to fm, and the output nodes (H, L) and other deep neural networks 410 and 610 of the present invention may be commonly applied to all nodes included.

노드(N)는 입력된 신호 x=[x1, x2, … , xn]에 가중치 w=[w1, w2, … , wn]를 적용한 후, 그 결과에 함수 F를 취한다. 여기서, 함수 F는 활성화 함수(activation function)이다. 이때, 입력이 동일한 경우에도, 출력은 가중치(W)에 따라 다른 값이 된다. 즉, 각 노드의 출력은 다음의 수학식 1과 같다. The node N is the input signal x=[x1, x2, ... , xn] with weight w=[w1, w2, … , wn], and take a function F on the result. Here, the function F is an activation function. At this time, even when the input is the same, the output has a different value according to the weight W. That is, the output of each node is as shown in Equation 1 below.

Figure 112020066801883-pat00040
Figure 112020066801883-pat00040

설명되지 않은 파라미터 중 b는 임계치이며, 이러한 임계치는 수학식 1에서

Figure 112020066801883-pat00041
의 값이 임계치 보다 작을 때 해당 노드가 활성화되지 않도록 하는 역할을 한다. 예를 들면, 노드(N)의 이전 계층의 노드가 3개라고 가정한다. 이에 따라, 해당 노드에 대해 3개의 입력(n=3) X1, X2, X3과 3개의 가중치 W1, W2, W3이 존재한다. 노드(N)는 3개의 입력 X1, X2, X3에 대응하는 가중치 W1, W2, W3을 곱한 값을 입력받고, 모두 합산한 후, 합산된 값을 전달 함수에 대입하여 출력을 산출한다. 구체적으로, 입력 [X1, X2, X3] = 0.5, -0.3, 0이라고 가정하고, 가중치 w=[W1, W2, W3] = 4, 5, 2라고 가정한다. 또한, 설명의 편의를 위하여 전달 함수는 'sgn()'이라고 가정하면, 다음과 같이 출력값이 산출된다. Among the unexplained parameters, b is a threshold, and this threshold is
Figure 112020066801883-pat00041
It serves to prevent the corresponding node from being activated when the value of is less than the threshold. For example, it is assumed that there are three nodes in the previous layer of the node N. Accordingly, there are three inputs (n=3) X1, X2, X3 and three weights W1, W2, and W3 for the corresponding node. The node N receives a value obtained by multiplying the weights W1, W2, and W3 corresponding to the three inputs X1, X2, and X3, sums them all up, and substitutes the summed value into the transfer function to calculate an output. Specifically, it is assumed that the input [X1, X2, X3] = 0.5, -0.3, 0, and the weight w = [W1, W2, W3] = 4, 5, 2. Also, for convenience of explanation, assuming that the transfer function is 'sgn()', the output value is calculated as follows.

x1 × w1 = 0.5 × 4 = 2 x1 × w1 = 0.5 × 4 = 2

x2 × w2 = - 0.3 × 5 = -1.5x2 × w2 = - 0.3 × 5 = -1.5

x3 × w3 = 0 × 2 = 0 x3 × w3 = 0 × 2 = 0

2 + (-1.5) + 0 = 0.5 2 + (-1.5) + 0 = 0.5

sgn(0.5) = 1 sgn(0.5) = 1

이와 같이, 본 발명의 심층신경망(310, 410, 610)의 어느 하나의 노드는 이전 계층의 노드값에 가중치를 적용한 값을 입력받고, 이를 합산하여 활성화함수에 따라 연산을 수행하고, 이러한 결과를 다음 계층의 입력으로 전달한다. In this way, any one node of the deep neural network (310, 410, 610) of the present invention receives a value obtained by applying a weight to the node value of the previous layer, sums them up, performs an operation according to the activation function, and displays these results It is passed as input to the next layer.

다시, 도 5를 참조하면, 출력층(GOL)의 2개의 출력노드(H, L) 각각은 근육에코발생도(Muscle echogenicity)의 고레벨(High grade) 및 저레벨(Low grade)에 대응한다. 즉, 제1 출력노드(H)는 고레벨(High grade)에 대응하며, 제2 출력노드(L)는 저레벨(Low grade)에 대응한다. 이에 따라, 제1 출력노드(H)의 출력값은 입력층(GIL)으로 입력된 초음파 영상의 골격근이 근육에코발생도의 고레벨(High grade)에 속할 확률이며, 제2 출력노드(L)의 출력값은 입력층(GIL)으로 입력된 초음파 영상의 골격근이 근육에코발생도의 저레벨(Low grade)에 속할 확률을 나타낸다. Referring again to FIG. 5 , each of the two output nodes H and L of the output layer GOL corresponds to a high grade and a low grade of muscle echogenicity. That is, the first output node H corresponds to a high grade, and the second output node L corresponds to a low grade. Accordingly, the output value of the first output node (H) is the probability that the skeletal muscle of the ultrasound image input to the input layer (GIL) belongs to a high grade of muscle echogenicity, and the output value of the second output node (L) denotes the probability that the skeletal muscle of the ultrasound image input to the input layer (GIL) belongs to the low grade of the muscle echogenicity.

이에 따라, 도 5 및 도 6을 참조하면, 출력노드(H, L) 각각은 복수의 제2 노드(g1 내지 gn)로부터 입력되는 복수의 제2 노드값 x=[x1, x2, … , xn]에 가중치 w=[w1, w2, … , wn]를 적용한 후, 그 결과에 활성화함수 F를 취하여 출력값을 산출한다. 예컨대, 제1 및 제2 출력노드(H, L) 각각의 출력값이 각각 0.789, 0.211이면, 골격근이 근육에코발생도의 고레벨(High grade)에 속할 확률이 79%이고, 저레벨(Low grade)에 속할 확률이 21%라는 것을 나타낸다. 이와 같이, 분류망(310)이 확률(0.789, 0.211)을 출력하면, 분류부(300)는 이러한 확률(0.789, 0.211)에 따라 해당 골격근이 근육에코발생도의 고레벨(High grade)에 속하는 것으로 판단할 수 있다. 그리고 분류부(300)는 분류망(310)의 출력인 확률과, 그 확률에 따른 판단, 즉, "근육에코발생도: 고레벨(High grade)"을 리포트부(700)에 제공한다. Accordingly, referring to FIGS. 5 and 6 , each of the output nodes H and L has a plurality of second node values x=[x1, x2, ... , xn] with weight w=[w1, w2, … , wn] is applied, and the output value is calculated by taking the activation function F on the result. For example, if the output values of the first and second output nodes H and L are 0.789 and 0.211, respectively, the probability that the skeletal muscle belongs to the high grade of the muscle echogenicity is 79%, and the low grade This indicates that there is a 21% chance of belonging. As such, when the classification network 310 outputs probabilities (0.789, 0.211), the classification unit 300 determines that the corresponding skeletal muscle belongs to a high grade of muscle echogenicity according to these probabilities (0.789, 0.211). can judge In addition, the classification unit 300 provides a probability that is the output of the classification network 310 and a determination according to the probability, that is, “muscle echogenicity: high grade” to the report unit 700 .

다음으로, 도 7 및 도 8을 참조하여 예측부(400)에 대해 설명한다. 전술한 바와 같이, 예측부(400)는 예측망(410)을 포함하며, 예측망(410)은 심층신경망 중 심층컨벌루션신경망(DCNN)이 될 수 있다. 본 발명의 예측망(410)은 복수의 계층을 포함한다. 즉, 예측망(410)은 입력층(EIL), 교번으로 반복되는 적어도 한 쌍의 컨벌루션층(ECL)과 풀링층(EPL), 적어도 하나의 완전연결층(EFL1, EFL2) 및 출력층(EOL)을 포함한다. 전처리부(200)의 영상처리부(220)로부터 골격근에 대한 초음파 영상을 수신하면, 예측망(410)을 통해 초음파 영상에 대해 픽셀 단위 연산을 수행하여 초음파 영상의 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 예측한다. Next, the prediction unit 400 will be described with reference to FIGS. 7 and 8 . As described above, the prediction unit 400 includes a prediction network 410 , and the prediction network 410 may be a deep convolutional neural network (DCNN) among deep neural networks. The prediction network 410 of the present invention includes a plurality of layers. That is, the prediction network 410 includes an input layer (EIL), at least a pair of alternately repeated convolutional layers (ECL) and a pooling layer (EPL), at least one fully connected layer (EFL1, EFL2), and an output layer (EOL). includes When an ultrasound image of skeletal muscle is received from the image processing unit 220 of the pre-processing unit 200, a pixel unit operation is performed on the ultrasound image through the prediction network 410 to cause sarcopenia in the skeletal muscle of the ultrasound image. It is calculated by probability, and whether or not sarcopenia progresses is predicted according to the calculated probability.

일 실시예에 따르면, 예측망(410)은 도 7에 도시된 바와 같이, 순차로 입력층(EIL), 컨볼루션층(ECL), 풀링층(EPL), 제1 완전연결계층(EFL1), 제2 완전연결계층(EFL2) 및 출력층(EOL)을 포함한다. According to an embodiment, the prediction network 410 sequentially includes an input layer (EIL), a convolution layer (ECL), a pooling layer (EPL), a first fully connected layer (EFL1), It includes a second fully connected layer (EFL2) and an output layer (EOL).

컨볼루션층(ECL) 및 풀링층(EPL)은 적어도 하나의 특징영상(FM: Feature Map)으로 구성된다. 특징영상(FM)은 이전 계층의 값에 대해 가중치가 적용된 연산 수행 결과로 생성된다. 이러한 가중치는 필터 혹은 커널(W)을 통해 적용된다. 본 발명의 실시예에서 컨볼루션층(ECL)의 컨벌루션 연산은 제1 필터(W1)가 사용되며, 풀링층(PL)의 풀링 연산은 제2 필터(W2)가 사용된다. The convolutional layer (ECL) and the pooling layer (EPL) are composed of at least one feature map (FM). The feature image (FM) is generated as a result of performing an operation in which weights are applied to the values of the previous layer. These weights are applied through a filter or kernel (W). In the embodiment of the present invention, the first filter W1 is used for the convolution operation of the convolutional layer ECL, and the second filter W2 is used for the pooling operation of the pooling layer PL.

입력층(EIL)에 전처리부(100)로부터 초음파 영상이 입력되면, 컨볼루션층(ECL)은 입력층(EIL)의 초음파 영상에 대해 제1 필터(W1)를 이용한 컨벌루션(convolution) 연산 및 활성화함수에 의한 연산을 수행하여 적어도 하나의 제1 특징 영상(FM1)을 도출한다. When an ultrasound image is input from the preprocessor 100 to the input layer EIL, the convolutional operation and activation of the convolutional layer ECL using the first filter W1 on the ultrasound image of the input layer EIL At least one first feature image FM1 is derived by performing an operation by a function.

이어서, 풀링층(EPL)은 컨볼루션층(ECL)의 적어도 하나의 제1 특징 영상(FM1)에 대해 제2 필터(W2)를 이용한 풀링(pooling 또는 sub-sampling) 연산을 수행하여 적어도 하나의 제2 특징 영상(FM2)을 도출한다. Subsequently, the pooling layer EPL performs a pooling or sub-sampling operation using the second filter W2 on at least one first feature image FM1 of the convolution layer ECL to obtain at least one A second feature image FM2 is derived.

제1 완결연결층(EFL1)은 복수의 제1 노드(f1 내지 fm)로 이루어진다. 제1 완결연결층(EFL1)의 복수의 제1 노드(f1 내지 fm)는 풀링층(EPL)의 적어도 하나의 제2 특징영상(FM2)에 대해 활성화함수에 의한 연산을 통해 복수의 제1 노드값을 산출한다. The first complete connection layer EFL1 includes a plurality of first nodes f1 to fm. The plurality of first nodes f1 to fm of the first complete connection layer EFL1 is formed by calculating the activation function with respect to at least one second feature image FM2 of the pooling layer EPL. Calculate the value.

제2 완결연결층(EFL2)은 복수의 제2 노드(g1 내지 gn)로 이루어진다. 제2 완결연결층(EFL2)의 복수의 제2 노드(f1 내지 fm)는 제1 완전연결층(EFL1)의 복수의 제1 노드값에 대해 활성화함수에 의한 연산을 통해 복수의 제2 노드값을 산출한다. The second complete connection layer EFL2 includes a plurality of second nodes g1 to gn. The plurality of second nodes f1 to fm of the second fully connected layer EFL2 has a plurality of second node values through an operation using an activation function on the plurality of first node values of the first fully connected layer EFL1. to calculate

전술한 컨벌루션층(ECL)과 제1 및 제2 완결연결층(EFL1, EFL2)에서 사용되는 활성화함수는 시그모이드(Sigmoid), 하이퍼볼릭탄젠트(tanh: Hyperbolic tangent), ELU(Exponential Linear Unit), ReLU(Rectified Linear Unit), Leakly ReLU, Maxout, Minout, Softmax 등을 예시할 수 있다. 이러한 활성화함수 중 어느 하나를 선택하여 사용할 수 있다. The activation functions used in the above-described convolutional layer (ECL) and the first and second complete connection layers (EFL1, EFL2) are sigmoid, hyperbolic tangent (tanh), and exponential linear unit (ELU). , ReLU (Rectified Linear Unit), Leakly ReLU, Maxout, Minout, Softmax, etc. may be exemplified. Any one of these activation functions can be selected and used.

도 7 및 도 8을 참조하면, 출력층(EOL)은 2개의 출력노드(P, N)를 포함한다. 도 8에 도시된 바와 같이, 제2 완결연결층(EFL2)의 복수의 제2 노드(g1 내지 gn) 각각은 가중치(W: weight)를 가지는 채널(점선으로 표시)로 출력층(OL)의 출력노드(P, N)와 연결된다. 다른 말로, 복수의 제2 노드(g1 내지 gn)의 복수의 제2 노드값은 가중치가 적용되어 출력노드(P, N)에 각각 입력된다. 이에 따라, 출력층(EOL)의 출력노드(P, N)는 제2 완결연결층(EFL2)의 가중치가 적용되는 복수의 제2 노드값에 대해 활성화함수에 의한 연산을 통해 입력층(EIL)으로 입력된 초음파 영상의 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 산출한다. 7 and 8 , the output layer EOL includes two output nodes P and N. As shown in FIG. 8 , each of the plurality of second nodes g1 to gn of the second complete connection layer EFL2 outputs the output layer OL to a channel (indicated by a dotted line) having a weight W: It is connected to the nodes (P, N). In other words, a weight is applied to the plurality of second node values of the plurality of second nodes g1 to gn and respectively input to the output nodes P and N. Accordingly, the output nodes P and N of the output layer EOL are converted to the input layer EIL through an activation function operation for a plurality of second node values to which the weight of the second finalized connection layer EFL2 is applied. Whether sarcopenia is progressing in the skeletal muscle of the input ultrasound image is calculated as a probability.

도 8을 참조하면, 출력층(EOL)의 2개의 출력노드(P, N) 각각은 근육감소증(Sarcopenia)에 대해 양성(Positive) 및 음성(Negative)에 대응한다. 즉, 제1 출력노드(P)는 양성(Positive)에 대응하며, 제2 출력노드(N)는 음성(Negative)에 대응한다. 이에 따라, 제1 출력노드(P)의 출력값은 입력층(EIL)으로 입력된 초음파 영상의 골격근에서 근육감소증이 있음을 나타내는 확률이고, 제2 출력노드(L)의 출력값은 입력층(EIL)으로 입력된 초음파 영상의 골격근이 근육감소증이 없음을 나타내는 확률이다. Referring to FIG. 8 , each of the two output nodes P and N of the output layer EOL corresponds to positive and negative for sarcopenia. That is, the first output node P corresponds to positive, and the second output node N corresponds to negative. Accordingly, the output value of the first output node P is a probability indicating that there is sarcopenia in the skeletal muscle of the ultrasound image input to the input layer EIL, and the output value of the second output node L is the input layer EIL. It is the probability that the skeletal muscle of the ultrasound image input as , indicates that there is no sarcopenia.

이에 따라, 도 7 및 도 8을 참조하면, 출력노드(P, N) 각각은 복수의 제2 노드(g1 내지 gn)로부터 입력되는 복수의 노드값에 가중치 w=[w1, w2, … , wn]를 적용한 후, 그 결과에 활성화함수 F를 취하여 출력값을 산출한다. 예컨대, 제1 및 제2 출력노드(P, N) 각각의 출력값이 각각 0.089, 0.911이면, 해당 골격근에 근육감소증이 있음을 나타내는 확률이 9%이고, 해당 골격근에 근육감소증이 없음을 나타내는 확률이 91%이다. 이와 같이, 예측망(410)이 확률(0.089, 0.911)을 출력하면, 예측부(400)는 이러한 확률(0.089, 0.911)에 따라 해당 골격근에 근육감소증이 없는 것으로 판단할 수 있다. 그리고 예측부(400)는 예측망(410)의 출력인 확률과, 그 확률에 따른 판단, 즉, "근육감소증: 음성"을 리포트부(700)에 제공한다. Accordingly, with reference to FIGS. 7 and 8 , each of the output nodes P and N has a weight w=[w1, w2, ... , wn] is applied, and the output value is calculated by taking the activation function F on the result. For example, if the output values of the first and second output nodes P and N are 0.089 and 0.911, respectively, the probability indicating that the skeletal muscle has sarcopenia is 9%, and the probability that the corresponding skeletal muscle does not have sarcopenia is 91%. As such, when the prediction network 410 outputs probabilities (0.089, 0.911), the prediction unit 400 may determine that there is no sarcopenia in the corresponding skeletal muscle according to these probabilities (0.089, 0.911). In addition, the prediction unit 400 provides a probability that is an output of the prediction network 410 and a determination according to the probability, that is, “sarcopenia: voice” to the report unit 700 .

한편, 도 2 및 도 7을 참조하면, 시각화부(500)는 예측망(410)을 이용하여 시각화 영상을 생성하기 위한 것이다. 시각화 영상은 입력되는 초음파 영상에서 근육감소증의 원인이 되는 영역을 도출하고, 이를 다른 영역과 구별되게 표시한 영상을 의미한다. 이러한 영상은 CAM(Class Activation Map), Grad-CAM 등이 될 수 있다. 이에 따라, 예측부(400)가 예측망(410)의 출력인 확률에 따라 근육감소증을 예측하면, 예측망(410)의 컨벌루션층(ECL)의 출력값 및 경사도(gradient)를 비교하여 초음파 영상에서 근육감소증의 원인이 되는 영역을 도출하고, 초음파 영상에서 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시한다. 이로써, 시각화 영상이 생성된다. 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시할 때, 색의 종류 및 밝기를 이용하여 표시할 수 있다. 이러한 방법에 대해 보다 상세하게 설명하기로 한다. Meanwhile, referring to FIGS. 2 and 7 , the visualization unit 500 is for generating a visualization image using the prediction network 410 . The visualization image refers to an image in which a region causing sarcopenia is derived from an input ultrasound image and displayed to be distinguished from other regions. Such an image may be a CAM (Class Activation Map), Grad-CAM, or the like. Accordingly, when the prediction unit 400 predicts sarcopenia according to the probability that it is the output of the prediction network 410 , the output value and the gradient of the convolutional layer (ECL) of the prediction network 410 are compared in the ultrasound image. A region causing sarcopenia is derived, and the region causing sarcopenia in the ultrasound image is displayed separately from other regions. Thereby, a visualization image is generated. When the region causing sarcopenia is displayed separately from other regions, it can be displayed using the type and brightness of the color. This method will be described in more detail.

시각화부(500)는 먼저, 예측망(410)의 연산의 결과를 수집한 후, 다음의 수학식 2에 따라 뉴런 중요도 가중치(neuron importance weight)를 산출한다. The visualization unit 500 first collects the results of the operation of the prediction network 410 and then calculates a neuron importance weight according to Equation 2 below.

Figure 112020066801883-pat00042
Figure 112020066801883-pat00042

여기서,

Figure 112020066801883-pat00043
는 뉴런 중요도 가중치이고,
Figure 112020066801883-pat00044
는 컨벌루션층(ECL)에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고,
Figure 112020066801883-pat00045
는 컨벌루션층(ECL)의 k번째 특징지도의 좌표(i, j)의 값을 의미한다. 특히,
Figure 112020066801883-pat00046
Figure 112020066801883-pat00047
Figure 112020066801883-pat00048
에 대해 가지는 경사도(gradient)를 의미한다. here,
Figure 112020066801883-pat00043
is the neuron importance weight,
Figure 112020066801883-pat00044
is the c-th value of the output before performing the convolution operation in the convolutional layer (ECL) and performing the operation by the activation function,
Figure 112020066801883-pat00045
is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer (ECL). Especially,
Figure 112020066801883-pat00046
Is
Figure 112020066801883-pat00047
go
Figure 112020066801883-pat00048
It means the gradient with respect to .

이어서, 시각화부(500)는 다음의 수학식 3에 따라 시각화 영상을 생성한다. Next, the visualization unit 500 generates a visualization image according to Equation 3 below.

Figure 112020066801883-pat00049
Figure 112020066801883-pat00049

여기서,

Figure 112020066801883-pat00050
는 시각화 영상의 좌표(i, j)의 값을 나타낸다. 또한,
Figure 112020066801883-pat00051
는 뉴런 중요도 가중치이고,
Figure 112020066801883-pat00052
는 컨벌루션층(ECL)의 k번째 특징영상(FM)의 좌표(i, j)의 값을 의미한다. 이와 같이, 시각화부(500)는
Figure 112020066801883-pat00053
를 산출한 후,
Figure 112020066801883-pat00054
값에 따라 색을 적용함으로써, 시각화 영상을 생성할 수 있다. here,
Figure 112020066801883-pat00050
denotes the value of the coordinates (i, j) of the visualization image. In addition,
Figure 112020066801883-pat00051
is the neuron importance weight,
Figure 112020066801883-pat00052
is the value of the coordinates (i, j) of the k-th feature image FM of the convolutional layer ECL. In this way, the visualization unit 500
Figure 112020066801883-pat00053
After calculating
Figure 112020066801883-pat00054
By applying a color according to a value, a visualization image can be generated.

도 9를 참조하여 식별부(600)에 대해 설명하기로 한다. 전술한 바와 같이, 식별부(600)는 식별망(610)을 포함하며, 식별망(610)은 심층신경망 중 다층퍼셉트론(MLP: Multilayer Perceptron)이 될 수 있다. 식별망(610)은 복수의 계층(PIL, PHL, POL)을 포함한다. 이러한 복수의 계층은 입력층(PIL), 복수의 은닉 계층(PHL: HL1 내지 HLk) 및 출력층(POL)을 포함한다. The identification unit 600 will be described with reference to FIG. 9 . As described above, the identification unit 600 includes an identification network 610, and the identification network 610 may be a multilayer perceptron (MLP) among deep neural networks. The identification network 610 includes a plurality of layers (PIL, PHL, POL). The plurality of layers includes an input layer (PIL), a plurality of hidden layers (PHL: HL1 to HLk), and an output layer (POL).

또한, 복수의 계층(PIL, PHL, POL) 각각은 복수의 노드를 포함한다. 예컨대, 도시된 바와 같이, 입력층(PIL)은 n개의 입력노드(i1 ~ in)를 포함하며, 출력층(OL)은 2개의 출력노드(P, N)를 포함할 수 있다. 또한, 은닉층(PHL) 중 제1 은닉계층(PHL1)은 a개의 노드(h11 ~ h1a)를 포함하고, 제2 은닉계층(PHL2)은 b개의 노드(h21 ~ h2b)를 포함하고, 제k 은닉계층(PHLk)은 c개의 노드(hk1 ~ hkc)를 포함할 수 있다. In addition, each of the plurality of layers (PIL, PHL, and POL) includes a plurality of nodes. For example, as illustrated, the input layer PIL may include n input nodes i1 to in, and the output layer OL may include two output nodes P and N. In addition, among the hidden layers (PHL), the first hidden layer (PHL1) includes a number of nodes (h11 to h1a), the second hidden layer (PHL2) includes b nodes (h21 to h2b), and the kth hidden layer (h21 to h2b) The layer PHLk may include c nodes hk1 to hkc.

복수의 계층의 복수의 노드 모두는 연산을 가진다. 특히, 서로 다른 계층의 복수의 노드는 가중치(W: weight)를 가지는 채널(점선으로 표시)로 연결된다. 다른 말로, 어느 하나의 노드의 연산 결과는 가중치가 적용되어 다음 계층 노드의 입력이 된다. 이러한 연결 관계에 대해 앞서 도 6을 참조하여 설명하였다. 즉, 식별망(610)의 어느 한 계층의 어느 하나의 노드는 이전 계층의 노드로부터의 입력에 가중치를 적용한 값을 입력받고, 이를 합산하여 전달 함수를 취하고, 이러한 결과를 다음 계층의 입력으로 전달한다. 이에 따라, 전처리부(200)의 특징처리부(230)로부터 골격근에 대한 초음파 영상으로부터 추출된 복수의 방사성 특징(Radiomics features)이 식별망(610)의 입력층(PIL)에 입력되면, 식별망(610)은 입력된 복수의 방사성 특징에 대해 복수의 계층(PIL, PHL, POL)의 가중치가 적용되는 복수의 연산을 수행하여 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 식별한다. All of the plurality of nodes of the plurality of layers have operations. In particular, a plurality of nodes of different layers are connected by a channel (indicated by a dotted line) having a weight (W). In other words, the calculation result of one node is weighted and becomes the input of the next layer node. This connection relationship has been described above with reference to FIG. 6 . That is, any one node of any one layer of the identification network 610 receives a value obtained by applying a weight to the input from the node of the previous layer, sums them up to obtain a transfer function, and transfers the result to the input of the next layer. do. Accordingly, when a plurality of radiomics features extracted from the ultrasound image of skeletal muscle from the feature processing unit 230 of the preprocessor 200 are input to the input layer PIL of the identification network 610, the identification network ( 610) calculates whether sarcopenia is in progress in skeletal muscle with a probability by performing a plurality of calculations in which weights of a plurality of layers (PIL, PHL, POL) are applied to a plurality of input radioactive features, According to the calculated probability, it is identified whether sarcopenia is progressing.

보다 구체적으로 설명하면, 복수의 방사성 특징(Radiomics features) 각각이 식별망(610)의 입력층(PIL)의 복수의 입력노드(i1 ~ in)에 입력되면, 제1 은닉층(PHL1)의 복수의 제1 은닉노드(h11 ~ h1a) 각각은 복수의 입력노드(i1 ~ in)의 복수의 방사성 특징(Radiomics features) 각각에 가중치가 적용된 값을 입력받고(점선으로 표시), 입력된 값을 모두 합산한 후, 합산된 값에 대해 활성화함수에 따른 연산을 수행하여 복수의 제1 은닉노드값을 산출한다. 이어서, 제2 은닉층(PHL2)의 복수의 제2 은닉노드(h21 ~ h2b) 각각은 복수의 제1 은닉노드(h11 ~ h1a)의 복수의 제1 은닉노드값 각각에 가중치가 적용된 값을 입력받고(점선으로 표시), 입력된 값을 모두 합산하고, 합산된 값에 대해 활성화함수에 따른 연산을 수행하여 복수의 제2 은닉노드값을 산출한다. 이와 같은 방식으로, 은닉층(PHL) 내에서 이전의 노드값이 가중치가 적용되어 전달되고, 연산을 통해 현재의 노드값이 산출된다. 이러한 과정을 반복하여, 제k 은닉계층(PHLk)의 복수의 제k 은닉노드(hk1 ~ hkc)의 복수의 제k 은닉노드값을 산출할 수 있다. More specifically, when each of a plurality of radiomics features is input to a plurality of input nodes i1 to in of the input layer PIL of the identification network 610, a plurality of radiomics features of the first hidden layer PHL1 Each of the first hidden nodes (h11 to h1a) receives a value in which a weight is applied to each of a plurality of radiomics features of a plurality of input nodes (i1 to in) (indicated by a dotted line), and sums all the input values Then, a plurality of first hidden node values are calculated by performing an operation according to the activation function on the summed values. Then, each of the plurality of second hidden nodes h21 to h2b of the second hidden layer PHL2 receives a value in which a weight is applied to each of the plurality of first hidden node values of the plurality of first hidden nodes h11 to h1a. (indicated by a dotted line), all input values are summed, and a plurality of second hidden node values are calculated by performing an operation according to the activation function on the summed values. In this way, a weight is applied to a previous node value and transmitted in the hidden layer (PHL), and a current node value is calculated through calculation. By repeating this process, a plurality of k-th hidden node values of the plurality of k-th hidden nodes hk1 to hkc of the k-th hidden layer PHLk may be calculated.

이에 따라, 도 9를 참조하면, 출력노드(P, N) 각각은 제k 은닉계층(PHLk)의 복수의 제k 은닉노드(hk1 ~ hkc)의 복수의 제k 은닉노드값에 가중치 w=[w1, w2, … , wc×2]가 적용된 값을 입력받고(점선으로 표시), 입력된 값을 모두 합산한 후, 합산된 값에 대해 활성화함수에 따른 연산을 수행하여 출력값을 산출한다. 출력층(POL)의 2개의 출력노드(P, N) 각각은 근육감소증(Sarcopenia)에 대해 양성(Positive) 및 음성(Negative)에 대응한다. 즉, 제1 출력노드(P)는 양성(Positive)에 대응하며, 제2 출력노드(N)는 음성(Negative)에 대응한다. 이에 따라, 제1 출력노드(P)의 출력값은 입력층(PIL)으로 입력된 방사성 특징(Radiomics features)을 가지는 골격근에 근육감소증이 있음을 나타내는 확률이고, 제2 출력노드(N)의 출력값은 입력층(PIL)으로 입력된 방사성 특징을 가지는 골격근이 근육감소증이 없음을 나타내는 확률이다. Accordingly, referring to FIG. 9 , each of the output nodes P and N has a weight w = [ w1, w2, … , wc×2] is received as input (indicated by a dotted line), and after summing all the input values, an operation is performed on the summed values according to the activation function to calculate an output value. Each of the two output nodes P and N of the output layer POL corresponds to positive and negative for Sarcopenia. That is, the first output node P corresponds to positive, and the second output node N corresponds to negative. Accordingly, the output value of the first output node (P) is a probability indicating that sarcopenia is present in skeletal muscle having radiomics features input to the input layer (PIL), and the output value of the second output node (N) is It is the probability that skeletal muscle with radioactive features input into the input layer (PIL) indicates the absence of sarcopenia.

예컨대, 제1 및 제2 출력노드(P, N) 각각의 출력값이 각각 0.089, 0.911이면, 해당 골격근에 근육감소증이 있음을 나타내는 확률이 9%이고, 해당 골격근에 근육감소증이 없음을 나타내는 확률이 91%이다. 이와 같이, 식별망(610)이 확률(0.089, 0.911)을 출력하면, 식별부(600)는 이러한 확률(0.089, 0.911)에 따라 해당 골격근에 근육감소증이 없는 것으로 판단할 수 있다. 그리고 식별부(600)는 식별망(610)의 출력인 확률과, 그 확률에 따른 판단, 예컨대, "근육감소증: 음성"을 리포트부(700)에 제공한다. For example, if the output values of the first and second output nodes P and N are 0.089 and 0.911, respectively, the probability indicating that the skeletal muscle has sarcopenia is 9%, and the probability that the corresponding skeletal muscle does not have sarcopenia is 91%. As such, when the identification network 610 outputs probabilities (0.089, 0.911), the identification unit 600 may determine that there is no sarcopenia in the corresponding skeletal muscle according to these probabilities (0.089, 0.911). And the identification unit 600 provides the output of the identification network 610 with a probability and a determination according to the probability, for example, “sarcopenia: voice” to the report unit 700 .

다시 도 2를 참조하면, 리포트부(700)는 분류부(300), 예측부(400), 시각화부(500) 및 식별부(600) 각각으로부터 정보를 제공받고, 제공 받은 정보를 종합하여 표시부(13)를 통해 표시할 수 있다. 구체적으로, 분류부(300)는 골격근이 근육에코발생도의 고레벨(High grade)에 속할 확률 및 저레벨(Low grade)에 속할 확률, 그 확률에 따른 판단, 즉, "근육에코발생도: 고레벨(High grade)" 혹은 "근육에코발생도: 저레벨(Low grade)"을 리포트부(700)에 제공한다. 예측부(400)는 골격근에 근육감소증이 있을 확률 및 근육감소증이 없을 확률과, 그 확률에 따른 판단, "근육감소증: 양성" 혹은 "근육감소증: 음성"을 리포트부(700)에 제공한다. 시각화부(500)는 시각화 영상을 리포트부(700)에 제공한다. 그리고 식별부(600)는 골격근에 근육감소증이 있을 확률 및 근육감소증이 없을 확률과, 그 확률에 따른 판단, "근육감소증: 양성" 혹은 "근육감소증: 음성"을 리포트부(700)에 제공할 수 있다. Referring back to FIG. 2 , the report unit 700 receives information from each of the classification unit 300 , the prediction unit 400 , the visualization unit 500 , and the identification unit 600 , and synthesizes the received information and displays the information. It can be indicated through (13). Specifically, the classification unit 300 determines the probability that the skeletal muscle belongs to a high grade and a low grade of the muscle echogenicity, that is, a determination according to the probability, that is, "muscle echogenicity: high level ( High grade)" or "muscle echogenicity: low grade" is provided to the report unit 700 . The prediction unit 400 provides the report unit 700 with a probability of having sarcopenia in the skeletal muscle and a probability of not having sarcopenia, and a determination based on the probability, "sarcopenia: positive" or "sarcopenia: negative". The visualization unit 500 provides a visualization image to the report unit 700 . In addition, the identification unit 600 provides the report unit 700 with a probability that skeletal muscle has sarcopenia and a probability that there is no sarcopenia, and a determination according to the probability, “sarcopenia: positive” or “sarcopenia: negative” to the report unit 700 . can

다음으로, 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 방법에 대해서 설명하기로 한다. 골격근의 상태를 진단하기 위해서는 심층신경망(310, 410, 610)의 학습이 요구된다. 이러한 방법에 대해서 설명하기로 한다. 먼저, 식별망(610)을 학습시키는 방법에 대해서 설명한다. 도 10은 본 발명의 실시예에 따른 골격근의 상태를 진단하도록 식별망을 학습시키는 방법을 설명하기 위한 흐름도이다. Next, a method for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention will be described. In order to diagnose the state of the skeletal muscle, learning of the deep neural networks 310 , 410 , and 610 is required. These methods will be described. First, a method for learning the identification network 610 will be described. 10 is a flowchart illustrating a method of training an identification network to diagnose a state of a skeletal muscle according to an embodiment of the present invention.

도 10을 참조하면, 학습부(100)는 S110 단계에서 근육감소증 여부가 알려진 환자의 초음파 영상으로부터 추출된 방사성 특징을 학습 데이터, 즉, 학습용 방사성 특징으로 마련한다. 그런 다음, 학습부(100)는 S120 단계에서 학습용 방사성 특징에 알려진 근육감소증 여부에 따라 레이블을 설정한다. 이때, 근육감소증이 있는 환자의 초음파 영상으로부터 추출된 방사성 특징에는 출력노드(P, N)에 대응하여 레이블 [1, 0]을 부여한다. 반면, 근육감소증이 없는 환자의 초음파 영상으로부터 추출된 방사성 특징에는 출력노드(P, N)에 대응하여 레이블 [0, 1]을 부여한다. Referring to FIG. 10 , the learning unit 100 prepares a radioactive feature extracted from an ultrasound image of a patient whose sarcopenia is known in step S110 as learning data, that is, a radioactive feature for learning. Then, the learning unit 100 sets a label according to whether or not sarcopenia known to the radioactive characteristic for learning in step S120. At this time, a label [1, 0] is given to the radioactive feature extracted from the ultrasound image of a patient with sarcopenia in correspondence with the output node (P, N). On the other hand, radioactive features extracted from ultrasound images of patients without sarcopenia are assigned labels [0, 1] corresponding to the output nodes (P, N).

다음으로, 학습부(100)는 S130 단계에서 해당하는 학습용 방사성 특징을 식별망(610)에 입력한다. 그러면, 식별망(610)은 S140 단계에서 입력된 학습용 방사성 특징에 대해 복수의 계층(PIL, PHL, POL)의 가중치가 적용되는 복수의 연산을 수행하여 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 나타내는 출력값을 산출한다. Next, the learning unit 100 inputs the radioactive characteristic for learning corresponding to the step S130 into the identification network 610 . Then, the identification network 610 performs a plurality of calculations in which the weights of a plurality of layers (PIL, PHL, POL) are applied to the radioactive features for learning input in step S140 to determine whether sarcopenia is in progress in the skeletal muscle. An output value indicating whether or not there is a probability is calculated.

그러면, 학습부(100)는 S150 단계에서 손실함수를 통해 출력값과 레이블의 차이인 이진 교차 엔트로피 손실이 최소가 되도록 식별망(610)의 파라미터, 즉, 가중치(w), 임계치(b) 등을 수정하는 최적화를 수행할 수 있다. Then, the learning unit 100 uses the loss function in step S150 to minimize the binary cross entropy loss that is the difference between the output value and the label. You can perform optimizations to correct it.

전술한 S120 단계 내지 S150 단계는 복수의 서로 다른 학습 데이터를 이용하여 반복하여 수행되며, 이러한 반복은 평가 지표를 통해 원하는 정확도에 도달할 때까지 이루어질 수 있다. Steps S120 to S150 described above are repeatedly performed using a plurality of different learning data, and such repetition may be performed until a desired accuracy is reached through an evaluation index.

다음으로, 본 발명의 실시예에 따른 골격근의 상태를 진단하도록 분류망 및 예측망을 학습시키는 방법을 설명하기 위한 흐름도이다. 도 11은 본 발명의 실시예에 따른 골격근의 상태를 진단하도록 분류망 및 예측망을 학습시키는 방법을 설명하기 위한 흐름도이다. Next, it is a flowchart for explaining a method of training a classification network and a prediction network to diagnose a state of a skeletal muscle according to an embodiment of the present invention. 11 is a flowchart illustrating a method for training a classification network and a prediction network to diagnose a state of a skeletal muscle according to an embodiment of the present invention.

도 11을 참조하면, 학습부(100)는 S210 단계에서 심층신경망인 분류망(310) 및 예측망(410)의 초기 모델을 마련한다. 초기 모델은 입력층(GIL/EIL), 교번으로 반복되는 적어도 한 쌍의 컨벌루션층(GCL/ECL) 및 풀링층(GPL/EPL), 적어도 하나의 완전연결층(GFL/EFL) 및 출력층(GOL/EOL)을 포함한다. 특히, 초기 모델은 출력층(GOL/EOL)의 출력노드의 활성화함수로 소프트맥스(softmax)가 선택된다. 즉, 초기 모델은 활성화함수가 소프트맥스인 출력노드를 포함한다. Referring to FIG. 11 , the learning unit 100 prepares initial models of a classification network 310 and a prediction network 410 that are deep neural networks in step S210 . The initial model included an input layer (GIL/EIL), at least one pair of alternating convolutional layers (GCL/ECL) and a pooling layer (GPL/EPL), at least one fully connected layer (GFL/EFL) and an output layer (GOL). /EOL). In particular, in the initial model, a softmax is selected as an activation function of an output node of the output layer (GOL/EOL). That is, the initial model includes an output node whose activation function is a softmax.

그런 다음, 학습부(100)는 S220 단계에서 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행한다. 여기서, 공개 데이터베이스에 저장된 영상은 이미지넷(ImageNet)의 데이터베이스의 이미지를 이용할 수 있다. 초기 학습은 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 레이블을 설정하고, 학습 데이터에 대한 초기 모델의 출력값을 산출하고, 손실함수에 따라 초기 모델의 출력값과 레이블과의 차이인 손실, 예컨대, 이진교차엔트로피 손실이 최소가 되도록 초기 모델의 파라미터를 수정하는 절차이다. Then, the learning unit 100 performs initial learning using the image stored in the public database for the initial model as training data in step S220. Here, the image stored in the public database may use an image of the database of ImageNet. Initial learning uses the image stored in the public database as training data to set a label, calculates the output value of the initial model for the training data, and according to the loss function, a loss that is the difference between the output value of the initial model and the label, e.g., binary It is a procedure to modify the parameters of the initial model so that the cross-entropy loss is minimized.

초기 학습이 완료되면, 학습부(100)는 S230 단계에서 초기 모델의 상기 예측망의 은닉층의 마지막층인 마지막 완전연결층을 초기화된 완전연결층(예컨대, GFL2/EFL2)으로 교체하고, 출력층의 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 초기 모델을 정규 모델로 변환한다. When the initial learning is completed, the learning unit 100 replaces the last fully connected layer, which is the last layer of the hidden layer of the prediction network of the initial model in step S230, with an initialized fully connected layer (eg, GFL2/EFL2), and The initial model is converted into a regular model by replacing the output node that performs the softmax operation with the output node that performs the sigmoid operation.

이어서, 학습부(100)는 S240 단계에서 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행한다. 정규 학습은 초음파 영상, 즉, B 모드 초음파 영상 혹은 SWE 초음파 영상으로부터 전처리부(100)를 통해 학습용 초음파 영상을 생성하고, 학습용 초음파 영상에 대해 레이블을 설정한 후, 학습용 초음파 영상에 대해 정규 모델의 출력값을 산출하고, 손실함수에 따라 초기 모델의 출력값과 레이블과의 차이인 손실이 최소가 되도록 정규 모델의 파라미터를 수정하는 절차이다. Next, the learning unit 100 performs regular learning using the ultrasound image as training data for the regular model in step S240 . In regular learning, an ultrasound image for training is generated from an ultrasound image, that is, a B-mode ultrasound image or a SWE ultrasound image through the preprocessor 100, a label is set for the ultrasound image for training, and the normal model is applied to the ultrasound image for training. It is a procedure of calculating the output value and modifying the parameters of the regular model so that the loss, which is the difference between the output value of the initial model and the label, is minimized according to the loss function.

그러면, 전술한 정규 모델에 대한 학습 방법에 대해 보다 상세하게 설명하기로 한다. 도 12는 본 발명의 실시예에 따른 골격근의 상태를 진단하도록 분류망 및 예측망에 대한 정규 학습을 수행하는 방법을 설명하기 위한 흐름도이다. 도 13은 본 발명의 실시예에 따른 골격근의 상태를 진단하도록 분류망 및 예측망에 대한 정규 학습을 수행하는 방법을 설명하기 위한 벡터공간의 예이다. Then, the learning method for the aforementioned regular model will be described in more detail. 12 is a flowchart illustrating a method of performing regular learning on a classification network and a prediction network to diagnose a state of a skeletal muscle according to an embodiment of the present invention. 13 is an example of a vector space for explaining a method of performing regular learning for a classification network and a prediction network to diagnose a state of a skeletal muscle according to an embodiment of the present invention.

도 12를 참조하면, 학습부(100)는 S310 단계에서 학습 데이터, 즉, 학습용 초음파 영상을 마련한다. 학습용 초음파 영상은 분류망(310)인 경우, 근육에코발생도의 등급이 알려진 환자의 초음파 영상이며, 예측망(410)인 경우, 근육감소증이 있는지 여부가 알려진 환자의 초음파 영상이다. 이와 같은, 학습용 초음파 영상이 마련되면, 학습부(100)는 다음의 수학식 4와 같은 손실함수를 이용하여 학습을 수행한다. Referring to FIG. 12 , the learning unit 100 prepares learning data, that is, an ultrasound image for learning in step S310 . In the case of the classification network 310 , the ultrasound image for learning is an ultrasound image of a patient whose level of muscle echogenicity is known, and in the case of the prediction network 410 , it is an ultrasound image of a patient whose sarcopenia is known. When such an ultrasound image for learning is prepared, the learning unit 100 performs learning using a loss function as shown in Equation 4 below.

Figure 112020066801883-pat00055
Figure 112020066801883-pat00055

수학식 4에서 E(y, f)는 손실함수를 나타낸다. yi는 출력층의 출력값이고, ei는 출력층의 출력값에 대응하는 분류 레이블이다. gij는 은닉층의 마지막층의 출력값이고, cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이다. i는 출력층의 노드에 대응하는 인덱스이고, j는 은닉층의 마지막층의 노드에 대응하는 인덱스이다. 그리고

Figure 112020066801883-pat00056
는 하이퍼파라미터를 나타낸다. In Equation 4, E(y, f) represents a loss function. yi is an output value of the output layer, and ei is a classification label corresponding to the output value of the output layer. gij is the output value of the last layer of the hidden layer, and cij is the hidden label corresponding to the output value of the last layer of the hidden layer. i is the index corresponding to the node of the output layer, and j is the index corresponding to the node of the last layer of the hidden layer. and
Figure 112020066801883-pat00056
is a hyperparameter.

학습부(100)는 S320 단계에서 분류 손실 최적화를 위한 설정을 수행한다. 이때, 학습부(100)는 하이퍼파라미터 및 분류 레이블을 설정한다. 이러한 S320 단계에서 학습부(100)는 하이퍼파라미터

Figure 112020066801883-pat00057
는 0으로 설정한다. 이에 따라, 분류 손실을 최적화하는 학습에서 수학식 4의 손실함수는 다음의 수학식 5와 같다. The learning unit 100 performs settings for optimizing classification loss in step S320. At this time, the learning unit 100 sets a hyperparameter and a classification label. In this step S320, the learning unit 100 is a hyperparameter
Figure 112020066801883-pat00057
is set to 0. Accordingly, in the learning to optimize the classification loss, the loss function of Equation 4 is the following Equation 5.

Figure 112020066801883-pat00058
Figure 112020066801883-pat00058

또한, S320 단계에서 학습부(100)는 분류 레이블을 설정한다. 분류망(310)인 경우, 근육에코발생도의 등급이 고레벨로 알려진 환자의 초음파 영상에는 출력노드(H, L)에 대응하여 분류 레이블 [1, 0]을 부여하고, 저레벨로 알려진 환자의 초음파 영상에는 출력노드(H, L)에 대응하여 분류 레이블 [0, 1]을 부여한다. 또한, 예측망(410)인 경우, 근육감소증이 있는 것으로 알려진 환자의 초음파 영상에는 출력노드(P, N)에 대응하여 분류 레이블 [1, 0]을 부여하고, 근육감소증이 없는 것으로 알려진 환자의 초음파 영상에는 출력노드(P, N)에 대응하여 분류 레이블 [0, 1]을 부여한다. Also, in step S320, the learning unit 100 sets a classification label. In the case of the classification network 310, a classification label [1, 0] is assigned to an ultrasound image of a patient known to have a high level of muscle echogenicity in response to the output nodes (H, L), and ultrasound of a patient known as a low level A classification label [0, 1] is given to an image corresponding to the output node (H, L). In addition, in the case of the prediction network 410, a classification label [1, 0] is assigned to an ultrasound image of a patient known to have sarcopenia in response to the output nodes (P, N), and A classification label [0, 1] is given to the ultrasound image corresponding to the output nodes (P, N).

전술한 바와 같이 하이퍼파라미터 및 분류 레이블이 설정되면, 학습부(100)는 S330 단계에서 분류 손실을 최적화하는 학습을 수행한다. 이때, 학습부(100)는 분류망(310) 및 식별망(410)의 정규 모델에 학습용 초음파 영상을 입력한다. 그러면, 정규 모델은 학습용 초음파 영상에 대해 복수의 계층 간 가중치가 적용되는 복수의 연산을 통해 출력값을 산출한다. 그러면, 학습부(100)는 하이퍼파라미터가 0으로 설정된 수학식 5와 같은 손실함수를 이용하여 앞서 설정된 분류 레이블과 출력층의 출력값의 차이를 나타내는 수학식 5의 손실함수의 분류 손실이 최소가 되도록 정규 모델의 파라미터를 수정하는 분류 손실 최적화를 수행한다. 이러한 분류 손실 최적화는 복수의 서로 다른 학습용 초음파 영상을 이용하여 반복하여 수행되며, 이러한 반복은 평가 지표를 통해 원하는 정확도에 도달할 때까지 이루어질 수 있다. When the hyperparameter and the classification label are set as described above, the learning unit 100 performs learning to optimize the classification loss in step S330 . At this time, the learning unit 100 inputs the training ultrasound image to the regular model of the classification network 310 and the identification network 410 . Then, the regular model calculates an output value through a plurality of operations in which a plurality of inter-layer weights are applied to the training ultrasound image. Then, the learning unit 100 is normalized so that the classification loss of the loss function of Equation 5 indicating the difference between the previously set classification label and the output value of the output layer is minimized by using the loss function as in Equation 5 in which the hyperparameter is set to 0. A classification loss optimization that modifies the parameters of the model is performed. This classification loss optimization is repeatedly performed using a plurality of different ultrasound images for training, and such repetition may be performed until a desired accuracy is reached through an evaluation index.

전술한 바와 같은 분류 손실 최적화에 따른 학습이 완료되면, 학습부(100)는 S340 단계에서 정규 모델의 분류 대상 클래스 별로 은닉벡터 중 대표벡터를 선택한다. 대표벡터를 선택하는 방법은 다음과 같다. 앞서(S330) 사용된 학습용 초음파 영상을 분류 손실 최적화가 이루어진 정규 모델에 다시 입력하여 정규 모델의 은닉층의 마지막층(예컨대, GFL2, EFL2)의 출력값을 도출한다. 도 5 및 도 8을 참조하면, 은닉층의 마지막층(예컨대, GFL2, EFL2)의 출력은 은닉벡터 G[g1, g2, g3, …, gn]이다. 학습부(100)는 복수의 은닉벡터를 소정의 벡터공간에 임베딩한다. 도 13 (A)에 소정의 벡터공간에 복수의 은닉벡터가 임베딩된 예를 도시하였다. 기본적으로, 분류 손실 최적화에 의해 벡터공간 상에서 분류선(Div)에 의해 복수의 은닉벡터 각 클래스별로 어느 정도의 구분이 이루어진다. 여기서, 분류 대상인 클래스는 분류망(310)인 경우, 근육에코발생도의 등급이 고레벨(H)인 클래스 및 저레벨(L)인 클래스를 포함하며, 예측망(410)인 경우, 근육감소증 양성(P)인 클래스 및 음성(N)인 클래스를 포함한다. 즉, 도 13의 (A)에서 제1 클래스는 분류망(310)인 경우, 근육에코발생도의 등급이 고레벨(H)인 클래스이고, 예측망(410)인 경우, 근육감소증이 양성(P)인 클래스일 수 있다. 또한, 도 13의 (A)에서 제2 클래스는 근육에코발생도의 등급이 저레벨(L)인 클래스이고, 예측망(410)인 경우, 근육감소증이 음성(N)인 클래스이다. 학습부(100)는 이러한 벡터 공간에 임베딩된 복수의 은닉벡터 중 각 클래스 별로 대표벡터를 선정할 수 있다. 대표벡터는 13의 (A)에 도시된 바와 같이, 동일한 클래스 내의 은닉벡터 중 중간값을 선정하는 것이 바람직하다. When the learning according to the classification loss optimization as described above is completed, the learning unit 100 selects a representative vector from among the hidden vectors for each classification target class of the regular model in step S340 . The method of selecting a representative vector is as follows. By re-inputting the training ultrasound image used above (S330) into the regular model with classification loss optimization, output values of the last layers (eg, GFL2, EFL2) of the hidden layer of the regular model are derived. 5 and 8, the output of the last layer of the hidden layer (eg, GFL2, EFL2) is the hidden vector G[g1, g2, g3, ... , gn]. The learning unit 100 embeds a plurality of hidden vectors in a predetermined vector space. 13A shows an example in which a plurality of hidden vectors are embedded in a predetermined vector space. Basically, some degree of classification is made for each class of a plurality of hidden vectors by the classification line Div in the vector space by classification loss optimization. Here, in the case of the classification network 310, the class to be classified includes a class having a high level (H) and a class having a low level (L), and in the case of the prediction network 410, sarcopenia positive ( It includes a class that is P) and a class that is negative (N). That is, in (A) of FIG. 13 , the first class is a class having a high level (H) of muscle echogenicity in the case of the classification network 310, and in the case of the prediction network 410, sarcopenia is positive (P). ) can be a class. In addition, in FIG. 13A , the second class is a class with a low level (L) of muscle echogenicity, and in the case of the prediction network 410, a class in which sarcopenia is negative (N). The learning unit 100 may select a representative vector for each class from among a plurality of hidden vectors embedded in such a vector space. As the representative vector, as shown in (A) of 13, it is preferable to select an intermediate value among the hidden vectors in the same class.

다음으로, 학습부(100)는 S350 단계에서 복합 손실 최적화를 위한 설정을 수행한다. 이때, 학습부(100)는 하이퍼파라미터, 분류 레이블 및 은닉 레이블을 설정한다. 이러한 S350 단계에서 학습부(100)는 하이퍼파라미터

Figure 112020066801883-pat00059
는 0.5로 설정한다. 이에 따라, 분류 손실을 최적화하는 학습에서 수학식 4의 손실함수는 다음의 수학식 6과 같다. Next, the learning unit 100 performs settings for optimizing the complex loss in step S350 . At this time, the learning unit 100 sets a hyperparameter, a classification label, and a hidden label. In this step S350, the learning unit 100 is a hyperparameter
Figure 112020066801883-pat00059
is set to 0.5. Accordingly, in the learning to optimize the classification loss, the loss function of Equation 4 is as shown in Equation 6 below.

Figure 112020066801883-pat00060
Figure 112020066801883-pat00060

또한, S350 단계에서 학습부(100)는 분류 레이블 및 은닉 레이블을 설정한다. 분류 레이블의 설정은, 분류망(310)인 경우, 근육에코발생도의 등급이 고레벨로 알려진 환자의 초음파 영상에는 출력노드(H, L)에 대응하여 분류 레이블 [1, 0]을 부여하고, 저레벨로 알려진 환자의 초음파 영상에는 출력노드(H, L)에 대응하여 분류 레이블 [0, 1]을 부여한다. 또한, 예측망(410)인 경우, 근육감소증이 있는 것으로 알려진 환자의 초음파 영상에는 출력노드(P, N)에 대응하여 분류 레이블 [1, 0]을 부여하고, 근육감소증이 없는 것으로 알려진 환자의 초음파 영상에는 출력노드(P, N)에 대응하여 분류 레이블 [0, 1]을 부여한다. Also, in step S350, the learning unit 100 sets a classification label and a hidden label. As for the classification label setting, in the case of the classification network 310, a classification label [1, 0] is given to the ultrasound image of a patient known to have a high level of muscle echogenicity in response to the output nodes (H, L), A classification label [0, 1] is assigned to an ultrasound image of a patient known as a low level corresponding to the output nodes (H, L). In addition, in the case of the prediction network 410, a classification label [1, 0] is assigned to an ultrasound image of a patient known to have sarcopenia in response to the output nodes (P, N), and A classification label [0, 1] is given to the ultrasound image corresponding to the output nodes (P, N).

특히, 은닉 레이블의 설정은, 앞서(S340) 선정한 대표벡터를 이용한다. 이에 따라, 각 클래스별로 대표벡터가 은닉 레이블로 설정된다. In particular, the setting of the hidden label uses the representative vector selected above (S340). Accordingly, a representative vector for each class is set as a hidden label.

전술한 바와 같이 하이퍼파라미터, 분류 레이블 및 은닉 레이블이 설정되면, 전술한 바와 같이 파라미터가 설정되면, 학습부(100)는 S360 단계에서 복합 손실을 최적화하는 학습을 수행한다. 이때, 학습부(100)는 분류망(310) 및 식별망(410)의 정규 모델에 학습용 초음파 영상을 입력한다. 그러면, 정규 모델은 학습용 초음파 영상에 대해 복수의 계층 간 가중치가 적용되는 복수의 연산을 통해 은닉벡터 및 출력값을 산출한다. 그러면, 학습부(100)는 하이퍼파라미터가 0.5로 설정된 수학식 6과 같은 손실함수를 이용하여 수학식 6의 손실함수의 복합 손실이 최소가 되도록 정규 모델의 파라미터를 수정하는 복합 손실 최적화를 수행한다. 이러한 복합 손실은 분류 레이블과 출력층의 출력값의 차이를 나타내는 분류 손실 및 은닉 벡터와 대표 벡터와의 차이를 나타내는 은닉 손실을 포함한다. 이러한 복합 손실 최적화는 복수의 서로 다른 학습용 초음파 영상을 이용하여 반복하여 수행되며, 이러한 반복은 평가 지표를 통해 원하는 정확도에 도달할 때까지 이루어질 수 있다. When the hyperparameter, classification label, and hidden label are set as described above, when the parameter is set as described above, the learning unit 100 performs learning to optimize the complex loss in step S360. At this time, the learning unit 100 inputs the training ultrasound image to the regular model of the classification network 310 and the identification network 410 . Then, the regular model calculates a hidden vector and an output value through a plurality of operations in which a plurality of inter-layer weights are applied to the training ultrasound image. Then, the learning unit 100 performs complex loss optimization by modifying the parameters of the regular model so that the complex loss of the loss function of Equation 6 is minimized by using the loss function as in Equation 6 in which the hyperparameter is set to 0.5. . The complex loss includes a classification loss indicating a difference between a classification label and an output value of an output layer, and a hiding loss indicating a difference between a hidden vector and a representative vector. The complex loss optimization is repeatedly performed using a plurality of different ultrasound images for training, and such repetition may be performed until a desired accuracy is reached through an evaluation index.

이러한 복합 손실에 의한 학습이 완료되면, 벡터 공간 상의 복수의 은닉벡터는 도 13의 (B)와 같이 분류될 수 있다. 즉, 복수의 은닉벡터는 벡터 공간 상에서 대표 벡터를 향하여 이동된다. 이에 따라, 클래스 별 구분이 보다 명확해짐을 알 수 있다. 즉, 복합 손실에 의한 학습에 의해 분류망(310) 및 예측망(410)의 분류 성능이 향상될 수 있다. When the learning by the complex loss is completed, a plurality of hidden vectors in the vector space may be classified as shown in FIG. 13B. That is, the plurality of hidden vectors are moved toward the representative vector in the vector space. Accordingly, it can be seen that the classification for each class becomes clearer. That is, the classification performance of the classification network 310 and the prediction network 410 may be improved by learning by complex loss.

다음으로, 전술한 바와 같은 방법으로 학습이 완료된 분류망(310), 예측망(410) 및 식별망(610)을 이용하여 골격근의 상태를 진단하는 방법에 대해서 설명하기로 한다. 도 14는 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 방법을 설명하기 위한 흐름도이다. Next, a method of diagnosing the state of skeletal muscle using the classification network 310 , the prediction network 410 , and the identification network 610 that has been trained in the manner described above will be described. 14 is a flowchart illustrating a method for diagnosing a state of skeletal muscle from an ultrasound image using a deep neural network according to an embodiment of the present invention.

도 14를 참조하면, 진단장치(10)의 전처리부(200)는 S410 단계에서 인터페이스부(11)를 통해 진단영상장치(20)로부터 초음파 영상을 입력받을 수 있다. 입력된 초음파 영상은 B 모드 초음파 영상 및 SWE 초음파 영상 중 어느 하나가 될 수 있다. Referring to FIG. 14 , the preprocessor 200 of the diagnostic apparatus 10 may receive an ultrasound image from the diagnostic imaging apparatus 20 through the interface 11 in operation S410 . The input ultrasound image may be any one of a B-mode ultrasound image and an SWE ultrasound image.

그러면, 전처리부(200)는 S420 단계에서 초음파 영상에 대한 전처리를 수행한다. 이때, 전처리부(200)는 초음파 영상의 크기를 조정하거나, 소정의 크기로 잘라낸 후, 이미지 파일 포맷으로 변환할 수 있다. 또한, 전처리부(200)는 초음파 영상으로부터 방사상 특징(Radiomics features)을 추출할 수 있다. 전처리부(200)는 전처리가 완료되면, 전처리된 초음파 영상 혹은 방사상 특징(Radiomics features)을 분류부(300), 예측부(400) 및 식별부(600)에 제공한다. Then, the preprocessor 200 performs preprocessing on the ultrasound image in step S420 . In this case, the preprocessor 200 may adjust the size of the ultrasound image or cut it to a predetermined size, and then convert it into an image file format. Also, the preprocessor 200 may extract radial features from the ultrasound image. When the preprocessing is completed, the preprocessing unit 200 provides the preprocessed ultrasound image or radiomics features to the classifying unit 300 , the predicting unit 400 , and the identifying unit 600 .

분류부(300), 예측부(400), 시각화부(500) 및 식별부(600)는 S430 단계에서 골격근에 대한 평가를 수행하여 평가 결과를 도출한다. 즉, 분류부(300)는 분류망(310)을 통해 초음파 영상에 대해 픽셀 단위 연산을 수행하여 골격근이 근육에코발생도(Muscle echogenicity)의 각 등급에 속할 확률을 산출하고, 산출된 확률에 따라 골격근의 근육에코발생도의 등급을 분류한다. 또한, 예측부(400)는 예측망(410)을 통해 초음파 영상에 대해 픽셀 단위 연산을 수행하여 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 예측한다. 그리고 시각화부(500)는 예측망(410)의 컨벌루션층(ECL)의 출력값 및 경사도(gradient)를 비교하여 초음파 영상에서 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시하는 시각화 영상을 생성한다. 또한, 식별부(600)는 식별망(610)을 통해 복수의 방사성 특징에 대한 연산을 수행하여 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 식별한다. The classification unit 300 , the prediction unit 400 , the visualization unit 500 , and the identification unit 600 perform evaluation on the skeletal muscle in step S430 to derive an evaluation result. That is, the classification unit 300 calculates the probability that the skeletal muscle belongs to each class of muscle echogenicity by performing pixel unit operation on the ultrasound image through the classification network 310 , and according to the calculated probability Classify the muscle echogenicity of skeletal muscle. In addition, the prediction unit 400 calculates with a probability whether sarcopenia is progressing in the skeletal muscle by performing a pixel unit operation on the ultrasound image through the prediction network 410 , and determines whether sarcopenia is progressing according to the calculated probability. predict And the visualization unit 500 compares the output value and the gradient of the convolutional layer (ECL) of the prediction network 410 to generate a visualization image that displays the region causing sarcopenia in the ultrasound image by distinguishing it from other regions. do. In addition, the identification unit 600 calculates with a probability whether sarcopenia is in progress in the skeletal muscle by performing an operation on a plurality of radioactive features through the identification network 610, and determines whether sarcopenia is in progress according to the calculated probability. identify

다음으로, 리포트부(700)는 S440 단계에서 분류부(300), 예측부(400), 시각화부(500) 및 식별부(600) 각각으로부터 정보를 제공받고, 제공 받은 정보를 종합한 종합 리포트를 표시부(13)를 통해 표시할 수 있다. 여기서, 종합 리포트는 근육에코발생도의 고레벨(High grade)에 속할 확률 및 저레벨(Low grade)에 속할 확률, 그 확률에 따른 판단, 즉, "근육에코발생도: 고레벨(High grade)" 혹은 "근육에코발생도: 저레벨(Low grade)"와, 초음파 영상을 기초로 예측되는 근육감소증이 있을 확률 및 근육감소증이 없을 확률과, 그 확률에 따른 판단, "근육감소증: 양성" 혹은 "근육감소증: 음성"과, 시각화 영상과, 방사상 특징을 기초로 식별된 근육감소증이 있을 확률 및 근육감소증이 없을 확률과, 그 확률에 따른 판단, "근육감소증: 양성" 혹은 "근육감소증: 음성"을 포함한다. Next, the report unit 700 receives information from each of the classification unit 300 , the prediction unit 400 , the visualization unit 500 , and the identification unit 600 in step S440 , and a comprehensive report synthesizing the received information may be displayed through the display unit 13 . Here, the comprehensive report indicates the probability of belonging to a high grade and a probability of belonging to a low grade of muscle echogenicity, and a judgment based on the probability, that is, "muscle echogenicity: high grade" or " Muscular echogenicity: Low grade, and the probability of having sarcopenia predicted based on ultrasound imaging and the probability of not having sarcopenia, and judging according to the probability, "Sarcopenia: benign" or "Sarcopenia: negative", the probability of having sarcopenia and the absence of sarcopenia identified on the basis of the visual image and the radial features, and the judgment according to the probability, "sarcopenia: positive" or "sarcopenia: negative" .

한편, 전술한 본 발명의 실시예에 따른 방법은 다양한 컴퓨터수단을 통하여 판독 가능한 프로그램 형태로 구현되어 컴퓨터로 판독 가능한 기록매체에 기록될 수 있다. 여기서, 기록매체는 프로그램 명령, 데이터 파일, 데이터구조 등을 단독으로 또는 조합하여 포함할 수 있다. 기록매체에 기록되는 프로그램 명령은 본 발명을 위하여 특별히 설계되고 구성된 것들이거나 컴퓨터 소프트웨어 당업자에게 공지되어 사용 가능한 것일 수도 있다. 예컨대 기록매체는 하드 디스크, 플로피 디스크 및 자기 테이프와 같은 자기 매체(magnetic media), CD-ROM, DVD와 같은 광 기록 매체(optical media), 플롭티컬 디스크(floptical disk)와 같은 자기-광 매체(magneto-optical media) 및 롬(ROM), 램(RAM), 플래시 메모리 등과 같은 프로그램 명령을 저장하고 수행하도록 특별히 구성된 하드웨어 장치를 포함한다. 프로그램 명령의 예에는 컴파일러에 의해 만들어지는 것과 같은 기계어 와이어뿐만 아니라 인터프리터 등을 사용해서 컴퓨터에 의해서 실행될 수 있는 고급 언어 와이어를 포함할 수 있다. 이러한 하드웨어 장치는 본 발명의 동작을 수행하기 위해 하나 이상의 소프트웨어 모듈로서 작동하도록 구성될 수 있으며, 그 역도 마찬가지이다. Meanwhile, the above-described method according to an embodiment of the present invention may be implemented in the form of a program readable by various computer means and recorded in a computer readable recording medium. Here, the recording medium may include a program command, a data file, a data structure, etc. alone or in combination. The program instructions recorded on the recording medium may be specially designed and configured for the present invention, or may be known and available to those skilled in the art of computer software. For example, the recording medium includes magnetic media such as hard disks, floppy disks and magnetic tapes, optical recording media such as CD-ROMs and DVDs, and magneto-optical media such as floppy disks ( magneto-optical media) and hardware devices specially configured to store and execute program instructions such as ROM, RAM, flash memory, and the like. Examples of program instructions may include not only machine language wires such as those generated by a compiler, but also high-level language wires that can be executed by a computer using an interpreter or the like. Such hardware devices may be configured to operate as one or more software modules to perform the operations of the present invention, and vice versa.

이상 본 발명을 몇 가지 바람직한 실시예를 사용하여 설명하였으나, 이들 실시예는 예시적인 것이며 한정적인 것이 아니다. 이와 같이, 본 발명이 속하는 기술분야에서 통상의 지식을 지닌 자라면 본 발명의 사상과 첨부된 특허청구범위에 제시된 권리범위에서 벗어나지 않으면서 균등론에 따라 다양한 변화와 수정을 가할 수 있음을 이해할 것이다. Although the present invention has been described above using several preferred embodiments, these examples are illustrative and not restrictive. As such, those of ordinary skill in the art to which the present invention pertains will understand that various changes and modifications can be made in accordance with the doctrine of equivalents without departing from the spirit of the present invention and the scope of rights set forth in the appended claims.

10: 관제장치 20: 진단영상장치
100: 학습부 200: 전처리부
210: 증강부 220: 영상처리부
230: 특징처리부 300: 분류부
310: 분류망 400: 예측부
410: 예측망 500: 시각화부
600: 식별부 610: 식별망
700: 리포트부
10: control device 20: diagnostic imaging device
100: learning unit 200: preprocessing unit
210: augmentation unit 220: image processing unit
230: feature processing unit 300: classification unit
310: classification network 400: prediction unit
410: prediction network 500: visualization unit
600: identification unit 610: identification network
700: report unit

Claims (19)

골격근의 상태를 진단하기 위한 장치에 있어서,
골격근에 대한 초음파 영상이 입력되면, 소정의 형식에 따라 상기 초음파 영상을 가공하는 전처리부; 및
입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하는 예측망을 통해
상기 가공된 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 산출하고,
산출된 확률에 따라 근육감소증 진행 여부를 예측하는 예측부;
를 포함하며,
상기 전처리부가 상기 골격근에 대한 초음파 영상으로부터 복수의 방사성 특징을 추출하면,
입력층, 은닉층 및 출력층을 포함하는 식별망을 통해 상기 복수의 방사성 특징에 대한 연산을 수행하여 상기 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 산출하고,
산출된 확률에 따라 근육감소증 진행 여부를 식별하는 식별부;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
In the apparatus for diagnosing the condition of skeletal muscle,
When an ultrasound image of skeletal muscle is input, a preprocessor processing the ultrasound image according to a predetermined format; and
Through a prediction network comprising an input layer, at least a pair of alternating convolutional and pooling layers, at least one fully connected layer, and an output layer
By performing a pixel unit operation on the processed ultrasound image, it is calculated as a probability whether Sarcopenia is in progress in the skeletal muscle,
a predictor for predicting whether sarcopenia progresses according to the calculated probability;
includes,
When the preprocessor extracts a plurality of radioactive features from the ultrasound image of the skeletal muscle,
Calculating with a probability whether Sarcopenia is in progress in the skeletal muscle by performing an operation on the plurality of radioactive features through an identification network including an input layer, a hidden layer, and an output layer,
an identification unit for identifying whether sarcopenia progresses according to the calculated probability;
characterized in that it further comprises
A device for diagnosing the condition of skeletal muscle.
제1항에 있어서,
상기 예측부가 근육감소증을 예측하면,
상기 예측망의 상기 컨벌루션층의 출력값 및 경사도(gradient)를 비교하여 상기 초음파 영상에서 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시하는 시각화 영상을 생성하는
시각화부;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
According to claim 1,
If the prediction unit predicts sarcopenia,
Comparing the output value and the gradient of the convolutional layer of the prediction network to generate a visualization image displaying a region causing sarcopenia in the ultrasound image by distinguishing it from other regions
visualization unit;
characterized in that it further comprises
A device for diagnosing the condition of skeletal muscle.
제2항에 있어서,
상기 시각화부는
상기 예측망의 연산의 결과를 수집한 후,
수학식
Figure 112020066801883-pat00061

에 따라 뉴런 중요도 가중치를 산출하고,
수학식
Figure 112020066801883-pat00062

에 따라 상기 시각화 영상을 생성하며,
상기
Figure 112020066801883-pat00063
는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고,
상기
Figure 112020066801883-pat00064
는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고,
상기
Figure 112020066801883-pat00065
Figure 112020066801883-pat00066
Figure 112020066801883-pat00067
에 대해 가지는 경사도(gradient)이고,
상기
Figure 112020066801883-pat00068
는 상기 뉴런 중요도 가중치이고,
상기
Figure 112020066801883-pat00069
는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
3. The method of claim 2,
The visualization unit
After collecting the results of the calculation of the prediction network,
formula
Figure 112020066801883-pat00061

Calculate the neuron importance weight according to
formula
Figure 112020066801883-pat00062

to generate the visualization image according to
remind
Figure 112020066801883-pat00063
is the c-th value of the output before performing the convolution operation in the convolutional layer and performing the operation by the activation function,
remind
Figure 112020066801883-pat00064
is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer,
remind
Figure 112020066801883-pat00065
Is
Figure 112020066801883-pat00066
go
Figure 112020066801883-pat00067
is the gradient with respect to
remind
Figure 112020066801883-pat00068
is the neuron importance weight,
remind
Figure 112020066801883-pat00069
is the value of the coordinates (i, j) of the visualization image
A device for diagnosing the condition of skeletal muscle.
삭제delete 제1항에 있어서,
상기 전처리부가 상기 골격근에 대한 초음파 영상을 상기 소정의 형식과 다른 형식으로 상기 초음파 영상을 가공하면,
입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하는 분류망을 통해
상기 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근이 근육에코발생도(Muscle echogenicity)의 각 등급에 속할 확률을 산출하고,
산출된 확률에 따라 상기 골격근의 근육에코발생도의 등급을 분류하는 분류부;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
According to claim 1,
When the preprocessor processes the ultrasound image of the skeletal muscle in a format different from the predetermined format,
Through a classification network comprising an input layer, at least a pair of alternating convolutional and pooling layers, at least one fully connected layer, and an output layer
calculating the probability that the skeletal muscle belongs to each class of muscle echogenicity by performing pixel unit operation on the ultrasound image,
a classification unit for classifying the degree of muscle echogenicity of the skeletal muscle according to the calculated probability;
characterized in that it further comprises
A device for diagnosing the condition of skeletal muscle.
제1항에 있어서,
상기 예측망의 상기 출력층이 소프트맥스 연산을 수행하는 출력노드로 이루어진 상기 예측망의 초기 모델을 마련하고, 상기 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행하며,
상기 초기 학습이 완료되면, 상기 예측망의 마지막 완전연결층을 초기화된 완전연결층으로 교체하고, 상기 출력층의 상기 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 상기 예측망의 정규 모델을 마련하고, 상기 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행하는 학습부;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
According to claim 1,
The output layer of the prediction network prepares an initial model of the prediction network consisting of an output node that performs a softmax operation, and performs initial learning using an image stored in a public database for the initial model as training data,
When the initial learning is completed, the last fully connected layer of the prediction network is replaced with an initialized fully connected layer, and the output node performing the softmax operation of the output layer is replaced with an output node performing the sigmoid operation. a learning unit that prepares a regular model of the prediction network and performs regular learning on the regular model by using an ultrasound image as training data;
characterized in that it further comprises
A device for diagnosing the condition of skeletal muscle.
제6항에 있어서,
상기 학습부는
상기 예측망에 대해 상기 정규 학습을 수행할 때,
손실함수
Figure 112022026989181-pat00070

를 이용하여
상기
Figure 112022026989181-pat00071
을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하고,
상기
Figure 112022026989181-pat00072
을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행하고,
상기 E(y, f)는 손실함수이고,
상기 yi는 출력층의 출력값이고,
상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고,
상기 gij는 은닉층의 마지막층의 출력값이고,
상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고,
상기 i는 출력층의 노드에 대응하는 인덱스이고,
상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고,
상기
Figure 112022026989181-pat00073
는 하이퍼파라미터인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
7. The method of claim 6,
the learning unit
When performing the regular learning on the prediction network,
loss function
Figure 112022026989181-pat00070

using
remind
Figure 112022026989181-pat00071
is set to 0, and classification loss optimization is performed to correct the parameters of the prediction network so that the classification loss, which is the difference between the output value of the prediction network and the classification label, is minimized,
remind
Figure 112022026989181-pat00072
is set to 0.5, and the composite loss including the classification loss, which is the difference between the output value of the prediction network and the classification label, and the concealment loss, which is the difference between the output value and the hidden label of the last layer of the hidden layer of the prediction network, is predicted to be minimal. performing complex loss optimization that modifies the parameters of the network;
The E(y, f) is a loss function,
Wherein yi is the output value of the output layer,
Wherein ei is a classification label corresponding to the output value of the output layer,
The gij is the output value of the last layer of the hidden layer,
Wherein cij is a hidden label corresponding to the output value of the last layer of the hidden layer,
wherein i is an index corresponding to the node of the output layer,
Where j is an index corresponding to the node of the last layer of the hidden layer,
remind
Figure 112022026989181-pat00073
is a hyperparameter
A device for diagnosing the condition of skeletal muscle.
골격근의 상태를 진단하기 위한 장치에 있어서,
골격근에 대한 초음파 영상이 입력되면, 서로 다른 2개의 형식에 따라 상기 초음파 영상을 가공하고, 상기 초음파 영상으로부터 복수의 방사성 특징을 추출하는 전처리부; 및
복수의 계층을 포함하는 분류망을 통해 상기 2개의 형식 중 어느 하나의 형식으로 가공된 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근의 근육에코발생도의 등급을 확률로 산출하고, 산출된 확률에 따라 상기 골격근의 근육에코발생도의 등급을 분류하는 분류부;
복수의 계층을 포함하는 예측망을 통해 상기 2개의 형식 중 다른 하나의 형식으로 가공된 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 예측하는 예측부; 및
복수의 계층을 포함하는 식별망을 통해 상기 복수의 방사성 특징에 대한 연산을 수행하여 상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 식별하는 식별부;
를 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
In the apparatus for diagnosing the condition of skeletal muscle,
When an ultrasound image of skeletal muscle is input, the preprocessor processes the ultrasound image according to two different formats and extracts a plurality of radioactive features from the ultrasound image; and
Through a classification network including a plurality of layers, a pixel unit operation is performed on an ultrasound image processed in any one of the two formats, and the grade of the muscle echogenicity of the skeletal muscle is calculated as a probability, and the calculated probability a classification unit for classifying the degree of muscle echogenicity of the skeletal muscle according to the
Through a prediction network including a plurality of layers, a pixel unit operation is performed on the ultrasound image processed in the other one of the two formats to calculate with a probability whether sarcopenia is in progress in the skeletal muscle, and the calculated a predictor for predicting whether sarcopenia progresses according to a probability; and
Identification for calculating whether sarcopenia is in progress in the skeletal muscle by performing calculations on the plurality of radioactive features through an identification network including a plurality of layers, and identifying whether sarcopenia is in progress according to the calculated probability wealth;
characterized in that it comprises
A device for diagnosing the condition of skeletal muscle.
제8항에 있어서,
상기 예측망은
입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하며,
상기 예측부가 근육감소증을 예측하면,
상기 예측망의 상기 컨벌루션층의 출력값 및 경사도(gradient)를 비교하여 상기 초음파 영상에서 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시하는 시각화 영상을 생성하는
시각화부;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
9. The method of claim 8,
The forecasting network is
an input layer, at least a pair of alternating convolutional and pooling layers, at least one fully connected layer and an output layer;
If the prediction unit predicts sarcopenia,
Comparing the output value and the gradient of the convolutional layer of the prediction network to generate a visualization image displaying a region causing sarcopenia in the ultrasound image by distinguishing it from other regions
visualization unit;
characterized in that it further comprises
A device for diagnosing the condition of skeletal muscle.
제9항에 있어서,
상기 시각화부는
상기 예측망의 연산의 결과를 수집한 후,
수학식
Figure 112020066801883-pat00074

에 따라 뉴런 중요도 가중치를 산출하고,
수학식
Figure 112020066801883-pat00075

에 따라 상기 시각화 영상을 생성하며,
상기
Figure 112020066801883-pat00076
는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고,
상기
Figure 112020066801883-pat00077
는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고,
상기
Figure 112020066801883-pat00078
Figure 112020066801883-pat00079
Figure 112020066801883-pat00080
에 대해 가지는 경사도(gradient)이고,
상기
Figure 112020066801883-pat00081
는 상기 뉴런 중요도 가중치이고,
상기
Figure 112020066801883-pat00082
는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
10. The method of claim 9,
The visualization unit
After collecting the results of the calculation of the prediction network,
formula
Figure 112020066801883-pat00074

Calculate the neuron importance weight according to
formula
Figure 112020066801883-pat00075

to generate the visualization image according to
remind
Figure 112020066801883-pat00076
is the c-th value of the output before performing the convolution operation in the convolutional layer and performing the operation by the activation function,
remind
Figure 112020066801883-pat00077
is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer,
remind
Figure 112020066801883-pat00078
Is
Figure 112020066801883-pat00079
go
Figure 112020066801883-pat00080
is the gradient with respect to
remind
Figure 112020066801883-pat00081
is the neuron importance weight,
remind
Figure 112020066801883-pat00082
is the value of the coordinates (i, j) of the visualization image
A device for diagnosing the condition of skeletal muscle.
제8항에 있어서,
상기 예측망의 출력층의 출력노드가 소프트맥스 연산을 수행하는 상기 예측망의 초기 모델을 마련하고, 상기 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행하며,
상기 초기 학습이 완료되면, 상기 분류망 및 상기 예측망의 은닉층의 마지막층을 초기화된 완전연결층으로 교체하고,
상기 출력층의 상기 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 상기 예측망의 정규 모델을 마련하고, 상기 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행하는 학습부;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
9. The method of claim 8,
An output node of the output layer of the prediction network prepares an initial model of the prediction network that performs softmax operation, and performs initial learning using an image stored in a public database for the initial model as training data,
When the initial learning is completed, the last layer of the hidden layer of the classification network and the prediction network is replaced with an initialized fully connected layer,
A regular model of the prediction network is prepared by replacing an output node that performs the softmax operation of the output layer with an output node that performs a sigmoid operation, and regular learning is performed using an ultrasound image as training data for the regular model. a learning unit to perform;
characterized in that it further comprises
A device for diagnosing the condition of skeletal muscle.
제11항에 있어서,
상기 학습부는
상기 예측망에 대해 상기 정규 학습을 수행할 때,
손실함수
Figure 112022026989181-pat00083

를 이용하여
상기
Figure 112022026989181-pat00084
을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하고,
상기
Figure 112022026989181-pat00085
을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행하고,
상기 E(y, f)는 손실함수이고,
상기 yi는 출력층의 출력값이고,
상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고,
상기 gij는 은닉층의 마지막층의 출력값이고,
상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고,
상기 i는 출력층의 노드에 대응하는 인덱스이고,
상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고,
상기
Figure 112022026989181-pat00086
는 하이퍼파라미터인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치.
12. The method of claim 11,
the learning unit
When performing the regular learning on the prediction network,
loss function
Figure 112022026989181-pat00083

using
remind
Figure 112022026989181-pat00084
is set to 0, and classification loss optimization is performed to correct the parameters of the prediction network so that the classification loss, which is the difference between the output value of the prediction network and the classification label, is minimized,
remind
Figure 112022026989181-pat00085
is set to 0.5, and the composite loss including the classification loss, which is the difference between the output value of the prediction network and the classification label, and the concealment loss, which is the difference between the output value and the hidden label of the last layer of the hidden layer of the prediction network, is predicted to be minimal. performing complex loss optimization that modifies the parameters of the network;
The E(y, f) is a loss function,
Wherein yi is the output value of the output layer,
Wherein ei is a classification label corresponding to the output value of the output layer,
The gij is the output value of the last layer of the hidden layer,
Wherein cij is a hidden label corresponding to the output value of the last layer of the hidden layer,
wherein i is an index corresponding to the node of the output layer,
Where j is an index corresponding to the node of the last layer of the hidden layer,
remind
Figure 112022026989181-pat00086
is a hyperparameter
A device for diagnosing the condition of skeletal muscle.
골격근의 상태를 진단하기 위한 방법에 있어서,
예측부가 골격근에 대한 초음파 영상이 입력되면, 입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하는 예측망을 통해 상기 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하는 단계; 및
상기 예측부가 산출된 확률에 따라 근육감소증 진행 여부를 예측하는 단계;
를 포함하며,
분류부가 상기 골격근에 대한 초음파 영상이 입력되면, 입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하는 분류망을 통해 상기 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근이 근육에코발생도의 각 등급에 속할 확률을 산출하는 단계;
상기 분류부가 상기 산출된 확률에 따라 상기 골격근의 근육에코발생도의 등급을 분류하는 단계;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 방법.
In the method for diagnosing the condition of skeletal muscle,
When the prediction unit receives an ultrasound image of skeletal muscle, the ultrasound image is processed in pixel units through a prediction network including an input layer, at least a pair of alternately repeated convolutional and pooling layers, at least one fully connected layer, and an output layer. calculating with a probability whether sarcopenia is in progress in the skeletal muscle by performing an operation; and
predicting whether sarcopenia progresses according to the probability calculated by the prediction unit;
includes,
When the classification unit receives the ultrasound image of the skeletal muscle, the ultrasound image is processed through a classification network including an input layer, at least a pair of alternately repeated convolutional and pooling layers, at least one fully connected layer, and an output layer. calculating a probability that the skeletal muscle belongs to each class of muscle echogenicity by performing a unit operation;
classifying, by the classification unit, the degree of muscle echogenicity of the skeletal muscle according to the calculated probability;
characterized in that it further comprises
A method for diagnosing a condition of skeletal muscle.
제13항에 있어서,
상기 예측부가 근육감소증을 예측하면, 시각화부가 상기 예측망의 상기 컨벌루션층의 출력값 및 경사도(gradient)를 비교하여 상기 초음파 영상에서 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시하는 시각화 영상을 생성하는 단계;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 방법.
14. The method of claim 13,
When the prediction unit predicts sarcopenia, the visualization unit compares the output value and the gradient of the convolutional layer of the prediction network to distinguish and display a region causing sarcopenia in the ultrasound image from other regions. generating;
characterized in that it further comprises
A method for diagnosing a condition of skeletal muscle.
제14항에 있어서,
상기 시각화 영상을 생성하는 단계는
상기 시각화부가 상기 예측망의 연산의 결과를 수집하는 단계;
상기 시각화부가
수학식
Figure 112020066801883-pat00087

에 따라 뉴런 중요도 가중치를 산출하는 단계;
상기 시각화부가
수학식
Figure 112020066801883-pat00088

에 따라 상기 시각화 영상을 생성하는 단계;
를 포함하며,
상기
Figure 112020066801883-pat00089
는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고,
상기
Figure 112020066801883-pat00090
는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고,
상기
Figure 112020066801883-pat00091
는 상기
Figure 112020066801883-pat00092
가 상기
Figure 112020066801883-pat00093
에 대해 가지는 경사도(gradient)이고,
상기
Figure 112020066801883-pat00094
는 상기 뉴런 중요도 가중치이고,
상기
Figure 112020066801883-pat00095
는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 방법.
15. The method of claim 14,
The step of generating the visualization image is
collecting, by the visualization unit, a result of the calculation of the prediction network;
the visualization unit
formula
Figure 112020066801883-pat00087

calculating a neuron importance weight according to
the visualization unit
formula
Figure 112020066801883-pat00088

generating the visualization image according to
includes,
remind
Figure 112020066801883-pat00089
is the c-th value of the output before performing the convolution operation in the convolutional layer and performing the operation by the activation function,
remind
Figure 112020066801883-pat00090
is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer,
remind
Figure 112020066801883-pat00091
is said
Figure 112020066801883-pat00092
is reminded
Figure 112020066801883-pat00093
is the gradient with respect to
remind
Figure 112020066801883-pat00094
is the neuron importance weight,
remind
Figure 112020066801883-pat00095
is the value of the coordinates (i, j) of the visualization image
A method for diagnosing a condition of skeletal muscle.
제13항에 있어서,
식별부가 상기 골격근에 대한 초음파 영상으로부터 추출된 복수의 방사성 특징이 입력되면, 입력층, 은닉층 및 출력층을 포함하는 식별망을 통해 상기 복수의 방사성 특징에 대한 연산을 수행하여 상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하는 단계; 및
상기 식별부가 산출된 확률에 따라 근육감소증 진행 여부를 식별하는 단계;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 방법.
14. The method of claim 13,
When the identification unit receives a plurality of radioactive features extracted from the ultrasound image of the skeletal muscle, sarcopenia progresses in the skeletal muscle by performing an operation on the plurality of radioactive features through an identification network including an input layer, a hidden layer, and an output layer. calculating whether or not it is being performed with a probability; and
identifying whether sarcopenia progresses according to the probability calculated by the identification unit;
characterized in that it further comprises
A method for diagnosing a condition of skeletal muscle.
삭제delete 제13항에 있어서,
상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하는 단계 전,
학습부가 상기 예측망의 상기 출력층이 소프트맥스 연산을 수행하는 출력노드로 이루어진 상기 예측망의 초기 모델을 마련하는 단계; 및
상기 학습부가 상기 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행하는 단계;
상기 학습부가 상기 예측망의 마지막 완전연결층을 초기화된 완전연결층으로 교체하고, 상기 출력층의 상기 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 상기 예측망의 정규 모델을 마련하는 단계; 및
상기 학습부가 상기 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행하는 단계;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 방법.
14. The method of claim 13,
Before the step of calculating with a probability whether sarcopenia is in progress in the skeletal muscle,
preparing, by a learning unit, an initial model of the prediction network including an output node in which the output layer of the prediction network performs a softmax operation; and
performing, by the learning unit, initial learning with respect to the initial model by using an image stored in a public database as training data;
The learning unit replaces the last fully connected layer of the prediction network with an initialized fully connected layer, and replaces the output node performing the softmax operation of the output layer with an output node performing the sigmoid operation, and preparing a canonical model; and
performing, by the learning unit, regular learning on the regular model using an ultrasound image as training data;
characterized in that it further comprises
A method for diagnosing a condition of skeletal muscle.
제18항에 있어서,
상기 정규 학습을 수행하는 단계는
상기 학습부가
손실함수
Figure 112022026989181-pat00096

를 이용하여
상기
Figure 112022026989181-pat00097
을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하는 단계;
상기 학습부가 상기
Figure 112022026989181-pat00098
을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행하는 단계;
를 포함하며,
상기 E(y, f)는 손실함수이고,
상기 yi는 출력층의 출력값이고,
상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고,
상기 gij는 은닉층의 마지막층의 출력값이고,
상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고,
상기 i는 출력층의 노드에 대응하는 인덱스이고,
상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고,
상기
Figure 112022026989181-pat00099
는 하이퍼파라미터인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 방법.
19. The method of claim 18,
The step of performing the regular learning is
the learning department
loss function
Figure 112022026989181-pat00096

using
remind
Figure 112022026989181-pat00097
after setting to 0, performing classification loss optimization by modifying parameters of the prediction network so that the classification loss, which is the difference between the output value of the prediction network and the classification label, is minimized;
the learning unit
Figure 112022026989181-pat00098
is set to 0.5, and the composite loss including the classification loss, which is the difference between the output value of the prediction network and the classification label, and the concealment loss, which is the difference between the output value and the hidden label of the last layer of the hidden layer of the prediction network, is predicted to be minimal. performing complex loss optimization to modify parameters of the network;
includes,
The E(y, f) is a loss function,
Wherein yi is the output value of the output layer,
Wherein ei is a classification label corresponding to the output value of the output layer,
The gij is the output value of the last layer of the hidden layer,
Wherein cij is a hidden label corresponding to the output value of the last layer of the hidden layer,
wherein i is an index corresponding to the node of the output layer,
Where j is an index corresponding to the node of the last layer of the hidden layer,
remind
Figure 112022026989181-pat00099
is a hyperparameter
A method for diagnosing a condition of skeletal muscle.
KR1020200079181A 2020-06-29 2020-06-29 Apparatus and method for diagnosing skeletal muscle condition from ultrasound images using deep neural network KR102436035B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020200079181A KR102436035B1 (en) 2020-06-29 2020-06-29 Apparatus and method for diagnosing skeletal muscle condition from ultrasound images using deep neural network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020200079181A KR102436035B1 (en) 2020-06-29 2020-06-29 Apparatus and method for diagnosing skeletal muscle condition from ultrasound images using deep neural network

Publications (2)

Publication Number Publication Date
KR20220001177A KR20220001177A (en) 2022-01-05
KR102436035B1 true KR102436035B1 (en) 2022-08-24

Family

ID=79348811

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020200079181A KR102436035B1 (en) 2020-06-29 2020-06-29 Apparatus and method for diagnosing skeletal muscle condition from ultrasound images using deep neural network

Country Status (1)

Country Link
KR (1) KR102436035B1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2808193C1 (en) * 2023-02-14 2023-11-24 Федеральное государственное бюджетное образовательное учреждение высшего образования "Смоленский государственный медицинский университет" министерства здравоохранения Российской Федерации Method for diagnosis of sarcopenia in elderly and senile patients using quantitative assessment of muscle tissue according to ultrasound study of rectum abdominal muscle
KR20240064107A (en) 2022-11-04 2024-05-13 경북대학교 산학협력단 Apparatus and method of classifying rotator cuf tears based on deep learning for ultrasonic wave images

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116721088A (en) * 2023-06-20 2023-09-08 北京积水潭医院 Artificial intelligent identification method and device for bone and meat tumor image based on deep learning

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180146947A1 (en) 2016-11-30 2018-05-31 MuscleSound, LLC Non-Invasive Determination of Muscle Tissue Quality and Intramuscular Fat
US20200129140A1 (en) 2017-06-21 2020-04-30 Charité - Universitätsmedizin Berlin System, method, and computer program product for determining cortical bone characteristics
JP2020130596A (en) 2019-02-19 2020-08-31 株式会社Cesデカルト Ultrasonic diagnostic apparatus, ultrasonic diagnostic program, and ultrasonic echo image analysis method
WO2021048776A1 (en) 2019-09-11 2021-03-18 Universita' Degli Studi Di Genova Method for evaluating the presence of sarcopenia by magnetic resonance or echography or computed tomography

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101870121B1 (en) * 2015-10-16 2018-06-25 재단법인 아산사회복지재단 System, method and program for analyzing blood flow by deep neural network
KR101961215B1 (en) * 2017-08-25 2019-03-27 남정림 Apparatus and method for estimating bone age based on image processing and machine learning
KR102140437B1 (en) * 2018-02-08 2020-08-03 주식회사 셀바스에이아이 Method for predicting cancer suspicious area and device for predicting cancer suspicious area using the same
KR102030533B1 (en) 2018-03-27 2019-10-10 울산대학교 산학협력단 Image processing apparatus for adopting human body morphometric based on artificial neural network for sarcopenia and image processing method using the same

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180146947A1 (en) 2016-11-30 2018-05-31 MuscleSound, LLC Non-Invasive Determination of Muscle Tissue Quality and Intramuscular Fat
US20200129140A1 (en) 2017-06-21 2020-04-30 Charité - Universitätsmedizin Berlin System, method, and computer program product for determining cortical bone characteristics
JP2020130596A (en) 2019-02-19 2020-08-31 株式会社Cesデカルト Ultrasonic diagnostic apparatus, ultrasonic diagnostic program, and ultrasonic echo image analysis method
WO2021048776A1 (en) 2019-09-11 2021-03-18 Universita' Degli Studi Di Genova Method for evaluating the presence of sarcopenia by magnetic resonance or echography or computed tomography

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Artificial intelligence in musculoskeletal ultrasound imaing", 대한초음파의학회, Vol.40, No.1, p.30~44, 2021
"The role of ultrasound as a diagnostic tool for sarcopenia", The Journal of Frailty & Aging, Vol. 7, No. 4, 2018

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20240064107A (en) 2022-11-04 2024-05-13 경북대학교 산학협력단 Apparatus and method of classifying rotator cuf tears based on deep learning for ultrasonic wave images
RU2808193C1 (en) * 2023-02-14 2023-11-24 Федеральное государственное бюджетное образовательное учреждение высшего образования "Смоленский государственный медицинский университет" министерства здравоохранения Российской Федерации Method for diagnosis of sarcopenia in elderly and senile patients using quantitative assessment of muscle tissue according to ultrasound study of rectum abdominal muscle

Also Published As

Publication number Publication date
KR20220001177A (en) 2022-01-05

Similar Documents

Publication Publication Date Title
US10706333B2 (en) Medical image analysis method, medical image analysis system and storage medium
US20190139641A1 (en) Artificial intelligence for physiological quantification in medical imaging
CN110070540B (en) Image generation method and device, computer equipment and storage medium
CN110475505A (en) Utilize the automatic segmentation of full convolutional network
US11350888B2 (en) Risk prediction for sudden cardiac death from image derived cardiac motion and structure features
JP6885517B1 (en) Diagnostic support device and model generation device
US10522253B2 (en) Machine-learnt prediction of uncertainty or sensitivity for hemodynamic quantification in medical imaging
JP2013537445A (en) Visualization adapted for direct use by physicians
KR102436035B1 (en) Apparatus and method for diagnosing skeletal muscle condition from ultrasound images using deep neural network
CN114494296A (en) Brain glioma segmentation method and system based on fusion of Unet and Transformer
CN110991254B (en) Ultrasonic image video classification prediction method and system
US20230230241A1 (en) System and method for detecting lung abnormalities
CN113554665A (en) Blood vessel segmentation method and device
CN114897780A (en) MIP sequence-based mesenteric artery blood vessel reconstruction method
CN113850753A (en) Medical image information calculation method and device, edge calculation equipment and storage medium
CN112633416A (en) Brain CT image classification method fusing multi-scale superpixels
CN112488971A (en) Medical image fusion method for generating countermeasure network based on spatial attention mechanism and depth convolution
US10909676B2 (en) Method and system for clinical decision support with local and remote analytics
CN113313728A (en) Intracranial artery segmentation method and system
EP3477551B1 (en) Machine-learnt prediction of uncertainty or sensitivity for hemodynamic quantification in medical imaging
CN115841476A (en) Method, device, equipment and medium for predicting life cycle of liver cancer patient
CN115018780A (en) Thyroid nodule segmentation method fusing global reasoning and MLP framework
CN115294023A (en) Liver tumor automatic segmentation method and device
CN116725560B (en) Low-dose DSA angiography method based on spatial diffusion model
KR102464422B1 (en) Apparatus, method for providing information needed, and recording medium for diagnosing lung damage

Legal Events

Date Code Title Description
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right