KR101092648B1 - Pet 검출기에서의 감마선 펄스 신호의 도달 시간 검출방법 - Google Patents
Pet 검출기에서의 감마선 펄스 신호의 도달 시간 검출방법 Download PDFInfo
- Publication number
- KR101092648B1 KR101092648B1 KR1020100039492A KR20100039492A KR101092648B1 KR 101092648 B1 KR101092648 B1 KR 101092648B1 KR 1020100039492 A KR1020100039492 A KR 1020100039492A KR 20100039492 A KR20100039492 A KR 20100039492A KR 101092648 B1 KR101092648 B1 KR 101092648B1
- Authority
- KR
- South Korea
- Prior art keywords
- signal
- pulse signal
- ray pulse
- detecting
- digital
- Prior art date
Links
- 230000005251 gamma ray Effects 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000000630 rising effect Effects 0.000 claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 20
- 230000003287 optical effect Effects 0.000 claims abstract description 13
- 239000013078 crystal Substances 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 9
- 230000003111 delayed effect Effects 0.000 claims description 5
- 238000002600 positron emission tomography Methods 0.000 description 52
- 238000001514 detection method Methods 0.000 description 6
- 238000002595 magnetic resonance imaging Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 102000006602 glyceraldehyde-3-phosphate dehydrogenase Human genes 0.000 description 4
- 108020004445 glyceraldehyde-3-phosphate dehydrogenase Proteins 0.000 description 4
- -1 Lutetium Aluminum Chemical compound 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- MARUHZGHZWCEQU-UHFFFAOYSA-N 5-phenyl-2h-tetrazole Chemical compound C1=CC=CC=C1C1=NNN=N1 MARUHZGHZWCEQU-UHFFFAOYSA-N 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- QWUZMTJBRUASOW-UHFFFAOYSA-N cadmium tellanylidenezinc Chemical compound [Zn].[Cd].[Te] QWUZMTJBRUASOW-UHFFFAOYSA-N 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- NZOCXFRGADJTKP-UHFFFAOYSA-K lutetium(3+);triiodide Chemical compound I[Lu](I)I NZOCXFRGADJTKP-UHFFFAOYSA-K 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 229910052688 Gadolinium Inorganic materials 0.000 description 1
- 229910014323 Lanthanum(III) bromide Inorganic materials 0.000 description 1
- 229910052765 Lutetium Inorganic materials 0.000 description 1
- OSGMVZPLTVJAFX-UHFFFAOYSA-N [Gd].[Lu] Chemical compound [Gd].[Lu] OSGMVZPLTVJAFX-UHFFFAOYSA-N 0.000 description 1
- ANDNPYOOQLLLIU-UHFFFAOYSA-N [Y].[Lu] Chemical compound [Y].[Lu] ANDNPYOOQLLLIU-UHFFFAOYSA-N 0.000 description 1
- 238000012984 biological imaging Methods 0.000 description 1
- 229910052797 bismuth Inorganic materials 0.000 description 1
- JCXGWMGPZLAOME-UHFFFAOYSA-N bismuth atom Chemical compound [Bi] JCXGWMGPZLAOME-UHFFFAOYSA-N 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- UIWYJDYFSGRHKR-UHFFFAOYSA-N gadolinium atom Chemical compound [Gd] UIWYJDYFSGRHKR-UHFFFAOYSA-N 0.000 description 1
- 239000002223 garnet Substances 0.000 description 1
- 229910052746 lanthanum Inorganic materials 0.000 description 1
- FZLIPJUXYLNCLC-UHFFFAOYSA-N lanthanum atom Chemical compound [La] FZLIPJUXYLNCLC-UHFFFAOYSA-N 0.000 description 1
- XKUYOJZZLGFZTC-UHFFFAOYSA-K lanthanum(iii) bromide Chemical compound Br[La](Br)Br XKUYOJZZLGFZTC-UHFFFAOYSA-K 0.000 description 1
- OHSVLFRHMCKCQY-UHFFFAOYSA-N lutetium atom Chemical compound [Lu] OHSVLFRHMCKCQY-UHFFFAOYSA-N 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- NJPPVKZQTLUDBO-UHFFFAOYSA-N novaluron Chemical compound C1=C(Cl)C(OC(F)(F)C(OC(F)(F)F)F)=CC=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F NJPPVKZQTLUDBO-UHFFFAOYSA-N 0.000 description 1
- 230000008506 pathogenesis Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004393 prognosis Methods 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000037152 sensory function Effects 0.000 description 1
- 238000007619 statistical method Methods 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
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Power Engineering (AREA)
- Optics & Photonics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
PET 검출기의 섬광결정에서 검출한 섬광을 전기적인 신호로 변환한 광센서로부터 전기적인 신호를 입력받아 아날로그-디지털 변환기(ADC)를 통해 변환된 감마선 펄스 신호에 대한 디지털 신호로부터 감마선 펄스 신호의 도달 시간을 검출하는 방법을 개시한다. 상기 방법은 (a) 소정의 시간 간격을 가진 상기 감마선 펄스 신호에 대한 디지털 신호의 전압값의 초기 상승선을 검출하는 단계, (b) 상기 초기 상승선 이전에 획득된 상기 소정의 시간 간격을 가진 상기 디지털 신호의 통계적 처리를 통하여 기준선을 검출하는 단계, (c) 상기 초기 상승선 및 상기 기준선과의 교차점을 계산하는 단계, 및 (d) 상기 교차점을 이용하여 상기 감마선 펄스 신호의 도달 시간을 검출하는 단계를 포함한다. PET 응용에 있어서 PET 검출기를 위한 간단하고 향상된 시간 분해능을 나타낸다.
Description
본 발명은 PET 검출기에서의 신호처리방법에 관한 것으로, 특히 PET 검출기에서의 감마선 펄스 신호의 도달 시간 검출방법에 관한 것이다.
세포, 전임상, 임상실험 및 환자의 진단을 위해 사용되는 의료 영상은 일반적으로 크게 구조적인 영상과 기능적인 영상으로 분류된다. 구조적인 영상은 인체의 구조 및 해부학적 영상을 의미하고, 기능적인 영상은 인체의 인지, 감각기능 등에 대한 기능 정보를 직접 또는 간접적인 방법으로 영상화하는 것이다. 구조 및 해부학적 영상 기술에는 컴퓨터단층촬영장치(Computed Tomography: CT), 자기공명단층촬영장치(Magnetic Resonance Imaging: MRI)등이 있고, 인체의 생리적 및 생화학적 작용을 관찰하여 기능 정보를 영상화하는 기술로서는 양전자방출단층촬영장치(Positron Emission Tomography :PET)가 널리 사용되고 있다.
PET은 비침습적으로 인체 기능을 계량화하는 강력한 생물학적 영상도구로서, 방사성 활성을 갖는 양전자 방출동위원소로 표지 된 생물학적 탐지자(probe) 분자를 체내에 주입한 후, 방사능의 분포를 단층촬영으로 재구성하여 영상화하여 인체의 각 장기 내의 생리적, 생화학적인 반응을 정량화할 수 있다. PET에 의해 제공되는 뇌, 장기 등의 인체 구조에 대한 기능적/분자학적 정보는 질병의 병인 연구, 진단 예후 판정 및 항암 치료 후 경과 관찰 등에 유용하게 이용할 수 있다.
도 1a는 PET 시스템을 도시한 사시도이고 도 1b는 PET 시스템 중 PET 검출기를 확대한 사시도이다.
도 1a 및 1b를 참조하면, PET 시스템(100)은 PET 검출기(120), PET 검출기(120)에 일측이 연결되는 케이블(130)과 케이블(130)의 타측에 연결되는 PET 회로부(140)를 포함한다. 환자를 이송하는 베드(111)는 받침대(113) 위에 이송가능하게 형성되고, PET 검출기(120)는 환자가 통과할 수 있도록 소정의 내경을 갖는 관 형상으로 이루어져 촬영구(도시되지 않음) 내부에 형성된다. 일측에는 다수개의 섬광결정(121)이 환형으로 배열되어 섬광결정 배열(122)를 이루고, 섬광결정 배열(122)는 PET 검출기(120)의 길이방향을 따라 다수개 배열되고, 섬광결정(121) 타측부는 광센서와 연결되어 PET 검출기(120)를 구성한다. PET 검출기(120)를 구성하는 광센서(도시되지 않음)는 섬광결정(121)에서 검출한 섬광을 전기적인 신호로 변환한다. 광센서로는 PMT(Photo-Multiplier Tube), PIN Diode(Positive-Intrinsic-Negative Diode), CdTe(Cadmium Telluride), CZT(Cadmium Zinc Telluride), APD(Avalanche Photo Diode) 또는 GAPD(Geiger mode Avalanche Photo Diode) 등이 이용될 수 있다.
PET 영상을 생성하기 위해서는 상기 광센서에서 출력되는 신호의 동시성을 측정해야하고, 동시성을 측정하기 위해서는 측정된 감마선의 도달 시간을 검출해야 한다. 이러한 PET 검출기에서의 감마선 펄스 신호의 도달 시간 검출방법으로는 디지털 상수 분수 차별법(Digital constant fraction discrimination: DCFD) 및 최대 상승 보간법(Maximum rise interpolation: MRI) 등이 있었다.
디지털 상수 분수 차별법(DCFD)은 아날로그 상수 분수 차별법의 디지털버전이며, 이는 다른 디지털 시간 측정법에 비해 더 나은 시간 분해능을 제공할 수 있다. 디지털 상수 분수 차별법은 디지털화된 감마선 펄스 신호는 지연되고 증폭되고나서 전도되며 그리고 나서 본래 펄스에 추가된다. 따라서, 단극 펄스는 양극 펄스로 변환된다. 펄스 도달 시간은 두 개의 제로-크로싱 샘플들 사이에 선형 보간법에 의해 계산된다. 이 방법의 정확성은 노이즈와 신호 펄스 영향에 의해 영향 받는다.
최대 상승 보간법(Maximum rise interpolation: MRI)은 기준선과 최대 상승선의 교차점을 계산하여 펄스 도달 시간을 검출한다. 최대 상승 보간법은 간단하고 빠르지만, 이 방법의 정확성은 ADC 샘플링 단계에 의해 영향받는다. 이는 최대 상승선의 기울기는 ADC(Analog-digital converter) 클럭의 샘플링 단계에 의해 다양해지기 때문이다. 통계적인 방법으로 성능을 향상시킬 수 있으나, 다량의 데이터가 후 처리(post-process)되어야 하는 단점이 있다. 따라서, PET 검출기에서 개선된 시간 분해능을 가지는 감마선 펄스 신호의 도달 시간 검출방법, 즉 디지털 시간 검출방법이 요구되는 실정이다.
본 발명이 이루고자 하는 기술적 과제는 PET 검출기에서 향상된 시간 분해능을 가지는 감마선 펄스 신호의 도달 시간 검출방법 및 그의 회로모듈을 제공하는 데 있다.
본 발명의 일 양태에 따르면, PET 검출기의 섬광결정에서 검출한 섬광을 전기적인 신호로 변환한 광센서로부터 전기적인 신호를 입력받아 아날로그-디지털 변환기(ADC)를 통해 변환된 감마선 펄스 신호에 대한 디지털 신호로부터 감마선 펄스 신호의 도달 시간을 검출하는 방법에 있어서, (a) 소정의 시간 간격을 가진 상기 감마선 펄스 신호에 대한 디지털 신호의 전압값의 초기 상승선을 검출하는 단계, (b) 상기 초기 상승선 이전에 획득된 상기 소정의 시간 간격을 가진 상기 디지털 신호의 통계적 처리를 통하여 기준선을 검출하는 단계, (c) 상기 초기 상승선 및 상기 기준선과의 교차점을 계산하는 단계, 및 (d) 상기 교차점을 이용하여 상기 감마선 펄스 신호의 도달 시간을 검출하는 단계를 포함하는 감마선 펄스 신호의 도달 시간 검출방법을 제공한다.
상기 (a) 단계는 (a-1) 최초의 상기 감마선 펄스 신호에서 상기 소정의 시간 간격을 가진 펄스 신호를 빼서 양의 피크 및 음의 피크를 가진 양극 펄스를 획득하는 단계, (a-2) 상기 양의 피크를 검출하는 단계,및 (a-3) 상기 초기 상승선을 위한 트리거링을 생성하는 단계를 포함할 수 있다.
상기 소정의 시간 간격은 상기 아날로그-디지털 변환기(ADC)의 샘플링 주파수에 따라 변화될 수 있다.
상기 감마선 펄스 신호의 도달 시간 검출방법은 상기 (a) 단계 이전에
상기 감마선 펄스 신호에 대한 디지털 신호를 증폭시키는 단계를 더 포함할 수 있다.
본 발명의 다른 양태에 따르면, 상기 감마선 펄스 신호의 도달 시간 검출방법을 수행하기 위해 디지털 처리 장치에 의해 실행될 수 있는 명령어들의 프로그램이 구현되어 있으며, 상기 디지털 처리 장치에 의해 판독될 수 있는 프로그램이 기록된 기록매체를 제공한다.
본 발명의 또 다른 양태에 따르면, PET 검출기의 섬광결정에서 검출한 섬광을 전기적인 신호로 변환한 광센서로부터 전기적인 신호를 입력받아 샘플링을 통하여 감마선 펄스 신호에 대한 디지털 신호를 출력하는 아날로그-디지털 변환부,및 상기 디지털 신호를 입력 받아 소정의 시간 간격을 가진 상기 디지털 신호의 전압값의 초기 상승선을 검출하고 상기 초기 상승선 이전에 획득된 상기 소정의 시간 간격을 가진 상기 디지털 신호의 통계적 처리를 통하여 기준선을 검출하여 상기 초기 상승선 및 상기 기준선과의 교차점으로부터 상기 감마선 펄스 신호의 도달 시간을 검출하는 신호처리부를 포함하는 PET 신호처리모듈을 제공한다.
본 실시예에 따른 PET 검출기에서의 감마선 펄스 신호의 도달 시간 검출방법은 PET 검출기에서 더 정확하게 감마선 펄스 신호의 도달 시간을 검출할 수 있다. PET 응용에 있어서 PET 검출기를 위한 간단하고 향상된 시간 분해능을 나타낸다. 본 실시예에 따른 PET 신호처리모듈이 적용된 PET 검출기는 향상된 성능의 PET 단층 영상을 나타낸다.
도 1a는 PET 시스템을 도시한 사시도이고 도 1b는 PET 시스템 중 PET 검출기를 확대한 사시도이다.
도 2는 PET 검출기의 신호를 보여주는 그래프이다.
도 3은 PET 검출기의 신호를 보여주는 그래프를 확대한 도면이다.
도 4는 PET 검출기의 감마선 펄스 신호를 보여주는 그래프이다.
도 5는 본 실시예에 따른 PET 검출기에서의 펄스 도달 시간을 획득하기 위한 방법을 설명하기 위한 그래프이다.
도 6은 본 실시예에 따른 감마선 펄스 신호의 도달 시간 검출방법을 설명하는 흐름도이다.
도 7은 본 실시예에 따른 초기 상승선을 획득하는 방법을 설명하는 흐름도이다.
도 8은 본 실시예에 따른 초기 상승선을 획득하는 방법을 설명하는 도면이다.
도 9는 본 실시예에 따라 기준선을 검출하기 위한 신호처리 단계의 일 예를 보여주는 도면이다.
도 10은 본 실시예에 따른 PET 신호처리모듈을 보여주는 블록도이다.
도 11은 본 실시예과 종래방법에 따른 시간 스펙트럼을 보여주는 그래프이다.
도 12는 본 실시예에 따른 PET 이미지를 보여주는 도면이다.
도 13은 도 12의 이미지의 프로파일을 보여주는 도면이다.
도 2는 PET 검출기의 신호를 보여주는 그래프이다.
도 3은 PET 검출기의 신호를 보여주는 그래프를 확대한 도면이다.
도 4는 PET 검출기의 감마선 펄스 신호를 보여주는 그래프이다.
도 5는 본 실시예에 따른 PET 검출기에서의 펄스 도달 시간을 획득하기 위한 방법을 설명하기 위한 그래프이다.
도 6은 본 실시예에 따른 감마선 펄스 신호의 도달 시간 검출방법을 설명하는 흐름도이다.
도 7은 본 실시예에 따른 초기 상승선을 획득하는 방법을 설명하는 흐름도이다.
도 8은 본 실시예에 따른 초기 상승선을 획득하는 방법을 설명하는 도면이다.
도 9는 본 실시예에 따라 기준선을 검출하기 위한 신호처리 단계의 일 예를 보여주는 도면이다.
도 10은 본 실시예에 따른 PET 신호처리모듈을 보여주는 블록도이다.
도 11은 본 실시예과 종래방법에 따른 시간 스펙트럼을 보여주는 그래프이다.
도 12는 본 실시예에 따른 PET 이미지를 보여주는 도면이다.
도 13은 도 12의 이미지의 프로파일을 보여주는 도면이다.
본 발명은 다양한 변경을 가할 수 있고 여러 가지 실시예를 가질 수 있는 바, 특정 실시예들을 도면에 예시하고 설명하고자 한다. 그러나, 이는 본 발명을 특정한 실시 형태에 대해 한정하려는 것이 아니며, 본 발명의 사상 및 기술 범위에 포함되는 모든 변경, 균등물 내지 대체물을 포함하는 것으로 이해되어야 한다.
어떤 구성요소가 다른 구성요소에 "연결되어" 있다거나 "접속되어" 있다고 언급된 때에는, 그 다른 구성요소에 직접적으로 연결되어 있거나 또는 접속되어 있을 수도 있지만, 중간에 다른 구성요소가 존재할 수도 있다고 이해되어야 할 것이다. 반면에, 어떤 구성요소가 다른 구성요소에 "직접 연결되어" 있다거나 "직접 접속되어" 있다고 언급된 때에는, 중간에 다른 구성요소가 존재하지 않는 것으로 이해되어야 할 것이다.
본 출원에서 사용한 용어는 단지 특정한 실시예를 설명하기 위해 사용된 것으로, 본 발명을 한정하려는 의도가 아니다. 단수의 표현은 문맥상 명백하게 다르게 뜻하지 않는 한, 복수의 표현을 포함한다. 본 출원에서, "포함하다" 또는 "가지다" 등의 용어는 명세서상에 기재된 특징, 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것이 존재함을 지정하려는 것이지, 하나 또는 그 이상의 다른 특징들이나 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것들의 존재 또는 부가 가능성을 미리 배제하지 않는 것으로 이해되어야 한다.
다르게 정의되지 않는 한, 기술적이거나 과학적인 용어를 포함해서 여기서 사용되는 모든 용어들은 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자에 의해 일반적으로 이해되는 것과 동일한 의미를 가지고 있다. 일반적으로 사용되는 사전에 정의되어 있는 것과 같은 용어들은 관련 기술의 문맥 상 가지는 의미와 일치하는 의미를 가지는 것으로 해석되어야 하며, 본 출원에서 명백하게 정의하지 않는 한, 이상적이거나 과도하게 형식적인 의미로 해석되지 않는다.
이하, 첨부된 도면을 참조하여 실시예를 상세히 설명하되, 도면 부호에 관계없이 동일하거나 대응하는 구성 요소는 동일한 참조 번호를 부여하고 이에 대한 중복되는 설명은 생략하기로 한다.
본 실시예에 따른 감마선 펄스 신호의 도달 시간 검출방법은 PET 검출기의 섬광결정에서 검출한 섬광을 전기적인 신호로 변환한 광센서로부터 전기적인 신호를 입력받아 아날로그-디지털 변환기(ADC)를 통해 변환된 감마선 펄스 신호에 대한 디지털 신호로부터 감마선 펄스 신호의 도달 시간을 검출한다.
본 실시예에 따른 PET 검출기에서의 감마선 펄스 신호의 도달 시간 검출방법에서 샘플링된 감마선 펄스 신호는 기준선(baseline)으로부터 초기 상승률로 전압이 상승하다가 최대 상승률로 전압이 상승하는 감마신호의 중요한 특성을 이용한다. 여기서, 초기 상승률(initial rise)은 최대 상승률(maximum rise) 보다 더 기울기가 작다. 또한, 초기 상승률(initial rise)과 최대 상승률(maximum rise) 사이에 예비 최대 상승률(pre-maximum rise)이 있을 수 있다. 이러한 예비 최대 상승률(pre-maximum rise)은 초기 상승률(initial rise) 보다 기울기가 크고 최대 상승률(maximum rise) 보다 기울기가 작을 수 있다.
본 실시예에서 펄스 도달 시간은 각 감마선 펄스 신호의 기준선과 초기 상승률을 가지는 초기 상승선의 교차점을 이용함으로써 획득할 수 있다. 또한, 본 실시예는 BGO(Bismuth Germanate), LSO(Lutetium Oxyorthosilicate), LYSO(Lutetium Yttrium Oxyorthosilicate), LuAP(Lutetium Aluminum Perovskite), LuYAP(Lutetium Yttrium Aluminum Perovskite), LaBr3 (Lanthanum Bromide) 및 LuI3(Lutetium Iodide), GSO(Gadolinium oxyorthosilicate), LGSO (lutetium gadolinium oxyorthosilicate), LuAG (Lutetium aluminum garnet) 등 중 적어도 하나를 이용한 광센서로 구성된 PET 검출기에서 이용할 수 있다.
도 2는 PET 검출기의 신호를 보여주는 그래프이고 도 3은 도 2를 확대한 그래프이다. 도 2 및 도 3을 참조하면, 시간에 따라 감마선 신호가 펄스 신호 형태로 출력되는 것을 확인할 수 있다.
도 4는 PET 검출기의 감마선 펄스 신호를 보여주는 그래프이고 도 5는 본 실시예에 따른 PET 검출기에서의 펄스 도달 시간을 획득하기 위한 방법을 설명하기 위한 그래프이다.
도 4 및 5를 함께 참조하면, 도 5는 도 4의 초기 50 ns부분을 확대해 본 그래프이다. 또한, 도 4는 도 3에서 감마선 펄스 신호(300)을 확대해 본 것이다. 도 5에서 포인트(a-2), 포인트(a-1), 포인트(a) 및 포인트(a+1)은 샘플링되어 디지털화된 감마선 펄스 신호의 각각 10 ns 이후의 포인트를 나타낸다. 포인트(a-2), 포인트(a-1), 포인트(a) 및 포인트(a+1)의 시간 간격은 ADC의 주파수에 따라 변경될 수 있다.여기서, 이러한 포인트들의 수는 다양할 수 있다. 예를 들어, 초기 상승선 전의 포인트로써 포인트(a-2) 외에 더 많은 포인트들이 존재할 수도 있다. V(a-2), V(a-1), V(a), V(a+1)은 각각의 포인트에서의 전압을 나타낸다. 본 실시예에서 사용한 구체적인 수치는 본 발명의 일 실시예를 설명하기 위한 것에 불과하므로 본 발명의 내용이 이러한 구체적인 수치에 한정되지 않는다. 따라서, 상기 10 ns 또는 50 ns는 다른 수치로 변경되어 적용될 수 있음은 당업자에게 자명할 것이다.
도 6은 본 실시예에 따른 감마선 펄스 신호의 도달 시간 검출방법을 설명하는 흐름도이다. 도 5 및 도 6을 함께 참조하면, PET 검출기에서의 감마선 펄스 신호를 샘플링을 통하여 디지털화한 다음, 소정의 시간 간격을 가진 감마선 펄스 신호에 대한 디지털 신호의 전압값의 초기 상승선을 검출한다(S610). 이러한 S610 이전에 감마선 펄스 신호에 대한 디지털 신호는 증폭될 수 있다. S610은 이하에서 자세히 설명한다.
도 7은 본 실시예에 따른 초기 상승선을 획득하는 방법을 설명하는 흐름도이고, 도 8은 본 실시예에 따른 초기 상승선을 획득하는 방법을 설명하는 도면이다.
도 7 및 도 8을 함께 참조하면, 먼저, 양극 펄스(bipolar pulse)를 생성한다(S611). 구체적으로는, 최초 펄스 신호로부터 10 ns 지연된 펄스 신호를 빼서 양의 피크 및 음의 피크를 가진 양극 펄스를 획득할 수 있다.
그런 다음, 양의 피크를 검출한다(S612). 구체적으로는, 최대 상승선(550)의 이차점 V(a+1)에 해당하는 시간점이 피크 센싱(peak sensing)에 의해 검출될 수 있다. 다음으로, 초기 상승선을 위한 트리거를 생성한다(S613). 상기 시간점을 이용하여 10 ns 지연된 펄스 및 20 ns 지연된 펄스의 트리거링(triggering)이 초기 상승선의 각 V(a) 및 V(a-1)를 저장할 수 있다.
다시 도 6을 참조하면, 기준선을 검출한다(S620).
도 9는 본 실시예에 따라 기준선을 검출하기 위한 신호처리 단계의 일 예를 보여주는 도면이다.
도 9를 참조하면, 초기 상승선 검출(S610)후, V(a-1) 및 V(a)가 감마선 펄스 신호의 도달 시간 계산을 위하여 저장될 수 있다. 초기 상승선 전에 획득된 적어도 한 개의 샘플값에 대한 평균값을 계산하여 각 펄스의 기준선 값을 검출한다.
다시 도 6을 참조하면, 각 펄스에 대해 초기 상승선(510)과 기준선(530)과의 교차점을 계산한다(S630). 이러한 초기 상승선(510)의 기울기는 최대 상승선(550)의 기울기에 비해 ADC 샘플링 단계에 의해 덜 영향을 받는다. 그런 다음, 펄스 도달 시간을 검출한다(S640). 본 실시예에 따른 펄스 도달 시간(TIRI)은 다음 수학식 1에 의해 획득될 수 있다.
여기서, Ta -1은 포인트(a-1)에서의 시간 및 △t1는 Ta -1과 초기 상승선(510)과 기준선(530)과의 교차점의 시간차를 나타내고, △t1은 다음 수학식 2에 의해 획득될 수 있다.
또한, 상기 △t1은 다음 수학식 3에 의해 획득될 수도 있다.
도 10은 본 실시예에 따른 PET 신호처리모듈을 보여주는 블록도이다.
도 10을 참조하면, PET 신호처리모듈(800)은 아날로그-디지털 변환부(ADC)(810), 신호처리부(820), 제어부(830) 및 저장부(840)를 포함한다.
아날로그-디지털 변환부(810)는 PET 검출기의 섬광결정에서 검출한 섬광을 전기적인 신호로 변환한 광센서로부터 전기적인 신호를 입력받아 샘플링을 통하여 감마선 펄스 신호에 대한 디지털 신호를 출력한다.
신호처리부(820)는 감마선 펄스 신호에 대한 디지털 신호를 입력받아 본 실시예에 따라 초기 상승선 및 기준선을 검출하고 초기 상승선 및 기준선의 교차점을 이용하여 펄스 도달 시간을 검출한다. 저장부(840)는 아날로그-디지털 변환부(810)로부터 감마선 펄스 신호에 대한 디지털 신호 또는 신호처리부(820)에서의 처리되는 데이터를 저장한다. 제어부(830)는 본 실시예에 따라 아날로그-디지털 변환부(810) 및 신호처리부(820)의 입력 및 출력을 제어할 수 있다.
이하에서는, 본 실시예의 효과에 대하여 설명한다.
표 1은 본 실시예과 종래방법(각각 30% 에너지 윈도우)에 따른 시간 분해능을 비교한 것이고 도 11은 본 실시예과 종래방법에 따른 시간 스펙트럼을 보여주는 그래프이다.
도달 시간 검출방법 | fast GAPD | slow GAPD |
본 실시예 | 0.7 ns | 1.5 ns |
DCFD | 1.5 ns | 2.2 ns |
MRI | 1.8 ns | 2.7 ns |
표 1 및 도 11을 함께 참조하면, fast GAPD에 대해서는 본 실시예의 시간 분해능는 DCFD 방법에 비해서는 53% 향상되고 MRI 방법에 비해서는 61% 향상되었다. slow GAPD에 대해서는 본 실시예의 시간 분해능는 DCFD 방법에 비해서는 32% 향상되고 MRI 방법에 비해서는 44% 향상되었다.
도 12는 본 실시예에 따른 PET 이미지를 보여주는 도면이고, 도 13은 도 12의 이미지의 프로파일을 보여주는 도면이다.
도 12 및 도 13를 참조하면, 획득된 이미지 공간 해상도(spatial resolution)는 3.1 mm FWHM임을 확인할 수 있으며, 향상된 성능의 PET 이미지를 나타냄을 확인할 수 있다.
본 실시예에서 사용되는 '~부'라는 용어는 소프트웨어 또는 FPGA(field-programmable gate array) 또는 ASIC과 같은 하드웨어 구성요소를 의미하며, '~부'는 어떤 역할들을 수행한다. 그렇지만 '~부'는 소프트웨어 또는 하드웨어에 한정되는 의미는 아니다. '~부'는 어드레싱할 수 있는 저장 매체에 있도록 구성될 수도 있고 하나 또는 그 이상의 프로세서들을 재생시키도록 구성될 수도 있다. 따라서, 일 예로서 '~부'는 소프트웨어 구성요소들, 객체지향 소프트웨어 구성요소들, 클래스 구성요소들 및 태스크 구성요소들과 같은 구성요소들과, 프로세스들, 함수들, 속성들, 프로시저들, 서브루틴들, 프로그램 코드의 세그먼트들, 드라이버들, 펌웨어, 마이크로코드, 회로, 데이터,데이터베이스, 데이터 구조들, 테이블들, 어레이들, 및 변수들을 포함한다. 구성요소들과 '~부'들 안에서 제공되는 기능은 더 작은 수의 구성요소들 및 '~부'들로 결합되거나 추가적인 구성요소들과 '~부'들로 더 분리될 수 있다. 뿐만 아니라, 구성요소들 및 '~부'들은 디바이스 또는 보안 멀티미디어카드 내의 하나 또는 그 이상의 CPU들을 재생시키도록 구현될 수도 있다.
이상에서 본 발명의 실시예에 대하여 상세하게 설명하였지만 본 발명의 권리범위는 이에 한정되는 것은 아니고 다음의 청구범위에서 정의하고 있는 본 발명의 기본 개념을 이용한 당업자의 여러 변형 및 개량 형태 또한 본 발명의 권리범위에 속하는 것이다.
또한, 프로그램은 기록 매체로서의 하드 디스크나 ROM(Read Only Memory)에 미리 기록하여 둘 수가 있거나 플렉시블 디스크, CD-ROM(Compact Disc Read Only Memory), MO(Magneto optical)디스크, DVD(Digital Versatile Disc), 자기 디스크, 반도체 메모리 등의 착탈 가능 기록 매체에, 일시 목표 또는 영속적으로 저장해 둘 수 있다. 이와 같은 착탈 가능 기록 매체는 패키지 소프트웨어로서 제공할 수 있다.
이상 본 발명에 대하여 실시예를 참조하여 설명하였지만, 해당 기술 분야의 통상의 지식을 가진 자는 본 발명의 기술적 사상 및 영역으로부터 벗어나지 않는 범위 내에서 본 발명을 다양하게 수정 및 변경시켜 실시할 수 있음을 이해할 수 있을 것이다. 따라서 상술한 실시예에 한정되지 않고, 본 발명은 이하의 특허청구범위의 범위 내의 모든 실시예들을 포함한다고 할 것이다.
800: PET 신호처리모듈 810: 아날로그-디지털 변환부(ADC)
820: 신호처리부 830: 제어부
840: 저장부
820: 신호처리부 830: 제어부
840: 저장부
Claims (7)
- PET 검출기의 섬광결정에서 검출한 섬광을 전기적인 신호로 변환한 광센서로부터 전기적인 신호를 입력받아 아날로그-디지털 변환기(ADC)를 통해 변환된 감마선 펄스 신호에 대한 디지털 신호로부터 감마선 펄스 신호의 도달 시간을 검출하는 방법에 있어서,
(a) 소정의 시간 간격을 가진 상기 감마선 펄스 신호에 대한 디지털 신호의 전압값의 초기 상승선을 검출하는 단계;
(b) 상기 초기 상승선 이전에 획득된 상기 소정의 시간 간격을 가진 상기 디지털 신호의 통계적 처리를 통하여 기준선을 검출하는 단계;
(c) 상기 초기 상승선 및 상기 기준선과의 교차점을 계산하는 단계; 및
(d) 상기 교차점을 이용하여 상기 감마선 펄스 신호의 도달 시간을 검출하는 단계를 포함하는 감마선 펄스 신호의 도달 시간 검출방법. - 제 1 항에 있어서,
상기 (a) 단계는
(a-1) 최초의 상기 감마선 펄스 신호로부터 소정 시간 지연된 펄스 신호를 빼서 양의 피크 및 음의 피크를 가진 양극 펄스를 획득하는 단계;
(a-2) 상기 양의 피크를 검출하는 단계; 및
(a-3) 상기 초기 상승선을 위한 트리거링을 생성하는 단계를 포함하는 것을 특징으로 하는 감마선 펄스 신호의 도달 시간 검출방법. - 제 1 항에 있어서,
상기 소정의 시간 간격은 상기 아날로그-디지털 변환기(ADC)의 샘플링 주파수에 따라 변화되는 것을 특징으로 하는 감마선 펄스 신호의 도달 시간 검출방법. - 제 1 항에 있어서,
상기 (a) 단계 이전에
상기 감마선 펄스 신호에 대한 디지털 신호를 증폭시키는 단계를 더 포함하는 것을 특징으로 하는 감마선 펄스 신호의 도달 시간 검출방법. - 제1항 내지 제4항 중 어느 한 항에 기재된 감마선 펄스 신호의 도달 시간 검출방법을 수행하기 위해 디지털 처리 장치에 의해 실행될 수 있는 명령어들의 프로그램이 구현되어 있으며, 상기 디지털 처리 장치에 의해 판독될 수 있는 프로그램이 기록된 기록매체.
- PET 검출기의 섬광결정에서 검출한 섬광을 전기적인 신호로 변환한 광센서로부터 전기적인 신호를 입력받아 샘플링을 통하여 감마선 펄스 신호에 대한 디지털 신호를 출력하는 아날로그-디지털 변환부;및
상기 디지털 신호를 입력 받아 소정의 시간 간격을 가진 상기 디지털 신호의 전압값의 초기 상승선을 검출하고 상기 초기 상승선 이전에 획득된 상기 소정의 시간 간격을 가진 상기 디지털 신호의 통계적 처리를 통하여 기준선을 검출하여 상기 초기 상승선 및 상기 기준선과의 교차점으로부터 상기 감마선 펄스 신호의 도달 시간을 검출하는 신호처리부를 포함하는 PET 신호처리모듈. - 제 6 항에 있어서,
상기 소정의 시간 간격은 상기 아날로그-디지털 변환부(ADC)의 주파수에 따라 변화되는 것을 특징으로 하는 PET 신호처리모듈.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020100039492A KR101092648B1 (ko) | 2010-04-28 | 2010-04-28 | Pet 검출기에서의 감마선 펄스 신호의 도달 시간 검출방법 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020100039492A KR101092648B1 (ko) | 2010-04-28 | 2010-04-28 | Pet 검출기에서의 감마선 펄스 신호의 도달 시간 검출방법 |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20110120015A KR20110120015A (ko) | 2011-11-03 |
KR101092648B1 true KR101092648B1 (ko) | 2011-12-13 |
Family
ID=45391302
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020100039492A KR101092648B1 (ko) | 2010-04-28 | 2010-04-28 | Pet 검출기에서의 감마선 펄스 신호의 도달 시간 검출방법 |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR101092648B1 (ko) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006029475A1 (en) | 2004-09-16 | 2006-03-23 | Southern Innovation International Pty Ltd | Method and apparatus for resolving individual signals in detector output data. |
CN104020500A (zh) | 2008-03-31 | 2014-09-03 | 南方创新国际股份有限公司 | 用于钻孔测井的方法和装置 |
JP2011519415A (ja) | 2008-03-31 | 2011-07-07 | サザン イノヴェーション インターナショナル プロプライアトリー リミテッド | 個別信号解像度を用いた放射線画像形成法 |
US8954300B2 (en) | 2008-03-31 | 2015-02-10 | Southern Innovation International Pty Ltd. | Screening method and apparatus |
JP5581333B2 (ja) | 2008-12-18 | 2014-08-27 | サザン イノヴェーション インターナショナル プロプライアトリー リミテッド | 数学的変換を使用することによってパイルアップされたパルスを分離する方法および装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005308738A (ja) | 2004-04-02 | 2005-11-04 | National Health Research Inst | デジタル化ポジトロン断層法(pet)放射イベント装置および使用法 |
JP2009097670A (ja) | 2007-10-18 | 2009-05-07 | Honda Motor Co Ltd | シフトインジケータ装置 |
-
2010
- 2010-04-28 KR KR1020100039492A patent/KR101092648B1/ko active IP Right Grant
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005308738A (ja) | 2004-04-02 | 2005-11-04 | National Health Research Inst | デジタル化ポジトロン断層法(pet)放射イベント装置および使用法 |
JP2009097670A (ja) | 2007-10-18 | 2009-05-07 | Honda Motor Co Ltd | シフトインジケータ装置 |
Also Published As
Publication number | Publication date |
---|---|
KR20110120015A (ko) | 2011-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Surti et al. | Advances in time-of-flight PET | |
JP5667446B2 (ja) | ポジトロン放出断層撮影法のパルス開始時間推定方法および一致ペア識別方法、ポジトロン放出断層撮影法スキャナ | |
KR101169708B1 (ko) | 큰 면적을 가진 마이크로셀로 구성된 gapd를 이용한 pet 검출기 모듈 | |
US6297506B1 (en) | System and method for reducing pile-up errors in multi-crystal gamma ray detector applications | |
Berg et al. | A combined time‐of‐flight and depth‐of‐interaction detector for total‐body positron emission tomography | |
US9606199B2 (en) | MR-compatible blood sampling system for PET imaging applications in combined PET/MR imaging system | |
US9903961B1 (en) | Photodetector array readout multiplexer having summing, pulse shaping, and dynamic-switching circuits | |
KR101092648B1 (ko) | Pet 검출기에서의 감마선 펄스 신호의 도달 시간 검출방법 | |
US8735831B2 (en) | Trigger methods in nuclear medical imaging | |
US20100065746A1 (en) | Systems and methods for calibrating a silicon photomultiplier-based positron emission tomography system | |
US20100084560A1 (en) | Timing response improvement in light-sharing detectors | |
KR102063828B1 (ko) | 방사선 영상 기기의 신호 검출 방법 및 그 방사선 영상 기기 | |
JP2006266996A (ja) | γ線の同時計数方法、画像作成方法及び核医学診断装置 | |
Nakanishi et al. | Comparison of noise equivalent count rates (NECRs) for the PET systems with different ring diameter and electronics | |
Wei et al. | A side-by-side LYSO/GAGG phoswich detector aiming for SPECT imaging | |
KR101085354B1 (ko) | 모델링을 이용한 pet―mri 융합장치에서의 pet 신호의 잡음제거방법 및 이를 이용한 pet―mri 융합장치에서의 pet 시스템 | |
KR101174485B1 (ko) | Pet 장치 및 pet 장치의 신호 처리 방법 | |
KR101595929B1 (ko) | 양전자방출단층촬영장치용 반응 위치 판별 방법 및 시스템 | |
KR101450806B1 (ko) | 양전자방출 단층촬영기기의 센서 이득 제어 방법 및 시스템 | |
Moses | Positron emission mammography imaging | |
US8754376B2 (en) | Systems and methods for determining a zero baseline value of a channel from a detector device | |
KR20110121081A (ko) | 필터링을 이용한 pet―mri 융합장치에서의 pet 신호의 잡음제거방법 및 이를 이용한 pet―mri 융합장치에서의 pet 시스템 | |
KR101646620B1 (ko) | 필터링을 이용한 pet-mri융합장치에서의 pet 신호의 잡음제거방법 및 이를 이용한 pet-mri 융합장치에서의 pet 시스템 | |
KR20060109380A (ko) | 3개 이상의 결정층이 형성된 섬광검출기 및 이를 이용한양전자 방출 단층촬영장치 | |
Cui et al. | Compact CdZnTe-based gamma camera for prostate cancer imaging |
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 | ||
FPAY | Annual fee payment |
Payment date: 20141205 Year of fee payment: 4 |
|
FPAY | Annual fee payment |
Payment date: 20151120 Year of fee payment: 5 |
|
FPAY | Annual fee payment |
Payment date: 20161123 Year of fee payment: 6 |
|
FPAY | Annual fee payment |
Payment date: 20171124 Year of fee payment: 7 |
|
FPAY | Annual fee payment |
Payment date: 20181226 Year of fee payment: 8 |