Dispersion-corrected extracorporeal arterial input functions in PET studies of mice: a comparison to intracorporeal microprobe measurements.

BACKGROUND: Kinetic modelling of dynamic PET typically requires knowledge of the arterial radiotracer concentration (arterial input function, AIF). Its accurate determination is very difficult in mice. AIF measurements in an extracorporeal shunt can be performed; however, this introduces catheter dispersion. We propose a framework for extracorporeal dispersion correction and validated it by comparison to invasively determined intracorporeal AIFs using implanted microprobes. RESULTS: The response of an extracorporeal radiation detector to radioactivity boxcar functions, characterised by a convolution-based dispersion model, gave best fits using double-gamma variate and single-gamma variate kernels compared to mono-exponential kernels for the investigated range of flow rates. Parametric deconvolution with the optimal kernels was performed on 9 mice that were injected with a bolus of 39 ± 25 MBq [18F]F-PSMA-1007 after application of an extracorporeal circulation for three different flow rates in order to correct for dispersion. Comparison with synchronous implantation of microprobes for invasive aortic AIF recordings showed favourable correspondence, with no significant difference in terms of area-under-curve after 300 s and 5000 s. One-tissue and two-tissue compartment model simulations were performed to investigate differences in kinetic parameters between intra- and extracorporeally measured AIFs. Results of the modelling study revealed kinetic parameters close to the chosen simulated values in all compartment models. CONCLUSION: The high correspondence of simultaneously intra- and extracorporeally determined AIFs and resulting model parameters establishes a feasible framework for extracorporeal dispersion correction. This should allow more precise and accurate kinetic modelling in small animal experiments.

Voxel-Based Analysis of the Relation of 3'-Deoxy-3'-[18F]fluorothymidine ([18F]FLT) PET and Diffusion-Weighted (DW) MR Signals in Subcutaneous Tumor Xenografts Does Not Reveal a Direct Spatial Relation of These Two Parameters.

PURPOSE: Multimodal molecular imaging allows a direct coregistration of different images, facilitating analysis of the spatial relation of various imaging parameters. Here, we further explored the relation of proliferation, as measured by [18F]FLT PET, and water diffusion, as an indicator of cellular density and cell death, as measured by diffusion-weighted (DW) MRI, in preclinical tumor models. We expected these parameters to be negatively related, as highly proliferative tissue should have a higher density of cells, hampering free water diffusion. PROCEDURES: Nude mice subcutaneously inoculated with either lung cancer cells (n = 11 A549 tumors, n = 20 H1975 tumors) or colorectal cancer cells (n = 13 Colo205 tumors) were imaged with [18F]FLT PET and DW-MRI using a multimodal bed, which was transferred from one instrument to the other within the same imaging session. Fiducial markers allowed coregistration of the images. An automatic post-processing was developed in MATLAB handling the spatial registration of DW-MRI (measured as apparent diffusion coefficient, ADC) and [18F]FLT image data and subsequent voxel-wise analysis of regions of interest (ROIs) in the tumor. RESULTS: Analyses were conducted on a total of 76 datasets, comprising a median of 2890 data points (ranging from 81 to 13,597). Scatterplots showing [18F]FLT vs. ADC values displayed various grades of relations (Pearson correlation coefficient (PCC) varied from - 0.58 to 0.49, median: -0.07). When relating PCC to tumor volume (median: 46 mm3, range: 3 mm3 to 584 mm3), lung tumors tended to have a more pronounced negative spatial relation of [18F]FLT and ADC with increasing tumor size. However, due to the low number of large tumors (> ~ 200 mm3), this conclusion has to be treated with caution. CONCLUSIONS: A spatial relation of water diffusion, as measured by DW-MRI, and cellular proliferation, as measured by [18F]FLT PET, cannot be detected in the experimental datasets investigated in this study.

Respiratory Motion Detection and Correction for MR Using the Pilot Tone: Applications for MR and Simultaneous PET/MR Examinations.

OBJECTIVES: The aim of this study was to develop a method for tracking respiratory motion throughout full MR or PET/MR studies that requires only minimal additional hardware and no modifications to the sequences. MATERIALS AND METHODS: Patient motion that is caused by respiration affects the quality of the signal of the individual radiofrequency receive coil elements. This effect can be detected as a modulation of a monofrequent signal that is emitted by a small portable transmitter placed inside the bore (Pilot Tone). The frequency is selected such that it is located outside of the frequency band of the actual MR readout experiment but well within the bandwidth of the radiofrequency receiver, that is, the oversampling area. Temporal variations of the detected signal indicate motion. After extraction of the signal from the raw data, principal component analysis was used to identify respiratory motion. The approach and potential applications during MR and PET/MR examinations that rely on a continuous respiratory signal were validated with an anthropomorphic, PET/MR-compatible motion phantom as well as in a volunteer study. RESULTS: Respiratory motion detection and correction were presented for MR and PET data in phantom and volunteer studies. The Pilot Tone successfully recovered the ground-truth respiratory signal provided by the phantom. CONCLUSIONS: The presented method provides reliable respiratory motion tracking during arbitrary imaging sequences throughout a full PET/MR study. All results can directly be transferred to MR-only applications as well.

Dynamic Cell Imaging in PET With Optimal Transport Regularization.

We propose a novel dynamic image reconstruction method from PET listmode data that could be particularly suited to tracking single or small numbers of cells. In contrast to conventional PET reconstruction our method combines the information from all detected events not only to reconstruct the dynamic evolution of the radionuclide distribution, but also to improve the reconstruction at each single time point by enforcing temporal consistency. This is achieved via optimal transport regularization where in principle, among all possible temporally evolving radionuclide distributions consistent with the PET measurement, the one is chosen with least kinetic motion energy. The reconstruction is found by convex optimization so that there is no dependence on the initialization of the method. We study its behaviour on simulated data of a human PET system and demonstrate its robustness even in settings with very low radioactivity. In contrast to previously reported cell tracking algorithms, our technique is oblivious to the number of tracked cells. Without any additional complexity one or multiple cells can be reconstructed, and the model automatically determines the number of particles. For instance, four radiolabelled cells moving at a velocity of 3.1 mm/s and a PET recorded count rate of 1.1 cps (for each cell) could be simultaneously tracked with a tracking accuracy of 5.3 mm inside a simulated human body.

Thymidine Metabolism as a Confounding Factor for 3'-Deoxy-3'-18F-Fluorothymidine Uptake After Therapy in a Colorectal Cancer Model.

Noninvasive monitoring of tumor therapy response helps in developing personalized treatment strategies. Here, we performed sequential PET and diffusion-weighted MRI to evaluate changes induced by a FOLFOX-like combination chemotherapy in colorectal cancer xenografts, to identify the cellular and molecular determinants of these imaging biomarkers. Methods: Tumor-bearing CD1 nude mice, engrafted with FOLFOX-sensitive Colo205 colorectal cancer xenografts, were treated with FOLFOX (5-fluorouracil, leucovorin, and oxaliplatin) weekly. On days 1, 2, 6, 9, and 13 of therapy, tumors were assessed by in vivo imaging and ex vivo analyses. In addition, HCT116 xenografts, which did not respond to the FOLFOX treatment, were imaged on day 1 of therapy. Results: In Colo205 xenografts, FOLFOX induced a profound increase in uptake of the proliferation PET tracer 3'-deoxy-3'-18F-fluorothymidine (18F-FLT) accompanied by increases in markers for proliferation (Ki-67, thymidine kinase 1) and for activated DNA damage response (γH2AX), whereas the effect on cell death was minimal. Because tracer uptake was unaltered in the HCT116 model, these changes appear to be specific for tumor response. Conclusion: We demonstrated that 18F-FLT PET can noninvasively monitor cancer treatment-induced molecular alterations, including thymidine metabolism and DNA damage response. The cellular or imaging changes may not, however, be directly related to therapy response as assessed by volumetric measurements.

PET attenuation correction for flexible MRI surface coils in hybrid PET/MRI using a 3D depth camera.

PET attenuation correction for flexible MRI radio frequency surface coils in hybrid PET/MRI is still a challenging task, as position and shape of these coils conform to large inter-patient variabilities. The purpose of this feasibility study is to develop a novel method for the incorporation of attenuation information about flexible surface coils in PET reconstruction using the Microsoft Kinect V2 depth camera. The depth information is used to determine a dense point cloud of the coil's surface representing the shape of the coil. From a CT template-acquired once in advance-surface information of the coil is extracted likewise and converted into a point cloud. The two point clouds are then registered using a combination of an iterative-closest-point (ICP) method and a partially rigid registration step. Using the transformation derived through the point clouds, the CT template is warped and thereby adapted to the PET/MRI scan setup. The transformed CT template is converted into an attenuation map from Hounsfield units into linear attenuation coefficients. The resulting fitted attenuation map is then integrated into the MRI-based patient-specific DIXON-based attenuation map of the actual PET/MRI scan. A reconstruction of phantom PET data acquired with the coil present in the field-of-view (FoV), but without the corresponding coil attenuation map, shows large artifacts in regions close to the coil. The overall count loss is determined to be around 13% compared to a PET scan without the coil present in the FoV. A reconstruction using the new µ-map resulted in strongly reduced artifacts as well as increased overall PET intensities with a remaining relative difference of about 1% to a PET scan without the coil in the FoV.

Anthropomorphic thorax phantom for cardio-respiratory motion simulation in tomographic imaging.

Patient motion during medical imaging using techniques such as computed tomography (CT), magnetic resonance imaging (MRI), positron emission tomography (PET), or single emission computed tomography (SPECT) is well known to degrade images, leading to blurring effects or severe artifacts. Motion correction methods try to overcome these degrading effects. However, they need to be validated under realistic conditions. In this work, a sophisticated anthropomorphic thorax phantom is presented that combines several aspects of a simulator for cardio-respiratory motion. The phantom allows us to simulate various types of cardio-respiratory motions inside a human-like thorax, including features such as inflatable lungs, beating left ventricular myocardium, respiration-induced motion of the left ventricle, moving lung lesions, and moving coronary artery plaques. The phantom is constructed to be MR-compatible. This means that we can not only perform studies in PET, SPECT and CT, but also inside an MRI system. The technical features of the anthropomorphic thorax phantom Wilhelm are presented with regard to simulating motion effects in hybrid emission tomography and radiotherapy. This is supplemented by a study on the detectability of small coronary plaque lesions in PET/CT under the influence of cardio-respiratory motion, and a study on the accuracy of left ventricular blood volumes.

Data-Driven Methods for the Determination of Anterior-Posterior Motion in PET.

Physiological motion combined with elongated scanning times in PET leads to image degradation and quantification errors. Correction approaches usually require 1-D signals that can be obtained with hardware-based or data-driven methods. Most of the latter are optimized or limited to capture internal motion along the superior-inferior (S-I) direction. In this work we present methods for also extracting anterior-posterior (A-P) motion from PET data and propose a set of novel weighting mechanisms that can be used to emphasize certain lines-of-response (LORs) for an increased sensitivity and better signal-to-noise ratio (SNR). The proper functioning of the methods was verified in a phantom experiment. Further, their application to clinical [18F]-FDG-PET data of 72 patients revealed that using the weighting mechanisms leads to signals with significantly higher spectral respiratory weights, i.e. signals with higher quality. Information about multi-dimensional motion is contained in PET data and can be derived with data-driven methods. Motion models or correction techniques such as respiratory gating might benefit from the proposed methods as they allow to describe the three-dimensional movements of PET-positive structures more precisely.

Sex differences in absolute myocardial perfusion. Non-invasive H2(15)O-PET in young healthy adults.

AIM: To investigate sex differences in myocardial perfusion especially in healthy individuals since former studies are rare and findings are controversial. Participants, methods: 26 subjects were enrolled: 16 healthy women (age: 34 ±7 years) were compared with 10 healthy men (age: 34 ± 3 years; p = ns). Myocardial blood flow (MBF) and coronary vascular resistance (CVR) were quantified at rest, during adenosine infusion and cold-pressor-testing, using positron emission tomography and radioactive-labelled water (H2(15)O-PET). RESULTS: Women showed higher MBF than men at rest (1.10 ± 0.18 vs. 0.85 ± 0.20 ml/min/ml; p = 0.003) and cold-stress (1.39 ± 0.38 vs. 1.06 ± 0.28 ml/min/ml; p = 0.026). Corrected for rate-pressure-product, baseline findings maintained significance (1.41 ± 0.33 vs. 1.16 ± 0.19 ml/min/ml; p = 0.024). CVR was lower in women at baseline (81 ± 14 vs. 107 ± 22 mmHg*ml(-1)*min*ml; p = 0.006) and during cold-pressor-testing (71 ± 17 vs. 91 ± 20 mmHg*ml(-1)*min*ml; p = 0.013). Under adenosine neither maximal MBF (4.06 ± 1.0 vs. 3.91 ± 0.88 ml/min/ml; p = ns) nor coronary flow reserve (3.07 ± 1.12 vs. 3.44 ± 0.92; p = ns) nor CVR (24 ± 8 vs. 24 ± 6 mmHg*ml(-1)*min*ml; p = ns) showed sex-related differences. CONCLUSION: Women show higher myocardial perfusion and lower coronary vascular resistance than men in physiologic states. Maximum perfusion and vasodilation under adenosine are not sex-specific.

Impact of Data-driven Respiratory Gating in Clinical PET.

Purpose To study the feasibility and impact of respiratory gating in positron emission tomographic (PET) imaging in a clinical trial comparing conventional hardware-based gating with a data-driven approach and to describe the distribution of determined parameters. Materials and Methods This prospective study was approved by the ethics committee of the University Hospital of Münster (AZ 2014-217-f-N). Seventy-four patients suspected of having abdominal or thoracic fluorine 18 fluorodeoxyglucose (FDG)-positive lesions underwent clinical whole-body FDG PET/computed tomographic (CT) examinations. Respiratory gating was performed by using a pressure-sensitive belt system (belt gating [BG]) and an automatic data-driven approach (data-driven gating [DDG]). PET images were analyzed for lesion uptake, metabolic volumes, respiratory shifts of lesions, and diagnostic image quality. Results Forty-eight patients had at least one lesion in the field of view, resulting in a total of 164 lesions analyzed (range of number of lesions per patient, one to 13). Both gating methods revealed respiratory shifts of lesions (4.4 mm ± 3.1 for BG vs 4.8 mm ± 3.6 for DDG, P = .76). Increase in uptake of the lesions compared with nongated values did not differ significantly between both methods (maximum standardized uptake value [SUVmax], +7% ± 13 for BG vs +8% ± 16 for DDG, P = .76). Similarly, gating significantly decreased metabolic lesion volumes with both methods (-6% ± 26 for BG vs -7% ± 21 for DDG, P = .44) compared with nongated reconstructions. Blinded reading revealed significant improvements in diagnostic image quality when using gating, without significant differences between the methods (DDG was judged to be inferior to BG in 22 cases, equal in 12 cases, and superior in 15 cases; P = .32). Conclusion Respiratory gating increases diagnostic image quality and uptake values and decreases metabolic volumes compared with nongated acquisitions. Data-driven approaches are clinically applicable alternatives to belt-based methods and might help establishing routine respiratory gating in clinical PET/CT. (©) RSNA, 2016 Online supplemental material is available for this article.

A contactless approach for respiratory gating in PET using continuous-wave radar.

PURPOSE: Respiratory gating is commonly used to reduce motion artifacts in positron emission tomography (PET). Clinically established methods for respiratory gating in PET require contact to the patient or a direct optical line between the sensor and the patient's torso and time consuming preparation. In this work, a contactless method for capturing a respiratory signal during PET is presented based on continuous-wave radar. METHODS: The proposed method relies on the principle of emitting an electromagnetic wave and detecting the phase shift of the reflected wave, modulated due to the respiratory movement of the patient's torso. A 24 GHz carrier frequency was chosen allowing wave propagation through plastic and clothing with high reflections at the skin surface. A detector module and signal processing algorithms were developed to extract a quantitative respiratory signal. The sensor was validated using a high precision linear table. During volunteer measurements and [(18)F] FDG PET scans, the radar sensor was positioned inside the scanner bore of a PET/computed tomography scanner. As reference, pressure belt (one volunteer), depth camera-based (two volunteers, two patients), and PET data-driven (six patients) signals were acquired simultaneously and the signal correlation was quantified. RESULTS: The developed system demonstrated a high measurement accuracy for movement detection within the submillimeter range. With the proposed method, small displacements of 25 µm could be detected, not considerably influenced by clothing or blankets. From the patient studies, the extracted respiratory radar signals revealed high correlation (Pearson correlation coefficient) to those derived from the external pressure belt and depth camera signals (r = 0.69-0.99) and moderate correlation to those of the internal data-driven signals (r = 0.53-0.70). In some cases, a cardiac signal could be visualized, due to the representation of the mechanical heart motion on the skin. CONCLUSIONS: Accurate respiratory signals were obtained successfully by the proposed method with high spatial and temporal resolution. By working without contact and passing through clothing and blankets, this approach minimizes preparation time and increases the convenience of the patient during the scan.

Correction of MRI-induced geometric distortions in whole-body small animal PET-MRI.

PURPOSE: The fusion of positron emission tomography (PET) and magnetic resonance imaging (MRI) data can be a challenging task in whole-body PET-MRI. The quality of the registration between these two modalities in large field-of-views (FOV) is often degraded by geometric distortions of the MRI data. The distortions at the edges of large FOVs mainly originate from MRI gradient nonlinearities. This work describes a method to measure and correct for these kind of geometric distortions in small animal MRI scanners to improve the registration accuracy of PET and MRI data. METHODS: The authors have developed a geometric phantom which allows the measurement of geometric distortions in all spatial axes via control points. These control points are detected semiautomatically in both PET and MRI data with a subpixel accuracy. The spatial transformation between PET and MRI data is determined with these control points via 3D thin-plate splines (3D TPS). The transformation derived from the 3D TPS is finally applied to real MRI mouse data, which were acquired with the same scan parameters used in the phantom data acquisitions. Additionally, the influence of the phantom material on the homogeneity of the magnetic field is determined via field mapping. RESULTS: The spatial shift according to the magnetic field homogeneity caused by the phantom material was determined to a mean of 0.1 mm. The results of the correction show that distortion with a maximum error of 4 mm could be reduced to less than 1 mm with the proposed correction method. Furthermore, the control point-based registration of PET and MRI data showed improved congruence after correction. CONCLUSIONS: The developed phantom has been shown to have no considerable negative effect on the homogeneity of the magnetic field. The proposed method yields an appropriate correction of the measured MRI distortion and is able to improve the PET and MRI registration. Furthermore, the method is applicable to whole-body small animal imaging routines including different standard MRI sequences.

A dual-Kinect approach to determine torso surface motion for respiratory motion correction in PET.

PURPOSE: Respiratory gating is commonly used to reduce blurring effects and attenuation correction artifacts in positron emission tomography (PET). Established clinically available methods that employ body-attached hardware for acquiring respiration signals rely on the assumption that external surface motion and internal organ motion are well correlated. In this paper, the authors present a markerless method comprising two Microsoft Kinects for determining the motion on the whole torso surface and aim to demonstrate its validity and usefulness-including the potential to study the external/internal correlation and to provide useful information for more advanced correction approaches. METHODS: The data of two Kinects are used to calculate 3D representations of a patient's torso surface with high spatial coverage. Motion signals can be obtained for any position by tracking the mean distance to a virtual camera with a view perpendicular to the surrounding surface. The authors have conducted validation experiments including volunteers and a moving high-precision platform to verify the method's suitability for providing meaningful data. In addition, the authors employed it during clinical (18)F-FDG-PET scans and exemplarily analyzed the acquired data of ten cancer patients. External signals of abdominal and thoracic regions as well as data-driven signals were used for gating and compared with respect to detected displacement of present lesions. Additionally, the authors quantified signal similarities and time shifts by analyzing cross-correlation sequences. RESULTS: The authors' results suggest a Kinect depth resolution of approximately 1 mm at 75 cm distance. Accordingly, valid signals could be obtained for surface movements with small amplitudes in the range of only few millimeters. In this small sample of ten patients, the abdominal signals were better suited for gating the PET data than the thoracic signals and the correlation of data-driven signals was found to be stronger with abdominal signals than with thoracic signals (average Pearson correlation coefficients of 0.74 ± 0.17 and 0.45 ± 0.23, respectively). In all cases, except one, the abdominal respiratory motion preceded the thoracic motion-a maximum delay of approximately 600 ms was detected. CONCLUSIONS: The method provides motion information with sufficiently high spatial and temporal resolution. Thus, it enables meaningful analysis in the form of comparisons between amplitudes and phase shifts of signals from different regions. In combination with a large field-of-view, as given by combining the data of two Kinect cameras, it yields surface representations that might be useful in the context of motion correction and motion modeling.

Motion correction of whole-body PET data with a joint PET-MRI registration functional.

Respiratory motion is known to degrade image quality in PET imaging. The necessary acquisition time of several minutes per bed position will inevitably lead to a blurring effect due to organ motion. A lot of research has been done with regards to motion correction of PET data. As full-body PET-MRI became available recently, the anatomical data provided by MRI is a promising source of motion information. Current PET-MRI-based motion correction approaches, however, do not take into account the available information provided by PET data. PET data, though, may add valuable additional information to increase motion estimation robustness and precision.

On transcending the impasse of respiratory motion correction applications in routine clinical imaging - a consideration of a fully automated data driven motion control framework.

Positron emission tomography (PET) is increasingly used for the detection, characterization, and follow-up of tumors located in the thorax. However, patient respiratory motion presents a unique limitation that hinders the application of high-resolution PET technology for this type of imaging. Efforts to transcend this limitation have been underway for more than a decade, yet PET remains for practical considerations a modality vulnerable to motion-induced image degradation. Respiratory motion control is not employed in routine clinical operations. In this article, we take an opportunity to highlight some of the recent advancements in data-driven motion control strategies and how they may form an underpinning for what we are presenting as a fully automated data-driven motion control framework. This framework represents an alternative direction for future endeavors in motion control and can conceptually connect individual focused studies with a strategy for addressing big picture challenges and goals.

A mass conservation-based optical flow method for cardiac motion correction in 3D-PET.

PURPOSE: Cardiac positron emission tomography (PET) images usually show two kinds of artifacts: the limited resolution of PET leads to partial volume effects and the motion of the heart induces blurring. These phenomena degrade the PET images and induce errors in the quantification. One method of reducing this problem is to use gated PET data. However, the reduction of information per phase leads to an increase in noise on the reconstructed images. Alternatively, the PET data have to be corrected for motion and partial volume effects. METHODS: Optical flow methods have been shown to accurately estimate the motion between PET image frames. These methods assume that the brightness of the objects remains constant between the frames. This condition is not fulfilled in cardiac PET data because the brightness of the cardiac muscle tissue (myocardium) is not accurately resolved due to the partial volume effect. Therefore, the use of a newly developed optical flow method based upon the conservation of mass condition is proposed to correct the cardiac PET data. Mass conservation is applicable to PET images as the total activity in the field of view may be assumed to remain almost constant, if the data are precorrected for radioactive decay. Two variants of the method using the quadratic and the nonquadratic penalization are presented. The methods were evaluated with respect to correlation coefficient, myocardial thickness and the blood pool activity in the left ventricle on phantom data and on 14 patient image volumes. RESULTS: The proposed methods showed that the cardiac motion can be efficiently corrected despite partial volume effects. The correlation coefficient between the image volumes increased from 0.87 to 0.98 on average. The change in myocardial thickness was reduced from 28% to 3%. The variation in blood pool activity was reduced from 80% to 8%. The algorithm needed only about 4 s for execution. CONCLUSIONS: A mass preserving optical flow method of cardiac motion correction in 3D PET data has been presented and tested on phantom as well as patient data. The results show that the motion was corrected for all datasets effectively.