KR101870890B1 - 콘빔 ct 영상의 금속으로 인한 아티팩트 보정방법 및 장치 - Google Patents

콘빔 ct 영상의 금속으로 인한 아티팩트 보정방법 및 장치 Download PDF

Info

Publication number
KR101870890B1
KR101870890B1 KR1020160156707A KR20160156707A KR101870890B1 KR 101870890 B1 KR101870890 B1 KR 101870890B1 KR 1020160156707 A KR1020160156707 A KR 1020160156707A KR 20160156707 A KR20160156707 A KR 20160156707A KR 101870890 B1 KR101870890 B1 KR 101870890B1
Authority
KR
South Korea
Prior art keywords
image
value
linear
calculating
linear combination
Prior art date
Application number
KR1020160156707A
Other languages
English (en)
Other versions
KR20180058062A (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 KR1020160156707A priority Critical patent/KR101870890B1/ko
Publication of KR20180058062A publication Critical patent/KR20180058062A/ko
Application granted granted Critical
Publication of KR101870890B1 publication Critical patent/KR101870890B1/ko

Links

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
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
    • A61B6/4085—Cone-beams
    • G06T5/005—
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/77—Retouching; Inpainting; Scratch removal

Landscapes

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

Abstract

본 발명은 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법 및 장치에 관한 것이다. 더 상세하게는 CT(컴퓨터단층촬영) 촬영된 영상의 이미지를 입력받고, 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출하여 금속으로 인한 손상영역을 추출하고, 추출된 손상영역을 선형보간(Linear Interpolation)하고, 보간된 손상영역데이터값을 제1메인영상으로 재구성하고, 재구성된 영상을 하나 또는 둘 이상의 이미지로 분할하고, 분할된 이미지를 기초로 선형조합값을 산출하고, 산출한 선형조합값과 영상이미지데이터값을 기초로 최종보정이미지 출력여부를 결정한다. 이를 통해 콘빔 CT 영상의 금속으로 인한 아티팩트 영역을 반복적으로 보정할 수 있으며, 점점 더 많은 CT 영상의 사용과 저선량 CT의 사용이 늘어나고 있음에 따라 미래 CT 산업 성장의 촉진제 역할을 할 수 있다.

Description

콘빔 CT 영상의 금속으로 인한 아티팩트 보정방법 및 장치 {Artifact correction method and device for metal of cone beam CT}
본 발명은 콘빔 CT 영상의 금속으로 인한 아티팩트 보정 방법 및 장치에 관한 것이다. 보다 상세하게는 여러 개의 프라이어(prior) 영상을 이용하여 콘빔 CT 영상의 금속으로 인한 아티팩트를 보정하는 방법 및 장치에 관한 기술이다.
도 1을 참조하면, CT(Computed Tomography : 컴퓨터 단층 촬영 영상)의 금속으로 인한 아티팩트는 특히 치아 CT(Computed Tomography : 컴퓨터 단층 촬영 영상)와 수술 후 고관절 CT(Computed Tomography : 컴퓨터 단층 촬영 영상)에서 많이 나타난다. 이를 보정하는 알고리즘을 적용하기 위해서는 첫 번째 전처리 단계로 금속 부분을 영상에서 분류해내는 작업이 필요하다. 일반적으로 지금까지는 주로 일정 임계값(threshold) 이상은 금속으로 분류하는 방법을 이용하였다. 하지만 이는 한 가지 종류의 금속일 경우에도 임계값에 따라 경계치를 제대로 추출하기 힘든 한계를 가진다. 뿐만 아니라 여러 개의 금속이 존재하는 경우 심한 금속 아티팩트가 존재한다. 도 1에서 보듯이 여러 개의 금속이 존재하는 경우는 금속 사이에 아티팩트로 인해 밝게 나타나는 부분이 생기므로 적정한 임계값을 찾는 것은 어렵다. 즉 서로 다른 종류의 금속이 존재할 경우 한 개의 임계값으로 복수 개의 금속을 한 번에 분류하기에 어려운 한계가 있다.
도 2는 종래에 정의된 하나의 프라이어(prior) 영상을 이용한 금속 인공물 감약 알고리즘 과정을 보여주는 도면이다. 금속 인공물 감약 알고리즘은 크게 두 가지로 분류할 수 있다. 첫 번째 방법은 보간 또는 근사를 이용한 방법, 두 번째 방법은 반복적 재구성 알고리즘을 이용한 방법이다.
먼저 첫 번째 방법은 금속영상을 앞으로 투사(forward projection)하여 나오는 사이노그램(sinogram)을 왜곡된 부분으로 보고, 사이노그램(sinogram)상에서 보간이나 근사 알고리즘을 사용하여 계산한다. 다만 사이노그램(sinogram)의 원리를 이용하지 않고 단순히 사이노그램(sinogram) 공간에서의 거리만을 고려한 보간 또는 근사 이용 방법은 새로운 인공물을 발생시킨다는 문제가 있다.
두 번째 방법은 반복적인 계산으로 인하여 효율적이지 않은 결과를 보이고 있다.
이에 따라 최근에는 첫 번째 방법과 두 번째 방법의 문제점을 보완하기 위해서 프라이어(prior) 영상이 투사(forward projection)된 사이노그램(sinogram)을 이용하여 왜곡된 부분을 복원하는 방법이 제시되고 있다. 다만 이 방법은 프라이어(prior) 영상 질에 따라 성능이 크게 영향을 받으므로, 좋은 프라이어(prior)를 찾는 것이 가장 결정적인 요인으로 작용한다. 현재까지 발표된 논문은 도 2에 제시된 바와 같이 임계값(threshold) 기반으로 정의된 프라이어(prior) 영상이다.
한국등록특허공보 제1591381호 한국공개특허공보 제2015-0095140호 미국등록특허공보 제8768027호
이에 본 발명은 상기와 같은 제반 사항을 고려하여 제안된 것으로, 현재 사용되고 있는 프라이어(prior) 영상을 기반으로 하면서 하나의 프라이어(prior)가 아닌 각 조직의 영상 값에 따른 여러 개의 서로 다른 프라이어(prior)를 이용하여 반복적으로 금속아티팩트 영역을 보정하는 콘빔 CT 영상의 금속으로 인한 아티팩트 보정 방법 및 장치를 제공하는 것을 목적으로 한다.
상기와 같은 목적을 달성하기 위하여 본 발명의 기술적 사상에 의한 콘빔 CT 영상의 금속아티팩트 영역 보정방법은, CT(컴퓨터단층촬영) 촬영된 영상의 이미지를 입력받는 영상이미지수신단계, 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출하여 금속으로 인한 손상영역을 추출하는 손상영역추출단계, 추출된 손상영역을 선형보간(Linear Interpolation)하는 제1선형보간단계, 보간된 손상영역데이터값을 제1메인영상으로 재구성하는 제1메인영상재구성단계, 재구성된 영상을 하나 또는 둘 이상의 이미지로 분할하는 이미지분할단계, 분할된 이미지를 기초로 선형조합값을 산출하는 선형조합값산출단계, 산출한 선형조합값과 상기 영상이미지데이터값을 기초로 최종보정이미지 출력여부를 결정하는 최종보정이미지출력결정단계를 포함하는 것을 특징으로 한다.
손상영역추출단계는 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출하는 영상이미지데이터값산출단계, 산출된 영상이미지데이터값을 서브영상으로 재구성하는 서브영상재구성단계, 재구성된 서브영상으로부터 금속으로 인한 손상영역의 데이터를 추출하는 손상영역데이터추출단계를 포함하는 것을 특징으로 한다.
선형조합값산출단계는 분할된 하나 또는 둘 이상의 이미지로부터 제1사이노그램(sinogram)값을 산출하는 제1사이노그램(sinogram)값산출단계, 산출된 하나 또는 둘 이상의 제1사이노그램(sinogram)값을 기저값으로 하여 비손상영역에 근사하도록 선형계수값을 계산하는 선형계수값계산단계, 계산된 선형계수값을 이용하여 선형조합값을 계산하는 선형조합값계산단계를 포함하는 것을 특징으로 한다. 이 때 선형계수값계산단계는 수학식 1로 정의되는 선형계수 계산식을 이용하여
Figure 112016114771226-pat00001
이 최대한 영상이미지데이터값(
Figure 112016114771226-pat00002
)과 근사한 선형계수값(
Figure 112016114771226-pat00003
)을 계산하는 것을 특징으로 한다.
[수학식 1]
Figure 112016114771226-pat00004
선형조합값계산단계는 수학식 2로 정의되는 선형조합 계산식을 이용하여 각각의 선형계수값(
Figure 112016114771226-pat00005
)와 기저값(
Figure 112016114771226-pat00006
) 곱의 합으로 계산하는 것을 특징으로 한다.
[수학식 2]
Figure 112016114771226-pat00007
최종보정이미지출력결정단계는 계산된 선형조합값에서 영상이미지데이터값을 뺀 값인 잔여오류값을 계산하는 잔여오류값계산단계, 계산된 잔여오류값의 손상영역을 선형보간(Linear Interpolation)하는 제2선형보간단계, 보간된 손상영역데이터값과 선형조합값을 더한 값인 보정사이노그램(sinogram)값을 계산하는 보정사이노그램(sinogram)값계산단계, 계산된 잔여오류값과 기설정한 임계값의 비교를 통해 최종보정이미지 출력여부를 판단하는 최종보정이미지출력여부판단단계를 포함하는 것을 특징으로 한다.
이 때 최종보정이미지출력여부판단단계는 잔여오류값이 기설정한 임계값보다 클 경우 제1선형보간단계 내지 최종보정이미지출력결정단계를 실행하고, 잔여오류값이 기설정한 임계값보다 작을 경우 상기 계산된 보정사이노그램(sinogram)값을 제2메인영상으로 재구성하고, 이를 출력하는 것을 특징으로 한다.
상기와 같은 목적을 달성하기 위하여 본 발명의 기술적 사상에 의한 콘빔 CT 영상의 금속아티팩트 영역 보정장치는, CT(컴퓨터단층촬영) 촬영된 영상의 이미지를 입력받는 영상이미지수신부, 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출하여 금속으로 인한 손상영역을 추출하는 손상영역추출부, 추출된 손상영역을 선형보간(Linear Interpolation)하는 제1선형보간부, 보간된 손상영역데이터값을 제1메인영상으로 재구성하는 제1메인영상재구성부, 재구성된 영상을 하나 또는 둘 이상의 이미지로 분할하는 이미지분할부, 분할된 이미지를 기초로 선형조합값을 산출하는 선형조합값산출부, 산출한 선형조합값과 상기 영상이미지데이터값을 기초로 최종보정이미지 출력여부를 결정하는 최종보정이미지출력결정부를 포함하는 것을 특징으로 한다.
손상영역추출부는 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출하는 영상이미지데이터값산출부, 산출된 영상이미지데이터값을 서브영상으로 재구성하는 서브영상재구성부, 재구성된 서브영상으로부터 금속으로 인한 손상영역의 데이터를 추출하는 손상영역데이터추출부를 포함하는 것을 특징으로 한다.
선형조합값산출부는 분할된 하나 또는 둘 이상의 이미지로부터 제1사이노그램(sinogram)값을 산출하는 제1사이노그램(sinogram)값산출부, 산출된 하나 또는 둘 이상의 제1사이노그램(sinogram)값을 기저값으로 하여 비손상영역에 근사하도록 선형계수값을 계산하는 선형계수값계산부, 계산된 선형계수값을 이용하여 선형조합값을 계산하는 선형조합값계산부를 포함하는 것을 특징으로 한다. 이 때 선형계수값계산부는 수학식 1로 정의되는 선형계수 계산식을 이용하여
Figure 112016114771226-pat00008
이 최대한 영상이미지데이터값(
Figure 112016114771226-pat00009
)과 근사한 선형계수값(
Figure 112016114771226-pat00010
)을 계산하는 것을 특징으로 한다.
[수학식 1]
Figure 112016114771226-pat00011
선형조합값계산부는 수학식 2로 정의되는 선형조합 계산식을 이용하여 각각의 선형계수값(
Figure 112016114771226-pat00012
)와 기저값(
Figure 112016114771226-pat00013
) 곱의 합으로 계산하는 것을 특징으로 한다.
[수학식 2]
Figure 112016114771226-pat00014
최종보정이미지출력결정부는 계산된 선형조합값에서 영상이미지데이터값을 뺀 값인 잔여오류값을 계산하는 잔여오류값계산부, 계산된 잔여오류값의 손상영역을 선형보간(Linear Interpolation)하는 제2선형보간부, 보간된 손상영역데이터값과 선형조합값을 더한 값인 보정사이노그램(sinogram)값을 계산하는 보정사이노그램(sinogram)값계산부, 계산된 잔여오류값과 기설정한 임계값의 비교를 통해 최종보정이미지 출력여부를 판단하는 최종보정이미지출력여부판단부를 포함하는 것을 특징으로 한다.
이 때 최종보정이미지출력여부판단부는 잔여오류값이 기설정한 임계값보다 클 경우 제1선형보간부 내지 최종보정이미지출력결정부를 실행하고, 잔여오류값이 기설정한 임계값보다 작을 경우 상기 계산된 보정사이노그램(sinogram)값을 제2메인영상으로 재구성하고, 이를 출력하는 것을 특징으로 한다.
이상에서 설명한 바와 같은 콘빔 CT 영상의 금속으로 인한 아티팩트 보정 방법 및 장치에 따르면,
첫째, 현재 사용되고 있는 프라이어(prior) 영상을 기반으로 하여 콘빔 CT 영상의 금속 아티팩트 영역을 반복적으로 보정할 수 있는 효과를 가진다.
둘째, 점점 더 많은 CT 영상의 사용과 저선량 CT의 사용이 늘어나고 있음에 따라 미래 CT 산업 성장의 촉진제가 될 수 있는 효과를 가진다.
도 1은 여러 개 존재하는 금속 인공물로 인해 밝게 나타나는 CT 영상을 보여주는 도면.
도 2는 종래에 정의된 하나의 프라이어(prior) 영상을 이용한 금속 인공물 감약 알고리즘 과정을 보여주는 도면.
도 3은 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법을 나타낸 순서도.
도 4는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치를 나타낸 구성도.
도 5는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치에서 손상영역추출부를 나타낸 구성도.
도 6은 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치에서 선형조합값산출부를 나타낸 구성도.
도 7은 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치에서 최종보정이미지출력결정부를 나타낸 구성도.
도 8은 본 발명의 일 실시예에 따른 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정 과정을 나타낸 도면.
도 9는 본 발명의 이미지분할과정에서 적용 가능한 요소기술 중 하나인 재귀적 액티브 컨투어(RAC) 알고리즘을 나타낸 도면.
도 10은 본 발명의 이미지분할과정에서 적용 가능한 요소기술 중 하나인 재귀적 액티브 컨투어(RAC) 알고리즘을 적용한 결과를 나타낸 도면.
도 11은 본 발명의 일 실시예에 따른 콘빔 CT 영상의 금속 아티팩트 영역 보정 결과를 나타낸 도면.
첨부한 도면을 참조하여 본 발명의 실시예들에 의한 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법 및 장치에 대하여 상세히 설명한다. 본 발명은 다양한 변경을 가할 수 있고 여러 가지 형태를 가질 수 있는바, 특정 실시예들을 도면에 예시하고 본문에 상세하게 설명하고자 한다. 그러나 이는 본 발명을 특정한 개시 형태에 대해 한정하려는 것이 아니며, 본 발명의 사상 및 기술범위에 포함되는 모든 영상을 포함하는 것으로 이해되어야 한다. 각 도면을 설명하면서 유사한 참조부호를 유사한 구성요소에 대해 사용하였다. 첨부된 도면에 있어서, 구조물들의 치수는 본 발명의 명확성을 기하기 위하여 실제보다 확대하거나, 개략적인 구성을 이해하기 위하여 실제보다 축소하여 도시한 것이다.
또한, 제1 및 제2 등의 용어는 다양한 구성요소들을 설명하는데 사용될 수 있지만, 상기 구성요소들은 상기 용어들에 의해 한정되어서는 안 된다. 상기 용어들은 하나의 구성요소를 다른 구성요소로부터 구별하는 목적으로만 사용된다. 예를 들어, 본 발명의 권리 범위를 벗어나지 않으면서 제1 구성요소는 제2 구성요소로 명명될 수 있고, 유사하게 제2 구성요소도 제1 구성요소로 명명될 수 있다. 한편, 다르게 정의되지 않는 한, 기술적이거나 과학적인 용어를 포함해서 여기서 사용되는 모든 용어들은 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자에 의해 일반적으로 이해되는 것과 동일한 의미를 가지고 있다. 일반적으로 사용되는 사전에 정의되어 있는 것과 같은 용어들은 관련 기술의 문맥 상 가지는 의미와 일치하는 의미를 가지는 것으로 해석되어야 하며, 본 출원에서 명백하게 정의하지 않는 한, 이상적이거나 과도하게 형식적인 의미로 해석되지 않는다.
도 3을 참조하면, 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법은 영상이미지수신단계(S100), 영상이미지데이터값산출단계(S210), 서브영상재구성단계(S230), 손상영역데이터추출단계(S250), 제1선형보간단계(S300), 제1메인영상재구성단계(S400), 이미지분할단계(S500), 제1사이노그램(sinogram)값산출단계(S610), 선형계수값계산단계(S630), 선형조합값계산단계(S650), 잔여오류값계산단계(S710), 제2선형보간단계(S730), 보정사이노그램(sinogram)값계산단계(S750), 최종보정이미지출력여부판단단계(S770)를 포함하여 구성된다.
영상이미지수신단계(S100)는 CT 촬영된 영상의 이미지를 입력받는다.
영상이미지데이터값산출단계(S210)는 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출한다.
서브영상재구성단계(S230)는 산출한 영상이미지데이터값을 서브영상으로 재구성한다.
손상영역데이터추출단계(S250)는 재구성한 서브영상으로부터 금속으로 인한 손상영역의 데이터를 추출한다.
제1선형보간단계(S300)는 추출된 손상영역을 선형보간한다.
제1메인영상재구성단계(S400)는 보간된 손상영역데이터값을 제1메인영상으로 재구성한다.
이미지분할단계(S500)는 재구성된 영상을 하나 또는 둘 이상의 이미지로 분할한다.
제1사이노그램(sinogram)값산출단계(S610)는 분할된 이미지로부터 1사이노그램(sinogram)값을 산출한다.
선형계수값계산단계(S630)는 사이노그램(sinogram)값을 기저값으로 하여 비손상영역에 근사하도록 선형계수값을 계산한다. 이 때 선형계수값계산단계(S630)는 수학식 1로 정의되는 선형계수 계산식을 이용하여
Figure 112016114771226-pat00015
이 최대한 영상이미지데이터값(
Figure 112016114771226-pat00016
)과 근사한 선형계수값(
Figure 112016114771226-pat00017
)을 계산한다.
[수학식 1]
Figure 112016114771226-pat00018
선형조합값계산단계(S650)는 계산한 선형계수값을 이용하여 선형조합값을 계산한다. 이 때 선형조합값계산단계(S650)는 수학식 2로 정의되는 선형조합 계산식을 이용하여 각각의 선형계수값(
Figure 112016114771226-pat00019
)와 기저값(
Figure 112016114771226-pat00020
) 곱의 합으로 계산한다.
[수학식 2]
Figure 112016114771226-pat00021
잔여오류값계산단계(S710)는 계산한 선형조합값에서 영상이미지데이터값을 뺀 값인 잔여오류값을 계산한다.
제2선형보간단계(S730)는 계산한 잔여오류값의 손상영역을 선형보간한다.
보정사이노그램(sinogram)값계산단계(S750)는 손상영역데이터값과 선형조합값을 더한 값인 보정사이노그램(sinogram)값을 계산한다.
최종보정이미지출력여부판단단계(S770)는 계산된 잔여오류값과 기설정한 임계값의 비교를 통해 최종보정이미지 출력여부를 판단한다.
잔여오류값이 기설정한 임계값보다 클 경우 제1선형보간단계(S300) 내지 최종보정이미지출력여부판단단계(S770)를 실행하고, 잔여오류값이 기설정한 임계값보다 작을 경우 보정사이노그램(sinogram)값을 제2메인영상으로 재구성하고, 이를 출력하는 것을 특징으로 한다.
도 4를 참조하면, 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치는 영상이미지수신부(100), 손상영역추출부(200), 제1선형보간부(300), 제1메인영상재구성부(400), 이미지분할부(500), 선형조합값산출부(600), 최종보정이미지출력결정부(700)를 포함하여 구성된다.
영상이미지수신부(100)는 CT 촬영된 영상의 이미지를 입력받는다.
손상영역추출부(200)는 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출하여 금속으로 인한 손상영역을 추출한다.
제1선형보간부(300)는 추출된 손상영역을 선형보간한다.
제1메인영상재구성부(400)는 선형보간된 손상영역데이터값을 제1메인영상으로 재구성한다.
이미지분할부(500)는 재구성된 영상을 하나 또는 둘 이상의 이미지로 분할한다.
선형조합값산출부(600)는 분할된 이미지를 기초로 선형조합값을 산출한다.
최종보정이미지출력결정부(700)는 산출한 선형조합값과 영상이미지데이터값을 기초로 최종보정이미지 출력여부를 결정한다.
도 5를 참조하면, 손상영역추출부(200)는 영상이미지데이터값산출부(210), 서브영상재구성부(230), 손상영역데이터추출부(250)를 포함하여 구성된다.
영상이미지데이터값산출부(210)는 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출한다.
서브영상재구성부(230)는 산출한 영상이미지데이터값을 서브영상으로 재구성한다.
손상영역데이터추출부(250)는 재구성한 서브영상으로부터 금속으로 인한 손상영역의 데이터를 추출한다.
도 6을 참조하면, 선형조합값산출부(600)는 제1사이노그램(sinogram)값산출부(610), 선형계수값계산부(630), 선형조합값계산부(650)를 포함하여 구성된다.
제1사이노그램(sinogram)값산출부(610)는 분할된 이미지로부터 1사이노그램(sinogram)값을 산출한다.
선형계수값계산부(630)는 사이노그램(sinogram)값을 기저값으로 하여 비손상영역에 근사하도록 선형계수값을 산출한다. 이 때 선형계수값계산부(S630)는 수학식 1로 정의되는 선형계수 계산식을 이용하여여
Figure 112016114771226-pat00022
이 최대한 영상이미지데이터값(
Figure 112016114771226-pat00023
)과 근사한 선형계수값(
Figure 112016114771226-pat00024
)을 계산한다.
[수학식 1]
Figure 112016114771226-pat00025
선형조합값계산부(650)는 계산한 선형계수값을 이용하여 선형조합값을 계산한다. 이 때 선형조합값계산부(650)는 수학식 2로 정의되는 선형조합 계산식을 이용하여 각각의 선형계수값(
Figure 112016114771226-pat00026
)와 기저값(
Figure 112016114771226-pat00027
) 곱의 합으로 계산한다.
[수학식 2]
Figure 112016114771226-pat00028
잔여오류값계산단계(710)는 계산한 선형조합값에서 영상이미지데이터값을 뺀 값인 잔여오류값을 계산한다.
도 7을 참조하면, 최종보정이미지출력결정부(700)는 잔여오류값계산부(710), 제2선형보간부(730), 보정사이노그램(sinogram)값계산부(750), 최종보정이미지출력여부판단부(770)를 포함하여 구성된다.
잔여오류값계산부(710)는 계산한 선형조합값에서 영상이미지데이터값을 뺀 값인 잔여오류값을 계산한다.
제2선형보간부(730)는 계산한 잔여오류값의 손상영역을 선형보간한다.
보정사이노그램(sinogram)값계산부(750)는 손상영역데이터값과 선형조합값을 더한 값인 보정사이노그램(sinogram)값을 계산한다.
최종보정이미지출력여부판단부(770)는 계산된 잔여오류값과 기설정한 임계값의 비교를 통해 최종보정이미지 출력여부를 판단한다. 이 때 잔여오류값이 기설정한 임계값보다 클 경우 제1선형보간부(300) 내지 최종보정이미지출력여부판단부(770)를 실행하고, 잔여오류값이 기설정한 임계값보다 작을 경우 보정사이노그램(sinogram)값을 제2메인영상으로 재구성하고, 이를 출력하는 것을 특징으로 한다.
도 8을 참조하면, 본 발명의 일 실시예에 따른 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정 과정을 구체적으로 제시한다.
1. 먼저 금속이 포함된 인체에 콘빔 CT를 투사한다. 콘빔 CT 투사로부터 얻어진 데이터인
Figure 112016114771226-pat00029
에는 금속으로 인하여 손상된 부분이 포함되어있다. 반복(iteration)을 위한 상수 k =1 으로 놓는다.
2. 얻어진
Figure 112016114771226-pat00030
을 보정 없이 재구성하면 금속 인공물을 포함한 영상이 재구성된다.
3. 2번 그림을 기초로 금속 부분을 추출한다. 이 때 사용될 수 있는 금속 추출 알고리즘은 임계값(threshold) 기반의 알고리즘과 액티브 컨투어(active contour) 기반의 알고리즘 등을 사용할 수 있다.
4. 추출한 금속 부분으로부터 금속 위치 정보를 이용하여 금속으로 인하여 손상되는 부분을 계산한다. 이 부분을 금속 트레이스(metal trace)라 한다.
5. 상기 얻은
Figure 112016114771226-pat00031
중에서 손상되지 않은 부분을 이용하여 금속 트레이스(metal trace)를 선형보간(Linear Interpolation)한다. 이 때 선형보간(Linear Interpolation)된 데이터를
Figure 112016114771226-pat00032
이라 한다.
6.
Figure 112016114771226-pat00033
을 영상으로 재구성한다. 이를
Figure 112016114771226-pat00034
라 하고,
Figure 112016114771226-pat00035
,
Figure 112016114771226-pat00036
이라고 한다.
7.
Figure 112016114771226-pat00037
영상을 분할하여 각기 다른 이미지 영역대로 나눈다. 이 때 적용될 수 있는 요소 기술로는 재귀적 액티브 컨투어(RAC) 알고리즘 등이 있다.
8. 서로 다른 값을 갖는 영역대로 분할한 이미지들을
Figure 112016114771226-pat00038
이라고 하면 이를 각각 콘빔 CT에 투사한 사이노그램(sinogram)들을
Figure 112016114771226-pat00039
이라고 한다.
Figure 112016114771226-pat00040
은 사이노그램(sinogram)상의 기저값으로 하여 비손상영역에 근사하도록 선형계수값을 계산한다. 이 때 수학식 1로 정의되는 선형계수 계산식을 이용하여
Figure 112016114771226-pat00041
이 최대한 영상이미지데이터값(
Figure 112016114771226-pat00042
)과 근사한 선형계수값(
Figure 112016114771226-pat00043
)을 계산한다.
[수학식 1]
Figure 112016114771226-pat00044
9. 산출된 선형계수값을 이용하여 선형조합값을 계산한다. 이 때 수학식 2로 정의되는 선형조합 계산식을 이용하여 각각의 선형계수값(
Figure 112016114771226-pat00045
)와 기저값(
Figure 112016114771226-pat00046
) 곱의 합으로 계산한다.
[수학식 2]
Figure 112016114771226-pat00047
10. 9번에서 얻어진
Figure 112016114771226-pat00048
과 상기 얻은
Figure 112016114771226-pat00049
사이의 차를 잔여오류값(residual)
Figure 112016114771226-pat00050
라고 한다. (
Figure 112016114771226-pat00051
)
11.
Figure 112016114771226-pat00052
에서의 금속 트레이스(metal trace)를 선형보간(Linear Interpolation)하고, 선형보간(Linear Interpolation)된
Figure 112016114771226-pat00053
을
Figure 112016114771226-pat00054
라고 한다.
12. k=k+1로 업데이트 하고, k 번째 보정된 사이노그램
Figure 112016114771226-pat00055
이라 한다.
13. 반복여부를 판단하기 위하여
Figure 112016114771226-pat00056
가 기설정한 임계값(threshold)보다 크면
Figure 112016114771226-pat00057
을 재구성한 영상 이라
Figure 112016114771226-pat00058
하고, 7번으로 돌아가서 7번 내지 12번을 반복한다. 반면
Figure 112016114771226-pat00059
가 기설정한 임계값(threshold)보다 작으면 상기 반복 과정을 멈추고 보정사이노그램(sinogram)(
Figure 112016114771226-pat00060
)을 제2메인영상(
Figure 112016114771226-pat00061
)으로 재구성하고, 이를 출력한다.
도 9를 참조하면, 본 발명의 이미지분할과정에서 적용 가능한 요소기술 중 하나인 재귀적 액티브 컨투어(RAC) 알고리즘을 구체적으로 제시한다.
먼저 분류하고자 하는 영상을
Figure 112016114771226-pat00062
이라 할 때,
Figure 112016114771226-pat00063
에 첫 번째 액티브 컨투어(active contour)를 적용하면 영상은 전경과 배경 부분으로 나뉜다(510). 이 때 흰색 영역은 전경, 검정색 영역은 배경이라 할 때, 전경 부분의 영상 값 중 가장 낮은 픽셀의 값으로 배경 부분의 영상 값을 변환하면
Figure 112016114771226-pat00064
이 된다. 전경 부분의 영상 값은
Figure 112016114771226-pat00065
이 해당되며, 그 중 가장 낮은 픽셀 값이 최소값이다. 따라서
Figure 112016114771226-pat00066
전경 부분의 최소값은
Figure 112016114771226-pat00067
가 된다.
Figure 112016114771226-pat00068
에 두 번째 액티브 컨투어(active contour)를 적용하면 영상은 전경과 배경 부분으로 나뉜다(520). 마찬가지로 흰색 영역은 전경, 검정색 영역은 배경이라 할 때, 전경 부분의 영상 값 중 가장 낮은 픽셀의 값으로 배경 부분의 영상 값을 변환하면
Figure 112016114771226-pat00069
가 된다. 전경 부분의 영상 값은
Figure 112016114771226-pat00070
가 해당되며, 그 중 가장 낮은 픽셀 값이 최소값이다. 따라서
Figure 112016114771226-pat00071
전경 부분의 최소값은
Figure 112016114771226-pat00072
가 된다.
Figure 112016114771226-pat00073
에 세 번째 액티브 컨투어(active contour)를 적용하면 영상은 전경과 배경 부분으로 나뉜다(530). 마찬가지로 흰색 영역은 전경, 검정색 영역은 배경이라 할 때, 전경 부분의 영상 값 중 가장 낮은 픽셀의 값으로 배경 부분의 영상 값을 변환하면
Figure 112016114771226-pat00074
가 된다. 전경 부분의 영상 값은
Figure 112016114771226-pat00075
이 해당되며,
Figure 112016114771226-pat00076
은
Figure 112016114771226-pat00077
전경 부분의 최소값이 된다.
Figure 112016114771226-pat00078
은 전경 부분의 영상 값이 0에 수렴하다고 판단된 상태이며
Figure 112016114771226-pat00079
의 영상 전체는 배경 부분이 되었다. 이 때 액티브 컨투어(active contour) 반복 과정을 멈추고, 첫 번째 액티브 컨투어(active contour)를 적용한 영상(510), 두 번째 액티브 컨투어(active contour)를 적용한 영상(520) 및 세 번째 액티브 컨투어(active contour)를 적용한 영상(530)을 합하면 최종 세그먼테이션(segmentation) 결과(
Figure 112016114771226-pat00080
)가 된다.
도 10을 참조하면, 본 발명의 이미지분할과정에서 적용 가능한 요소기술 중 하나인 재귀적 액티브 컨투어(RAC) 알고리즘을 적용한 결과이다. 오른쪽 그림과 같이 RAC가 적용된 결과 빨간색과 파란색 영역으로 구분된 경계가 왼쪽 그림보다 뚜렷한 것이 확인 가능하다.
도 11을 참조하면, 본 발명의 일 실시예에 따른 콘빔 CT 영상의 금속 아티팩트 영역 보정 결과로서 여러 개의 프라이어(prior)를 앞으로 투사하여 만든 기저 사이노그램(sinogram)을 이용하였다. 왼쪽 그림과 같이 금속 인공물로 인해 밝게 나타나는 부분이 생기므로 적정한 임계값을 찾는 것은 어려웠다. 하지만 본 발명에서 제시한 바와 같이 여러 개의 프라이어(prior)를 앞으로 투사하여 만든 기저 사이노그램(sinogram)을 이용한 결과, 오른쪽 그림과 같이 금속 인공물로 인한 아티팩트가 효과적으로 감소한 것을 볼 수 있다.
이상에서 본 발명의 바람직한 실시예를 설명하였으나, 본 발명은 다양한 변화와 변경 및 균등물을 사용할 수 있다. 본 발명은 상기 실시예를 적절히 변형하여 동일하게 응용할 수 있음이 명확하다. 따라서 상기 기재 내용은 하기 특허청구범위의 한계에 의해 정해지는 본 발명의 범위를 한정하는 것이 아니다.
100 : 영상이미지수신부 200 : 손상영역추출부
210 : 영상이미지데이터값산출부 230 : 서브영상재구성부
250 : 손상영역데이터추출부 300 : 제1선형보간부
400 : 제1메인영상재구성부 500 : 이미지분할부
510 :
Figure 112016114771226-pat00081
에 첫 번째 액티브 컨투어를 적용한 영상
520 :
Figure 112016114771226-pat00082
에 두 번째 액티브 컨투어를 적용한 영상
530 :
Figure 112016114771226-pat00083
에 세 번째 액티브 컨투어를 적용한 영상
600 : 선형조합값산출부 610 : 제1사이노그램(sinogram)값산출부
630 : 선형계수값계산부 650 : 선형조합값계산부
700 : 최종보정이미지출력결정부 710 : 잔여오류값계산부
730 : 제2선형보간부 750 : 보정사이노그램(sinogram)값계산부
770 : 최종보정이미지출력여부판단부

Claims (14)

  1. 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치에 의해 수행되는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법에 있어서,
    CT(컴퓨터단층촬영) 촬영된 영상의 이미지를 영상이미지수신부에서 입력받는 영상이미지수신단계;
    상기 입력받은 영상의 이미지로부터 손상영역추출부에서 영상이미지데이터값을 산출하여 금속으로 인한 손상영역을 추출하는 손상영역추출단계;
    상기 추출된 손상영역을 제1선형보간부에서 선형보간(Linear Interpolation)하는 제1선형보간단계;
    상기 보간된 손상영역데이터값을 제1메인영상재구성부에서 제1메인영상으로 재구성하는 제1메인영상재구성단계;
    상기 재구성된 영상을 이미지분할부에서 하나 또는 둘 이상의 이미지로 분할하는 이미지분할단계;
    상기 분할된 이미지를 기초로 선형조합값산출부에서 선형조합값을 산출하는 선형조합값산출단계;
    상기 산출한 선형조합값과 상기 영상이미지데이터값을 기초로 최종보정이미지출력결정부에서 최종보정이미지 출력여부를 결정하는 최종보정이미지출력결정단계;를 포함하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법.
  2. 제 1항에 있어서,
    상기 선형조합값산출단계는,
    상기 분할된 하나 또는 둘 이상의 이미지로부터 제1사이노그램(sinogram)값을 산출하는 제1사이노그램(sinogram)값산출단계;
    상기 산출된 하나 또는 둘 이상의 제1사이노그램(sinogram)값을 기저값으로 하여 비손상영역에 근사하도록 선형계수값을 계산하는 선형계수값계산단계; 및
    상기 계산된 선형계수값을 이용하여 선형조합값을 계산하는 선형조합값계산단계;를 포함하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법.
  3. 제 2항에 있어서,
    상기 선형계수값계산단계는,
    수학식 1로 정의되는 선형계수 계산식을 이용하여
    Figure 112016114771226-pat00084
    이 최대한 영상이미지데이터값(
    Figure 112016114771226-pat00085
    )과 근사한 선형계수값(
    Figure 112016114771226-pat00086
    )을 계산하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법.
    [수학식 1]
    Figure 112016114771226-pat00087
  4. 제 2항에 있어서,
    상기 선형조합값 계산단계는,
    수학식 2로 정의되는 선형조합 계산식을 이용하여 각각의 선형계수값(
    Figure 112016114771226-pat00088
    )와 기저값(
    Figure 112016114771226-pat00089
    ) 곱의 합으로 계산하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법.
    [수학식 2]
    Figure 112016114771226-pat00090
  5. 제 1항에 있어서,
    상기 손상영역추출단계는,
    상기 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출하는 영상이미지데이터값산출단계;
    상기 산출된 영상이미지데이터값을 서브영상으로 재구성하는 서브영상재구성단계; 및
    상기 재구성된 서브영상으로부터 금속으로 인한 손상영역의 데이터을 추출하는 손상영역데이터추출단계;를 포함하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법.
  6. 제 1항에 있어서,
    상기 최종보정이미지출력결정단계는,
    상기 계산된 선형조합값에서 상기 영상이미지데이터값을 뺀 값인 잔여오류값을 계산하는 잔여오류값계산단계;
    상기 계산된 잔여오류값의 손상영역을 선형보간(Linear Interpolation)하는 제2선형보간단계;
    상기 보간된 손상영역데이터값과 상기 선형조합값을 더한 값인 보정사이노그램(sinogram)값을 계산하는 보정사이노그램(sinogram)값계산단계; 및
    상기 계산된 잔여오류값과 기설정한 임계값의 비교를 통해 최종보정이미지 출력여부를 판단하는 최종보정이미지출력여부판단단계;를 포함하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법.
  7. 제 6항에 있어서,
    상기 최종보정이미지출력여부판단단계는,
    상기 잔여오류값이 기설정한 임계값보다 클 경우 상기 제1선형보간단계 내지 상기 최종보정이미지출력결정단계를 실행하고,
    상기 잔여오류값이 기설정한 임계값보다 작을 경우 상기 계산된 보정사이노그램(sinogram)값을 제2메인영상으로 재구성하고, 이를 출력하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정방법.
  8. 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치에 있어서,
    CT(컴퓨터단층촬영) 촬영된 영상의 이미지를 입력받는 영상이미지수신부;
    상기 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출하여 금속으로 인한 손상영역을 추출하는 손상영역추출부;
    상기 추출된 손상영역을 선형보간(Linear Interpolation)하는 제1선형보간부;
    상기 보간된 손상영역데이터값을 제1메인영상으로 재구성하는 제1메인영상재구성부;
    상기 재구성된 영상을 하나 또는 둘 이상의 이미지로 분할하는 이미지분할부;
    상기 분할된 이미지를 기초로 선형조합값을 산출하는 선형조합값산출부;
    상기 산출한 선형조합값과 상기 영상이미지데이터값을 기초로 최종보정이미지 출력여부를 결정하는 최종보정이미지출력결정부;를 포함하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치.
  9. 제 8항에 있어서,
    상기 선형조합값산출부는,
    상기 분할된 하나 또는 둘 이상의 이미지로부터 제1사이노그램(sinogram)값을 산출하는 제1사이노그램(sinogram)값산출부;
    상기 산출된 하나 또는 둘 이상의 제1사이노그램(sinogram)값을 기저값으로 하여 비손상영역에 근사하도록 선형계수값을 계산하는 선형계수값계산부; 및
    상기 계산된 선형계수값을 이용하여 선형조합값을 계산하는 선형조합값계산부;를 포함하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치.
  10. 제 9항에 있어서,
    상기 선형계수값계산부는,
    수학식 1로 정의되는 선형계수 계산식을 이용하여
    Figure 112016114771226-pat00091
    이 최대한 영상이미지데이터값(
    Figure 112016114771226-pat00092
    )과 근사한 선형계수값(
    Figure 112016114771226-pat00093
    )을 계산하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치.
    [수학식 1]
    Figure 112016114771226-pat00094
  11. 제 9항에 있어서,
    상기 선형조합값 계산부는,
    수학식 2로 정의되는 선형조합 계산식을 이용하여 각각의 선형계수값(
    Figure 112016114771226-pat00095
    )와 기저값(
    Figure 112016114771226-pat00096
    ) 곱의 합으로 계산하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치.
    [수학식 2]
    Figure 112016114771226-pat00097
  12. 제 8항에 있어서,
    상기 손상영역추출부는,
    상기 입력받은 영상의 이미지로부터 영상이미지데이터값을 산출하는 영상이미지데이터값산출부;
    상기 산출된 영상이미지데이터값을 서브영상으로 재구성하는 서브영상재구성부; 및
    상기 재구성된 서브영상으로부터 금속으로 인한 손상영역의 데이터을 추출하는 손상영역데이터추출부;를 포함하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치.
  13. 제 8항에 있어서,
    상기 최종보정이미지출력결정부는,
    상기 계산된 선형조합값에서 상기 영상이미지데이터값을 뺀 값인 잔여오류값을 계산하는 잔여오류값계산부;
    상기 계산된 잔여오류값의 손상영역을 선형보간(Linear Interpolation)하는 제2선형보간부;
    상기 보간된 손상영역데이터값과 상기 선형조합값을 더한 값인 보정사이노그램(sinogram)값을 계산하는 보정사이노그램(sinogram)값계산부; 및
    상기 계산된 잔여오류값과 기설정한 임계값의 비교를 통해 최종보정이미지 출력여부를 판단하는 최종보정이미지출력여부판단부;를 포함하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치.
  14. 제 13항에 있어서,
    상기 최종보정이미지출력여부판단부는,
    상기 잔여오류값이 기설정한 임계값보다 클 경우 상기 제1선형보간부 내지 상기 최종보정이미지출력결정부를 실행하고,
    상기 잔여오류값이 기설정한 임계값보다 작을 경우 상기 계산된 보정사이노그램(sinogram)값을 제2메인영상으로 재구성하고, 이를 출력하는 것을 특징으로 하는 콘빔 CT 영상의 금속으로 인한 아티팩트 영역 보정장치.
KR1020160156707A 2016-11-23 2016-11-23 콘빔 ct 영상의 금속으로 인한 아티팩트 보정방법 및 장치 KR101870890B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020160156707A KR101870890B1 (ko) 2016-11-23 2016-11-23 콘빔 ct 영상의 금속으로 인한 아티팩트 보정방법 및 장치

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020160156707A KR101870890B1 (ko) 2016-11-23 2016-11-23 콘빔 ct 영상의 금속으로 인한 아티팩트 보정방법 및 장치

Publications (2)

Publication Number Publication Date
KR20180058062A KR20180058062A (ko) 2018-05-31
KR101870890B1 true KR101870890B1 (ko) 2018-06-25

Family

ID=62454672

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020160156707A KR101870890B1 (ko) 2016-11-23 2016-11-23 콘빔 ct 영상의 금속으로 인한 아티팩트 보정방법 및 장치

Country Status (1)

Country Link
KR (1) KR101870890B1 (ko)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102248191B1 (ko) * 2019-04-19 2021-05-04 연세대학교 산학협력단 콘빔 컴퓨터 단층촬영 영상에서 콘빔아티팩트를 보정하는 방법 및 장치, 그것을 포함하는 콘빔 컴퓨터 단층촬영 장치

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20110019264A (ko) * 2009-08-19 2011-02-25 (주)바텍이우홀딩스 압축 센싱 이론을 이용한 ｘ-선 전산 단층 촬영 영상의 금속성 음영 제거 방법
KR101018051B1 (ko) * 2008-10-23 2011-03-02 한국과학기술원 엑스선 촬영 영상에서 메탈 아티팩트를 감소시키는 방법
JP2011156302A (ja) * 2010-02-03 2011-08-18 Kyoto Univ Ｘ線ｃｔ画像処理方法，ｘ線ｃｔプログラムおよび該プログラムが搭載されたｘ線ｃｔ装置
CN102567958B (zh) * 2010-12-31 2014-06-25 飞利浦(中国)投资有限公司 一种去除ct图像中金属伪影的图像后处理方法
JP2014530709A (ja) * 2011-10-24 2014-11-20 コーニンクレッカ フィリップス エヌ ヴェ Ｃｔスライス画像のための動作補償されたセカンドパス金属アーティファクト補正

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101018051B1 (ko) * 2008-10-23 2011-03-02 한국과학기술원 엑스선 촬영 영상에서 메탈 아티팩트를 감소시키는 방법
KR20110019264A (ko) * 2009-08-19 2011-02-25 (주)바텍이우홀딩스 압축 센싱 이론을 이용한 ｘ-선 전산 단층 촬영 영상의 금속성 음영 제거 방법
JP2011156302A (ja) * 2010-02-03 2011-08-18 Kyoto Univ Ｘ線ｃｔ画像処理方法，ｘ線ｃｔプログラムおよび該プログラムが搭載されたｘ線ｃｔ装置
CN102567958B (zh) * 2010-12-31 2014-06-25 飞利浦(中国)投资有限公司 一种去除ct图像中金属伪影的图像后处理方法
JP2014530709A (ja) * 2011-10-24 2014-11-20 コーニンクレッカ フィリップス エヌ ヴェ Ｃｔスライス画像のための動作補償されたセカンドパス金属アーティファクト補正

Also Published As

Publication number Publication date
KR20180058062A (ko) 2018-05-31

Similar Documents

Publication Publication Date Title
CN106683144B (zh) 一种图像迭代重建方法及装置
US10092266B2 (en) Metal artifact correction and noise-adjustment of CT scan image
US7831097B2 (en) System and method for image reconstruction
US10395396B2 (en) X-ray CT apparatus, reconfiguration arithmetic apparatus, and X-ray CT image generation method
US7929746B2 (en) System and method for processing imaging data
JP6026214B2 (ja) 連続マルチスケール再構成において詳細画像を補うｘ線コンピュータ断層撮像装置（ｘ線ｃｔ装置）、医用画像処理装置及び医用画像処理方法
EP3438922B1 (en) Method for processing at least one x-ray image
CN103026379A (zh) 推算图像噪音水平的方法
US9514549B2 (en) Method for reducing metal artifact in computed tomography
CN111223156B (zh) 一种用于牙科锥束ct系统的金属伪影消除方法
CN110533738B (zh) 重建数据处理方法、装置、医学成像系统及存储介质
US8768045B2 (en) Method for acquiring a 3D image dataset freed of traces of a metal object
US20070071296A1 (en) Radiographic image processing apparatus for processing radiographic image taken with radiation, method of radiographic image processing, and computer program product therefor
CN116485809B (zh) 基于自注意力与感受野调整的牙齿实例分割方法及系统
WO2021165053A1 (en) Out-of-distribution detection of input instances to a model
KR101870890B1 (ko) 콘빔 ct 영상의 금속으로 인한 아티팩트 보정방법 및 장치
US11631160B2 (en) Image processing apparatus, image processing method, and X-ray CT apparatus
CN113793272A (zh) 图像降噪方法及装置、存储介质、终端
US10932743B2 (en) Determining image values in marked pixels of at least one projection image
US10969771B2 (en) Computed tomography for non-destructive evaluation of manufactured parts
Kim et al. Blind single image super resolution with low computational complexity
US20100061656A1 (en) Noise reduction of an image signal
CN111247555B (zh) 用于校正ct图像的方法、计算机程序产品和计算机系统
KR20220071554A (ko) 의료 영상 융합 시스템
US8180124B2 (en) Registration method for projections in dual energy

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