KR100942699B1 - 요골원위부 골간단 경계추출 방법 및 장치 - Google Patents
요골원위부 골간단 경계추출 방법 및 장치 Download PDFInfo
- Publication number
- KR100942699B1 KR100942699B1 KR1020070127883A KR20070127883A KR100942699B1 KR 100942699 B1 KR100942699 B1 KR 100942699B1 KR 1020070127883 A KR1020070127883 A KR 1020070127883A KR 20070127883 A KR20070127883 A KR 20070127883A KR 100942699 B1 KR100942699 B1 KR 100942699B1
- Authority
- KR
- South Korea
- Prior art keywords
- boundary
- distal
- landmark
- distal radius
- pixel
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 70
- 238000005381 potential energy Methods 0.000 claims abstract description 41
- 238000000137 annealing Methods 0.000 claims abstract description 26
- 238000000605 extraction Methods 0.000 claims abstract description 24
- 238000005452 bending Methods 0.000 claims description 10
- 210000002303 tibia Anatomy 0.000 claims description 9
- 210000001185 bone marrow Anatomy 0.000 claims description 5
- 210000004872 soft tissue Anatomy 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims 1
- 238000007781 pre-processing Methods 0.000 claims 1
- 210000000988 bone and bone Anatomy 0.000 abstract description 43
- 230000008569 process Effects 0.000 description 8
- 230000012010 growth Effects 0.000 description 7
- 238000005259 measurement Methods 0.000 description 7
- 241000270295 Serpentes Species 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 230000008468 bone growth Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000035800 maturation Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000002745 epiphysis Anatomy 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 210000000623 ulna Anatomy 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Geometry (AREA)
- Optics & Photonics (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
본 발명에 의한 요골원위부 골간단 경계를 추출할 수 있는 방법 및 장치는, 엑스선 영상에서 요골원위부를 포함하는 관심영역을 설정하고, 그레이 레벨의 그래디언트를 이용하여 상기 관심 영역에 대한 포텐셜 에너지 분포를 설정하고, 상기 관심영역에서 요골원위부 골간단의 양 측면과 골단쪽에 근접한 임시경계를 설정한 후, 메트로폴리스 어닐링 방법을 이용하여 상기 설정된 임시 경계를 조정하여, 최소의 에너지를 갖는 요골원위부 골간단 상세 경계를 추출하도록 구현된다.
골연령, 엑스선 영상, 메트로폴리스 어닐링, 그래디언트, 경계 추출
Description
본 발명은 엑스선 영상으로부터 어린이의 성인 신장 예측에 이용되는 요골원위부의 골간단 경계를 추출하는 방법 및 장치에 관한 것이다.
최근 어린이 성장에 관심이 많아지면서, 어린이의 성인 신장을 예측할 수 있는 임상적 방법 개발과, 성인 신장을 증가시킬 수 있는 치료 방법의 개발은 소아 성장연구 분야의 주요 관심사가 되었다.
이 중에서, 어린이의 성인 신장 예측은 해당 어린이의 골연령 측정에 기반하여 이루어진다. 여기서 골연령이란 해당 어린이의 물리적 연령과 다른 생리학적인 골성숙 연령을 의미하는 것으로, 물리적 연령과 반드시 일치하는 것은 아니며, 골성숙이 완료된 상태의 골연령은 누구나 같도록 되어있다.
따라서 골연령은 해당 어린이의 성장 잔여 시간에 대한 척도로 사용되어, 골연령이 물리적 연령보다 많으면 성장 잔여 시간이 짧아서 향후 해당 아동의 성장 여력은 작은 것으로 판단한다. 그리고 골연령 측정을 통하여 성인 신장이 작을 것 으로 예측되는 어린이는 골성숙이 완료되기 전에 성장 치료를 통하여 성인 신장을 증가시킬 수 있다.
현재 대표적으로 사용되는 골연령 측정 방법으로는, 아틀라스에 기반한 Greulich-Pyle (GP) 방법과, 골의 성숙도 정도를 점수로 표기하는 기법인 Tanner-Whitehouse (TW) 방법이 있다. 이 중에서 TW 방법이 좀 더 객관적이어서 현재 널리 활용되고 있으며, 상기 TW 방법은 TW2로 개선되고, 2000년에는 TW3가 발표되었다.
TW3 방법에서는 엑스선 영상에 나타난 골 성장점의 성숙 점수를 합하여 RUS(radius, ulna, short bones) 점수와 손목뼈(carpal) 점수를 얻는다. 각 골 성장점 성숙 단계는 보통 A-I 의 9단계로 이루어지며, 각 골성숙 단계에는 상응하는 성숙 점수가 부여되어 있다.
보통 임상 의사는 TW3 판정 방법에 기반하여, 손의 엑스선 사진으로부터 골의 성숙 단계를 판독하고, 이를 성숙 점수로 환산하고, 각 골의 성숙 점수를 합하여 RUS 점수나 carpal 점수를 계산한 후 해당 아동의 골연령을 측정한다.
기존의 골연령 측정 방법은, 임상 의사의 주관적인 판단에 의하여 이루어지기 때문에, 임상 의사 간에 판독 오차가 발생할 수 있으며, 동일한 의사도 시간에 따라 상이한 판독 결과를 제공할 수 있다.
따라서 컴퓨터 알고리즘을 통하여 골성숙 단계를 객관적으로 계측하는 장비의 개발이 요구되었으며, 현재까지 몇몇 컴퓨터 알고리즘을 이용한 골성숙 단계 계측 장비가 개발되어 활용되고 있다.
그러나, 현재 임상에서 활용되고 있는 골연령 측정 관련 장비의 대부분은, 획득된 엑스선 영상으로부터 의사가 TW3 골연령 측정을 하는 과정을 용이하게 도와주는 도구만 제공할 뿐 스스로 골성숙 단계를 계측할 수 없거나, 골성숙 단계를 계측하는 기능이 포함되어 있더라도 의사의 주관적 개입이 필요하거나, 골성숙 계측 과정이 자동화한 되어 있더라도 계측 결과의 정확성을 보장할 수 없는 경우가 많다.
더 구체적으로 설명하면, 골성숙 단계의 판정은, 골 성장점을 형성하는 골단골(epiphysis)과 골간단(metaphysis)의 형상 변화와 두 골의 융합 상태를 적절히 분류하고 평가하여 이루어진다. 따라서 정확한 골성숙 단계의 판정을 위해서는, 골 성장점을 형성하는 골단골과 골간단의 경계를 정확히 추출하여야 한다.
그런데, 기존의 방법에 따르면, 엑스선 영상에서 골단골과 골간단의 경계를 추출하는 것은, 다양한 경계 검출(edge detection)용 마스크를 적용하거나 그레이레벨 문턱치(thresholding) 설정 방법에 의해 이루어졌다. 그러나 실제 엑스선 영상을 이용하여 상기와 같은 방법으로 경계를 검출할 경우 많은 잡신호가 섞여 있어 그 검출 결과에 대한 신뢰도가 떨어진다는 문제점이 있다.
이외에, 엑스선 영상에서 잡신호의 영향을 줄이기 위하여 active shape model (ASM)이나 snake 와 같은 방법을 사용하여 골 경계선을 검출하기도 하는데, 이 경우는 많은 표본 영상을 사전에 획득하여 임상 의사의 수작업을 통하여 골경계 데이터베이스를 구축해야 하거나, snake에 곡률 특성을 주는 방법에 따라서 골경계 검출 성능이 달라지지며, 실제로 다양한 골경계에 어울리는 snake를 만들기도 어렵다는 문제점이 있다.
상기와 같이, 기존의 엑스선 영상으로부터 골연령 측정에 중요한 기준이 되는 요골원위부의 골간단 경계를 추출하는 방법은, 추출결과에 대한 정확도가 떨어져, 골연령 측정에 대한 신뢰도가 저하된다는 문제점이 있었다.
이에 본 발명은, 요골원위부 엑스선 영상에서 골간단의 경계를 정확하고 간단하게 추출하기 위한 것으로서, 그레이 레벨의 그래디언트(gradiant) 및 메트로폴리스 어닐링 방법을 이용하여, 잡신호가 많은 엑스선 영상으로부터 가장 신뢰성 있는 요골원위부 골간단 경계를 추출할 수 있는 방법 및 장치를 제공하고자 한다.
상기 과제를 해결하기 위한 수단으로서, 본 발명의 일 실시 형태에 의한 요골원위부 골간단 경계추출 방법은, 엑스선 영상에서 요골원위부를 포함하는 관심영역을 설정하는 단계; 그레이레벨의 그래디언트를 이용하여 상기 관심 영역에 대한 포텐셜에너지 분포를 설정하는 단계; 상기 관심영역에서 요골원위부 골간단의 양 측면과 골단쪽에 근접한 임시경계를 설정하는 단계; 및 메트로폴리스 어닐링 방법을 이용하여 상기 설정된 임시 경계를 조정하여, 최소의 에너지를 갖는 요골원위부 골간단 상세 경계를 추출하는 단계를 포함하여 이루어진다.
더하여, 본 발명의 다른 실시 형태에 의한 요골원위부 골간단 경계추출 장치는, 엑스선 영상에서 요골원위부를 포함하는 관심영역을 설정하는 관심영역 설정 부; 상기 관심영역에서 그레이레벨의 그래디언트로부터 포텐셜에너지 분포를 설정하는 포텐셜에너지 설정부; 상기 관심영역에서 요골원위부 골간단의 임시경계를 설정하는 전처리부; 및 메트로폴리스 어닐링 방법으로 상기 설정한 임시 경계를 조정하여 요골원위부 골간단 상세 경계를 추출하는 경계 추출부를 포함하여 구성된다.
본 발명에 따르면 손의 엑스선 영상을 이용하여 요골 원위부 골간단 경계를 추출하고자 할 때, 그레이 레벨의 그래디언트 및 메트로폴리스 어닐링 방법에 의하여 요골 원위부 골간단 경계를 간단하고, 정확하게 추출할 수 있으며, 그 결과, 어린이 성장 예측에 중요한 지표가 되는 골연령을 보다 정확히 측정할 수 있는 우수한 효과가 있다.
이하 첨부된 도면을 참조하여 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자가 본 발명을 용이하게 실시할 수 있는 바람직한 실시 예를 상세히 설명한다. 다만, 본 발명의 바람직한 실시 예에 대한 동작 원리를 상세하게 설명함에 있어 관련된 공지 기능 또는 구성에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우에는 그 상세한 설명을 생략한다.
또한, 도면 전체에 걸쳐 유사한 기능 및 작용을 하는 부분에 대해서는 동일한 도면 부호를 사용한다.
덧붙여, 명세서 전체에서, 어떤 부분이 다른 부분과 '연결'되어 있다고 할때, 이는 '직접적으로 연결'되어 있는 경우뿐만 아니라, 그 중간에 다른 소자를 사이에 두고 '간접적으로 연결'되어 있는 경우도 포함한다. 또한 어떤 구성 요소를 '포함'한다는 것은, 특별히 반대되는 기재가 없는 한 다른 구성요소를 제외하는 것이 아니라, 다른 구성요소를 더 포함할 수 있는 것을 의미한다.
도 1은 본 발명의 바람직한 실시 예에 따른 요골원위부 골간단 경계 추출 방법을 도시하는 흐름도이다.
도 1을 참조하면, 본 발명에 따른 요골 원위부 골간단 경계 추출 방법은, 먼저 단계 11에서, 엑스선 영상을 입력받아, 상기 엑스선 영상으로부터 요골원위부를 포함하는 관심영역을 설정한다. 상기 엑스선 영상은 일반적으로 이용되는 아날로그 필름 방식의 엑스선 촬영 장치나 디지털 센서 방식의 엑스선 촬영 장치로 획득할 수 있다. 본 발명에 따른 요골 원위부 골간단 경계추출 방법은 엑스선 영상에 나타난 요골 원위부와 인근 연부 조직의 밝기 차이를 이용한 것으로서, 상기 엑스선 촬영 장치를 이용한 엑스선 촬영은 요골 원위부와 인근 연부조직의 대조도가 커지도록 실시하는 것이 바람직하다.
도 2는 엑스선 촬영 영상(21)과 요골원위부 주위에 설정한 관심영역(22)의 일 예를 나타낸다. 상기 도 2에 예시된 바와 같이, 요골(20)이 영상의 한 축에 나란히 나타나도록 엑스선 촬영하는 것이 바람직하며, 관심 영역(22)의 중심은 대략 상기 요골(20)의 원위부 골간단 중앙부에 놓이도록 설정하는 것이 바람직하다. 한편, 상기 도 2의 예시에서, 엑스선 영상(21)의 공간 해상도는 200PPI(pixels per inch)로 하였으며, 그레이 레벨은 256(8-비트) 레벨로 하였다. 관심영역은 236ㅧ236 픽셀로 하였다.
다시 도 1을 참조하면, 본 발명에 따른 요골원위부 골간단 경계추출 방법은, 단계 12에서, 상기 설정한 관심영역에 대하여 그레이 레벨의 그래디언트로부터 각 픽셀의 포텐셜 에너지를 설정하여, 상기 관심영역에서의 포텐셜에너지 분포를 구한다.
여기서, 그래디언트의 계산은 3ㅧ3 Sobel 필터(S)를 이용하는 것이 바람직하다. 또한 잡신호를 더욱 줄이기 위하여 크기를 키운 Sobel 필터를 이용하는 것도 가능하다. 상기 Sobel 필터(S)를 이용함에 있어서, 각 픽셀에서의 그래디언트의 크기는, x 방향 Sobel 필터값(Sx)과 y 방향 Sobel 필터값(Sy)의 조합 으로 설정하고, 상기 픽셀에서의 포텐셜에너지는 로 설정한다. 이렇게 하면 그래디언트가 큰 픽셀(즉, 요골과 연부 조직의 경계)에서 포텐셜에너지가 낮아지는 역학계를 만들 수 있다.
도 3은 상기 단계 12에 의해서, 관심영역(22)에서 설정한 포텐셜 에너지가 나타난 영상을 예시한다. 상기와 같이 그래디언트가 큰 픽셀에서 포텐셜 에너지가 낮아지도록 설정한 영상으로서, 영상 내에서 밝기가 어두워질수록 포텐셜에너지가 낮아진다.
다시 도 1을 참조하면, 본 발명에 따른 요골원위부 골간단 경계추출 방법, 그 다음 단계 13에서, 상기 관심영역의 포텐셜 에너지 영상으로부터 요골원위부 골간단의 양 측면 경계와 골단쪽 임시경계를 설정한다.
도 4는 본 발명의 바람직한 일 실시 예에 따른, 상기 단계 13에서의 임시 경계를 설정하는 과정에 대한 상세 흐름도이다.
도 4를 참조하면, 상기 단계 13에서의 전처리 과정은, 먼저 단계 131을 통해, 상기 관심영역 내에서 요골이 배치된 축(예를 들어, X축)에 평행한 관심영역의 양측 경계에 두 개의 직선을 설정하고, 상기 골단이 배치되는 축(예를 들어 Y축)에 평행한 관심영역의 경계선에 또 다른 직선을 설정하고, 그 다음 단계 132에서, 상기 세 직선의 중심을 관심영역의 중심 방향으로 평행이동하면서, 동시에 각 직선을 회전시켜, 상기 세 직선 각각에 대하여 각 직선상에 놓인 픽셀들의 평균 포텐셜에너지가 최소가 되는 위치와, 회전각을 찾는다. 즉, 배경과 요골의 경계에서 포텐셜 에너지가 최소가 되므로, 상기 세 직선이 각각 요골원위부의 경계에 가장 근접하게 놓일 수 록, 그 평균 포텐셜 에너지는 최소가 된다. 따라서, 단계 133에서, 평균 포텐셜 에너지가 최소가 되는 위치와 회전각을 갖는 세 직선으로 둘러싸인 경계를 요골원위부 골간단 임시경계로 설정한다.
도 5는 상기 단계 13에 의한 전처리 과정의 바람직한 실시 예를 나타낸 것으로서, 부호 51a, 51b, 51c는 단계 131에서 각각의 초기 위치에 놓인 세 직선을 가 리키고, 부호 52a, 52b, 52c는 단계 132에 의해서 이동 및 회전되어 최소의 평균 포텐션 에너지를 갖는 위치에 놓인 세 직선을 가리키며, 53은 상기 세 직선(52a,52b,52c)에 의해 만들어진 임시 경계를 가리킨다.
다시 도 1을 참조하면, 본 발명에 따른 요골원위부 골간단 경계추출 방법은, 단계 14에서, 상기 설정한 임시경계를 메트로폴리스 어닐링 방법으로 처리하여, 요골원위부 골간단의 상세 경계를 추출한다.
도 6은 상기 메트로폴리스 어닐링 방법에 의한 요골원위부 골간단 경계의 상세 추출 과정의 바람직한 실시 예를 도시하는 흐름도이다.
도 6을 참조하면, 본 발명은 요골원부 골간단 상세 경계를 추출하기 위하여, 먼저 단계 141에서, 상기 설정한 임시경계를 분할하여 다수의 랜드마크 픽셀을 관심영역에 설정한다. 도 7은 다수(N)의 랜드마크 픽셀의 설정 예를 도시한 것으로서, 포텐셜 에너지 분포를 나타낸 영상에 설정된 임시 경계를 일정 간격으로 나누어, 상기 임시 경계 상의 해당 픽셀을 랜드마크 픽셀(71)로 설정하였다. 편의상 임시 경계의 양 끝에 위치한 두 개의 랜드마크(1번 및 N번)는 관심영역의 경계선을 따라서 이동하도록 제한을 가한다.
그리고, 단계 142에서, 상기 각 랜드마크 픽셀의 에너지(E)를 설정한다. 여기서, 랜드마크 픽셀의 에너지를 설정하는 방법은 다양하게 이루어질 수 있는데, 예를 들면, 해당 랜드마크 픽셀의 포텐셜에너지와 벤딩에너지의 합으로 설정한다. n번째 랜드마크 픽셀의 벤딩에너지는 전후(n-1,n+1) 랜드마크 픽셀들이 형성하는 각도(θ)를 이용하여 로 설정할 수 있다. 이 벤딩에너지는 포텐셜에너지의 감소가 없는 상태에서 랜드마크 픽셀에 굴곡이 생기는 것을 방지하기 위한 일종의 패널티이다. n번째 랜드마크 픽셀의 벤딩에너지의 계수()는 상수로 설정할 수도 있고, 전후에 위치한 n-1, n, n+1 번째 랜드마크 픽셀의 평균 포텐셜에너지에 비례하도록 설정할 수도 있다. 도 7의 실시 예에서는 후자의 방법을 택한다. 그리고 나서 단계 143에서 메트로폴리스 어닐링을 시작하기 위한 최고온도()를 설정한다. 상기 는 최저 포텐셜에너지보다 큰 값으로 하는 것이 바람직하며, 본 실시 예에서는 편의상 최고 포텐셜에너지의 5배로 설정한다.
메트로폴리스 어닐링은 에너지가 주어진 물리계에서 온도를 서서히 낮추면서 그 온도에서 가장 보편적인 물리적 상태를 만들어 해당 물리계의 낮은 에너지 상태의 구조를 찾는 방법이다. 이 방법은 해당 물리계가 국소점에 빠져 실질적 최소 상태에 이르지 못하는 것을 방지하며, 어떤 온도에서 가장 보편적인 물리적 상태를 만들어 낸다. 메트로폴리스 방법에서는 온도 T인 물리계의 상태를 변화시키려 할 때, 에너지 변화가 ΔE이면 max[,1]의 확률로 상태 변화를 확정하여 가장 보편적인 물리적 상태를 만들어낸다.
따라서, 상기와 같이 최고 온도 가 설정되면, 단계 144에서 메트로폴 리스 어닐링 온도 T를 로 한 상태에서, 무작위 혹은 순서대로 선택한 랜드마크 픽셀의 위치를 인근으로 이동시키되, 시스템의 에너지 변화를 계산하여 메트로폴리스 방법으로 이동을 확정한다. 그리고 단계 145에서, 상기 단계 144의 과정을 설정된 횟수(총 랜드마크 픽셀수(N)의 다수 배)만큼 실시하도록 반복시키고, 상기 설정횟수만큼 이동이 이루어진 후에는, 단계 146에서 온도 T를 사전에 설정된 단위 온도(ΔT) 만큼 낮춘 후에, 단계 147에서 현재의 메트로폴리스 어닐링 온도 T를 확인하여, 최저 온도가 될 때까지 온도를 낮춰가면서, 상기 단계 144를 반복시킨다.
상기 단계 144에서 랜드마크의 위치를 인근으로 이동시킬 때 전후 두 개의 랜드마크가 만드는 직선의 수직방향(도 7의 화살표 방향)으로 이동하되, 이동 위치는 현재 랜드마크 위치에서 사전에 설정한 거리 이내의 픽셀을 무작위 선택한다.
상기 단계 144를 반복실행함에 의하여, 상기 설정된 각 랜드마크 픽셀은 최저 에너지를 갖는 픽셀로 이동되어, 가장 유사한 요골원위부 골간단의 경계를 추출할 수 있게 된다.
도 8은, 상술한 메트로폴리스 어닐링 방법에 의하여 상기 단계 13에서 설정된 임시 경계를 조정하여, 요골원위부 골간단 상세 경계를 추출한 결과의 일 예를 도시한 것이다.
앞서 언급한 바와 같이 기존의 ASM을 이용하여 요골원위부 골간단 경계를 추출하려는 경우에는, 상당한 양의 요골원위부 골간단 경계추출 데이터베이스를 갖추고 있어야 하며, snake 등을 이용하여 요골원위부 골간단을 추출할 경우에는 snake 의 특성에 따라 골간단 경계 검출 성능이 달라질 수 있었으나, 본 발명에 의한 요골원위부 골간단 경계추출 방법에 의하면, 별도의 데이터베이스없이 엑스선 영상에서 물리적으로 가장 그럴법한 경계선을 찾도록 함으로써, 간편하면서도 뛰어난 성능의 골간단 경계추출 성능을 보여준다.
도 9는 본 발명의 바람직한 실시 예에 따른 요골원위부 골간단 경계추출 장치를 도시하는 블록도이다.
도 9를 참조하면, 본 발명에 따른 요골원위부 골간단 경계추출 장치는, 손의 엑스선 영상으로부터 요골원위부를 포함하는 관심영역을 설정하는 관심영역 설정부(201)와, 상기 관심영역의 각 픽셀의 포텐셜 에너지를 그레이레벨의 그래디언트로부터 설정하여 포텐셜에너지 분포 영상을 만드는 포텐셜에너지 설정부(202)와, 상기 관심영역의 포텐셜 에너지 분포 영상에서 요골원위부 골간단의 임시 경계를 설정하는 전처리부(203)와, 상기 포텐셜 에너지 분포 영상의 임시 경계를 메트로폴리스 어닐링 방법으로 조정하여 에너지가 가장 낮은 부분으로 이루어지는 요골원위부 골간단 상세 경계를 추출하는 경계 추출부(204)를 포함한다.
상술한 요골원위부 골간단 경계 추출 장치에서의 동작은, 앞서 설명한 요골원위부 골간단 경계 추출 방법에 따라 이루어지는 것으로서, 불필요한 중복 설명을 피하기 위하여 상기 동작의 상세 설명은 생략한다.
본 발명에 의한 요골원위부 골간단 경계추출 장치는, 컴퓨터로 읽을 수 있는 기록매체에 저장되는 컴퓨터가 읽을 수 있는 코드로 구현될 수 있다. 즉, 본 발명 의 요골원위부 골간단 경계 추출 장치는, 컴퓨터 프로그램으로 구현되어, 손의 엑스선 촬영 영상 데이터로부터 골연령 측정에 필요한 요골원위부 골간단 경계를 비교적 간단하고 정확하게 추출할 수 있다. 상기에서, 컴퓨터가 읽을 수 있는 기록매체는 컴퓨터 시스템에 의하여 읽혀질 수 있는 데이터가 저장되는 모든 종류의 기록장치를 포함한다. 컴퓨터가 읽을 수 있는 기록매체의 예로는 ROM, RAM, CD_ROM, 자기테이프, 플로피디스크 및 광 데이터 저장장치 등이 있으며, 또한 캐리어 웨이브(예를 들어 인터넷을 통한 전송)의 형태로 구현되는 것도 포함한다. 또한 컴퓨터가 읽을 수 있는 기록매체는 네트워크로 연결된 컴퓨터 시스템에 분산되어, 분산방식으로 컴퓨터가 읽을 수 있는 코드로 저장되고 실행될 수 있다.
이상에서 설명한 본 발명은 전술한 실시 예 및 첨부된 도면에 의해 한정되는 것이 아니고, 본 발명의 기술적 사상을 벗어나지 않는 범위 내에서 여러 가지 치환, 변형 및 변경할 수 있다는 것은 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 당업자에게 있어 명백할 것이다.
도 1은 본 발명의 바람직한 실시 형태에 따른 요골원위부 골간단 경계 추출 방법을 나타낸 흐름도,
도 2는 본 발명에 의한 요골원위부 골간단 경계 추출에 이용되는 엑스선 영상 및 관심영역 설정 예를 보인 예시도,
도 3은 본 발명에 의한 포텐셜 에너지 분포 설정 예를 보인 예시도,
도 4는 본 발명의 바람직한 실시 형태에 따른 요골원위부 골간단 경계 추출 방법에 있어서, 임시 경계를 설정하는 상세 과정을 나타낸 흐름도,
도 5는 본 발명에 의한 임시 경계 설정 예를 보인 예시도,
도 6은 본 발명의 바람직한 실시 형태에 따른 요골원위부 골간단 경계 추출 방법에 있어서, 요골원위부 골간단 상세 경계의 설정 과정을 구체적으로 나타낸 흐름도,
도 7은 본 발명에 의한 요골원위부 골간단 상세 경계의 설정에 있어서, 랜드마크 픽셀의 설정 예를 도시한 예시도,
도 9은 본 발명에 의한 요골원위부 골간단 상세 경계 추출 결과의 일 예를 도시한 예시도, 그리고
도 9는 본 발명의 바람직한 실시 예에 따른 요골원위부 골간단 경계 추출 장치를 도시하는 블록도이다.
Claims (17)
- 엑스선 영상에서 요골원위부를 포함하는 관심영역을 설정하는 단계;그레이레벨의 그래디언트를 이용하여 상기 관심 영역에 대한 포텐셜에너지 분포를 설정하는 단계;상기 관심영역에서 요골원위부 골간단의 양 측면과 골단쪽에 근접한 임시경계를 설정하는 단계; 및메트로폴리스 어닐링 방법을 이용하여 상기 설정된 임시 경계를 조정하여, 최소의 에너지를 갖는 요골원위부 골간단 상세 경계를 추출하는 단계를 포함하는 요골원위부 골간단 경계추출 방법.
- 제 1 항에 있어서,상기 엑스선 영상은요골 원위부와 인근 연부조직의 밝기 차가 크도록 촬영된 것임을 특징으로 하는 요골원위부 골간단 경계추출 방법.
- 제 1 항에 있어서,상기 관심 영역에 대한 포텐셜 에너지 분포를 설정하는 단계는,상기 관심 영역 내 각 픽셀의 포텐셜 에너지는 그레이 레벨의 그래디언트의 절대값에 -1을 곱한 값으로 설정함에 의해 이루어지는 것을 특징으로 하는 요골원 위부 골간단 경계추출 방법.
- 제 1 항에 있어서, 상기 임시 경계를 설정하는 단계는,상기 관심영역의 경계선중 요골원위부 양 측면의 위치한 2개의 경계선과, 골단쪽에 위치한 경계선을 세 개의 직선으로 설정하는 단계;상기 설정한 세 개의 직선 각각에 대하여, 직선의 중심을 상기 관심영역의 중심 방향으로 평행이동하면서 동시에 직선을 회전시켜 직선상에 위치한 픽셀들의 포텐셜에너지 평균값을 구하여, 직선의 포센셜 에너지 평균값이 최소가 되는 위치와 회전각을 찾는 단계; 및상기 위치와 회전각을 찾는 단계를 통하여 찾은 위치 및 회전각을 갖는 세 개의 직선으로 둘러싸인 부분을 요골원위부 골간단에 대한 임시경계로 설정하는 단계를 포함하는 요골원위부 골간단 경계추출 방법.
- 삭제
- 제 1 항에 있어서, 상기 요골원위부 골간단 상세 경계를 추출하는 단계는,상기 임시 경계를 분할하여, 임시 경계 선상에 위치한 다수의 랜드마크 픽셀을 설정하는 단계;상기 각 랜드마크 픽셀의 에너지를 설정하는 단계;메트로폴리스 어닐링을 시작하기 위한 최고온도를 설정하는 단계;무작위 혹은 순서대로 복수의 랜드마크 픽셀을 선택하여, 메트로폴리스 어닐링 방법에 의해 이동시키는 단계;사전에 설정된 횟수가 도달되었는 지를 확인하여, 설정 횟수에 도달할때까지 상기 복수의 랜드마크 픽셀의 선택 및 메트로폴리스 어닐링 방법에 의한 이동을 반복하여 수행하도록 하는 단계;상기 사전에 설정된 횟수가 도달하면, 상기 메트로폴리스 어닐링 방법을 통해서 구해진 매트로폴리스 어닐링 온도를 단위 온도만큼 감소시키는 단계; 및상기 메트로폴리스 어닐링 방법을 통해서 구해진 메트로폴리스 어닐링 온도가 사전에 설정한 최저온도 미만인지를 확인하여, 최저 온도 미만이 아니면, 상기 다수의 랜드마크 픽셀을 선택하여, 메트로폴리스 어닐링 방법에 의해 이동시키는 단계부터 반복 수행하고, 최저 온도 미만이면 현재 상태의 랜드마크 픽셀들을 연결한 경계를 요골원위부 골간단 상세 경계로 추출하는 단계를 포함하는 요골원위부 골간단 경계추출 방법.
- 제 6 항에 있어서,상기 다수의 랜드마크 픽셀에서, 임시 경계의 양 끝에 위치한 두 개의 랜드마크 픽셀은 관심영역의 경계선을 따라서 이동하도록 설정하는 것을 특징으로 하는 요골원위부 골간단 경계추출 방법.
- 제 6 항에 있어서,상기 최고온도는 최저 포텐셜에너지보다 큰 값으로 설정하는 것을 특징으로 하는 요골원위부 골간단 경계추출 방법.
- 제 6 항에 있어서,상기 다수의 랜드마크 픽셀을 선택하여, 메트로폴리스 어닐링 방법에 의해 이동시키는 단계는,선택한 랜드마크 픽셀의 전후 랜드마크 픽셀이 만드는 직선의 수직방향으로 상기 선택한 랜드마크 픽셀을 이동시키되, 현재 위치에서 사전에 설정한 거리 이내의 픽셀을 무작위로 선택하는 것을 특징으로 하는 요골원위부 골간단 경계추출 방법.
- 엑스선 영상에서 요골원위부를 포함하는 관심영역을 설정하는 관심영역 설정부;상기 관심영역에서 그레이레벨의 그래디언트로부터 포텐셜에너지 분포를 설정하는 포텐셜에너지 설정부;상기 관심영역에서 요골원위부 골간단의 임시경계를 설정하는 전처리부; 및메트로폴리스 어닐링 방법으로 상기 설정한 임시 경계를 조정하여 요골원위부 골간단 상세 경계를 추출하는 경계 추출부를 포함하는 요골원위부 골간단 경계추출 장치.
- 제 13 항에 있어서, 상기 포텐셜에너지 설정부는,상기 관심 영역 내 각 픽셀의 포텐셜 에너지를, 상기 픽셀의 그레이 레벨의 그래디언트의 절대값에 -1을 곱한 값으로 설정하여 포텐셜에너지 분포를 구하는 것을 특징으로 하는 요골원위부 골간단 경계추출 장치.
- 제 13 항에 있어서, 상기 전처리부는,상기 관심영역에서 영역을 표시하는 경계선 4개 중 요골원위부 양 측면의 위치한 두 개의 경계선과 골단쪽 경계선을 세 개의 직선으로 설정하고, 설정된 세 개의 직선을 상기 관심영역의 중심 방향으로 평행이동하면서 동시에 회전시키면서, 설정된 세 개의 직선 각각에 위치한 픽셀들의 포텐셜에너지 평균값을 구하여, 설정된 세 개의 직선 각각의 포센셜 에너지 평균값이 최소가 될 때, 설정된 세 개의 직선으로 둘러싸인 부분을 임시 경계로 설정하는 것을 특징으로 하는 요골원위부 골간단 경계추출 장치.
- 제 13 항에 있어서, 상기 경계추출부는상기 임시 경계를 분할하여, 임시 경계 선상에 위치한 다수의 랜드마크 픽셀을 설정하고, 상기 랜드마크 픽셀의 에너지를 설정한 후, 사전에 설정한 최고온도에서 최저 온도의 범위에서 단위온도만큼 감소시키면서 메트로폴리스 어닐링에 의해 랜드마크 픽셀을 이동시켜, 요골원위부 요골원위부 골간단 상세 경계를 추출하는 것을 특징으로 하는 요골원위부 골간단 경계추출 장치.
- 제 16 항에 있어서, 상기 경계 추출부는상기 다수의 랜드마크 픽셀 각각의 전후 랜드마크 픽셀이 만드는 직선의 수직방향으로 상기 다수의 랜드마크 픽셀 각각을 이동시키되, 현재 위치에서 사전에 설정한 거리 이내의 픽셀을 무작위로 선택하는 것을 특징으로 하는 요골원위부 골간단 경계추출 장치.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020070127883A KR100942699B1 (ko) | 2007-12-10 | 2007-12-10 | 요골원위부 골간단 경계추출 방법 및 장치 |
US12/111,320 US7983461B2 (en) | 2007-12-10 | 2008-04-29 | Method and system for extracting distal radius metaphysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020070127883A KR100942699B1 (ko) | 2007-12-10 | 2007-12-10 | 요골원위부 골간단 경계추출 방법 및 장치 |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20090060900A KR20090060900A (ko) | 2009-06-15 |
KR100942699B1 true KR100942699B1 (ko) | 2010-02-17 |
Family
ID=40721726
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020070127883A KR100942699B1 (ko) | 2007-12-10 | 2007-12-10 | 요골원위부 골간단 경계추출 방법 및 장치 |
Country Status (2)
Country | Link |
---|---|
US (1) | US7983461B2 (ko) |
KR (1) | KR100942699B1 (ko) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201117804D0 (en) * | 2011-10-14 | 2011-11-30 | Siemens Medical Solutions | Automatic local contrast quantification tool |
US10867374B2 (en) * | 2017-02-13 | 2020-12-15 | Sensors Unlimited, Inc. | Auto-focusing system and method by determining contrast difference between adjacent pixels using sobel filter |
CN107563440B (zh) * | 2017-09-01 | 2020-07-28 | 京东方科技集团股份有限公司 | 一种图像处理方法及装置 |
US20210142477A1 (en) * | 2018-08-01 | 2021-05-13 | China Medical University Hospital | Bone Age Assessment And Height Prediction Model, System Thereof And Prediction Method Thereof |
US11741586B2 (en) * | 2018-09-04 | 2023-08-29 | Bonewise Inc. | Apparatus, method and recording medium storing program for assessing bone age |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07168944A (ja) * | 1993-05-28 | 1995-07-04 | South Birmingham District Health Authority | 画像処理方法およびその装置 |
KR20060008849A (ko) * | 2005-12-23 | 2006-01-27 | 학교법인 한양학원 | Tw2 영상의 기하학적 특성값 정의를 이용한 뼈 나이 측정방법 및 장치 |
KR20070077319A (ko) * | 2006-01-23 | 2007-07-26 | 길호석 | 자동적으로 rus 점수를 산출하여 성장을 예측하는 방법 |
KR20080052090A (ko) * | 2006-12-07 | 2008-06-11 | 한국전자통신연구원 | 수지골 경계 추출 영상 획득 방법 및 그 시스템 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AT390183B (de) * | 1985-07-16 | 1990-03-26 | Boehler Gmbh | Knochenimplantat fuer endoprothesen |
EP0228785B1 (en) * | 1985-11-11 | 1993-01-13 | Teijin Limited | Bone evaluation method |
US5931780A (en) * | 1993-11-29 | 1999-08-03 | Arch Development Corporation | Method and system for the computerized radiographic analysis of bone |
DE69513300T2 (de) * | 1994-08-29 | 2000-03-23 | Torsana A/S, Skodsborg | Bestimmungsverfahren |
WO1997042602A1 (en) * | 1996-05-06 | 1997-11-13 | Torsana A/S | A method of estimating skeletal status |
JP2004535882A (ja) | 2001-07-24 | 2004-12-02 | サンライト メディカル リミテッド | 超音波を用いた骨年齢評価方法 |
EP1293925A1 (en) * | 2001-09-18 | 2003-03-19 | Agfa-Gevaert | Radiographic scoring method |
US7769214B2 (en) * | 2002-12-05 | 2010-08-03 | The Trustees Of The University Of Pennsylvania | Method for measuring structural thickness from low-resolution digital images |
KR100686534B1 (ko) | 2004-11-03 | 2007-02-23 | 윤지현 | 증감지 비투과방식의 구조를 이용한 골밀도 진단장치 |
US7848893B2 (en) * | 2005-06-02 | 2010-12-07 | Hans Henrik Thodberg | Method for skeletal maturity determination |
-
2007
- 2007-12-10 KR KR1020070127883A patent/KR100942699B1/ko not_active IP Right Cessation
-
2008
- 2008-04-29 US US12/111,320 patent/US7983461B2/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07168944A (ja) * | 1993-05-28 | 1995-07-04 | South Birmingham District Health Authority | 画像処理方法およびその装置 |
KR20060008849A (ko) * | 2005-12-23 | 2006-01-27 | 학교법인 한양학원 | Tw2 영상의 기하학적 특성값 정의를 이용한 뼈 나이 측정방법 및 장치 |
KR20070077319A (ko) * | 2006-01-23 | 2007-07-26 | 길호석 | 자동적으로 rus 점수를 산출하여 성장을 예측하는 방법 |
KR20080052090A (ko) * | 2006-12-07 | 2008-06-11 | 한국전자통신연구원 | 수지골 경계 추출 영상 획득 방법 및 그 시스템 |
Also Published As
Publication number | Publication date |
---|---|
US7983461B2 (en) | 2011-07-19 |
US20090148022A1 (en) | 2009-06-11 |
KR20090060900A (ko) | 2009-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8965083B2 (en) | Automatic detection of vertebrae boundaries in spine images | |
US9594975B2 (en) | Image processing device and spinal canal evaluation method | |
NL1027538C2 (nl) | Werkwijzen en apparaten voor verwerken van beelddata om te helpen bij detecteren van ziekte. | |
US8411927B2 (en) | Marker detection in X-ray images | |
EP1881453A2 (en) | A medical image-processing apparatus and a method for processing medical images | |
US8649843B2 (en) | Automated calcium scoring of the aorta | |
WO2007117506A2 (en) | System and method for automatic detection of internal structures in medical images | |
KR100942699B1 (ko) | 요골원위부 골간단 경계추출 방법 및 장치 | |
US9940537B2 (en) | Medical image processing apparatus, medical imaging apparatus and medical image processing method | |
EP3047455B1 (en) | Method and system for spine position detection | |
CN105556567B (zh) | 用于脊椎位置探测的方法和系统 | |
CN104507392A (zh) | 图像处理装置及图像处理方法 | |
CN113826140B (zh) | 血管造影数据分析 | |
KR100853187B1 (ko) | 수지골 경계 추출 영상 획득 방법 및 그 시스템 | |
KR20180051267A (ko) | 골 등급 측정 장치 및 방법 | |
US20150254850A1 (en) | System for detecting blood vessel structures in medical images | |
Sha et al. | The improved faster-RCNN for spinal fracture lesions detection | |
CN104915989A (zh) | 基于ct影像的血管三维分割方法 | |
CN112184623A (zh) | 脊柱椎体的椎间隙分析方法、设备和存储介质 | |
CN109615656A (zh) | 一种基于模型搜索的脊柱定位方法 | |
CN105956587A (zh) | 一种基于形状约束的膝关节磁共振图像序列半月板自动提取方法 | |
CN113538352B (zh) | 脑卒中组织窗评价值的获取方法及装置、存储介质 | |
KR101750173B1 (ko) | 방사선 유효 피폭량 자동 계산 시스템 및 방법 | |
JP3639283B2 (ja) | 画像処理装置及び画像処理方法 | |
US20230103920A1 (en) | Device and method of angiography for cerebrovascular obliteration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A201 | Request for examination | ||
E902 | Notification of reason for refusal | ||
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant | ||
LAPS | Lapse due to unpaid annual fee |