1
|
Yu Z, Rahman MA, Schindler T, Laforest R, Jha AK. A physics and learning-based transmission-less attenuation compensation method for SPECT. PROCEEDINGS OF SPIE--THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING 2021; 11595. [PMID: 34658480 DOI: 10.1117/12.2582350] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Attenuation compensation (AC) is a pre-requisite for reliable quantification and beneficial for visual interpretation tasks in single-photon emission computed tomography (SPECT). Typical AC methods require the availability of an attenuation map, which is obtained using a transmission scan, such as a CT scan. This has several disadvantages such as increased radiation dose, higher costs, and possible misalignment between SPECT and CT scans. Also, often a CT scan is unavailable. In this context, we and others are showing that scattered photons in SPECT contain information to estimate the attenuation distribution. To exploit this observation, we propose a physics and learning-based method that uses the SPECT emission data in the photopeak and scatter windows to perform transmission-less AC in SPECT. The proposed method uses data acquired in the scatter window to reconstruct an initial estimate of the attenuation map using a physics-based approach. A convolutional neural network is then trained to segment this initial estimate into different regions. Pre-defined attenuation coefficients are assigned to these regions, yielding the reconstructed attenuation map, which is then used to reconstruct the activity distribution using an ordered subsets expectation maximization (OSEM)-based reconstruction approach. We objectively evaluated the performance of this method using highly realistic simulation studies conducted on the clinically relevant task of detecting perfusion defects in myocardial perfusion SPECT. Our results showed no statistically significant differences between the performance achieved using the proposed method and that with the true attenuation maps. Visually, the images reconstructed using the proposed method looked similar to those with the true attenuation map. Overall, these results provide evidence of the capability of the proposed method to perform transmission-less AC and motivate further evaluation.
Collapse
Affiliation(s)
- Zitong Yu
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO, USA, 63130
| | - Md Ashequr Rahman
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO, USA, 63130
| | - Thomas Schindler
- Mallinckrodt Institute of Radiology, Washington University in St. Louis, St. Louis, MO, USA, 63110
| | - Richard Laforest
- Mallinckrodt Institute of Radiology, Washington University in St. Louis, St. Louis, MO, USA, 63110
| | - Abhinav K Jha
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO, USA, 63130.,Mallinckrodt Institute of Radiology, Washington University in St. Louis, St. Louis, MO, USA, 63110
| |
Collapse
|
2
|
Rahman A, Zhu Y, Clarkson E, Kupinski MA, Frey EC, Jha AK. Fisher information analysis of list-mode SPECT emission data for joint estimation of activity and attenuation distribution. INVERSE PROBLEMS 2020; 36:084002. [PMID: 33071423 PMCID: PMC7561050 DOI: 10.1088/1361-6420/ab958b] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
The potential to perform attenuation and scatter compensation (ASC) in single-photon emission computed tomography (SPECT) imaging without a separate transmission scan is highly significant. In this context, attenuation in SPECT is primarily due to Compton scattering, where the probability of Compton scatter is proportional to the attenuation coefficient of the tissue and the energy of the scattered photon and the scattering angle are related. Based on this premise, we investigated whether the SPECT scattered-photon data acquired in list-mode (LM) format and including the energy information can be used to estimate the attenuation map. For this purpose, we propose a Fisher-information-based method that yields the Cramer-Rao bound (CRB) for the task of jointly estimating the activity and attenuation distribution using only the SPECT emission data. In the process, a path-based formalism to process the LM SPECT emission data, including the scattered-photon data, is proposed. The Fisher information method was implemented on NVIDIA graphics processing units (GPU) for acceleration. The method was applied to analyze the information content of SPECT LM emission data, which contains up to first-order scattered events, in a simulated SPECT system with parameters modeling a clinical system using realistic computational studies with 2-D digital synthetic and anthropomorphic phantoms. The method was also applied to LM data containing up to second-order scatter for a synthetic phantom. Experiments with anthropomorphic phantoms simulated myocardial perfusion and dopamine transporter (DaT)-Scan SPECT studies. The results show that the CRB obtained for the attenuation and activity coefficients was typically much lower than the true value of these coefficients. An increase in the number of detected photons yielded lower CRB for both the attenuation and activity coefficients. Further, we observed that systems with better energy resolution yielded a lower CRB for the attenuation coefficient. Overall, the results provide evidence that LM SPECT emission data, including the scattered photons, contains information to jointly estimate the activity and attenuation coefficients.
Collapse
Affiliation(s)
- Ashequr Rahman
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO, USA
- Mallinckrodt Institute of Radiology, Washington University in St. Louis, St. Louis, MO, USA
| | - Yansong Zhu
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, MD, USA
- Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, MD, USA
- Department of Physics & Astronomy, University of British Columbia, Canada
| | - Eric Clarkson
- College of Optical Sciences, University of Arizona, Tucson AZ, USA
| | | | - Eric C Frey
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, MD, USA
| | - Abhinav K Jha
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO, USA
- Mallinckrodt Institute of Radiology, Washington University in St. Louis, St. Louis, MO, USA
| |
Collapse
|
3
|
Protonotarios NE, Spyrou GM, Kastis GA. Automatic cumulative sums contour detection of FBP-reconstructed multi-object nuclear medicine images. Comput Biol Med 2017; 85:43-52. [PMID: 28433871 DOI: 10.1016/j.compbiomed.2017.04.010] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2017] [Revised: 03/27/2017] [Accepted: 04/12/2017] [Indexed: 11/26/2022]
Abstract
The problem of determining the contours of objects in nuclear medicine images has been studied extensively in the past, however most of the analysis has focused on a single object as opposed to multiple objects. The aim of this work is to develop an automated method for determining the contour of multiple objects in positron emission tomography (PET) and single photon emission computed tomography (SPECT) filtered backprojection (FBP) reconstructed images. These contours can be used for computing body edges for attenuation correction in PET and SPECT, as well as for eliminating streak artifacts outside the objects, which could be useful in compressive sensing reconstruction. Contour detection has been accomplished by applying a modified cumulative sums (CUSUM) scheme in the sinogram. Our approach automatically detects all objects in the image, without requiring a priori knowledge of the number of distinct objects in the reconstructed image. This method has been tested in simulated phantoms, such as an image-quality (IQ) phantom and two digital multi-object phantoms, as well as a real NEMA phantom and a clinical thoracic study. For this purpose, a GE Discovery PET scanner was employed. The detected contours achieved root mean square accuracy of 1.14 pixels, 1.69 pixels and 3.28 pixels and a Hausdorff distance of 3.13, 3.12 and 4.50 pixels, for the simulated image-quality phantom PET study, the real NEMA phantom and the clinical thoracic study, respectively. These results correspond to a significant improvement over recent results obtained in similar studies. Furthermore, we obtained an optimal sub-pattern assignment (OSPA) localization error of 0.94 and 1.48, for the two-objects and three-objects simulated phantoms, respectively. Our method performs efficiently for sets of convex objects and hence it provides a robust tool for automatic contour determination with precise results.
Collapse
Affiliation(s)
- Nicholas E Protonotarios
- Research Center of Mathematics, Academy of Athens, Athens 11527, Greece; Department of Mathematics, National Technical University of Athens, Zografou Campus, Athens 15780, Greece.
| | - George M Spyrou
- Center of Systems Biology, Biomedical Research Foundation of the Academy of Athens, Athens 11527, Greece; Bioinformatics ERA Chair, The Cyprus Institute of Neurology and Genetics, Ayios Dometios, 2370 Nicosia, Cyprus
| | - George A Kastis
- Research Center of Mathematics, Academy of Athens, Athens 11527, Greece
| |
Collapse
|
4
|
Berker Y, Li Y. Attenuation correction in emission tomography using the emission data--A review. Med Phys 2016; 43:807-32. [PMID: 26843243 PMCID: PMC4715007 DOI: 10.1118/1.4938264] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2015] [Revised: 11/19/2015] [Accepted: 11/25/2015] [Indexed: 11/07/2022] Open
Abstract
The problem of attenuation correction (AC) for quantitative positron emission tomography (PET) had been considered solved to a large extent after the commercial availability of devices combining PET with computed tomography (CT) in 2001; single photon emission computed tomography (SPECT) has seen a similar development. However, stimulated in particular by technical advances toward clinical systems combining PET and magnetic resonance imaging (MRI), research interest in alternative approaches for PET AC has grown substantially in the last years. In this comprehensive literature review, the authors first present theoretical results with relevance to simultaneous reconstruction of attenuation and activity. The authors then look back at the early history of this research area especially in PET; since this history is closely interwoven with that of similar approaches in SPECT, these will also be covered. We then review algorithmic advances in PET, including analytic and iterative algorithms. The analytic approaches are either based on the Helgason-Ludwig data consistency conditions of the Radon transform, or generalizations of John's partial differential equation; with respect to iterative methods, we discuss maximum likelihood reconstruction of attenuation and activity (MLAA), the maximum likelihood attenuation correction factors (MLACF) algorithm, and their offspring. The description of methods is followed by a structured account of applications for simultaneous reconstruction techniques: this discussion covers organ-specific applications, applications specific to PET/MRI, applications using supplemental transmission information, and motion-aware applications. After briefly summarizing SPECT applications, we consider recent developments using emission data other than unscattered photons. In summary, developments using time-of-flight (TOF) PET emission data for AC have shown promising advances and open a wide range of applications. These techniques may both remedy deficiencies of purely MRI-based AC approaches in PET/MRI and improve standalone PET imaging.
Collapse
Affiliation(s)
- Yannick Berker
- Department of Radiology, University of Pennsylvania, 3620 Hamilton Walk, Philadelphia, Pennsylvania 19104
| | - Yusheng Li
- Department of Radiology, University of Pennsylvania, 3620 Hamilton Walk, Philadelphia, Pennsylvania 19104
| |
Collapse
|
5
|
Angelis GI, Kyme AZ, Ryder WJ, Fulton RR, Meikle SR. Attenuation correction for freely moving small animal brain PET studies based on a virtual scanner geometry. Phys Med Biol 2014; 59:5651-66. [DOI: 10.1088/0031-9155/59/19/5651] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
6
|
Chang T, Diab RH, Clark JW, Mawlawi OR. Investigating the use of nonattenuation corrected PET images for the attenuation correction of PET data. Med Phys 2014; 40:082508. [PMID: 23927353 DOI: 10.1118/1.4816304] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
PURPOSE The aim of this study is to investigate the feasibility of using the nonattenuated PET images (PET-NAC) as a means for the AC of PET data. METHODS A three-step iterative segmentation process is proposed. In step 1, a patient's body contour is segmented from the PET-NAC using an active contour algorithm. Voxels inside the contour are then assigned a value of 0.096 cm(-1) to represent the attenuation coefficient of soft tissue at 511 keV. This segmented attenuation map is then used to correct for attenuation the raw PET data and the resulting PET images are used as the input to Step 2 of the process. In step 2, the lung region is segmented using an optimal thresholding approach and the corresponding voxels are assigned a value of 0.024 cm(-1) representing the attenuation coefficients of lung tissue at 511 keV. The updated attenuation map is then used for a second time to correct for attenuation the raw PET data, and the resulting PET images are used as the input to step 3. The purpose of Step 3 is to delineate parts of the heart and liver in the lung contour using a region growing approach since these parts were unavoidably excluded in the lung contour in step 2. These parts are then corrected by using a value of 0.096 cm(-1) in the attenuation map. Finally the attenuation coefficients of the bed are included based on CT images to eliminate the impact of the couch on the accuracy of AC. The final attenuation map is then used to AC the raw PET data and generates the final PET image, which we name iterative AC PET (PET-IAC). To assess the proposed segmentation approach, a phantom and 14 patients (with a total of 55 lesions including bone) were scanned on a GE Discovery-RX PET∕CT scanner. PET-IAC images were generated using the proposed process and compared to those of CT-AC PET (PET-CTAC). Visual inspection, lesion SUV, and voxel by voxel histograms between PET-IAC and PET-CTAC for phantom and patient studies were performed to assess the accuracy of image quantification. RESULTS Visual inspection showed a small difference in lung parenchyma between the PET-IAC and PET-CTAC. Tumor SUV based on PET-IAC were on average different by 3%±9% (6%±7%) compared to the SUVs from the PET-CTAC in the phantom (patient) studies. For bone lesions only, the average difference was 3%±6%. The histogram comparing PET-CTAC and PET-IAC resulted in an average regression line of y=(1.08±0.07)x+(0.00007±0.0013), with R2=0.978±0.0057. CONCLUSIONS Preliminary results suggest that PET-NAC for the AC of PET images is feasible. Such an approach can potentially be used for dedicated PET or PET∕MR hybrid systems while minimizing scan time or potential image artifacts, respectively.
Collapse
Affiliation(s)
- Tingting Chang
- Department of Electrical and Computer Engineering, Rice University, 6100 Main Street, MS-366, Houston, Texas 77005, USA
| | | | | | | |
Collapse
|
7
|
Jha AK, Clarkson E, Kupinski MA, Barrett HH. Joint reconstruction of activity and attenuation map using LM SPECT emission data. PROCEEDINGS OF SPIE--THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING 2013; 8668. [PMID: 26236067 DOI: 10.1117/12.2008111] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Attenuation and scatter correction in single photon emission computed tomography (SPECT) imaging often requires a computed tomography (CT) scan to compute the attenuation map of the patient. This results in increased radiation dose for the patient, and also has other disadvantages such as increased costs and hardware complexity. Attenuation in SPECT is a direct consequence of Compton scattering, and therefore, if the scattered photon data can give information about the attenuation map, then the CT scan may not be required. In this paper, we investigate the possibility of joint reconstruction of the activity and attenuation map using list-mode (LM) SPECT emission data, including the scattered-photon data. We propose a path-based formalism to process scattered-photon data. Following this, we derive analytic expressions to compute the Cramér-Rao bound (CRB) of the activity and attenuation map estimates, using which, we can explore the fundamental limit of information-retrieval capacity from LM SPECT emission data. We then suggest a maximum-likelihood (ML) scheme that uses the LM emission data to jointly reconstruct the activity and attenuation map. We also propose an expectation-maximization (EM) algorithm to compute the ML solution.
Collapse
Affiliation(s)
- Abhinav K Jha
- College of Optical Sciences, University of Arizona, Tucson, AZ, USA
| | - Eric Clarkson
- College of Optical Sciences, University of Arizona, Tucson, AZ, USA ; Department of Medical Imaging, University of Arizona, Tucson, AZ, USA
| | - Matthew A Kupinski
- College of Optical Sciences, University of Arizona, Tucson, AZ, USA ; Department of Medical Imaging, University of Arizona, Tucson, AZ, USA
| | - Harrison H Barrett
- College of Optical Sciences, University of Arizona, Tucson, AZ, USA ; Department of Medical Imaging, University of Arizona, Tucson, AZ, USA
| |
Collapse
|
8
|
Nuyts J, Bal G, Kehren F, Fenchel M, Michel C, Watson C. Completion of a truncated attenuation image from the attenuated PET emission data. IEEE TRANSACTIONS ON MEDICAL IMAGING 2013; 32:237-246. [PMID: 23014717 DOI: 10.1109/tmi.2012.2220376] [Citation(s) in RCA: 91] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
Positron emission tomographs (PETs) are currently almost exclusively designed as hybrid systems. The current standard is the PET/CT combination, while prototype PET/MRI systems are being studied by several research groups. One problem in these systems is that the transaxial field-of-view of the second system is smaller than that of the PET camera and does not provide complete attenuation data. Because this second system provides the image for PET attenuation and scatter correction, the smaller FOV causes truncation of the attenuation map, producing bias in the attenuation corrected activity image. In this paper, we propose a maximum-a-posteriori algorithm for estimating the missing part of the attenuation map from the PET emission data. The method is evaluated on five artificially truncated 18F-FDGPET/CT studies, where it reduced the error on the reconstructed PET activities from 20% to less than 7%. The results on a PET/MRI patient study with 18F-FDG are presented as well.
Collapse
Affiliation(s)
- Johan Nuyts
- Nuclear Medicine, KU Leuven, B-3000 Leuven, Belgium.
| | | | | | | | | | | |
Collapse
|
9
|
Martinez-Möller A, Nekolla SG. Attenuation correction for PET/MR: Problems, novel approaches and practical solutions. Z Med Phys 2012; 22:299-310. [DOI: 10.1016/j.zemedi.2012.08.003] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2012] [Revised: 07/12/2012] [Accepted: 08/08/2012] [Indexed: 10/27/2022]
|
10
|
Catana C, van der Kouwe A, Benner T, Michel CJ, Hamm M, Fenchel M, Fischl B, Rosen B, Schmand M, Sorensen AG. Toward implementing an MRI-based PET attenuation-correction method for neurologic studies on the MR-PET brain prototype. J Nucl Med 2010; 51:1431-8. [PMID: 20810759 DOI: 10.2967/jnumed.109.069112] [Citation(s) in RCA: 256] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
UNLABELLED Several factors have to be considered for implementing an accurate attenuation-correction (AC) method in a combined MR-PET scanner. In this work, some of these challenges were investigated, and an AC method based entirely on the MRI data obtained with a single dedicated sequence was developed and used for neurologic studies performed with the MR-PET human brain scanner prototype. METHODS The focus was on the problem of bone-air segmentation, selection of the linear attenuation coefficient for bone, and positioning of the radiofrequency coil. The impact of these factors on PET data quantification was studied in simulations and experimental measurements performed on the combined MR-PET scanner. A novel dual-echo ultrashort echo time (DUTE) MRI sequence was proposed for head imaging. Simultaneous MR-PET data were acquired, and the PET images reconstructed using the proposed DUTE MRI-based AC method were compared with the PET images that had been reconstructed using a CT-based AC method. RESULTS Our data suggest that incorrectly accounting for the bone tissue attenuation can lead to large underestimations (>20%) of the radiotracer concentration in the cortex. Assigning a linear attenuation coefficient of 0.143 or 0.151 cm(-1) to bone tissue appears to give the best trade-off between bias and variability in the resulting images. Not identifying the internal air cavities introduces large overestimations (>20%) in adjacent structures. On the basis of these results, the segmented CT AC method was established as the silver standard for the segmented MRI-based AC method. For an integrated MR-PET scanner, in particular, ignoring the radiofrequency coil attenuation can cause large underestimations (i.e., <or=50%) in the reconstructed images. Furthermore, the coil location in the PET field of view has to be accurately known. High-quality bone-air segmentation can be performed using the DUTE data. The PET images obtained using the DUTE MRI- and CT-based AC methods compare favorably in most of the brain structures. CONCLUSION A DUTE MRI-based AC method considering all these factors was implemented. Preliminary results suggest that this method could potentially be as accurate as the segmented CT method and could be used for quantitative neurologic MR-PET studies.
Collapse
Affiliation(s)
- Ciprian Catana
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts 02129, USA.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
11
|
Lemmens C, Montandon ML, Nuyts J, Ratib O, Dupont P, Zaidi H. Impact of metal artefacts due to EEG electrodes in brain PET/CT imaging. Phys Med Biol 2008; 53:4417-29. [PMID: 18670051 DOI: 10.1088/0031-9155/53/16/013] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The goal of this study is to investigate the impact of electroencephalogram (EEG) electrodes on the visual quality and quantification of (18)F-FDG PET images in neurological PET/CT examinations. For this purpose, the scans of 20 epilepsy patients with EEG monitoring were used. The CT data were reconstructed with filtered backprojection (FBP) and with a metal artefact reduction (MAR) algorithm. Both data sets were used for CT-based attenuation correction (AC) of the PET data. Also, a calculated AC (CALC) technique was considered. A volume of interest (VOI)-based analysis and a voxel-based quantitative analysis were performed to compare the different AC methods. Images were also evaluated visually by two observers. It was shown with simulations and phantom measurements that from the considered AC methods, the MAR-AC can be used as the reference in this setting. The visual assessment of PET images showed local hot spots outside the brain corresponding to the locations of the electrodes when using FBP-AC. In the brain, no abnormalities were observed. The quantitative analysis showed a very good correlation between PET-FBP-AC and PET-MAR-AC, with a statistically significant positive bias in the PET-FBP-AC images of about 5-7% in most brain voxels. There was also good correlation between PET-CALC-AC and PET-MAR-AC, but in the PET-CALC-AC images, regions with both a significant positive and negative bias were observed. EEG electrodes give rise to local hot spots outside the brain and a positive quantification bias in the brain. However, when diagnosis is made by mere visual assessment, the presence of EEG electrodes does not seem to alter the diagnosis. When quantification is performed, the bias becomes an issue especially when comparing brain images with and without EEG monitoring.
Collapse
Affiliation(s)
- Catherine Lemmens
- Department of Nuclear Medicine and Medical Imaging Center, University Hospital Gasthuisberg and Katholieke Universiteit Leuven, Leuven, Belgium
| | | | | | | | | | | |
Collapse
|
12
|
Generation of planar images from lung ventilation/perfusion SPECT. Ann Nucl Med 2008; 22:437-45. [DOI: 10.1007/s12149-007-0120-1] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2007] [Accepted: 12/19/2007] [Indexed: 10/21/2022]
|
13
|
Zaidi H, Montandon ML, Meikle S. Strategies for attenuation compensation in neurological PET studies. Neuroimage 2007; 34:518-41. [PMID: 17113312 DOI: 10.1016/j.neuroimage.2006.10.002] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2006] [Revised: 09/29/2006] [Accepted: 10/03/2006] [Indexed: 11/29/2022] Open
Abstract
Molecular brain imaging using positron emission tomography (PET) has evolved into a vigorous academic field and is progressively gaining importance in the clinical arena. Significant progress has been made in the design of high-resolution three-dimensional (3-D) PET units dedicated to brain research and the development of quantitative imaging protocols incorporating accurate image correction techniques and sophisticated image reconstruction algorithms. However, emerging clinical and research applications of molecular brain imaging demand even greater levels of accuracy and precision and therefore impose more constraints with respect to the quantitative capability of PET. It has long been recognized that photon attenuation in tissues is the most important physical factor degrading PET image quality and quantitative accuracy. Quantitative PET image reconstruction requires an accurate attenuation map of the object under study for the purpose of attenuation compensation. Several methods have been devised to correct for photon attenuation in neurological PET studies. Significant attention has been devoted to optimizing computational performance and to balancing conflicting requirements. Approximate methods suitable for clinical routine applications and more complicated approaches for research applications, where there is greater emphasis on accurate quantitative measurements, have been proposed. The number of scientific contributions related to this subject has been increasing steadily, which motivated the writing of this review as a snapshot of the dynamically changing field of attenuation correction in cerebral 3D PET. This paper presents the physical and methodological basis of photon attenuation and summarizes state of the art developments in algorithms used to derive the attenuation map aiming at accurate attenuation compensation of brain PET data. Future prospects, research trends and challenges are identified and directions for future research are discussed.
Collapse
Affiliation(s)
- Habib Zaidi
- Division of Nuclear Medicine, Geneva University Hospital, CH-1211 Geneva 4, Switzerland.
| | | | | |
Collapse
|
14
|
Barnden LR, Dickson J, Hutton BF. Detection and validation of the body edge in low count emission tomography images. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2006; 84:153-61. [PMID: 16978729 DOI: 10.1016/j.cmpb.2006.08.001] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/26/2006] [Revised: 07/14/2006] [Accepted: 08/03/2006] [Indexed: 05/11/2023]
Abstract
Segmentation of the body edge in tomographic images with low count, noisy edges is needed for both PET and SPECT respiratory motion correction and brain SPECT attenuation correction. To reduce noise we re-projected tomographic images and searched for edges in the projection count profiles or their spatial derivatives. Edge location versus projection angle was fitted with cosine basis functions after rejecting outliers and including information about edges of previous sections. We processed 10 s duration FDG PET of the thorax, HMPAO brain, DAT brain and lung perfusion SPECT. Stable edges for all four types of scan were achieved but with different profiles. Edges were validated against edges of coregistered CT or MRI. The best root mean square (rms) accuracy was 8.2 mm (PET) and 5.2 mm (brain SPECT). Inter-scan variability (standard deviation) in the estimated-to-control edge distance for 17 PET scans was 0.4 mm, for 25 ordered subset expectation maximisation (OSEM) reconstructed brain SPECT 1.0 mm and for 18 filtered back-projection (FBP) reconstructed brain SPECT 1.4 mm.
Collapse
Affiliation(s)
- Leighton R Barnden
- Department of Nuclear Medicine, The Queen Elizabeth Hospital, 28 Woodville Road, Woodville, SA 5011, Australia.
| | | | | |
Collapse
|
15
|
Zaidi H, Montandon ML, Slosman DO. Magnetic resonance imaging-guided attenuation and scatter corrections in three-dimensional brain positron emission tomography. Med Phys 2003; 30:937-48. [PMID: 12773003 DOI: 10.1118/1.1569270] [Citation(s) in RCA: 179] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
Reliable attenuation correction represents an essential component of the long chain of modules required for the reconstruction of artifact-free, quantitative brain positron emission tomography (PET) images. In this work we demonstrate the proof of principle of segmented magnetic resonance imaging (MRI)-guided attenuation and scatter corrections in three-dimensional (3D) brain PET. We have developed a method for attenuation correction based on registered T1-weighted MRI, eliminating the need of an additional transmission (TX) scan. The MR images were realigned to preliminary reconstructions of PET data using an automatic algorithm and then segmented by means of a fuzzy clustering technique which identifies tissues of significantly different density and composition. The voxels belonging to different regions were classified into air, skull, brain tissue and nasal sinuses. These voxels were then assigned theoretical tissue-dependent attenuation coefficients as reported in the ICRU 44 report followed by Gaussian smoothing and addition of a good statistics bed image. The MRI-derived attenuation map was then forward projected to generate attenuation correction factors (ACFs) to be used for correcting the emission (EM) data. The method was evaluated and validated on 10 patient data where TX and MRI brain images were available. Qualitative and quantitative assessment of differences between TX-guided and segmented MRI-guided 3D reconstructions were performed by visual assessment and by estimating parameters of clinical interest. The results indicated a small but noticeable improvement in image quality as a consequence of the reduction of noise propagation from TX into EM data. Considering the difficulties associated with preinjection TX-based attenuation correction and the limitations of current calculated attenuation correction, MRI-based attenuation correction in 3D brain PET would likely be the method of choice for the foreseeable future as a second best approach in a busy nuclear medicine center and could be applied to other functional brain imaging modalities such as SPECT.
Collapse
Affiliation(s)
- Habib Zaidi
- Division of Nuclear Medicine, Geneva University Hospital, CH-1211 Geneva 4, Switzerland.
| | | | | |
Collapse
|
16
|
De Volder AG. Functional brain imaging of childhood clinical disorders with PET and SPECT. Dev Sci 2002. [DOI: 10.1111/1467-7687.00374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
17
|
Vandermeeren Y, Olivier E, Sébire G, Cosnard G, Bol A, Sibomana M, Michel C, De Volder AG. Increased FDG uptake in the ipsilesional sensorimotor cortex in congenital hemiplegia. Neuroimage 2002; 15:949-60. [PMID: 11906235 DOI: 10.1006/nimg.2001.1041] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
The resting brain metabolism was estimated in six children suffering from a right congenital hemiplegia (CH) of subcortical origin. This estimate was based on the 18F-labeled 2-deoxy-2-fluoro-d-glucose (FDG) uptake measured by means of positron emission tomography and compared, using statistical parametric mapping (SPM99), with that of six control subjects. The contrast [CH children - Controls] showed that CH children had two loci of relatively higher FDG uptake. The larger voxel cluster was found in the ipsilesional hemisphere and comprised the primary motor and somatosensory cortices and left inferior parietal lobule. The other cluster was located in the contralesional hemisphere and encompassed the primary motor cortex, callosomarginal sulcus, and cingulate gyrus. The reverse contrast [Controls - CH children] showed that control subjects had a relatively higher FDG uptake bilaterally in the temporal and hippocampal gyri, the rostral part of the brain stem, the thalami, the putamen, and the superior frontal gyri. A crossed cerebellar diaschisis was not observed in CH children. This relatively higher FDG uptake in the ipsi- and contralesional motor areas of CH children stands out in contrast to the hypometabolism (diaschisis) frequently observed in adult stroke patients with a subcortical lesion. This increased FDG uptake in the disconnected ipsilesional motor areas may reflect a long-term adaptation leading, for example, to an increased synaptic density and/or activity or to a change in the density of glucose transporters.
Collapse
Affiliation(s)
- Yves Vandermeeren
- Laboratory of Neurophysiology, Université Catholique de Louvain School of Medicine, B-1200 Brussels, Belgium
| | | | | | | | | | | | | | | |
Collapse
|
18
|
Nuyts J, Dupont P, Stroobants S, Benninck R, Mortelmans L, Suetens P. Simultaneous maximum a posteriori reconstruction of attenuation and activity distributions from emission sinograms. IEEE TRANSACTIONS ON MEDICAL IMAGING 1999; 18:393-403. [PMID: 10416801 DOI: 10.1109/42.774167] [Citation(s) in RCA: 180] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
In order to perform attenuation correction in emission tomography an attenuation map is required. We propose a new method to compute this map directly from the emission sinogram, eliminating the transmission scan from the acquisition protocol. The problem is formulated as an optimization task where the objective function is a combination of the likelihood and an a priori probability. The latter uses a Gibbs prior distribution to encourage local smoothness and a multimodal distribution for the attenuation coefficients. Since the attenuation process is different in positron emission tomography (PET) and single photon emission tomography (SPECT), a separate algorithm for each case is derived. The method has been tested on mathematical phantoms and on a few clinical studies. For PET, good agreement was found between the images obtained with transmission measurements and those produced by the new algorithm in an abdominal study. For SPECT, promising simulation results have been obtained for nonhomogeneous attenuation due to the presence of the lungs.
Collapse
Affiliation(s)
- J Nuyts
- Department of Nuclear Medicine, KU Leuven, Belgium.
| | | | | | | | | | | |
Collapse
|
19
|
De Volder AG, Catalan-Ahumada M, Robert A, Bol A, Labar D, Coppens A, Michel C, Veraart C. Changes in occipital cortex activity in early blind humans using a sensory substitution device. Brain Res 1999; 826:128-34. [PMID: 10216204 DOI: 10.1016/s0006-8993(99)01275-5] [Citation(s) in RCA: 62] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The purpose of this study was to investigate the neural networks involved when using an ultrasonic echolocation device, which is a substitution prosthesis for blindness through audition. Using positron emission tomography with fluorodeoxyglucose, regional brain glucose metabolism was measured in the occipital cortex of early blind subjects and blindfolded controls who were trained to use this prosthesis. All subjects were studied under two different activation conditions: (i) during an auditory control task, (ii) using the ultrasonic echolocation device in a spatial distance and direction evaluation task. Results showed that the abnormally high metabolism already observed in early blind occipital cortex at rest [C. Veraart, A.G. De Volder, M.C. Wanet-Defalque, A. Bol, C. Michel, A.M. Goffinet, Glucose utilization in human visual cortex is, respectively elevated and decreased in early versus late blindness, Brain Res. 510 (1990) 115-121.] was also present during the control task and showed a trend to further increase during the use of the ultrasonic echolocation device. This specific difference in occipital cortex activity between the two tasks was not observed in control subjects. The metabolic recruitment of the occipital cortex in early blind subjects using a substitution prosthesis could reflect a concurrent stimulation of functional cross-modal sensory connections. Given the unfamiliarity of the task, it could be interpreted as a prolonged plasticity in the occipital cortex early deprived of visual afferences.
Collapse
Affiliation(s)
- A G De Volder
- Neural Rehabilitation Engineering Laboratory, University of Louvain, Avenue Hippocrate 54, UCL-54.46, B-1200, Brussels, Belgium
| | | | | | | | | | | | | | | |
Collapse
|
20
|
Bentourkia M, Michel C, Ferriere G, Bol A, Coppens A, Sibomana M, Bausart R, Labar D, De Volder AG. Evolution of brain glucose metabolism with age in epileptic infants, children and adolescents. Brain Dev 1998; 20:524-9. [PMID: 9840673 DOI: 10.1016/s0387-7604(98)00040-0] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
During the first years of life, the human brain undergoes repetitive modifications in its anatomical, functional, and synaptic construction to reach the complex functional organization of the adult central nervous system. As an attempt to gain further insight in those maturation processes, the evolution of cerebral metabolic activity was investigated as a function of age in epileptic infants, children and adolescents. The regional cerebral metabolic rates for glucose (rCMRGlc) were measured with positron emission tomography (PET) in 60 patients aged from 6 weeks to 19 years, who were affected by complex partial epilepsy. They were scanned at rest, without premedication, in similar conditions to 20 epileptic adults and in 49 adult controls. The distribution of brain metabolic activity successively extended from sensorimotor areas and thalamus in epileptic newborns to temporo-parietal and frontal cortices and reached the adult pattern after 1 year of age. The measured rCMRGlc in the cerebral cortex, excluding the epileptic lesions, increased from low values in infants to a maximum between 4 and 12 years, before it declined to stabilize at the end of the second decade of life. Similar age-related changes in glucose metabolic rates were not observed in the adult groups. Despite the use of medications, the observed variations of rCMRGlc with age in young epileptic humans confirm those previously described in pediatric subjects. These metabolic changes are in full agreement with the current knowledge of the synaptic density evolution in the human brain.
Collapse
Affiliation(s)
- M Bentourkia
- Positron Tomography Laboratory, University of Louvain, School of Medicine, Louvain-La-Neuve, Belgium
| | | | | | | | | | | | | | | | | |
Collapse
|
21
|
van Elmbt LR, Keyeux A, Demeure R. Validation of automated brain contour determination in normal and abnormal cerebral single-photon emission tomography. EUROPEAN JOURNAL OF NUCLEAR MEDICINE 1995; 22:537-42. [PMID: 7556299 DOI: 10.1007/bf00817278] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
A contour detection algorithm for cerebral studies, using the method of Tomitani, has been implemented on a single-photon emission tomographic (SPET) system. It is based on the detetion by threshold of the brain edge in the sinogram and does not depend on the reconstruction algorithm. Thirteen normal subjects underwent an examination on both computed tomography (CT) and SPET using a head holder to ensure the reproducibility of the positioning. The CT scan contour of the brain was drawn manually according to the brain parenchyma limits. The SPET brain contour was obtained by use of the Tomitani algorithm after the threshold had been determined on an active cylindrical phantom. Using a threshold of 37% of the maximum uptake, the length of the contour as well as the area obtained with SPET and CT were not found to be statistically different. The method of Tomitani, which is simpler and faster then previous methods, provides contours which superimpose very well with CT scan images. Application to patients with unilateral pathological defects is possible by requiring that the contour is symmetrical.
Collapse
Affiliation(s)
- L R van Elmbt
- Center of Nuclear Medecine, Université Catholique de Louvain, School of Medicine, Brussels, Belgium
| | | | | |
Collapse
|
22
|
Bettinardi V, Gilardi MC, Cargnel S, Rizzo G, Teräs M, Striano G, Fazio F. A hybrid method of attenuation correction for positron emission tomography brain studies. EUROPEAN JOURNAL OF NUCLEAR MEDICINE 1994; 21:1279-84. [PMID: 7875164 DOI: 10.1007/bf02426690] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
A hybrid method for attenuation correction (HAC) in positron emission tomography (PET) brain studies is proposed. The technique requires the acquisition of two short (1 min) transmission scans immediately before or after the emission study, with the patient and the head fixation system in place and after removing the patient from the scanner with the head fixation system alone. The method combines a uniform map of attenuation coefficients for the patient's head with measured attenuation coefficients for the head fixation system to generate a hybrid attenuation map. The HAC method was calibrated on 30 PET cerebral studies for comparison with the conventional measured attenuation correction method by ROI analysis. Average differences of less than 3% were found for cortical and subcortical regions. The HAC technique is particularly suitable in a PET clinical environment, allowing a reduction of the total study time, greater comfort for patients and an increase in patient throughput.
Collapse
Affiliation(s)
- V Bettinardi
- INB-CNR, Scientific Institute H San Raffaele, Univeristy of Milan, Italy
| | | | | | | | | | | | | |
Collapse
|
23
|
De Volder AG, Gadisseux JF, Michel CJ, Maloteaux JM, Bol AC, Grandin CB, Duprez TP, Evrard P. Brain glucose utilization in band heterotopia: synaptic activity of "double cortex". Pediatr Neurol 1994; 11:290-4. [PMID: 7702687 DOI: 10.1016/0887-8994(94)90003-5] [Citation(s) in RCA: 44] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Regional brain glucose utilization was investigated with positron emission tomography and fluorodeoxyglucose in 2 patients with a seizure disorder associated with diffuse band heterotopia, a condition known as "double cortex." Although 1 patient was examined shortly after the onset of the first seizures, the other had a long history of intractable epilepsy before examination. Magnetic resonance imaging revealed a symmetric and generalized band of ectopic gray matter and an overlying normal-looking cortex, without focal abnormality. Metabolic studies yielded comparable results in both patients, with similar and even higher glucose uptake in the layer of gray matter heterotopia compared to the normal cortex. These data suggest the persistence of some synaptic activity in the heterotopic neurons, which seems unaffected by age or by the time-course of epilepsy.
Collapse
Affiliation(s)
- A G De Volder
- Department of Pediatric Neurology, University of Louvain, Brussels, Belgium
| | | | | | | | | | | | | | | |
Collapse
|
24
|
De Volder AG, Michel C, Thauvoy C, Willems G, Ferrière G. Brain glucose utilisation in acquired childhood aphasia associated with a sylvian arachnoid cyst: recovery after shunting as demonstrated by PET. J Neurol Neurosurg Psychiatry 1994; 57:296-300. [PMID: 7512624 PMCID: PMC1072817 DOI: 10.1136/jnnp.57.3.296] [Citation(s) in RCA: 59] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Regional brain glucose utilisation was investigated with PET and fluorodeoxyglucose (FDG) in a case of epileptic aphasia (Landau-Kleffner syndrome) associated with a left sylvian arachnoid cyst. CT and MRI had failed to disclose any mass effect of the cyst on surrounding brain structures. Sequential metabolic measurements showed a comparable pronounced hypometabolism in cortical regions around the cyst, involving speech areas, and suggested mild but chronic compression of the developing brain. After placement of a cyst-peritoneal shunt system, significant metabolic improvement occurred in all cortical regions, especially the inferior frontal gyrus and the perisylvian area, with predominant residual deficit in the left superior temporal gyrus. These findings were correlated with a pronounced increase in word fluency and slower progress in verbal auditory comprehension. This report suggests that PET is able to evaluate the functional disturbances associated with expanding arachnoid cysts, and to follow the neurological improvement after drainage.
Collapse
Affiliation(s)
- A G De Volder
- Positron Tomography Laboratory, University of Louvain, Belgium
| | | | | | | | | |
Collapse
|
25
|
Catalán-Ahumada M, Deggouj N, De Volder A, Melin J, Michel C, Veraart C. High metabolic activity demonstrated by positron emission tomography in human auditory cortex in case of deafness of early onset. Brain Res 1993; 623:287-92. [PMID: 8221110 DOI: 10.1016/0006-8993(93)91439-y] [Citation(s) in RCA: 42] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Glucose metabolism has been studied in the auditory cortex of human subjects with deafness of early onset, and compared to normal subjects with ears plugged. The metabolism in the auditory cortex and in the association auditory cortex was higher in deaf subjects than in normal subjects. This result is compared to similar observations that we made previously in the visual cortex of human subjects with blindness of early onset.
Collapse
Affiliation(s)
- M Catalán-Ahumada
- Neural Rehabilitation Engineering Laboratory, University of Louvain, Brussels, Belgium
| | | | | | | | | | | |
Collapse
|
26
|
Wienhard K. The FDG model and its application in clinical PET studies. JOURNAL OF NEURAL TRANSMISSION. SUPPLEMENTUM 1992; 37:39-52. [PMID: 1402851 DOI: 10.1007/978-3-7091-9209-2_4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The FDG method, as it is applied in clinical PET studies is reviewed. The influence of different implementations of the method and instrumental inaccuracies on the values of cerebral metabolic rate of glucose is discussed. For the comparison of the results between different groups standardized procedures are recommended.
Collapse
Affiliation(s)
- K Wienhard
- Max-Planck-Institut für Neurologische Forschung, Köln (Lindenthal), Federal Republic of Germany
| |
Collapse
|
27
|
Bol A, Vanmelckenbeke P, Michel C, Cogneau M, Goffinet AM. Measurement of cerebral blood flow with a bolus of oxygen-15-labelled water: comparison of dynamic and integral methods. EUROPEAN JOURNAL OF NUCLEAR MEDICINE 1990; 17:234-41. [PMID: 2083558 DOI: 10.1007/bf00812363] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
A method is presented for the measurement of cerebral blood flow (CBF) with a bolus of water labelled with oxygen 15. The method, which has been evaluated in normal volunteers, is based on Kety's model, with two additional parameters to account for the difference in the time of tracer arrival in the radial and carotid arteries ("delay") and for dispersion of the tracer in the body and/or blood counting systems. It combines the advantages of: (i) dynamic data collection for estimation of delay and dispersion; (ii) robustness and linearity of CBF estimates with an integral method; and (iii) simplicity of continuous external monitoring of arterial blood radioactivity, particularly with repeated measurements. An optimized protocol is proposed for routine applications in neurological and neurophysiological studies.
Collapse
Affiliation(s)
- A Bol
- Positron Tomography Laboratory, University of Louvain, Belgium
| | | | | | | | | |
Collapse
|