KR20190013161A - 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법 - Google Patents

비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법 Download PDF

Info

Publication number
KR20190013161A
KR20190013161A KR1020170097322A KR20170097322A KR20190013161A KR 20190013161 A KR20190013161 A KR 20190013161A KR 1020170097322 A KR1020170097322 A KR 1020170097322A KR 20170097322 A KR20170097322 A KR 20170097322A KR 20190013161 A KR20190013161 A KR 20190013161A
Authority
KR
South Korea
Prior art keywords
positron emission
pet
emission tomography
image
data
Prior art date
Application number
KR1020170097322A
Other languages
English (en)
Other versions
KR101958098B1 (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 KR1020170097322A priority Critical patent/KR101958098B1/ko
Publication of KR20190013161A publication Critical patent/KR20190013161A/ko
Application granted granted Critical
Publication of KR101958098B1 publication Critical patent/KR101958098B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • 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/037—Emission tomography
    • 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/10072—Tomographic images
    • G06T2207/10104—Positron emission tomography [PET]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/416—Exact reconstruction

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Biophysics (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)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine (AREA)

Abstract

본 발명 비스듬히 배치한 섬광결정(120)을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법은, 섬광결정(120)을 비스듬히 배치한 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)로 획득한 list mode 데이터를 재배열하는 단계와; 재배열한 데이터를 통해 영상재구성하는 단계; 로 이루어지는 것으로,
본 발명은 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법을 통하여, 섬광결정을 비스듬히 배치한 양전자 방출 단층 촬영기(PET)와, 상기 양전자 방출 단층 촬영기(PET)로 획득한 데이터의 재배열 및 영상재구성 방법을 사용함으로써, 민감도를 증가시키고, Sinogram에서의 갭의 발생 가능성을 줄이고, 중심축방향의 해상도 향상 및 전반적인 영상재구성 향상을 제공하는 현저한 효과가 있다.

Description

비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법{image reconstruction method of positron emission tomography with scintillation crystal disposed obliquely}
본 발명은 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법에 관한 것으로, 더욱 자세하게는 섬광결정을 비스듬히 배치한 양전자 방출 단층 촬영기(PET)와, 상기 양전자 방출 단층 촬영기(PET)로 획득한 데이터의 재배열 및 영상재구성 방법을 사용함으로써, 민감도를 증가시키고, Sinogram에서의 갭의 발생 가능성을 줄이고, 중심축방향의 해상도 향상 및 전반적인 영상재구성 향상을 제공하는 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법에 관한 것이다.
양전자 방출 단층 촬영기(PET, Positron Emission Tomography)는 인체에 주사 또는 흡입된 방사성의약품의 체내분포를 영상의 형태로 표현하는 의료기기로, 방사성동위원소에서 양전자가 방출 직후에 전자와 결합하여 발생하는 쌍소멸에 의해 거의 180도 반대방향으로 방출되는 두 개의 감마선을 검출하여 얻은 정보를 이용하여 영상을 재구성하는 장치이다.
예를 들어, 2개의 검출기인 b1,b2가 각각 t1,t2 시간에 동일한 양전자방출에 의한 감마선을 검출하였다면, 2개의 검출기를 연결하는 어느 중간부분에 양전자가 하나 존재했다는 것을 의미한다. 특히, b1의 위치, b2의 위치, t1 및 t2를 포함하는 관측정보가 나열된 데이터를 list mode 데이터라 하며, 상기 list mode 데이터에서 검출기로 수집한 데이터는 Raw Data이되,
Figure pat00001
로 표현하고, 상기
Figure pat00002
는 2개의 검출기인 b1,b2을 잇는 반응영역에 분포한 양전자의 개수이다.
상기 2개의 검출기인 에서 동시에 관측된 감마선을 초래한 양전자가 위치하는 반응영역은 대응선(LOR, Line-Of-Response)과 대응관(TOR, Tube-Of-Response)를 통해 관측현상을 모델화하여 결정한다. 상기 반응영역에서 관측된 데이터는 Raw Data이되,
Figure pat00003
으로 표현하며, 대응선(LOR)을 통해 모델화된 모델을 사용해서 데이터를 획득하되, 검출기의 배치와는 무관하게 미리 정해진 각도 θ와, 중심축으로부터의 거리 u와, 높이 h와, 중심축에 대해 직각에서 벗어난 정도 d에 따라 재배열한 데이터를 Sinogram이라 하고, Yθ,u,h,d로 표현한다.
예를 들어, 종래의 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)의 경우에는 개의 섬광결정을 포함하고, 상기 섬광결정이 검출기의 중심축에서 직각방향으로 배치된 것으로, 다수 개의 Sinogram을 가지게 되는 것이며, 상기 Sinogram는 r2개의 θ,u를 인덱스로 갖는 것이다. 특히, 다수 개의 Sinogram 중에서, r개의 Sinogram에 속하는 데이터는 이론적으로 중심축에 대해 직각인 평면상에서 방출된 것에 대응하는 것이다.
상기 Sinogram을 사용하는 포맷은 검출기의 배치와는 무관하게 미리 정해진 규칙에 따라 재배열할 수 있는 장점이 있으나, 재배열 과정에서 오류가 발생할 여지가 있다.
상기 관측영역에서 관측된 데이터를 대응선(LOR)을 통해 모델화된 모델을 사용해서 물리적 현상을 분석할 수 있으며, 상기 물리적 현상은 Positron Range, Photon non-collinearity, Depth-of-Interaction 등이 있다.
상기 Positron Range는 양전자의 원래의 위치와 쌍소멸이 벌어지는 위치의 차이인 것으로, 사용한 방사선의약품의 종류에 따라 다르지만 약 1mm정도이다. 상기 Positron Range가 궁극적으로 양전자 방출 단층 촬영기(PET)의 공간분해능에 대한 한계를 결정한다.
상기 Photon non-collinearity는 양전자의 전자의 쌍소멸은 충돌 각도에 따라 쌍소멸의 결과인 감마선의 방출이 180도에서 약간 다르게 됨으로써, 데이터 분석 시 부정적인 현상을 띄는 현상이다.
상기 Depth-of-Interaction은 검출기의 민감도를 높이려는 과정에서 등장한다. 상기 검출기의 민감도를 높이는 방법으로 표면에 비해 길이가 상대적으로 긴 검출기를 주로 사용하고 있으며, 광자가 검출기의 옆면을 통과하여 관측되는 경우, 정확한 반응깊이(Depth-of-Interaction, DOI)를 추정하는 것이 어렵게 됨으로써, 공간해상도의 저하로 나타난다.
또한, 상기 양전자 방출 단층 촬영기(PET)로 획득한 데이터를 통해 영상을 재구성하는 알고리즘은 list mode 데이터를 어떻게 분석하느냐에 따라 사용하는 데이터가 구분되는 것이며, 상기 데이터는 Raw Data
Figure pat00004
또는 재배열한 Sinogram Yθ,u,h,d를 사용할 수 있는 것이다.
상기 Raw Data
Figure pat00005
를 사용하는 경우에는 관측된 list mode 데이터를 Raw Data
Figure pat00006
로 재배열하고 TOR 모델을 사용하여 분석한다.
즉, 상기 Raw Data
Figure pat00007
과 재구성하고자 하는 방사성의약품의 체내분포영상
Figure pat00008
와의 선형관계를 정의하는 시스템 행렬
Figure pat00009
을 미리 계산한 후,
Figure pat00010
,
Figure pat00011
에 반복적 수치해석방법을 적용하여 방사성의약품의 체내분포영상
Figure pat00012
를 반복적으로 적용한다. 이때, v는 화소 인덱스이며,
Figure pat00013
는 화소 v에서의 방사성의약품의 체내분포밀도이다.
상기 반복적 수치해석방법으로 가장 널리 사용되고 있는 방법은 최대우도 기댓값-최대화(ML-EM, Maximum Likelihood-Expectation Maximization)이다.
상기 Raw Data에 대응하는 시스템 행렬을 반복법에 적용하는 방법은 관측된 데이터와 재구성하고자 하는 영상과의 관계를 정확히 구현할 수 있다는 장점이 있으나, 관측된 데이터와 재구성하고자 하는 영상과의 관계를 결정하는 양전자 방출 단층 촬영기(PET)의 설계적 특징과, 물리적 현상을 많이 고려할수록 영상재구성 과정에서 필요한 투사 및 역투사 과정에 많은 계산이 필요하며, 또한, 시스템 행렬
Figure pat00014
를 미리 계산하여 저장해 두기위해 많은 메모리를 필요로 한다.
상기 재배열한 Sinogram Yθ,u,h,d를 사용하는 경우에는 관측된 list mode 데이터를 LOR 데이터로 이용하여 재배열한 Sinogram Yθ,u,h,d을 얻고, 상기 재배열된 데이터
Figure pat00015
와 재구성하고자 하는 방사성의약품의 체내분포영상
Figure pat00016
와의 선형관계
Figure pat00017
에 반복적 수치해석방법을 적용하여 방사성의약품의 체내분포영상
Figure pat00018
를 반복적으로 적용한다.
즉, 상기 재배열된 데이터
Figure pat00019
와 재구성하고자 하는 방사성의약품의 체내분포영상
Figure pat00020
와의 선형관계를 정의하는 시스템 행렬
Figure pat00021
는 미리 계산하여 저장된 것을 사용하기 보다는 필요할 때 마다 계산하여 사용하는 방식을 이용한다. 상기와 같은 재배열과정은 행렬이용방식이 가능하게 하는 형태로 진행해야 한다.
상기 재배열된 Sinogram을 이용하는 방법은 Raw Data를 사용하는 방법에 비해 재배열 과정에서 오차가 발생하나, 시스템 행렬
Figure pat00022
를 미리 계산하여 저장해 둘 필요가 없고, 영상재구성에 상대적으로 적은 계산을 필요로 하므로, 영상재구성 방법은 재배열된 Sinogram을 이용할 수 있는 것이다.
그리고 상기 재배열된 Sinogram을 이용하는 방법의 또 다른 장점은 양전자 방출 단층 촬영기(PET)의 영상재구성에 널리 사용되는 여현역투사(FBP, Filtered Backprojection)방법에 적용가능하다는 것이다. 상기 여현역투사(FBP)는 동일 평면상의 얻은 Sinogram Yθ,u,h,d에 대해 ‘1차원 푸리에 변환’, ‘램프필터링’, ‘1차원 푸리에 역변환’, ‘2차원 역투사’를 순차적으로 적용하여 2차원 영상을 재구성하는 과정으로, 도식화하면
Figure pat00023
와 같다. 이때, 상기
Figure pat00024
는 높이 h에서의 f의 단층영상이며, 상기 FFT는 Fast Fourier Transform이며, 상기 IFFT는 Inverse Fast Fourier Transform이다.
상기 여현역투사(FBP)를 이용한 전체 3차원 영상재구성은 2차원 단층영상을 순차적으로 재구성함으로써 얻는 것이다.
예를 들어, 종래의 양전자 방출 단층 촬영기(PET)의 경우에는 가로[r]*세로[c]개의 섬광결정을 포함하고, 상기 섬광결정이 검출기의 중심축에서 직각방향으로 배치된 것으로, r개의 동일 평면상에서 발생한 얻은 Sinogram 집합 Yθ,u,h,d에 여현역투사(FBP)를 이용하여 r개의 단층영상을 재구성한다.
또한, 동일 평면상에 있지 않은 r2-r개의 Sinogram에서 Yθ,u,h,d, d≠0 를 이용하여, 동일 평면상의 r개의 Sinogram에서 Yθ,u,h,d 중간에 있는 평면 Sinogram을 r-1개만큼 보간법을 이용하여 생성한 후, 단층영상을 재구성한다.
즉, 상기 여현역투사(FBP)는 2r-1개의 평행단층내의 Sinogram인 Direct Sinogram을 이용하여 2r-1개의 단층영상을 재구성하여 3차원 영상재구성을 수행한다.
상기 여현역투사(FBP)를 하는 방법은 고속 푸리에 변환을 이용하여 빠른 영상재구성을 제공하나, 주어진 데이터가 모든 각도와 위치에 충분히 존재해야 한다는 제한조건이 있고, 중심축에 대해 직각평면에서 많이 벗어난 데이터를 이용하기가 어렵고, 검출기간의 간격이 큰 경우, 특히, 검출기 링을 여러 겹 쌓아 만든 양전자 방출 단층 촬영기(PET)가 검출기 링간의 간격이 큰 경우에는 보간법에 의한 평면에서의 Sinogram이 부정확하여 재구성 영상에 왜곡현상이 발생한다.
상기 양전자 방출 단층 촬영기(PET)의 데이터는 여러 종류의 관측오차를 필연적으로 수반하게 된다. 상기 관측오차의 원인으로는 감마선이 인체내에서 흡수되어 관측하지 못하는 현상인 감쇠와, 감마선의 산란현상에 의해 실제 위치가 아닌 다른 위치에서의 관측되는 현상인 산란과 함께, 양전자의 방출, 양전자-전자의 쌍소멸, 감마선의 발생의 과정에 내재한 임의성(Randomness)을 들 수 있다. 기술한 관측의 임의성은 포아송 분포(Poisson distribution)를 따르는 것으로, 상기 임의성에 의한 관측오차를 줄이는 가장 효과적 방법은 관측횟수, 즉 민감도를 증가시키는 것이다.
상기 양전자 방출 단층 촬영기(PET)의 민감도를 증가시키는 방법의 하나로 논문 [1] (Yang et.a al, “Tapered LSO Arrays for Small Animal PET”, Physics in Medicine and Biology, 56(1), 139-153, 2011)과 논문 [2](James and Cherry, "Simulation study of spatial resolution and sensitivity for the tapered depth of interaction PET detectors for small animal imaging", Physics in Medicine and Biology, 55, N63-N74, 2010)에서는 사다리꼴 모양의 섬광결정을 사용하는 것을 제시하였다. 그러나 상기 사다리꼴 모양의 섬광결정은 직사각형 모양의 섬광결정에 비해 제작단가가 높다.
상기 양전자 방출 단층 촬영기(PET)의 시차오차를 감소시키는 방법으로 논문 [3] (“N.-Y. Lee and Y. Choi, "Simulation studies on depth of interaction effect correction using a Monte Carlo computed system matrix for brain positron emission tomography”, Computer Methods and Programs in Biomedicine, 108(2), 820-831, 2012)에서 TOR 관측모델을 이용한 시스템 행렬을 생성과 이를 이용한 DOI 보정을 제시하였으나, TOR 관측모델을 이용한 방법은 상기에서 언급한 바와 같이 많은 계산을 필요로 하다.
상기 양전자 방출 단층 촬영기(PET)의 시차오차를 감소시키는 방법의 하나는 다층 섬광결정을 이용하는 것이다. 살기 다층 섬광결정을 이용하는 경우에는 TOR 관측모델에 의한 시스템 행렬을 이용하지 않고도 상당한 정도의 DOI 보정을 할 수 있는데, 검출기모듈을 중심축에 대해 직각으로 배치하는 기존의 PET에서는 다층 섬광결정의 사용이 중심축방향의 해상도 향상에는 큰 도움을 주지 않는다.
상기와 같이 양전자 방출 단층 촬영을 하기위한 종래기술로써, 등록특허공보 등록번호 제10-1207710호의 양전자 단층촬영에서 영상의 해상도를 향상시키는 해상도 향상 장치 및 방법에는, 영상 의료 기기의 해상도 향상 장치에 있어서, 측정 대상체에 조사된 방사선에 응답하여 반응선을 검출하는 반응선 검출부; 상기 검출된 반응선으로부터 여현 곡선(sinogram)을 추출하는 여현 곡선 추출부; 및 상기 영상 의료 기기에 관한 블러 커널을 추정하고, 상기 추출된 여현 곡선을 고해상도의 여현 곡선으로 변환하는 고해상도 변환부를 포함하는 것을 특징으로 하는 해상도 향상 장치라고 기재되어 있다.
다른 종래기술로써, 등록특허공보 등록번호 제10-1356881호의 고해상도 양전자 방출 단층 촬영에서 병렬 처리를 위해 영상을 재구성하는 방법 및 장치에는, 측정 대상체에 조사된 방사선에 응답하여 검출기들로부터 검출된 반응선((Line-Of-Response; LOR)들을 콘 빔(cone beam) 형태의 여현 곡선 포맷(sinogram format)으로 변환하는 단계; 상기 변환된 반응선들을 기초로 역투영(back projection)를 수행하는 단계; 및 상기 역투영의 수행 결과를 이용하여 영상을 재구성하는 단계를 포함하는 병렬 처리를 위해 영상을 재구성하는 방법이라고 기재되어 있다.
그러나 상기와 같은 종래의 기술은 섬광결정이 직각으로 배치되거나, 또는 방사형으로 배치된 것으로, 민감도가 낮고, Sinogram에서의 갭의 발생 가능성이 높고, 중심축방향의 해상도 및 전반적인 영상재구성이 좋지 않은 단점이 있었다.
본 발명 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법에 있어서, 섬광결정을 비스듬히 배치한 양전자 방출 단층 촬영기(PET)와, 상기 양전자 방출 단층 촬영기(PET)로 획득한 데이터의 재배열 및 영상재구성 방법을 사용함으로써, 민감도를 증가시키고, Sinogram에서의 갭의 발생 가능성을 줄이고, 중심축방향의 해상도 향상 및 전반적인 영상재구성 향상을 제공하는 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법을 제공하고자 한다.
본 발명 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법은, 섬광결정을 비스듬히 배치한 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)로 획득한 list mode 데이터를 재배열하는 단계와; 재배열한 데이터를 통해 영상재구성하는 단계; 로 이루어지는 것을 특징으로 한다.
본 발명은 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법을 통하여, 섬광결정을 비스듬히 배치한 양전자 방출 단층 촬영기(PET)와, 상기 양전자 방출 단층 촬영기(PET)로 획득한 데이터의 재배열 및 영상재구성 방법을 사용함으로써, 민감도를 증가시키고, Sinogram에서의 갭의 발생 가능성을 줄이고, 중심축방향의 해상도 향상 및 전반적인 영상재구성 향상을 제공하는 현저한 효과가 있다.
도 1은 섬광결정의 집합으로 구성된 검출기모듈을 중심축으로부터의 방사방향에 대해 비스듬히 배치한 검출기 링이 설치된 양전자 방출 단층 촬영기(PET)의 평면도
도 2는 섬광결정의 집합으로 구성된 검출기모듈을 중심축의 직각에 대해 비스듬히 배치한 검출기 링이 설치된 양전자 방출 단층 촬영기의 측면도
도 3은 검출기 링을 구성하는 검출기모듈을 도시한 개념도
본 발명 비스듬히 배치한 섬광결정(120)을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법은, 섬광결정(120)을 비스듬히 배치한 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)로 획득한 list mode 데이터를 재배열하는 단계와; 재배열한 데이터를 통해 영상재구성하는 단계; 로 이루어지는 것을 특징으로 한다.
또한, 상기 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)는 섬광결정(120)의 집합으로 구성된 검출기모듈(110)을 중심축에 대해 상방으로 비스듬히 배치하는 것을 특징으로 한다.
또는, 상기 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)는 섬광결정(120)의 집합으로 구성된 검출기모듈(110)을 중심축에 대해 원주방향으로 비스듬히 배치한 것을 특징으로 한다.
본 발명을 첨부도면에 의해 상세히 설명하면 다음과 같다.
도 1은 섬광결정(120)의 집합으로 구성된 검출기모듈(110)을 중심축으로부터의 방사방향에 대해 비스듬히 배치한 검출기 링(100)이 설치된 양전자 방출 단층 촬영기(PET)의 평면도, 도 2는 섬광결정(120)의 집합으로 구성된 검출기모듈(110)을 중심축의 직각에 대해 비스듬히 배치한 검출기 링(100)이 설치된 양전자 방출 단층 촬영기의 측면도, 도 3은 검출기 링(100)을 구성하는 검출기모듈(110)을 도시한 개념도이다.
본 발명에 대해 구체적으로 기술하면, 본 발명 비스듬히 배치한 섬광결정(120)을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법은, 섬광결정(120)을 비스듬히 배치한 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)로 획득한 list mode 데이터를 재배열하는 단계와; 재배열한 데이터를 통해 영상재구성하는 단계; 로 이루어지는 것이다.
도 1 내지 2는 검출기 링(100)을 구성하는 검출기모듈(110)을 도시한 것이다.
상기 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)는 섬광결정(120)의 집합으로 구성된 검출기모듈(110)을 중심축으로부터의 방사방향에 대해 (예를들면, 원주방향으로) 비스듬히 배치한 것이다.
더욱 상세하게는, 도 1에 도시된 바와 같이, 상기 양전자 방출 단층 촬영기(PET)는 다수 개의 검출기모듈(110)을 중심축으로부터의 방사방향에 대해 각도
Figure pat00025
만큼 비스듬히 배치한 것이다.
이때, 상기 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)의 중심축으로부터의 방사방향으로 검출기모듈(110)을 각도
Figure pat00026
만큼 비스듬히 배치할 시, 상기 각도
Figure pat00027
는 섬광결정(120)의 크기, 검출기 링(100)의 반지름, 검출기 링(100)에 설치할 검출기모듈(110)의 수 등에 따라 결정되되, 가변적으로 결정할 수 있는 것이다.
실시예로써, 상기 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)는 16개의 검출기모듈(110)을 반경이 165mm인 검출기 링(100)에 배치하고, 상기 검출기모듈(110)은 크기가 3x3x20㎣인 섬광결정(120)을 4x4 array 형태로 집합시킬 때, 각도
Figure pat00028
는 20°로 하는 것이다.
상기 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)는 섬광결정(120)의 집합으로 구성된 검출기모듈(110)을 중심축의 직각에 대해 (상방으로) 비스듬히 배치하는 것이다.
더욱 상세하게는, 도 2에 도시된 바와 같이, 상기 양전자 방출 단층 촬영기(PET)는 검출기모듈(110)을 중심축의 직각방향에 대해 각도
Figure pat00029
만큼 비스듬히 배치한 것이다.
이때, 상기 양전자 방출 단층 촬영기(PET, Positron Emission Tomography) 중심축의 직각으로부터 검출기모듈(110)을 각도
Figure pat00030
만큼 비스듬히 배치할 시, 상기 각도
Figure pat00031
는 섬광결정(120)의 크기, 검출기 링(100) 간의 간격, 검출기 링(100)의 반지름 등에 따라 결정하되, 가변적으로 결정할 수 있는 것이다.
실시예로써, 상기 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)는 검출기모듈(110)을 반경이 165mm인 검출기 링(100)에 배치하고, 상기 검출기모듈(110)은 크기가 3x3x20㎣인 섬광결정(120)을 4x4 array 형태로 집합시킨 것이며, 검출기 링(100) 간의 간격이 2mm일 때, 각도
Figure pat00032
는 10°로 하는 것이다.
한편, 상기 양전자 방출 단층 촬영기(PET)에서 관측된 list mode 데이터는 사용하고자 하는 영상재구성 방법에 따라 3가지 데이터 포맷으로 재배열 할 수 있는 것으로, 상기 3가지 데이터 포맷은 Raw Data 포맷과, Rebinned Sinogram 포맷과, Direct Sinogram 포맷으로 이루어진다.
상기 Raw Data 포맷은 List mode 데이터의 각 관측에 대응하는 2개의 검출기의 위치인 b1,b2에 따라 데이터를 재배열하는 것으로, 상기 데이터는 ML-EM을 비롯한 다양한 반복적 영상재구성 방법에 사용될 수 있다.
상기 Rebinned Sinogram 포맷은 List mode 데이터의 각 관측에 대응하는 2개의 검출기의 위치인 b1,b2에 따라 해당 관측에 해당하는 대응선(LOR)의 높이 h와, 중심축에 대해 직각에서 벗어난 정도 d와, 상기 대응선(LOR)을 xy 평면에 투사한 직선의 각도 θ와, 중심축으로부터의 거리 u를 기준으로 데이터를 재배열하여 Yθ,u,h,d을 계산함으로써, 상기 데이터는 ML-EM을 비롯한 다양한 반복적 영상재구성 방법과 함께 여현역투사(FBP)에 사용될 수 있다.
상기 Direct Sinogram 포맷은 재구성하고자 하는 영상을 중심이 z1,z2,…,zx이고, 두께가 △z인 단층들의 2차원 영상재구성을 통해 얻으려는 경우에 사용할 수 있으며, List mode 데이터의 각 관측에 대해 대응하는 2개의 검출기의 위치인 b1,b2에 따라 해당 관측에 해당하는 대응선(LOR)이 어느 단층면에 포함되는지 계산하고. 만약 중심이 zi에 해당된다면, 해당 대응선(LOR)을 xy 평면에 투사한 직선의 각도 θ와 중심축으로부터의 거리 u와, 높이 zi를 기준으로 재배열하여
Figure pat00033
을 계산하는 것으로, 여현역투사(FBP) 또는 ML-EM을 비롯한 다양한 반복적 영상재구성 방법에 사용될 수 있다.
상기 ML-EM을 비롯한 다양한 반복적 영상재구성 방법에 필요한 시스템 행렬은 Raw Data를 사용하는 경우, 양전자 방출 단층 촬영기(PET)의 기하적 및 물리적 특성에 대한 몬테카를로(Monte Carlo) 모사방법을 이용하여 계산할 수 있는 것으로, Sinogram을 ML-EM 등의 반복적 영상재구성 방법에 적용하는 경우에는 선추적방식(Ray-Tracing Method)에 의해 시스템 행렬을 계산할 수 있는 것이며, 이 때, 시스템 행렬을 미리 계산하여 저장된 것을 사용하거나, 또는 필요할 때 마다 계산하여 사용하는 것이다.
즉, 본 발명은 섬광결정(120)을 비스듬히 배치한 양전자 방출 단층 촬영기(PET)와, 상기 양전자 방출 단층 촬영기(PET)로 획득한 데이터의 재배열 및 영상재구성 방법을 사용함으로써, 민감도를 증가시키고, Sinogram에서의 갭의 발생 가능성을 줄이고, 중심축방향의 해상도 향상 및 전반적인 영상재구성 향상을 제공하는 것이다.
도 3에 도시된 바와 같이, 상기 검출기모듈(110)은 섬광결정(120)과, 광전자 증배관(130)으로 이루어진다.
상기 섬광결정(120)은 방사선이 입사되면 섬광을 발하게 되는 것이며, 반응 깊이를 정확히 측정하기 위해 다수 개의 층으로 이루어 질 수 있다.
그리고 상기 광전자증배관(130)은 섬광결정(120)에서 발생하는 섬광을 전기 신호로 변환하는 것이다.
상기 양전자 방출 단층 촬영기(PET)에서 관측된 list mode 데이터를 사용하여 TOR 관측모델기반 반복법, LOR 관측모델기반 반복법, LOR 관측모델기반 FBP 중에 하나의 방법을 사용하여 영상재구성을 한다.
특히, 상기 list mode 데이터는 각 관측에 대해 해당 관측을 검출한 한 쌍의 검출기의 위치 정보 b1,b2와, 위치 정보가 관측된 시간 t1,t2을 포함하고 있다.
상기 TOR 관측모델기반 반복법은 Raw Data에 ML-EM 등의 반복적 영상재구성 방법을 적용하는 것으로, 시스템 행렬
Figure pat00034
의 획득, Raw Data의 재배열, 반복적 영상재구성을 순차적으로 진행한다.
이때, 상기 반복적 영상재구성에는 Raw Data
Figure pat00035
과, 재구성하고자 하는 영상 f의 화소 v에서의 값
Figure pat00036
의 선형관계를 규정하는 시스템 행렬
Figure pat00037
가 필요하며, 시스템 행렬
Figure pat00038
를 획득해야 한다.
상기 시스템 행렬
Figure pat00039
의 획득할 시에는, 상기 양전자 방출 단층 촬영기(PET)를 제작한 후, 각 화소 v로 이동시키면서 각 화소에서 N개의 양전자가 방출될 때 각각의 b에서 발생하는 감마선의 수를 측정하여 획득하거나, 또는, Monte Carlo 모사실험을 이용하여 계산하여 획득하는 것이다.
상기 Monte Carlo 모사실험에 의한 시스템 행렬 생성방법은 논문 [4] “Rafecas M, Mosler B, Dietz M, Pogl M, Stamatakis A,McElroy M P, and Ziegler S I (2004) Use of a Monte-Carlo based probability matrixfor 3D iterative reconstruction of MADPET-II data, IEEE Trans. Nucl. Sci. 51(5)2597-2605”에 기재된 바와 같이, 종래에 널리 알려진 기술임으로 자세한 기재는 생략한다.
상기 Raw Data을 재배열할 시에는, 상기 Raw Data에서 2개의 검출기의 위치 b1,b2을 인덱스로 갖는 벡터
Figure pat00040
를 모든 b에 대해
Figure pat00041
으로 초기화한 후에, list mode 데이터를 순차적으로 읽으면서 각 관측에 대응하는 2개의 검출기의 위치 b1,b2에 대해
Figure pat00042
(
Figure pat00043
의 값을 1 증가)하여 재배열되는 것이다.
상기 반복적 영상재구성을 할 시 에는, Raw Data
Figure pat00044
와 시스템 행렬
Figure pat00045
를 반복법을 이용하여 영상으로 재구성할 수 있다.
상기 영상인 f를 반복적으로 재구성하는 ML-EM는 다음의 수식을 따른다.
for
Figure pat00046
상기 LOR 관측모델기반 반복법은 LOR의 결정, 스케일 보정, Sinogram으로의 재배열을 순차적으로 진행하는 것으로, List mode 데이터를 Sinogram으로 재배열하는 과정에서는 두 검출기의 어느 두 점을 연결하는 선분을 LOR로 해야 할지 결정해야 한다.
상기 LOR의 결정을 할 시에는, 다음의 절차를 시행하여 결정하는 것이다.
먼저, Raw Data 는 2개의 검출기의 위치
Figure pat00047
에 대해서만 1이고 다른 검출기쌍에 대해서는 0인 것이며, 상기 Raw Data
Figure pat00048
에 대해 TOR(Tube-Of-Response) 관측모델을 기반으로 생성한 시스템 행렬
Figure pat00049
의 역투사를 계산하는 것으로,
Figure pat00050
로 계산하고, 상기 계산 결과는 2개의 검출기의 위치 b1,b2를 잇는 튜브(Tube)모양의 3차원 영상
Figure pat00051
이 된다.
이때, 상기 3차원 영상
Figure pat00052
에 대해, b1과 b2의 x 좌표가 다르다면 b1과 b2사이를 지나는 평면 x는 a일 때, 가능한 모든 a에 대해
Figure pat00053
의 화소 값이 가장 큰 화소
Figure pat00054
를 획득한 후, 획득한
Figure pat00055
를 직선으로 연결하여 LOR 선분
Figure pat00056
을 구할 수 있다.
반대로, 상기 3차원 영상
Figure pat00057
에 대해, b1과 b2의 x 좌표가 같다면 b1과 b2 사이를 지나는 평면 y는 a일 때, 가능한 모든 a에 대해
Figure pat00058
의 화소 값이 가장 큰 화소
Figure pat00059
를 획득한 후, 획득한
Figure pat00060
를 직선으로 연결하여 LOR 선분
Figure pat00061
을 구할 수 있다.
상기 스케일 보정을 할 시에는, 상기 LOR 관측모델에 의한 실제 데이터를 표현할 경우, TOR 관측모델에 비해 더 많은 오차를 수반한다. 상기 수반되는 오차를 보상하기 위해 먼저 List mode 데이터를 Raw Data
Figure pat00062
로 재배열한 후에 다음의 수식과 같이 스케일 보상(scale compensation)을 수행한다.
Figure pat00063
상기
Figure pat00064
는 관측된 list mode 데이터를 TOR 관측모델기반으로 재배열한 것이고, 상기
Figure pat00065
는 관측영역의 모든 곳에서 동일한 화소 값을 갖는 참조팬텀(Reference Phantom)의 list mode 데이터를 TOR 관측모델기반으로 재배열한 것이고, 상기
Figure pat00066
에 대한 참조팬텀에 대해 선분
Figure pat00067
상에서의 선적분값이며, 상기 선적분값은 Ray Tracing 방법에 의해 계산 가능하다.
상기 Rebinned Sinogram으로의 재배열을 할 시에는, List mode 데이터를 Sinogram으로 재배열하는 것이며, 스케일 보정에 따른
Figure pat00068
에 대해, 2개의 검출기의 위치 b1,b2의 중간높이 h와, 상기 중간높이 h에서 중심축을을 수직으로 통과하는 평면과 2개의 검출기를 있는 LOR(Line-Of-Response)가 이루는 각도 d와, 해당 LOR을 xy 평면에 투사한 직선의 각도 θ와 중심축으로부터의 거리 u, 즉, 4개의 인덱스를 갖는 벡터 Yθ,u,h,d로 재배열하는 것이다.
구체적으로 벡터 Yθ,u,h,d로 재배열하는 방법은, 모든 Yθ,u,h,d에 대해 Yθ,u,h,d=0로 초기화한 후에, 스케일 보상된
Figure pat00069
의 2개의 검출기의 위치 b1,b2에 대응하는 θ,u,h,d를 계산하고,
Figure pat00070
로 재배열 하는 것이며, 모든 b에 재배열을 하는 것이다.
이후, 각 θ,u,h,d에 재배열이 이루어진 횟수만큼 를 나누는 것이다.
상기 시스템 행렬
Figure pat00071
의 계산을 할 시에는, 재배열된 데이터
Figure pat00072
와 재구성하고자 하는 영상 f의 화소 v에서의 값
Figure pat00073
와의 선형관계를 규정하는 시스템 행렬
Figure pat00074
가 필요하다.
상기 시스템 행렬
Figure pat00075
는 θ,u,h,d가 규정하는 직선이 화소 v내의 길이로 계산할 수 있으며, 상기 시스템 행렬
Figure pat00076
는 미리 계산하여 저장된 것을 사용하기 보다는 필요할 때 마다 계산하여 사용하는 방식을 이용한다.
상기 반복적 영상재구성을 할 시 에는, 재배열된 데이터
Figure pat00077
와 시스템 행렬
Figure pat00078
에 반복법을 적용하여 영상을 재구성할 수 있다.
상기 영상 f를 반복적으로 재구성하는 실시예 ML-EM는 다음의 수식을 따른다.
for
Figure pat00079
상기 LOR 관측모델기반 FBP은 Direct Sinogram으로의 재배열, FBP를 이용한 영상재구성을 순차적으로 진행하는 것이다.
상기 Direct Sinogram으로의 재배열을 할 시에는, 관측된 list mode 데이터를 Raw Data로의 재배열, 스케일 보정, Sinogram으로의 재배열의 과정을 순차적으로 적용하여 Sinogarm
Figure pat00080
으로 재배열하는 것이다.
특히, 상기 재배열하고자 하는 영상을 중심이 z1,z2,…,zn이며 두께가 △z인 단층들로 이루어진 2차원 영상재구성을 통해 얻으려는 경우, 계산된 Sinogarm
Figure pat00081
에 대해서 주어진 zi 단층에 속하는 θ,u,h,d에 해당하는
Figure pat00082
를 모아서 Direct Sinogram
Figure pat00083
로 재배열하는 것이다.
상기 FBP를 이용한 영상재구성을 할 시에는, 재배열된 Direct Sinogram
Figure pat00084
에 FBP를 적용하여 zi-단층에서의 영상
Figure pat00085
를 재구성하여 전체영상을 재구성한다.
상기 영상
Figure pat00086
는 아래의 수식을 통해 계산할 수 있다.
Figure pat00087
상기
Figure pat00088
는 높이 h에서의 f의 단층영상이며, FFT는 Fast Fourier Transform이며, IFFT는 Inverse Fast Fourier Transform이다.
본 발명의 구체적인 구성 요소 등과 같은 특정 사항들과 한정된 실시예 및 도면에 의해 설명하였다. 그러나, 이는 본 발명의 전반적인 이해를 돕기 위해서 제공된 것일 뿐, 본 발명은 상기의 실시예에 한정되는 것은 아니다. 본 발명이 속하는 분야에서 통상적인 지식을 바탕으로 이러한 기재로부터 다양한 수정 및 변형이 가능하므로, 본 발명의 사상은 설명된 실시예에 국한되어 정해져서는 아니되며, 후술하는 특허 청구범위뿐 아니라 이 특허 청구범위와 균등하거나 등가적 변형이 있는 모든 것들은 본 발명 사상의 범주에 속한다고 할 것이다.
도 1과 도 2에서 유추해 볼 수 있듯이 비스듬히 배치된 검출기모듈(110)로 인하여 PET 내부의 관측영역에서 방출된 감마선이 검출될 가능성이 높아진다. 즉, 민감도를 높일 수 있다.
그리고 섬광결정(120)이 다수의 층인 경우에 비스듬히 배치된 검출기모듈(110)은 Rebinned Sinogram 또는 Direct Sinogram으로의 재배열과정에서 해당 영역에 대한 재배열데이터가 하나도 없는 경우인 gap이 발생할 가능성을 줄여준다.
그리고 섬광결정(120)이 한 층인 경우에도 비스듬히 배치된 검출기모듈(110)은 TOR 관측모델기반 반복적 영상재구성을 통해 종축, 횡축 방향으로 향상된 해상도를 제공한다.
그러므로 본 발명은 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법을 통하여, 섬광결정을 비스듬히 배치한 양전자 방출 단층 촬영기(PET)와, 상기 양전자 방출 단층 촬영기(PET)로 획득한 데이터의 재배열 및 영상재구성 방법을 사용함으로써, 민감도를 증가시키고, Sinogram에서의 갭의 발생 가능성을 줄이고, 중심축방향의 해상도 향상 및 전반적인 영상재구성 향상을 제공하는 현저한 효과가 있다.
100 : 검출기 링 110 : 검출기모듈
120 : 섬광결정 130 : 광전자 증배관

Claims (3)

  1. 섬광결정(120)을 비스듬히 배치한 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)로 획득한 list mode 데이터를 재배열하는 단계와; 재배열한 데이터를 통해 영상재구성하는 단계; 로 이루어지는 것을 특징으로 하는 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법
  2. 제 1항에 있어서, 상기 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)는 섬광결정(120)의 집합으로 구성된 검출기모듈(110)을 중심축에 대해 상방으로 비스듬히 배치하는 것을 특징으로 하는 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법
  3. 제 1항에 있어서, 상기 양전자 방출 단층 촬영기(PET, Positron Emission Tomography)는 섬광결정(120)의 집합으로 구성된 검출기모듈(110)을 중심축에 대해 원주방향으로 비스듬히 배치한 것을 특징으로 하는 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법
KR1020170097322A 2017-07-31 2017-07-31 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법 KR101958098B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020170097322A KR101958098B1 (ko) 2017-07-31 2017-07-31 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020170097322A KR101958098B1 (ko) 2017-07-31 2017-07-31 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법

Publications (2)

Publication Number Publication Date
KR20190013161A true KR20190013161A (ko) 2019-02-11
KR101958098B1 KR101958098B1 (ko) 2019-03-13

Family

ID=65370055

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020170097322A KR101958098B1 (ko) 2017-07-31 2017-07-31 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법

Country Status (1)

Country Link
KR (1) KR101958098B1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20210110236A (ko) * 2020-02-28 2021-09-07 가천대학교 산학협력단 Pet-mri 다중 영상 정합 시스템 및 방법

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20070041754A (ko) * 2004-08-20 2007-04-19 시로나 덴탈 시스템스 게엠베하 기울어진 기하학적 구조를 이용한 3ｄ 재구성
KR20090126850A (ko) * 2008-06-05 2009-12-09 연세대학교 산학협력단 반응깊이측정을 위한 양전자방출단층촬영기의 준블록배열형 검출기모듈 및 그 검출방법
JP2010101666A (ja) * 2008-10-22 2010-05-06 Hitachi Medical Corp 核医学診断装置
KR20120077836A (ko) * 2010-12-31 2012-07-10 중원대학교 산학협력단 이중 검출기 구조를 이용한 방사선 검출기 및 이를 이용한 방사선 영상장치
KR20120122158A (ko) * 2011-04-28 2012-11-07 서강대학교산학협력단 한층 섬광결정을 이용하여 영상을 재구성하는 방법
KR20160001253A (ko) * 2014-06-27 2016-01-06 동서대학교산학협력단 양방향 방사선 검출기
CN106821402A (zh) * 2016-12-14 2017-06-13 赛诺联合医疗科技（北京）有限公司 构建pet图像的方法和装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20070041754A (ko) * 2004-08-20 2007-04-19 시로나 덴탈 시스템스 게엠베하 기울어진 기하학적 구조를 이용한 3ｄ 재구성
KR20090126850A (ko) * 2008-06-05 2009-12-09 연세대학교 산학협력단 반응깊이측정을 위한 양전자방출단층촬영기의 준블록배열형 검출기모듈 및 그 검출방법
JP2010101666A (ja) * 2008-10-22 2010-05-06 Hitachi Medical Corp 核医学診断装置
KR20120077836A (ko) * 2010-12-31 2012-07-10 중원대학교 산학협력단 이중 검출기 구조를 이용한 방사선 검출기 및 이를 이용한 방사선 영상장치
KR20120122158A (ko) * 2011-04-28 2012-11-07 서강대학교산학협력단 한층 섬광결정을 이용하여 영상을 재구성하는 방법
KR20160001253A (ko) * 2014-06-27 2016-01-06 동서대학교산학협력단 양방향 방사선 검출기
CN106821402A (zh) * 2016-12-14 2017-06-13 赛诺联合医疗科技（北京）有限公司 构建pet图像的方法和装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20210110236A (ko) * 2020-02-28 2021-09-07 가천대학교 산학협력단 Pet-mri 다중 영상 정합 시스템 및 방법

Also Published As

Publication number Publication date
KR101958098B1 (ko) 2019-03-13

Similar Documents

Publication Publication Date Title
Tong et al. Image reconstruction for PET/CT scanners: past achievements and future challenges
EP3067864B1 (en) Iterative reconstruction with enhanced noise control filtering
JP6440230B2 (ja) マルチモダリティ撮像のシステム及び方法
Gaudin et al. Performance simulation of an ultrahigh resolution brain PET scanner using 1.2-mm pixel detectors
Pratx et al. Online detector response calculations for high-resolution PET image reconstruction
WO2019000389A1 (en) SYSTEM AND METHOD FOR POSITRON EMISSION TOMOGRAPHY
US20100067758A1 (en) Point Spread Function Radial Component Implementation in Joseph's Forward Projector
Zhang et al. Development and evaluation of a LOR-based image reconstruction with 3D system response modeling for a PET insert with dual-layer offset crystal design
US10482634B2 (en) Systems and methods for imaging with anisotropic voxels
WO2012120791A1 (ja) ポジトロンｃｔ装置および再構成方法
Kang et al. Initial results of a mouse brain PET insert with a staggered 3-layer DOI detector
US9495770B2 (en) Practical model based CT construction
CN103908278A (zh) 图像重建方法及装置、医疗成像系统
KR101233984B1 (ko) 한층 섬광결정을 이용하여 영상을 재구성하는 방법
KR101958098B1 (ko) 비스듬히 배치한 섬광결정을 이용한 양전자 방출 단층 촬영기의 영상재구성 방법
Kao et al. Accurate image reconstruction using DOI information and its implications for the development of compact PET systems
Li et al. LOR-interleaving image reconstruction for PET imaging with fractional-crystal collimation
KR101493683B1 (ko) 콘-빔 기반 반응선 재구성을 이용한 초고해상도 ｐｅｔ 영상 재구성 장치 및 그 방법
KR102283454B1 (ko) 양전자 방출 단층 촬영기의 비행시간 정보를 이용하여 관측된 데이터에 포함된 반응깊이 정보를 보정하여 영상을 재구성하는 영상재구성방법
JP6052425B2 (ja) 輪郭画像生成装置および核医学診断装置
Lee et al. Precise system models using crystal penetration error compensation for iterative image reconstruction of preclinical quad-head PET
US7890282B2 (en) Estimation of crystal efficiency with axially compressed sinogram
JP2011002306A (ja) Ｐｅｔ装置の逐次近似画像再構成法
Bai et al. PET image reconstruction: methodology and quantitative accuracy
Protonotarios et al. Automatic cumulative sums contour detection of FBP-reconstructed multi-object nuclear medicine images

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