RESUMEN
Improving sensitivity and spatial resolution in small animal Positron Emission Tomography imaging instrumentation constitutes one of the main goals of nuclear imaging research. These parameters are degraded by the presence of gaps between the detectors. The present manuscript experimentally validates our prototype of an edge-less pre-clinical PET system based on a single LYSO:Ce annulus with an inner diameter of 62 mm and 10 outer facets of 26 × 52 mm2. Scintillation light is read out by arrays of 8 × 8 SiPMs coupled to the facets, using a projection readout of the rows and columns signals. The readout provides accurate Depth of Interaction (DOI). We have implemented a calibration that mitigates the DOI-dependency of the transaxial and axial impact coordinates, and the energy photopeak gain. An energy resolution of 23.4 ± 1.8% was determined. Average spatial resolution of 1.4 ± 0.2 and 1.3 ± 0.4 mm FWHM were achieved for the radial and axial directions, respectively. We found a peak sensitivity of 3.8% at the system center, and a maximum NECR at 40.6 kcps for 0.27 mCi. The image quality was evaluated using reconstructed images of an array of sources and the NEMA image quality phantom was also studied.
RESUMEN
PURPOSE: Human ageing is associated with a regional reduction in cerebral neuronal activity as assessed by numerous studies on brain glucose metabolism and perfusion, grey matter (GM) density and white matter (WM) integrity. As glucose metabolism may impact energetics to maintain myelin integrity, but changes in functional connectivity may also alter regional metabolism, we conducted a cross-sectional simultaneous FDG PET/MR study in a large cohort of healthy volunteers with a wide age range, to directly assess the underlying associations between reduced glucose metabolism, GM atrophy and decreased WM integrity in a single ageing cohort. METHODS: In 94 healthy subjects between 19.9 and 82.5 years (mean 50.1 ± 17.1; 47 M/47F, MMSE ≥ 28), simultaneous FDG-PET, structural MR and diffusion tensor imaging (DTI) were performed. Voxel-wise associations between age and grey matter (GM) density, RBV partial-volume corrected (PVC) glucose metabolism, white matter (WM) fractional anisotropy (FA) and mean diffusivity (MD), and age were assessed. Clusters representing changes in glucose metabolism correlating significantly with ageing were used as seed regions for tractography. Both linear and quadratic ageing models were investigated. RESULTS: An expected age-related reduction in GM density was observed bilaterally in the frontal, lateral and medial temporal cortex, striatum and cerebellum. After PVC, relative FDG uptake was negatively correlated with age in the inferior and midfrontal, cingulate and parietal cortex and subcortical regions, bilaterally. FA decreased with age throughout the entire brain WM. Four white matter tracts were identified connecting brain regions with declining glucose metabolism with age. Within these, relative FDG uptake in both origin and target clusters correlated positively with FA (0.32 ≤ r ≤ 0.71) and negatively with MD (- 0.75 ≤ r ≤ - 0.41). CONCLUSION: After appropriate PVC, we demonstrated that regional cerebral glucose metabolic declines with age and that these changes are related to microstructural changes in the interconnecting WM tracts. The temporal course and potential causality between ageing effects on glucose metabolism and WM integrity should be further investigated in longitudinal cohort PET/MR studies.
Asunto(s)
Envejecimiento , Glucosa , Sustancia Blanca , Envejecimiento/metabolismo , Encéfalo/diagnóstico por imagen , Encéfalo/metabolismo , Estudios Transversales , Imagen de Difusión Tensora/métodos , Glucosa/metabolismo , Sustancia Gris/diagnóstico por imagen , Humanos , Sustancia Blanca/diagnóstico por imagen , Sustancia Blanca/metabolismo , Sustancia Blanca/patologíaRESUMEN
Several research groups are studying organ-dedicated limited angle positron emission tomography (PET) systems to optimize performance-cost ratio, sensitivity, access to the patient and/or flexibility. Often open systems are considered, typically consisting of two detector panels of various sizes. Such systems provide incomplete sampling due to limited angular coverage and/or truncation, which leads to artefacts in the reconstructed activity images. In addition, these organ-dedicated PET systems are usually stand-alone systems, and as a result, no attenuation information can be obtained from anatomical images acquired in the same imaging session. It has been shown that the use of time-of-flight information reduces incomplete data artefacts and enables the joint estimation of the activity and the attenuation factors. In this work, we explore with simple 2D simulations the performance and stability of a joint reconstruction algorithm, for imaging with a limited angle PET system. The reconstruction is based on the so-called MLACF (Maximum Likelihood Attenuation Correction Factors) algorithm and uses linear attenuation coefficients in a known-tissue-class region to obtain absolute quantification. Different panel sizes and different time-of-flight (TOF) resolutions are considered. The noise propagation is compared to that of MLEM reconstruction with exact attenuation correction (AC) for the same PET system. The results show that with good TOF resolution, images of good visual quality can be obtained. If also a good scatter correction can be implemented, quantitative PET imaging will be possible. Further research, in particular on scatter correction, is required.
RESUMEN
A data-driven method is proposed for rigid motion estimation directly from time-of-flight (TOF)-positron emission tomography (PET) emission data. Rigid motion parameters (translations and rotations) are estimated from the first and second moments of the emission data masked in a spherical volume. The accuracy of the method is analyzed on 3D analytical simulations of the PET-SORTEO brain phantom, and subsequently tested on18F-FDG as well as11C-PIB brain datasets acquired on a TOF-PET/CT scanner. The estimated inertia-based motion is later compared to rigid motion parameters obtained by directly registering the short frame backprojections. We find that the method provides sub mm/degree accuracies for the estimated rigid motion parameters for counts corresponding to typical 0.5 s, 1 s, and 2 s18F-FDG brain scans, with the current TOF resolutions clinically available. The method provides robust motion estimation for different types of patient motion, most notably for a continuous patient motion case where conventional frame-based approaches which rely on little to no intra-frame motion of short time intervals could fail. The method relies on the detection of stable eigenvectors for accurate motion estimation, and a monitoring of this condition can reveal time-frames where the motion estimation is less accurate, such as in dynamic PET studies.
Asunto(s)
Procesamiento de Imagen Asistido por Computador , Tomografía Computarizada por Tomografía de Emisión de Positrones , Algoritmos , Encéfalo/diagnóstico por imagen , Humanos , Movimiento (Física) , Movimiento , Fantasmas de Imagen , Tomografía de Emisión de PositronesRESUMEN
In the last two decades, it has been shown that anatomically-guided PET reconstruction can lead to improved bias-noise characteristics in brain PET imaging. However, despite promising results in simulations and first studies, anatomically-guided PET reconstructions are not yet available for use in routine clinical because of several reasons. In light of this, we investigate whether the improvements of anatomically-guided PET reconstruction methods can be achieved entirely in the image domain with a convolutional neural network (CNN). An entirely image-based CNN post-reconstruction approach has the advantage that no access to PET raw data is needed and, moreover, that the prediction times of trained CNNs are extremely fast on state of the art GPUs which will substantially facilitate the evaluation, fine-tuning and application of anatomically-guided PET reconstruction in real-world clinical settings. In this work, we demonstrate that anatomically-guided PET reconstruction using the asymmetric Bowsher prior can be well-approximated by a purely shift-invariant convolutional neural network in image space allowing the generation of anatomically-guided PET images in almost real-time. We show that by applying dedicated data augmentation techniques in the training phase, in which 16 [18F]FDG and 10 [18F]PE2I data sets were used, lead to a CNN that is robust against the used PET tracer, the noise level of the input PET images and the input MRI contrast. A detailed analysis of our CNN in 36 [18F]FDG, 18 [18F]PE2I, and 7 [18F]FET test data sets demonstrates that the image quality of our trained CNN is very close to the one of the target reconstructions in terms of regional mean recovery and regional structural similarity.
Asunto(s)
Encéfalo/diagnóstico por imagen , Procesamiento de Imagen Asistido por Computador/métodos , Redes Neurales de la Computación , Tomografía de Emisión de Positrones , Fluorodesoxiglucosa F18 , Humanos , Imagen por Resonancia Magnética , Nortropanos , Radiofármacos , Tirosina/análogos & derivadosRESUMEN
Positron emission tomography (PET) plays an increasingly important role in research and clinical applications, catalysed by remarkable technical advances and a growing appreciation of the need for reliable, sensitive biomarkers of human function in health and disease. Over the last 30 years, a large amount of the physics and engineering effort in PET has been motivated by the dominant clinical application during that period, oncology. This has led to important developments such as PET/CT, whole-body PET, 3D PET, accelerated statistical image reconstruction, and time-of-flight PET. Despite impressive improvements in image quality as a result of these advances, the emphasis on static, semi-quantitative 'hot spot' imaging for oncologic applications has meant that the capability of PET to quantify biologically relevant parameters based on tracer kinetics has not been fully exploited. More recent advances, such as PET/MR and total-body PET, have opened up the ability to address a vast range of new research questions, from which a future expansion of applications and radiotracers appears highly likely. Many of these new applications and tracers will, at least initially, require quantitative analyses that more fully exploit the exquisite sensitivity of PET and the tracer principle on which it is based. It is also expected that they will require more sophisticated quantitative analysis methods than those that are currently available. At the same time, artificial intelligence is revolutionizing data analysis and impacting the relationship between the statistical quality of the acquired data and the information we can extract from the data. In this roadmap, leaders of the key sub-disciplines of the field identify the challenges and opportunities to be addressed over the next ten years that will enable PET to realise its full quantitative potential, initially in research laboratories and, ultimately, in clinical practice.
Asunto(s)
Inteligencia Artificial , Neoplasias/diagnóstico por imagen , Tomografía Computarizada por Tomografía de Emisión de Positrones/métodos , Tomografía Computarizada por Tomografía de Emisión de Positrones/tendencias , Tomografía de Emisión de Positrones/métodos , Tomografía de Emisión de Positrones/tendencias , Historia del Siglo XX , Historia del Siglo XXI , Humanos , Procesamiento de Imagen Asistido por Computador/métodos , Imagenología Tridimensional , Cinética , Oncología Médica/métodos , Oncología Médica/tendencias , Tomografía Computarizada por Tomografía de Emisión de Positrones/historia , Pronóstico , Radiofármacos , Biología de Sistemas , Tomografía Computarizada por Rayos XRESUMEN
This study evaluates the performance of the Bruker positron emission tomograph (PET) insert combined with a BioSpec 70/30 USR magnetic resonance imaging (MRI) scanner using the manufacturer acceptance protocol and the NEMA NU 4-2008 for small animal PET. The PET insert is made of 3 rings of 8 monolithic LYSO crystals (50 × 50 × 10 mm3) coupled to silicon photomultipliers (SiPM) arrays, conferring an axial and transaxial FOV of 15 cm and 8 cm. The MRI performance was evaluated with and without the insert for the following radiofrequency noise, magnetic field homogeneity and image quality. For the PET performance, we extended the NEMA protocol featuring system sensitivity, count rates, spatial resolution and image quality to homogeneity and accuracy for quantification using several MRI sequences (RARE, FLASH, EPI and UTE). The PET insert does not show any adverse effect on the MRI performances. The MR field homogeneity is well preserved (Diameter Spherical Volume, for 20 mm of 1.98 ± 4.78 without and -0.96 ± 5.16 Hz with the PET insert). The PET insert has no major effect on the radiofrequency field. The signal-to-noise ratio measurements also do not show major differences. Image ghosting is well within the manufacturer specifications (<2.5%) and no RF noise is visible. Maximum sensitivity of the PET insert is 11.0% at the center of the FOV even with simultaneous acquisition of EPI and RARE. PET MLEM resolution is 0.87 mm (FWHM) at 5 mm off-center of the FOV and 0.97 mm at 25 mm radial offset. The peaks for true/noise equivalent count rates are 410/240 and 628/486 kcps for the rat and mouse phantoms, and are reached at 30.34/22.85 and 27.94/22.58 MBq. PET image quality is minimally altered by the different MRI sequences. The Bruker PET insert shows no adverse effect on the MRI performance and demonstrated a high sensitivity, sub-millimeter resolution and good image quality even during simultaneous MRI acquisition.
Asunto(s)
Imagen por Resonancia Magnética/instrumentación , Tomografía de Emisión de Positrones/instrumentación , Animales , Diseño de Equipo , Modelos Lineales , Ratones , Fantasmas de Imagen , Ratas , Relación Señal-RuidoRESUMEN
BACKGROUND: Yoga is increasingly popular worldwide with several physical and mental benefits, but the underlying neurobiology remains unclear. Whereas many studies have focused on pure meditational aspects, the triad of yoga includes meditation, postures, and breathing. We conducted a cross-sectional study comparing experienced yoga practitioners to yoga-naive healthy subjects using a multiparametric 2 × 2 design with simultaneous positron emission tomography/magnetic resonance (PET/MR) imaging. METHODS: 18F-FDG PET, morphometric and diffusion tensor imaging, resting state fMRI, and MR spectroscopy were acquired in 10 experienced (4.8 ± 2.3 years of regular yoga experience) yoga practitioners and 15 matched controls in rest and after a single practice (yoga practice and physical exercise, respectively). RESULTS: In rest, decreased regional glucose metabolism in the medial temporal cortex, striatum, and brainstem was observed in yoga practitioners compared to controls (p < 0.0001), with a significant inverse correlation of resting parahippocampal and brainstem metabolism with years of regular yoga practice (ρ < - 0.63, p < 0.05). A single yoga practice resulted in significant hypermetabolism in the cerebellum (p < 0.0001). None of the MR measures differed, both at rest and after intervention. CONCLUSIONS: Experienced yoga practitioners show regional long-term decreases in glucose metabolism related to years of practice. To elucidate a potential causality, a prospective longitudinal study in yoga-naive individuals is warranted.
RESUMEN
Importance: During a time with the potential for novel treatment strategies, early detection of disease manifestations at an individual level in presymptomatic carriers of a hexanucleotide repeat expansion in the C9orf72 gene (preSxC9) is becoming increasingly relevant. Objectives: To evaluate changes in glucose metabolism before symptom onset of amyotrophic lateral sclerosis or frontotemporal dementia in preSxC9 using simultaneous fluorine 18-labeled fluorodeoxyglucose ([18F]FDG positron emission tomographic (PET) and magnetic resonance imaging as well as the mutation's association with clinical and fluid biomarkers. Design, Setting, and Participants: A prospective, case-control study enrolled 46 participants from November 30, 2015, until December 11, 2018. The study was conducted at the neuromuscular reference center of the University Hospitals Leuven, Leuven, Belgium. Main Outcomes and Measures: Neuroimaging data were spatially normalized and analyzed at the voxel level at a height threshold of P < .001, cluster-level familywise error-corrected threshold of P < .05, and statistical significance was set at P < .05 for the volume-of-interest level analysis, using Benjamini-Hochberg correction for multiple correction. W-score maps were computed using the individuals serving as controls as a reference to quantify the degree of [18F]FDG PET abnormality. The threshold for abnormality on the W-score maps was designated as an absolute W-score greater than or equal to 1.96. Neurofilament levels and performance on cognitive and neurologic examinations were determined. All hypothesis tests were 1-sided. Results: Of the 42 included participants, there were 17 with the preSxC9 mutation (12 women [71%]; mean [SD] age, 51 [9] years) and 25 healthy controls (12 women [48%]; mean [SD] age, 47 [10] years). Compared with control participants, significant clusters of relative hypometabolism were found in frontotemporal regions, basal ganglia, and thalami of preSxC9 participants and relative hypermetabolism in the peri-Rolandic region, superior frontal gyrus, and precuneus cortex. W-score frequency maps revealed reduced glucose metabolism with local maxima in the insular cortices, central opercular cortex, and thalami in up to 82% of preSxC9 participants and increased glucose metabolism in the precentral gyrus and precuneus cortex in up to 71% of preSxC9 participants. Other findings in the preSxC9 group were upper motor neuron involvement in 10 participants (59%), cognitive abnormalities in 5 participants (29%), and elevated neurofilament levels in 3 of 16 individuals (19%) who underwent lumbar puncture. Conclusions and Relevance: The results suggest that [18F]FDG PET can identify glucose metabolic changes in preSxC9 at an individual level, preceding significantly elevated neurofilament levels and onset of symptoms.
Asunto(s)
Esclerosis Amiotrófica Lateral/diagnóstico , Proteína C9orf72/genética , Corteza Cerebral/diagnóstico por imagen , Demencia Frontotemporal/diagnóstico , Tomografía de Emisión de Positrones/normas , Síntomas Prodrómicos , Adulto , Anciano , Esclerosis Amiotrófica Lateral/genética , Esclerosis Amiotrófica Lateral/metabolismo , Esclerosis Amiotrófica Lateral/patología , Biomarcadores , Estudios de Casos y Controles , Corteza Cerebral/metabolismo , Corteza Cerebral/patología , Expansión de las Repeticiones de ADN , Femenino , Fluorodesoxiglucosa F18 , Demencia Frontotemporal/genética , Demencia Frontotemporal/metabolismo , Demencia Frontotemporal/patología , Heterocigoto , Humanos , Imagen por Resonancia Magnética , Masculino , Persona de Mediana Edad , Imagen Multimodal , Estudios ProspectivosRESUMEN
Neuroinflammation has been implicated in amyotrophic lateral sclerosis (ALS) and can be visualized using translocator protein (TSPO) radioligands. To become a reliable pharmacodynamic biomarker for ALS multicenter trials, TSPO radioligands have some challenges to overcome. We aimed to investigate whether multicenter data pooling of different TSPO tracers (11C-PBR28 and 18F-DPA714) is feasible, after validation of an established 11C-PBR28 PET pseudo reference analysis technique for 18F-DPA714. Methods: Seven ALS patients from Belgium (58.9 ± 6.7 y old, 5 men and 2 women), 8 healthy volunteers from Belgium (52.1 ± 15.2 y old, 3 men and 5 women), 7 ALS patients from the United States (53.4 ± 9.8 y old, 5 men and 2 women), and 7 healthy volunteers from the United States (54.6 ± 9.6 y old, 4 men and 3 women) from a previously published study underwent dynamic 18F-DPA714 (Leuven, Belgium) or 11C-PBR28 (Boston, Massachusetts) PET/MRI. For 18F-DPA714, maps of total volume of distribution (VT) were compared with SUV ratio (SUVR) images from 40 to 60 min after injection (SUVR40-60) calculated using the pseudo reference regions cerebellum, occipital cortex, and whole brain (WB) without ventricles. For 11C-PBR28, SUVR images from 60 to 90 min after injection using the WB without ventricles were calculated. Results: In line with previous studies, increased 18F-DPA714 uptake (17.0% ± 5.6%) in primary motor cortices was observed in ALS subjects, as measured by both VT and SUVR40-60 approaches. The highest sensitivity was found for SUVR calculated using the WB without ventricles (average cluster, 21.6% ± 0.1%). 18F-DPA714 VT ratio was highly correlated with the SUVR40-60 (r > 0.8, P < 0.001). A similar pattern of increased uptake (average cluster, 20.5% ± 0.5%) in the primary motor cortices was observed in ALS subjects for 11C-PBR28 SUVR calculated using the WB without ventricles. Analysis of the 18F-DPA714 and 11C-PBR28 data together resulted in a more extensive pattern of significantly increased glial activation bilaterally in the primary motor cortices. Conclusion: The same pseudo reference region analysis technique for 11C-PBR28 PET can be extended toward 18F-DPA714 PET. Therefore, in ALS, standardized analysis across these 2 tracers enables pooling of TSPO PET data across multiple centers and increases the power of TSPO as a biomarker for future therapeutic trials.
Asunto(s)
Acetamidas/farmacocinética , Esclerosis Amiotrófica Lateral/diagnóstico por imagen , Radioisótopos de Carbono/farmacocinética , Radioisótopos de Flúor/farmacocinética , Tomografía de Emisión de Positrones/métodos , Pirazoles/farmacocinética , Piridinas/farmacocinética , Pirimidinas/farmacocinética , Receptores de GABA/metabolismo , Adulto , Anciano , Esclerosis Amiotrófica Lateral/tratamiento farmacológico , Ensayos Clínicos como Asunto , Estudios de Factibilidad , Femenino , Humanos , Masculino , Persona de Mediana EdadRESUMEN
With increasing improvements in the time of flight (TOF) resolution of positron emission tomography (PET) scanners, an accurate model of the TOF measurements is becoming increasingly important. This work considers two parameters of the TOF kernel; the relative positioning of the timing data-bins and the timing resolution along each line of response (LOR). Similar to an existing data-driven method, we assume that any shifts of data-bins along lines of response can be modelled as differences between crystal timing offsets. Inspired by this, timing resolutions of all LORs are modelled as the hypotenuse of timing resolutions of the crystal-pairs in coincidence. Furthermore, in order to mitigate the influence of potential inaccuracies of detector-pair sensitivities on crystal timing resolutions, relative LOR sensitivities are modelled as the product of efficiency factors for the two crystals in coincidence. We validate estimating maps of crystal timing offsets, timing resolutions and efficiencies from the emission data using noisy simulations of a brain phantom. Results are shown for phantom and patient data scanned on clinically available TOF-PET scanners. We find that the estimation of crystal timing resolutions is more sensitive to the data statistics than the estimation of crystal timing offsets. As a result, estimation of crystal timing properties could either be limited to high count emission data, or be obtained utilizing additional regularizations on the estimates. Using a more accurate model of the TOF acquisition, improvements are observed in standard activity reconstructions as well as joint reconstructions of activity and attenuation.
Asunto(s)
Encéfalo/diagnóstico por imagen , Encéfalo/fisiología , Procesamiento de Imagen Asistido por Computador/métodos , Tomografía de Emisión de Positrones/métodos , Humanos , Modelos Estadísticos , Fantasmas de ImagenRESUMEN
Time-of-flight (TOF) PET data provide an effective means for attenuation correction (AC) when no (or incomplete or inaccurate) attenuation information is available. Since MR scanners provide little information on photon attenuation of different tissue types, AC in hybrid PET/MR scanners has always been challenging. In this contribution, we aim at validating the activity reconstructions of the maximum-likelihood ordered-subsets activity and attenuation (OSAA) reconstruction algorithm on a patient brain data set. We present a quantitative comparison of joint reconstructions with the current clinical gold standard-ordered-subsets expectation maximization-using CT-based AC in PET/CT, as well as the current state of the art in PET/MR, that is, zero time echo (ZTE)-based AC. Methods: The TOF PET emission data were initially used in a preprocessing stage to estimate crystal maps of efficiencies, timing offsets, and timing resolutions. Applying these additional corrections during reconstructions, OSAA, ZTE-based, and the vendor-provided atlas-based AC techniques were analyzed and compared with CT-based AC. In our initial study, we used the CT-based estimate of the expected scatter and later used the ZTE-based and OSAA attenuation estimates to compute the expected scatter contribution of the data during reconstructions. In all reconstructions, a maximum-likelihood scaling of the single-scatter simulation estimate to the emission data was used for scatter correction. The reconstruction results were analyzed in the 86 segmented regions of interest of the Hammers atlas. Results: Our quantitative analysis showed that, in practice, a tracer activity difference of +0.5% (±2.1%) and +0.1% (±2.3%) could be expected for the state-of-the-art ZTE-based and OSAA AC methods, respectively, in PET/MR compared with the clinical gold standard in PET/CT. Conclusion: Joint activity and attenuation estimation methods can provide an effective solution to the challenging AC problem for brain studies in hybrid TOF PET/MR scanners. With an accurate TOF-based (timing offsets and timing resolutions) calibration, and similar to the results of the state-of-the-art method in PET/MR, regional errors of joint TOF PET reconstructions are within a few percentage points.
Asunto(s)
Encéfalo/diagnóstico por imagen , Procesamiento de Imagen Asistido por Computador/métodos , Imagen por Resonancia Magnética , Imagen Multimodal , Tomografía de Emisión de Positrones , Humanos , Factores de TiempoRESUMEN
There are drawbacks with using a Positron Emission Tomography (PET) scanner design employing the traditional arrangement of multiple detectors in an array format. Typically PET systems are constructed with many regular gaps between the detector modules in a ring or box configuration, with additional axial gaps between the rings. Although this has been significantly reduced with the use of the compact high granularity SiPM photodetector technology, such a scanner design leads to a decrease in the number of annihilation photons that are detected causing lower scanner sensitivity. Moreover, the ability to precisely determine the line of response (LOR) along which the positron annihilated is diminished closer to the detector edges because the spatial resolution there is degraded due to edge effects. This happens for both monolithic based designs, caused by the truncation of the scintillation light distribution, but also for detector blocks that use crystal arrays with a number of elements that are larger than the number of photosensors and, therefore, make use of the light sharing principle. In this report we present a design for a small-animal PET scanner based on a single monolithic annulus-like scintillator that can be used as a PET insert in high-field Magnetic Resonance systems. We provide real data showing the performance improvement when edge-less modules are used. We also describe the specific proposed design for a rodent scanner that employs facetted outside faces in a single LYSO tube. In a further step, in order to support and prove the proposed edgeless geometry, simulations of that scanner have been performed and lately reconstructed showing the advantages of the design.
RESUMEN
This paper presents new data driven methods for the time of flight (TOF) calibration of positron emission tomography (PET) scanners. These methods are derived from the consistency condition for TOF PET, they can be applied to data measured with an arbitrary tracer distribution and are numerically efficient because they do not require a preliminary image reconstruction from the non-TOF data. Two-dimensional simulations are presented for one of the methods, which only involves the two first moments of the data with respect to the TOF variable. The numerical results show that this method estimates the detector timing offsets with errors that are larger than those obtained via an initial non-TOF reconstruction, but remain smaller than [Formula: see text] of the TOF resolution and thereby have a limited impact on the quantitative accuracy of the activity image estimated with standard maximum likelihood reconstruction algorithms.
Asunto(s)
Algoritmos , Procesamiento de Imagen Asistido por Computador/métodos , Tomografía de Emisión de Positrones/métodos , Calibración , Humanos , Factores de TiempoRESUMEN
Methods for joint activity reconstruction and attenuation reconstruction of time-of-flight (TOF) PET data provide an effective solution to attenuation correction when no (or incomplete or inaccurate) information on attenuation is available. One of the main barriers limiting use of these methods in clinical practice is their lack of validation in a relatively large patient database. In this contribution, we aim to validate reconstruction performed with maximum-likelihood activity reconstruction and attenuation registration (MLRR) in a whole-body patient dataset. Furthermore, a partial validation (because the scale problem of the algorithm is avoided for now) of reconstruction performed with maximum-likelihood activity and attenuation (MLAA) is also provided. We present a quantitative comparison between these 2 methods of joint reconstruction and the current clinical gold standard, maximum-likelihood expectation maximization (MLEM) with CT-based attenuation correction. Methods: The whole-body TOF PET emission data of each patient dataset were processed as a whole to reconstruct an activity volume covering all the acquired bed positions, helping reduce the problem of a scale per bed position in MLAA to a global scale for the entire activity volume. Three reconstruction algorithms were used: MLEM, MLRR, and MLAA. A maximum-likelihood scaling of the single-scatter simulation estimate to the emission data was used for scatter correction. The reconstruction results for various regions of interest were then analyzed. Results: The joint reconstructions of the whole-body patient dataset provided better quantification than the gold standard in cases of PET and CT misalignment caused by patient or organ motion. Our quantitative analysis showed a difference of -4.2% ± 2.3% between MLRR and MLEM and a difference of -7.5% ± 4.6% between MLAA and MLEM, averaged over all regions of interest. Conclusion: Joint reconstruction of activity and attenuation provides a useful means to estimate tracer distribution when CT-based-attenuation images are subject to misalignment or are not available. With an accurate estimate of the scatter contribution in the emission measurements, the joint reconstructions of TOF PET data are within clinically acceptable accuracy.
Asunto(s)
Procesamiento de Imagen Asistido por Computador/métodos , Tomografía de Emisión de Positrones , Algoritmos , Fluorodesoxiglucosa F18 , Humanos , Factores de Tiempo , Imagen de Cuerpo EnteroRESUMEN
In this article, we evaluate Parallel Level Sets (PLS) and Bowsher's method as segmentation-free anatomical priors for regularized brain positron emission tomography (PET) reconstruction. We derive the proximity operators for two PLS priors and use the EM-TV algorithm in combination with the first order primal-dual algorithm by Chambolle and Pock to solve the non-smooth optimization problem for PET reconstruction with PLS regularization. In addition, we compare the performance of two PLS versions against the symmetric and asymmetric Bowsher priors with quadratic and relative difference penalty function. For this aim, we first evaluate reconstructions of 30 noise realizations of simulated PET data derived from a real hybrid positron emission tomography/magnetic resonance imaging (PET/MR) acquisition in terms of regional bias and noise. Second, we evaluate reconstructions of a real brain PET/MR data set acquired on a GE Signa time-of-flight PET/MR in a similar way. The reconstructions of simulated and real 3D PET/MR data show that all priors were superior to post-smoothed maximum likelihood expectation maximization with ordered subsets (OSEM) in terms of bias-noise characteristics in different regions of interest where the PET uptake follows anatomical boundaries. Our implementation of the asymmetric Bowsher prior showed slightly superior performance compared with the two versions of PLS and the symmetric Bowsher prior. At very high regularization weights, all investigated anatomical priors suffer from the transfer of non-shared gradients.
Asunto(s)
Encéfalo/diagnóstico por imagen , Procesamiento de Imagen Asistido por Computador/métodos , Tomografía de Emisión de Positrones/métodos , Algoritmos , Humanos , Imagen por Resonancia Magnética , Fantasmas de ImagenRESUMEN
Scatter correction is typically done using a simulation of the single scatter, which is then scaled to account for multiple scatters and other possible model mismatches. This scaling factor is determined by fitting the simulated scatter sinogram to the measured sinogram, using only counts measured along LORs that do not intersect the patient body, i.e. 'scatter-tails'. Extending previous work, we propose to scale the scatter with a plane dependent factor, which is determined as an additional unknown in the maximum likelihood (ML) reconstructions, using counts in the entire sinogram rather than only the 'scatter-tails'. The ML-scaled scatter estimates are validated using a Monte-Carlo simulation of a NEMA-like phantom, a phantom scan with typical contrast ratios of a 68Ga-PSMA scan, and 23 whole-body 18F-FDG patient scans. On average, we observe a 12.2% change in the total amount of tracer activity of the MLEM reconstructions of our whole-body patient database when the proposed ML scatter scales are used. Furthermore, reconstructions using the ML-scaled scatter estimates are found to eliminate the typical 'halo' artifacts that are often observed in the vicinity of high focal uptake regions.
Asunto(s)
Imagenología Tridimensional/métodos , Método de Montecarlo , Fantasmas de Imagen , Tomografía Computarizada por Tomografía de Emisión de Positrones/métodos , Radiofármacos , Dispersión de Radiación , Recuento Corporal Total/métodos , Fluorodesoxiglucosa F18 , Humanos , Procesamiento de Imagen Asistido por ComputadorRESUMEN
For quantitative tracer distribution in positron emission tomography, attenuation correction is essential. In a hybrid PET/CT system the CT images serve as a basis for generation of the attenuation map, but in PET/MR, the MR images do not have a similarly simple relationship with the attenuation map. Hence attenuation correction in PET/MR systems is more challenging. Typically either of two MR sequences are used: the Dixon or the ultra-short time echo (UTE) techniques. However these sequences have some well-known limitations. In this study, a reconstruction technique based on a modified and optimized non-TOF MLAA is proposed for PET/MR brain imaging. The idea is to tune the parameters of the MLTR applying some information from an attenuation image computed from the UTE sequences and a T1w MR image. In this MLTR algorithm, an [Formula: see text] parameter is introduced and optimized in order to drive the algorithm to a final attenuation map most consistent with the emission data. Because the non-TOF MLAA is used, a technique to reduce the cross-talk effect is proposed. In this study, the proposed algorithm is compared to the common reconstruction methods such as OSEM using a CT attenuation map, considered as the reference, and OSEM using the Dixon and UTE attenuation maps. To show the robustness and the reproducibility of the proposed algorithm, a set of 204 [18F]FDG patients, 35 [11C]PiB patients and 1 [18F]FET patient are used. The results show that by choosing an optimized value of [Formula: see text] in MLTR, the proposed algorithm improves the results compared to the standard MR-based attenuation correction methods (i.e. OSEM using the Dixon or the UTE attenuation maps), and the cross-talk and the scale problem are limited.
Asunto(s)
Encéfalo/diagnóstico por imagen , Procesamiento de Imagen Asistido por Computador/métodos , Imagen por Resonancia Magnética , Imagen Multimodal , Tomografía de Emisión de Positrones , Algoritmos , Fluorodesoxiglucosa F18 , Humanos , Reproducibilidad de los Resultados , Factores de TiempoRESUMEN
Previously, maximum-likelihood methods have been proposed to jointly estimate the activity image and the attenuation image or the attenuation sinogram from time-of-flight (TOF) positron emission tomography (PET) data. In this contribution, we propose a method that addresses the possible alignment problem of the TOF-PET emission data and the computed tomography (CT) attenuation data, by combining reconstruction and registration. The method, called MLRR, iteratively reconstructs the activity image while registering the available CT-based attenuation image, so that the pair of activity and attenuation images maximise the likelihood of the TOF emission sinogram. The algorithm is slow to converge, but some acceleration could be achieved by using Nesterov's momentum method and by applying a multi-resolution scheme for the non-rigid displacement estimation. The latter also helps to avoid local optima, although convergence to the global optimum cannot be guaranteed. The results are evaluated on 2D and 3D simulations as well as a respiratory gated clinical scan. Our experiments indicate that the proposed method is able to correct for possible misalignment of the CT-based attenuation image, and is therefore a very promising approach to suppressing attenuation artefacts in clinical PET/CT. When applied to respiratory gated data of a patient scan, it produced deformations that are compatible with breathing motion and which reduced the well known attenuation artefact near the dome of the liver. Since the method makes use of the energy-converted CT attenuation image, the scale problem of joint reconstruction is automatically solved.