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 PDFInfo
- 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
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/042—Knowledge-based neural networks; Logical representations of neural networks
-
- G06N3/0427—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20084—Artificial 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
본 발명은 골격근 상태 진단 기술에 관한 것으로, 보다 상세하게는, 심층신경망(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).
본 발명의 목적은 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치 및 이를 위한 방법을 제공함에 있다. 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.
상기 시각화부는 상기 예측망의 연산의 결과를 수집한 후, 수학식 에 따라 뉴런 중요도 가중치를 산출하고, 수학식 에 따라 상기 시각화 영상을 생성하며, 상기 는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고, 상기 는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고, 상기 는 가 에 대해 가지는 경사도(gradient)이고, 상기 는 상기 뉴런 중요도 가중치이고, 상기 는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 한다. After the visualization unit collects the results of the calculation of the prediction network, Calculate the neuron importance weight according to the formula to generate the visualization image according to the 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, is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer, Is go is the gradient with respect to , and is the neuron importance weight, and 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.
상기 학습부는 상기 예측망에 대해 상기 정규 학습을 수행할 때, 손실함수 를 이용하여 상기 을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하고, 상기 을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행하고, 상기 E(y, f)는 손실함수이고, 상기 yi는 출력층의 출력값이고, 상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고, 상기 gij는 은닉층의 마지막층의 출력값이고, 상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고, 상기 i는 출력층의 노드에 대응하는 인덱스이고, 상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고, 상기 는 하이퍼파라미터인 것을 특징으로 한다. When the learning unit performs the regular learning on the prediction network, a loss function above using 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 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 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.
상기 시각화부는 상기 예측망의 연산의 결과를 수집한 후, 수학식 에 따라 뉴런 중요도 가중치를 산출하고, 수학식 에 따라 상기 시각화 영상을 생성하며, 상기 는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고, 상기 는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고, 상기 는 가 에 대해 가지는 경사도(gradient)이고, 상기 는 상기 뉴런 중요도 가중치이고, 상기 는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 한다. After the visualization unit collects the results of the calculation of the prediction network, Calculate the neuron importance weight according to the formula to generate the visualization image according to the 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, is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer, Is go is the gradient with respect to , and is the neuron importance weight, and 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.
상기 학습부는 상기 예측망에 대해 상기 정규 학습을 수행할 때, 손실함수 를 이용하여 상기 을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하고, 상기 을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행한다. When the learning unit performs the regular learning on the prediction network, a loss function above using 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 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는 은닉층의 마지막층의 노드에 대응하는 인덱스이고, 상기 는 하이퍼파라미터인 것을 특징으로 한다. 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 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.
상기 시각화 영상을 생성하는 단계는 상기 시각화부가 상기 예측망의 연산의 결과를 수집하는 단계와, 상기 시각화부가 수학식 에 따라 뉴런 중요도 가중치를 산출하는 단계와, 상기 시각화부가 수학식 에 따라 상기 시각화 영상을 생성하는 단계를 포함한다. 여기서, 상기 는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고, 상기 는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고, 상기 는 상기 가 상기 에 대해 가지는 경사도(gradient)이고, 상기 는 상기 뉴런 중요도 가중치이고, 상기 는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 한다. The generating of the visualization image includes: the visualization unit collecting the results of the calculation of the prediction network; calculating the neuron importance weight according to and generating the visualization image according to the Here, the 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, is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer, is said is reminded is the gradient with respect to , and is the neuron importance weight, and 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.
상기 정규 학습을 수행하는 단계는 상기 학습부가 손실함수 를 이용하여 상기 을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하는 단계와, 상기 학습부가 상기 을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행하는 단계를 포함한다. 여기서, 상기 E(y, f)는 손실함수이고, 상기 yi는 출력층의 출력값이고, 상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고, 상기 gij는 은닉층의 마지막층의 출력값이고, 상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고, 상기 i는 출력층의 노드에 대응하는 인덱스이고, 상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고, 상기 는 하이퍼파라미터인 것을 특징으로 한다. In the step of performing the regular learning, the learning unit is a loss function above using 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; 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 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
진단영상장치(20)는 기본적으로, 초음파(USG: ultrasonography) 영상을 촬영하기 위한 것이다. 본 발명의 실시예에서 진단영상장치(20)는 B 모드 초음파 영상 및 SWE 초음파 영상 양자 모두를 촬영할 수 있는 하나의 장치로 설명될 것이지만, 진단영상장치(20)는 B 모드 초음파 영상을 촬영하는 장치와 SWE 초음파 영상을 촬영하는 장치로 구분되어 2개의 장치로 구현될 수 있다. The
진단장치(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
인터페이스부(11)는 진단영상장치(20)로부터 초음파 영상을 수신하기 위한 것이다. 인터페이스부(11)는 필요에 따라 무선 통신, 유선 통신 등을 통해 진단영상장치(20)와 연결되어 초음파 영상을 수신할 수 있다. The
입력부(12)는 진단장치(10)의 각 종 기능, 동작 등을 제어하기 위한 사용자의 키 조작을 입력받고 입력 신호를 생성하여 제어부(15)에 전달한다. 입력부(12)는 키보드, 마우스 등을 예시할 수 있다. 입력부(12)는 전원 on/off를 위한 전원키, 문자키, 숫자키, 방향키 등을 포함할 수 있다. 입력부(12)의 기능은 표시부(13)가 터치스크린으로 구현된 경우, 표시부(13)에서 이루어질 수 있으며, 표시부(13)만으로 모든 기능을 수행할 수 있는 경우, 입력부(12)는 생략될 수도 있다. The
표시부(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
저장부(14)는 진단장치(10)의 동작에 필요한 각 종 데이터, 애플리케이션, 진단장치(10)의 동작에 따라 구성되는 각 종 데이터를 저장하는 역할을 수행한다. 이러한 저장부(14)는 크게 프로그램 영역과 데이터 영역을 포함할 수 있다. 프로그램 영역은 진단장치(10)의 부팅(booting) 및 운영(operation)을 위한 운영체제(OS, Operating System), 본 발명의 실시예에 따른 애플리케이션 등을 저장할 수 있다. 데이터 영역은 본 발명의 실시예에 따른 진단을 위한 각 종 데이터를 저장할 수 있다. 저장부(14)에 저장되는 각 종 데이터는 사용자의 조작에 따라, 삭제, 변경, 추가될 수 있다. The
제어부(15)는 진단장치(10)의 전반적인 동작 및 진단장치(10)의 내부 블록들 간 신호 흐름을 제어하고, 데이터를 처리하는 데이터 처리 기능을 수행할 수 있다. 이러한 제어부(15)는 중앙처리장치(CPU: Central Processing Unit), 디지털신호처리기(DSP: Digital Signal Processor) 등이 될 수 있다. 또한, 제어부(15)는 추가로 이미지 프로세서(Image processor) 혹은 GPU(Graphic Processing Unit)를 더 구비할 수 있다. 이러한 제어부(15)의 구체적인 동작에 대해서는 아래에서 더 상세하게 설명될 것이다. The
다음으로, 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 장치의 구성에 대해 보다 상세하게 설명하기로 한다. 도 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
학습부(100)는 본 발명의 실시예에 따른 심층신경망(DNN: Deep Neural Network)을 학습시키기 위한 일련의 프로세스를 수행한다. 한편, 분류부(300), 예측부(400) 및 식별부(600)는 각각 심층신경망으로 분류망(310), 예측망(410) 및 식별망(610)을 포함한다. 학습부(100)는 이러한 분류망(310), 예측망(410) 및 식별망(610)을 학습시킨다. 이러한 학습 프로세스에 대해서는 아래에서 더 상세하게 설명하기로 한다. The
도 3을 참조하면, 전처리부(200)는 진단영상장치(20)로부터 인터페이스부(11)를 통해 초음파 영상을 수신하면, 분류부(300)의 분류망(310), 예측부(400)의 예측망(410) 및 식별부(600)의 식별망(610) 중 적어도 하나에 부합하는 형식에 맞춰 수신된 초음파 영상을 가공하거나, 정보를 추출하여 가공된 초음파 영상 혹은 정보를 제공한다. 이를 위하여 전처리부(200)는 증강부(210), 영상처리부(220) 및 특징처리부(230)를 포함한다. Referring to FIG. 3 , when the
증강부(210)는 본 발명의 실시예에 따른 학습 데이터의 수를 증강시키기 위한 것이다. 일 실시예에 따르면, 증강부(210)는 학습용으로 마련된 초음파 영상이 입력되면, 플립핑(flipping), 시프트(shift), 쉬어(shear), 확대, 축소 및 회전 등의 변형을 통해 학습 데이터의 수를 증강시킬 수 있다. The
영상처리부(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
특징처리부(230)는 초음파 영상에서 방사성 특징을 추출하기 위한 것이다. 먼저, 특징처리부(230)는 B-모드 초음파 영상에서 모든 복수의 후보 방사성 특징(Radiomics features)을 추출한다. 그런 다음, 특징처리부(230)는 추출된 복수의 후보 방사성 특징으로부터 LASSO(least absolute shrinkage and selection operator) 알고리즘을 통해 루소 계수(LASSO regression coefficients)가 높은 순으로 식별망의 입력층(PIL)의 크기에 맞춰, 소정 개수, 예컨대, 30개의 방사성 특징을 도출하였다. 이러한 방사성 특징은 다음의 표 1과 같다. The
다음으로, 도 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
컨볼루션층(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
이어서, 풀링층(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
노드(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
설명되지 않은 파라미터 중 b는 임계치이며, 이러한 임계치는 수학식 1에서 의 값이 임계치 보다 작을 때 해당 노드가 활성화되지 않도록 하는 역할을 한다. 예를 들면, 노드(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 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
다음으로, 도 7 및 도 8을 참조하여 예측부(400)에 대해 설명한다. 전술한 바와 같이, 예측부(400)는 예측망(410)을 포함하며, 예측망(410)은 심층신경망 중 심층컨벌루션신경망(DCNN)이 될 수 있다. 본 발명의 예측망(410)은 복수의 계층을 포함한다. 즉, 예측망(410)은 입력층(EIL), 교번으로 반복되는 적어도 한 쌍의 컨벌루션층(ECL)과 풀링층(EPL), 적어도 하나의 완전연결층(EFL1, EFL2) 및 출력층(EOL)을 포함한다. 전처리부(200)의 영상처리부(220)로부터 골격근에 대한 초음파 영상을 수신하면, 예측망(410)을 통해 초음파 영상에 대해 픽셀 단위 연산을 수행하여 초음파 영상의 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 예측한다. Next, the
일 실시예에 따르면, 예측망(410)은 도 7에 도시된 바와 같이, 순차로 입력층(EIL), 컨볼루션층(ECL), 풀링층(EPL), 제1 완전연결계층(EFL1), 제2 완전연결계층(EFL2) 및 출력층(EOL)을 포함한다. According to an embodiment, the
컨볼루션층(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
이어서, 풀링층(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
한편, 도 2 및 도 7을 참조하면, 시각화부(500)는 예측망(410)을 이용하여 시각화 영상을 생성하기 위한 것이다. 시각화 영상은 입력되는 초음파 영상에서 근육감소증의 원인이 되는 영역을 도출하고, 이를 다른 영역과 구별되게 표시한 영상을 의미한다. 이러한 영상은 CAM(Class Activation Map), Grad-CAM 등이 될 수 있다. 이에 따라, 예측부(400)가 예측망(410)의 출력인 확률에 따라 근육감소증을 예측하면, 예측망(410)의 컨벌루션층(ECL)의 출력값 및 경사도(gradient)를 비교하여 초음파 영상에서 근육감소증의 원인이 되는 영역을 도출하고, 초음파 영상에서 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시한다. 이로써, 시각화 영상이 생성된다. 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시할 때, 색의 종류 및 밝기를 이용하여 표시할 수 있다. 이러한 방법에 대해 보다 상세하게 설명하기로 한다. Meanwhile, referring to FIGS. 2 and 7 , the
시각화부(500)는 먼저, 예측망(410)의 연산의 결과를 수집한 후, 다음의 수학식 2에 따라 뉴런 중요도 가중치(neuron importance weight)를 산출한다. The
여기서, 는 뉴런 중요도 가중치이고, 는 컨벌루션층(ECL)에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고, 는 컨벌루션층(ECL)의 k번째 특징지도의 좌표(i, j)의 값을 의미한다. 특히, 는 가 에 대해 가지는 경사도(gradient)를 의미한다. here, is the neuron importance weight, 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, is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer (ECL). Especially, Is go It means the gradient with respect to .
이어서, 시각화부(500)는 다음의 수학식 3에 따라 시각화 영상을 생성한다. Next, the
여기서, 는 시각화 영상의 좌표(i, j)의 값을 나타낸다. 또한, 는 뉴런 중요도 가중치이고, 는 컨벌루션층(ECL)의 k번째 특징영상(FM)의 좌표(i, j)의 값을 의미한다. 이와 같이, 시각화부(500)는 를 산출한 후, 값에 따라 색을 적용함으로써, 시각화 영상을 생성할 수 있다. here, denotes the value of the coordinates (i, j) of the visualization image. In addition, is the neuron importance weight, is the value of the coordinates (i, j) of the k-th feature image FM of the convolutional layer ECL. In this way, the
도 9를 참조하여 식별부(600)에 대해 설명하기로 한다. 전술한 바와 같이, 식별부(600)는 식별망(610)을 포함하며, 식별망(610)은 심층신경망 중 다층퍼셉트론(MLP: Multilayer Perceptron)이 될 수 있다. 식별망(610)은 복수의 계층(PIL, PHL, POL)을 포함한다. 이러한 복수의 계층은 입력층(PIL), 복수의 은닉 계층(PHL: HL1 내지 HLk) 및 출력층(POL)을 포함한다. The
또한, 복수의 계층(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
보다 구체적으로 설명하면, 복수의 방사성 특징(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
이에 따라, 도 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
다시 도 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
다음으로, 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 방법에 대해서 설명하기로 한다. 골격근의 상태를 진단하기 위해서는 심층신경망(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
도 10을 참조하면, 학습부(100)는 S110 단계에서 근육감소증 여부가 알려진 환자의 초음파 영상으로부터 추출된 방사성 특징을 학습 데이터, 즉, 학습용 방사성 특징으로 마련한다. 그런 다음, 학습부(100)는 S120 단계에서 학습용 방사성 특징에 알려진 근육감소증 여부에 따라 레이블을 설정한다. 이때, 근육감소증이 있는 환자의 초음파 영상으로부터 추출된 방사성 특징에는 출력노드(P, N)에 대응하여 레이블 [1, 0]을 부여한다. 반면, 근육감소증이 없는 환자의 초음파 영상으로부터 추출된 방사성 특징에는 출력노드(P, N)에 대응하여 레이블 [0, 1]을 부여한다. Referring to FIG. 10 , the
다음으로, 학습부(100)는 S130 단계에서 해당하는 학습용 방사성 특징을 식별망(610)에 입력한다. 그러면, 식별망(610)은 S140 단계에서 입력된 학습용 방사성 특징에 대해 복수의 계층(PIL, PHL, POL)의 가중치가 적용되는 복수의 연산을 수행하여 골격근에 근육감소증(Sarcopenia)이 진행되고 있는지 여부를 확률로 나타내는 출력값을 산출한다. Next, the
그러면, 학습부(100)는 S150 단계에서 손실함수를 통해 출력값과 레이블의 차이인 이진 교차 엔트로피 손실이 최소가 되도록 식별망(610)의 파라미터, 즉, 가중치(w), 임계치(b) 등을 수정하는 최적화를 수행할 수 있다. Then, the
전술한 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
그런 다음, 학습부(100)는 S220 단계에서 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행한다. 여기서, 공개 데이터베이스에 저장된 영상은 이미지넷(ImageNet)의 데이터베이스의 이미지를 이용할 수 있다. 초기 학습은 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 레이블을 설정하고, 학습 데이터에 대한 초기 모델의 출력값을 산출하고, 손실함수에 따라 초기 모델의 출력값과 레이블과의 차이인 손실, 예컨대, 이진교차엔트로피 손실이 최소가 되도록 초기 모델의 파라미터를 수정하는 절차이다. Then, the
초기 학습이 완료되면, 학습부(100)는 S230 단계에서 초기 모델의 상기 예측망의 은닉층의 마지막층인 마지막 완전연결층을 초기화된 완전연결층(예컨대, GFL2/EFL2)으로 교체하고, 출력층의 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 초기 모델을 정규 모델로 변환한다. When the initial learning is completed, the
이어서, 학습부(100)는 S240 단계에서 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행한다. 정규 학습은 초음파 영상, 즉, B 모드 초음파 영상 혹은 SWE 초음파 영상으로부터 전처리부(100)를 통해 학습용 초음파 영상을 생성하고, 학습용 초음파 영상에 대해 레이블을 설정한 후, 학습용 초음파 영상에 대해 정규 모델의 출력값을 산출하고, 손실함수에 따라 초기 모델의 출력값과 레이블과의 차이인 손실이 최소가 되도록 정규 모델의 파라미터를 수정하는 절차이다. Next, the
그러면, 전술한 정규 모델에 대한 학습 방법에 대해 보다 상세하게 설명하기로 한다. 도 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
수학식 4에서 E(y, f)는 손실함수를 나타낸다. yi는 출력층의 출력값이고, ei는 출력층의 출력값에 대응하는 분류 레이블이다. gij는 은닉층의 마지막층의 출력값이고, cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이다. i는 출력층의 노드에 대응하는 인덱스이고, j는 은닉층의 마지막층의 노드에 대응하는 인덱스이다. 그리고 는 하이퍼파라미터를 나타낸다. 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 is a hyperparameter.
학습부(100)는 S320 단계에서 분류 손실 최적화를 위한 설정을 수행한다. 이때, 학습부(100)는 하이퍼파라미터 및 분류 레이블을 설정한다. 이러한 S320 단계에서 학습부(100)는 하이퍼파라미터 는 0으로 설정한다. 이에 따라, 분류 손실을 최적화하는 학습에서 수학식 4의 손실함수는 다음의 수학식 5와 같다. The
또한, 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
전술한 바와 같이 하이퍼파라미터 및 분류 레이블이 설정되면, 학습부(100)는 S330 단계에서 분류 손실을 최적화하는 학습을 수행한다. 이때, 학습부(100)는 분류망(310) 및 식별망(410)의 정규 모델에 학습용 초음파 영상을 입력한다. 그러면, 정규 모델은 학습용 초음파 영상에 대해 복수의 계층 간 가중치가 적용되는 복수의 연산을 통해 출력값을 산출한다. 그러면, 학습부(100)는 하이퍼파라미터가 0으로 설정된 수학식 5와 같은 손실함수를 이용하여 앞서 설정된 분류 레이블과 출력층의 출력값의 차이를 나타내는 수학식 5의 손실함수의 분류 손실이 최소가 되도록 정규 모델의 파라미터를 수정하는 분류 손실 최적화를 수행한다. 이러한 분류 손실 최적화는 복수의 서로 다른 학습용 초음파 영상을 이용하여 반복하여 수행되며, 이러한 반복은 평가 지표를 통해 원하는 정확도에 도달할 때까지 이루어질 수 있다. When the hyperparameter and the classification label are set as described above, the
전술한 바와 같은 분류 손실 최적화에 따른 학습이 완료되면, 학습부(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
다음으로, 학습부(100)는 S350 단계에서 복합 손실 최적화를 위한 설정을 수행한다. 이때, 학습부(100)는 하이퍼파라미터, 분류 레이블 및 은닉 레이블을 설정한다. 이러한 S350 단계에서 학습부(100)는 하이퍼파라미터 는 0.5로 설정한다. 이에 따라, 분류 손실을 최적화하는 학습에서 수학식 4의 손실함수는 다음의 수학식 6과 같다. Next, the
또한, 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
특히, 은닉 레이블의 설정은, 앞서(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
이러한 복합 손실에 의한 학습이 완료되면, 벡터 공간 상의 복수의 은닉벡터는 도 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
다음으로, 전술한 바와 같은 방법으로 학습이 완료된 분류망(310), 예측망(410) 및 식별망(610)을 이용하여 골격근의 상태를 진단하는 방법에 대해서 설명하기로 한다. 도 14는 본 발명의 실시예에 따른 심층신경망을 이용하여 초음파 영상으로부터 골격근의 상태를 진단하기 위한 방법을 설명하기 위한 흐름도이다. Next, a method of diagnosing the state of skeletal muscle using the
도 14를 참조하면, 진단장치(10)의 전처리부(200)는 S410 단계에서 인터페이스부(11)를 통해 진단영상장치(20)로부터 초음파 영상을 입력받을 수 있다. 입력된 초음파 영상은 B 모드 초음파 영상 및 SWE 초음파 영상 중 어느 하나가 될 수 있다. Referring to FIG. 14 , the
그러면, 전처리부(200)는 S420 단계에서 초음파 영상에 대한 전처리를 수행한다. 이때, 전처리부(200)는 초음파 영상의 크기를 조정하거나, 소정의 크기로 잘라낸 후, 이미지 파일 포맷으로 변환할 수 있다. 또한, 전처리부(200)는 초음파 영상으로부터 방사상 특징(Radiomics features)을 추출할 수 있다. 전처리부(200)는 전처리가 완료되면, 전처리된 초음파 영상 혹은 방사상 특징(Radiomics features)을 분류부(300), 예측부(400) 및 식별부(600)에 제공한다. Then, the
분류부(300), 예측부(400), 시각화부(500) 및 식별부(600)는 S430 단계에서 골격근에 대한 평가를 수행하여 평가 결과를 도출한다. 즉, 분류부(300)는 분류망(310)을 통해 초음파 영상에 대해 픽셀 단위 연산을 수행하여 골격근이 근육에코발생도(Muscle echogenicity)의 각 등급에 속할 확률을 산출하고, 산출된 확률에 따라 골격근의 근육에코발생도의 등급을 분류한다. 또한, 예측부(400)는 예측망(410)을 통해 초음파 영상에 대해 픽셀 단위 연산을 수행하여 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 예측한다. 그리고 시각화부(500)는 예측망(410)의 컨벌루션층(ECL)의 출력값 및 경사도(gradient)를 비교하여 초음파 영상에서 근육감소증의 원인이 되는 영역을 다른 영역과 구분하여 표시하는 시각화 영상을 생성한다. 또한, 식별부(600)는 식별망(610)을 통해 복수의 방사성 특징에 대한 연산을 수행하여 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하고, 산출된 확률에 따라 근육감소증 진행 여부를 식별한다. The
다음으로, 리포트부(700)는 S440 단계에서 분류부(300), 예측부(400), 시각화부(500) 및 식별부(600) 각각으로부터 정보를 제공받고, 제공 받은 정보를 종합한 종합 리포트를 표시부(13)를 통해 표시할 수 있다. 여기서, 종합 리포트는 근육에코발생도의 고레벨(High grade)에 속할 확률 및 저레벨(Low grade)에 속할 확률, 그 확률에 따른 판단, 즉, "근육에코발생도: 고레벨(High grade)" 혹은 "근육에코발생도: 저레벨(Low grade)"와, 초음파 영상을 기초로 예측되는 근육감소증이 있을 확률 및 근육감소증이 없을 확률과, 그 확률에 따른 판단, "근육감소증: 양성" 혹은 "근육감소증: 음성"과, 시각화 영상과, 방사상 특징을 기초로 식별된 근육감소증이 있을 확률 및 근육감소증이 없을 확률과, 그 확률에 따른 판단, "근육감소증: 양성" 혹은 "근육감소증: 음성"을 포함한다. Next, the
한편, 전술한 본 발명의 실시예에 따른 방법은 다양한 컴퓨터수단을 통하여 판독 가능한 프로그램 형태로 구현되어 컴퓨터로 판독 가능한 기록매체에 기록될 수 있다. 여기서, 기록매체는 프로그램 명령, 데이터 파일, 데이터구조 등을 단독으로 또는 조합하여 포함할 수 있다. 기록매체에 기록되는 프로그램 명령은 본 발명을 위하여 특별히 설계되고 구성된 것들이거나 컴퓨터 소프트웨어 당업자에게 공지되어 사용 가능한 것일 수도 있다. 예컨대 기록매체는 하드 디스크, 플로피 디스크 및 자기 테이프와 같은 자기 매체(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.
상기 예측부가 근육감소증을 예측하면,
상기 예측망의 상기 컨벌루션층의 출력값 및 경사도(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.
상기 시각화부는
상기 예측망의 연산의 결과를 수집한 후,
수학식
에 따라 뉴런 중요도 가중치를 산출하고,
수학식
에 따라 상기 시각화 영상을 생성하며,
상기 는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고,
상기 는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고,
상기 는 가 에 대해 가지는 경사도(gradient)이고,
상기 는 상기 뉴런 중요도 가중치이고,
상기 는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치. 3. The method of claim 2,
The visualization unit
After collecting the results of the calculation of the prediction network,
formula
Calculate the neuron importance weight according to
formula
to generate the visualization image according to
remind 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 is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer,
remind Is go is the gradient with respect to
remind is the neuron importance weight,
remind is the value of the coordinates (i, j) of the visualization image
A device for diagnosing the condition of skeletal muscle.
상기 전처리부가 상기 골격근에 대한 초음파 영상을 상기 소정의 형식과 다른 형식으로 상기 초음파 영상을 가공하면,
입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하는 분류망을 통해
상기 초음파 영상에 대해 픽셀 단위 연산을 수행하여 상기 골격근이 근육에코발생도(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.
상기 예측망의 상기 출력층이 소프트맥스 연산을 수행하는 출력노드로 이루어진 상기 예측망의 초기 모델을 마련하고, 상기 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행하며,
상기 초기 학습이 완료되면, 상기 예측망의 마지막 완전연결층을 초기화된 완전연결층으로 교체하고, 상기 출력층의 상기 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 상기 예측망의 정규 모델을 마련하고, 상기 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행하는 학습부;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치. 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.
상기 학습부는
상기 예측망에 대해 상기 정규 학습을 수행할 때,
손실함수
를 이용하여
상기 을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하고,
상기 을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행하고,
상기 E(y, f)는 손실함수이고,
상기 yi는 출력층의 출력값이고,
상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고,
상기 gij는 은닉층의 마지막층의 출력값이고,
상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고,
상기 i는 출력층의 노드에 대응하는 인덱스이고,
상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고,
상기 는 하이퍼파라미터인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치. 7. The method of claim 6,
the learning unit
When performing the regular learning on the prediction network,
loss function
using
remind 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 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 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.
상기 예측망은
입력층, 교번으로 반복되는 적어도 한 쌍의 컨벌루션층 및 풀링층, 적어도 하나의 완전연결층 및 출력층을 포함하며,
상기 예측부가 근육감소증을 예측하면,
상기 예측망의 상기 컨벌루션층의 출력값 및 경사도(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.
상기 시각화부는
상기 예측망의 연산의 결과를 수집한 후,
수학식
에 따라 뉴런 중요도 가중치를 산출하고,
수학식
에 따라 상기 시각화 영상을 생성하며,
상기 는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고,
상기 는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고,
상기 는 가 에 대해 가지는 경사도(gradient)이고,
상기 는 상기 뉴런 중요도 가중치이고,
상기 는 상기 시각화 영상의 좌표(i, j)의 값인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치. 10. The method of claim 9,
The visualization unit
After collecting the results of the calculation of the prediction network,
formula
Calculate the neuron importance weight according to
formula
to generate the visualization image according to
remind 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 is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer,
remind Is go is the gradient with respect to
remind is the neuron importance weight,
remind is the value of the coordinates (i, j) of the visualization image
A device for diagnosing the condition of skeletal muscle.
상기 예측망의 출력층의 출력노드가 소프트맥스 연산을 수행하는 상기 예측망의 초기 모델을 마련하고, 상기 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행하며,
상기 초기 학습이 완료되면, 상기 분류망 및 상기 예측망의 은닉층의 마지막층을 초기화된 완전연결층으로 교체하고,
상기 출력층의 상기 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 상기 예측망의 정규 모델을 마련하고, 상기 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행하는 학습부;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치. 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.
상기 학습부는
상기 예측망에 대해 상기 정규 학습을 수행할 때,
손실함수
를 이용하여
상기 을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하고,
상기 을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행하고,
상기 E(y, f)는 손실함수이고,
상기 yi는 출력층의 출력값이고,
상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고,
상기 gij는 은닉층의 마지막층의 출력값이고,
상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고,
상기 i는 출력층의 노드에 대응하는 인덱스이고,
상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고,
상기 는 하이퍼파라미터인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 장치. 12. The method of claim 11,
the learning unit
When performing the regular learning on the prediction network,
loss function
using
remind 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 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 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.
상기 예측부가 근육감소증을 예측하면, 시각화부가 상기 예측망의 상기 컨벌루션층의 출력값 및 경사도(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.
상기 시각화 영상을 생성하는 단계는
상기 시각화부가 상기 예측망의 연산의 결과를 수집하는 단계;
상기 시각화부가
수학식
에 따라 뉴런 중요도 가중치를 산출하는 단계;
상기 시각화부가
수학식
에 따라 상기 시각화 영상을 생성하는 단계;
를 포함하며,
상기 는 컨벌루션층에서 컨벌루션 연산을 수행하고 활성화함수에 의한 연산을 수행하기 전, 출력의 c번째 값이고,
상기 는 컨벌루션층의 k번째 특징지도의 좌표(i, j)의 값이고,
상기 는 상기 가 상기 에 대해 가지는 경사도(gradient)이고,
상기 는 상기 뉴런 중요도 가중치이고,
상기 는 상기 시각화 영상의 좌표(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
calculating a neuron importance weight according to
the visualization unit
formula
generating the visualization image according to
includes,
remind 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 is the value of the coordinates (i, j) of the k-th feature map of the convolutional layer,
remind is said is reminded is the gradient with respect to
remind is the neuron importance weight,
remind is the value of the coordinates (i, j) of the visualization image
A method for diagnosing a condition of skeletal muscle.
식별부가 상기 골격근에 대한 초음파 영상으로부터 추출된 복수의 방사성 특징이 입력되면, 입력층, 은닉층 및 출력층을 포함하는 식별망을 통해 상기 복수의 방사성 특징에 대한 연산을 수행하여 상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하는 단계; 및
상기 식별부가 산출된 확률에 따라 근육감소증 진행 여부를 식별하는 단계;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 방법. 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.
상기 골격근에 근육감소증이 진행되고 있는지 여부를 확률로 산출하는 단계 전,
학습부가 상기 예측망의 상기 출력층이 소프트맥스 연산을 수행하는 출력노드로 이루어진 상기 예측망의 초기 모델을 마련하는 단계; 및
상기 학습부가 상기 초기 모델에 대해 공개 데이터베이스에 저장된 영상을 학습 데이터로 이용하여 초기 학습을 수행하는 단계;
상기 학습부가 상기 예측망의 마지막 완전연결층을 초기화된 완전연결층으로 교체하고, 상기 출력층의 상기 소프트맥스 연산을 수행하는 출력노드를 시그모이드 연산을 수행하는 출력노드로 교체하여 상기 예측망의 정규 모델을 마련하는 단계; 및
상기 학습부가 상기 정규 모델에 대해 학습 데이터로 초음파 영상을 이용하여 정규 학습을 수행하는 단계;
를 더 포함하는 것을 특징으로 하는
골격근의 상태를 진단하기 위한 방법. 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.
상기 정규 학습을 수행하는 단계는
상기 학습부가
손실함수
를 이용하여
상기 을 0으로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실이 최소가 되도록 예측망의 파라미터를 수정하는 분류 손실 최적화를 수행하는 단계;
상기 학습부가 상기 을 0.5로 설정한 후, 상기 예측망의 출력값과 분류 레이블의 차이인 분류 손실과 상기 예측망의 은닉층의 마지막층의 출력값과 은닉 레이블과의 차이인 은닉 손실을 포함하는 복합 손실이 최소가 되도록 예측망의 파라미터를 수정하는 복합 손실 최적화를 수행하는 단계;
를 포함하며,
상기 E(y, f)는 손실함수이고,
상기 yi는 출력층의 출력값이고,
상기 ei는 출력층의 출력값에 대응하는 분류 레이블이고,
상기 gij는 은닉층의 마지막층의 출력값이고,
상기 cij는 은닉층의 마지막층의 출력값에 대응하는 은닉 레이블이고,
상기 i는 출력층의 노드에 대응하는 인덱스이고,
상기 j는 은닉층의 마지막층의 노드에 대응하는 인덱스이고,
상기 는 하이퍼파라미터인 것을 특징으로 하는
골격근의 상태를 진단하기 위한 방법. 19. The method of claim 18,
The step of performing the regular learning is
the learning department
loss function
using
remind 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 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 is a hyperparameter
A method for diagnosing a condition of skeletal muscle.
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)
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)
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)
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)
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 |
-
2020
- 2020-06-29 KR KR1020200079181A patent/KR102436035B1/en active IP Right Grant
Patent Citations (4)
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)
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)
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 |