KR20200032477A - Apparatus and method for estimating bone structure using ultrasonic nonlinearity - Google Patents
Apparatus and method for estimating bone structure using ultrasonic nonlinearity Download PDFInfo
- Publication number
- KR20200032477A KR20200032477A KR1020180111553A KR20180111553A KR20200032477A KR 20200032477 A KR20200032477 A KR 20200032477A KR 1020180111553 A KR1020180111553 A KR 1020180111553A KR 20180111553 A KR20180111553 A KR 20180111553A KR 20200032477 A KR20200032477 A KR 20200032477A
- Authority
- KR
- South Korea
- Prior art keywords
- bone
- frequency
- bone structure
- power spectrum
- ultrasound
- 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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0875—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of bone
-
- 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/54—Control of the diagnostic device
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Rheumatology (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
본 발명은 초음파의 비선형 특성을 이용한 골구조 예측 장치 및 방법에 관한 것이다.The present invention relates to an apparatus and method for predicting bone structure using nonlinear characteristics of ultrasound.
일반적으로, 골다공증은 골량의 감소와 골구조의 파괴로 인해 작은 충격에도 골절이 쉽게 발생하는 전신적 골질환으로 정의되며, 고령 인구가 증가하고 있는 현대사회에서 골다공증은 당뇨병 및 심혈관 질환과 함께 가장 심각한 노인성 질환으로 인식되고 있다.In general, osteoporosis is defined as a systemic bone disease in which fractures easily occur even with a small impact due to a decrease in bone mass and destruction of the bone structure. Osteoporosis is the most serious senile age in the modern society, where the elderly population is increasing, along with diabetes and cardiovascular disease. It is recognized as a disease.
골다공증은 발병하고 나면 정상 상태로 되돌릴 수 없는 비가역성 질환이라는 점에서 조기 진단 및 진단에 의한 예방이 중요시되는 질병으로, 종래에 골다공증 진단을 위한 여러 가지 진단방법 중 대표적으로 이중 에너지 X-선 흡수계측법(dual energy X-ray absorptiometry; DEXA)을 이용한 골밀도 측정이 있다.Osteoporosis is an irreversible disease that is irreversible and cannot be returned to its normal state after onset. As a result, early diagnosis and prevention by diagnosis are important. Among the various diagnostic methods for diagnosing osteoporosis, dual energy X-ray absorption measurement method ( Bone density measurement using dual energy X-ray absorptiometry (DEXA).
그러나, 골다공증 진단을 위한 골밀도 측정은 골다공증의 진단 뿐 아니라 치료반응의 평가를 위하여 반복적인 측정이 이루어져야 하기 때문에 방사선을 이용하는 이중 에너지 X-선 흡수계측법은 방사선에 반복적으로 노출되어야 하는 위험성이 존재하였으며, 특히 임산부의 골밀도 측정에는 어려움이 있었다.However, since the measurement of bone density for the diagnosis of osteoporosis needs to be repeated for the diagnosis of osteoporosis as well as for the evaluation of treatment response, the dual energy X-ray absorption measurement method using radiation has a risk of being repeatedly exposed to radiation. In particular, it was difficult to measure bone density in pregnant women.
그에 따라, 인체에 무해한 골다공증 진단 방법에 대한 연구가 다각도로 이루어지고 있으며, 이러한 연구의 일환으로 초음파를 이용하여 골다공증을 진단하는 정량적 초음파(quantitative ultrasound; QUS) 기술이 제시되었다. 정량적 초음파 기술은 인체의 발 뒤꿈치 뼈 즉, 종골(calcaneus)에 초음파를 투과시키고, 종골의 음속(speed of sound; SOS) 및 광대역 초음파 감쇠(normalized broadband ultrasound attenuation; nBUA)와 같은 초음파 특성을 측정하여 전신의 골밀도(bone mineral density; BMD)를 간접적으로 예측하여 골다공증을 진단하는 방법으로써, 방사선을 이용한 골다공증 진단 방법에 비해 반복적인 검사에도 방사선 노출의 위험성이 전혀 없고, 진단기기의 가격 및 검진비용이 상대적으로 저렴하다.Accordingly, research on a method of diagnosing osteoporosis harmless to the human body has been conducted at various angles, and as part of such research, a quantitative ultrasound (QUS) technique for diagnosing osteoporosis using ultrasound has been proposed. Quantitative ultrasound technology transmits ultrasound to the heel bone of the human body, that is, the calcaneus, and measures ultrasound characteristics such as speed of sound (SOS) and normalized broadband ultrasound attenuation (nBUA). As a method of diagnosing osteoporosis by indirectly predicting the bone mineral density (BMD) of the whole body, there is no risk of radiation exposure even after repeated examinations compared to the method of diagnosing osteoporosis using radiation. It is relatively inexpensive.
그런데, 동일한 골밀도를 갖더라도 연령이 높은 경우나 골절상의 병력 및 스테로이드류의 약물 치료를 받은 경우는 골절률이 더 높으며, 골다공증 약제의 골절률 감소 효과는 골밀도의 변화만으로는 설명하기가 어렵다. 따라서, 골밀도뿐만 아니라 골구조의 변화도 함께 고려하여 종합적으로 골강도를 예측해야만 골다공증의 원인을 규명하고 보다 정확하게 골다공증을 진단할 수 있다.However, even if they have the same bone density, they have a higher fracture rate when they are older or have a history of fractures and are treated with steroids, and the effect of reducing the fracture rate of osteoporotic drugs is difficult to explain by only changing the bone density. Therefore, it is necessary to comprehensively predict bone strength by considering not only bone density but also changes in bone structure, so that the cause of osteoporosis can be identified and the osteoporosis diagnosed more accurately.
최근에는 초음파 측정 방식을 이용하여 골구조를 예측하는 기술이 연구되고 있으며, 대한민국 공개특허공보 제10-2018-0037351호(이하, '종래기술'이라 함)에는 초음파 비선형 변수를 이용한 골구조 예측 장치 및 방법에 관한 기술이 제시된 바 있다. 종래기술은 특정 수식을 이용하여 해면질골의 비선형 변수를 구한 후, 이를 이용하여 해면질골의 골구조를 예측하는 방식이었으나, 비선형 변수를 도출하기 위한 특정 수식 내에는 많은 변수들이 존재하고 수식이 복잡함에 따라 수식을 계산하는 과정이 복잡해지고, 측정 오차가 발생할 우려가 있으며, 골구조의 예측시 정밀도가 저하될 수 있으므로 종래기술보다 계산을 간소화하면서도 보다 정밀하게 골구조를 예측하는 기술의 개발이 시급한 실정이다. Recently, a technique for predicting a bone structure using an ultrasonic measurement method has been researched, and in Korean Patent Publication No. 10-2018-0037351 (hereinafter referred to as 'prior art'), a bone structure prediction apparatus using ultrasonic nonlinear variables And methods have been presented. The prior art was a method of obtaining a nonlinear variable of the spongy bone using a specific formula and then predicting the bone structure of the spongy bone using this, but there are many variables in the specific formula for deriving the nonlinear variable and the formula is complicated. Accordingly, the process of calculating the formula is complicated, there is a risk of measurement error, and the precision may be lowered when predicting the bone structure. Therefore, it is urgent to develop a technology that more accurately predicts the bone structure while simplifying the calculation than the prior art. to be.
본 발명은 상술한 문제점을 해결하기 위한 것으로, 골구조 예측시 사용되는 변수의 개수를 줄이고, 수식을 간소화하여 간단한 계산으로도 손쉽게 골구조를 예측할 수 있고, 측정 오차를 줄여 골구조 예측 정밀도를 높이는 기술을 제공하는데 목적이 있다.The present invention is to solve the above-mentioned problems, by reducing the number of variables used in predicting the bone structure, by simplifying the formula can easily predict the bone structure with a simple calculation, reducing the measurement error to increase the accuracy of bone structure prediction The purpose is to provide technology.
본 발명이 해결하려는 과제는 전술한 과제로 제한되지 아니하며, 언급되지 아니한 또 다른 기술적 과제들은 후술할 내용으로부터 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 명확하게 이해될 수 있을 것이다.The problems to be solved by the present invention are not limited to the above-described problems, and other technical problems that are not mentioned will be clearly understood by a person having ordinary knowledge in the technical field to which the present invention pertains from contents to be described later.
이러한 목적을 달성하기 위하여 본 발명의 일 실시형태로서, 초음파의 비선형 특성을 이용한 골구조 예측 장치는 인가된 전기적 신호를 초음파로 변환하여 해면질골에 입사시키는 초음파 송신부; 상기 해면질골을 통과한 초음파를 수신하고, 이를 전기적 신호로 변환하는 초음파 수신부; 상기 해면질골에 초음파가 입사되도록 제1주파수를 갖는 송신 초음파에 대한 전기적 신호를 생성하여 상기 초음파 송신부로 전송하는 파형 발생부; 상기 초음파 수신부에서 변환된 전기적 신호를 검출하는 신호처리부; 상기 신호처리부에서 검출된 전기적 신호를 변환함으로써, 상기 제1주파수 및 상기 제1주파수의 2배인 주파수인 제2주파수에 대한 파워 스펙트럼 레벨을 도출하는 변환부; 상기 제1주파수의 파워 스펙트럼 레벨과 상기 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 산출하는 연산부; 및 상기 연산부에서 산출된 값과 골구조 간의 상관관계를 이용하여 상기 해면질골의 골구조를 예측하는 골구조 예측부;를 포함할 수 있다.In order to achieve this object, as an embodiment of the present invention, an apparatus for predicting bone structure using nonlinear characteristics of ultrasonic waves includes: an ultrasonic transmitter that converts an applied electrical signal into ultrasonic waves and enters the spongy bone; An ultrasonic wave receiving unit that receives ultrasonic waves that have passed through the spongy bone and converts them into electrical signals; A waveform generator that generates an electrical signal for transmitting ultrasonic waves having a first frequency so that ultrasonic waves are incident on the spongy bone and transmits them to the ultrasonic transmitter; A signal processor that detects the electrical signal converted by the ultrasonic receiver; A conversion unit that derives a power spectrum level for the first frequency and a second frequency that is twice the frequency of the first frequency by converting the electrical signal detected by the signal processing unit; A calculating unit calculating a difference value between the power spectrum level of the first frequency and the power spectrum level of the second frequency; And a bone structure prediction unit for predicting the bone structure of the spongy bone by using a correlation between the value calculated by the operation unit and the bone structure.
그리고, 해면질골의 골구조는 골소주 선밀도를 포함할 수 있다.In addition, the bone structure of the spongy bone may include bone shochu linear density.
또한, 해면질골의 골구조는 골 용적비율을 포함할 수 있다.In addition, the bone structure of the spongy bone may include a bone volume ratio.
아울러, 골구조 예측부는 상기 연산부에서 산출된 값을 1차 선형 방정식에 대입함으로써, 상기 골소주 선밀도를 예측하되, 상기 1차 선형 방정식은, 해면질골 샘플에 초음파를 투과시켜 도출한 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 독립변수로 설정하고, 상기 해면질골 샘플의 골소주 선밀도를 종속변수로 설정하며, 상기 독립변수와 종속변수를 선형회귀분석함으로써 도출될 수 있다.In addition, the bone structure prediction unit predicts the linear density of the bone shochu by substituting the value calculated by the calculation unit into a linear linear equation, wherein the primary linear equation is the power of the first frequency derived by transmitting ultrasound to the spongy bone sample. It can be derived by setting the difference between the spectral level and the power spectral level of the second frequency as an independent variable, setting the bone mineral density of the spongy bone sample as a dependent variable, and linear regression analysis of the independent variable and the dependent variable.
그리고, 골구조 예측부는 상기 연산부에서 산출된 값을 1차 선형 방정식에 대입함으로써, 상기 골 용적비율을 예측하되, 상기 1차 선형 방정식은, 해면질골 샘플에 초음파를 투과시켜 도출한 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 독립변수로 설정하고, 상기 해면질골 샘플의 골 용적비율을 종속변수로 설정하며, 상기 독립변수와 종속변수를 선형회귀분석함으로써 도출될 수 있다.In addition, the bone structure prediction unit predicts the bone volume ratio by substituting a value calculated by the calculation unit into a first-order linear equation, wherein the first linear equation has a first frequency derived by transmitting ultrasound to a spongy bone sample. It can be derived by setting the difference between the power spectrum level and the power spectrum level of the second frequency as an independent variable, setting the bone volume ratio of the spongy bone sample as a dependent variable, and linear regression analysis of the independent variable and the dependent variable. have.
또한, 해면질골은 상기 초음파 송신부 및 초음파 수신부의 사이에 위치하되, 전달 매질의 특성에 의해 발생하는 신호 왜곡을 최소화하도록 상기 해면질골과 초음파 수신부 간의 직선거리는 상기 해면질골과 초음파 송신부 간의 직선거리에 비해 상대적으로 짧게 형성될 수 있다.In addition, the spongy bone is located between the ultrasonic transmitter and the ultrasonic receiver, the linear distance between the spongy bone and the ultrasonic receiver to minimize the signal distortion caused by the characteristics of the transmission medium compared to the linear distance between the spongy bone and the ultrasonic transmitter It can be formed relatively short.
이러한 목적을 달성하기 위하여 본 발명의 다른 실시형태로서, 초음파의 비선형 특성을 이용한 골구조 예측 방법은 파형 발생부가 제1주파수를 갖는 송신 초음파에 대한 전기적 신호를 생성하여 초음파 송신부로 전송하면, 상기 초음파 송신부가 해면질골에 초음파를 입사시키는 초음파 입사 단계; 초음파 수신부가 상기 해면질골을 통과한 초음파를 수신하고, 이를 전기적 신호로 변환하는 신호 수신 단계; 신호처리부가 상기 초음파 수신부에서 변환된 전기적 신호를 검출하는 신호 검출 단계; 변환부가 상기 신호처리부에서 검출된 전기적 신호를 변환함으로써, 상기 제1주파수 및 상기 제1주파수의 2배인 주파수인 제2주파수에 대한 파워 스펙트럼 레벨을 도출하는 변환 단계; 연산부가 상기 제1주파수의 파워 스펙트럼 레벨과 상기 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 산출하는 연산 단계; 및 골구조 예측부가 상기 연산부에서 산출된 값과 골구조 간의 상관관계를 이용하여 상기 해면질골의 골구조를 예측하는 골구조 예측 단계;를 포함할 수 있다.In order to achieve this object, as another embodiment of the present invention, a method for predicting bone structure using a nonlinear characteristic of ultrasonic waves is generated when the waveform generator generates an electrical signal for transmitting ultrasonic waves having a first frequency and transmits them to the ultrasonic transmitter. An ultrasonic incident step in which the transmitter injects ultrasound into the spongy bone; A signal receiving step of receiving an ultrasonic wave through the spongy bone and converting it into an electrical signal; A signal detection step in which a signal processing unit detects an electrical signal converted by the ultrasonic receiving unit; A conversion step of converting an electrical signal detected by the signal processing unit to derive a power spectrum level for the first frequency and a second frequency that is twice the frequency of the first frequency; A calculating step of calculating a difference between the power spectrum level of the first frequency and the power spectrum level of the second frequency; And a bone structure prediction step in which the bone structure prediction unit predicts the bone structure of the spongy bone by using a correlation between the value calculated by the operation unit and the bone structure.
그리고, 골구조 예측 단계에서 상기 해면질골의 골구조는 골소주 선밀도를 포함하고, 상기 골구조 예측부는 상기 연산부에서 산출된 값을 1차 선형 방정식에 대입함으로써, 상기 골소주 선밀도를 예측하되, 상기 1차 선형 방정식은, 해면질골 샘플에 초음파를 투과시켜 도출한 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 독립변수로 설정하고, 상기 해면질골 샘플의 골소주 선밀도를 종속변수로 설정하며, 상기 독립변수와 종속변수를 선형회귀분석함으로써 도출될 수 있다.And, in the bone structure prediction step, the bone structure of the spongy bone includes bone suso linear density, and the bone structure predictor predicts the bone suso linear density by substituting a value calculated by the calculation unit into a linear linear equation. In the linear equation, the difference value between the power spectrum level of the first frequency and the power spectrum level of the second frequency derived by transmitting ultrasound to the spongy bone sample is set as an independent variable, and the linear density of the bone shochu of the spongy bone sample is a dependent variable. It can be derived by linear regression analysis of the independent and dependent variables.
또한, 골구조 예측 단계에서 상기 해면질골의 골구조는 골 용적비율을 포함하고, 상기 골구조 예측부는 상기 연산부에서 산출된 값을 1차 선형 방정식에 대입함으로써, 상기 골 용적비율을 예측하되, 상기 1차 선형 방정식은, 해면질골 샘플에 초음파를 투과시켜 도출한 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 독립변수로 설정하고, 상기 해면질골 샘플의 골 용적비율을 종속변수로 설정하며, 상기 독립변수와 종속변수를 선형회귀분석함으로써 도출될 수 있다.In addition, in the bone structure prediction step, the bone structure of the spongy bone includes a bone volume ratio, and the bone structure prediction unit predicts the bone volume ratio by substituting a value calculated by the calculation unit into a linear equation. The first linear equation sets the difference value between the power spectrum level of the first frequency and the power spectrum level of the second frequency derived by transmitting ultrasound to the spongy bone sample as an independent variable, and sets the bone volume ratio of the spongy bone sample. It is set as a dependent variable, and can be derived by linear regression analysis of the independent variable and the dependent variable.
그리고, 초음파 입사 단계 및 신호 수신 단계에서 상기 해면질골은 상기 초음파 송신부 및 초음파 수신부의 사이에 위치하되, 전달 매질의 특성에 의해 발생하는 신호 왜곡을 최소화하도록 상기 해면질골과 초음파 수신부 간의 직선거리는 상기 해면질골과 초음파 송신부 간의 직선거리에 비해 상대적으로 짧게 형성될 수 있다.Further, in the ultrasonic incidence step and the signal receiving step, the spongy bone is located between the ultrasound transmitter and the ultrasound receiver, but the linear distance between the spongy bone and the ultrasound receiver is minimized to minimize signal distortion caused by characteristics of the transmission medium. It may be formed relatively short compared to the straight line distance between the bone and the ultrasound transmitter.
상술한 과제의 해결 수단은 단지 예시적인 것으로서, 본 발명을 제한하려는 의도로 해석되지 않아야 한다. 상술한 예시적인 실시예 외에도, 도면 및 발명의 상세한 설명에 기재된 추가적인 실시예가 존재할 수 있다.The means for solving the above-described problems are merely exemplary and should not be construed as limiting the present invention. In addition to the exemplary embodiments described above, there may be additional embodiments described in the drawings and detailed description of the invention.
이상에서 설명한 바와 같이, 본 발명의 일 실시예에 의하면 골구조의 예측시 사용되는 변수의 개수를 최소로 줄이고, 수식을 간소화시킴으로써 간단한 계산으로도 골구조를 예측할 수 있으므로 많은 변수와 복잡한 수식을 이용하여 골구조를 예측하는 기술에 비해 편리하고 손쉽게 골구조를 예측할 수 있다.As described above, according to an embodiment of the present invention, the number of variables used in the prediction of the bone structure is reduced to a minimum, and by simplifying the equation, the bone structure can be predicted by simple calculation, so many variables and complex equations are used. Therefore, it is more convenient and easier to predict the bone structure than the technology for predicting the bone structure.
또한, 일 실시예에서는 골구조의 예측시 사용되는 변수로서, 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨만을 이용하며, 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 통해 간단하게 골구조를 예측할 수 있으므로 골구조의 예측시 사용되는 변수의 개수가 기존에 비해 대폭 감소되고, 복잡한 수식 또한 간소화되어 기존에 비해 측정 오차를 줄일 수 있고, 결과적으로 골구조 예측 정밀도가 향상되는 효과가 있다.In addition, in one embodiment, as a variable used in the prediction of the bone structure, only the power spectrum level of the first frequency and the power spectrum level of the second frequency are used, and the power spectrum level of the first frequency and the power spectrum level of the second frequency are used. Since the bone structure can be easily predicted through the difference value between, the number of variables used in the prediction of the bone structure is greatly reduced compared to the previous one, and the complex formula is also simplified, reducing the measurement error compared to the previous one, and consequently the bone structure. There is an effect that the prediction precision is improved.
본 발명의 다양한 실시예에 따른 효과들은 이상에서 언급한 효과들로 제한되지 않으며, 언급되지 않은 또 다른 효과들은 청구범위의 기재로부터 통상의 기술자에게 명확하게 이해될 수 있을 것이다.Effects according to various embodiments of the present invention are not limited to the above-mentioned effects, and other effects not mentioned will be clearly understood by those skilled in the art from the claims.
도1은 본 발명의 일 실시예에 따른 초음파의 비선형 특성을 이용한 골구조 예측 장치를 개략적으로 도시한 개념도이다.
도2는 해면질골 샘플들로부터 획득한 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값 및 골 용적비율 간의 상관관계를 도시한 그래프이다.
도3은 해면질골 샘플들로부터 획득한 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값 및 골소주 선밀도 간의 상관관계를 도시한 그래프이다.
도4는 본 발명의 일 실시예에 따른 초음파의 비선형 특성을 이용한 골구조 예측 방법을 도시한 흐름도이다.1 is a conceptual diagram schematically illustrating an apparatus for predicting bone structure using nonlinear characteristics of ultrasound according to an embodiment of the present invention.
2 is a graph showing the correlation between the difference between the first power spectrum level and the second power spectrum level obtained from the spongy bone samples and the bone volume ratio.
FIG. 3 is a graph showing the correlation between the difference between the first power spectrum level and the second power spectrum level obtained from the spongy bone samples and the line density of the bone shochu.
4 is a flowchart illustrating a method for predicting bone structure using nonlinear characteristics of ultrasound according to an embodiment of the present invention.
본 발명의 바람직한 실시예에 대하여 첨부된 도면을 참조하여 더 구체적으로 설명하되, 이미 주지되어진 기술적 부분에 대해서는 설명의 간결함을 위해 생략하거나 압축하기로 한다.With reference to the accompanying drawings, a preferred embodiment of the present invention will be described in more detail, but for the brevity of description, the technical parts that are already well-known will be omitted or compressed.
본 명세서에서 본 발명의 "일" 또는 "하나의" 실시예에 대한 언급들은 반드시 동일한 실시예에 대한 것은 아니며, 이들은 적어도 하나를 의미한다는 것에 유의해야 한다.It should be noted that references herein to “one” or “one” embodiments of the present invention are not necessarily referring to the same embodiment, and they mean at least one.
이하의 실시예에서, 제1, 제2 등의 용어는 한정적인 의미가 아니라 하나의 구성 요소를 다른 구성 요소와 구별하는 목적으로 사용되었다.In the following examples, terms such as first and second are not used in a limiting sense, but for the purpose of distinguishing one component from other components.
이하의 실시예에서, 단수의 표현은 문맥상 명백하게 다른 의미를 뜻하지 않는 한, 복수의 표현을 포함한다.In the following embodiments, the singular expressions include the plural expressions unless the context clearly indicates different meanings.
이하의 실시예에서, 포함하다 또는 가지다 등의 용어는 명세서상에 기재된 특징 또는 구성요소가 존재함을 의미하는 것이고, 하나 이상의 다른 특징들 또는 구성요소가 부가될 가능성을 미리 배제하는 것은 아니다.In the following embodiments, terms such as include or have mean that the features or components described in the specification exist, and do not preclude the possibility that one or more other features or components may be added.
도면에서는 설명의 편의를 위하여 구성 요소들이 그 크기가 과장 또는 축소될 수 있다. 예컨대, 도면에서 나타난 각 구성의 크기 및 두께는 설명의 편의를 위해 임의로 나타내었으므로, 본 발명이 반드시 도시된 바에 한정되지 않는다.In the drawings, the size of components may be exaggerated or reduced for convenience of description. For example, since the size and thickness of each component shown in the drawings are arbitrarily shown for convenience of description, the present invention is not necessarily limited to what is shown.
<초음파의 비선형 특성을 이용한 골구조 예측 장치에 대한 설명><Description of a device for predicting bone structure using nonlinear characteristics of ultrasonic waves>
도1은 본 발명의 일 실시예에 따른 초음파의 비선형 특성을 이용한 골구조 예측 장치를 개략적으로 도시한 개념도이다. 도1을 참조하면, 초음파의 비선형 특성을 이용한 골구조 예측 장치(10)는 하우징(100), 초음파 송신부(200), 초음파 수신부(300), 파형 발생부(400), 신호처리부(500), 변환부(600), 연산부(700), 골구조 예측부(800) 및 출력부(900)를 포함할 수 있다.1 is a conceptual diagram schematically illustrating an apparatus for predicting bone structure using nonlinear characteristics of ultrasound according to an embodiment of the present invention. Referring to Figure 1, the bone
하우징(100)은 매질(110)을 수용할 수 있는 내부 공간을 가지며, 상부가 개방된 형태로 적용될 수 있다. 일 실시예에서 매질(110)은 초음파를 전파하는 역할을 하며, 인체의 연부조직과 유사한 음향특성을 갖는 물(일 예로, 증류수)이 초음파 전파 매질로서 이용될 수 있다. The
초음파 송신부(200)는 후술될 파형 발생부(400)로부터 전기적 신호를 수신하면, 인가된 전기적 신호를 초음파로 변환하여 해면질골(20)에 입사시킬 수 있다. 그리고, 초음파 수신부(300)는 해면질골(20)을 통과한 초음파를 수신하고, 이를 전기적 신호로 변환할 수 있다. 이때, 초음파 송신부(200), 초음파 수신부(300) 및 해면질골(20)은 매질(110)로 채워진 상태의 하우징(100)의 내부 공간에 위치할 수 있다.When the
그리고, 일 실시예에서 해면질골(20)은 초음파 송신부(200) 및 초음파 수신부(300)의 사이에 위치하되, 해면질골(20) 일면과 초음파 수신부(300) 간의 직선거리(D2)는 해면질골(20) 타면과 초음파 송신부(200) 간의 직선거리(D1)에 비해 상대적으로 짧게 형성됨으로써, 전달 매질의 특성에 의해 발생하는 신호 왜곡을 최소화시킬 수 있다. 즉, 초음파 송신부(200)로부터 조사된 초음파는 해면질골(20)을 투과한 후에 초음파 수신부(300)에 도달하게 되는데, 해면질골(20)을 통과한 초음파가 초음파 수신부(300)로 이동하는 과정에서 전달 매질의 비선형 특성에 의해 초음파 신호의 진폭이나 파형이 미묘하게 달라질 수 있기 때문에 해면질골(20) 일면과 초음파 수신부(300) 간의 직선거리(D2)는 해면질골(20) 타면과 초음파 송신부(200) 간의 직선거리(D1)에 비해 상대적으로 짧게 형성시켜 물에서의 비선형 특성을 최소화하여 초음파 측정의 정확도를 높일 수 있다.In addition, in one embodiment, the
한편, 실시하기에 따라, 하우징(100) 내에는 초음파 송신부(200), 초음파 수신부(300) 및 해면질골(20)의 위치를 고정할 수 있는 별도의 고정장치(미도시)가 더 포함될 수도 있다.On the other hand, according to the implementation, the
파형 발생부(400)는 해면질골(20)에 초음파가 입사되도록 제1주파수를 갖는 송신 초음파에 대한 전기적 신호를 생성하여 초음파 송신부(200)로 전송할 수 있다. 여기서, 제1주파수는 송신 초음파가 갖는 기본 주파수(일 예로, 0.5MHz)를 의미한다.The
신호처리부(500)는 초음파 수신부(300)에서 변환된 전기적 신호를 검출할 수 있다. 초음파 송신부(200)에서 기본 주파수를 갖는 송신 초음파를 조사 대상물체인 해면질골(20)에 입사시키면, 입사된 초음파는 초음파 수신부(300)에 수신된다. 그리고, 수신된 초음파는 초음파 수신부(300)에서 전기적 신호로 변환된다. 이때, 변환된 전기적 신호로부터 기본 주파수 및 기본 주파수의 2배, 3배 등의 주파수 성분이 함께 측정되는 현상을 초음파의 비선형 특성이라고 한다.The
변환부(600)는 신호처리부(500)에서 검출된 전기적 신호를 변환함으로써, 주파수 함수에 대한 파워 스펙트럼 레벨을 구할 수 있다. 즉, 변환부(600)는 제1주파수에 대한 파워 스펙트럼 레벨(이하, '제1파워 스펙트럼 레벨'이라 함)과 제1주파수의 2배인 주파수인 제2주파수(일 예로, 1.0MHz)에 대한 파워 스펙트럼 레벨(이하, '제2파워 스펙트럼 레벨'이라 함)을 도출할 수 있다. The
예를 들어, 변환부(600)는 신호처리부(500)가 검출한 전기적 신호를 고속 푸리에 변환함으로써, 0.5MHz인 제1주파수의 파워 스펙트럼 레벨은 12dB이고, 1.0MHz인 제2주파수의 파워 스펙트럼 레벨은 -35dB이라고 도출할 수 있다.For example, the
연산부(700)는 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값을 산출할 수 있다. 즉, 연산부(700)는 제1파워 스펙트럼 레벨에서 제2파워 스펙트럼 레벨을 감한 값을 산출할 수 있다. 일 예로, 제1파워 스펙트럼 레벨이 12dB이고, 제2파워 스펙트럼 레벨이 -35dB일 때, 연산부(700)는 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값이 47dB라고 산출하게 된다.The
골구조 예측부(800)는 연산부(700)에서 산출된 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값 및 골구조 간의 상관관계를 이용하여 해면질골(20)의 골구조를 예측할 수 있다. 여기서, 해면질골(20)의 골구조는 골소주 선밀도에 대한 데이터 및 골 용적비율에 대한 데이터 중 적어도 어느 하나를 포함할 수 있다.The bone structure prediction unit 800 may predict the bone structure of the
이때, 골소주 선밀도는 단위부피당 골소주의 양을 의미하며, 골 용적비율 값을 골소주의 두께 값으로 나눈 것이다. 그리고, 골 용적비율은 기둥 형상의 골소주(bone trabecula)를 포함하는 조직 체적에서 골소주가 차지하는 골체적의 비율을 의미한다.At this time, the line density of bone shochu means the amount of bone shochu per unit volume, and is the bone volume ratio value divided by the thickness value of bone shochu. In addition, the bone volume ratio refers to the ratio of the bone volume occupied by the bone shochu in the tissue volume including the pillar-shaped bone trabecula.
본 발명자는 골 용적비율이 서로 다른 해면질골 샘플에 각각 초음파를 조사했을 때, 해면질골 샘플마다 제1파워 스펙트럼 레벨값과 제2파워 스펙트럼 레벨값이 상이하다는 것을 확인하였고, 해면질골 샘플에 초음파를 조사하여 획득한 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값 및 해면질골 샘플의 골구조(예를 들어, 골소주 선밀도, 골 용적비율) 간에는 긴밀한 상관관계가 있다는 것을 밝혀냄으로써, 이를 이용해 실제 환자의 해면질골의 골구조를 예측할 수 있는 기술을 완성하였다. The present inventors confirmed that the first power spectrum level value and the second power spectrum level value were different for each spongy bone sample when ultrasound was applied to the spongy bone samples having different bone volume ratios, and the ultrasound was applied to the spongy bone samples. It is revealed that there is a close correlation between the difference between the first power spectral level and the second power spectral level obtained by investigation and the bone structure of the spongy bone sample (eg, bone somatic density, bone volume ratio). The technology to predict the bone structure of the spongy bone of the patient was completed.
즉, 환자의 해면질골에 초음파를 조사하여 실제 해면질골의 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값을 구하고, 구한 차이값을 선형 방정식에 대입하면 환자의 골구조(골 용량비율 또는 골소주 선밀도)를 쉽게 예측할 수 있다.That is, when the patient's spongy bone is irradiated with ultrasound to obtain a difference value between the first power spectrum level and the second power spectrum level of the actual spongy bone, and substituting the difference value into the linear equation, the patient's bone structure (bone capacity ratio or Bone shochu linear density) can be easily predicted.
이에 대해 보다 구체적으로 설명하면 다음과 같다. 본 발명자는 소의 대퇴골 근위부를 이용하여 제작된 32개의 해면질골 샘플(두께 15mm의 직육면체 형태)을 대상으로 각각 초음파를 조사하였다. 그로부터 각 샘플의 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨을 구하고, 제1파워 스펙트럼 레벨 값에서 제2파워 스펙트럼 레벨 값을 감한 차이값을 도출하였다. 그리고, 마이크로 전산화 단층 촬영(micro-computed tomography)을 이용하여 각 해면질골 샘플을 촬영하고, 촬영된 이미지를 소정의 프로그램(일 예로, skyscan CTAn 소프트웨어)으로 분석하여 각 해면질골 샘플의 골 용적비율, 골소주의 두께, 골소주 선밀도 값 등을 도출하였다.This will be described in more detail as follows. The present inventors respectively irradiated ultrasound on 32 spongy bone samples (rectangular cuboids having a thickness of 15 mm) produced using the proximal portion of the femur. From this, the first power spectrum level and the second power spectrum level of each sample were obtained, and a difference value obtained by subtracting the second power spectrum level value from the first power spectrum level value was derived. Then, using a micro-computed tomography (micro-computed tomography) to take a sample of each spongy bone, by analyzing the image taken with a predetermined program (eg, skyscan CTAn software), the bone volume ratio of each spongy bone sample, The thickness of bone shochu and the line density of bone shochu were derived.
그 후, 도출된 해면질골 샘플의 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값을 독립변수로 설정하고, 해면질골 샘플의 골 용적비율을 종속변수로 설정하며, 독립변수와 종속변수를 이용하여 선형회귀분석함으로써 도2의 1차 선형 방정식을 도출하였다.Then, the difference value between the first power spectrum level and the second power spectrum level of the derived spongy bone sample is set as an independent variable, the bone volume ratio of the spongy bone sample is set as a dependent variable, and the independent variable and the dependent variable are set. The linear linear equation of FIG. 2 was derived by using linear regression analysis.
도2는 32개의 해면질골 샘플들로부터 획득한 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값 및 골 용적비율 간의 상관관계를 도시한 그래프이다. 도2를 참조하면, 선형회귀분석을 통해 얻어진 피어슨(Pearson) 상관계수 R이 -0.85이므로 제1파워 스펙트럼 레벨과 제1파워 스펙트럼 레벨 간의 차이값은 골 용적비율과 강한 음의 상관관계를 갖는다는 것을 확인할 수 있다. 여기서, 상관관계는 독립변수와 종속변수 간의 관계의 강도를 의미한다. 또한, 도2의 그래프에 도시된 1차 선형 방정식(즉, 골 용적비율 도출용 방정식)은 하기 수학식 1로 표현될 수 있다.2 is a graph showing the correlation between the difference between the first power spectrum level and the second power spectrum level and the bone volume ratio obtained from 32 spongy bone samples. Referring to FIG. 2, since the Pearson correlation coefficient R obtained through the linear regression analysis is -0.85, the difference value between the first power spectrum level and the first power spectrum level has a strong negative correlation with the bone volume ratio. You can confirm that. Here, the correlation means the strength of the relationship between the independent variable and the dependent variable. In addition, the first linear equation (ie, the equation for deriving the bone volume ratio) illustrated in the graph of FIG. 2 may be expressed by Equation 1 below.
[수학식 1][Equation 1]
Y=A1X+B1 Y = A 1 X + B 1
(여기서, Y는 종속변수인 해면질골 샘플의 골 용적비율, X는 독립변수인 해면질골 샘플의 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값, A1은 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값이 증가함에 따른 골 용적비율의 변화량, B1은 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값이 20일 때 골 용적비율의 값)(Where Y is the bone volume ratio of the spongy bone sample as a dependent variable, X is the difference between the first power spectrum level and the second power spectrum level of the spongy bone sample as an independent variable, A 1 is the first power spectrum level and the The amount of change in the bone volume ratio as the difference between the two power spectral levels increases, B 1 is the value of the bone volume ratio when the difference between the first power spectral level and the second power spectral level is 20)
한편, 도출된 해면질골 샘플의 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값을 독립변수로 설정하고, 해면질골 샘플의 골소주 선밀도를 종속변수로 설정하며, 독립변수와 종속변수를 선형회귀분석함으로써 도3의 1차 선형 방정식을 도출하였다.On the other hand, the difference between the first power spectrum level and the second power spectrum level of the derived spongy bone sample is set as an independent variable, the bone density of the spongy bone sample is set as a dependent variable, and the independent variable and the dependent variable are linearly regressed. By analysis, the linear equation of FIG. 3 was derived.
도3은 해면질골 샘플들로부터 획득한 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값 및 골소주 선밀도(골소주 선밀도 단위는 mm-1) 간의 상관관계를 도시한 그래프이다. 도3을 참조하면, 선형회귀분석을 통해 얻어진 피어슨 상관계수 R이 -0.71이므로 제1파워 스펙트럼 레벨과 제1파워 스펙트럼 레벨 간의 차이값은 골소주 선밀도와 강한 음의 상관관계를 갖는다는 것을 확인할 수 있다. 또한, 도3의 그래프에 도시된 1차 선형 방정식(즉, 골소주 선밀도 도출용 방정식)은 하기 수학식 2로 표현될 수 있다.FIG. 3 is a graph showing the correlation between the difference between the first power spectrum level and the second power spectrum level obtained from the spongy bone samples and the bone strain linear density (the bone strain linear density unit is mm −1 ). Referring to FIG. 3, since the Pearson correlation coefficient R obtained through the linear regression analysis is -0.71, it can be confirmed that the difference between the first power spectrum level and the first power spectrum level has a strong negative correlation with bone line density. . In addition, the first-order linear equation shown in the graph of FIG. 3 (that is, the equation for deriving bone somatic linear density) may be expressed by Equation 2 below.
[수학식 2][Equation 2]
Y=A2X+B2 Y = A 2 X + B 2
(여기서, Y는 종속변수인 해면질골 샘플의 골소주 선밀도, X는 독립변수인 해면질골 샘플의 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값, A2는 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값이 증가함에 따른 골소주 선밀도의 변화량, B2는 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값이 20일 때 골소주 선밀도의 값)(Where Y is the bone density of the bone strain of the spongy bone sample, which is a dependent variable, X is the difference between the first power spectrum level and the second power spectrum level of the spongy bone sample, which is an independent variable, A 2 is the first power spectrum level and the second The amount of change in bone mineral density as the difference between power spectral levels increases, B 2 is the value of bone mineral density when the difference between the first power spectrum level and the second power spectrum level is 20)
따라서, 수학식 1과 수학식 2를 통해 이미 1차 선형 방정식이 확보된 상태이므로, 골구조 예측부(800)는 연산부(700)에서 산출된 차이값을 수학식 1이나 수학식 2의 1차 선형 방정식의 X값에 대입함으로써, 실제 환자의 해면질골의 골 용적비율 또는 골소주 선밀도를 구할 수 있다. 물론, 이때 얻어지는 골 용적비율이나 골소주 선밀도 결과값은 해면질골 샘플의 골구조 데이터에 근거하여 환자의 골구조를 예측한 결과값이다. Therefore, since the first linear equation is already secured through Equation 1 and Equation 2, the bone structure predicting unit 800 may calculate the difference value calculated by the
즉, 해면질골 샘플의 제1파워 스펙트럼 레벨과 제2파워 스펙트럼 레벨 간의 차이값은 해면질골 샘플의 실제 골 용적비율이나 골소주 선밀도와 밀접한 상관관계를 가지고 있기 때문에 이러한 상관관계를 이용하여 실제 환자의 해면질골의 골구조를 예측할 수 있는 것이다.That is, since the difference between the first power spectral level and the second power spectral level of the spongy bone sample has a close correlation with the actual bone volume ratio or bone density of the spongy bone sample, the correlation of the actual patient's spongy quality using this correlation It is possible to predict the bone structure of the bone.
출력부(900)는 디스플레이 장치나 스피커 등으로 구현될 수 있다. 출력부(900)는 골구조 예측부(800)에서 예측된 해면질골(20)의 골구조(즉, 골 용적비율과 골소주 선밀도) 결과값을 화면에 출력하거나 음성으로 안내할 수 있다.The output unit 900 may be implemented as a display device or a speaker. The output unit 900 may output a result of the bone structure (ie, bone volume ratio and bone strain linear density) of the
<초음파의 비선형 특성을 이용한 골구조 예측 방법에 대한 설명><Explanation of bone structure prediction method using nonlinear characteristics of ultrasound>
본 발명의 일 실시예에 따른 초음파의 비선형 특성을 이용한 골구조 예측 방법에 대하여 도4에 도시된 흐름도를 따라 설명하고, 도1 내지 도3에 도시된 도면을 참조하여 설명하되, 편의상 순서를 붙여 설명하기로 한다.A method for predicting bone structure using nonlinear characteristics of ultrasonic waves according to an embodiment of the present invention will be described with reference to the flowchart shown in FIG. 4, and described with reference to the drawings shown in FIGS. I will explain.
1. 초음파 입사 단계<S401>1. Ultrasonic incidence step <S401>
본 단계에서는 파형 발생부(400)가 제1주파수를 갖는 송신 초음파에 대한 전기적 신호를 생성하여 초음파 송신부(200)로 전송하면, 초음파 송신부(200)가 해면질골(20)에 초음파를 입사시키는 과정이 이루어질 수 있다. 여기서, 초음파 송신부(200)는 제1주파수와 동일한 중심 주파수(일 예로, 0.5MHz)를 갖는 비집속형 트랜스듀서로 적용될 수 있다.In this step, when the
2. 신호 수신 단계<S402>2. Signal reception step <S402>
본 단계에서는 초음파 수신부(300)가 해면질골(20)을 통과한 초음파를 수신하고, 이를 전기적 신호로 변환하는 과정이 진행될 수 있다. 여기서, 초음파 수신부(300)는 초음파 송신부(200)의 중심 주파수의 2배인 중심 주파수(일 예로, 1.0MHz)를 갖는 비집속형 트랜스듀서로 적용될 수 있다. 만일, 초음파 수신부(300)의 중심 주파수가 초음파 송신부(200)와 동일할 경우에는 초음파의 비선형 특성에 의해 발생하는 제2주파수를 민감하게 감지하기가 어려우므로 비선형 특성을 이용한 골구조 예측 정밀도가 저하될 우려가 있다. 따라서, 초음파 수신부(300)는 초음파 송신부(200)의 중심 주파수에 비해 2배 큰 중심 주파수를 갖는 트랜스듀서로 적용하는 것이 바람직하다.In this step, the
3. 신호 검출 단계<S403>3. Signal detection step <S403>
본 단계에서 신호처리부(500)는 초음파 수신부(300)로부터 변환된 전기적 신호를 검출할 수 있다.In this step, the
4. 변환 단계<S404>4. Conversion step <S404>
본 단계에서는 변환부(600)가 신호처리부(500)에서 검출된 전기적 신호를 고속 푸리에 변환함으로써, 제1주파수 및 제1주파수의 2배인 주파수인 제2주파수에 대한 파워 스펙트럼 레벨을 도출하는 과정이 진행될 수 있다.In this step, the
5. 연산 단계<S405>5. Operation step <S405>
본 단계에서는 연산부(700)가 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 산출하는 과정이 이루어질 수 있다. 즉, 연산부(700)는 제1파워 스펙트럼 레벨에서 제2파워 스펙트럼 레벨을 감한 값을 산출할 수 있다.In this step, a process of calculating the difference value between the power spectrum level of the first frequency and the power spectrum level of the second frequency may be performed by the
6. 골구조 예측 단계<S406>6. Bone structure prediction step <S406>
본 단계에서 골구조 예측부(800)는 연산부(700)에서 산출된 값과 골구조 간의 상관관계를 이용하여 해면질골(20)의 골구조를 예측하는 과정이 진행될 수 있다. 즉, 골구조 예측부(800)는 단계 S405에서 연산부(700)가 산출한 차이값을 전술한 수학식 1의 X값에 대입함으로써, 실제 환자의 해면질골의 골 용적비율을 예측할 수 있다. 아울러, 골구조 예측부(800)는 단계 S405에서 연산부(700)가 산출한 차이값을 전술한 수학식 2의 X값에 대입함으로써, 실제 환자의 해면질골의 골소주 선밀도를 예측할 수 있다.In this step, the bone structure prediction unit 800 may proceed with a process of predicting the bone structure of the
아울러, 골구조 예측부(800)는 예측된 골소주 선밀도 또는 골 용적비율 값을 기설정된 기준값과 비교하여 환자의 골다공증 여부를 진단할 수 있다. 예를 들어, 환자의 골소주 선밀도 값이나 골 용적비율 값이 기준값 미만일 때, 골구조 예측부(800)는 환자의 골다공증이 진행 중이라고 판단할 수 있다.In addition, the bone structure prediction unit 800 may diagnose the patient's osteoporosis by comparing the predicted bone somatic linear density or bone volume ratio value with a predetermined reference value. For example, when the patient's bone mineral density or bone volume ratio value is less than the reference value, the bone structure prediction unit 800 may determine that the patient's osteoporosis is in progress.
7. 출력 단계<S407>7. Output step <S407>
골구조 예측부(800)에서 예측된 해면질골(20)의 골구조(즉, 골 용적비율과 골소주 선밀도) 결과값 또는 예측된 골구조 결과값에 따른 골다공증 진단 결과를 본 단계에서 출력부(900)가 화면에 출력하거나 사용자에게 음성으로 안내하는 과정이 이루어질 수 있다.The bone structure prediction unit 800 predicts the osteoporosis diagnosis result according to the bone structure (ie, bone volume ratio and bone strain linear density) result value or the predicted bone structure result value predicted by the
상술한 바와 같이, 본 발명의 일 실시예에 따른 초음파의 비선형 특성을 이용한 골구조 예측 장치 및 방법은 골구조의 예측시 사용되는 변수의 개수를 최소로 줄이고, 수식을 간소화시킴으로써 간단한 계산으로도 골구조를 예측할 수 있으므로 많은 변수와 복잡한 수식을 이용하여 골구조를 예측하는 기술에 비해 편리하고 손쉽게 골구조를 예측할 수 있다.As described above, the apparatus and method for predicting bone structure using nonlinear characteristics of ultrasonic waves according to an embodiment of the present invention can reduce the number of variables used in the prediction of the bone structure to a minimum and simplify the equation to achieve the goal even with simple calculation. Since the structure can be predicted, it is more convenient and easier to predict the bone structure than the technology for predicting the bone structure using many variables and complex formulas.
또한, 일 실시예에서는 골구조의 예측시 사용되는 변수로서, 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨만을 이용하며, 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 통해 간단하게 골구조를 예측할 수 있으므로 골구조의 예측시 사용되는 변수의 개수가 기존에 비해 대폭 감소되고, 복잡한 수식 또한 간소화되어 기존에 비해 측정 오차를 줄일 수 있고, 결과적으로 골구조 예측 정밀도가 향상되는 효과가 있다.In addition, in one embodiment, as a variable used in the prediction of the bone structure, only the power spectrum level of the first frequency and the power spectrum level of the second frequency are used, and the power spectrum level of the first frequency and the power spectrum level of the second frequency are used. Since the bone structure can be easily predicted through the difference value between, the number of variables used in the prediction of the bone structure is greatly reduced compared to the previous one, and the complex formula is also simplified, reducing the measurement error compared to the previous one, and consequently the bone structure. There is an effect that the prediction precision is improved.
아울러, 일 실시예에 따르면, 환자의 해면질골의 골구조를 파악하기 위해 마이크로 전산화 단층 촬영을 실시하지 않아도 해면질골의 골소주 선밀도나 골 용적비율 등의 골구조를 쉽게 예측할 수 있으므로 사용상 편의성이 높고, 골다공증의 진단이 신속하게 이루어지는 효과가 있다.In addition, according to one embodiment, it is easy to use because it is easy to predict the bone structure, such as the linear density of the spongy bone or bone volume ratio of the spongy bone, even without micro-computed tomography to determine the bone structure of the spongy bone of the patient, high convenience in use It is effective in quickly diagnosing osteoporosis.
위에서 설명한 바와 같이 본 발명에 대한 구체적인 설명은 첨부된 도면을 참조한 실시예에 의해서 이루어졌지만, 상술한 실시예는 본 발명의 바람직한 예를 들어 설명하였을 뿐이기 때문에, 본 발명이 상기의 실시예에만 국한되는 것으로 이해되어져서는 아니 되며, 본 발명의 권리범위는 후술하는 청구범위 및 그 균등개념으로 이해되어져야 할 것이다.As described above, the detailed description of the present invention has been made by the embodiments with reference to the accompanying drawings. However, since the above-described embodiments are merely described as preferred examples of the present invention, the present invention is limited to the above-described embodiments only. It should not be understood as being, and the scope of the present invention should be understood as the claims and the equivalent concept described later.
10 : 초음파의 비선형 특성을 이용한 골구조 예측 장치
100 : 하우징
110 : 매질
200 : 초음파 송신부
300 : 초음파 수신부
400 : 파형 발생부
500 : 신호처리부
600 : 변환부
700 : 연산부
800 : 골구조 예측부
900 : 출력부
D1 : 해면질골과 초음파 송신부 간의 직선거리
D2 : 해면질골과 초음파 수신부 간의 직선거리
20 : 해면질골10: Bone structure prediction device using nonlinear characteristics of ultrasound
100: housing
110: medium
200: ultrasonic transmitter
300: ultrasonic receiver
400: waveform generator
500: signal processing unit
600: conversion unit
700: operation unit
800: bone structure prediction unit
900: output unit
D1: straight line distance between the spongy bone and the ultrasound transmitter
D2: straight line distance between the spongy bone and the ultrasound receiver
20: spongy bone
Claims (10)
상기 해면질골을 통과한 초음파를 수신하고, 이를 전기적 신호로 변환하는 초음파 수신부;
상기 해면질골에 초음파가 입사되도록 제1주파수를 갖는 송신 초음파에 대한 전기적 신호를 생성하여 상기 초음파 송신부로 전송하는 파형 발생부;
상기 초음파 수신부에서 변환된 전기적 신호를 검출하는 신호처리부;
상기 신호처리부에서 검출된 전기적 신호를 변환함으로써, 상기 제1주파수 및 상기 제1주파수의 2배인 주파수인 제2주파수에 대한 파워 스펙트럼 레벨을 도출하는 변환부;
상기 제1주파수의 파워 스펙트럼 레벨과 상기 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 산출하는 연산부; 및
상기 연산부에서 산출된 값과 골구조 간의 상관관계를 이용하여 상기 해면질골의 골구조를 예측하는 골구조 예측부;를 포함하는 것을 특징으로 하는
초음파의 비선형 특성을 이용한 골구조 예측 장치.An ultrasonic transmitter that converts the applied electrical signal into ultrasonic waves and enters the spongy bone;
An ultrasonic wave receiving unit that receives ultrasonic waves that have passed through the spongy bone and converts them into electrical signals;
A waveform generator for generating an electrical signal for transmitting ultrasonic waves having a first frequency so that ultrasonic waves are incident on the spongy bone and transmitting them to the ultrasonic transmitter;
A signal processor that detects the electrical signal converted by the ultrasonic receiver;
A conversion unit that derives a power spectrum level for the first frequency and a second frequency that is twice the frequency of the first frequency by converting the electrical signal detected by the signal processing unit;
A calculating unit calculating a difference value between the power spectrum level of the first frequency and the power spectrum level of the second frequency; And
And a bone structure prediction unit for predicting the bone structure of the spongy bone by using a correlation between the value calculated by the operation unit and the bone structure.
A device for predicting bone structure using nonlinear characteristics of ultrasound.
상기 해면질골의 골구조는 골소주 선밀도를 포함하는 것을 특징으로 하는
초음파의 비선형 특성을 이용한 골구조 예측 장치.According to claim 1,
The bone structure of the spongy bone is characterized in that it comprises the bone density of the shochu
A device for predicting bone structure using nonlinear characteristics of ultrasound.
상기 해면질골의 골구조는 골 용적비율을 포함하는 것을 특징으로 하는
초음파의 비선형 특성을 이용한 골구조 예측 장치.According to claim 2,
The bone structure of the spongy bone is characterized in that it comprises a bone volume ratio
A device for predicting bone structure using nonlinear characteristics of ultrasound.
상기 골구조 예측부는 상기 연산부에서 산출된 값을 1차 선형 방정식에 대입함으로써, 상기 골소주 선밀도를 예측하되,
상기 1차 선형 방정식은, 해면질골 샘플에 초음파를 투과시켜 도출한 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 독립변수로 설정하고, 상기 해면질골 샘플의 골소주 선밀도를 종속변수로 설정하며, 상기 독립변수와 종속변수를 선형회귀분석함으로써 도출된 것을 특징으로 하는
초음파의 비선형 특성을 이용한 골구조 예측 장치.According to claim 2,
The bone structure prediction unit predicts the bone strain linear density by substituting a value calculated by the calculation unit into a linear linear equation,
In the first linear equation, the difference value between the power spectrum level of the first frequency and the power spectrum level of the second frequency derived by transmitting ultrasound to the spongy bone sample is set as an independent variable, and the linear density of the bone shochu of the spongy bone sample is determined. Set as a dependent variable, characterized by being derived by linear regression analysis of the independent variable and the dependent variable
A device for predicting bone structure using nonlinear characteristics of ultrasound.
상기 골구조 예측부는 상기 연산부에서 산출된 값을 1차 선형 방정식에 대입함으로써, 상기 골 용적비율을 예측하되,
상기 1차 선형 방정식은, 해면질골 샘플에 초음파를 투과시켜 도출한 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 독립변수로 설정하고, 상기 해면질골 샘플의 골 용적비율을 종속변수로 설정하며, 상기 독립변수와 종속변수를 선형회귀분석함으로써 도출된 것을 특징으로 하는
초음파의 비선형 특성을 이용한 골구조 예측 장치.According to claim 3,
The bone structure prediction unit predicts the bone volume ratio by substituting a value calculated by the calculation unit into a linear linear equation,
In the first linear equation, the difference value between the power spectrum level of the first frequency and the power spectrum level of the second frequency derived by transmitting ultrasound to the spongy bone sample is set as an independent variable, and the bone volume ratio of the spongy bone sample is Is set as a dependent variable, and the independent variable and the dependent variable are derived by linear regression analysis.
A device for predicting bone structure using nonlinear characteristics of ultrasound.
상기 해면질골은 상기 초음파 송신부 및 초음파 수신부의 사이에 위치하되, 전달 매질의 특성에 의해 발생하는 신호 왜곡을 최소화하도록 상기 해면질골과 초음파 수신부 간의 직선거리는 상기 해면질골과 초음파 송신부 간의 직선거리에 비해 상대적으로 짧게 형성된 것을 특징으로 하는
초음파의 비선형 특성을 이용한 골구조 예측 장치.According to claim 1,
The cavernous bone is positioned between the ultrasonic transmitter and the ultrasonic receiver, but the linear distance between the cavernous bone and the ultrasonic receiver is minimized to minimize signal distortion caused by the characteristics of the transmission medium, compared to the linear distance between the cavernous bone and the ultrasonic transmitter. Characterized in that the short formed
A device for predicting bone structure using nonlinear characteristics of ultrasound.
초음파 수신부가 상기 해면질골을 통과한 초음파를 수신하고, 이를 전기적 신호로 변환하는 신호 수신 단계;
신호처리부가 상기 초음파 수신부에서 변환된 전기적 신호를 검출하는 신호 검출 단계;
변환부가 상기 신호처리부에서 검출된 전기적 신호를 변환함으로써, 상기 제1주파수 및 상기 제1주파수의 2배인 주파수인 제2주파수에 대한 파워 스펙트럼 레벨을 도출하는 변환 단계;
연산부가 상기 제1주파수의 파워 스펙트럼 레벨과 상기 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 산출하는 연산 단계; 및
골구조 예측부가 상기 연산부에서 산출된 값과 골구조 간의 상관관계를 이용하여 상기 해면질골의 골구조를 예측하는 골구조 예측 단계;를 포함하는 것을 특징으로 하는
초음파의 비선형 특성을 이용한 골구조 예측 방법.When the waveform generator generates an electrical signal for transmitting ultrasound having a first frequency and transmits it to the ultrasonic transmitting unit, an ultrasonic incident step in which the ultrasonic transmitting unit injects ultrasonic waves into the spongy bone;
A signal receiving step of receiving an ultrasonic wave through the spongy bone and converting it into an electrical signal;
A signal detection step in which a signal processing unit detects an electrical signal converted by the ultrasonic receiving unit;
A conversion step of converting an electrical signal detected by the signal processing unit to derive a power spectrum level for the first frequency and a second frequency that is twice the frequency of the first frequency;
A calculating step of calculating a difference between the power spectrum level of the first frequency and the power spectrum level of the second frequency; And
And a bone structure prediction step of predicting the bone structure of the spongy bone by using a correlation between the bone structure prediction unit and the value calculated by the operation unit.
A method for predicting bone structure using nonlinear characteristics of ultrasound.
상기 골구조 예측 단계에서 상기 해면질골의 골구조는 골소주 선밀도를 포함하고, 상기 골구조 예측부는 상기 연산부에서 산출된 값을 1차 선형 방정식에 대입함으로써, 상기 골소주 선밀도를 예측하되,
상기 1차 선형 방정식은, 해면질골 샘플에 초음파를 투과시켜 도출한 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 독립변수로 설정하고, 상기 해면질골 샘플의 골소주 선밀도를 종속변수로 설정하며, 상기 독립변수와 종속변수를 선형회귀분석함으로써 도출된 것을 특징으로 하는
초음파의 비선형 특성을 이용한 골구조 예측 방법.The method of claim 7,
In the bone structure prediction step, the bone structure of the spongy bone includes bone soju line density, and the bone structure predictor predicts the bone soju line density by substituting a value calculated by the calculation unit into a linear linear equation,
In the first linear equation, the difference value between the power spectrum level of the first frequency and the power spectrum level of the second frequency derived by transmitting ultrasound to the spongy bone sample is set as an independent variable, and the linear density of the bone shochu of the spongy bone sample is determined. Set as a dependent variable, characterized by being derived by linear regression analysis of the independent variable and the dependent variable
A method for predicting bone structure using nonlinear characteristics of ultrasound.
상기 골구조 예측 단계에서 상기 해면질골의 골구조는 골 용적비율을 포함하고, 상기 골구조 예측부는 상기 연산부에서 산출된 값을 1차 선형 방정식에 대입함으로써, 상기 골 용적비율을 예측하되,
상기 1차 선형 방정식은, 해면질골 샘플에 초음파를 투과시켜 도출한 제1주파수의 파워 스펙트럼 레벨과 제2주파수의 파워 스펙트럼 레벨 간의 차이값을 독립변수로 설정하고, 상기 해면질골 샘플의 골 용적비율을 종속변수로 설정하며, 상기 독립변수와 종속변수를 선형회귀분석함으로써 도출된 것을 특징으로 하는
초음파의 비선형 특성을 이용한 골구조 예측 방법.The method of claim 7,
In the bone structure prediction step, the bone structure of the spongy bone includes a bone volume ratio, and the bone structure prediction unit predicts the bone volume ratio by substituting a value calculated by the calculation unit into a linear linear equation,
In the first linear equation, the difference value between the power spectrum level of the first frequency and the power spectrum level of the second frequency derived by transmitting ultrasound to the spongy bone sample is set as an independent variable, and the bone volume ratio of the spongy bone sample is Is set as a dependent variable, and the independent variable and the dependent variable are derived by linear regression analysis.
A method for predicting bone structure using nonlinear characteristics of ultrasound.
상기 초음파 입사 단계 및 신호 수신 단계에서 상기 해면질골은 상기 초음파 송신부 및 초음파 수신부의 사이에 위치하되, 전달 매질의 특성에 의해 발생하는 신호 왜곡을 최소화하도록 상기 해면질골과 초음파 수신부 간의 직선거리는 상기 해면질골과 초음파 송신부 간의 직선거리에 비해 상대적으로 짧게 형성된 것을 특징으로 하는
초음파의 비선형 특성을 이용한 골구조 예측 방법.The method of claim 7,
In the ultrasonic incidence step and the signal receiving step, the spongy bone is located between the ultrasound transmitter and the ultrasound receiver, but the straight line distance between the spongy bone and the ultrasound receiver is minimized to minimize signal distortion caused by characteristics of the transmission medium. It is characterized in that it is formed relatively short compared to the linear distance between the ultrasonic transmitter and the
A method for predicting bone structure using nonlinear characteristics of ultrasound.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020180111553A KR102159012B1 (en) | 2018-09-18 | 2018-09-18 | Apparatus and method for estimating bone structure using ultrasonic nonlinearity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020180111553A KR102159012B1 (en) | 2018-09-18 | 2018-09-18 | Apparatus and method for estimating bone structure using ultrasonic nonlinearity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| KR20200032477A true KR20200032477A (en) | 2020-03-26 |
| KR102159012B1 KR102159012B1 (en) | 2020-09-24 |
Family
ID=69958604
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| KR1020180111553A KR102159012B1 (en) | 2018-09-18 | 2018-09-18 | Apparatus and method for estimating bone structure using ultrasonic nonlinearity |
Country Status (1)
| Country | Link |
|---|---|
| KR (1) | KR102159012B1 (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20050042542A (en) * | 2003-11-03 | 2005-05-10 | 학교법인 성균관대학 | Nondestructive acoustic evaluation device and method by using nonlinear acoustic responses |
| KR101102784B1 (en) * | 2010-07-07 | 2012-01-05 | 강원대학교산학협력단 | Estimation method of bone mineral density |
| KR20140045124A (en) * | 2012-10-08 | 2014-04-16 | 강원대학교산학협력단 | Method and apparatus for estimating bone mineral density using ultrasonic nonlinear parameter |
| KR20160049642A (en) * | 2014-10-28 | 2016-05-10 | 강원대학교산학협력단 | Method and apparatus for estimating bone mineral density using a time reversal acoustic focusing technique |
| JP6189227B2 (en) * | 2013-02-20 | 2017-08-30 | 株式会社東芝 | Ultrasonic flaw detector and its evaluation method |
| KR101840349B1 (en) * | 2016-11-15 | 2018-03-21 | 강원대학교산학협력단 | Apparatus and method for estimating bone mineral density using ultrasonic sum frequency component |
| KR20180037351A (en) * | 2016-10-04 | 2018-04-12 | 강원대학교산학협력단 | Apparatus and method for estimating bone structure using ultrasonic nonlinear parameter |
-
2018
- 2018-09-18 KR KR1020180111553A patent/KR102159012B1/en active IP Right Grant
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20050042542A (en) * | 2003-11-03 | 2005-05-10 | 학교법인 성균관대학 | Nondestructive acoustic evaluation device and method by using nonlinear acoustic responses |
| KR101102784B1 (en) * | 2010-07-07 | 2012-01-05 | 강원대학교산학협력단 | Estimation method of bone mineral density |
| KR20140045124A (en) * | 2012-10-08 | 2014-04-16 | 강원대학교산학협력단 | Method and apparatus for estimating bone mineral density using ultrasonic nonlinear parameter |
| JP6189227B2 (en) * | 2013-02-20 | 2017-08-30 | 株式会社東芝 | Ultrasonic flaw detector and its evaluation method |
| KR20160049642A (en) * | 2014-10-28 | 2016-05-10 | 강원대학교산학협력단 | Method and apparatus for estimating bone mineral density using a time reversal acoustic focusing technique |
| KR20180037351A (en) * | 2016-10-04 | 2018-04-12 | 강원대학교산학협력단 | Apparatus and method for estimating bone structure using ultrasonic nonlinear parameter |
| KR101840349B1 (en) * | 2016-11-15 | 2018-03-21 | 강원대학교산학협력단 | Apparatus and method for estimating bone mineral density using ultrasonic sum frequency component |
Also Published As
| Publication number | Publication date |
|---|---|
| KR102159012B1 (en) | 2020-09-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4043181A (en) | Ultrasonic pulse-echo apparatus | |
| Vallet et al. | In vivo characterization of cortical bone using guided waves measured by axial transmission | |
| US8419643B2 (en) | Ultrasonic method and apparatus for assessment of bone | |
| JPH01503199A (en) | Ultrasonic evaluation method and device for evaluating bone condition in living organisms | |
| KR20140035932A (en) | Ultrasound apparatus for assessing the quality of a patient's bone tissue | |
| CN104013419A (en) | Adaptive acoustic pressure estimation in medical ultrasound | |
| KR20180037351A (en) | Apparatus and method for estimating bone structure using ultrasonic nonlinear parameter | |
| JPH08280671A (en) | Ultrasonic bon evaluating apparatus | |
| KR102140538B1 (en) | Apparatus for estimating bone mineral density and bone structure using frequency dependence of ultrasonic backscatter coefficient | |
| KR101412785B1 (en) | Method and apparatus for estimating bone mineral density using ultrasonic nonlinear parameter | |
| KR101345870B1 (en) | Method for estimating bone mineral density using ultrasonic attenuation | |
| Litniewski et al. | Ultrasonic scanner for in vivo measurement of cancellous bone properties from backscattered data | |
| KR101432871B1 (en) | Measuring method and device of bone density by using dispersion rate of ultrasonic phase velocity | |
| KR101102784B1 (en) | Estimation method of bone mineral density | |
| Litniewski et al. | Detection of bone disease with ultrasound—comparison with bone densitometry | |
| KR101840349B1 (en) | Apparatus and method for estimating bone mineral density using ultrasonic sum frequency component | |
| KR102159012B1 (en) | Apparatus and method for estimating bone structure using ultrasonic nonlinearity | |
| KR102303922B1 (en) | Method for estimating bone mineral density and bone structure using ultrasonic attenuation coefficient and phase velocity | |
| KR101510525B1 (en) | Method for Estimating Bone Structure of Proximal Femur by Using Ultrasonic Attenuation and Backscatter Coefficients | |
| KR102364309B1 (en) | Method for estimating bone mineral density and bone structure using ultrasonic attenuation coefficient and backscatter coefficient | |
| CN112244780A (en) | Bone density measuring device and method based on photoacoustic signals | |
| KR100581229B1 (en) | Method for Measuring the Density of Shinbone by Using Lamb Wave | |
| JP3472376B2 (en) | Ultrasonic bone evaluation device | |
| KR20020082378A (en) | A osteoporosis apparatus and method by ultrasound longitudinal transmission | |
| KR101643101B1 (en) | Method And Apparatus For Estimating Bone Mineral Density Using Bidirectional Ultrasound Transverse Transmission Method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| E902 | Notification of reason for refusal | ||
| E701 | Decision to grant or registration of patent right |