KR20200102028A - Method for converting scan information of computed tomography scanner into bone parameters - Google Patents
Method for converting scan information of computed tomography scanner into bone parameters Download PDFInfo
- Publication number
- KR20200102028A KR20200102028A KR1020190019767A KR20190019767A KR20200102028A KR 20200102028 A KR20200102028 A KR 20200102028A KR 1020190019767 A KR1020190019767 A KR 1020190019767A KR 20190019767 A KR20190019767 A KR 20190019767A KR 20200102028 A KR20200102028 A KR 20200102028A
- Authority
- KR
- South Korea
- Prior art keywords
- computed tomography
- scan information
- tomography apparatus
- phantoms
- subject
- Prior art date
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
- A61B6/505—Clinical applications involving diagnosis of bone
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
본 발명은 피검체의 골격 파라미터 취득방법에 관한 것으로서, 특히 컴퓨터 단층 촬영장치를 통한 스캔 정보를 골격 파라미터로 전환하는 방법에 관한 것이다.The present invention relates to a method of acquiring skeleton parameters of a subject, and in particular, to a method of converting scan information through a computed tomography apparatus into skeleton parameters.
골밀도(bone mineral density) 검사는 통상적으로 특수한 장치를 통해 실시된다. 예를 들어 이중에너지 X선 골밀도 검사(Dual-energy X-ray absorptiometry, DEXA) 및 정량적 컴퓨터 단층 촬영 검사(Quantitative computed tomography, QCT) 등이 있다. Bone mineral density testing is usually carried out using a special device. Examples include dual-energy X-ray absorptiometry (DEXA) and quantitative computed tomography (QCT).
이중에너지 X선 골밀도 검사 기술은 두 종류의 상이한 에너지의 X선을 이용하여 피검 조직을 스캔하고, 피검 조직을 관통한 X선이 섬광 검출기에 의해 수신된 후 컴퓨터 처리를 거쳐 골밀도를 획득하는 방식이다.The dual-energy X-ray bone density test technology is a method of scanning a tissue to be examined using two types of X-rays of different energies, and obtaining bone density through computer processing after X-rays penetrating the tissue are received by a scintillation detector. .
정량적 컴퓨터 단층 촬영 검사는 실무적으로 인체 모형을 포함하거나 또는 인체 모형을 포함하지 않는 두 가지 유형으로 구분된다. 모형을 이용한 QCT는 피검 조직을 밀도가 이미 알려진 팬텀(phantom)과 동기적으로 스캔하여, 팬텀의 스캔 결과를 참고로, 피검조직의 스캔 정보를 팬텀의 스캔 정보와 비교하여 골밀도를 획득한다. 모형을 포함하지 않는 QCT는 근육과 지방 조직의 스캔 정보를 참고로 하여, 피검 조직의 스캔 정보와 비교함으로써 골밀도를 획득한다. QCT의 단점은 방사량이 크고 스캔 시간이 비교적 길며, 측정의 정확성이 상대적으로 낮다는데 있다.Quantitative computed tomography examinations are practically divided into two types, with or without manikin. The QCT using the model scans a tissue to be examined synchronously with a phantom whose density is known, and obtains bone density by comparing scan information of the tissue to be examined with the scan information of the phantom with reference to the scan result of the phantom. QCT without a model obtains bone density by referring to scan information of muscle and adipose tissue and comparing it with scan information of the tissue to be examined. The disadvantages of QCT are that the radiation amount is large, the scan time is relatively long, and the measurement accuracy is relatively low.
이로써 알 수 있듯이, 골밀도는 현재 모두 특수 또는 전용 장치를 이용하여 검사가 실시된다.As can be seen from this, all bone mineral density is currently tested using a special or dedicated device.
또한, 비록 컴퓨터 단층 촬영장치(Computed Tomography, CT)는 이미 생체 조직 검사에 광범위하게 응용되고 있으나, 현재 컴퓨터 단층 촬영장치는 피검 조직의 삼차원 영상만 구축할 수 있고 골격 파라미터는 판별할 수 없다.In addition, although computed tomography (CT) is already widely applied to biological tissue examination, the current computed tomography can only construct a three-dimensional image of a tissue to be examined and cannot determine skeletal parameters.
상기 결점을 감안하여, 본 발명은 상용(시판)되는 컴퓨터 단층 촬영장치로 촬영한 스캔 정보를 통해, 컴퓨터 단층 촬영장치 자체의 에너지 감쇠계수를 획득하여 스캔 정보를 피검체가 갖고 있는 골격 파라미터로 전환함으로써 골격 검사 능력을 향상시킬 수 있다.In view of the above drawbacks, the present invention obtains the energy attenuation coefficient of the computed tomography apparatus itself through scan information taken by a commercial (commercially available) computed tomography apparatus, and converts the scan information into a skeleton parameter possessed by the subject. By doing so, it is possible to improve the skeletal examination ability.
상기 목적을 달성하기 위하여, 본 발명인 컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법은 컴퓨터 단층 촬영장치를 제공하는 단계; 피검체 및 성분이 이미 알려진 두 팬텀을 제공하는 단계; 컴퓨터 단층 촬영장치를 통해 피검체 및 두 팬텀을 촬영하여 대응되는 피검체 스캔 정보 및 두 팬텀의 스캔 정보를 획득하는 단계; 계산장치를 통해 피검체의 스캔 정보 및 두 팬텀의 스캔 정보를 수신하는 단계; 계산장치가 두 팬텀에 대응하며 두 팬텀의 기지의 성분 및 두 팬텀의 스캔 정보를 포함하는 물리함수모델을 통해 컴퓨터 단층 촬영장치의 에너지 감쇠계수를 계산하는 단계; 두 팬텀의 기지의 성분의 에너지 감쇠계수 및 이상적인 피검체의 에너지 감쇠계수와 관련이 있는 에너지 보정계수를 계산장치로 제공하는 단계; 및 계산장치가 컴퓨터 단층 촬영장치의 에너지 감쇠계수, 피검체 스캔 정보 및 에너지 보정 계수를 포함하는 진실관계함수를 통해 피검체의 골격 파라미터를 획득하는 단계를 포함한다.In order to achieve the above object, a method of converting scan information of a computed tomography apparatus according to the present invention into a skeleton parameter comprises: providing a computed tomography apparatus; Providing two phantoms in which the subject and components are known; Acquiring scan information of the subject and scan information of the two phantoms by photographing the subject and the two phantoms through a computed tomography apparatus; Receiving scan information of a subject and scan information of two phantoms through a calculation device; Calculating, by the computing device, an energy attenuation coefficient of the computed tomography apparatus through a physical function model corresponding to the two phantoms and including a known component of the two phantoms and scan information of the two phantoms; Providing an energy attenuation coefficient of a known component of the two phantoms and an energy correction coefficient related to an energy attenuation coefficient of an ideal subject to a calculation device; And obtaining, by the computing device, a skeleton parameter of the subject through a truth relationship function including an energy attenuation coefficient of the computed tomography apparatus, scan information of the subject, and an energy correction coefficient.
이와 같이, 본 발명은 성분이 이미 알려진 두 팬텀 및 대응되는 두 팬텀의 스캔 정보를 통해 컴퓨터 단층 촬영장치의 에너지 감쇠계수를 획득하고, 에너지 감쇠계수와 에너지 보정계수를 통해 피검체 스캔 정보의 골격 파라미터를 획득할 수 있다.In this way, the present invention obtains the energy attenuation coefficient of the computed tomography apparatus through scan information of two phantoms with known components and corresponding two phantoms, and the skeleton parameter of the subject scan information through the energy attenuation coefficient and the energy correction coefficient. Can be obtained.
본 발명이 제공하는 컴퓨터 단층촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법의 조성, 특징 또는 방법에 관하여, 후속되는 실시방식에서 상세히 설명할 것이다. 그러나 본 발명 분야에서 통상적인 지식을 갖춘 자라면 상기 상세한 설명 및 본 발명에 열거된 특정 실시예의 실시는 단지 본 발명을 설명하기 위한 것일 뿐, 본 발명의 특허출원범위를 제한하기 위한 것이 아님을 이해할 수 있을 것이다.The composition, feature or method of a method for converting scan information of a computed tomography apparatus provided by the present invention into a skeleton parameter will be described in detail in a subsequent implementation manner. However, those with ordinary knowledge in the field of the present invention will understand that the detailed description and implementation of the specific embodiments listed in the present invention are only for describing the present invention and are not intended to limit the scope of the present invention. I will be able to.
도 1은 본 발명인 컴퓨터 단층촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법의 단계 흐름도이다.
도 2는 도 1을 실행하는 전환 시스템의 설명도이다.1 is a step flow diagram of a method for converting scan information of a computed tomography apparatus according to the present invention into a skeleton parameter.
Fig. 2 is an explanatory diagram of a switching system for executing Fig. 1;
이하, 각 도면을 결합하여 대응하는 바람직한 실시예를 열거함으로써 본 발명인 컴퓨터 단층 촬영장치의 스캔정보를 골격 파라미터로 전환하는 방법의 기술특징 및 달성되는 효과에 대해 설명한다. 그러나 각 도면 중 컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법의 단계 및 단계의 수량은 단지 본 발명의 기술특징을 설명하기 위한 것일 뿐, 본 발명에 대해 제한을 구성하지 않는다.Hereinafter, technical features and achieved effects of a method of converting scan information of the computed tomography apparatus of the present invention into skeleton parameters by enumerating corresponding preferred embodiments by combining the drawings will be described. However, the steps and the number of steps in the method of converting scan information of a computed tomography apparatus into a skeleton parameter in each drawing are only for describing the technical features of the present invention, and do not constitute a limitation on the present invention.
도 1에 도시된 바와 같이, 본 발명의 컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법(100)은 상기 도면 중의 단계를 통해 구현되며, 단 기타 실시예에서, 단계의 수량은 더 많거나 또는 더 적을 수도 있고, 단계의 순서는 조정 가능하며, 따라서 단계의 수량 및 순서는 본 실시예로 한정되지 않는다.As shown in Fig. 1, the
단계 110은 두 팬텀을 제공하는 단계이고, 단계 111은 두 팬텀을 분석하는 단계이다. 그 중, 분석단계는 기지(旣知)의 성분 분석 방식을 통해 팬텀과 피검체의 성분을 획득할 수 있다. 단계 113은 피검체를 제공하는 단계이며, 본 실시예에서, 피검체는 생체의 골격(cortical and marrow)이다. 본 실시예의 기지의 성분의 상기 파라미터 또는 비율은 설명용이며, 본 발명을 제한하기 위한 것이 아니다.
본 실시예에서, 팬텀의 성분은 인산수소칼륨(K2HPO4) 및 물(H2O)의 혼합물이다. 팬텀의 제작은 특정 비율의 인산수소칼륨과 물을 혼합한 후 몰드캐비티에 주입하고, 경화 후 다분절(단) 골격 팬텀을 형성한다. 팬텀의 각 분절(단)에 충전되는 용액은 당량 팬텀 밀도가 다름을 나타내며, 당량 팬텀 밀도 범위는 0.1 내지 0.9g/cm3 사이이다. 팬텀은 다수의 상이한 영역의 세퍼레이터를 구비하며, 각 세퍼레이터 내의 팬텀 파라미터 비율은 상이하고, 또한 동일한 세퍼레이터 내의 재질 비율은 평균치 상태이어야 한다.In this example, the component of the phantom is a mixture of potassium hydrogen phosphate (K 2 HPO 4 ) and water (H 2 O). In the manufacture of the phantom, a certain ratio of potassium hydrogen phosphate and water is mixed, then injected into the mold cavity, and after curing, a multi-segmented (single) skeleton phantom is formed. The solution charged in each segment (end) of the phantom indicates that the equivalent phantom density is different, and the equivalent phantom density range is between 0.1 and 0.9 g/cm 3 . The phantom includes a plurality of separators in different regions, the phantom parameter ratio in each separator is different, and the material ratio in the same separator should be in an average state.
상기 두 팬텀의 비율은 상이하며, 이어서, 두 팬텀을 각각 제1 팬텀 및 제2 팬텀으로 칭한다. 제1 팬텀과 제2 팬텀의 비율은 하기 표 1과 같으며, V1과 V2는 각각 부피백분율의 수치를 나타내고, 또한 양자는 상이하다. W1과 W2는 각각 특정 중량백분율의 수치를 나타내며, 또한 양자는 상이하다. 과 역시 각각 상이한 밀도의 제1 팬텀 및 제2 팬텀을 나타내며, 따라서 팬텀은 표 1에 따라 상이한 성분과 비율의 팬텀으로 제조될 수 있다.The ratios of the two phantoms are different, and then, the two phantoms are referred to as a first phantom and a second phantom, respectively. The ratios of the first phantom and the second phantom are shown in Table 1 below, and V1 and V2 each represent a value of a volume percentage, and both are different. W1 and W2 each represent a value of a specific weight percentage, and both are different. and Also, each of the first and second phantoms of different densities is represented, and thus the phantom may be manufactured with different components and ratios of phantoms according to Table 1.
표 2는 표 1의 설명에 따라 획득된 제1 팬텀, 제2 팬텀 및 물의 성분의 부피백분율, 중량백분율, 밀도 및 스캔 정보의 관련 수치이며, 상기 표 중의 수치는 단지 설명용일 뿐, 본 발명을 제한하는 것이 아니다.Table 2 is the related values of the volume percentage, weight percentage, density and scan information of the components of the first phantom, the second phantom, and water obtained according to the description of Table 1, and the numerical values in the table are for illustration only, and the present invention is It is not limiting.
팬텀 분석은 표 3과 같이, 원소주기율표를 통해 K2HPO4 및 H2O의 각 원소(K, H, P 및 O)의 원자번호()와 원자량()을 획득할 수 있으며, 이후 표 1을 통해 팬텀에서 차지하는 각 원소의 중량백분율() 및 밀도()를 획득한다.As shown in Table 3, phantom analysis is performed using K 2 HPO 4 through the periodic table of the elements. And the atomic number of each element (K, H, P and O) of H 2 O ( ) And atomic weight ( ) Can be obtained, and the weight percentage of each element occupied by the phantom ( ) And density ( ).
단계 130은 컴퓨터 단층 촬영장치(Computed Tomography; 약칭 CT)를 제공하는 단계로서, 컴퓨터 단층 촬영장치는 디지털 기하 계산을 거친 후 입체 방사선 의료 영상을 구축한다. 컴퓨터 단층 촬영장치는 광자 콜리메이터를 구비한 임의의 장치일 수 있으며, 예를 들어 멀티디텍터 CT와 이중에너지 CT 등이다. 본 실시예에서는 단일에너지 계수를 이용하여 교정하므로, 스캔 횟수 및 방사선 조사량을 감소시킬 수 있으나, 단 기타 실시예에서, 컴퓨터 단층 촬영장치는 이중에너지 조사를 채택할 수도 있다.In
단계 131은 제1 팬텀, 제2 팬텀 및 피검체를 촬영하여 대응되는 제1 팬텀 스캔 정보(CT number), 제2 팬텀 스캔 정보 및 피검체 스캔 정보를 획득하는 단계이다. 스캔 정보는 CT값(hounsfield unit, HU)이라고도 칭하며, CT값()은 생체(피검체)의 국부 조직을 측정 시의 계량단위로서, 광감쇠값을 나타내기도 한다.
단계 140은 피검체의 스캔 정보, 제1 팬텀의 스캔 정보 및 제2 팬텀의 스캔 정보를 수신하는 단계이다.
단계 150은 컴퓨터 단층 촬영장치의 에너지 감쇠계수를 계산하는 단계로서, 에너지 감쇠계수는 제1 팬텀과 제2 팬텀의 성분, 제1 팬텀의 스캔 정보, 제2 팬텀의 스캔 정보, 제1 팬텀에 대응하는 물리함수모델 및 제2 팬텀에 대응하는 물리함수모델을 근거로 계산하여 획득된다. In
단계 150에서 지칭하는 물리함수모델은 수학식(1)과 수학식(2)을 포함한다. 그 중 수학식(1)은 광감쇠값 을 계산하기 위한 것이고, 또한 의 은 임의의 숫자 또는 부호를 통해 각각 제1 팬텀, 제2 팬텀 및 물 등을 나타내는 것이며, 는 아보가르도 상수이고, 은 밀도이며, 는 중량백분율이고, 는 원자량이며, , , 는 원자번호 상관계수이고, , 는 각각 컴퓨터 단층 촬영장치의 에너지 감쇠계수이며, 에너지 감쇠계수는 레일리산란계수인 및 광전흡수계수 를 포함한다. 아보가르도 상수는 6.022E23이다.The physical function model referred to in
수학식(2)의 등호 좌변은 CT값을 나타내고, 우변은 팬텀(피검)체의 광감쇠값을 물로 나눈 광감쇠값 이다. 수학식(2)의 등호 좌변이 나타내는 CT값은 각각 컴퓨터 단층 촬영장치의 촬영을 통해 획득된 각 팬텀의 CT값()이고, CT값은 수학식(2)를 통해 전환된 후 팬텀(피검)체의 광감쇠값 을 물로 나눈 광감쇠값 과 동등한 상수값이 된다.The left side of the equal sign in Equation (2) represents the CT value, and the right side is the light attenuation value of the phantom (subject). Light attenuation value divided by water to be. The CT values indicated by the left side of the equal sign in Equation (2) are the CT values of each phantom obtained through imaging by a computed tomography apparatus ( ), and the CT value is the light attenuation value of the phantom (subject) after being converted through Equation (2) Light attenuation value divided by water It becomes a constant value equal to
수학식(1) 중 은 i번째 원소의 원자번호로서, 예를 들어 H(수소)는 1이고, 산소(O)는 8이며, 참고문헌 Attix, F. H. (2008). Introduction to radiological physics and radiation dosimetry. John Wiley & Sons.과 같이, 는 광전흡수의 유효 원자번호이고, 는 레일리산란의 유효 원자번호이다. 이밖에, 광전흡수의 유효원자번호 중의 와 레일리산란의 유효원자번호의 는 Rutherford, R. A., Pullan, B. R., & Isherwood, I. (1976). Measurement of effective atomic number and electron density using an EMI scanner. Neuroradiology, 11(1), 15-21. Rutherford를 참조할 수 있다. 본 실시예에서, 는 1.86이고, 는 3.62이다. 산소를 예로 들면, 는 14868.79384이고, 는 382.6814077이다. 관련 파라미터는 표 4와 같으며, 상기 파라미터들은 설명용 일뿐, 본 발명에 대해 제한을 구성하지 않는다.In Equation (1) Is the atomic number of the i-th element, for example, H (hydrogen) is 1, oxygen (O) is 8, and references Attix, FH (2008). Introduction to radiological physics and radiation dosimetry. Like John Wiley & Sons., Is the effective atomic number of photoelectric absorption, Is the effective atomic number of Rayleigh scattering. In addition, among the effective atomic numbers of photoelectric absorption, And the effective atomic number of Rayleigh scattering Rutherford, RA, Pullan, BR, & Isherwood, I. (1976). Measurement of effective atomic number and electron density using an EMI scanner. Neuroradiology, 11(1), 15-21. See Rutherford. In this example, Is 1.86, Is 3.62. Take oxygen as an example, Is 14868.79384, Is 382.6814077. Related parameters are shown in Table 4, and the above parameters are for illustrative purposes only and do not constitute a limitation on the present invention.
표 4에서의 각 파라미터값을 수학식(1)에 대입하여 제1 팬텀의 광감쇠값 , 제2 팬텀의 광감쇠값 및 물의 광감쇠값 을 계산하며, 수학식(1) 중의 , , , , 원자번호(), 원자량() 및 중량백분율()은 물리성분을 통해 대응되는 상수값을 획득할 수 있다. 따라서, 수학식(1) 중 에너지 감쇠계수 , 만 미지수이다.The light attenuation value of the first phantom by substituting each parameter value in Table 4 into Equation (1) , Light attenuation value of the second phantom And light attenuation value of water And in Equation (1) , , , , Atomic number ( ), atomic weight ( ) And weight percentage ( ) Can obtain a corresponding constant value through the physical component. Therefore, the energy attenuation coefficient in Equation (1) , Only is unknown.
이후, 단계 150에서의 수학식(1)로 획득된 제1 팬텀의 광감쇠값 , 제2 팬텀의 광감쇠값 및 물의 광감쇠값 을 수학식(2)에 대입하고, 표 2에서의 제1 팬텀 및 제2 팬텀의 CT값을 매칭시켜 제1 팬텀의 광감쇠값 , 제2 팬텀의 광감쇠값 및 물의 광감쇠값 의 2개의 비율 수학식을 구축한다. 따라서, 두 비율 수학식의 해를 구하면 컴퓨터 단층 촬영장치의 에너지 감쇠계수 , 를 획득할 수 있으며, 본 실시예에서, 에너지 감쇠계수 , 는 각각 -4.960584907와 0.010816298이다.Thereafter, the light attenuation value of the first phantom obtained by Equation (1) in
따라서, 단계150의 설명을 통해 알 수 있듯이, 팬텀의 수량은 최소 2개면 되나, 단 팬텀의 수량은 2개로 한정되지 않고 더 많을 수도 있다.Therefore, as can be seen through the description of
이어서, 단계 170은 에너지 보정계수 X를 제공하는 단계로서, 에너지 보정계수 X는 이상적인 골격과 팬텀의 특정 성분의 에너지 의존성 차이를 통해 보정한다. 이상적인 골격의 에너지 감쇠계수는 미국 국립 표준기술연구소(국립표준기술연구소, NIST)의 광자 단면적 데이터베이스(XCOM:Photon Cross Sections Database)와 국제 방사선 단위 측정 위원회(International Commission on Radiation Units and Measurements, ICRU) 중의 조직 대체물 방사량 측정 및 측정 보고서(Tissue Substitutes in Radiation Dosimetry and Measurement Report-44)로부터 관련 물리량, 방사선 단위와 방사성의 국제적 허용 건의, 및 통합 계량 프로그램과 물리 데이터를 조회할 수 있다. 에너지 보정계수 X는 팬텀의 주 성분(인산수소칼륨)의 에너지 감쇠계수를 이상적인 골격으로 나눈 에너지 감쇠계수의 비율이다.Subsequently, in
단계 190은 진실관계함수를 구축하는 단계 191과, 피검체 골격 파라미터를 획득하는 단계 193을 포함한다. 다시 말해 단계 190은 진실관계함수를 통해 상기 피검체의 골격 파라미터를 획득한다. 진실관계함수는 수학식(3)과 같이, 피검체의 광감쇠값과 컴퓨터 단층 촬영장치로 촬영한 CT값()의 대응관계를 획득하며, 수학식(3)은 컴퓨터 단층 촬영장치의 에너지 감쇠계수 , , 피검체의 스캔정보() 및 에너지 보정계수 X와 관련이 있다. 수학식(3)이 수학식(2)와 다른 점은 수학식(3)은 에너지 보정계수를 통해 보정을 수행함으로써 피검체의 광감쇠값 획득의 정확성을 높인다는데 있다. 이와 같이, 광감쇠값 는 수학식(1)을 통해 골격 파라미터, 예를 들어 골격 밀도 를 획득할 수 있다.Step 190 includes
이와 같이, 상용(시판)되는 컴퓨터 단층촬영장치의 에너지 감쇠계수는 모두 차이가 있고, 또한 고정상수가 아니기 때문에, 단계 150을 통해 컴퓨터 단층 촬영장치의 에너지 감쇠계수 , 를 획득할 수 있으며, 따라서, 상용(시판) 컴퓨터 단층 촬영장치는 상기 단계를 통해 상관 에너지 감쇠계수를 획득할 수 있다. 이후, 다시 단계 170을 통해 에너지 보정계수 X를 획득함으로써, 단계 190에서 구축된 진실관계함수가 더욱 정확해질 수 있다.As described above, since the energy attenuation coefficients of commercial (commercially available) computed tomography apparatuses are all different and are not fixed constants, the energy attenuation coefficient of the computed tomography apparatuses through
본 발명은 컴퓨터 단층 촬영장치를 통해 촬영한 CT값을 상기 단계를 통해 전환하여 생체 골격의 밀도 파라미터를 획득할 수 있으며, 따라서 수술 계획 및 품질 향상에 도움이 된다. 예를 들어 골나사 작업 시, 의료인이 골격의 밀도 파라미터를 통해 골격의 각 분절(단)의 밀도 파라미터를 이해할 수 있어, 골밀도가 비교적 양호한 고정위치를 용이하게 선택하여 골나사를 고정시킬 수 있다. The present invention can obtain a density parameter of a living body skeleton by converting a CT value photographed through a computed tomography apparatus through the above steps, thus helping to improve operation planning and quality. For example, during bone screw work, a medical practitioner can understand the density parameter of each segment (end) of the skeleton through the density parameter of the skeleton, so that the bone screw can be fixed by easily selecting a fixed position with relatively good bone density.
단계 150의 계산은 컴퓨터 단층 촬영장치에 내장된 계산장치 또는 외부의 계산장치를 통해 수행되며, 다시 말해, 컴퓨터 단층 촬영장치에 내장된 계산장치 또는 외부의 계산장치를 막론하고, 수학식(1)과 수학식(2)는 계산장치가 처리 및 연산을 수행하도록 계산장치 내에 구축될 수 있다.The calculation of
단계 170은 입력 또는 미리 설정된 프로그램을 통해 에너지 보정계수 X를 계산장치로 제공하는 단계이다. 단계 190의 진실관계함수 구축은 진실관계함수를 컴퓨터 단층 촬영장치에 내장된 계산장치 또는 외부의 계산장치에 구축하는 것이며, 이와 같이, 컴퓨터 단층 촬영장치의 에너지 감쇠계수 , 에 대응되는 컴퓨터 단층 촬영장치가 촬영한 목표 골격의 CT값(정보)을 이미 알고 있다면, 대응되는 진실관계함수를 통해 목표 골격의 골격밀도 파라미터로 전환할 수 있다.Step 170 is a step of providing the energy correction coefficient X to the calculation device through input or a preset program. In
계산장치는 컴퓨터 단층 촬영장치가 스캔한 CT값(의료영상 또는 스캔정보라고도 칭함)을 피검체의 골격 파라미터(예를 들어 밀도)를 포함하는 의료영상 또는 스캔정보로 전환한다.The computing device converts the CT value (also referred to as a medical image or scan information) scanned by the computed tomography apparatus into a medical image or scan information including a skeleton parameter (eg, density) of a subject.
도 2에 도시된 바와 같이, 전환 시스템(300)은 컴퓨터 단층 촬영장치(310) 및 계산장치(330)를 포함한다. 계산장치(330)는 에너지 감쇠계수 계산모듈(331) 및 관계함수 구축모듈(333)을 포함하며, 관계함수 구축모듈(333)에 에너지 감쇠계수 계산모듈(331)이 연결된다.As shown in FIG. 2, the
컴퓨터 단층 촬영장치(310)에 계산장치(330)가 접속되며, 계산장치가 단계 140을 실행할 수 있도록, 접속은 유선 연결, 무선 연결이거나 또는 저장매체가 스캔정보를 계산장치(330)로 제공하는 것일 수 있다. 또한, 컴퓨터 단층 촬영장치(310)는 한 대 이상, 예를 들어 두 대, 세 대 또는 그 이상일 수 있다.The
에너지 감쇠계수 계산모듈(331)은 단계 150을 실행하며, 에너지 감쇠계수 계산모듈(331)은 스캔정보를 수신할 수 있고(단계 131), 또한 스캔정보에 의거하여 물리함수모델을 실행함으로써(단계 150), 컴퓨터 단층 촬영장치에 대응되는 에너지 감쇠계수 , 를 획득할 수 있다.The energy attenuation
관계함수 구축모듈(333)은 단계 170 및 단계 190을 실행하며, 단계 170의 에너지 보정계수 X는 관계함수 구축모듈(333) 내에 저장 또는 구축되고, 단계 190의 진실관계함수 역시 관계함수 구축모듈(333) 내에 구축된다.The relationship
이와 같이, 본 발명의 계산장치 및 그의 전환방법은 컴퓨터 단층 촬영장치를 통해 촬영한 스캔정보 및 기지의 팬텀 성분을 컴퓨터 단층 촬영장치의 에너지 감쇠계수로 전환하여 획득하고, 에너지 감쇠계수를 이용하여 피검체(골격)의 골격 파라미터를 획득할 수 있으며, 골격 파라미터는 밀도 및 분포를 포함한다. 이와 같이, 수술 계획 시 골격의 상태를 더욱 정확하게 획득할 수 있어 효율 및 의료 품질이 향상된다.As described above, the calculating apparatus and its conversion method of the present invention are obtained by converting the scan information and known phantom components photographed through the computed tomography apparatus into the energy attenuation coefficient of the computed tomography apparatus, and the energy attenuation coefficient is used. Skeletal parameters of the specimen (skeleton) can be obtained, and the skeleton parameters include density and distribution. In this way, it is possible to more accurately obtain the state of the skeleton when planning an operation, thereby improving efficiency and medical quality.
기타 실시예에서, 에너지 감쇠계수 계산모듈(31) 및 관계함수 구축모듈(33)은 각각 상이한 계산장치를 통해 실행될 수 있으며, 이는 에너지 감쇠계수 계산모듈은 계산장치 내에 구축되고, 관계함수 구축모듈은 즉 다른 계산장치에 구축됨을 나타낸다. 예를 들어 그 중 하나의 계산장치는 상기 단계 150을 실행하고, 다른 계산장치는 이미 획득한 컴퓨터 단층 촬영장치의 에너지 감쇠계수를 수신하여, 상기 단계 170 및 단계 190을 통해 피검체의 골격 파라미터를 획득할 수 있다. 따라서, 계산장치는 하나로 한정되지 않는다.In other embodiments, the energy attenuation coefficient calculation module 31 and the relationship function building module 33 may be executed through different calculation devices, respectively, which the energy attenuation factor calculation module is built in the calculation device, and the relationship function building module is In other words, it indicates that it is built on another computing device. For example, one of the calculation devices executes the
마지막으로, 재차 강조할 점은, 본 발명의 상기 실시예에 공개된 단계 순서와 구성소자는 단지 예를 들어 설명한 것일 뿐, 본 발명의 범위를 제한하기 위한 것이 아니며, 기타 단계 순서의 변경, 등가 소자의 교체 또는 변화 역시 본원의 특허 출원 범위에 포함되어야 한다.Finally, it should be emphasized again that the sequence of steps and components disclosed in the above embodiments of the present invention are only described as examples, and are not intended to limit the scope of the present invention, and other step sequence changes, equivalents. Replacement or change of the device should also be included in the scope of the patent application of the present application.
100: 방법
110-193: 단계
300: 전환시스템
310: 컴퓨터 단층 촬영장치
330: 계산장치
331: 에너지 감쇠계수 계산모듈
333: 관계함수 구축모듈100: Method 110-193: Steps
300: conversion system 310: computed tomography apparatus
330: calculation device 331: energy attenuation coefficient calculation module
333: relation function building module
Claims (6)
컴퓨터 단층 촬영장치를 제공하는 단계;
피검체 및 성분이 이미 알려진 두 팬텀을 제공하는 단계;
상기 컴퓨터 단층 촬영장치를 통해 상기 피검체 및 상기 두 팬텀을 촬영하여 대응되는 피검체 스캔 정보 및 두 팬텀의 스캔 정보를 획득하는 단계;
계산장치를 통해 상기 피검체의 스캔 정보 및 상기 두 팬텀의 스캔 정보를 수신하는 단계;
상기 계산장치가 상기 두 팬텀에 대응하며 상기 두 팬텀의 기지(旣知)의 성분 및 상기 두 팬텀의 스캔 정보를 포함하는 물리함수모델을 통해 상기 컴퓨터 단층 촬영장치의 에너지 감쇠계수를 계산하는 단계;
상기 두 팬텀의 기지의 성분의 에너지 감쇠계수 및 이상적인 피검체의 에너지 감쇠계수와 관련된 에너지 보정계수를 상기 계산장치로 제공하는 단계; 및
상기 계산장치가 상기 컴퓨터 단층 촬영장치의 에너지 감쇠계수, 상기 피검체의 스캔 정보 및 상기 에너지 보정 계수를 포함하는 진실관계함수를 통해 상기 피검체의 골격 파라미터를 획득하는 단계를 포함하는
컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법.In a method for converting scan information of a computed tomography apparatus into a skeleton parameter,
Providing a computed tomography apparatus;
Providing two phantoms in which the subject and components are known;
Capturing the subject and the two phantoms through the computed tomography apparatus to obtain corresponding scan information of the subject and scan information of the two phantoms;
Receiving scan information of the subject and scan information of the two phantoms through a calculation device;
Calculating, by the computing device, an energy attenuation coefficient of the computed tomography apparatus through a physical function model corresponding to the two phantoms and including a known component of the two phantoms and scan information of the two phantoms;
Providing an energy attenuation coefficient of the known components of the two phantoms and an energy correction coefficient related to the energy attenuation coefficient of an ideal subject to the calculation device; And
And obtaining, by the calculation device, a skeleton parameter of the subject through a truth relationship function including an energy attenuation coefficient of the computed tomography apparatus, scan information of the subject, and the energy correction coefficient.
A method of converting scan information of a computed tomography apparatus into skeleton parameters.
상기 두 팬텀의 기지의 성분은 인산수소칼륨과 물의 원소의 밀도, 중량 백분율, 원자량 및 원자번호를 포함하는, 컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법.The method of claim 1,
Known components of the two phantoms include the density, weight percentage, atomic weight, and atomic number of elements of potassium hydrogen phosphate and water, a method of converting scan information of a computed tomography apparatus into a skeleton parameter.
상기 에너지 감쇠계수는 광전흡수계수 및 레일리산란계수를 포함하는, 컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법.The method of claim 1,
The energy attenuation coefficient includes a photoelectric absorption coefficient and a Rayleigh scattering coefficient. A method of converting scan information of a computed tomography apparatus into a skeleton parameter.
상기 계산장치는 에너지 감쇠계수 계산모듈 및 관계함수 구축모듈을 포함하고, 상기 에너지 감쇠계수 계산모듈은 상기 물리함수모델을 포함하며, 상기 관계함수 구축모듈은 상기 에너지 보정계수 및 상기 진실관계함수를 포함하는, 컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법.The method of claim 1,
The calculation device includes an energy attenuation coefficient calculation module and a relationship function building module, the energy attenuation factor calculation module includes the physical function model, and the relationship function building module includes the energy correction factor and the truth relationship function. A method of converting scan information of a computed tomography apparatus into skeleton parameters.
상기 물리함수모델은
을 포함하며,
은 두 팬텀의 광감쇠값이고, 는 물의 광감쇠값이며, 은 아보가르도상수이고, 는 밀도이며, 는 중량백분율이고, 는 원자량이며, , , 는 원자번호 상관계수이고, 는 상기 컴퓨터 단층 촬영장치의 레일리산란계수이고, 는 상기 컴퓨터 단층 촬영장치의 광전흡수계수이며, 는 상기 컴퓨터 단층 촬영장치의 스캔정보인, 컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법.The method of claim 1,
The physical function model is
Including,
Is the light attenuation value of the two phantoms, Is the light attenuation value of water, Is the Avogardo constant, Is the density, Is the weight percentage, Is the atomic weight, , , Is the atomic number correlation coefficient, Is the Rayleigh scattering coefficient of the computed tomography apparatus, Is the photoelectric absorption coefficient of the computed tomography apparatus, Is a method of converting scan information of the computed tomography apparatus, which is the scan information of the computed tomography apparatus, into a skeleton parameter.
상기 진실관계함수는
를 포함하며,
는 상기 피검체의 스캔 정보이고, 는 상기 피검체의 광감쇠값이며, X는 상기 에너지 보정계수인, 컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법.The method of claim 1,
The truth relationship function is
Including,
Is the scan information of the subject, Is a light attenuation value of the subject, and X is the energy correction factor, a method of converting scan information of a computed tomography apparatus into a skeleton parameter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020190019767A KR102167859B1 (en) | 2019-02-20 | 2019-02-20 | Method for converting scan information of computed tomography scanner into bone parameters |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020190019767A KR102167859B1 (en) | 2019-02-20 | 2019-02-20 | Method for converting scan information of computed tomography scanner into bone parameters |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20200102028A true KR20200102028A (en) | 2020-08-31 |
KR102167859B1 KR102167859B1 (en) | 2020-10-21 |
Family
ID=72234375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020190019767A KR102167859B1 (en) | 2019-02-20 | 2019-02-20 | Method for converting scan information of computed tomography scanner into bone parameters |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR102167859B1 (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004147863A (en) * | 2002-10-30 | 2004-05-27 | Aloka Co Ltd | X-ray bone density measuring apparatus |
KR20040051938A (en) * | 2002-12-13 | 2004-06-19 | 한국전자통신연구원 | Method for calibrating bone mineral density index variation and storage medium for storing program of executing the same |
JP2011239830A (en) * | 2010-05-14 | 2011-12-01 | Gunma Univ | Multipurpose phantom and using method thereof |
KR20180059326A (en) * | 2016-11-25 | 2018-06-04 | 서울대학교산학협력단 | Quantitative system and method for alveolar bone dnesity using computed tomography |
JP2018514363A (en) * | 2015-05-07 | 2018-06-07 | アンバロール | Method for estimating the distribution of bone density in at least part of a human skeleton |
-
2019
- 2019-02-20 KR KR1020190019767A patent/KR102167859B1/en active IP Right Grant
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004147863A (en) * | 2002-10-30 | 2004-05-27 | Aloka Co Ltd | X-ray bone density measuring apparatus |
KR20040051938A (en) * | 2002-12-13 | 2004-06-19 | 한국전자통신연구원 | Method for calibrating bone mineral density index variation and storage medium for storing program of executing the same |
JP2011239830A (en) * | 2010-05-14 | 2011-12-01 | Gunma Univ | Multipurpose phantom and using method thereof |
JP2018514363A (en) * | 2015-05-07 | 2018-06-07 | アンバロール | Method for estimating the distribution of bone density in at least part of a human skeleton |
KR20180059326A (en) * | 2016-11-25 | 2018-06-04 | 서울대학교산학협력단 | Quantitative system and method for alveolar bone dnesity using computed tomography |
Also Published As
Publication number | Publication date |
---|---|
KR102167859B1 (en) | 2020-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8290232B2 (en) | System and method for quantitative imaging of chemical composition to decompose more than two materials | |
US6904118B2 (en) | Method and apparatus for generating a density map using dual-energy CT | |
Byng et al. | X-ray characterization of breast phantom materials | |
US7885373B2 (en) | System and method for quantitative imaging of chemical composition to decompose multiple materials | |
Koyama et al. | Radiation dose evaluation in tomosynthesis and C‐arm cone‐beam CT examinations with an anthropomorphic phantom | |
US6909771B2 (en) | Three component x-ray bone densitometry | |
Hashim et al. | Evaluation of organ doses and specific k effective dose of 64-slice CT thorax examination using an adult anthropomorphic phantom | |
Spampinato et al. | Dosimetric changes with computed tomography automatic tube-current modulation techniques | |
Gulam et al. | Bone mineral measurement of phalanges: comparison of radiographic absorptiometry and area dual X-ray absorptiometry | |
Homolka et al. | Optimization of the composition of phantom materials for computed tomography | |
US11517278B2 (en) | System and method for basis material decomposition with general physical constraint for multi-energy computed tomography | |
KR102167859B1 (en) | Method for converting scan information of computed tomography scanner into bone parameters | |
US10846855B2 (en) | Method for converting scan information of a computed tomography scanner into bone parameters | |
Dinten et al. | Dual-energy x-ray absorptiometry using 2D digital radiography detector: application to bone densitometry | |
TWI695387B (en) | Method for converting scanning information of computer tomography scanner into bone parameters | |
Hintenlang et al. | Physical phantoms for experimental radiation dosimetry | |
JP6751791B1 (en) | How to convert computed tomography scan information into skeletal parameters | |
CN111493914A (en) | Method for converting scanning information of computer tomography scanner into bone parameters | |
Frimaio et al. | X-ray spectrometry applied for determination of linear attenuation coefficient of tissue-equivalent materials | |
Tang et al. | X-ray projection simulation based on physical imaging model | |
US11864942B2 (en) | Methods for calibrating a spectral X-ray imaging system to perform material decomposition | |
DE102019110648A1 (en) | Method for converting scan information from computed tomography scanners into bone parameters | |
Lai et al. | A curve fitting approach using ann for converting ct number to linear attenuation coefficient for ct-based pet attenuation correction | |
Ali et al. | Low-cost chest paediatric phantom for dose optimisation: construction and validation | |
Brevitt et al. | Avoiding Defensive Medicine among Pandemic |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant |