We present a standalone, scalable and high-throughput software platform for PET image reconstruction and analysis. We focus on high fidelity modelling of the acquisition processes to provide high ...accuracy and precision quantitative imaging, especially for large axial field of view scanners. All the core routines are implemented using parallel computing available from within the Python package
NiftyPET
, enabling easy access, manipulation and visualisation of data at any processing stage. The pipeline of the platform starts from MR and raw PET input data and is divided into the following processing stages: (1) list-mode data processing; (2) accurate attenuation coefficient map generation; (3) detector normalisation; (4) exact forward and back projection between sinogram and image space; (5) estimation of reduced-variance random events; (6) high accuracy fully 3D estimation of scatter events; (7) voxel-based partial volume correction; (8) region- and voxel-level image analysis. We demonstrate the advantages of this platform using an amyloid brain scan where all the processing is executed from a single and uniform computational environment in Python. The high accuracy acquisition modelling is achieved through span-1 (no axial compression) ray tracing for true, random and scatter events. Furthermore, the platform offers uncertainty estimation of any image derived statistic to facilitate robust tracking of subtle physiological changes in longitudinal studies. The platform also supports the development of new reconstruction and analysis algorithms through restricting the axial field of view to any set of rings covering a region of interest and thus performing fully 3D reconstruction and corrections using real data significantly faster. All the software is available as open source with the accompanying wiki-page and test data.
Patient respiratory motion during PET image acquisition leads to blurring in the reconstructed images and may cause significant artifacts, resulting in decreased lesion detectability, inaccurate ...standard uptake value calculation and incorrect treatment planning in radiation therapy. To reduce these effects data can be regrouped into (nearly) 'motion-free' gates prior to reconstruction by selecting the events with respect to the breathing phase. This gating procedure therefore needs a respiratory signal: on current scanners it is obtained from an external device, whereas with data driven (DD) methods it can be directly obtained from the raw PET data. DD methods thus eliminate the use of external equipment, which is often expensive, needs prior setup and can cause patient discomfort, and they could also potentially provide increased fidelity to the internal movement. DD methods have been recently applied on PET data showing promising results. However, many methods provide signals whose direction with respect to the physical motion is uncertain (i.e. their sign is arbitrary), therefore a maximum in the signal could refer either to the end-inspiration or end-expiration phase, possibly causing inaccurate motion correction. In this work we propose two novel methods, CorrWeights and CorrSino, to detect the correct direction of the motion represented by the DD signal, that is obtained by applying principal component analysis (PCA) on the acquired data. They only require the PET raw data, and they rely on the assumption that one of the major causes of change in the acquired data related to the chest is respiratory motion in the axial direction, that generates a cranio-caudal motion of the internal organs. We also implemented two versions of a published registration-based method, that require image reconstruction. The methods were first applied on XCAT simulations, and later evaluated on cancer patient datasets monitored by the Varian Real-time Position ManagementTM (RPM) device, selecting the lower chest bed positions. For each patient different time intervals were evaluated ranging from 50 to 300 s in duration. The novel methods proved to be generally more accurate than the registration-based ones in detecting the correct sign of the respiratory signal, and their failure rates are lower than 3% when the DD signal is highly correlated with the RPM. They also have the advantage of faster computation time, avoiding reconstruction. Moreover, CorrWeights is not specifically related to PCA and considering its simple implementation, it could easily be applied together with any DD method in clinical practice.
Patient motion due to respiration can lead to artefacts and blurring in positron emission tomography (PET) images, in addition to quantification errors. The integration of PET with magnetic resonance ...(MR) imaging in PET-MR scanners provides complementary clinical information, and allows the use of high spatial resolution and high contrast MR images to monitor and correct motion-corrupted PET data. In this paper we build on previous work to form a methodology for respiratory motion correction of PET data, and show it can improve PET image quality whilst having minimal impact on clinical PET-MR protocols. We introduce a joint PET-MR motion model, using only 1 min per PET bed position of simultaneously acquired PET and MR data to provide a respiratory motion correspondence model that captures inter-cycle and intra-cycle breathing variations. In the model setup, 2D multi-slice MR provides the dynamic imaging component, and PET data, via low spatial resolution framing and principal component analysis, provides the model surrogate. We evaluate different motion models (1D and 2D linear, and 1D and 2D polynomial) by computing model-fit and model-prediction errors on dynamic MR images on a data set of 45 patients. Finally we apply the motion model methodology to 5 clinical PET-MR oncology patient datasets. Qualitative PET reconstruction improvements and artefact reduction are assessed with visual analysis, and quantitative improvements are calculated using standardised uptake value (SUVpeak and SUVmax) changes in avid lesions. We demonstrate the capability of a joint PET-MR motion model to predict respiratory motion by showing significantly improved image quality of PET data acquired before the motion model data. The method can be used to incorporate motion into the reconstruction of any length of PET acquisition, with only 1 min of extra scan time, and with no external hardware required.
The last decade has seen a great proliferation of supervised learning pipelines for individual diagnosis and prognosis in Alzheimer's disease. As more pipelines are developed and evaluated in the ...search for greater performance, only those results that are relatively impressive will be selected for publication. We present an empirical study to evaluate the potential for optimistic bias in classification performance results as a result of this selection. This is achieved using a novel, resampling-based experiment design that effectively simulates the optimisation of pipeline specifications by individuals or collectives of researchers using cross validation with limited data. Our findings indicate that bias can plausibly account for an appreciable fraction (often greater than half) of the apparent performance improvement associated with the pipeline optimisation, particularly in small samples. We discuss the consistency of our findings with patterns observed in the literature and consider strategies for bias reduction and mitigation.
Purpose
Patients with idiopathic pulmonary fibrosis (IPF) show increased PET signal at sites of morphological abnormality on high-resolution computed tomography (HRCT). The purpose of this ...investigation was to investigate the PET signal at sites of normal-appearing lung on HRCT in IPF.
Methods
Consecutive IPF patients (22 men, 3 women) were prospectively recruited. The patients underwent
18
F-FDG PET/HRCT. The pulmonary imaging findings in the IPF patients were compared to the findings in a control population. Pulmonary uptake of
18
F-FDG (mean SUV) was quantified at sites of morphologically normal parenchyma on HRCT. SUVs were also corrected for tissue fraction (TF). The mean SUV in IPF patients was compared with that in 25 controls (patients with lymphoma in remission or suspected paraneoplastic syndrome with normal PET/CT appearances).
Results
The pulmonary SUV (mean ± SD) uncorrected for TF in the controls was 0.48 ± 0.14 and 0.78 ± 0.24 taken from normal lung regions in IPF patients (
p
< 0.001). The TF-corrected mean SUV in the controls was 2.24 ± 0.29 and 3.24 ± 0.84 in IPF patients (
p
< 0.001).
Conclusion
IPF patients have increased pulmonary uptake of
18
F-FDG on PET in areas of lung with a normal morphological appearance on HRCT. This may have implications for determining disease mechanisms and treatment monitoring.
Pharmacokinetic modelling on dynamic positron emission tomography (PET) data is a quantitative technique. However, the long acquisition time is prohibitive for routine clinical use. Instead, the ...semi-quantitative standardised uptake value ratio (SUVR) from a shorter static acquisition is used, despite its sensitivity to blood flow confounding longitudinal analysis. A method has been proposed to reduce the dynamic acquisition time for quantification by incorporating cerebral blood flow (CBF) information from arterial spin labelling (ASL) magnetic resonance imaging (MRI) into the pharmacokinetic modelling. In this work, we optimise and validate this framework for a study of ageing and preclinical Alzheimer's disease. This methodology adapts the simplified reference tissue model (SRTM) for a reduced acquisition time (RT-SRTM) and is applied to 18F-florbetapir PET data for amyloid-β quantification. Evaluation shows that the optimised RT-SRTM can achieve amyloid burden estimation from a 30-min PET/MR acquisition which is comparable with the gold standard SRTM applied to 60 min of PET data. Conversely, SUVR showed a significantly higher error and bias, and a statistically significant correlation with tracer delivery due to the influence of blood flow. The optimised RT-SRTM produced amyloid burden estimates which were uncorrelated with tracer delivery indicating its suitability for longitudinal studies.
In this study, we explore the use of a spatially-variant penalty strength in penalized image reconstruction using anatomical priors to reduce the dependence of lesion contrast on surrounding activity ...and lesion location. This work builds on a previous method to make the local perturbation response (LPR) approximately spatially invariant. While the dependence of lesion contrast on the local properties introduced by the anatomical penalty is intentional, the method aims to reduce the influence from surroundings lying along the lines of response (LORs) but not in the penalty neighborhood structure. The method is evaluated using simulated data, assuming that the anatomical information is absent or well-aligned with the corresponding activity images. Since the parallel level sets (PLS) penalty is convex and has shown promising results in the literature, it is chosen as the representative anatomical penalty and incorporated into the previously proposed preconditioned algorithm (L-BFGS-B-PC) for achieving good image quality and fast convergence rate. A 2D disc phantom with a feature at the center and a 3D XCAT thorax phantom with lesions inserted in different slices are used to study how surrounding activity and lesion location affect the visual appearance and quantitative consistency. A bias and noise analysis is also performed with the 2D disc phantom. The consistency of the algorithm convergence rate with respect to different data noise and background levels is also investigated using the XCAT phantom. Finally, an example of reconstruction for a patient dataset with inserted pseudo lesions is used as a demonstration in a clinical context. We show that applying the spatially-variant penalization with PLS can reduce the dependence of the lesion contrast on the surrounding activity and lesion location. It does not affect the bias and noise trade-off curves for matched local resolution. Moreover, when using the proposed penalization, significant improvement in algorithm convergence rate and convergence consistency is observed.
Abstract
Statistical iterative reconstruction is now widely used in clinical practice and has contributed to significant improvement in image quality in recent years. Although primarily used for ...reconstruction in emission tomography (both single photon emission computed tomography (SPECT) and positron emission tomography (PET)) there is increasing interest in also applying similar algorithms to x-ray computed tomography (CT). There is increasing complexity in the factors that are included in the reconstruction, a demonstration of the versatility of the approach. Research continues with exploration of methods for further improving reconstruction quality with effective correction for various sources of artefact.
This work is an extension of our recent work on joint activity reconstruction/motion estimation (JRM) from positron emission tomography (PET) data. We performed JRM by maximization of the penalized ...log-likelihood in which the probabilistic model assumes that the same motion field affects both the activity distribution and the attenuation map. Our previous results showed that JRM can successfully reconstruct the activity distribution when the attenuation map is misaligned with the PET data, but converges slowly due to the significant cross-talk in the likelihood. In this paper, we utilize time-of-flight PET for JRM and demonstrate that the convergence speed is significantly improved compared to JRM with conventional PET data.
This paper reports on the feasibility of using a quasi-Newton optimization algorithm, limited-memory Broyden-Fletcher-Goldfarb-Shanno with boundary constraints (L-BFGS-B), for penalized image ...reconstruction problems in emission tomography (ET). For further acceleration, an additional preconditioning technique based on a diagonal approximation of the Hessian was introduced. The convergence rate of L-BFGS-B and the proposed preconditioned algorithm (L-BFGS-B-PC) was evaluated with simulated data with various factors, such as the noise level, penalty type, penalty strength and background level. Data of three 18 F-FDG patient acquisitions were also reconstructed. Results showed that the proposed L-BFGS-B-PC outperforms L-BFGS-B in convergence rate for all simulated conditions and the patient data. Based on these results, L-BFGS-B-PC shows promise for clinical application.