KR102054210B1 - 치아 ct 영상을 이용한 자동 파노라마 영상 생성 방법 및 시스템 - Google Patents

치아 ct 영상을 이용한 자동 파노라마 영상 생성 방법 및 시스템 Download PDF

Info

Publication number
KR102054210B1
KR102054210B1 KR1020180026260A KR20180026260A KR102054210B1 KR 102054210 B1 KR102054210 B1 KR 102054210B1 KR 1020180026260 A KR1020180026260 A KR 1020180026260A KR 20180026260 A KR20180026260 A KR 20180026260A KR 102054210 B1 KR102054210 B1 KR 102054210B1
Authority
KR
South Korea
Prior art keywords
region
slide
image
generating
curve
Prior art date
Application number
KR1020180026260A
Other languages
English (en)
Other versions
KR20190105763A (ko
Inventor
채옥삼
류병용
김재면
김병우
정소현
Original Assignee
경희대학교 산학협력단
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 경희대학교 산학협력단 filed Critical 경희대학교 산학협력단
Priority to KR1020180026260A priority Critical patent/KR102054210B1/ko
Publication of KR20190105763A publication Critical patent/KR20190105763A/ko
Application granted granted Critical
Publication of KR102054210B1 publication Critical patent/KR102054210B1/ko

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • A61B6/14
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/51Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for dentistry

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

본 발명의 파노라마 영상 생성 시스템에서, 치아의 CT 영상을 이용하여 자동으로 상기 치아의 파노라마 영상을 생성하는 방법에 있어서, 상기 CT 영상에서 치아 영역을 추출하는 단계, 상기 치아 영역에 포함되는 치열 궁 곡선을 생성하는 단계, 상기 치아 영역과 상기 치열 궁 곡선의 최대 거리를 이웃 화소 범위로 설정하는 단계, 상기 치열 궁 곡선에 포함된 각 점에서 상기 이웃 화소 범위만큼 이격된 점으로 구성된 폐곡선을 생성하는 단계, 상기 CT 영상에 포함된 제1 슬라이드에 폐곡선과 상응하는 위치의 폐곡선 영역을 추출하는 a 단계, 상기 폐곡선 영역에서 상기 치열 궁 곡선을 축으로 휘어짐을 보정하여 직사각 영역을 생성하는 b 단계, 상기 직사각 영역에 포함된 픽셀에서 동일한 x 값을 갖는 하나 이상의 픽셀을 이용하여 상기 x 값에 대응되는 대표 픽셀을 생성하는 c 단계, 각 x 값에 대응되는 상기 대표 픽셀을 조합하여 상기 제1 슬라이드의 픽셀 라인을 생성하는 d 단계, 상기 CT 영상에 포함된 슬라이드 각각에 대하여 상기 a 내지 상기 d 단계를 수행하고, 각 슬라이드의 픽셀 라인을 합성하여 상기 파노라마 영상을 생성하는 단계를 포함하는 것을 특징으로 한다.

Description

치아 CT 영상을 이용한 자동 파노라마 영상 생성 방법 및 시스템{METHOD AND SYSTEM FOR AUTOMATIC GENERATION OF PANORAMIC IMAGE USING DENTAL CT}
본 발명은 치아 파노라마 영상 생성 방법 및 시스템에 관한 것으로, 보다 구체적으로 치아 CT 영상을 이용하여 자동으로 파노라마 영상을 생성하는 방법 및 시스템에 관한 것이다.
최근 치아를 치료하기 위해 사용되는 영상 기술이 빠르게 발전하고 있다. 의료 기술 분야는 사람의 신체에 직접적인 영향을 미치기 때문에 그 어떤 분야보다 보수적으로 개발됨에도 불구하고 치과 분야는 3D 스캐너, 3D 조각기, 프린터 및 치아 CT(Computed Tomography, 전산화 단층 촬영) 등의 기술을 활용한 디지털 기술을 도입하여 빠르게 발전하고 있다. 이러한 디지털 기술 중 치아 CT는 구강 내부의 정확한 3D 정보를 제공하기 때문에 치과 진료용 영상으로 사용되고 있다. 또한, 치아 CT는 치아 진단을 위한 2D 및 3D 시각화, 치아 교정치료, 임플란트, 보철 치료 등을 위한 CAD/CAM 생산 등 다양한 분야에도 널리 사용되고 있다.
나아가 치과 진료에 있어서 치아 CT의 보급과 가시화 소프트웨어의 발달로 2D 및 3D의 CT 영상 가시화를 통해 환자의 상태를 진단하고 문제가 있는 치아 부분을 설명하는 과정이 일반화되었다. 치아 CT 영상 가시화 소프트웨어를 사용한 치과 진료 과정에 있어서 가장 빈번하게 활용되는 기술은 치아 파노라마 영상을 재구성하는 것이다. 이 때, 치아 CT에서 획득한 3D 영상은 가시화 소프트웨어를 이용하여 X-Ray 파노라마 영상과 유사한 형태를 갖는 파노라마 영상으로 재구성될 수 있다.
치아 CT에서 파노라마 영상을 재구성하기 위해 사람마다 각기 다른 치열 궁(Dental Arch)의 위치를 수동으로 입력 받아야 한다. 치열 궁의 위치를 자동으로 입력 받기 위해 다양한 치열 궁을 수학적으로 표현하고 검출하는 연구가 진행되고 있다. 관련된 종래의 특허에는 대한민국 공개특허공보 제10-2016-0083788호의 '치아 파노라마 영상 생성방법, 이를 위한 장치 및 이를 기록한 기록매체'가 있다. 그러나 기존의 자동 파노라마 영상 생성 방법은 일반적인 상황에서는 적합한 결과를 도출하지만 부정교합이 심한 치아에서는 특정 치아가 파노라마 영상에 노출되지 않는다는 단점이 있다. 또한 자동으로 파노라마 영상을 생성함에 있어서, 치열 궁의 이웃 화소를 어느 범위까지 수용해야 하는 지에 대한 점이 전혀 고려되지 않는다는 단점이 있다. 이웃 화소를 적게 수용하면 치아가 뚜렷하게 보여지지만 수용되는 화소가 적기 때문에 치아의 충분한 구조를 보여주기 어려운 반면 이웃 화소를 많이 수용하면 치아가 다른 정보와 섞여 뚜렷하게 보이지 않지만 충분한 구조를 보여줄 수 있기 때문에 적절한 이웃 화소 범위를 예측하는 기술이 필요하다.
본 발명은 전술한 문제를 해결하기 위한 것으로, 치아 CT 영상에서 파노라마 영상을 자동으로 재구성하는 것을 일 목적으로 한다.
또한 본 발명은 치아 CT 영상에 포함된 슬라이드에서 영상에 따라 설정되는 임계 값을 이용하는 Otsu 이진화 알고리즘을 적용하여 치아 영역을 추출하는 것을 일 목적으로 한다.
또한 본 발명은 적합한 이웃 화소 범위를 설정하여 모든 치아가 정확하게 나타날 수 있도록 하는 파노라마 영상을 재구성하는 것을 일 목적으로 한다.
또한 본 발명은 매개변수 곡선을 이용하여 적절한 치열 궁 곡선을 생성하는 것을 일 목적으로 한다.
이러한 목적을 달성하기 위한 본 발명은 파노라마 영상 생성 시스템에서 치아의 CT 영상을 이용하여 자동으로 상기 치아의 파노라마 영상을 생성하는 방법에 있어서, 상기 CT 영상에서 치아 영역을 추출하는 단계, 상기 치아 영역에 포함되는 치열 궁 곡선을 생성하는 단계, 상기 치아 영역과 상기 치열 궁 곡선의 최대 거리를 이웃 화소 범위로 설정하는 단계, 상기 치열 궁 곡선에 포함된 각 점에서 상기 이웃 화소 범위만큼 이격된 점으로 구성된 폐곡선을 생성하는 단계, 상기 CT 영상에 포함된 제1 슬라이드에 폐곡선과 상응하는 위치의 폐곡선 영역을 추출하는 a 단계, 상기 폐곡선 영역에서 상기 치열 궁 곡선을 축으로 휘어짐을 보정하여 직사각 영역을 생성하는 b 단계, 상기 직사각 영역에 포함된 픽셀에서 동일한 x 값을 갖는 하나 이상의 픽셀을 이용하여 상기 x 값에 대응되는 대표 픽셀을 생성하는 c 단계, 각 x 값에 대응되는 상기 대표 픽셀을 조합하여 상기 제1 슬라이드의 픽셀 라인을 생성하는 d 단계, 상기 CT 영상에 포함된 슬라이드 각각에 대하여 상기 a 내지 상기 d 단계를 수행하고, 각 슬라이드의 픽셀 라인을 합성하여 상기 파노라마 영상을 생성하는 단계를 포함하는 것을 일 특징으로 한다.
또한 상기 치아 영역을 추출하는 단계는, 상기 CT 영상에 포함된 슬라이드 각각에 대하여 얼굴의 최하단을 기준으로 기 설정된 크기의 후보 영역을 추출하는 단계, 상기 후보 영역에 포함된 픽셀의 밝기 값을 이용하여 경조직을 추출하는 것을 일 특징으로 한다.
나아가 상기 슬라이드에서 추출된 상기 경조직에 대하여 동일한 위치에 포함된 픽셀의 밝기 값에 대한 평균 값을 연산하여 누적 슬라이드를 생성하는 단계, 상기 누적 슬라이드에 포함된 경조직의 영역 중 크기가 가장 큰 영역을 치아 영역으로 설정하는 단계를 포함하는 것을 일 특징으로 한다.
또한 상기 경조직을 추출하는 단계는 상기 후보 영역의 밝기 값을 Otsu 알고리즘을 이용하여 이진화하는 것을 특징으로 하는 것을 일 특징으로 한다.
나아가 상기 누적 슬라이드를 생성하는 단계는, 상기 슬라이드 각각에서 추출된 경조직의 상응하는 위치에 포함된 픽셀의 밝기 값에 대한 평균 값을 연산하여 누적 슬라이드를 생성하는 단계, 상기 누적 슬라이드를 Otsu 알고리즘을 이용하여 이진화하는 단계를 포함하는 것을 일 특징으로 한다.
또한 상기 대표 픽셀을 생성하는 단계는 상기 직사각 영역을 픽셀 단위로 분할하여 동일한 x 값을 갖는 하나 이상의 픽셀의 평균 값을 연산하는 단계, 상기 평균 값을 상기 대표 픽셀의 밝기 값으로 설정하는 단계를 포함하는 것을 일 특징으로 한다.
나아가 본 발명은 치아의 CT 영상을 이용하여 자동으로 상기 치아의 파노라마 영상을 생성하는 시스템에 있어서, 상기 CT 영상에서 치아 영역을 추출하는 영역 추출부, 상기 치아 영역에 포함되는 치열 궁 곡선을 생성하고, 상기 치아 영역과 상기 치열 궁 곡선 사이의 최대 거리를 이웃 화소 범위로 설정하고, 상기 치열 궁 곡선에 포함된 각 점에서 상기 이웃 화소 범위만큼 이격된 점으로 구성된 폐곡선을 생성하는 곡선 생성부, 상기 CT 영상에 포함된 제1 슬라이드에 폐곡선과 상응하는 위치의 폐곡선 영역을 추출하고, 상기 폐곡선 영역에서 상기 치열 궁 곡선을 축으로 휘어짐을 보정하여 직사각 영역을 생성하고, 상기 직사각 영역에 포함된 픽셀에서 동일한 x 값에 대응되는 대표 픽셀을 연산하고, 각 x 값에 대응되는 상기 대표 픽셀을 조합하여 상기 제1 슬라이드의 픽셀 라인을 생성하고, 상기 CT 영상에 포함된 슬라이드 각각에 대하여 반복되는 픽셀 라인 생성부, 각 슬라이드의 픽셀 라인을 합성하여 상기 파노라마 영상을 생성하는 영상 생성부를 포함하는 것을 일 특징으로 한다.
또한 상기 영역 추출부는 상기 CT 영상에 포함된 슬라이드 각각에 대하여 얼굴의 최하단을 기준으로 기 설정된 크기의 후보 영역을 추출하고, 상기 후보 영역에 포함된 픽셀의 밝기 값을 이용하여 경조직을 추출하고, 상기 슬라이드에서 추출된 상기 경조직에 대하여 동일한 위치에 포함된 픽셀의 밝기 값에 대한 평균 값을 연산하여 누적 슬라이드를 생성하고, 상기 누적 슬라이드에 포함된 경조직의 영역 중 크기가 가장 큰 영역을 치아 영역으로 설정하는 것을 일 특징으로 한다.
나아가 상기 픽셀 라인 생성부는 상기 직사각 영역을 픽셀 단위로 분할하여 동일한 x 값을 갖는 하나 이상의 픽셀의 평균 값을 연산하고, 상기 평균 값을 상기 대표 픽셀의 밝기 값으로 설정하는 것을 일 특징으로 한다.
전술한 바와 같은 본 발명에 의하면, 치아 CT 영상에서 파노라마 영상을 자동으로 재구성할 수 있다.
또한 본 발명은 치아 CT 영상에 포함된 슬라이드에서 영상에 따라 설정되는 임계 값을 이용하는 Otsu 이진화 알고리즘을 적용하여 치아 영역을 추출할 수 있다.
또한 본 발명은 적합한 이웃 화소 범위를 설정하여 모든 치아가 정확하게 나타날 수 있도록 하는 파노라마 영상을 재구성할 수 있다.
또한 본 발명은 매개변수 곡선을 이용하여 적절한 치열 궁 곡선을 생성할 수 있다.
도1은 본 발명의 일 실시 예에 의한 자동 파노라마 영상 생성 시스템의 구성을 도시한 도면이다.
도2는 본 발명의 일 실시 예에 의한 자동 파노라마 영상 생성 방법을 설명하기 위한 도면이다.
도3은 본 발명의 일 실시 예에 의한 치아 영역을 추출하는 방법을 설명하기 위한 도면이다.
도4는 본 발명의 일 실시 예에 의한 픽셀 라인을 추출하는 방법을 설명하기 위한 도면이다.
도5는 본 발명의 일 실시 예에 의한 부정 교합을 가진 환자의 CT 영상을 도시한 도면이다.
도6은 본 발명의 일 실시 예에 의한 3mm 크기를 갖는 이웃 화소 범위를 갖는 파노라마 영상을 도시한 도면이다.
도7은 본 발명의 일 실시 예에 의한 20mm 크기를 갖는 이웃 화소 범위를 갖는 파노라마 영상을 도시한 도면이다.
도8은 본 발명의 일 실시 예에 의한 경조직을 추출하는 방법을 설명하기 위한 도면이다.
도9는 본 발명의 일 실시 예에 의한 누적 영상을 도시한 도면이다.
도10은 본 발명의 일 실시 예에 의한 누적 영상에서 치아 영역을 추출하는 방법을 설명하기 위한 도면이다.
도11은 본 발명의 일 실시 예에 의한 치열 궁 곡선을 도시한 도면이다.
도12는 본 발명의 일 실시 예에 의한 폐곡선을 생성하는 방법을 설명하기 위한 도면이다.
도13은 본 발명의 일 실시 예에 의한 픽셀 라인을 추출하는 방법을 설명하기 위한 도면이다.
도14는 본 발명의 일 실시 예에 의한 자동 생성된 파노라마 영상을 도시한 도면이다.
도15는 본 발명의 일 실시 예에 의한 수동 생성된 파노라마 영상을 도시한 도면이다.
전술한 목적, 특징 및 장점은 첨부된 도면을 참조하여 상세하게 후술되며, 이에 따라 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자가 본 발명의 기술적 사상을 용이하게 실시할 수 있을 것이다. 본 발명을 설명함에 있어서 본 발명과 관련된 공지 기술에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우에는 상세한 설명을 생략한다.
도면에서 동일한 참조부호는 동일 또는 유사한 구성요소를 가리키는 것으로 사용되며, 명세서 및 특허청구의 범위에 기재된 모든 조합은 임의의 방식으로 조합될 수 있다. 그리고 다른 식으로 규정하지 않는 한, 단수에 대한 언급은 하나 이상을 포함할 수 있고, 단수 표현에 대한 언급은 또한 복수 표현을 포함할 수 있음이 이해되어야 한다.
본 명세서에서 사용되는 용어는 단지 특정 예시적 실시 예들을 설명할 목적을 가지고 있으며 한정할 의도로 사용되는 것이 아니다. 본 명세서에서 사용된 바와 같은 단수적 표현들은 또한, 해당 문장에서 명확하게 달리 표시하지 않는 한, 복수의 의미를 포함하도록 의도될 수 있다. 용어 "및/또는," "그리고/또는"은 그 관련되어 나열되는 항목들의 모든 조합들 및 어느 하나를 포함한다. 용어 "포함한다", "포함하는", "포함하고 있는", "구비하는", "갖는", "가지고 있는" 등은 내포적 의미를 갖는 바, 이에 따라 이러한 용어들은 그 기재된 특징, 정수, 단계, 동작, 요소, 및/또는 컴포넌트를 특정하며, 하나 이상의 다른 특징, 정수, 단계, 동작, 요소, 컴포넌트, 및/또는 이들의 그룹의 존재 혹은 추가를 배제하지 않는다. 본 명세서에서 설명되는 방법의 단계들, 프로세스들, 동작들은, 구체적으로 그 수행 순서가 확정되는 경우가 아니라면, 이들의 수행을 논의된 혹은 예시된 그러한 특정 순서로 반드시 해야 하는 것으로 해석돼서는 안 된다. 추가적인 혹은 대안적인 단계들이 사용될 수 있음을 또한 이해해야 한다.
또한, 각각의 구성요소는 각각 하드웨어 프로세서로 구현될 수 있고, 위 구성요소들이 통합되어 하나의 하드웨어 프로세서로 구현될 수 있으며, 또는 위 구성요소들이 서로 조합되어 복수 개의 하드웨어 프로세서로 구현될 수도 있다.
이하, 첨부된 도면을 참조하여 본 발명에 따른 바람직한 실시 예를 상세히 설명하기로 한다.
도1은 본 발명의 일 실시 예에 의한 자동 파노라마 영상 생성 시스템의 구성을 도시한 도면이다.
도1을 참조하면 자동 파노라마 영상 생성 시스템은 치아 CT 영상을 이용하여 자동으로 치아의 파노라마 영상을 생성할 수 있다.
치아 CT 영상을 생성하는 장치는 보편적으로 DICOM(Digital Imaging and Communications in Medicine, 의료용 디지털 영상 및 통신) 형식을 사용하여 CT 영상을 생성하고 제공할 수 있다. DICOM은 의료용 기기에서 디지털 영상 표현 또는 통신에 사용되는 여러 가지 표준을 의미한다. DICOM 형식은 RCS(Reference Coordinates System, 기준 좌표계)을 사용할 수 있다. RCS를 이용하면 치아 CT 영상의 각 슬라이드에서 환자의 오른쪽 측면부터 왼쪽 측면으로 향하는 축을 x축, 전면부터 후면으로 향하는 축을 y축, 발부터 머리로 향하는 축을 z축으로 정의할 수 있다. DICOM 형식의 헤더인 (0020, 0032)를 참조하면, 치아 CT 영상 좌측 상단에 기재된 좌표 값(x, y, z) 값은 mm 단위로 제공될 수 있다. 예를 들어, z 값이 작으면 발에 가까운 것을 의미하고, z 값이 크면 머리에 가까운 것을 의미한다.
자동 파노라마 영상 생성 시스템은 영역 추출부(100), 곡선 생성부(400), 픽셀 라인 추출부(200), 영상 생성부(300)를 포함할 수 있다.
영역 추출부(100)는 치아 CT 영상의 각 슬라이드에서 치아 영역을 추출할 수 있다. 영역 추출부(100)는 치아 CT 영상의 슬라이드에서 RCS 좌표의 z값이 가장 작은 슬라이드를 기준으로 치아 영역을 추출할 수 있다. 보다 구체적으로 영역 추출부(100)는 치아를 모두 포함할 수 있도록 하는 범위를 설정하여 얼굴의 최하단을 기준으로 해당 범위의 크기를 갖는 후보 영역을 추출할 수 있다. 일반적으로 사람의 치아는 턱 부분부터 70 mm 이내에 포함되어 있기 때문에 영역 추출부(100)는 치아를 포함할 수 있도록 하는 범위를 70mm 내외로 설정할 수 있다.
영역 추출부(100)는 치아를 포함하는 후보 영역을 추출하면, 치아를 포함하는 경조직을 추출할 수 있다. 치아 CT 영상은 공기질, 연조직 그리고 경조직을 포함할 수 있기 때문에 치아에 대한 보다 정확한 정보를 얻기 위해서는 경조직만을 추출할 필요가 있다.
경조직을 추출하기 위해 영상 추출부(100)는 CT 값을 이용할 수 있다. CT 값은 CT 영상에 포함된 화소의 흡수 정도를 의미하는 단위로, 0의 값을 갖는 물을 중심으로 공기는 -500, 뼈는 +500으로 설정되어 있다. CT 값은 장치의 종류나 촬영 조건에 따라 달라질 수 있다. 따라서 영상 추출부(100)는 CT 값이 음수인 공기질을 제거하기 위하여 CT 값에 대한 임계 값을 0으로 설정하여 CT 값이 0 이상인 부분을 추출할 수 있다.
CT 값이 0 이상인 부분을 추출한 후, 영상 추출부(100)는 연조직과 경조직을 분리하기 위하여 이진화 알고리즘을 적용할 수 있다. 본 발명에서 사용되는 이진화 알고리즘은 Otsu 이진화 알고리즘을 사용할 수 있으나 그 외의 이진화 알고리즘을 사용할 수 있다. Otsu 이진화 알고리즘은 영상의 CT 값에 대한 히스토그램의 형태가 쌍봉형이라고 가정하였을 때 그 문턱값을 임계 값으로 설정하여 영상을 이진화하는 알고리즘이다. 즉 Otsu 이진화 알고리즘은 모든 영역에 대한 분산을 최대로 설정하여 적절한 임계 값을 자동으로 설정할 수 있다. 도8을 참조하면, 본 발명은 Otsu 이진화 알고리즘을 적용하여 경조직에 대한 밝기 값은 255로, 연조직에 대한 밝기 값은 0의 값으로 설정하여 경조직을 분리한 슬라이드를 생성할 수 있다. 영상 추출부(100)는 분리된 경조직의 영역을 이용하여 치아 영역을 추출할 수 있다.
도5를 참조하면, 기존의 치아 영역을 추출하여 파노라마 영상을 생성하면 부정 교합과 같은 치아를 가진 환자의 경우 모든 치아를 포함하기 어렵다는 어려움이 있다. 따라서 CT 영상에 도시된 모든 치아를 포함하는 파노라마 영상을 생성하기 위하여, 영역 추출부(100)는 CT 영상에 포함된 모든 슬라이드에서 추출된 경조직의 영역을 누적할 수 있다. 영역 추출부(100)는 누적된 경조직에 대하여 동일한 위치에 포함된 픽셀의 밝기 값에 대한 평균 값을 연산하여 누적 슬라이드를 생성할 수 있다. 생성된 누적 슬라이드는 도9에 도시된 바와 같은 형상을 가질 수 있으며 경조직의 범위를 정확하게 인식할 수 있도록 한다.
그러나 누적 슬라이드에는 치아 뿐 아니라 뼈와 같은 경조직도 함께 추출되기 때문에 치아 영역만을 추출할 필요가 있다. 영역 추출부(100)는 치아 영역을 추출하기 위하여 Otsu 이진화 알고리즘을 이용하여 경조직과 기타 영역을 뚜렷하게 구분할 수 있다. 영역 추출부(100)는 구분된 경조직에서 치아 영역을 추출하기 위하여 추출된 각 경조직의 영역에 대한 크기를 연산하여 가장 큰 크기를 갖는 영역을 치아 영역으로 설정하고 그 외의 영역을 모두 제거할 수 있다.
영역 추출부(100)가 누적 슬라이드에서 치아 영역을 추출하는 과정은 도10을 참조할 수 있다. 도10을 참조하면 영역 추출부(100)는 누적 슬라이드(1001)에 Otsu 알고리즘을 적용하여 경조직을 뚜렷하게 구분하고(1002), 구분된 경조직 중 크기가 가장 큰 영역을 치아 영역(1003)으로 설정할 수 있다.
이와 같이 치아 영역을 생성함으로써 본원 발명은 기존의 방법으로 추출된 치아 영역보다 치아에 대한 정보를 풍부하게 가질 수 있다는 장점이 있다.
곡선 생성부(400)는 추출된 치아 영역에 포함되는 치열 궁 곡선을 생성하고, 치아 영역과 치열 궁 곡선 사이 사이의 최대 거리를 이웃 화소 범위로 설정하고, 치열 궁 곡선에 포함된 각 점에서 이웃 화소 범위만큼 이격된 점으로 구성된 폐곡선을 생성할 수 있다. 보다 구체적으로 곡선 생성부(400)는 치열 궁 곡선 생성부(410), 이웃 화소 범위 생성부(420) 그리고 폐곡선 생성부(470)를 포함할 수 있다.
치아 영상을 이용하여 파노라마 영상을 생성하기 위해서는 치열 궁 곡선을 필요로 한다. 따라서 치열 궁 곡선 생성부(410)는 치아 영역의 중앙을 관통하는 4차선 곡선인 치열 궁 곡선을 생성할 수 있다. 치열 궁 곡선은 수학식 1과 같이 표현될 수 있다.
Figure 112018022491382-pat00001
수학식 1의
Figure 112018022491382-pat00002
는 4차식 곡선의 계수를 의미한다. 누적 슬라이드에서 검출된 치아 영역을 이용하여 4차식 곡선의 계수를 연산하기 위하여 치열 궁 곡선 생성부(410)는 치아 영역에 포함된 픽셀의 위치를
Figure 112018022491382-pat00003
로 정의하고 수학식 2와 같은 행렬식을 생성할 수 있다.
Figure 112018022491382-pat00004
수학식 2의 N은 치아 영역에 포함된 픽셀의 수를 의미한다. 치열 궁 곡선 생성부(410)는 수학식 2를 연산하기 위해 최소 자승법을 사용할 수 있다. 최소 자승법은 근사 값과 실제 값의 오차의 제곱의 합이 최소가 되는 근사 값을 검출하여 함수을 추정하는 방법으로, 오류의 발생이 적다는 장점이 있다. 치열 궁 곡선 생성부(410)는 최소 자승법의 연산에 있어서 대수적 방법을 사용할 수 있다. 대수적 방법은 문제를 행렬식 형태로 표현한 후 선형대수학을 적용하는 방법이다. 수학식 2는
Figure 112018022491382-pat00005
의 형태를 갖는 행렬식으로,
Figure 112018022491382-pat00006
는 역행렬을 가질 수 없기 때문에 의사 역행렬(Pseudo Inverse)를 이용하여 연산할 수 있다. 치열 궁 곡선 생성부(410)는 의사 역행렬을 이용하여 수학식 2의 행렬식을
Figure 112018022491382-pat00007
와 같이 연산할 수 있다.
치열 궁 곡선 생성부(410)는 연산된 4차식 곡선의 계수를 이용하여 4차식 곡선을 매개변수 곡선의 형태를 갖도록 하기 위해 치아 영역에 포함되는 픽셀 중 가장 작은 x값인
Figure 112018022491382-pat00008
과 가장 큰 x값인
Figure 112018022491382-pat00009
를 연산할 수 있다. 치열 궁 곡선 생성부(410)는
Figure 112018022491382-pat00010
Figure 112018022491382-pat00011
를 이용하여 매개변수 곡선을 수학식 3 내지 수학식 4와 같이 생성할 수 있다.
Figure 112018022491382-pat00012
Figure 112018022491382-pat00013
수학식 3 내지 수학식 4에서
Figure 112018022491382-pat00014
Figure 112018022491382-pat00015
는 치열 궁의 위치를 나타내는 매개변수 곡선이고,
Figure 112018022491382-pat00016
의 범위를 가질 수 있다. 즉 치열 궁 곡선 생성부(410)는
Figure 112018022491382-pat00017
Figure 112018022491382-pat00018
를 치열 궁 곡선으로 설정할 수 있다. 도11에 도시된 예시와 같이 치열 궁 곡선 생성부(410)는 4차식 곡선의 계수를 이용하여
Figure 112018022491382-pat00019
Figure 112018022491382-pat00020
를 치아 영역에 도시할 수 있다. 도11에 도시된 치열 궁 곡선은 치아 영역 내부에 적절하게 위치하고 있는 것을 확인할 수 있다. 본 발명의 치열 궁 곡선은 모든 치아 정보가 포함된 누적 슬라이드에서 추출되었기 때문에 모든 치아 정보를 포함할 수 있다.
이웃 화소 범위 생성부(420)는 모든 치아를 포함하는 파노라마 영상을 생성하기 위해 적절한 이웃 화소 범위를 연산할 수 있다. 이웃 화소의 값을 너무 작게 설정하면 치아가 뚜렷하게 나타나는 반면 치아의 구조를 충분히 포함할 수 없다. 예를 들어 도6에 도시된 바와 같이 이웃 화소 범위를 3 mm의 크기로 설정하였을 경우 치아는 뚜렷하게 확인할 수 있으나 누락된 치아가 존재한다. 또한 이웃 화소의 값을 너무 크게 설정하면 치아는 뚜렷하지 않게 나타나는 반면 치아의 구조를 충분히 포함할 수 있다. 예를 들어 도7에 도시된 바와 같이 이웃 화소 범위를 20 mm의 크기로 설정하였을 경우 모든 치아를 확인할 수 있으나 명확한 구조를 확인할 수 없다.
보다 구체적으로 이웃 화소 범위 생성부(420)는 누적 슬라이드를 이용하여 이웃 화소 범위를 예측할 수 있다. 이웃 화소 범위 생성부(420)는 치열 궁 곡선
Figure 112018022491382-pat00021
Figure 112018022491382-pat00022
를 지나는 N개의 점인
Figure 112018022491382-pat00023
,
Figure 112018022491382-pat00024
를 수학식 5와 같이 생성할 수 있다.
Figure 112018022491382-pat00025
이웃 화소 범위 생성부(420)는 이웃 화소 범위를 설정하기 위해 각 점에 해당되는 치열 궁 곡선의 수직 벡터
Figure 112018022491382-pat00026
를 이용할 수 있다. 이웃 화소 범위 생성부(420)는 수직 벡터가 치아 영역을 벗어나지 않도록 하는 길이를 연산할 수 있다. 예를 들어 도12의 1201은 N이 100일 경우 수직 벡터의 길이를 도시한 도면으로, 이웃 화소 범위 생성부(420)는 생성된 수직 벡터가 치아 영역에서 벗어나지 않도록 1202와 같이 수직 벡터 라인을 설정할 수 있다. 다시 말해서 이웃 화소 범위 생성부(420)는 치열 궁 곡선에 포함된 N개의 점에서 생성된 수직 벡터가 치아 영역을 벗어나지 않도록 하는 길이를 연산할 수 있다.
이웃 화소 범위 생성부(420)는 N개의 점에 대해 연산된 수직 벡터의 길이 중 최대 값을 이웃 화소 범위로 설정할 수 있다.
폐곡선 생성부(430)는 치열 궁 곡선에 포함된 각 점에서 이웃 화소 범위만큼 이격된 점으로 구성된 폐곡선을 생성할 수 있다. 예를 들어 생성된 폐곡선은 도12의 1203과 같은 형태를 가질 수 있다.
픽셀 라인 추출부(200)는 CT 영상에 포함된 각 슬라이드에 폐곡선과 상응하는 위치의 폐곡선 영역을 추출하고, 치열 궁 곡선을 축으로 휘어짐을 보정하여 직사각 영역을 생성하고, 직사각 영역에 포함된 픽셀에서 동일한 x 값을 갖는 하나 이상의 픽셀을 이용하여 x 값에 대응되는 대표 픽셀을 생성하고, 각 x 값에 대응되는 대표 픽셀을 조합하여 각 슬라이드 별로 픽셀 라인을 추출할 수 있다. 치열 궁 곡선과 이웃 화소 범위를 통해 생성된 폐곡선을 이용하여 픽셀 라인 추출부(200)는 보다 정확한 파노라마 영상을 생성할 수 있도록 한다.
보다 구체적으로 픽셀 라인 추출부(200)는 DICOM 형식의 헤더인 (0028, 0030)에서 제공하는 CT 영상에 포함된 픽셀의 실제 크기(mm 단위) 비율
Figure 112018022491382-pat00027
을 이용하여 치열 궁 곡선의 실제 길이(mm 단위)의 근사 값을 연산할 수 있다. 픽셀 라인 추출부(200)가 근사 값을 연산하는 과정은 수학식 6과 같다.
Figure 112018022491382-pat00028
수학식 6의
Figure 112018022491382-pat00029
는 치열 궁 곡선의 mm 단위 길이를 의미하고,
Figure 112018022491382-pat00030
는 픽셀의 가로 길이와 실제 길이의 비율을,
Figure 112018022491382-pat00031
는 픽셀의 세로 길이와 실제 길이의 비율을 의미한다. 픽셀 라인 추출부(200)는
Figure 112018022491382-pat00032
를 이용하여 생성되는 파노라마 영상의 가로 길이를 연산할 수 있다. 보다 구체적으로 픽셀 라인 추출부(200)는 각 픽셀의 실제 길이 비율이
Figure 112018022491382-pat00033
와 같은 값을 갖도록 수학식 7과 같이 파노라마 영상의 가로 길이를 연산할 수 있다.
Figure 112018022491382-pat00034
수학식 7에서
Figure 112018022491382-pat00035
는 파노라마 영상을 생성하기 위해 사용되는 치열 궁 곡선의 샘플의 숫자를 의미한다. 도13을 참조하면, 픽셀 라인 추출부(200)는
Figure 112018022491382-pat00036
개로 나누어진 폐곡선 영역(1301)의 휘어짐을 보정하여 직사각 영역(1302)을 생성할 수 있다.
픽셀 라인 추출부(200)는 직사각 영역을 픽셀 단위로 분할하여 동일한 x 값을 갖는 하나 이상의 픽셀을 이용하여 x 값에 대응되는 대표 픽셀을 생성할 수 있다. 보다 구체적으로 픽셀 라인 추출부(200)는 동일한 x 값을 갖는 하나 이상의 픽셀에 대한 평균 값을 연산하여 대표 픽셀을 생성할 수 있다. 픽셀 라인 추출부(200)는 x 값에 대응되는 대표 픽셀을 조합하여 슬라이드 별로 하나의 픽셀 라인(1303)을 생성할 수 있다. 즉 픽셀 라인은 CT 영상에 포함된 모든 후보 슬라이드에 생성될 수 있다.
영상 생성부(300)는 슬라이드 별로 생성된 픽셀 라인을 합성하여 파노라마 영상을 생성할 수 있다. 영상 생성부(300)는 하나 이상의 픽셀 라인을 합성하여 파노라마 영상을 생성할 수 있다. 이 때, 영상 생성부(300)는 RCS 좌표의 z값이 작은 슬라이드가 아래 쪽으로 가도록 픽셀 라인을 순서대로 합성할 수 있다. 생성된 파노라마 영상은 도14에 도시된 형태를 가질 수 있다. 도14에 도시된 파노라마 영상은 부정 교합을 갖고 있는 환자의 치아 CT 영상을 통해 생성한 것으로, 동일 환자의 치아에 대하여 기존의 방법으로 파노라마 영상을 생성하면 도15와 같은 형태를 가질 수 있다. 도15의 파노라마 영상은 부정 교합이 발생한 치아의 정보를 나타낼 수 없는 반면 본원 발명에 의해 생성된 도14의 파노라마 영상은 모든 치아의 정보를 나타낼 수 있다.
이하에서는 도2 내지 도4를 참조하여 본 발명의 일 실시 예에 의한 자동 파노라마 영상 생성 방법을 설명한다. 자동 파노라마 영상 생성 방법에 관한 설명에서 전술한 자동 파노라마 영상 생성 시스템과 중복되는 세부 실시 예는 생략될 수 있다.
도2는 본 발명의 일 실시 예에 의한 자동 파노라마 영상 생성 방법을 설명하기 위한 도면이다.
도2를 참조하면, 파노라마 영상 생성 시스템은 치아 CT 영상에 포함된 모든 슬라이드에 대하여 Otsu 이진화 알고리즘을 이용하여 경조직과 기타 영역을 구분하여 경조직에 포함된 치아 영역을 추출할 수 있다(S100).
파노라마 영상 생성 시스템은 치아 영역에 대한 파노라마 영상을 생성하기 위해 치아 영역의 중심을 관통하는 치열 궁 곡선을 생성할 수 있다(S200).
파노라마 영상 생성 시스템은 치열 궁 곡선에 포함된 N개의 점에서 치아 영역을 벗어나지 않도록 하는 수직 벡터를 생성할 수 있다. 파노라마 영상 생성 시스템은 생성된 수직 벡터의 길이에 대한 최대 값을 이웃 화소 범위로 설정할 수 있다(S300).
파노라마 영상 생성 시스템은 치열 궁 곡선에 포함된 각 점에서 이웃 화소 범위만큼 이격된 점으로 구성된 폐곡선을 생성할 수 있다(S400).
폐곡선이 생성되면, 파노라마 영상 생성 시스템은 치아 CT 영상에 포함된 슬라이드의 실제 크기 비율을 이용하여 치열 궁 곡선의 근사 값을 연산하고, 생성된 폐곡선에 대하여 치열 궁 곡선을 축으로 하도록 휘어짐을 보정할 수 있다. 또한, 휘어짐이 보정된 폐곡선을 픽셀 단위로 분할하여 동일한 y축에 존재하는 픽셀의 평균 값을 연산할 수 있다. 파노라마 영상 생성 시스템은 각 슬라이드별로 연산된 평균 값을 이용하여 픽셀 라인을 생성할 수 있다(S500).
파노라마 영상 생성 시스템은 슬라이드 별로 생성된 평균 값을 합성하여 파노라마 영상을 생성할 수 있다(S600).
도3은 본 발명의 일 실시 예에 의한 치아 영역을 추출하는 방법을 설명하기 위한 도면이다.
도3을 참조하면, 파노라마 영상 생성 시스템은 치아 CT 영상의 슬라이드에서 모든 치아를 포함할 수 있도록 하는 범위를 설정하여 후보 슬라이드를 추출할 수 있다(S110).
파노라마 영상 생성 시스템은 치아를 보다 뚜렷하게 나타내기 위해 Otsu 이진화 알고리즘을 이용하여 각 후보 슬라이드를 이진화할 수 있다(S120). 이 후, 파노라마 영상 생성 시스템은 이진화된 각 후보 슬라이드에서 밝기 값이 255로 설정된 경조직을 추출하고(S130), 각 후보 슬라이드의 상응하는 위치에 존재하는 픽셀의 픽셀 값 대한 평균 값을 연산(S140)하여 누적 슬라이드를 생성할 수 있다.
누적 슬라이드는 뚜렷한 형태를 갖지 않기 때문에 파노라마 영상 시스템은 슬라이드에서 치아 또는 뼈를 포함하는 경조직을 보다 명확하게 확인할 수 있도록 하기 위하여 Otsu 이진화 알고리즘을 이용하여 누적 슬라이드를 이진화할 수 있다. 파노라마 영상 생성 시스템은 누적 슬라이드를 이진화하면, 이진화된 누적 슬라이드에 나타난 하나 이상의 경조직의 크기를 연산하여 가장 큰 크기를 갖는 영역을 치아 영역으로 추출할 수 있다(S150).
도4는 본 발명의 일 실시 예에 의한 픽셀 라인을 추출하는 방법을 설명하기 위한 도면이다.
도4를 참조하면, 파노라마 영상 생성 시스템은 치아 CT 영상에 포함된 픽셀의 실제 길이의 근사 값을 연산할 수 있다(S510).
파노라마 영상 생성 시스템은 폐곡선에 포함된 치열 궁 곡선을 축으로 하여 폐곡선의 휘어짐을 보정할 수 있다.
파노라마 영상 생성 시스템은 보정된 폐곡선을 픽셀 단위로 분할하여(S520) 동일한 y 축에 존재하는 픽셀의 평균 값을 연산할 수 있다(S530).
평균 값을 연산한 파노라마 영상 시스템은 누적 슬라이드 별로 하나의 픽셀 라인을 생성할 수 있다(S540).
본 명세서와 도면에 개시된 본 발명의 실시 예들은 본 발명의 기술 내용을 쉽게 설명하고 본 발명의 이해를 돕기 위해 특정 예를 제시한 것뿐이며, 본 발명의 범위를 한정하고자 하는 것은 아니다. 여기에 개시된 실시 예들 이외에도 본 발명의 기술적 사상에 바탕을 둔 다른 변형 예들이 실시 가능하다는 것은 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자에게 자명한 것이다.

Claims (8)

  1. 파노라마 영상 생성 시스템에서 치아의 CT 영상을 이용하여 자동으로 상기 치아의 파노라마 영상을 생성하는 방법에 있어서,
    상기 CT 영상에서 치아 영역을 추출하는 단계;
    상기 치아 영역에 포함되는 치열 궁 곡선을 생성하는 단계;
    상기 치아 영역과 상기 치열 궁 곡선의 최대 거리를 이웃 화소 범위로 설정하는 단계;
    상기 치열 궁 곡선에 포함된 각 점에서 상기 이웃 화소 범위만큼 이격된 점으로 구성된 폐곡선을 생성하는 단계;
    상기 CT 영상에 포함된 제1 슬라이드에 폐곡선과 상응하는 위치의 폐곡선 영역을 추출하는 a 단계;
    상기 폐곡선 영역에서 상기 치열 궁 곡선을 축으로 휘어짐을 보정하여 직사각 영역을 생성하는 b 단계;
    상기 직사각 영역에 포함된 픽셀에서 동일한 x 값을 갖는 하나 이상의 픽셀을 이용하여 상기 x 값에 대응되는 대표 픽셀을 생성하는 c 단계;
    각 x 값에 대응되는 상기 대표 픽셀을 조합하여 상기 제1 슬라이드의 픽셀 라인을 생성하는 d 단계;
    상기 CT 영상에 포함된 슬라이드 각각에 대하여 상기 a 내지 상기 d 단계를 수행하고, 각 슬라이드의 픽셀 라인을 합성하여 상기 파노라마 영상을 생성하는 단계를 포함하는 자동 파노라마 영상 생성 방법.
  2. 제1항에 있어서,
    상기 치아 영역을 추출하는 단계는,
    상기 CT 영상에 포함된 슬라이드 각각에 대하여 얼굴의 최하단을 기준으로 기 설정된 크기의 후보 영역을 추출하는 단계;
    상기 후보 영역에 포함된 픽셀의 밝기 값을 이용하여 경조직을 추출하는 단계;
    상기 슬라이드에서 추출된 상기 경조직에 대하여 동일한 위치에 포함된 픽셀의 밝기 값에 대한 평균 값을 연산하여 누적 슬라이드를 생성하는 단계;
    상기 누적 슬라이드에 포함된 경조직의 영역 중 크기가 가장 큰 영역을 치아 영역으로 설정하는 단계를 포함하는 자동 파노라마 영상 생성 방법.
  3. 제2항에 있어서,
    상기 경조직을 추출하는 단계는
    상기 후보 영역의 밝기 값을 Otsu 알고리즘을 이용하여 이진화하는 것을 특징으로 하는 자동 파노라마 영상 생성 방법.
  4. 제2항에 있어서,
    상기 누적 슬라이드를 생성하는 단계는,
    상기 슬라이드 각각에서 추출된 경조직의 상응하는 위치에 포함된 픽셀의 밝기 값에 대한 평균 값을 연산하여 누적 슬라이드를 생성하는 단계;
    상기 누적 슬라이드를 Otsu 알고리즘을 이용하여 이진화하는 단계를 포함하는 자동 파노라마 영상 생성 방법.
  5. 제1항에 있어서,
    상기 대표 픽셀을 생성하는 단계는
    상기 직사각 영역을 픽셀 단위로 분할하여 동일한 x 값을 갖는 하나 이상의 픽셀의 평균 값을 연산하는 단계;
    상기 평균 값을 상기 대표 픽셀의 밝기 값으로 설정하는 단계를 포함하는 자동 파노라마 영상 생성 방법.
  6. 치아의 CT 영상을 이용하여 자동으로 상기 치아의 파노라마 영상을 생성하는 시스템에 있어서,
    상기 CT 영상에서 치아 영역을 추출하는 영역 추출부;
    상기 치아 영역에 포함되는 치열 궁 곡선을 생성하고, 상기 치아 영역과 상기 치열 궁 곡선 사이의 최대 거리를 이웃 화소 범위로 설정하고, 상기 치열 궁 곡선에 포함된 각 점에서 상기 이웃 화소 범위만큼 이격된 점으로 구성된 폐곡선을 생성하는 곡선 생성부;
    상기 CT 영상에 포함된 제1 슬라이드에 폐곡선과 상응하는 위치의 폐곡선 영역을 추출하고, 상기 폐곡선 영역에서 상기 치열 궁 곡선을 축으로 휘어짐을 보정하여 직사각 영역을 생성하고, 상기 직사각 영역에 포함된 픽셀에서 동일한 x 값에 대응되는 대표 픽셀을 연산하고, 각 x 값에 대응되는 상기 대표 픽셀을 조합하여 상기 제1 슬라이드의 픽셀 라인을 생성하고, 상기 CT 영상에 포함된 슬라이드 각각에 대하여 반복되는 픽셀 라인 생성부;
    각 슬라이드의 픽셀 라인을 합성하여 상기 파노라마 영상을 생성하는 영상 생성부를 포함하는 것을 특징으로 하는 자동 파노라마 영상 생성 시스템.
  7. 제6항에 있어서,
    상기 영역 추출부는 상기 CT 영상에 포함된 슬라이드 각각에 대하여 얼굴의 최하단을 기준으로 기 설정된 크기의 후보 영역을 추출하고, 상기 후보 영역에 포함된 픽셀의 밝기 값을 이용하여 경조직을 추출하고, 상기 슬라이드에서 추출된 상기 경조직에 대하여 동일한 위치에 포함된 픽셀의 밝기 값에 대한 평균 값을 연산하여 누적 슬라이드를 생성하고, 상기 누적 슬라이드에 포함된 경조직의 영역 중 크기가 가장 큰 영역을 치아 영역으로 설정하는 것을 특징으로 하는 자동 파노라마 영상 생성 시스템.
  8. 제6항에 있어서,
    상기 픽셀 라인 생성부는
    상기 직사각 영역을 픽셀 단위로 분할하여 동일한 x 값을 갖는 하나 이상의 픽셀의 평균 값을 연산하고, 상기 평균 값을 상기 대표 픽셀의 밝기 값으로 설정하는 것을 특징으로 하는 자동 파노라마 영상 생성 시스템.
KR1020180026260A 2018-03-06 2018-03-06 치아 ct 영상을 이용한 자동 파노라마 영상 생성 방법 및 시스템 KR102054210B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020180026260A KR102054210B1 (ko) 2018-03-06 2018-03-06 치아 ct 영상을 이용한 자동 파노라마 영상 생성 방법 및 시스템

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020180026260A KR102054210B1 (ko) 2018-03-06 2018-03-06 치아 ct 영상을 이용한 자동 파노라마 영상 생성 방법 및 시스템

Publications (2)

Publication Number Publication Date
KR20190105763A KR20190105763A (ko) 2019-09-18
KR102054210B1 true KR102054210B1 (ko) 2019-12-10

Family

ID=68070647

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020180026260A KR102054210B1 (ko) 2018-03-06 2018-03-06 치아 ct 영상을 이용한 자동 파노라마 영상 생성 방법 및 시스템

Country Status (1)

Country Link
KR (1) KR102054210B1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20220138919A (ko) 2021-04-06 2022-10-14 (주)레비스톤 3차원 구강내 열화상 영상 생성 시스템 및 방법

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102251651B1 (ko) * 2019-03-19 2021-05-14 오스템임플란트 주식회사 파노라믹 영상 생성을 위한 커브 라인 자동 생성방법 및 이를 위한 치과영상 처리장치
KR102302245B1 (ko) 2019-11-07 2021-09-14 부산대학교 산학협력단 딥러닝을 이용한 파노라마 이중상 및 허상 제거를 위한 장치 및 방법
KR102313749B1 (ko) * 2019-12-23 2021-10-18 주식회사 메가젠임플란트 인공지능 기반 자동 구강 ct 색상변환장치 및 그 장치의 구동방법
KR102346199B1 (ko) * 2020-02-05 2022-01-03 오스템임플란트 주식회사 파노라믹 영상 생성 방법 및 이를 위한 영상 처리장치
KR102566311B1 (ko) * 2022-06-17 2023-08-10 김정란 신경치료를 위한 엔도 가이드 생성 방법 및 시스템

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090052617A1 (en) 2007-02-22 2009-02-26 J. Morita Manufacturing Corporation Image processing method, image display method, image processing program, storage medium, image processing apparatus and X-ray imaging apparatus
US20140153799A1 (en) 2011-04-25 2014-06-05 Vatech Co., Ltd. Method and apparatus for obtaining panoramic image
CN105608747A (zh) 2015-12-17 2016-05-25 首都师范大学 一种从三维牙科锥形束ct数据中提取全景图的方法
KR101827354B1 (ko) 2016-09-23 2018-02-08 오스템임플란트 주식회사 치아 영상 처리 방법, 이를 위한 장치, 및 이를 기록한 기록 매체

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101669178B1 (ko) * 2014-12-31 2016-10-27 오스템임플란트 주식회사 치아 파노라마 영상 생성방법, 이를 위한 장치 및 이를 기록한 기록매체

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090052617A1 (en) 2007-02-22 2009-02-26 J. Morita Manufacturing Corporation Image processing method, image display method, image processing program, storage medium, image processing apparatus and X-ray imaging apparatus
US20140153799A1 (en) 2011-04-25 2014-06-05 Vatech Co., Ltd. Method and apparatus for obtaining panoramic image
CN105608747A (zh) 2015-12-17 2016-05-25 首都师范大学 一种从三维牙科锥形束ct数据中提取全景图的方法
KR101827354B1 (ko) 2016-09-23 2018-02-08 오스템임플란트 주식회사 치아 영상 처리 방법, 이를 위한 장치, 및 이를 기록한 기록 매체

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20220138919A (ko) 2021-04-06 2022-10-14 (주)레비스톤 3차원 구강내 열화상 영상 생성 시스템 및 방법

Also Published As

Publication number Publication date
KR20190105763A (ko) 2019-09-18

Similar Documents

Publication Publication Date Title
KR102054210B1 (ko) 치아 ct 영상을 이용한 자동 파노라마 영상 생성 방법 및 시스템
US11464467B2 (en) Automated tooth localization, enumeration, and diagnostic system and method
US11443423B2 (en) System and method for constructing elements of interest (EoI)-focused panoramas of an oral complex
EP2560572B1 (en) Reduction and removal of artifacts from a three-dimensional dental x-ray data set using surface scan information
Kumar et al. Descriptive analysis of dental X-ray images using various practical methods: A review
CN106846307B (zh) 基于锥形束计算机体层摄影的图像处理方法及装置
Hosntalab et al. Segmentation of teeth in CT volumetric dataset by panoramic projection and variational level set
US20150254816A1 (en) Alignment of mixed-modality data sets for reduction and removal of imaging artifacts
KR101105494B1 (ko) 환자 맞춤형 3차원 인체 뼈 모델 재구성 방법
CN111388125B (zh) 一种口腔正畸治疗前后牙齿移动量计算方法和装置
Kim et al. Automatic extraction of inferior alveolar nerve canal using feature-enhancing panoramic volume rendering
US20220012888A1 (en) Methods and system for autonomous volumetric dental image segmentation
US20210217170A1 (en) System and Method for Classifying a Tooth Condition Based on Landmarked Anthropomorphic Measurements.
CN113164153A (zh) 用于生成2d全景图像的方法和装置
Gao et al. Automatic tooth region separation for dental CT images
CN117769391A (zh) 在自然头位上拍摄的三维cbct影像中基于机器学习的用于牙齿矫正诊断的头部计测参数导出方法
KR102302249B1 (ko) 영상처리와 cnn을 이용한 자동 3차원 세팔로메트리 장치 및 방법
JP2022530298A (ja) X線歯科ボリュームトモグラフィにおける金属アーチファクト低減の方法
US20220358740A1 (en) System and Method for Alignment of Volumetric and Surface Scan Images
US20070053570A1 (en) Image processing method, and computer-readable recording medium in which image processing program is recorded
CN107408301B (zh) 使用通道检测对图像数据中的对象的分割
JP2007014435A (ja) 画像処理装置、画像処理方法およびそのプログラム
KR102231334B1 (ko) 치아 배열 분석 방법, 이를 이용한 파노라마 이미지 생성 방법 및 그 장치
JP2004152043A (ja) 差分画像の補正方法および画像処理装置
JP2006175057A (ja) 画像処理装置、そのプログラムおよび方法

Legal Events

Date Code Title Description
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant