KR20040051938A - 골밀도 지수 보정 방법 및 이 프로그램을 기록한 컴퓨터로읽을 수 있는 기록매체 - Google Patents
골밀도 지수 보정 방법 및 이 프로그램을 기록한 컴퓨터로읽을 수 있는 기록매체 Download PDFInfo
- Publication number
- KR20040051938A KR20040051938A KR1020020079637A KR20020079637A KR20040051938A KR 20040051938 A KR20040051938 A KR 20040051938A KR 1020020079637 A KR1020020079637 A KR 1020020079637A KR 20020079637 A KR20020079637 A KR 20020079637A KR 20040051938 A KR20040051938 A KR 20040051938A
- Authority
- KR
- South Korea
- Prior art keywords
- bone density
- phantom
- density index
- bmd
- gray level
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 210000000988 bone and bone Anatomy 0.000 title claims description 29
- 229910052500 inorganic mineral Inorganic materials 0.000 title claims description 6
- 239000011707 mineral Substances 0.000 title claims description 6
- 238000003860 storage Methods 0.000 title description 4
- 230000037182 bone density Effects 0.000 claims abstract description 66
- 230000008859 change Effects 0.000 claims description 33
- 238000003384 imaging method Methods 0.000 claims description 20
- 238000012937 correction Methods 0.000 claims description 18
- 210000004872 soft tissue Anatomy 0.000 claims description 17
- 229920000642 polymer Polymers 0.000 claims description 12
- PPBRXRYQALVLMV-UHFFFAOYSA-N Styrene Chemical compound C=CC1=CC=CC=C1 PPBRXRYQALVLMV-UHFFFAOYSA-N 0.000 claims description 8
- -1 polyethylene Polymers 0.000 claims description 7
- 239000004698 Polyethylene Substances 0.000 claims description 3
- 239000004743 Polypropylene Substances 0.000 claims description 3
- 229920000728 polyester Polymers 0.000 claims description 3
- 229920000573 polyethylene Polymers 0.000 claims description 3
- 229920001155 polypropylene Polymers 0.000 claims description 3
- 229920002635 polyurethane Polymers 0.000 claims description 3
- 239000004814 polyurethane Substances 0.000 claims description 3
- 239000004952 Polyamide Substances 0.000 claims description 2
- 229920000058 polyacrylate Polymers 0.000 claims description 2
- 229920002647 polyamide Polymers 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 9
- 238000007689 inspection Methods 0.000 abstract 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 24
- 229910052782 aluminium Inorganic materials 0.000 description 24
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 23
- 230000006870 function Effects 0.000 description 23
- 210000000707 wrist Anatomy 0.000 description 13
- 238000001739 density measurement Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 208000001132 Osteoporosis Diseases 0.000 description 6
- 238000010521 absorption reaction Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 238000002835 absorbance Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000009547 dual-energy X-ray absorptiometry Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000001054 cortical effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000006386 Bone Resorption Diseases 0.000 description 1
- 206010017076 Fracture Diseases 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 208000001164 Osteoporotic Fractures Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000024279 bone resorption Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 229940011871 estrogen Drugs 0.000 description 1
- 239000000262 estrogen Substances 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000028327 secretion Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 210000000623 ulna Anatomy 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
본 발명은 엑스선 영상을 이용하여 골밀도를 측정하고자 할 때, 골밀도 지수 보정 방법과 상기 방법을 실현시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체에 관한 것으로, 대상물과 함께 촬영한 서로 다른 2개 이상의 두께를 갖는 영역을 구비하는 팬텀 엑스선 영상을 획득하고, 이들 팬텀을 이용하여 대상물의 골밀도 지수를 보정하는 방법을 제공한다. 이를 통해서 엑스선 촬영조건의 변화에 기인한 골밀도 지수 변화를 보정할 수 있게 된다.
Description
본 발명은 엑스선 영상을 이용하여 골밀도를 측정하고자 할 때, 엑스선 촬영조건의 변화에 기인한 골밀도 지수 변화 보정 방법과 상기 방법을 실현시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체에 관한 것이다.
골다공증은 비정상적으로 골량이 감소된 상태로 동통, 척추 및 대퇴부 등의 골절, 골의 변형 등을 동반하는 병적인 상태이다. 일반적으로 젊은 사람은 골의 생성과 흡수가 균형을 이루고 있다. 그러나, 골 생성 기능이 저하되는 노인이나 에스트로젠의 분비가 저하되어 골 흡수 기능이 향상되는 폐경기 이후의 여성은 골다공증의 발생 빈도가 아주 높다. 그러나, 일단 감소된 골량을 효과적으로 회복시키는 것은 쉽지 않다. 따라서, 골다공증의 발생방지와 조기치료가 무엇보다 중요하다. 이를 위해서는 일상적으로 널리 사용할 수 있는 골밀도 측정 방법의 개발이 매우 중요하다.
정량적인 골밀도 측정 방법으로는 QCT(Quantitative Computed Tomography),DEXA(Dual Energy X-ray Absorptiometry), 초음파 골밀도 계측기 등이 있다.
QCT는 뼈의 3차원 영상을 얻을 수 있어서, 피질골(cortical bone)과 해면골(trabecular bone)의 골밀도를 분리하여 계측할 수 있을 뿐만 아니라, 골의 구조적 안정성까지 예측할 수 있는 장점이 있다. 이러한 장점에도 불구하고 QCT 장비의 가격이 아주 높고 방사선 조사량이 단순 엑스선 촬영의 수백 배에 이르기 때문에 일상적인 골다공증 검사 방법으로 사용하기에는 많은 제약이 따른다.
DEXA는 QCT에 비하여 장비 가격이 저렴하고, 방사선 조사량이 적고, 골밀도 계측의 정확도 및 재현도가 높아서 2차원 골밀도 측정 및 골다공증 치료의 경과판정에 널리 사용되고 있다.
초음파 골밀도 계측기는 매질에 따른 초음파 전달속도의 변화나 반사계수의 변화 등을 이용하여 골밀도 및 골의 구조강도를 계측하는 장비인데, 재현성과 정확성이 높지 않은 단점이 있다.
상기의 QCT, DEXA, 초음파 골밀도 계측기 등의 골밀도 측정 방법이 있음에도 불구하고 엑스선 영상을 이용하여 골다공증 진단 및 골밀도를 측정하는 방법의 연구는 꾸준히 진행되고 있다. 그 이유는 대부분의 병원이 엑스선 촬영장비를 기본으로 갖추고 있어 골밀도 계측장비를 구매하기 위한 추가 비용이 소요되지 않으며, 엑스선 영상의 해상도가 높아 골밀도 측정 및 해면골의 패턴분석이 가능하기 때문이다. 특히, 해면골 패턴은 골다공증에 의한 골절 가능성에 관한 유용한 정보를 제공하는 것으로 추정되고 있어, 많은 연구자들이 엑스선 해면골 패턴을 연구하고 있다.
현재까지, 엑스선 영상을 이용하여 골밀도를 계측하기 위한 몇 가지 방법이 제시되었는데, 임상에서 널리 사용되지 않는 실정이다. 이것은 엑스선 영상의 획득과정에 수반하는 각종 불안정성(엑스선 촬영조건 및 필름 현상과정의 불안정성)이 상당한 골밀도 오차를 일으킬 수 있기 때문이다. 특히, 엑스선 촬영조건의 불안정성은 필름 현상과정의 불안정성과는 달리 사용자의 주의만으로는 제거하기 어렵다.
따라서, 본 발명은 상기 문제점을 해결하기 위하여 제안된 것으로, 엑스선 영상을 이용하여 골밀도를 측정하고자 할 때, 엑스선 촬영조건의 변화에 기인한 골밀도 지수 변화를 보상하기 위하여 팬텀을 이용하는 방법과 상기 방법을 실현시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체를 제공하는데 그 목적이 있다.
도 1은 본 발명의 바람직한 실시예에 따른 골밀도 지수 변화 보상 방법을 실현시키기 위한 하드웨어 시스템을 도시한 구성 블록도이다.
도 2는 아크릴 팬텀을 포함하는 엑스선 영상을 이용한 골밀도 지수 보정 방법에 대한 일실시예 흐름도이다.
도 3은 아크릴 팬텀 구조의 일예를 도시한 도면이다.
도 4는 아크릴 팬텀을 포함하는 엑스선 영상의 일예를 도시한 도면이다.
도 5는 알루미늄 팬텀 구조의 일예를 도시한 도면이다.
도 6은 아크릴 팬텀 및 알루미늄 팬텀을 포함하는 엑스선 영상의 일예를 도시한 도면이다.
도 7은 요골 골밀도 지수 측정 방법에 대한 일실시예 흐름도이다.
도 8은 알루미늄 팬텀의 두께와 평균 그레이 레벨의 상관관계를 도시한 그래프이다.
도 9는 알루미늄 팬텀을 사용하여 보정한 손목 엑스선 영상 및 요골 골밀도 지수를 측정하기 위하여 선택한 직사각형 관심영역의 일예시도이다.
도 10은 도 9의 한 횡단선l에서의 픽셀좌표에 따른 그레이 레벨 프로파일을 도시한 그래프이다.
* 도면의 주요 부분에 대한 부호의 설명
11 : 입출력 장치 12 : 주 기억장치
13 : 보조 기억장치 14 : 마이크로 프로세서
상술한 문제점을 해결하기 위한 수단으로서, 본 발명은 골밀도 지수를 보정하는 방법에 있어서, (a) 대상물과 2개 이상의 서로 다른 두께를 갖는 영역을 구비하는 팬텀을 동시에 촬영한 영상을 얻는 단계와, (b) 촬영된 영상으로부터 대상물의 골밀도 지수와 팬텀의 각 영역에서의 평균 그레이 레벨을 계산하는 단계와, (c) (a) 및 (b) 단계를 N회 반복 측정하여, 골밀도 지수와 팬텀 각각 영역에서의 평균 그레이 레벨 사이의 상관 관계식을 추출하는 단계와, (d) 상관 관계식을 통해서, 대상물의 엑스선 촬영시의 조건 변화에 기인한 요골 골밀도 지수 변화를 보정하는 방법을 제공한다.
한편, 팬텀은 M(2 이상)개의 서로 다른 두께를 갖는 영역을 가지며, 상기 상관관계식은 H(G,0,0,0,...)=0을 만족하는 연속함수를 이용하여 다음 식에 의해 표시된다.
η=H(G, A-A0,B-B0,C-C0, ...), (여기서, η는 대상물의 골밀도 지수 변화량, G는 대상물의 골밀도 지수, A, B, C...는 팬텀 영역의 서로 다른 M개의 두께 영역에서의 평균 그레이 레벨, A0,B0,C0...는 엑스선 촬영을 N회 실시한 후 팬텀 영역의 서로 다른 M개의 두께 영역 각각의 영상에서 획득한 A, B, C.. 를 평균한 값).
"대상물"이라 함은 골밀도를 측정하고자 하는 뼈의 부위를 의미하는 것으로 인체 또는 동물의 모든 뼈의 부위를 포함가능하다.
"팬텀"은 피부조직과 유사한 물질을 이용하는 것이 바람직하며, 예를 들어 아크릴계 폴리머(acrylic polymer), 스타이렌계 폴리머(styrene polymer), 폴리에틸렌(polyethylene), 폴리프로필렌(polypropylene), 폴리에스터계 폴리머(polyester polymer), 폴리아미드계 폴리머(polyamide polymer), 폴리우레탄계 폴리머(polyurethane polymer) 등이 가능하다.
한편, 상술한 문제점을 해결하기 위한 수단으로서 본 발명의 다른 측면은 상기 골밀도 지수 보정 방법을 실현시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체를 제공한다.
이하, 첨부된 도면들을 참조하면서 본 발명의 바람직한 실시예에 따른 골밀도 지수 변화 보상 방법에 대해 상세히 설명한다. 본 실시예에서는 2개의 서로다른 두께 영역을 갖는 아크릴 팬텀을 이용한 경우를 설명하고 있으나, 본 발명이 이에 한정되는 것이 아님은 자명하다.
도 1 은 본 발명의 바람직한 실시예에 따른 골밀도 지수 변화 보상 방법을 실현시키기 위한 하드웨어 시스템을 도시한 구성 블록도이다.
하드웨어 시스템은 외부의 사용자와 데이터를 입/출력하기 위한 입출력 장치(11), 엑스선 영상을 이용하여 골밀도를 지수를 측정하는 과정에서 필요한 각종 데이터들을 저장하기 위한 주/보조 기억장치(12,13) 및 상기 주/보조 기억장치(12,13)와 입/출력 장치(11)를 제어하고, 엑스선 영상을 이용한 골밀도 지수를 측정하며, 아크릴 팬텀을 이용한 골밀도 지수 변화 보상 방법을 수행하기 위한 제반 연산처리를 수행하는 마이크로 프로세서(14)를 구비한다.
상기 입/출력 장치(11)는 모니터, 프린터, 엑스선 필름을 디지털 영상화하기 위한 엑스선 필름스캐너 등이 가능하고, 엑스선 필름 대신에 디지털 이미지 센서[고체촬상소자(CCD) 또는 씨모스(CMOS) 센서]를 사용할 경우에는 디지털 이미지 센서가 상기 입/출력 장치(11)에 포함될 수 있다.
주/보조 기억장치(12,13)은 본 발명의 골밀도 지수 변화 보정방법을 구현하는 프로그램을 컴퓨터로 읽을 수 있는 형태로 각종 기록매체인 씨디롬, 램, 롬, 플로피 디스크, 하드 디스크, 광자기 디스크 등이 가능하다.
상기와 같은 시스템을 통하여, 골밀도 지수 변화 보상 방법을 실행시키기 위한 프로그램을 실행하게 되는데, 이와 같은 처리 과정을 포함하고 있는 프로그램을 상기 마이크로 프로세서(14)에 내장시킨 상태에서 엑스선 영상을 입출력 장치(11)로 입력하여 상기 프로그램을 실행시키면, 프로그램이 골밀도 지수를 측정하고, 골밀도 지수 변화 보상 방법을 수행한다.
이하, 도 2 내지 도 4를 참조하여 아크릴 팬텀을 포함하는 엑스선 영상을 이용한 골밀도 지수 보정 방법에 대한 일실시예를 설명한다.
먼저, 도 3에 도시한 아크릴 팬텀과 대상물(본 실시예에서는, 요골)을 함께 촬영한 엑스선 영상을 획득한다(S101). 도 3은 본 발명에 따른 아크릴 팬텀을 도시한 도면이고 본 실시예에서는 2단 아크릴 팬텀의 일예를 도시하고 있다. 아크릴의 밑면은 60mm x 30mm 이고, 높이는 각각 60mm와 40mm 이다. 도 4는 아크릴 팬텀과 손목이 함께 촬영된 엑스선 영상의 일예를 도시하고 있다.
엑스선 영상을 얻을 때, 엑스선 발생장치의 관전압(kVp)은 일정한 값을 유지한다. 예를 들어 엑스선 발생장치의 관전압은 50kVp(이하, 표준 관전압)를 유지한다. 엑스선 영상은 단순 엑스선 촬영을 통하여 얻어진 엑스선 필름을 엑스선 필름스캐너로 디지털화 함으로써 획득할 수 있다. 디지털 이미지 센서를 사용할 경우에는 필름스캐너를 이용하는 과정 없이 직접 엑스선 영상을 얻는다. 엑스선 영상을 얻을 때, 공간 해상도는 예컨대 200PPI(Pixels Per Inch), 각 픽셀은 8-비트 깊이의 256색 그레이레벨로 할 수 있다.
이후, 획득한 손목 엑스선 영상의 요골 골밀도 지수(G)와 아크릴의 평균 그레이레벨을 계산한다(S103). 아크릴의 평균 그레이레벨은 각각 40mm 및 60mm 두께 영역에서 각각 측정한 값이다. 이 경우, 손목 엑스선 영상의 요골 골밀도 지수(G)의 측정은 통상적인 방법으로 가능하고, 손목 엑스선 영상에 본 발명의 보정방법과 무관한 다른 보정 방법을 1차로 실시하고 난 다음, 본 발명을 실시하는 것도 가능하다. 예를 들어, 동출원인에 의한 한국특허출원제2001-45123호에서는 손목 엑스선 영상과 알루미늄 팬텀 영상을 함께 촬영하여 엑스선 영상을 보정하는 방법을 상세하게 개시하고 있는바, 본 발명에 의한 아크릴 팬텀을 이용한 보정 방법 이전에 미리 알루미늄 팬텀을 이용한 보정을 별개로 실시하는 것도 가능하다. 이에 대해서는 상세히 후술한다.
이후, 요골 골밀도 지수 변화량을 계산한다(S105). 이 방식을 상세히 설명하면, 60mm 아크릴의 보정 후 평균 그레이레벨을 A, 40mm 아크릴의 보정 후 평균 그레이레벨을 B라고 가정하자. 또한, G0 ,A0및 B0는 각각 정확한 표준 관전압(50kVp)에서의 G, A 및 B라고 가정하자. 만약, 엑스선 발생 조건이 정확하다면, G, A 및 B는 각각 G0 ,A0및 B0와 동일하다. 그러나, 실제 임상에서 사용하는 엑스선 촬영장비는 엑스선 발생 조건이 변화할 수 있어서 G, A 및 B는 각각 G0 ,A0및 B0와 일치하지 않는다. 그런데, 엑스선 촬영실험을 하면, 골밀도 지수 변화량(η=G-G0), G, A-A0및 B-B0사이에는 근사적 함수관계 있음을 알 수 있다. 따라서, η는 일반적으로 다음의 식 (1)과 같이 표현할 수 있다.
η=H(G, A-A0, B-B0)) (1)
H는 H(G, 0, 0)을 만족하는 연속함수이다. 다만, 본 실시예에서는 아크릴 팬텀이 2개의 서로 다른 두께를 갖는 영역이 포함하는 경우를 예를 들어 설명하고 있으나, M개로 그 범위를 확장하는 것도 가능하다. 이 경우, 상기 식(1)는 다음식 (2)로 표현할 수 있다.
η=H(G, A-A0, B-B0, C-C0,...) (2)
한편, 엑스선 실험에 의하면, 식(2)의 서로 다른 2개의 두께 영역을 갖는 아크릴인 경우 B-B0가 증가하면 η은 감소하고, A-A0가 증가하면 η은 증가하는 현상이 있었으며, η는 G에 비례하는 경향이 발견되었다. 이러한 현상을 통해서 실험식을 만들어 보면, 식(3)와 같은 관계식이 가능하다.
η=C1G(C2(A-A0)-(B-B0)) (3)
여기에서 c1, c2는 상수이다. 상기 식(3)이 엑스선 촬영조건 변화에 기인한 요골 골밀도 지수 변화량이다.
이후, 요골 골밀도 지수(G)에서 상기 요골 골밀도 지수 변화량(η)을 차감하여 요골 골밀도 지수 변화를 보상한다(S107). 요골 골밀도 지수의 변화량을 아크릴의 그레이 레벨을 이용하여 식(3)으로 표현했으므로, 표준 관전압에서의 요골 골밀도 지수는 다음의 식(4)와 같이 표현된다.
G0=G-C1G(c2(A-A0)-(B-B0)) (4)
상기 식(4)가 엑스선 촬영조건 변화에 기인한 요골 골밀도 지수 변화를 40mm 및 60mm 아크릴을 이용하여 보상한 결과이다. 한편, 수학식 4를 실제 임상에서 사용하기 위해서는 상수 c1, c2, A0및 B0를 결정해야 한다. 그러나, A0및 B0는 표준 관전압의 안정성이 완벽하지 못할 수 있으므로 다수의 측정결과를 평균하여 설정한다. 예를 들면, 표준 관전압(50kVp)으로 엑스선 촬영을 10회 실시한 후, 각각의 영상에서 획득한 A 및 B를 평균하여 A0및 B0로 설정한다. 한편, 상수 c1, c2도 다수의 피험자에 대한 엑스선 촬영을 통하여 결정한다. 각 피험자에 대하여 10회의 엑스선 촬영을 실시하고 c1, c2를 변화시키면서 동일 피험자 내에서 η의 표준편차(ηS)를 계산한다. 이후, 상수 c1, c2는 전체 피험자의 ηS의 합이 최소가 되는 값으로 설정한다. 피험자 10인을 대상으로 한 손목 엑스선 촬영 실험에서 설정된 최적 c1는 0.005, c2는 0.72이다. 한편, 최적 c2는 엑스선 장비에 따라서 0.6~0.8사이의 변화를 보인다.
(손목 엑스선 영상의 요골 골밀도 지수(G)의 측정)
이하, 전술한 바와 같이, 한국특허출원 제2001-45123호에 개시되어 있는 손목 엑스선 영상과 알루미늄 팬텀 영상을 함께 촬영하여 엑스선 영상을 보정하는 방법과 본 발명의 아크릴 팬텀을 이용한 골밀도 지수 보정 방법을 함께 실시하는 방법에 대해서 설명한다. 이 경우는 알루미늄 팬텀과 아크릴 팬텀을 손목 영상과 함께 촬영한다(도 6 참조). 알루미늄 팬텀은 엑스선 영상을 보정(calibration)하기 위하여 사용되는 것으로, 엑스선 발생장치의 특성, 스크린 및 필름 종류, 필름 현상조건, 및 엑스선 필름 디지타이저 특성(디지털 이미지 센서의 경우에는 고체촬상소자 또는 씨모스 센서의 특성) 등에 의한 영상특성 변화를 보정할 수 있다. 도 5에서는 이와 같은 알루미늄 팬텀의 구조의 일예를 도시하고 있다. 도 7은 알루미늄 팬텀과 함께 손목을 촬영한 엑스선 영상을 이용하여 요골 골밀도 지수(G)의 측정 과정을 도시한 흐름도이다.
먼저, 촬영한 영상으로부터 알루미늄 팬텀과 손목의 영상정보를 얻는다(S301).
손목 엑스선 영상을 보정하여 각 픽셀에서의 엑스선 흡수도를 정량화하기 위하여 도 5 에 도시한 바와 같이 다수의 단계로 이루어진 알루미늄 팬텀을 사용한다. 도 5 에 예시한 알루미늄 팬텀은 밑면이 40㎜ x 40㎜인 정사각형이고 두께가 12㎜인 알루미늄 판을 8개의 단계를 갖는 방사상 계단으로 가공하여 만든다. 상기 알루미늄 팬텀의 각 단계의 높이는 낮은 것으로부터 1.5, 3.0, 4.5, 6.0, 7.5, 9.0, 10.5, 12.0㎜이다.
엑스선을 이용한 골밀도 측정은 일반적으로 뼈에 의한 엑스선 흡수도를 이용한다. 따라서, 알루미늄 팬텀의 각 단계별 두께와 평균 그레이레벨 사이의 관계를 이용하여 전체 엑스선 영상을 보정한다(S303). 이렇게 하면, 보정한 엑스선 영상은 알루미늄의 엑스선 흡수도에 상대적인 수치로 전환되어 인체에 의한 엑스선 흡수도를 정량화한 것이 된다.
알루미늄 팬텀을 이용하여 전체 엑스선 영상을 보정하기 위해서, 먼저 임의의 두께에서의 그레이레벨 값을 계산해야 한다. 임의의 알루미늄 두께에서의 그레이레벨을 계산하기 위하여, 도 8 의 알루미늄 팬텀의 두께와 평균 그레이 레벨의 상관관계를 도시한 데이터를 적절한 함수를 이용하여 피팅(fitting)할 수 있다. 도 8을 참고하면, 알루미늄 두께가 커지면서 처음에는 그레이레벨의 기울기가 점차 증가하다가 나중에는 다시 감소한다. 이러한 특징을 보이는 대표적인 함수가 탄젠트 하이퍼볼릭(tangent hyperbolic) 함수이다. 따라서, f(t) = a + b*tanh(c*t+d) 형태의 피팅함수를 사용하는 것이 바람직하다. 여기서, f(t)는 그레이레벨, t는 알루미늄의 두께(㎜단위), a, b, c, d는 각각 피팅 파라미터이다. 여기서, tanh 함수는 t=-d/c인 점을 중심으로 완전한 대칭이다. 그러나, 도 8 에 도시한 것처럼, 알루미늄 팬텀의 그레이레벨 프로파일은 완전한 대칭이 아니다. 따라서, 하나의 tanh 함수를 이용하여 피팅하는 것은 상당한 피팅오차를 야기할 수 있기 때문에, 데이터를 두 개의 영역으로 나누어 피팅을 실시한다. 상기의 피팅함수는 4개의 피팅 파라미터가 있어서 한 피팅 영역에 4개 이상의 데이터를 필요로 한다.
첫번째 영역은 높이가 0㎜ 단계로부터 7.5㎜인 단계까지 6개의 데이터로 구성하고, 두번째 영역은 높이가 4.5㎜인 단계로부터 12.0㎜인 단계까지 6개의 데이터로 구성한다. 여기서, 두 개의 영역을 서로 중복되게 한 것은 두 영역에서의 피팅 결과를 부드럽게 연결하기 위한 이행영역을 만들기 위함이다. 이 후, 각 영역에서 상기 피팅함수로 피팅을 실시한다. 피팅은 레벤버그-마르쿼트(Levenberg-Marquardt) 피팅 방법을 이용할 수 있다.
첫번째 영역에서의 피팅결과를 f1(t), 두 번째 영역에서의 피팅결과를 f2(t)라 하고, 그 결과를 도 8에 각각 실선과 점선으로 표시하였다. 이제 두 개의 피팅함수를 다음의 방법에 의하여 하나의 최종 피팅 함수 F(t)로 만든다. 먼저, t <= 4.5 구간에서는 F(t) = f1(t)로 하고, t >= 7.5 구간에서는 F(t) = f2(t)로 한다. 4.5 < t < 7.5 구간에서는 F(t) = x*f1(t) + (1-x)f2(t)로 한다. 여기서, x = (7.5-t)/3으로 정의한다. 함수 F(t)는 단순 증가함수이기 때문에 역함수 F-1(g)는 일의적으로 결정된다. 여기서, g는 그레이 레벨이다. 이제, 엑스선 영상의 각 픽셀의 그레이 레벨 g를 최종 피팅함수 F(t)를 이용하여 다음과 같이 보정한다. g>=F(12)이면 보정값을 255로, g<= F(0)이면 보정값을 0으로, 그 외에는 보정값을 256*F-1(g)/12의 정수부분으로 한다.
도 9는 상술한 방식에 의하여 보정한 영상의 일예시도이다. 도 9의 각 픽셀은 인체에 의한 엑스선 흡수도와 직결된 그레이 레벨로 이루어져 있다. 인체에 의한 엑스선 흡수도는 뼈와 연부조직의 밀도 및 두께 등에 의하여 결정된다. 그러나, 엑스선 골영상은 중첩된 연부조직에 의한 엑스선 흡수도를 동시에 포함하고 있다. 따라서, 엑스선 골영상의 그레이레벨을 이용하여 뼈만의 밀도를 측정하기 위해서는 엑스선 골영상에 포함된 연부조직에 의한 엑스선 흡수도를 제거하여야 한다.
이후, 엑스선 골영상에 포함된 연부조직에 의한 엑스선 흡수도를 제거하기 위하여 요골 영역에 관심영역을 설정한다(S305). 도 9에 도시한 관심영역의 크기는예컨대 350x300 픽셀이고, 요골 영역 좌우에 연부조직 영역이 포함되도록 한다.
이후, 요골 영역에서 연부조직에 의한 백그라운드 트렌드를 계산하기 위하여 피팅함수를 선택한다(S307). 도 10에 도 9의 관심영역내 하나의 횡단선 l에서의 그레이레벨 프로파일을 실선으로 도시하였다. 도 10에서 가로는 픽셀의 좌표이고 세로는 그레이레벨이다. 도 10에서 구간 al~bl,bl~cl,cl~dl은 각각 연부조직 구간, 요골 구간, 요골과 척골(ulna) 사이의 연부조직 구간이다. 도 10의 요골 구간에서, 연부조직에 의한 백그라운드 트렌드를 정확하게 계산하는 것은 불가능하다. 따라서, 백그라운드 트렌드를 계산하기 위한 근사적 방법을 사용한다. 여기서, 사용하는 방법은 연부조직 al~bl, cl~dl구간과의 그레이레벨 프로파일을 요골 구간 bl~cl으로 내삽(interpolation)하여 백그라운드 트렌드로 설정하는 것이다. 내삽을 위해서는 먼저 미분가능한 피팅 함수를 선택한다. 일반적으로, 피팅 함수로서 다항식이 적당하지만, 여기서는 좀더 제한을 가하여 4차 이하의 다항식[P(x) = C0+ C1x + C2x2+ C3x3+ C4x4]을 피팅함수로 사용한다. 여기서 C0, C1, C2, C3, C4는 피팅 파라미터이다.
이후, 연부조직 영역의 그레이레벨 프로 파일을 상기 피팅함수를 사용하여 요골 영역으로 내삽하여 백그라운드 트렌드를 계산한다(S309). 이 과정을 좀 더 상세히 설명하면, 요골에 인접한 연부조직 구간 al~bl, cl~dl에 있는 그레이레벨 프로파일을 상기 피팅함수를 사용하여 요골 영역으로 내삽한다. 한편, 내삽은 레벤버그-마르쿼트(Levenberg-Marquardt) 피팅 방법으로 실시한다. 요골 영역의 내삽결과가 연부조직에 의한 백그라운 트렌드이다. 도 10 의 Pl은 상기 피팅함수로 내삽하여 설정한 백그라운드 트렌드이다. 이제, 상기의 횡단선l을 관심영역 안의 모든 줄(row)로 이동하면서 연부조직에 의한 백그라운드 트렌드를 계산한다.
이와 같이 연부조직에 의한 백그라운드 트렌드가 설정되었으면, 요골 영역의 그레이레벨에서 백그라운드 트렌드를 제거한 후 요골 골밀도 지수(G)를 계산한다(S311). 이 과정을 좀 더 상세히 살펴보면, 상기 골 영역의 그레이레벨에서 상기 설정된 백그라운드 트렌드를 제거하는 과정 및 상기 요골 영역에서 상기 백그라운드 트렌드를 제거한 그레이레벨의 평균을 요골 골밀도 지수로 설정하는 과정을 포함한다. 따라서, 요골 골밀도 지수는 다음의 식(5)과 같이 표현된다.
(5)
상기 수학식 5에서 Gln은 연부조직에 의한 백그라운 트렌드가 제거된 그레이레벨 프로파일이고, n은 픽셀의 인덱스이고, A()는 골 영역의 면적이다. 여기서, 골밀도 지수 G가 바로 골밀도의 척도이다.
이와 같은 방식으로 측정된 골밀도 지수(G)는 도 2의 S103에 골밀도 지수로 이용될 수 있다. 다음으로 전술한 S103단계 내지 S107단계를 수행하여 본 발명에 의한 골밀도 지수 보정 방법을 수행한다.
상술한 바와 같이, 본 발명의 구체적인 예를 보였는데, 여기서 사용하는 수치들이나 영상들은 본 발명에 의한 방법의 성능 향상을 위하여 달라질 수 있다. 즉, 기존의 엑스선 영상을 이용한 골밀도 측정에서는 엑스선 촬영조건 등의 변화에 기인한 골밀도 변화를 보상하지 않았기 때문에 엑스선 촬영장비의 불안정성에 의한 측정오차가 발생하였다. 본 발명에서는 엑스선 촬영조건의 변화에 기인한 골밀도 변화를 보상하기 위하여 40mm 및 60mm 아크릴을 이용하였다.
이상에서 설명한 본 발명은 전술한 실시예 및 첨부된 도면에 의해 한정되는 것이 아니고, 본 발명의 기술적 사상을 벗어나지 않는 범위 내에서 여러 가지 치환, 변형 및 변경이 가능하다는 것이 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 있어 명백할 것이다.
상기와 같은 본 발명은, 엑스선 영상을 이용하여 골밀도를 측정하고자 할 때, 엑스선 영상의 획득과정에 수반하는 각종 불안정성 (엑스선 촬영조건 및 필름 현상과정의 불안정성)등에 기인한 골밀도 지수 변화를 보정하여 골밀도 측정을 더욱 정확하게 하는 효과가 있다.
Claims (7)
- 골밀도 지수를 보정하는 방법에 있어서,(a) 대상물과 2개 이상의 서로 다른 두께를 갖는 영역을 구비하는 팬텀을 동시에 촬영한 영상을 얻는 단계;(b) 촬영된 영상으로부터 상기 대상물의 골밀도 지수와 상기 팬텀의 각 영역에서의 평균 그레이 레벨을 계산하는 단계;(c) 상기 (a) 및 (b) 단계를 N회 반복 측정하여, 골밀도 지수와 상기 팬텀 각각 영역에서의 평균 그레이 레벨 사이의 상관 관계식을 추출하는 단계; 및(d) 상기 상관 관계식을 통해서, 상기 대상물의 엑스선 촬영시의 조건 변화에 기인한 요골 골밀도 지수 변화를 보정하는 것을 특징으로 하는 골밀도 지수 보정 방법.
- 제 1 항에 있어서,상기 팬텀은 M(2 이상)개의 서로 다른 두께를 갖는 영역을 가지며, 상기 상관관계식은 H(G,0,0)=0을 만족하는 연속함수를 이용하여 식(6)에 의해 표시되는 것을 특징으로 하는 골밀도 지수 보정 방법.η=H(G, A-A0, B-B0, C-C0, ...)) (6)(여기서, η는 대상물의 골밀도 지수 변화량, G는 대상물의 골밀도 지수, A,B, C...는 팬텀 영역의 서로 다른 M개의 두께 영역에서의 평균 그레이 레벨, A0,B0,C0...는 엑스선 촬영을 N회 실시한 후 팬텀 영역의 서로 다른 M개의 두께 영역 각각의 영상에서 획득한 A, B, C.. 를 평균한 값)
- 제 2 항에 있어서,상기 M은 2이며, 상기 연속함수는 식(7)으로 표시된 함수인 것을 특징으로 하는 골밀도 지수 보정 방법.η=C1G(C2(A-A0)-(B-B0)) (여기서, c1, c2는 상수) (7)
- 제 1 항에 있어서,상기 팬텀은 아크릴계 폴리머(acrylic polymer), 스타이렌계 폴리머(styrene polymer), 폴리에틸렌(polyethylene), 폴리프로필렌(polypropylene), 폴리에스터계 폴리머(polyester polymer), 폴리아미드계 폴리머(polyamide polymer) 또는 폴리우레탄계 폴리머(polyurethane polymer)로 만들어 진 것을 특징으로 하는 골밀도 지수 보정 방법.
- 제 1 항에 있어서,상기 (b) 단계의 골밀도 지수는,상기 대상물 인근 연부조직에 의한 그레이레벨을 상기 대상물 영역으로 내삽하여 백그라운드 트렌드로 설정하는 단계, 상기 대상물 영역의 그레이레벨에서 상기 백그라운드 트렌드를 제거하는 단계, 상기 대상물 영역에서 상기 백그라운드 트렌드가 제거된 그레이레벨의 평균을 골밀도 지수로 설정하는 단계를 포함하는 것을 특징으로 하는 골밀도 지수 보정 방법.
- 상기 제 5 항에 있어서,상기 골밀도 지수는 식(8)로 표현되는 것을 특징으로 하는 골밀도 지수 보정 방법.(8)(여기서, Gln은 백그라운 트렌드가 제거된 그레이레벨 프로파일이고, n은 픽셀의 인덱스이고, A()는 대상물 영역의 면적, G는 골밀도 지수)
- 제 1 내지 제 6 항 중 어느 하나의 항에 의한 골밀도 지수 보정 방법을 실현시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020020079637A KR20040051938A (ko) | 2002-12-13 | 2002-12-13 | 골밀도 지수 보정 방법 및 이 프로그램을 기록한 컴퓨터로읽을 수 있는 기록매체 |
US10/672,726 US20040114726A1 (en) | 2002-12-13 | 2003-09-25 | Method for Calibrating bone mineral density index variation and recording medium for storing program for executing the same |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020020079637A KR20040051938A (ko) | 2002-12-13 | 2002-12-13 | 골밀도 지수 보정 방법 및 이 프로그램을 기록한 컴퓨터로읽을 수 있는 기록매체 |
Publications (1)
Publication Number | Publication Date |
---|---|
KR20040051938A true KR20040051938A (ko) | 2004-06-19 |
Family
ID=32501406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020020079637A KR20040051938A (ko) | 2002-12-13 | 2002-12-13 | 골밀도 지수 보정 방법 및 이 프로그램을 기록한 컴퓨터로읽을 수 있는 기록매체 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20040114726A1 (ko) |
KR (1) | KR20040051938A (ko) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20200102028A (ko) * | 2019-02-20 | 2020-08-31 | 하이윈 테크놀로지스 코포레이션 | 컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7471761B2 (en) * | 2005-09-15 | 2008-12-30 | Schick Technologies, Inc. | System and method for computing oral bone mineral density with a panoramic x-ray system |
US7734080B2 (en) * | 2005-12-08 | 2010-06-08 | Electronics And Telecommunications Research Institute | Method and device for measuring pulmonary nodule volume using dual-energy substraction image |
ES2396772B1 (es) * | 2010-05-20 | 2014-01-27 | Universitat De Barcelona | Método y sistema de recopilación de datos para su uso en la detección y/o diagnóstico de las lesiones óseas. |
GB201107385D0 (en) * | 2011-05-04 | 2011-06-15 | Materialise Nv | Medical imaging calibration device |
CN103892856B (zh) * | 2012-12-31 | 2017-10-20 | 深圳先进技术研究院 | 一种获取人体骨密度值的方法及系统 |
EP2964093A4 (en) * | 2013-03-06 | 2016-12-07 | Marika Pty Ltd | EVALUATION OF OPTICAL DENSITY GRADIENTS AND VARIATIONS |
FR3004099B1 (fr) * | 2013-04-03 | 2015-04-03 | Designers Developers Distributors Associates D3A Medical Systems | Procede et systeme de caracterisation d'un tissu osseux. |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0542130A (ja) * | 1991-08-09 | 1993-02-23 | Matsushita Electric Ind Co Ltd | 校正用フアントムおよびデータ校正法 |
JPH07148143A (ja) * | 1993-11-26 | 1995-06-13 | Teijin Ltd | 骨計測方法及び装置 |
JP2001128967A (ja) * | 1999-11-01 | 2001-05-15 | Aloka Co Ltd | 超音波診断装置用ファントム、ファントムを用いた誤差修正システム及び誤差修正方法 |
KR20010055220A (ko) * | 1999-12-10 | 2001-07-04 | 오길록 | 톱니 모양의 래크를 이용한 골소주 지표 보정방법 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5335260A (en) * | 1992-11-25 | 1994-08-02 | Arnold Ben A | Calibration phantom and improved method of quantifying calcium and bone density using same |
-
2002
- 2002-12-13 KR KR1020020079637A patent/KR20040051938A/ko not_active Application Discontinuation
-
2003
- 2003-09-25 US US10/672,726 patent/US20040114726A1/en not_active Abandoned
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0542130A (ja) * | 1991-08-09 | 1993-02-23 | Matsushita Electric Ind Co Ltd | 校正用フアントムおよびデータ校正法 |
JPH07148143A (ja) * | 1993-11-26 | 1995-06-13 | Teijin Ltd | 骨計測方法及び装置 |
JP2001128967A (ja) * | 1999-11-01 | 2001-05-15 | Aloka Co Ltd | 超音波診断装置用ファントム、ファントムを用いた誤差修正システム及び誤差修正方法 |
KR20010055220A (ko) * | 1999-12-10 | 2001-07-04 | 오길록 | 톱니 모양의 래크를 이용한 골소주 지표 보정방법 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20200102028A (ko) * | 2019-02-20 | 2020-08-31 | 하이윈 테크놀로지스 코포레이션 | 컴퓨터 단층 촬영장치의 스캔 정보를 골격 파라미터로 전환하는 방법 |
Also Published As
Publication number | Publication date |
---|---|
US20040114726A1 (en) | 2004-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR100343777B1 (ko) | 톱니 모양의 래크를 이용한 골소주 지표 보정방법 | |
Ortman et al. | Subtraction radiography and computer assisted densitometric analyses of standardized radiographs: a comparison study with 125I absorptiometry | |
Blake et al. | Applications of bone densitometry for osteoporosis | |
Rosholm et al. | Estimation of bone mineral density by digital X-ray radiogrammetry: theoretical background and clinical testing | |
Steiger et al. | Morphometric X-ray absorptiometry of the spine: correlation in vivo with morphometric radiography | |
JPH0331061B2 (ko) | ||
US7623618B2 (en) | Method for scattered radiation correction in X-ray imaging, and X-ray imaging system for this purpose | |
CN107485405B (zh) | 一种利用参考模块测量骨密度的装置 | |
KR100479207B1 (ko) | 엑스선 영상을 이용한 골밀도 측정 방법 | |
KR101087137B1 (ko) | 악골밀도 측정 방법 및 악골밀도 측정을 위한 보조장치 | |
KR20040051938A (ko) | 골밀도 지수 보정 방법 및 이 프로그램을 기록한 컴퓨터로읽을 수 있는 기록매체 | |
Maloul et al. | The impact of voxel size-based inaccuracies on the mechanical behavior of thin bone structures | |
Böttcher et al. | Influence of image-capturing parameters on digital X-ray radiogrammetry | |
KR20020046453A (ko) | 엑스선 영상을 이용한 해면골 발달정도 측정 방법 | |
US7046834B2 (en) | Method for measuring bone mineral density by using X-ray image | |
Atmaca et al. | Axial plane analysis of Lenke 1A adolescent idiopathic scoliosis as an aid to identify curve characteristics | |
Horiguchi et al. | Coronary artery calcium scoring on low-dose prospective electrocardiographically-triggered 64-slice CT | |
US9271690B2 (en) | Method and system to estimate visceral adipose tissue by restricting subtraction of subcutaneous adipose tissue to coelom projection region | |
US7734080B2 (en) | Method and device for measuring pulmonary nodule volume using dual-energy substraction image | |
Van Rijn et al. | Radiographic absorptiometry of the middle phalanx (digit II) in a Caucasian pediatric population: normative data | |
Tung et al. | Skeletal and dental relationships in vertical/non-vertical growers using CBCT | |
US9254101B2 (en) | Method and system to improve visceral adipose tissue estimate by measuring and correcting for subcutaneous adipose tissue composition | |
CN114782375B (zh) | 骨质密度测量方法、装置和设备 | |
JPS639462B2 (ko) | ||
KR20040058150A (ko) | 골 충실도를 이용한 골밀도 측정 방법 및 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A201 | Request for examination | ||
E902 | Notification of reason for refusal | ||
E601 | Decision to refuse application |