Uniform acquisition modelling across PET imaging systems: Unified scatter modelling

© 2016 IEEE. PET imaging is an important tool commonly used for studying disease by research consortia which implement multi-centre studies to improve the statistical power of findings. The UK government launched the Dementias Platform UK to facilitate one of the world's largest dementia population study involving national centres equipped with state-of-the-art PET/MR scanners from two major vendors. However, the difference in PET detector technology between the two scanners involved makes the standardisation of data acquisition and image reconstruction necessary. We propose a new approach to PET acquisition system modelling across different PET systems and technologies, focusing in particular on unified scatter estimation across TOF (time-of-flight) and non-TOF PET systems. The proposed scatter modelling is fully 3D and voxel based, as opposed to the popular line-of-response driven methods. This means that for each emitting voxel an independent 3D scatter estimate is found, inherently preserving the necessary information for TOF calculations as well as accounting for the large axial field of view. With adequate sampling of the input images, the non-TOF scatter estimate is identical to the summed TOF estimates across TOF bins, without an additional computational cost used for the TOF estimation. The model is implemented using the latest NVIDA GPU CUDA platform, allowing finer sampling of image space which is more essential for accurate TOF modelling. The high accuracy of the proposed scatter model is validated using Monte Carlo simulations. The model is deployed in our stand-alone image reconstruction pipeline for the Biograph mMR scanner, demonstrating accurate 3D scatter estimates resulting in uniform reconstruction for a high statistics phantom scan

Uniform acquisition modelling across PET imaging systems: unified scatter modelling

RIN factor of all samples used for Illumina sequencing. (PDF 225 kb

Rapid processing of PET list-mode data for efficient uncertainty estimation and data analysis

In this technical note we propose a rapid and scalable software solution for the processing of PET list-mode data, which allows the efficient integration of list mode data processing into the workflow of image reconstruction and analysis. All processing is performed on the graphics processing unit (GPU), making use of streamed and concurrent kernel execution together with data transfers between disk and CPU memory as well as CPU and GPU memory. This approach leads to fast generation of multiple bootstrap realisations, and when combined with fast image reconstruction and analysis, it enables assessment of uncertainties of any image statistic and of any component of the image generation process (e.g., random correction, image processing) within reasonable time frames (e.g., within five minutes per realisation). This is of particular value when handling complex chains of image generation and processing. The software outputs the following: (1) estimate of expected random event data for noise reduction; (2) dynamic prompt and random sinograms of span-1 and span-11 and (3) variance estimates based on multiple bootstrap realisations of (1) and (2) assuming reasonable count levels for acceptable accuracy. In addition, the software produces statistics and visualisations for immediate quality control and crude motion detection, such as: (1) count rate curves; (2) centre of mass plots of the radiodistribution for motion detection; (3) video of dynamic projection views for fast visual list-mode skimming and inspection; (4) full normalisation factor sinograms. To demonstrate the software, we present an example of the above processing for fast uncertainty estimation of regional SUVR (standard uptake value ratio) calculation for a single PET scan of ¹⁸F-florbetapir using the Siemens Biograph mMR scanner

Short acquisition time PET quantification using MRI-based pharmacokinetic parameter synthesis

Positron Emission Tomography (PET) with pharmacokinetic (PK) modelling is a quantitative molecular imaging technique, however the long data acquisition time is prohibitive in clinical practice. An approach has been proposed to incorporate blood flow information from Arterial Spin Labelling (ASL) Magnetic Resonance Imaging (MRI) into PET PK modelling to reduce the acquisition time. This requires the conversion of cerebral blood flow (CBF) maps, measured by ASL, into the relative tracer delivery parameter (R 1 ) used in the PET PK model. This was performed regionally using linear regression between population R 1 and ASL values. In this paper we propose a novel technique to synthesise R 1 maps from ASL data using a database with both R 1 and CBF maps. The local similarity between the candidate ASL image and those in the database is used to weight the propagation of R 1 values to obtain the optimal patient specific R 1 map. Structural MRI data is also included to provide information within common regions of artefact in ASL data. This methodology is compared to the linear regression technique using leave one out analysis on 32 subjects. The proposed method significantly improves regional R 1 estimation (p < 0.001), reducing the error in the pharmacokinetic modelling. Furthermore, it allows this technique to be extended to a voxel level, increasing the clinical utility of the images

NiftyPET: A high-throughput software platform for high quantitative accuracy and precision PET imaging and analysis

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 coeffi- cient 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

Tau accumulation in autosomal dominant Alzheimer’s disease: a longitudinal [18F]flortaucipir study

Cortical tau accumulation is a key pathological event that partly defines Alzheimer’s disease (AD) onset and is associated with cognitive decline and future disease progression. However, an improved understanding of the timing and pattern of early tau deposition in AD and how this may be tracked in vivo is needed. Data from 59 participants involved in two longitudinal cohort studies of autosomal dominant AD (ADAD) were used to investigate whether tau PET can detect and track presymptomatic change; seven participants were symptomatic, and 52 were asymptomatic but at a 50% risk of carrying a pathogenic mutation. All had baseline flortaucipir (FTP) PET, MRI and clinical assessments; 26 individuals had more than one FTP PET scan. Standardised uptake value ratios (SUVRs) in prespecified regions of interest (ROIs) were obtained using inferior cerebellar grey matter as the reference region. We compared the changes in FTP SUVRs between presymptomatic carriers, symptomatic carriers and non-carriers, adjusting for age, sex and study site. We also investigated the relationship between regional FTP SUVRs and estimated years to/from symptom onset (EYO). Compared to both non-carriers and presymptomatic carriers, FTP SUVRs were significantly higher in symptomatic carriers in all ROIs tested (p 0.05), although increased FTP signal uptake was seen posteriorly in some individuals around the time of expected symptom onset. When we examined the relationship of FTP SUVR with respect to EYO, the earliest significant regional difference between mutation carriers and non-carriers was detected within the precuneus prior to estimated symptom onset in some cases. This study supports preliminary studies suggesting that presymptomatic tau tracer uptake is rare in ADAD. In cases where early uptake was seen, there was often a predilection for posterior regions (the precuneus and post-cingulate) as opposed to the medial temporal lobe, highlighting the importance of examining in vivo tau uptake beyond the confines of traditional Braak staging

Topological Surface States Protected From Backscattering by Chiral Spin Texture

Topological insulators are a new class of insulators in which a bulk gap for electronic excitations is generated by strong spin orbit coupling. These novel materials are distinguished from ordinary insulators by the presence of gapless metallic boundary states, akin to the chiral edge modes in quantum Hall systems, but with unconventional spin textures. Recently, experiments and theoretical efforts have provided strong evidence for both two- and three-dimensional topological insulators and their novel edge and surface states in semiconductor quantum well structures and several Bi-based compounds. A key characteristic of these spin-textured boundary states is their insensitivity to spin-independent scattering, which protects them from backscattering and localization. These chiral states are potentially useful for spin-based electronics, in which long spin coherence is critical, and also for quantum computing applications, where topological protection can enable fault-tolerant information processing. Here we use a scanning tunneling microscope (STM) to visualize the gapless surface states of the three-dimensional topological insulator BiSb and to examine their scattering behavior from disorder caused by random alloying in this compound. Combining STM and angle-resolved photoemission spectroscopy, we show that despite strong atomic scale disorder, backscattering between states of opposite momentum and opposite spin is absent. Our observation of spin-selective scattering demonstrates that the chiral nature of these states protects the spin of the carriers; they therefore have the potential to be used for coherent spin transport in spintronic devices.Comment: to be appear in Nature on August 9, 200

An Automated and High Precision Quantitative Analysis of the ACR Phantom

A novel phantom-imaging platform for automated and high precision imaging of the American College of Radiology (ACR) PET phantom is proposed. The platform facilitates the generation of an accurate μ-map for PET/MR systems with a robust alignment based on two-stage image registration using specifically designed PET templates. The automated analysis of PET images uses a set of granular composite volume of interest (VOI) templates in a 0.5 mm resolution grid for sampling of the system response to the insert step functions. The impact of the activity outside the field of view (FOV) was evaluated using two acquisitions of 30 minutes each, with and without the activity outside the FOV. Iterative image reconstruction was employed with and without modelled shift-invariant point spread function (PSF) and varying ordered subsets expectation maximisation (OSEM) iterations. Uncertainty analysis of all image-derived statistics was performed using bootstrap resampling of the list-mode data. We found that the activity outside the FOV can adversely affect the imaging planes close to the edge of the axial FOV, reducing the contrast, background uniformity and overall quantitative accuracy. The PSF had a positive impact on contrast recovery (although it slows convergence). The proposed platform may be helpful in a more informative evaluation of PET systems and image reconstruction methods

Multi-tracer model for staging cortical amyloid deposition using PET imaging

OBJECTIVE: To develop and evaluate a model for staging cortical amyloid deposition using PET with high generalizability. METHODS: 3027 subjects (1763 Cognitively Unimpaired (CU), 658 Impaired, 467 Alzheimer's disease (AD) dementia, 111 non-AD dementia, and 28 with missing diagnosis) from six cohorts (EMIF-AD, ALFA, ABIDE, ADC, OASIS-3, ADNI) who underwent amyloid PET were retrospectively included; 1049 subjects had follow-up scans. Applying dataset-specific cut-offs to global Standard Uptake Value ratio (SUVr) values from 27 regions, single-tracer and pooled multi-tracer regional rankings were constructed from the frequency of abnormality across 400 CU subjects (100 per tracer). The pooled multi-tracer ranking was used to create a staging model consisting of four clusters of regions as it displayed a high and consistent correlation with each single-tracer ranking. Relationships between amyloid stage, clinical variables and longitudinal cognitive decline were investigated. RESULTS: SUVr abnormality was most frequently observed in cingulate, followed by orbitofrontal, precuneal, and insular cortices, then the associative, temporal and occipital regions. Abnormal amyloid levels based on binary global SUVr classification were observed in 1.0%, 5.5%, 17.9%, 90.0%, and 100.0% of stage 0-4 subjects, respectively. Baseline stage predicted decline in MMSE (ADNI: N=867, F=67.37, p3000 subjects across cohorts and radiotracers, and detects pre-global amyloid burden and distinct risk profiles of cognitive decline within globally amyloid-positive subjects

The amyloid imaging for the prevention of Alzheimer's disease consortium: A European collaboration with global impact

Background: Amyloid-β (Aβ) accumulation is considered the earliest pathological change in Alzheimer's disease (AD). The Amyloid Imaging to Prevent Alzheimer's Disease (AMYPAD) consortium is a collaborative European framework across European Federation of Pharmaceutical Industries Associations (EFPIA), academic, and ‘Small and Medium-sized enterprises’ (SME) partners aiming to provide evidence on the clinical utility and cost-effectiveness of Positron Emission Tomography (PET) imaging in diagnostic work-up of AD and to support clinical trial design by developing optimal quantitative methodology in an early AD population. The AMYPAD studies: In the Diagnostic and Patient Management Study (DPMS), 844 participants from eight centres across three clinical subgroups (245 subjective cognitive decline, 342 mild cognitive impairment, and 258 dementia) were included. The Prognostic and Natural History Study (PNHS) recruited pre-dementia subjects across 11 European parent cohorts (PCs). Approximately 1600 unique subjects with historical and prospective data were collected within this study. PET acquisition with [18F]flutemetamol or [18F]florbetaben radiotracers was performed and quantified using the Centiloid (CL) method. Results: AMYPAD has significantly contributed to the AD field by furthering our understanding of amyloid deposition in the brain and the optimal methodology to measure this process. Main contributions so far include the validation of the dual-time window acquisition protocol to derive the fully quantitative non-displaceable binding potential (BPND), assess the value of this metric in the context of clinical trials, improve PET-sensitivity to emerging Aβ burden and utilize its available regional information, establish the quantitative accuracy of the Centiloid method across tracers and support implementation of quantitative amyloid-PET measures in the clinical routine. Future steps: The AMYPAD consortium has succeeded in recruiting and following a large number of prospective subjects and setting up a collaborative framework to integrate data across European PCs. Efforts are currently ongoing in collaboration with ARIDHIA and ADDI to harmonize, integrate, and curate all available clinical data from the PNHS PCs, which will become openly accessible to the wider scientific community