Resolution enhancement, noise suppression, and joint T2* decay estimation in dual-echo sodium-23 MR imaging using anatomically guided reconstruction

PURPOSE

Sodium MRI is challenging because of the low tissue concentration of the 23 Na nucleus and its extremely fast biexponential transverse relaxation rate. In this article, we present an iterative reconstruction framework using dual-echo 23 Na data and exploiting anatomical prior information (AGR) from high-resolution, low-noise, 1 H MR images. This framework enables the estimation and modeling of the spatially varying signal decay due to transverse relaxation during readout (AGRdm), which leads to images of better resolution and reduced noise resulting in improved quantification of the reconstructed 23 Na images.

METHODS

The proposed framework was evaluated using reconstructions of 30 noise realizations of realistic simulations of dual echo twisted projection imaging (TPI) 23 Na data. Moreover, three dual echo 23 Na TPI brain datasets of healthy controls acquired on a 3T Siemens Prisma system were reconstructed using conventional reconstruction, AGR and AGRdm.

RESULTS

Our simulations show that compared to conventional reconstructions, AGR and AGRdm show improved bias-noise characteristics in several regions of the brain. Moreover, AGR and AGRdm images show more anatomical detail and less noise in the reconstructions of the experimental data sets. Compared to AGR and the conventional reconstruction, AGRdm shows higher contrast in the sodium concentration ratio between gray and white matter and between gray matter and the brain stem.

CONCLUSION

AGR and AGRdm generate 23 Na images with high resolution, high levels of anatomical detail, and low levels of noise, potentially enabling high-quality 23 Na MR imaging at 3T.

The SNR of time-of-flight positron emission tomography data for joint reconstruction of the activity and attenuation images

Objective.Measurement of the time-of-flight (TOF) difference of each coincident pair of photons increases the effective sensitivity of positron emission tomography (PET). Many authors have analyzed the benefit of TOF for quantification and hot spot detection in the reconstructed activity images. However, TOF not only improves the effective sensitivity, it also enables the joint reconstruction of the tracer concentration and attenuation images. This can be used to correct for errors in CT- or MR-derived attenuation maps, or to apply attenuation correction without the help of a second modality. This paper presents an analysis of the effect of TOF on the variance of the jointly reconstructed attenuation and (attenuation corrected) tracer concentration images.Approach.The analysis is performed for PET systems that have a distribution of possibly non-Gaussian TOF-kernels, and includes the conventional Gaussian TOF-kernel as a special case. Non-Gaussian TOF-kernels are often observed in novel detector designs, which make use of two (or more) different mechanisms to convert the incoming 511 keV photon to optical photons. The analytical result is validated with a simple 2D simulation.Main results.We show that if two different TOF-kernels are equivalent for image reconstruction with known attenuation, then they are also equivalent for joint reconstruction of the activity and the attenuation images. The variance increase in the activity, caused by also jointly reconstructing the attenuation image, vanishes when the TOF-resolution approaches perfection.Significance.These results are of interest for PET detector development and for the development of stand-alone PET systems.

The SNR of Positron Emission Data With Gaussian and Non-Gaussian Time-of-Flight Kernels, With Application to Prompt Photon Coincidence

It is well known that measurement of the time-of-flight (TOF) increases the information provided by coincident events in positron emission tomography (PET). This information increase propagates through the reconstruction and improves the signal-to-noise ratio in the reconstructed images. Takehiro Tomitani has analytically computed the gain in variance in the reconstructed image, provided by a particular TOF resolution, for the center of a uniform disk and for a Gaussian TOF kernel. In this paper we extend this result, by computing the signal-to-noise ratio (SNR) contributed by individual coincidence events for two different tasks. One task is the detection of a hot spot in the center of a uniform cylinder. The second one is the same as that considered by Tomitani, i.e. the reconstruction of the central voxel in the image of a uniform cylinder. In addition, we extend the computation to non-Gaussian TOF kernels. It is found that a modification of the TOF-kernel changes the SNR for both tasks in almost exactly the same way. The proposed method can be used to compare TOF-systems with different and possibly event-dependent TOF-kernels, as encountered when prompt photons, such as Cherenkov photons are present, or when the detector is composed of different scintillators. The method is validated with simple 2D simulations and illustrated by applying it to PET detectors producing optical photons with event-dependent timing characteristics.

An optically-guided cochlear implant sheath for real-time monitoring of electrode insertion into the human cochlea

In cochlear implant surgery, insertion of perimodiolar electrode arrays into the scala tympani can be complicated by trauma or even accidental translocation of the electrode array within the cochlea. In patients with partial hearing loss, cochlear trauma can not only negatively affect implant performance, but also reduce residual hearing function. These events have been related to suboptimal positioning of the cochlear implant electrode array with respect to critical cochlear walls of the scala tympani (modiolar wall, osseous spiral lamina and basilar membrane). Currently, the position of the electrode array in relation to these walls cannot be assessed during the insertion and the surgeon depends on tactile feedback, which is unreliable and often comes too late. This study presents an image-guided cochlear implant device with an integrated, fiber-optic imaging probe that provides real-time feedback using optical coherence tomography during insertion into the human cochlea. This novel device enables the surgeon to accurately detect and identify the cochlear walls ahead and to adjust the insertion trajectory, avoiding collision and trauma. The functionality of this prototype has been demonstrated in a series of insertion experiments, conducted by experienced cochlear implant surgeons on fresh-frozen human cadaveric cochleae.

Low-dose CT allows for accurate proton therapy dose calculation and plan optimization

Objective. Protons offer a more conformal dose delivery compared to photons, yet they are sensitive to anatomical changes over the course of treatment. To minimize range uncertainties due to anatomical variations, a new CT acquisition at every treatment session would be paramount to enable daily dose calculation and subsequent plan adaptation. However, the series of CT scans results in an additional accumulated patient dose. Reducing CT radiation dose and thereby decreasing the potential risk of radiation exposure to patients is desirable, however, lowering the CT dose results in a lower signal-to-noise ratio and therefore in a reduced quality image. We hypothesized that the signal-to-noise ratio provided by conventional CT protocols is higher than needed for proton dose distribution estimation. In this study, we aim to investigate the effect of CT imaging dose reduction on proton therapy dose calculations and plan optimization. Approach. To verify our hypothesis, a CT dose reduction simulation tool has been developed and validated to simulate lower-dose CT scans from an existing standard-dose scan. The simulated lower-dose CTs were then used for proton dose calculation and plan optimization and the results were compared with those of the standard-dose scan. The same strategy was adopted to investigate the effect of CT dose reduction on water equivalent thickness (WET) calculation to quantify CT noise accumulation during integration along the beam. Main results. The similarity between the dose distributions acquired from the low-dose and standard-dose CTs was evaluated by the dose-volume histogram and the 3D Gamma analysis. The results on an anthropomorphic head phantom and three patient cases indicate that CT imaging dose reduction up to 90% does not have a significant effect on proton dose calculation and plan optimization. The relative error was employed to evaluate the similarity between WET maps and was found to be less than 1% after reducing the CT imaging dose by 90%. Significance. The results suggest the possibility of using low-dose CT for proton therapy dose estimation, since the dose distributions acquired from the standard-dose and low-dose CTs are clinically equivalent.

Experimental validation of a rodent PET scanner prototype based on a single LYSO crystal tube

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.

Can nuclear imaging accurately detect scar in ischemic cardiac resynchronization therapy candidates?

BACKGROUND

Accurate scar assessment is crucial in cardiac resynchronization therapy (CRT) candidates, since its presence is a negative predictor for CRT response. Therefore, we assessed the performance of different PET parameters to detect scar in CRT candidates.

METHODS

Twenty-nine CRT candidates underwent 18F-fluorodeoxyglucose (18F-FDG)-PET/computed tomography (CT), resting 13N-NH3-PET/CT and cardiac magnetic resonance (CMR) prior to CRT implantation. Segmental 18F-FDG uptake, late 13N-NH3 uptake and absolute myocardial blood flow (MBF) were evaluated for scar detection using late gadolinium enhancement (LGE) CMR as reference. A receiver operator characteristic (ROC) area under the curve (AUC) ≥0.8 indicated a good accuracy of the methods evaluated.

RESULTS

Scar was present in 111 of 464 segments. None of the approaches could reliably identify segments with nontransmural scar, except for 18F-FDG uptake in the lateral wall (AUC 0.83). Segmental transmural scars could be detected with all methods (AUC ≥ 0.8), except for septal 18F-FDG uptake and MBF in the inferior wall (AUC < 0.8). Late 13N-NH3 uptake was the best parameter for transmural scar detection, independent of its location, with a sensitivity of 80% and specificity of 92% using a cutoff of 66% of the maximum tracer activity.

CONCLUSIONS

Late 13N-NH3 uptake is superior to 13N-NH3 MBF and 18F-FDG in detecting transmural scar, independently of its location. However, none of the tested PET parameters was able to accurately detect nontransmural scar.

Design and modeling of a high resolution and high sensitivity PET brain scanner with double-ended readout

In the wake of recent advancements in scintillator, photodetector, and low-noise fast electronics technologies, as well as in fast reconstruction software, positron emission tomography (PET) scanners have seen considerable improvements in spatial resolution, time resolution, and absolute sensitivity. To continue this trend, we present a helmet type PET brain scanner design that combines high solid angle coverage and double-ended readout of 30 mm-thick scintillator crystals to achieve excellent absolute sensitivity, depth of interaction resolution, and time resolution. This scanner comprises 598 detector arrays, each with 8×8 Lu_1.8Y_0.2SiO₅:Ce (LYSO:Ce) crystals with dimensions 3.005×3.005×30 mm³one-to-one coupled on either end to silicon photomultipliers (SiPMs). Our Monte Carlo simulations based in the platform Geant4 predict that this scanner would attain an absolute sensitivity to a 35 cm line source placed at the center of the radial field of view of (17.1±0.1) %, a depth of interaction resolution of (3.99±0.05) mm, and a coincidence time resolution of (198±5) ps. Our simulations also predict radial, tangential, and axial spatial resolutions at the center of the field of view of 3.3 mm, 3.1 mm, and 3.3 mm, respectively. As this set of simultaneous parameters compares favorably to today's most advanced clinical PET scanners and other proposed designs, this scanner has a good chance of becoming a preferred tool for high quality brain imaging.

Segmentation-guided multi-modal registration of liver images for dose estimation in SIRT

Purpose

Selective internal radiation therapy (SIRT) requires a good liver registration of multi-modality images to obtain precise dose prediction and measurement. This study investigated the feasibility of liver registration of CT and MR images, guided by segmentation of the liver and its landmarks. The influence of the resulting lesion registration on dose estimation was evaluated.

Methods

The liver segmentation was done with a convolutional neural network (CNN), and the landmarks were segmented manually. Our image-based registration software and its liver-segmentation-guided extension (CNN-guided) were tuned and evaluated with 49 CT and 26 MR images from 20 SIRT patients. Each liver registration was evaluated by the root mean square distance (RMSD) of mean surface distance between manually delineated liver contours and mass center distance between manually delineated landmarks (lesions, clips, etc.). The root mean square of RMSDs (RRMSD) was used to evaluate all liver registrations. The CNN-guided registration was further extended by incorporating landmark segmentations (CNN&LM-guided) to assess the value of additional landmark guidance. To evaluate the influence of segmentation-guided registration on dose estimation, mean dose and volume percentages receiving at least 70 Gy (V70) estimated on the ^99mTc-labeled macro-aggregated albumin (^99mTc-MAA) SPECT were computed, either based on lesions from the reference ^99mTc-MAA CT (reference lesions) or from the registered floating CT or MR images (registered lesions) using the CNN- or CNN&LM-guided algorithms.

Results

The RRMSD decreased for the floating CTs and MRs by 1.0 mm (11%) and 3.4 mm (34%) using CNN guidance for the image-based registration and by 2.1 mm (26%) and 1.4 mm (21%) using landmark guidance for the CNN-guided registration. The quartiles for the relative mean dose difference (the V70 difference) between the reference and registered lesions and their correlations [25th, 75th; r] are as follows: [− 5.5% (− 1.3%), 5.6% (3.4%); 0.97 (0.95)] and [− 12.3% (− 2.1%), 14.8% (2.9%); 0.96 (0.97)] for the CNN&LM- and CNN-guided CT to CT registrations, [− 7.7% (− 6.6%), 7.0% (3.1%); 0.97 (0.90)] and [− 15.1% (− 11.3%), 2.4% (2.5%); 0.91 (0.78)] for the CNN&LM- and CNN-guided MR to CT registrations.

Conclusion

Guidance by CNN liver segmentations and landmarks markedly improves the performance of the image-based registration. The small mean dose change between the reference and registered lesions demonstrates the feasibility of applying the CNN&LM- or CNN-guided registration to volume-level dose prediction. The CNN&LM- and CNN-guided registrations for CTs can be applied to voxel-level dose prediction according to their small V70 change for most lesions. The CNN-guided MR to CT registration still needs to incorporate landmark guidance for smaller change of voxel-level dose estimation.

Rigid motion tracking using moments of inertia in TOF-PET brain studies

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 on¹⁸F-FDG as well as¹¹C-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 s¹⁸F-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.

Time of Flight in Perspective: Instrumental and Computational Aspects of Time Resolution in Positron Emission Tomography

The first time-of-flight positron emission tomography (TOF-PET) scanners were developed as early as in the 1980s. However, the poor light output and low detection efficiency of TOF-capable detectors available at the time limited any gain in image quality achieved with these TOF-PET scanners over the traditional non-TOF PET scanners. The discovery of LSO and other Lu-based scintillators revived interest in TOF-PET and led to the development of a second generation of scanners with high sensitivity and spatial resolution in the mid-2000s. The introduction of the silicon photomultiplier (SiPM) has recently yielded a third generation of TOF-PET systems with unprecedented imaging performance. Parallel to these instrumentation developments, much progress has been made in the development of image reconstruction algorithms that better utilize the additional information provided by TOF. Overall, the benefits range from a reduction in image variance (SNR increase), through allowing joint estimation of activity and attenuation, to better reconstructing data from limited angle systems. In this work, we review these developments, focusing on three broad areas: 1) timing theory and factors affecting the time resolution of a TOF-PET system; 2) utilization of TOF information for improved image reconstruction; and 3) quantification of the benefits of TOF compared to non-TOF PET. Finally, we offer a brief outlook on the TOF-PET developments anticipated in the short and longer term. Throughout this work, we aim to maintain a clinically driven perspective, treating TOF as one of multiple (and sometimes competitive) factors that can aid in the optimization of PET imaging performance.

2D feasibility study of joint reconstruction of attenuation and activity in limited angle TOF-PET

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.

Impact of left bundle branch block on myocardial perfusion and metabolism: A positron emission tomography study

BACKGROUND

Better understanding of pathophysiological changes, induced by left bundle branch block (LBBB), may improve patient selection for cardiac resynchronization therapy (CRT). Therefore, we assessed the effect of LBBB on regional glucose metabolism, 13N-NH3-derived absolute and semiquantitative myocardial blood flow (MBF), and their relation in non-ischemic CRT candidates.

METHODS

Twenty-five consecutive non-ischemic patients with LBBB underwent 18F-FDG and resting dynamic 13N-NH3 PET/CT prior to CRT implantation. Regional 18F-FDG uptake, absolute MBF, and late 13N-NH3 uptake were analyzed and corresponding septal-to-lateral wall ratios (SLR) were calculated. Segmental analysis was performed to evaluate "reverse mismatch," "mismatch," and "match" patterns, based on late 13N-NH3/18F-FDG uptake ratios.

RESULTS

A significantly lower 18F-FDG uptake was observed in the septum compared to the lateral wall (SLR 0.53 ± 0.17). A similar pattern was observed for MBF (SLR 0.68 ± 0.18), whereas late 13N-NH3 uptake showed a homogeneous distribution (SLR 0.96 ± 0.13). 13N-NH3/18F-FDG "mismatch" and "reverse mismatch" segments were predominantly present in the lateral (52%) and septal wall (61%), respectively.

CONCLUSIONS

Non-ischemic CRT candidates with LBBB demonstrate lower glucose uptake and absolute MBF in the septum compared to the lateral wall. However, late static 13N-NH3 uptake showed a homogenous distribution, reflecting a composite measure of altered regional MBF and metabolism, induced by LBBB.

A motion correction approach for oral and maxillofacial cone-beam CT imaging

Patient movement affects image quality in oral and maxillofacial cone-beam computed tomography imaging. While many efforts are made to minimize the possibility of motion during a scan, relatively little attention has been given to motion correction after acquisition. We propose a novel method which can improve the image quality after an oral and maxillofacial scan. The proposed method is based on our previous work and is a retrospective motion estimation and motion compensation (ME/MC) approach that iteratively estimates and compensates for rigid pose change over time. During motion estimation, image update and motion update are performed alternately in a multi-resolution scheme to obtain the motion. We propose use of a feature-based motion update and patch-based image update in the iterative estimation process, to alleviate the effect of limited scan field of view on estimation. During motion compensation, a fine-resolution image reconstruction was performed with compensation for the estimated motion. The proposed ME/MC method was evaluated with simulations, phantom and patient studies. Two experts in dentomaxillofacial radiology assessed the diagnostic importance of the resulting motion artifact suppression. The quality of the reconstructed images was improved after motion compensation, and most of the image artifacts were eliminated. Quantitative analysis by comparison to a reference image and by calculation of a sharpness metric agreed with the qualitative observation. The results are promising, and further evaluation is required to assess the clinical value of the proposed method.

Quantitative PET in the 2020s: a roadmap

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.

Approximating anatomically-guided PET reconstruction in image space using a convolutional neural network

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.

Characterization of a preclinical PET insert in a 7 tesla MRI scanner: beyond NEMA testing

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 mm³) 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.

Roadmap toward the 10 ps time-of-flight PET challenge

Since the seventies, positron emission tomography (PET) has become an invaluable medical molecular imaging modality with an unprecedented sensitivity at the picomolar level, especially for cancer diagnosis and the monitoring of its response to therapy. More recently, its combination with x-ray computed tomography (CT) or magnetic resonance (MR) has added high precision anatomic information in fused PET/CT and PET/MR images, thus compensating for the modest intrinsic spatial resolution of PET. Nevertheless, a number of medical challenges call for further improvements in PET sensitivity. These concern in particular new treatment opportunities in the context personalized (also called precision) medicine, such as the need to dynamically track a small number of cells in cancer immunotherapy or stem cells for tissue repair procedures. A better signal-to-noise ratio (SNR) in the image would allow detecting smaller size tumours together with a better staging of the patients, thus increasing the chances of putting cancer in complete remission. Moreover, there is an increasing demand for reducing the radioactive doses injected to the patients without impairing image quality. There are three ways to improve PET scanner sensitivity: improving detector efficiency, increasing geometrical acceptance of the imaging device and pushing the timing performance of the detectors. Currently, some pre-localization of the electron-positron annihilation along a line-of-response (LOR) given by the detection of a pair of annihilation photons is provided by the detection of the time difference between the two photons, also known as the time-of-flight (TOF) difference of the photons, whose accuracy is given by the coincidence time resolution (CTR). A CTR of about 10 picoseconds FWHM will ultimately allow to obtain a direct 3D volume representation of the activity distribution of a positron emitting radiopharmaceutical, at the millimetre level, thus introducing a quantum leap in PET imaging and quantification and fostering more frequent use of 11C radiopharmaceuticals. The present roadmap article toward the advent of 10 ps TOF-PET addresses the status and current/future challenges along the development of TOF-PET with the objective to reach this mythic 10 ps frontier that will open the door to real-time volume imaging virtually without tomographic inversion. The medical impact and prospects to achieve this technological revolution from the detection and image reconstruction point-of-views, together with a few perspectives beyond the TOF-PET application are discussed.

Quantitative comparison of pre-treatment predictive and post-treatment measured dosimetry for selective internal radiation therapy using cone-beam CT for tumor and liver perfusion territory definition

Background

Selective internal radiation therapy (SIRT) is a promising treatment for unresectable hepatic malignancies. Predictive dose calculation based on a simulation using ^99mTc-labeled macro-aggregated albumin (^99mTc-MAA) before the treatment is considered as a potential tool for patient-specific treatment planning. Post-treatment dose measurement is mainly performed to confirm the planned absorbed dose to the tumor and non-tumor liver volumes. This study compared the predicted and measured absorbed dose distributions.

Methods

Thirty-one patients (67 tumors) treated by SIRT with resin microspheres were analyzed. Predicted and delivered absorbed dose was calculated using ^99mTc-MAA-SPECT and ⁹⁰Y-TOF-PET imaging. The voxel-level dose distribution was derived using the local deposition model. Liver perfusion territories and tumors have been delineated on contrast-enhanced CBCT images, which have been acquired during the ^99mTc-MAA work-up. Several dose-volume histogram (DVH) parameters together with the mean dose for liver perfusion territories and non-tumoral and tumoral compartments were evaluated.

Results

A strong correlation between the predicted and measured mean dose for non-tumoral volume was observed (r = 0.937). The ratio of measured and predicted mean dose to this volume has a first, second, and third interquartile range of 0.83, 1.05, and 1.25. The difference between the measured and predicted mean dose did not exceed 11 Gy. The correlation between predicted and measured mean dose to the tumor was moderate (r = 0.623) with a mean difference of − 9.3 Gy. The ratio of measured and predicted tumor mean dose had a median of 1.01 with the first and third interquartile ranges of 0.58 and 1.59, respectively. Our results suggest that ^99mTc-MAA-based dosimetry could predict under or over dosing of the non-tumoral liver parenchyma for almost all cases. For more than two thirds of the tumors, a predictive absorbed dose correctly indicated either good tumor dose coverage or under-dosing of the tumor.

Conclusion

Our results highlight the predictive value of ^99mTc-MAA-based dose estimation to predict non-tumor liver irradiation, which can be applied to prescribe an optimized activity aiming at avoiding liver toxicity. Compared to non-tumoral tissue, a poorer agreement between predicted and measured absorbed dose is observed for tumors.

Estimation of Crystal Timing Properties and Efficiencies for the Improvement of (Joint) Maximum-Likelihood Reconstructions in TOF-PET

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.

Analysis of partial volume correction on quantification and regional heterogeneity in cardiac PET

BACKGROUND

The partial volume correction (PVC) of cardiac PET datasets using anatomical side information during reconstruction is appealing but not straightforward. Other techniques, which do not make use of additional anatomical information, could be equally effective in improving the reconstructed myocardial activity.

METHODS

Resolution modeling in combination with different noise suppressing priors was evaluated as a means to perform PVC. Anatomical priors based on a high-resolution CT are compared to non-anatomical, edge-preserving priors (relative difference and total variation prior). The study is conducted on ex vivo datasets from ovine hearts. A simulation study additionally clarifies the relationship between prior effectiveness and myocardial wall thickness.

RESULTS

Simple resolution modeling during data reconstruction resulted in over- and underestimation of activity, which hampers the absolute left ventricular quantification when compared to the ground truth. Both the edge-preserving and the anatomy-based PVC techniques improve the absolute quantification, with comparable results (Student t-test, P = .17). The relative tracer distribution was preserved with any reconstruction technique (repeated ANOVA, P = .98).

CONCLUSIONS

The use of edge-preserving priors emerged as optimal choice for quantification of tracer uptake in the left ventricular wall of the available datasets. Anatomical priors visually outperformed edge-preserving priors when the thinnest structures were of interest.

Benefits of Using a Spatially-Variant Penalty Strength With Anatomical Priors in PET Reconstruction

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.

Low septal to lateral wall 18F-FDG ratio is highly associated with mechanical dyssynchrony in non-ischemic CRT candidates

Background

In order to better understand the concept of mechanical dyssynchrony, a promising hallmark of cardiac resynchronization therapy (CRT) response, we investigated its effect on regional myocardial metabolism and myocardial blood flow (MBF) in non-ischemic CRT candidates.

Results

Thirty consecutive non-ischemic CRT eligible patients underwent static ¹⁸F-FDG and resting dynamic ¹³N-NH₃ PET/CT. ¹⁸F-FDG uptake and MBF for septal and lateral wall were analysed and septal-to-lateral wall ratios (SLR) were calculated. Based on the presence of mechanical dyssynchrony (septal flash and/or apical rocking) on echocardiography, patients were divided into 2 groups, with (n = 23) and without (n = 7) mechanical dyssynchrony.

Patients with mechanical dyssynchrony had significantly lower ¹⁸F-FDG SUVmean in the septum compared with the lateral wall (5.58 ± 2.65 vs 11.19 ± 4.10, p < 0.0001), while patients without mechanical dyssynchrony had a more homogeneous ¹⁸F-FDG distribution (7.33 ± 2.88 vs 8.31 ± 2.50, respectively, p = 0.30). Similarly, MBF was significantly different between the septal and lateral wall in the dyssynchrony group (0.57 ± 0.11 ml/g/min vs 0.92 ± 0.23 ml/g/min, respectively, p < 0.0001), whereas no difference was observed in the non-dyssynchrony group (0.61 ± 0.23 ml/g/min vs 0.77 ± 0.21 ml/g/min, respectively, p = 0.16). ¹⁸F-FDG SLR, but not MBF SLR, was associated with the presence of mechanical dyssynchrony and showed a significant inverse correlation with volumetric reverse remodeling after CRT (r = − 0.62, p = 0.001).

Conclusions

Non-ischemic heart failure patients with mechanical dyssynchrony demonstrate heterogeneous regional metabolism and MBF compared with patients without dyssynchrony. However, only ¹⁸F-FDG SLR appeared to be highly associated with the presence of mechanical dyssynchrony.

Trial registration

Clinicaltrials, NCT02537782. Registered 2 September 2015.

Partial volume and motion correction in cardiac PET: First results from an in vs ex vivo comparison using animal datasets

BACKGROUND

In a previous study on ex vivo, static cardiac datasets, we investigated the benefits of performing partial volume correction (PVC) in cardiac 18F-Fluorodeoxyglucose(FDG) PET datasets. In the present study, we extend the analysis to in vivo cardiac datasets, with the aim of defining which reconstruction technique maximizes quantitative accuracy and, ultimately, makes PET a better diagnostic tool for cardiac pathologies.

METHODS

In vivo sheep datasets were acquired and reconstructed with/without motion correction and using several reconstruction algorithms (with/without resolution modeling, with/without non-anatomical priors). Corresponding ex vivo scans of the excised sheep hearts were performed on a small-animal PET scanner (Siemens Focus 220, microPET) to provide high-resolution reference data unaffected by respiratory and cardiac motion. A comparison between the in vivo cardiac reconstructions and the corresponding ex vivo ground truth was performed.

RESULTS

The use of an edge-preserving prior (Total Variation (TV) prior in this work) in combination with motion correction reduces the bias in absolute quantification when compared to the standard clinical reconstructions (- 0.83 vs - 3.74 SUV units), when the end-systolic gate is considered. At end-diastole, motion correction improves absolute quantification but the PVC with priors does not improve the similarity to the ground truth more than a regular iterative reconstruction with motion correction and without priors. Relative quantification was not influenced much by the chosen reconstruction algorithm.

CONCLUSIONS

The relative ranking of the algorithms suggests superiority of the PVC reconstructions with dual gating in terms of overall absolute quantification and noise properties. A well-tuned edge-preserving prior, such as TV, enhances the noise properties of the resulting images of the heart. The end-systolic gate yields the most accurate quantification of cardiac datasets.

Quantification and discriminative power of 18F-FE-PE2I PET in patients with Parkinson's disease

RATIONALE

Dopamine transporter (DAT) imaging is an important adjunct in the diagnostic workup of patients with Parkinsonism. ¹⁸F-FE-PE2I is a suitable PET radioligand for DAT quantification and imaging with good pharmacokinetics. The aim of this study was to determine a clinical optimal simplified reference tissue-based image acquisition protocol and to compare the discriminatory value and effect size for ¹⁸F-FE-PE2I to that for ¹²³I-FP-CIT scan currently used in clinical practice.

METHODS

Nine patients with early Parkinson's disease (PD, 64.3 ± 6.8 years, 3M), who had previously undergone a ¹²³I-FP-CIT scan as part of their diagnostic workup, and 34 healthy volunteers (HV, 47.7 ± 16.8 years, 13M) underwent a 60-min dynamic ¹⁸F-FE-PE2I PET-MR scan on a GE Signa 3T PET-MR. Based on dynamic data and MR-based VOI delineation, BP_ND, semi-quantitative uptake ratio and SUVR_[t1-t2] images were calculated using either occipital cortex or cerebellum as reference region. For start-and-end time of the SUVR interval, three time frames [t₁-t₂] were investigated: [15-40] min, [40-60] min, and [50-60] min postinjection. Data for putamen (PUT) and caudate nucleus-putamen ratio (CPR) were compared in terms of quantification bias versus BP_ND and discriminative power.

RESULTS

Using occipital cortex as reference region resulted in smaller bias of SUVR with respect to BP_ND + 1 and higher correlation between SUVR and BP_ND + 1 compared with using cerebellum, irrespective of SUVR [t₁-t₂] interval. Smallest bias was observed with the [15-40]-min time window, in accordance with previous literature. The correlation between BP_ND + 1 and SUVR was slightly better for the late time windows. Discriminant analysis between PD and HV using both PUT and CPR SUVRs showed an accuracy of ≥ 90%, for both reference regions and all studied time windows. Semi-quantitative ¹²³I-FP-CIT and ¹⁸F-FE-PE2I values and relative decrease in the striatum for patients were highly correlated, with a higher effect size for ¹⁸F-FE-PE2I for PUT and CPR SUVR.

CONCLUSION

¹⁸F-FE-PE2I is a suitable radioligand for in vivo DAT imaging with high discriminative power between early PD and healthy controls. Whereas a [15-40]-min window has lowest bias with respect to BP_ND, a [50-60]-min time window at pseudoequilibrium can be advocated in terms of clinical feasibility with optimal discriminative power. The occipital cortex may be slightly preferable as reference region because of the higher time stability, stronger correlation of SUVR with BP_ND + 1, and lower bias. Moreover, the data suggest that the diagnostic accuracy of a 10-min static ¹⁸F-FE-PE2I scan is non-inferior compared with ¹²³I-FP-CIT scan used in standard clinical practice.

P619Stress-strain loop area better represents regional myocardial work than pressure-strain loop area in the dyssynchronous and remodelled left ventricle

Background

Pressure-strain-loops (PSL) have been suggested as surrogate parameter of regional myocardial work. However, in left ventricles (LV) with inhomogeneous remodelling, e.g. due to left bundle branch block (LBBB), wall stress may be unevenly distributed. Stress-strain loops (SSL) include information on both regional wall thickness and curvature, and may therefore provide a better surrogate.

Study plan

We therefore compared the correlation of segmental myocardial work estimated through both PSL and SSL to segmental myocardial glucose metabolism as a gold standard, in an animal model of pacing-induced LV remodelling.

Methods

Twelve sheep developed LV dilatation, thinned septum and thickened lateral wall, due to eight weeks of rapid right-atrial and right-ventricular free wall (DDD) pacing (at 180 bpm), causing a LBBB-like dyssynchrony. Invasive LV pressure and echocardiographic speckle tracking based circumferential strain were used to construct PSL. SSL were calculated by considering in addition dynamic changes in segmental myocardial wall thickness and curvature using the formula of Laplace. 18F-fluorodeoxyglucose (FDG)-uptake was measured by positron emission tomography (PET) in absolute values as standardised uptake ratio (SUR). Spatial resolution of PET was improved by ECG- and breathing-gating and using anatomical priors. All imaging was performed during dyssynchronous DDD-pacing and synchronous AAI-pacing (right-atrial pacing only), at baseline (n=3 animals), and after eight weeks of pacing induced remodelling (n=12 animals).

Results

Both at baseline (Fig. A+B) and after 8 weeks (Fig. D+E), switching between AAI and DDD-pacing caused an acute re-distribution of regional myocardial work as measured by both PSL and SSL. In contrast to PSL, however, SSL identified more regional differences among walls in remodelled hearts and showed clearer regional changes when switching between AAI and DDD-pacing. The correlation between regional work, assessed by PSL and SSL, and metabolism by PET, was comparable at baseline (r=0.65 and r=0.64, respectively) (Fig. C). In remodelled hearts after 8 weeks, however, the correlation of regional work assessed by SSL and glucose uptake by PET was significantly higher compared to PSL (r=0.73 vs. r=0.59, respectively; p<0.05) (Fig. F).

Conclusions

Regional myocardial work assessed by stress-strain loops correlates significantly better to regional metabolism as measured by PET glucose uptake, particularly after remodelling. Our findings therefore suggest that integrating information on wall thickness and curvature is essential for the reliable assessment of regional myocardial work, especially in dyssynchronous and remodelled left ventricles.

A Quantitative Evaluation of Joint Activity and Attenuation Reconstruction in TOF PET/MR Brain Imaging

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.

Multi-modal image analysis for semi-automatic segmentation of the total liver and liver arterial perfusion territories for radioembolization

Purpose

We have developed a multi-modal imaging approach for SIRT, combining ^99mTc-MAA SPECT/CT and/or ⁹⁰Y PET, ¹⁸F-FDG PET/CT, and contrast-enhanced CBCT for voxel-based dosimetry, as a tool for treatment planning and verification. For radiation dose prediction calculations, a segmentation of the total liver volume and of the liver perfusion territories is required.

Method

In this paper, we proposed a procedure for multi-modal image analysis to assist SIRT treatment planning. The pre-treatment ¹⁸F-FDG PET/CT, ^99mTc-MAA SPECT/CT, and contrast-enhanced CBCT images were registered to a common space using an initial rigid, followed by a deformable registration. The registration was scored by an expert using Likert scores. The total liver was segmented semi-automatically based on the PET/CT and SPECT/CT images, and the liver perfusion territories were determined based on the CBCT images. The segmentations of the liver and liver lobes were compared to the manual segmentations by an expert on a CT image.

Result

Our methodology showed that multi-modal image analysis can be used for determination of the liver and perfusion territories using CBCT in SIRT using all pre-treatment studies. The results for image registration showed acceptable alignment with limited impact on dosimetry.

The image registration performs well according to the expert reviewer (scored as perfect or with little misalignment in 94% of the cases). The semi-automatic liver segmentation agreed well with manual liver segmentation (dice coefficient of 0.92 and an average Hausdorff distance of 3.04 mm). The results showed disagreement between lobe segmentation using CBCT images compared to lobe segmentation based on CT images (average Hausdorff distance of 14.18 mm), with a high impact on the dosimetry (differences up to 9 Gy for right and 21 Gy for the left liver lobe).

Conclusion

This methodology can be used for pre-treatment dosimetry and for SIRT planning including the determination of the activity that should be administered to achieve the therapeutical goal. The inclusion of perfusion CBCT enables perfusion-based definition of the liver lobes, which was shown to be markedly different from the anatomical definition in some of the patients.

Feasibility Study of a Small Animal PET Insert Based on a Single LYSO Monolithic Tube

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.

Joint Reconstruction of Activity and Attenuation in Time-of-Flight PET: A Quantitative Analysis

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.

Time-of-flight PET time calibration using data consistency

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.

Combined PET/MRI: Global Warming-Summary Report of the 6th International Workshop on PET/MRI, March 27-29, 2017, Tübingen, Germany

The 6th annual meeting to address key issues in positron emission tomography (PET)/magnetic resonance imaging (MRI) was held again in Tübingen, Germany, from March 27 to 29, 2017. Over three days of invited plenary lectures, round table discussions and dialogue board deliberations, participants critically assessed the current state of PET/MRI, both clinically and as a research tool, and attempted to chart future directions. The meeting addressed the use of PET/MRI and workflows in oncology, neurosciences, infection, inflammation and chronic pain syndromes, as well as deeper discussions about how best to characterise the tumour microenvironment, optimise the complementary information available from PET and MRI, and how advanced data mining and bioinformatics, as well as information from liquid biomarkers (circulating tumour cells and nucleic acids) and pathology, can be integrated to give a more complete characterisation of disease phenotype. Some issues that have dominated previous meetings, such as the accuracy of MR-based attenuation correction (AC) of the PET scan, were finally put to rest as having been adequately addressed for the majority of clinical situations. Likewise, the ability to standardise PET systems for use in multicentre trials was confirmed, thus removing a perceived barrier to larger clinical imaging trials. The meeting openly questioned whether PET/MRI should, in all cases, be used as a whole-body imaging modality or whether in many circumstances it would best be employed to give an in-depth study of previously identified disease in a single organ or region. The meeting concluded that there is still much work to be done in the integration of data from different fields and in developing a common language for all stakeholders involved. In addition, the participants advocated joint training and education for individuals who engage in routine PET/MRI. It was agreed that PET/MRI can enhance our understanding of normal and disrupted biology, and we are in a position to describe the in vivo nature of disease processes, metabolism, evolution of cancer and the monitoring of response to pharmacological interventions and therapies. As such, PET/MRI is a key to advancing medicine and patient care.

Evaluation of Parallel Level Sets and Bowsher's Method as Segmentation-Free Anatomical Priors for Time-of-Flight PET Reconstruction

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.

Plane-dependent ML scatter scaling: 3D extension of the 2D simulated single scatter (SSS) estimate

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 ⁶⁸Ga-PSMA scan, and 23 whole-body ¹⁸F-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.

Optimized MLAA for quantitative non-TOF PET/MR of the brain

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 [¹⁸F]FDG patients, 35 [¹¹C]PiB patients and 1 [¹⁸F]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.

The effect of isoflurane on 18F-FDG uptake in the rat brain: a fully conscious dynamic PET study using motion compensation

BACKGROUND

In preclinical positron emission tomography (PET) studies an anaesthetic is used to ensure that the animal does not move during the scan. However, anaesthesia may have confounding effects on the drug or tracer kinetics under study, and the nature of these effects is usually not known.

METHOD

We have implemented a protocol for tracking the rigid motion of the head of a fully conscious rat during a PET scan and performing a motion compensated list-mode reconstruction of the data. Using this technique we have conducted eight rat studies to investigate the effect of isoflurane on the uptake of ¹⁸F-FDG in the brain, by comparing conscious and unconscious scans.

RESULTS

Our results indicate that isoflurane significantly decreases the whole brain uptake, as well as decreasing the relative regional FDG uptake in the cortex, diencephalon, and inferior colliculi, while increasing it in the vestibular nuclei. No statistically significant changes in FDG uptake were observed in the cerebellum and striata.

CONCLUSION

The applied event-based motion compensation technique allowed for the investigation of the effect of isoflurane on FDG uptake in the rat brain using fully awake and unrestrained rats, scanned dynamically from the moment of injection. A significant effect of the anaesthesia was observed in various regions of the brain.

Combined PET/MRI: from Status Quo to Status Go. Summary Report of the Fifth International Workshop on PET/MR Imaging; February 15-19, 2016; Tübingen, Germany

This article provides a collaborative perspective of the discussions and conclusions from the fifth international workshop of combined positron emission tomorgraphy (PET)/magnetic resonance imaging (MRI) that was held in Tübingen, Germany, from February 15 to 19, 2016. Specifically, we summarise the second part of the workshop made up of invited presentations from active researchers in the field of PET/MRI and associated fields augmented by round table discussions and dialogue boards with specific topics. This year, this included practical advice as to possible approaches to moving PET/MRI into clinical routine, the use of PET/MRI in brain receptor imaging, in assessing cardiovascular diseases, cancer, infection, and inflammatory diseases. To address perceived challenges still remaining to innovatively integrate PET and MRI system technologies, a dedicated round table session brought together key representatives from industry and academia who were engaged with either the conceptualisation or early adoption of hybrid PET/MRI systems. Discussions during the workshop highlighted that emerging unique applications of PET/MRI such as the ability to provide multi-parametric quantitative and visual information which will enable not only overall disease detection but also disease characterisation would eventually be regarded as compelling arguments for the adoption of PET/MR. However, as indicated by previous workshops, evidence in favour of this observation is only growing slowly, mainly due to the ongoing inability to pool data cohorts from independent trials as well as different systems and sites. The participants emphasised that moving from status quo to status go entails the need to adopt standardised imaging procedures and the readiness to act together prospectively across multiple PET/MRI sites and vendors.