Self-gated MRI motion modeling for respiratory motion compensation in integrated PET/MRI

Abdominal MR Multitasking for radiotherapy treatment planning: A motion-resolved and multicontrast 3D imaging approach

PURPOSE

Radiotherapy treatment planning (RTP) using MR has been used increasingly for the abdominal site. Multiple contrast weightings and motion-resolved imaging are desired for accurate delineation of the target and various organs-at-risk and patient-tailored planning. Current MR protocols achieve these through multiple scans with distinct contrast and variable respiratory motion management strategies and acquisition parameters, leading to a complex and inaccurate planning process. This study presents a standalone MR Multitasking (MT)-based technique to produce volumetric, motion-resolved, multicontrast images for abdominal radiotherapy treatment planning.

METHODS

The MT technique resolves motion and provides a wide range of contrast weightings by repeating a magnetization-prepared (saturation recovery and T2 preparations) spoiled gradient-echo readout series and adopting the MT image reconstruction framework. The performance of the technique was assessed through digital phantom simulations and in vivo studies of both healthy volunteers and patients with liver tumors.

RESULTS

In the digital phantom study, the MT technique presented structural details and motion in excellent agreement with the digital ground truth. The in vivo studies showed that the motion range was highly correlated (R2 = 0.82) between MT and 2D cine imaging. MT allowed for a flexible contrast-weighting selection for better visualization. Initial clinical testing with interobserver analysis demonstrated acceptable target delineation quality (Dice coefficient = 0.85 ± 0.05, Hausdorff distance = 3.3 ± 0.72 mm).

CONCLUSION

The developed MT-based, abdomen-dedicated technique is capable of providing motion-resolved, multicontrast volumetric images in a single scan, which may facilitate abdominal radiotherapy treatment planning.

Development of quantitative PET/MR imaging for measurements of hepatic portal vein input function: a phantom study

BACKGROUND

Accurate pharmacokinetic modelling in PET necessitates measurements of an input function, which ideally is acquired non-invasively from image data. For hepatic pharmacokinetic modelling two input functions need to be considered, to account for the blood supply from the hepatic artery and portal vein. Image-derived measurements at the portal vein are challenging due to its small size and image artifacts caused by respiratory motion. In this work we seek to demonstrate, using phantom experiments, how a dedicated PET/MR protocol can tackle these challenges and potentially provide input function measurements of the portal vein in a clinical setup.

METHODS

A custom 3D printed PET/MR phantom was constructed to mimic the liver and portal vein. PET/MR acquisitions were made with emulated respiratory motion. The PET/MR imaging protocol consisted of high-resolution anatomical MR imaging of the portal vein, followed by a PET acquisition in parallel to a dedicated motion-tracking MR sequence. Motion tracking and deformation information were extracted from PET data and subsequently used in PET reconstruction to produce dynamic series of motion-free PET images. Anatomical MR images were used post PET reconstruction for partial volume correction of the input function measurements.

RESULTS

Reconstruction of dynamic PET data with motion-compensation provided nearly motion-free series of PET frame data, suitable for image derived input function measurements of the portal vein. After partial volume correction, the individual input function measurements were within a 16.1% error range from the true activity in the portal vein compartment at the time of PET acquisition.

CONCLUSION

The proposed protocol demonstrates clinically feasible PET/MR imaging of the liver for pharmacokinetic studies with accurate quantification of the portal vein input function, including correction for respiratory motion and partial volume effects.

Free-breathing high-resolution respiratory-gated radial stack-of-stars magnetic resonance imaging of the upper abdomen at 7 T

Ultrahigh field magnetic resonance imaging (MRI) (≥ 7 T) has the potential to provide superior spatial resolution and unique image contrast. Apart from radiofrequency transmit inhomogeneities in the body at this field strength, imaging of the upper abdomen faces additional challenges associated with motion-induced ghosting artifacts. To address these challenges, the goal of this work was to develop a technique for high-resolution free-breathing upper abdominal MRI at 7 T with a large field of view. Free-breathing 3D gradient-recalled echo (GRE) water-excited radial stack-of-stars data were acquired in seven healthy volunteers (five males/two females, body mass index: 19.6-24.8 kg/m2) at 7 T using an eight-channel transceive array coil. Two volunteers were also examined at 3 T. In each volunteer, the liver and kidney regions were scanned in two separate acquisitions. To homogenize signal excitation, the time-interleaved acquisition of modes (TIAMO) method was used with personalized pairs of B1 shims, based on a 23-s Cartesian fast low angle shot (FLASH) acquisition. Utilizing free-induction decay navigator signals, respiratory-gated images were reconstructed at a spatial resolution of 0.8 × 0.8 × 1.0 mm3. Two experienced radiologists rated the image quality and the impact of B1 inhomogeneity and motion-related artifacts on multipoint scales. The images of all volunteers showcased effective water excitation and were accurately corrected for respiratory motion. The impact of B1 inhomogeneity on image quality was minimal, underscoring the efficacy of the multitransmit TIAMO shim. The high spatial resolution allowed excellent depiction of small structures such as the adrenal glands, the proximal ureter, the diaphragm, and small blood vessels, although some streaking artifacts persisted in liver image data. In direct comparisons with 3 T performed for two volunteers, 7-T acquisitions demonstrated increases in signal-to-noise ratio of 77% and 58%. Overall, this work demonstrates the feasibility of free-breathing MRI in the upper abdomen at submillimeter spatial resolution at a magnetic field strength of 7 T.

Data-driven system matrix manipulation enabling fast functional imaging and intra-image nonrigid motion correction in tomography

Tomographic imaging modalities are described by large system matrices. Sparse sampling and tissue motion degrade system matrix and image quality. Various existing techniques improve the image quality without correcting the system matrices. Here, we compress the system matrices to improve computational efficiency (e.g., 42 times) using singular value decomposition and fast Fourier transform. Enabled by the efficiency, we propose (1) fast sparsely sampling functional imaging by incorporating a densely sampled prior image into the system matrix, which maintains the critical linearity while mitigating artifacts and (2) intra-image nonrigid motion correction by incorporating the motion as subdomain translations into the system matrix and reconstructing the translations together with the image iteratively. We demonstrate the methods in 3D photoacoustic computed tomography with significantly improved image qualities and clarify their applicability to X-ray CT and MRI or other types of imperfections due to the similarities in system matrices.

Multiparametric Oncologic Hybrid Imaging: Machine Learning Challenges and Opportunities

BACKGROUND

Machine learning (ML) is considered an important technology for future data analysis in health care.

METHODS

The inherently technology-driven fields of diagnostic radiology and nuclear medicine will both benefit from ML in terms of image acquisition and reconstruction. Within the next few years, this will lead to accelerated image acquisition, improved image quality, a reduction of motion artifacts and - for PET imaging - reduced radiation exposure and new approaches for attenuation correction. Furthermore, ML has the potential to support decision making by a combined analysis of data derived from different modalities, especially in oncology. In this context, we see great potential for ML in multiparametric hybrid imaging and the development of imaging biomarkers.

RESULTS AND CONCLUSION

In this review, we will describe the basics of ML, present approaches in hybrid imaging of MRI, CT, and PET, and discuss the specific challenges associated with it and the steps ahead to make ML a diagnostic and clinical tool in the future.

KEY POINTS

· ML provides a viable clinical solution for the reconstruction, processing, and analysis of hybrid imaging obtained from MRI, CT, and PET..

Respiratory motion management using a single rapid MRI scan for a 0.35 T MRI-Linac system

BACKGROUND

MRI has a rapidly growing role in radiation therapy (RT) for treatment planning, real-time image guidance, and beam gating (e.g., MRI-Linac). Free-breathing 4D-MRI is desirable in respiratory motion management for therapy. Moreover, high-quality 3D-MRIs without motion artifacts are needed to delineate lesions. Existing MRI methods require multiple scans with lengthy acquisition times or are limited by low spatial resolution, contrast, and signal-to-noise ratio.

PURPOSE

We developed a novel method to obtain motion-resolved 4D-MRIs and motion-integrated 3D-MRI reconstruction using a single rapid (35-45 s scan on a 0.35 T MRI-Linac.

METHODS

Golden-angle radial stack-of-stars MRI scans were acquired from a respiratory motion phantom and 12 healthy volunteers (n = 12) on a 0.35 T MRI-Linac. A self-navigated method was employed to detect respiratory motion using 2000 (acquisition time = 5-7 min) and the first 200 spokes (acquisition time = 35-45 s). Multi-coil non-uniform fast Fourier transform (MCNUFFT), compressed sensing (CS), and deep-learning Phase2Phase (P2P) methods were employed to reconstruct motion-resolved 4D-MRI using 2000 spokes (MCNUFFT2000) and 200 spokes (CS200 and P2P200). Deformable motion vector fields (MVFs) were computed from the 4D-MRIs and used to reconstruct motion-corrected 3D-MRIs with the MOtion Transformation Integrated forward-Fourier (MOTIF) method. Image quality was evaluated quantitatively using the structural similarity index measure (SSIM) and the root mean square error (RMSE), and qualitatively in a blinded radiological review.

RESULTS

Evaluation using the respiratory motion phantom experiment showed that the proposed method reversed the effects of motion blurring and restored edge sharpness. In the human study, P2P200 had smaller inaccuracy in MVFs estimation than CS200. P2P200 had significantly greater SSIMs (p < 0.0001) and smaller RMSEs (p < 0.001) than CS200 in motion-resolved 4D-MRI and motion-corrected 3D-MRI. The radiological review found that MOTIF 3D-MRIs using MCNUFFT2000 exhibited the highest image quality (scoring > 8 out of 10), followed by P2P200 (scoring > 5 out of 10), and then motion-uncorrected (scoring < 3 out of 10) in sharpness, contrast, and artifact-freeness.

CONCLUSIONS

We have successfully demonstrated a method for respiratory motion management for MRI-guided RT. The method integrated self-navigated respiratory motion detection, deep-learning P2P 4D-MRI reconstruction, and a motion integrated reconstruction (MOTIF) for 3D-MRI using a single rapid MRI scan (35-45 s) on a 0.35 T MRI-Linac system.

Quality assurance for magnetic resonance angiography of the chest in patients suspected of pulmonary embolism during iodinated contrast shortage in the emergency department setting

Objectives

COVID-19 lockdowns resulted in a global shortage of iodinated contrast media. Therefore, alternative imaging protocols were devised to evaluate patients arriving to the emergency department (ED) with suspicion of pulmonary embolism (PE). This quality assurance (QA) aims to compare diagnostic potential between alternative magnetic resonance angiography (MRA) protocol over the gold standard computed tomography angiography (CTA) by evaluating MRA imaging quality, scanner type/imaging sequence, and any risk of misdiagnosis in patients with symptoms of PE.

Material and Methods

This retrospective study compromised of 55 patients who arrived to ED and underwent MRA of the chest for suspicion of PE during the months of May to June 2022. Data regarding their chief complaints, imaging sequence, and MRA results were collected. Two fellowship-trained faculty radiologists reviewed the MRA scans of the patients and scored the quality using a Likert scale.

Results

Two patients were positive for PE and 53 patients showed negative results. Regarding the scan quality issues, motion was noted in 80% of the 55 studies that we reviewed. Significant associations (P < 0.009) between Likert scale scores and initial complaint category were found. The characteristic symptoms associated with suspicion of PE, namely, shortness of breath, chest pain, and cough were distributed among the 1 and 2 categories, reflecting the most optimal vessel opacification scores. We found no risk of misdiagnosis after reviewing the electronic medical record for follow-up appointments within 6 months of ED visit.

Conclusion

Patients were screened for PE with MRA as an alternative imaging tool during times of contrast shortage. Further, evaluation of MRA with CTA, side by side, in a larger patient population is required to increase the validity of our QA study.

Anatomical validation of automatic respiratory motion correction for coronary 18F-sodium fluoride positron emission tomography by expert measurements from four-dimensional computed tomography

BACKGROUND

Respiratory motion correction is of importance in studies of coronary plaques employing ¹⁸ F-NaF; however, the validation of motion correction techniques mainly relies on indirect measures such as test-retest repeatability assessments. In this study, we aim to compare and, thus, validate the respiratory motion vector fields obtained from the positron emission tomography (PET) images directly to the respiratory motion observed during four-dimensional cine-computed tomography (CT) by an expert observer.

PURPOSE

To investigate the accuracy of the motion correction employed in a software (FusionQuant) used for evaluation of ¹⁸ F-NaF PET studies by comparing the respiratory motion of the coronary plaques observed in PET to the respiratory motion observed in 4D cine-CT images.

METHODS

This study included 23 patients who undertook thoracic PET scans for the assessment of coronary plaques using ¹⁸ F-sodium fluoride (¹⁸ F-NaF). All patients underwent a 5-s cine-CT (4D-CT), a coronary CT angiography (CTA), and ¹⁸ F-NaF PET. The 4D-CT and PET scan were reconstructed into 10 phases. Respiratory motion was estimated for the non-contrast visible coronary plaques using diffeomorphic registrations (PET) and compared to respiratory motion observed on 4D-CT. We report the PET motion vector fields obtained in the three principal axes in addition to the 3D motion. Statistical differences were examined using paired t-tests. Signal-to-noise ratios (SNR) are reported for the single-phase images (end-expiratory phase) and for the motion-corrected image-series (employing the motion vector fields extracted during the diffeomorphic registrations).

RESULTS

In total, 19 coronary plaques were identified in 16 patients. No statistical differences were observed for the maximum respiratory motion observed in x, y, and the 3D motion fields (magnitude and direction) between the CT and PET (X direction: 4D CT = 2.5 ± 1.5 mm, PET = 2.4 ± 3.2 mm; Y direction: 4D CT = 2.3 ± 1.9 mm, PET = 0.7 ± 2.9 mm, 3D motion: 4D CT = 6.6 ± 3.1 mm, PET = 5.7 ± 2.6 mm, all p ≥ 0.05). Significant differences in respiratory motion were observed in the systems' Z direction: 4D CT = 4.9 ± 3.4 mm, PET = 2.3 ± 3.2 mm, p = 0.04. Significantly improved SNR is reported for the motion corrected images compared to the end-expiratory phase images (end-expiratory phase = 6.8±4.8, motion corrected = 12.2±4.5, p = 0.001).

CONCLUSION

Similar respiratory motion was observed in two directions and 3D for coronary plaques on 4D CT as detected by automatic respiratory motion correction of coronary PET using FusionQuant. The respiratory motion correction technique significantly improved the SNR in the images.

Deep-learning synthesized pseudo-CT for MR high-resolution pediatric cranial bone imaging (MR-HiPCB)

PURPOSE

CT is routinely used to detect cranial abnormalities in pediatric patients with head trauma or craniosynostosis. This study aimed to develop a deep learning method to synthesize pseudo-CT (pCT) images for MR high-resolution pediatric cranial bone imaging to eliminating ionizing radiation from CT.

METHODS

3D golden-angle stack-of-stars MRI were obtained from 44 pediatric participants. Two patch-based residual UNets were trained using paired MR and CT patches randomly selected from the whole head (NetWH) or in the vicinity of bone, fractures/sutures, or air (NetBA) to synthesize pCT. A third residual UNet was trained to generate a binary brain mask using only MRI. The pCT images from NetWH (pCTNetWH ) in the brain area and NetBA (pCTNetBA ) in the nonbrain area were combined to generate pCTCom . A manual processing method using inverted MR images was also employed for comparison.

RESULTS

pCTCom (68.01 ± 14.83 HU) had significantly smaller mean absolute errors (MAEs) than pCTNetWH (82.58 ± 16.98 HU, P < 0.0001) and pCTNetBA (91.32 ± 17.2 HU, P < 0.0001) in the whole head. Within cranial bone, the MAE of pCTCom (227.92 ± 46.88 HU) was significantly lower than pCTNetWH (287.85 ± 59.46 HU, P < 0.0001) but similar to pCTNetBA (230.20 ± 46.17 HU). Dice similarity coefficient of the segmented bone was significantly higher in pCTCom (0.90 ± 0.02) than in pCTNetWH (0.86 ± 0.04, P < 0.0001), pCTNetBA (0.88 ± 0.03, P < 0.0001), and inverted MR (0.71 ± 0.09, P < 0.0001). Dice similarity coefficient from pCTCom demonstrated significantly reduced age dependence than inverted MRI. Furthermore, pCTCom provided excellent suture and fracture visibility comparable to CT.

CONCLUSION

MR high-resolution pediatric cranial bone imaging may facilitate the clinical translation of a radiation-free MR cranial bone imaging method for pediatric patients.

MR-assisted PET respiratory motion correction using deep-learning based short-scan motion fields

PURPOSE

We evaluated the impact of PET respiratory motion correction (MoCo) in a phantom and patients. Moreover, we proposed and examined a PET MoCo approach using motion vector fields (MVFs) from a deep-learning reconstructed short MRI scan.

METHODS

The evaluation of PET MoCo was performed in a respiratory motion phantom study with varying lesion sizes and tumor to background ratios (TBRs) using a static scan as the ground truth. MRI-based MVFs were derived from either 2000 spokes (MoCo2000 , 5-6 min acquisition time) using a Fourier transform reconstruction or 200 spokes (MoCoP2P200 , 30-40 s acquisition time) using a deep-learning Phase2Phase (P2P) reconstruction and then incorporated into PET MoCo reconstruction. For six patients with hepatic lesions, the performance of PET MoCo was evaluated using quantitative metrics (SUVmax , SUVpeak , SUVmean , lesion volume) and a blinded radiological review on lesion conspicuity.

RESULTS

MRI-assisted PET MoCo methods provided similar results to static scans across most lesions with varying TBRs in the phantom. Both MoCo2000 and MoCoP2P200 PET images had significantly higher SUVmax , SUVpeak , SUVmean and significantly lower lesion volume than non-motion-corrected (non-MoCo) PET images. There was no statistical difference between MoCo2000 and MoCoP2P200 PET images for SUVmax , SUVpeak , SUVmean or lesion volume. Both radiological reviewers found that MoCo2000 and MoCoP2P200 PET significantly improved lesion conspicuity.

CONCLUSION

An MRI-assisted PET MoCo method was evaluated using the ground truth in a phantom study. In patients with hepatic lesions, PET MoCo images improved quantitative and qualitative metrics based on only 30-40 s of MRI motion modeling data.

Cranial vault imaging for pediatric head trauma using a radial VIBE MRI sequence

OBJECTIVE

Head trauma is the most common indication for a CT scan. In this pilot study, the authors assess the feasibility of a 5-minute high-resolution 3D golden-angle (GA) stack-of-stars radial volumetric interpolated breath-hold examination (VIBE) MRI sequence (GA-VIBE) to obtain clinically acceptable cranial bone images and identify cranial vault fractures compared to CT.

METHODS

Patients younger than 18 years of age presenting after head trauma were eligible for the study. Three clinicians reviewed and assessed 1) slice-by-slice volumetric CT and inverted MR images, and 2) 3D reconstructions obtained from inverted MR images and the gold standard (CT). For each image set, reviewers noted on 5-point Likert scales whether they recommended that a repeat scan be performed and the presence or absence of cranial vault fractures.

RESULTS

Thirty-one patients completed MRI after a clinical head CT scan was performed. Based on CT imaging, 8 of 31 patients had cranial fractures. Two of 31 patients were sedated as part of their clinical MRI scan. In 30 (97%) of 31 MRI reviews, clinicians agreed (or strongly agreed) that the image quality was acceptable for clinical diagnosis. Overall, comparing MRI to acceptable gold-standard CT, sensitivity and specificity of fracture detection were 100%. Furthermore, there were no discrepancies between CT and MRI in classification of fracture type or location.

CONCLUSIONS

When compared with the gold standard (CT), the volumetric and 3D reconstructed images using the GA-VIBE sequence were able to produce clinically acceptable cranial images with excellent ability to detect cranial vault fractures.

Evaluating different methods of MR-based motion correction in simultaneous PET/MR using a head phantom moved by a robotic system

Background

Due to comparatively long measurement times in simultaneous positron emission tomography and magnetic resonance (PET/MR) imaging, patient movement during the measurement can be challenging. This leads to artifacts which have a negative impact on the visual assessment and quantitative validity of the image data and, in the worst case, can lead to misinterpretations. Simultaneous PET/MR systems allow the MR-based registration of movements and enable correction of the PET data. To assess the effectiveness of motion correction methods, it is necessary to carry out measurements on phantoms that are moved in a reproducible way. This study explores the possibility of using such a phantom-based setup to evaluate motion correction strategies in PET/MR of the human head.

Method

An MR-compatible robotic system was used to generate rigid movements of a head-like phantom. Different tools, either from the manufacturer or open-source software, were used to estimate and correct for motion based on the PET data itself (SIRF with SPM and NiftyReg) and MR data acquired simultaneously (e.g. MCLFIRT, BrainCompass). Different motion estimates were compared using data acquired during robot-induced motion. The effectiveness of motion correction of PET data was evaluated by determining the segmented volume of an activity-filled flask inside the phantom. In addition, the segmented volume was used to determine the centre-of-mass and the change in maximum activity concentration.

Results

The results showed a volume increase between 2.7 and 36.3% could be induced by the experimental setup depending on the motion pattern. Both, BrainCompass and MCFLIRT, produced corrected PET images, by reducing the volume increase to 0.7–4.7% (BrainCompass) and to -2.8–0.4% (MCFLIRT). The same was observed for example for the centre-of-mass, where the results show that MCFLIRT (0.2–0.6 mm after motion correction) had a smaller deviation from the reference position than BrainCompass (0.5–1.8 mm) for all displacements.

Conclusions

The experimental setup is suitable for the reproducible generation of movement patterns. Using open-source software for motion correction is a viable alternative to the vendor-provided motion-correction software.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-022-00442-6.

Cardiac MR: From Theory to Practice

Cardiovascular disease (CVD) is the leading single cause of morbidity and mortality, causing over 17. 9 million deaths worldwide per year with associated costs of over $800 billion. Improving prevention, diagnosis, and treatment of CVD is therefore a global priority. Cardiovascular magnetic resonance (CMR) has emerged as a clinically important technique for the assessment of cardiovascular anatomy, function, perfusion, and viability. However, diversity and complexity of imaging, reconstruction and analysis methods pose some limitations to the widespread use of CMR. Especially in view of recent developments in the field of machine learning that provide novel solutions to address existing problems, it is necessary to bridge the gap between the clinical and scientific communities. This review covers five essential aspects of CMR to provide a comprehensive overview ranging from CVDs to CMR pulse sequence design, acquisition protocols, motion handling, image reconstruction and quantitative analysis of the obtained data. (1) The basic MR physics of CMR is introduced. Basic pulse sequence building blocks that are commonly used in CMR imaging are presented. Sequences containing these building blocks are formed for parametric mapping and functional imaging techniques. Commonly perceived artifacts and potential countermeasures are discussed for these methods. (2) CMR methods for identifying CVDs are illustrated. Basic anatomy and functional processes are described to understand the cardiac pathologies and how they can be captured by CMR imaging. (3) The planning and conduct of a complete CMR exam which is targeted for the respective pathology is shown. Building blocks are illustrated to create an efficient and patient-centered workflow. Further strategies to cope with challenging patients are discussed. (4) Imaging acceleration and reconstruction techniques are presented that enable acquisition of spatial, temporal, and parametric dynamics of the cardiac cycle. The handling of respiratory and cardiac motion strategies as well as their integration into the reconstruction processes is showcased. (5) Recent advances on deep learning-based reconstructions for this purpose are summarized. Furthermore, an overview of novel deep learning image segmentation and analysis methods is provided with a focus on automatic, fast and reliable extraction of biomarkers and parameters of clinical relevance.

Multiparametric Oncologic Hybrid Imaging: Machine Learning Challenges and Opportunities

BACKGROUND

 Machine learning (ML) is considered an important technology for future data analysis in health care.

METHODS

 The inherently technology-driven fields of diagnostic radiology and nuclear medicine will both benefit from ML in terms of image acquisition and reconstruction. Within the next few years, this will lead to accelerated image acquisition, improved image quality, a reduction of motion artifacts and - for PET imaging - reduced radiation exposure and new approaches for attenuation correction. Furthermore, ML has the potential to support decision making by a combined analysis of data derived from different modalities, especially in oncology. In this context, we see great potential for ML in multiparametric hybrid imaging and the development of imaging biomarkers.

RESULTS AND CONCLUSION

 In this review, we will describe the basics of ML, present approaches in hybrid imaging of MRI, CT, and PET, and discuss the specific challenges associated with it and the steps ahead to make ML a diagnostic and clinical tool in the future.

KEY POINTS

  · ML provides a viable clinical solution for the reconstruction, processing, and analysis of hybrid imaging obtained from MRI, CT, and PET..

CITATION FORMAT

· Küstner T, Hepp T, Seith F. Multiparametric Oncologic Hybrid Imaging: Machine Learning Challenges and Opportunities. Fortschr Röntgenstr 2022; DOI: 10.1055/a-1718-4128.

Tumour-to-normal tissue (T/N) dosimetry ratios role in assessment of 90Y selective internal radiation therapy (SIRT)

OBJECTIVE

The purpose of our work is to assess the role of tumour-to-normal tissue (T/N) dosimetry ratios for predicting response in patients undergoing locoregional therapy to the liver with ⁹⁰Y microspheres.

METHODS

A total of 39 patients (7 female:32 male, mean age 68.3 ± 7.6 years), underwent positron emission tomography (PET)/CT imaging after treatment with ⁹⁰Y microspheres. For attenuation correction and localization of the ⁹⁰Y microspheres, the low-dose, non-diagnostic CT images from PET/CT were used. The acquisition took 15 min and the reconstruction matrix size was 200 × 200 × 75 mm and voxel size of 4.07 × 4.07 × 3.00 mm. For dosimetry calculations, the local deposition method with known activity of ⁹⁰Y was used. For each patient, regions of interest for tumour(s) and whole liver were manually created; the normal tissue region of interest was created automatically. mRECIST criteria on MRI done at 1 month post-treatment and subsequently every 3 months after ⁹⁰Y treatment, were used to assess response.

RESULTS

For 39 patients, the mean liver, tumour and normal tissue doses (mean ± SD) were, 55.17 ± 26.04 Gy, 911.87 ± 866.54 Gy and 47.79 ± 20.47 Gy, respectively. Among these patients, 31 (79%) showed complete response (CR) and 8 (21%) showed progression of disease (PD). For patients with CR, the mean T/N dose ratio obtained was 24.91 (range 3.09-80.12) and for patients with PD, the mean T/N dose ratio was significantly lower, at 6.69 (range 0.36-14.75).

CONCLUSION

Our data show that patients with CR have a statistically higher T/N dose ratio than those with PD. Because, the number of PD cases was limited and partial volume effect was not considered, further investigation is warranted.

ADVANCES IN KNOWLEDGE

T/N dosimetry ratios can be used for assessing response in patients undergoing locoregional therapy to the liver with ⁹⁰Y microspheres.

Phase2Phase: Respiratory Motion-Resolved Reconstruction of Free-Breathing Magnetic Resonance Imaging Using Deep Learning Without a Ground Truth for Improved Liver Imaging

OBJECTIVES

Respiratory binning of free-breathing magnetic resonance imaging data reduces motion blurring; however, it exacerbates noise and introduces severe artifacts due to undersampling. Deep neural networks can remove artifacts and noise but usually require high-quality ground truth images for training. This study aimed to develop a network that can be trained without this requirement.

MATERIALS AND METHODS

This retrospective study was conducted on 33 participants enrolled between November 2016 and June 2019. Free-breathing magnetic resonance imaging was performed using a radial acquisition. Self-navigation was used to bin the k-space data into 10 respiratory phases. To simulate short acquisitions, subsets of radial spokes were used in reconstructing images with multicoil nonuniform fast Fourier transform (MCNUFFT), compressed sensing (CS), and 2 deep learning methods: UNet3DPhase and Phase2Phase (P2P). UNet3DPhase was trained using a high-quality ground truth, whereas P2P was trained using noisy images with streaking artifacts. Two radiologists blinded to the reconstruction methods independently reviewed the sharpness, contrast, and artifact-freeness of the end-expiration images reconstructed from data collected at 16% of the Nyquist sampling rate. The generalized estimating equation method was used for statistical comparison. Motion vector fields were derived to examine the respiratory motion range of 4-dimensional images reconstructed using different methods.

RESULTS

A total of 15 healthy participants and 18 patients with hepatic malignancy (50 ± 15 years, 6 women) were enrolled. Both reviewers found that the UNet3DPhase and P2P images had higher contrast (P < 0.01) and fewer artifacts (P < 0.01) than the CS images. The UNet3DPhase and P2P images were reported to be sharper than the CS images by 1 reviewer (P < 0.01) but not by the other reviewer (P = 0.22, P = 0.18). UNet3DPhase and P2P were similar in sharpness and contrast, whereas UNet3DPhase had fewer artifacts (P < 0.01). The motion vector lengths for the MCNUFFT800 and P2P800 images were comparable (10.5 ± 4.2 mm and 9.9 ± 4.0 mm, respectively), whereas both were significantly larger than CS2000 (7.0 ± 3.9 mm; P < 0.0001) and UNnet3DPhase800 (6.9 ± 3.2; P < 0.0001) images.

CONCLUSIONS

Without a ground truth, P2P can reconstruct sharp, artifact-free, and high-contrast respiratory motion-resolved images from highly undersampled data. Unlike the CS and UNet3DPhase methods, P2P did not artificially reduce the respiratory motion range.

Synergistic motion compensation strategies for positron emission tomography when acquired simultaneously with magnetic resonance imaging

Subject motion in positron emission tomography (PET) is a key factor that degrades image resolution and quality, limiting its potential capabilities. Correcting for it is complicated due to the lack of sufficient measured PET data from each position. This poses a significant barrier in calculating the amount of motion occurring during a scan. Motion correction can be implemented at different stages of data processing either during or after image reconstruction, and once applied accurately can substantially improve image quality and information accuracy. With the development of integrated PET-MRI (magnetic resonance imaging) scanners, internal organ motion can be measured concurrently with both PET and MRI. In this review paper, we explore the synergistic use of PET and MRI data to correct for any motion that affects the PET images. Different types of motion that can occur during PET-MRI acquisitions are presented and the associated motion detection, estimation and correction methods are reviewed. Finally, some highlights from recent literature in selected human and animal imaging applications are presented and the importance of motion correction for accurate kinetic modelling in dynamic PET-MRI is emphasized. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'.

Motion estimation and correction for simultaneous PET/MR using SIRF and CIL

SIRF is a powerful PET/MR image reconstruction research tool for processing data and developing new algorithms. In this research, new developments to SIRF are presented, with focus on motion estimation and correction. SIRF's recent inclusion of the adjoint of the resampling operator allows gradient propagation through resampling, enabling the MCIR technique. Another enhancement enabled registering and resampling of complex images, suitable for MRI. Furthermore, SIRF's integration with the optimization library CIL enables the use of novel algorithms. Finally, SPM is now supported, in addition to NiftyReg, for registration. Results of MR and PET MCIR reconstructions are presented, using FISTA and PDHG, respectively. These demonstrate the advantages of incorporating motion correction and variational and structural priors. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'.

A decade of multi-modality PET and MR imaging in abdominal oncology

OBJECTIVES

To investigate trends observed in a decade of published research on multimodality PET(/CT)+MR imaging in abdominal oncology, and to explore how these trends are reflected by the use of multimodality imaging performed at our institution.

METHODS

First, we performed a literature search (2009-2018) including all papers published on the multimodality combination of PET(/CT) and MRI in abdominal oncology. Retrieved papers were categorized according to a structured labelling system, including study design and outcome, cancer and lesion type under investigation and PET-tracer type. Results were analysed using descriptive statistics and evolutions over time were plotted graphically. Second, we performed a descriptive analysis of the numbers of MRI, PET/CT and multimodality PET/CT+MRI combinations (performed within a ≤14 days interval) performed during a similar time span at our institution.

RESULTS

Published research papers involving multimodality PET(/CT)+MRI combinations showed an impressive increase in numbers, both for retrospective combinations of PET/CT and MRI, as well as hybrid PET/MRI. Main areas of research included new PET-tracers, visual PET(/CT)+MRI assessment for staging, and (semi-)quantitative analysis of PET-parameters compared to or combined with MRI-parameters as predictive biomarkers. In line with literature, we also observed a vast increase in numbers of multimodality PET/CT+MRI imaging in our institutional data.

CONCLUSIONS

The tremendous increase in published literature on multimodality imaging, reflected by our institutional data, shows the continuously growing interest in comprehensive multivariable imaging evaluations to guide oncological practice.

ADVANCES IN KNOWLEDGE

The role of multimodality imaging in oncology is rapidly evolving. This paper summarizes the main applications and recent developments in multimodality imaging, with a specific focus on the combination of PET+MRI in abdominal oncology.

Pitfalls on PET/MRI

A decade of PET/MRI clinical imaging has passed and many of the pitfalls are similar to those on earlier studies. However, techniques to overcome them have emerged and continue to develop. Although clinically significant lung nodules are demonstrable, smaller nodules may be detected using ultrashort/zero echo-time (TE) lung MRI. Fast reconstruction ultrashort TE sequences have also been used to achieve high-resolution lung MRI even with free-breathing. The introduction and improvement of time-of-flight scanners and increasing the axial length of the PET detector arrays have more than doubled the sensitivity of the PET part of the system. MRI for attenuation correction has provided many potential pitfalls, including misclassification of tissue classes based on MRI information for attenuation correction. Although the use of short echo times have helped to address these pitfalls, one of the most exciting developments has been the use of deep learning algorithms and computational neural networks to rapidly provide soft tissue, fat, bone and air information for the attenuation correction as a supplement to the attenuation correction information from fat-water imaging. Challenges with motion correction, particularly respiratory and cardiac remain but are being addressed with respiratory monitors and using PET data. In order to address truncation artefacts, the system manufacturers have developed methods to extend the MR field-of-view for the purpose of the attenuation and scatter corrections. General pitfalls like stitching of body sections for individual studies, optimum delivery of images for viewing and reporting, and resource implications for the sheer volume of data generated remain Methods to overcome these pitfalls serve as a strong foundation for the future of PET/MRI. Advances in the underlying technology with significant evolution in hard-ware and software and the exiting developments in use of deep learning algorithms and computational neural networks will drive the next decade of PET/MRI imaging.

Quantitative comparison of data-driven gating and external hardware gating for 18F-FDG PET-MRI in patients with esophageal tumors

Background

Respiratory motion during PET imaging reduces image quality. Data-driven gating (DDG) based on principal component analysis (PCA) can be used to identify respiratory signals. The use of DDG, without need for external devices, would greatly increase the feasibility of using respiratory gating in a routine clinical setting. The objective of this study was to evaluate data-driven gating in relation to external hardware gating and regular static image acquisition on PET-MRI data with respect to SUV_max and lesion volumes.

Methods

Sixteen patients with esophageal or gastroesophageal cancer (Siewert I and II) underwent a 6-min PET scan on a Signa PET-MRI system (GE Healthcare) 1.5–2 h after injection of 4 MBq/kg ¹⁸F-FDG. External hardware gating was done using a respiratory bellow device, and DDG was performed using MotionFree (GE Healthcare). The DDG raw data files and the external hardware-gating raw files were created on a Matlab-based toolbox from the whole 6-min scan LIST-file. For comparison, two 3-min static raw files were created for each patient. Images were reconstructed using TF-OSEM with resolution recovery with 2 iterations, 28 subsets, and 3-mm post filter. SUV_max and lesion volume were measured in all visible lesions, and noise level was measured in the liver. Paired t-test, linear regression, Pearson correlation, and Bland-Altman analysis were used to investigate difference, correlation, and agreement between the methods.

Results

A total number of 30 lesions were included in the study. No significant differences between DDG and external hardware-gating SUV_max or lesion volumes were found, but the noise level was significantly reduced in the DDG images. Both DDG and external hardware gating demonstrated significantly higher SUV_max (9.4% for DDG, 10.3% for external hardware gating) and smaller lesion volume (− 5.4% for DDG, − 6.6% for external gating) in comparison with non-gated static images.

Conclusions

Data-driven gating with MotionFree for PET-MRI performed similar to external device gating for esophageal lesions with respect to SUV_max and lesion volume. Both gating methods significantly increased the SUV_max and reduced the lesion volume in comparison with non-gated static acquisition. DDG resulted in reduced image noise compared to external device gating and static images.

Adaptive Weighting Landmark-Based Group-Wise Registration on Lung DCE-MRI Images

Image registration of lung dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) is challenging because the rapid changes in intensity lead to non-realistic deformations of intensity-based registration methods. To address this problem, we propose a novel landmark-based registration framework by incorporating landmark information into a group-wise registration. Robust principal component analysis is used to separate motion from intensity changes caused by a contrast agent. Landmark pairs are detected on the resulting motion components and then incorporated into an intensity-based registration through a constraint term. To reduce the negative effect of inaccurate landmark pairs on registration, an adaptive weighting landmark constraint is proposed. The method for calculating landmark weights is based on an assumption that the displacement of a good matched landmark is consistent with those of its neighbors. The proposed method was tested on 20 clinical lung DCE-MRI image series. Both visual inspection and quantitative assessment are used for the evaluation. Experimental results show that the proposed method effectively reduces the non-realistic deformations in registration and improves the registration performance compared with several state-of-the-art registration methods.

Motion correction for PET data using subspace-based real-time MR imaging in simultaneous PET/MR

Image quality of positron emission tomography (PET) reconstructions is degraded by subject motion occurring during the acquisition. Magnetic resonance (MR)-based motion correction approaches have been studied for PET/MR scanners and have been successful at capturing regular motion patterns, when used in conjunction with surrogate signals (e.g. navigators) to detect motion. However, handling irregular respiratory motion and bulk motion remains challenging. In this work, we propose an MR-based motion correction method relying on subspace-based real-time MR imaging to estimate motion fields used to correct PET reconstructions. We take advantage of the low-rank characteristics of dynamic MR images to reconstruct high-resolution MR images at high frame rates from highly undersampled k-space data. Reconstructed dynamic MR images are used to determine motion phases for PET reconstruction and estimate phase-to-phase nonrigid motion fields able to capture complex motion patterns such as irregular respiratory and bulk motion. MR-derived binning and motion fields are used for PET reconstruction to generate motion-corrected PET images. The proposed method was evaluated on in vivo data with irregular motion patterns. MR reconstructions accurately captured motion, outperforming state-of-the-art dynamic MR reconstruction techniques. Evaluation of PET reconstructions demonstrated the benefits of the proposed method in terms of motion artifacts reduction, improving the contrast-to-noise ratio by up to a factor 3 and achieveing a target-to-background ratio up to 90% superior compared to standard/uncorrected methods. The proposed method can improve the image quality of motion-corrected PET reconstructions in clinical applications.

Data-driven, projection-based respiratory motion compensation of PET data for cardiac PET/CT and PET/MR imaging

BACKGROUND

Respiratory patient motion causes blurring of the PET images that may impact accurate quantification of perfusion and infarction extents in PET myocardial viability studies. In this study, we investigate the feasibility of correcting for respiratory motion directly in the PET-listmode data prior to image reconstruction using a data-driven, projection-based, respiratory motion compensation (DPR-MoCo) technique.

METHODS

The DPR-MoCo method was validated using simulations of a XCAT phantom (Biograph mMR PET/MR) as well as experimental phantom acquisitions (Biograph mCT PET/CT). Seven patient studies following a dual-tracer (18F-FDG/13N-NH3) imaging-protocol using a PET/MR-system were also evaluated. The performance of the DPR-MoCo method was compared against reconstructions of the acquired data (No-MoCo), a reference gate method (gated) and an image-based MoCo method using the standard reconstruction-transform-average (RTA-MoCo) approach. The target-to-background ratio (TBRLV) in the myocardium and the noise in the liver (CoVliver) were evaluated for all acquisitions. For all patients, the clinical effect of the DPR-MoCo was assessed based on the end-systolic (ESV), the end-diastolic volumes (EDV) and the left ventricular ejection fraction (EF) which were compared to functional values obtained from the cardiac MR.

RESULTS

The DPR-MoCo and the No-MoCo images presented with similar noise-properties (CoV) (P = .12), while the RTA-MoCo and reference-gate images showed increased noise levels (P = .05). TBRLV values increased for the motion limited reconstructions when compared to the No-MoCo reconstructions (P > .05). DPR-MoCo results showed higher correlation with the functional values obtained from the cardiac MR than the No-MoCo results, though non-significant (P > .05).

CONCLUSION

The projection-based DPR-MoCo method helps to improve PET image quality of the myocardium without the need for external devices for motion tracking.

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.

3D pediatric cranial bone imaging using high-resolution MRI for visualizing cranial sutures: a pilot study

OBJECTIVE

There is an unmet need to perform imaging in young children and obtain CT-equivalent cranial bone images without subjecting the patients to radiation. In this study, the authors propose using a high-resolution fast low-angle shot golden-angle 3D stack-of-stars radial volumetric interpolated breath-hold examination (GA-VIBE) MRI sequence that is intrinsically robust to motion and has enhanced bone versus soft-tissue contrast.

METHODS

Patients younger than 11 years of age, who underwent clinical head CT scanning for craniosynostosis or other cranial malformations, were eligible for the study. 3D reconstructed images created from the GA-VIBE MRI sequence and the gold-standard CT scan were randomized and presented to 3 blinded reviewers. For all image sets, each reviewer noted the presence or absence of the 6 primary cranial sutures and recorded on 5-point Likert scales whether they recommended a second scan be performed.

RESULTS

Eleven patients (median age 1.8 years) underwent MRI after clinical head CT scanning was performed. Five of the 11 patients were sedated. Three clinicians reviewed the images, and there were no cases, either with CT scans or MR images, in which a reviewer agreed a repeat scan was required for diagnosis or surgical planning. The reviewers reported clear imaging of the regions of interest on 99% of the CT reviews and 96% of the MRI reviews. With CT as the standard, the sensitivity and specificity of the GA-VIBE MRI sequence to detect suture closure were 97% and 96%, respectively (n = 198 sutures read).

CONCLUSIONS

The 3D reconstructed images using the GA-VIBE sequence in comparison to the CT scans created clinically acceptable cranial images capable of detecting cranial sutures. Future directions include reducing the scan time, improving motion correction, and automating postprocessing for clinical utility.

Impact of respiratory motion correction on lesion visibility and quantification in thoracic PET/MR imaging

The impact of a method for MR-based respiratory motion correction of PET data on lesion visibility and quantification in patients with oncologic findings in the lung was evaluated. Twenty patients with one or more lesions in the lung were included. Hybrid imaging was performed on an integrated PET/MR system using 18F-FDG as radiotracer. The standard thoracic imaging protocol was extended by a free-breathing self-gated acquisition of MR data for motion modelling. PET data was acquired simultaneously in list-mode for 5-10 mins. One experienced radiologist and one experienced nuclear medicine specialist evaluated and compared the post-processed data in consensus regarding lesion visibility (scores 1-4, 4 being best), image noise levels (scores 1-3, 3 being lowest noise), SUVmean and SUVmax. Motion-corrected (MoCo) images were additionally compared with gated images. Non-motion-corrected free-breathing data served as standard of reference in this study. Motion correction generally improved lesion visibility (3.19 ± 0.63) and noise ratings (2.95 ± 0.22) compared to uncorrected (2.81 ± 0.66 and 2.95 ± 0.22, respectively) or gated PET data (2.47 ± 0.93 and 1.30 ± 0.47, respectively). Furthermore, SUVs (mean and max) were compared for all methods to estimate their respective impact on the quantification. Deviations of SUVmax were smallest between the uncorrected and the MoCo lesion data (average increase of 9.1% of MoCo SUVs), while SUVmean agreed best for gated and MoCo reconstructions (MoCo SUVs increased by 1.2%). The studied method for MR-based respiratory motion correction of PET data combines increased lesion sharpness and improved lesion activity quantification with high signal-to-noise ratio in a clinical setting. In particular, the detection of small lesions in moving organs such as the lung and liver may thus be facilitated. These advantages justify the extension of the PET/MR imaging protocol by 5-10 minutes for motion correction.

The impact of atlas-based MR attenuation correction on the diagnosis of FDG-PET/MR for Alzheimer's diseases- A simulation study combining multi-center data and ADNI-data

Background

The purpose of this study was to assess the impact of vendor-provided atlas-based MRAC on FDG PET/MR for the evaluation of Alzheimer’s disease (AD) by using simulated images.

Methods

We recruited 47 patients, from two institutions, who underwent PET/CT and PET/MR (GE SIGNA) examination for oncological staging. From the PET raw data acquired on PET/MR, two FDG-PET series were generated, using vendor-provided MRAC (atlas-based) and CTAC. The following simulation steps were performed in MNI space: After spatial normalization and smoothing of the PET datasets, we calculated the error map for each patient, PET_MRAC/PET_CTAC. We multiplied each of these 47 error maps with each of the 203 Alzheimer’s Disease Neuroimaging Initiative (ADNI) cases after the identical normalization and smoothing. This resulted in 203*47 = 9541 datasets. To evaluate the probability of AD in each resulting image, a cumulative t-value was calculated automatically using commercially-available software (PMOD PALZ) which has been used in multiple large cohort studies. The diagnostic accuracy for the discrimination of AD and predicting progression from mild cognitive impairment (MCI) to AD were evaluated in simulated images compared with ADNI original images.

Results

The accuracy and specificity for the discrimination of AD-patients from normal controls were not substantially impaired, but sensitivity was slightly impaired in 5 out of 47 datasets (original vs. error; 83.2% [CI 75.0%-89.0%], 83.3% [CI 74.2%-89.8%] and 83.1% [CI 75.6%-88.3%] vs. 82.7% [range 80.4–85.0%], 78.5% [range 72.9–83.3%,] and 86.1% [range 81.4–89.8%]). The accuracy, sensitivity and specificity for predicting progression from MCI to AD during 2-year follow-up was not impaired (original vs. error; 62.5% [CI 53.3%-69.3%], 78.8% [CI 65.4%-88.6%] and 54.0% [CI 47.0%-69.1%] vs. 64.8% [range 61.5–66.7%], 75.7% [range 66.7–81.8%,] and 59.0% [range 50.8–63.5%]). The worst 3 error maps show a tendency towards underestimation of PET scores.

Conclusion

FDG-PET/MR based on atlas-based MR attenuation correction showed similar diagnostic accuracy to the CT-based method for the diagnosis of AD and the prediction of progression of MCI to AD using commercially-available software, although with a minor reduction in sensitivity.

Free-breathing fat and R2 * quantification in the liver using a stack-of-stars multi-echo acquisition with respiratory-resolved model-based reconstruction

Purpose

To develop a free‐breathing hepatic fat and R2∗ quantification method by extending a previously described stack‐of‐stars model‐based fat‐water separation technique with additional modeling of the transverse relaxation rate R2∗.

Methods

The proposed technique combines motion‐robust radial sampling using a stack‐of‐stars bipolar multi‐echo 3D GRE acquisition with iterative model‐based fat‐water separation. Parallel‐Imaging and Compressed‐Sensing principles are incorporated through modeling of the coil‐sensitivity profiles and enforcement of total‐variation (TV) sparsity on estimated water, fat, and R2∗ parameter maps. Water and fat signals are used to estimate the confounder‐corrected proton‐density fat fraction (PDFF). Two strategies for handling respiratory motion are described: motion‐averaged and motion‐resolved reconstruction. Both techniques were evaluated in patients (n = 14) undergoing a hepatobiliary research protocol at 3T. PDFF and R2∗ parameter maps were compared to a breath‐holding Cartesian reference approach.

Results

Linear regression analyses demonstrated strong (r > 0.96) and significant (P ≪ .01) correlations between radial and Cartesian PDFF measurements for both the motion‐averaged reconstruction (slope: 0.90; intercept: 0.07%) and the motion‐resolved reconstruction (slope: 0.90; intercept: 0.11%). The motion‐averaged technique overestimated hepatic R2∗ values (slope: 0.35; intercept: 30.2 1/s) compared to the Cartesian reference. However, performing a respiratory‐resolved reconstruction led to better R2∗ value consistency (slope: 0.77; intercept: 7.5 1/s).

Conclusions

The proposed techniques are promising alternatives to conventional Cartesian imaging for fat and R2∗ quantification in patients with limited breath‐holding capabilities. For accurate R2∗ estimation, respiratory‐resolved reconstruction should be used.

MR-based cardiac and respiratory motion correction of PET: application to static and dynamic cardiac 18F-FDG imaging

Motion of the myocardium deteriorates the quality and quantitative accuracy of cardiac PET images. We present a method for MR-based cardiac and respiratory motion correction of cardiac PET data and evaluate its impact on estimation of activity and kinetic parameters in human subjects. Three healthy subjects underwent simultaneous dynamic 18F-FDG PET and MRI on a hybrid PET/MR scanner. A cardiorespiratory motion field was determined for each subject using navigator, tagging and golden-angle radial MR acquisitions. Acquired coincidence events were binned into cardiac and respiratory phases using electrocardiogram and list mode-driven signals, respectively. Dynamic PET images were reconstructed with MR-based motion correction (MC) and without motion correction (NMC). Parametric images of 18F-FDG consumption rates (Ki) were estimated using Patlak's method for both MC and NMC images. MC alleviated motion artifacts in PET images, resulting in improved spatial resolution, improved recovery of activity in the myocardium wall and reduced spillover from the myocardium to the left ventricle cavity. Significantly higher myocardium contrast-to-noise ratio and lower apparent wall thickness were obtained in MC versus NMC images. Likewise, parametric images of Ki calculated with MC data had improved spatial resolution as compared to those obtained with NMC. Consistent with an increase in reconstructed activity concentration in the frames used during kinetic analyses, MC led to the estimation of higher Ki values almost everywhere in the myocardium, with up to 18% increase (mean across subjects) in the septum as compared to NMC. This study shows that MR-based motion correction of cardiac PET results in improved image quality that can benefit both static and dynamic studies.

Higher-order singular value decomposition-based lung parcellation for breathing motion management

Positron emission tomography (PET) imaging of the lungs is confounded by respiratory motion-induced blurring artifacts that degrade quantitative accuracy. Gating and motion-compensated image reconstruction are frequently used to correct these motion artifacts in PET. In the absence of voxel-by-voxel deformation measures, surrogate signals from external markers are used to track internal motion and generate gated PET images. The objective of our work is to develop a group-level parcellation framework for the lungs to guide the placement of markers depending on the location of the internal target region. We present a data-driven framework based on higher-order singular value decomposition (HOSVD) of deformation tensors that enables identification of synchronous areas inside the torso and on the skin surface. Four-dimensional (4-D) magnetic resonance (MR) imaging based on a specialized radial pulse sequence with a one-dimensional slice-projection navigator was used for motion capture under free-breathing conditions. The deformation tensors were computed by nonrigidly registering the gated MR images. Group-level motion signatures obtained via HOSVD were used to cluster the voxels both inside the volume and on the surface. To characterize the parcellation result, we computed correlation measures across the different regions of interest (ROIs). To assess the robustness of the parcellation technique, leave-one-out cross-validation was performed over the subject cohort, and the dependence of the result on varying numbers of gates and singular value thresholds was examined. Overall, the parcellation results were largely consistent across these test cases with Jaccard indices reflecting high degrees of overlap. Finally, a PET simulation study was performed which showed that, depending on the location of the lesion, the selection of a synchronous ROI may lead to noticeable gains in the recovery coefficient. Accurate quantitative interpretation of PET images is important for lung cancer management. Therefore, a guided motion monitoring approach is of utmost importance in the context of pulmonary PET imaging.

Gating Approaches in Cardiac PET Imaging

Cardiac PET provides high sensitivity and high negative predictive value in the diagnosis of coronary artery disease and cardiomyopathies. Cardiac, respiratory as well as bulk patient motion have detrimental effects on thoracic PET imaging, in particular on cardiovascular PET imaging where the motion can affect the PET images quantitatively as well as qualitatively. Gating can ameliorate the unfavorable impact of motion additionally enabling evaluation of left ventricular systolic function. In this article, the authors review the recent advances in gating approaches and highlight the advances in data-driven approaches, which hold promise in motion detection without the need for complex hardware setup.

Quantitative comparison of pre-therapy 99mTc-macroaggregated albumin SPECT/CT and post-therapy PET/MR studies of patients who have received intra-arterial radioembolization therapy with 90Y microspheres

OBJECTIVE

The aim of our study was to compare yttrium -90 (⁹⁰Y) dosimetry obtained from pre-therapy ^99mTc-macroaggregated albumin (MAA) SPECT/CT versus post-therapy PET/MRI imaging among patients with primary or metastatic hepatic tumors.

MATERIALS AND METHODS

Prior to ⁹⁰Y radioembolization (RE), 32 patients underwent a scan using MAA mimicking ⁹⁰Y distribution. After RE with ⁹⁰Y microspheres, the patients were imaged on a PET/MRI system. Reconstructed images were transferred to a common platform and used to calculate ⁹⁰Y dosimetry. The Passing-Bablok regression scatter diagram and the Bland and Altman method were used to analyze the difference between dosimetry values.

RESULTS

For MAA and PET/MRI modalities, the mean liver doses for all 32 subjects were 43.0 ± 20.9 Gy and 46.5 ± 22.7 Gy, respectively, with a mean difference of 3.4 ± 6.2 Gy. The repeatibility coefficient was 12.1 (27.0% of the mean). The Spearman rank correlation coefficient was high (ρ = 0.92). Although, there was a substantial difference in the maximum doses to the liver between the modalities, the mean liver doses were relatively close, with a difference of 24.0% or less.

CONCLUSIONS

The two main contributors to the difference between dosimetry calculations using MAA versus ⁹⁰Y PET/MRI can be attributed to the changes in catheter positioning as well as the liver ROIs used for the calculations. In spite of these differences, our results demonstrate that the dosimetry values obtained from pre-therapy MAA SPECT/CT scans and PET/MRI post-therapy ⁹⁰Y studies were not significantly different.

Super-resolution T2-weighted 4D MRI for image guided radiotherapy

BACKGROUND AND PURPOSE

The superior soft-tissue contrast of 4D-T2w MRI motivates its use for delineation in radiotherapy treatment planning. We address current limitations of slice-selective implementations, including thick slices and artefacts originating from data incompleteness and variable breathing.

MATERIALS AND METHODS

A method was developed to calculate midposition and 4D-T2w images of the whole thorax from continuously acquired axial and sagittal 2D-T2w MRI (1.5 × 1.5 × 5.0 mm3). The method employed image-derived respiratory surrogates, deformable image registration and super-resolution reconstruction. Volunteer imaging and a respiratory motion phantom were used for validation. The minimum number of dynamic acquisitions needed to calculate a representative midposition image was investigated by retrospectively subsampling the data (10-30 dynamic acquisitions).

RESULTS

Super-resolution 4D-T2w MRI (1.0 × 1.0 × 1.0 mm3, 8 respiratory phases) did not suffer from data incompleteness and exhibited reduced stitching artefacts compared to sorted multi-slice MRI. Experiments using a respiratory motion phantom and colour-intensity projection images demonstrated a minor underestimation of the motion range. Midposition diaphragm differences in retrospectively subsampled acquisitions were <1.1 mm compared to the full dataset. 10 dynamic acquisitions were found sufficient to generate midposition MRI.

CONCLUSIONS

A motion-modelling and super-resolution method was developed to calculate high quality 4D/midposition T2w MRI from orthogonal 2D-T2w MRI.

Hybrid positron emission tomography-magnetic resonance of the heart: current state of the art and future applications

Hybrid positron emission tomography-magnetic resonance (PET-MR) imaging is a novel imaging modality with emerging applications for cardiovascular disease. PET-MR aims to combine the high-spatial resolution morphological and functional assessment afforded by magnetic resonance imaging (MRI) with the ability of positron emission tomography (PET) for quantification of metabolism, perfusion, and inflammation. The fusion of these two modalities into a single imaging platform not only represents an opportunity to acquire complementary information from a single scan, but also allows motion correction for PET with reduction in ionising radiation. This article presents a brief overview of PET-MR technology followed by a review of the published literature on the clinical cardio-vascular applications of PET and MRI performed separately and with hybrid PET-MR.

Abdominal DCE-MRI reconstruction with deformable motion correction for liver perfusion quantification

PURPOSE

Abdominal dynamic contrast-enhanced (DCE) MRI suffers from motion-induced artifacts that can blur images and distort contrast-agent uptake curves. For liver perfusion analysis, image reconstruction with rigid-body motion correction (RMC) can restore distorted portal-venous input functions (PVIF) to higher peak amplitudes. However, RMC cannot correct for liver deformation during breathing. We present a reconstruction algorithm with deformable motion correction (DMC) that enables correction of breathing-induced deformation in the whole abdomen.

METHODS

Raw data from a golden-angle stack-of-stars gradient-echo sequence were collected for 54 DCE-MRI examinations of 31 patients. For each examination, a respiratory motion signal was extracted from the data and used to reconstruct 21 breathing states from inhale to exhale. The states were aligned with deformable image registration to the end-exhale state. Resulting deformation fields were used to correct back-projection images before reconstruction with view sharing. Images with DMC were compared to uncorrected images and images with RMC.

RESULTS

DMC significantly increased the PVIF peak amplitude compared to uncorrected images (P << 0.01, mean increase: 8%) but not compared to RMC. The increased PVIF peak amplitude significantly decreased estimated portal-venous perfusion in the liver (P << 0.01, mean decrease: 8 ml/(100 ml·min)). DMC also removed artifacts in perfusion maps at the liver edge and reduced blurring of liver tumors for some patients.

CONCLUSIONS

DCE-MRI reconstruction with DMC can restore motion-distorted uptake curves in the abdomen and remove motion artifacts from reconstructed images and parameter maps but does not significantly improve perfusion quantification in the liver compared to RMC.

Comparison of PET/CT and PET/MR imaging and dosimetry of yttrium-90 (90Y) in patients with unresectable hepatic tumors who have received intra-arterial radioembolization therapy with 90Y microspheres

BACKGROUND

The aim of our study was to compare 90Y dosimetry obtained from PET/MRI versus PET/CT post-therapy imaging among patients with primary or metastatic hepatic tumors. First, a water-filled Jaszczak phantom containing fillable sphere with 90Y-chloride was acquired on both the PET/CT and PET/MRI systems, in order to check the cross-calibration of the modalities. Following selective internal radiation therapy (SIRT) with 90Y microspheres, 32 patients were imaged on a PET/CT system, immediately followed by a PET/MRI study. Reconstructed images were transferred to a common platform and used to calculate 90Y dosimetry. A Passing-Bablok regression scatter diagram and the Bland and Altman method were used to analyze the difference between the dosimetry values.

RESULTS

The phantom study showed that both modalities were calibrated with less than 1% error. The mean liver doses for the 32 subjects calculated from PET/CT and PET/MRI were 51.6 ± 24.7 Gy and 46.5 ± 22.7 Gy, respectively, with a mean difference of 5.1 ± 5.0 Gy. The repeatability coefficient was 9.0 (18.5% of the mean). The Spearman rank correlation coefficient was very high, ρ = 0.97. Although the maximum dose to the liver can be significantly different (up to 40%), mean liver doses from each modalities were relatively close, with a difference of 18.5% or less.

CONCLUSIONS

The two main contributors to the difference in 90Y dosimetry calculations using PET/CT versus PET/MRI can be attributed to the differences in regions of interest (ROIs) and differences attributed to attenuation correction. Due to the superior soft-tissue contrast of MRI, liver contours are usually better seen than in CT images. However, PET/CT provides better quantification of PET images, due to better attenuation correction. In spite of these differences, our results demonstrate that the dosimetry values obtained from PET/MRI and PET/CT in post-therapy 90Y studies were similar.

4D dose calculation for pencil beam scanning proton therapy of pancreatic cancer using repeated 4DMRI datasets

4D magnetic resonance imaging (4DMRI) has a high potential for pancreatic cancer treatments using proton therapy, by providing time-resolved volumetric images with a high soft-tissue contrast without exposing the patient to any additional imaging dose. In this study, we aim to show the feasibility of 4D treatment planning for pencil beam scanning (PBS) proton therapy of pancreatic cancer, based on five repeated 4DMRI datasets and 4D dose calculations (4DDC) for one pancreatic cancer patient. To investigate the dosimetric impacts of organ motion, deformation vector fields were extracted from 4DMRI, which were then used to warp a static CT of the patient, so as to generate synthetic 4DCT (4DCT-MRI). CTV motion amplitudes  <15 mm were observed for this patient. The results from 4DDC show pronounced interplay effects in the CTV with dose homogeneity d5/d95 and dose coverage v95 being 1.14 and 91%, respectively, after a single fraction of the treatment. An averaging effect was further observed when increasing the number of fractions. Motion effects can become less dominant and dose homogeneity d5/d95  =  1.03 and dose coverage v95  =  [Formula: see text] within the CTV can be achieved after 28 fractions. The observed inter-fractional organ and tumor motion variations underline the importance of 4D imaging before and during PBS proton therapy.

Improved Detection of Small Pulmonary Nodules Through Simultaneous MR/PET Imaging

Magnetic resonance (MR)/PET scanners provide an imaging platform that enables simultaneous acquisition of MR and PET data in perfect spatial and temporal registration. This feature allows improving image quality for the MR and PET images obtained during the course of an examination. In this work the authors demonstrate the use of prospective MR-based motion tracking information for removing motion blur in MR/PET images of small pulmonary nodules. The theoretical basis for the algorithms is presented alongside clinical examples of its use.

CAPTURE: Consistently Acquired Projections for Tuned and Robust Estimation: A Self-Navigated Respiratory Motion Correction Approach

OBJECTIVES

In this study, we present a fully automated and robust self-navigated approach to obtain 4-dimensional (4-D) motion-resolved images during free breathing.

MATERIALS AND METHODS

The proposed method, Consistently Acquired Projections for Tuned and Robust Estimation (CAPTURE), is a variant of the stack-of-stars gradient-echo sequence. A 1-D navigator was consistently acquired at a fixed azimuthal angle for all stacks of spokes to reduce nonphysiological signal contamination due to system imperfections. The resulting projections were then "tuned" using complex phase rotation to adapt to scan-to-scan variations, followed by the detection of the respiratory curve. Four-dimensional motion-corrected and uncorrected images were then reconstructed via respiratory and temporal binning, respectively.This Health Insurance Portability and Accountability Act-compliant study was performed with Institutional Review Board approval. A phantom experiment was performed using a custom-made deformable motion phantom with an adjustable frequency and amplitude. For in vivo experiments, 10 healthy participants and 12 liver tumor patients provided informed consent and were imaged with the CAPTURE sequence.Two radiologists, blinded to which images were motion-corrected and which were not, independently reviewed the images and scored the image quality using a 5-point Likert scale.

RESULTS

In the respiratory motion phantom experiment, CAPTURE reversed the effects of motion blurring and restored edge sharpness from 36% to 78% of that observed in the images from the static scan.Despite large intra- and intersubject variability in respiration patterns, CAPTURE successfully detected the respiratory motion signal in all participants and significantly improved the image quality according to the subjective radiological assessments of 2 raters (P < 0.05 for both raters) with a 1 to 2-point improvement in the median Likert scores across the whole set of participants. Small lesions (<1 cm in size) which might otherwise be missed on uncorrected images because of motion blurring were more clearly depicted on the CAPTURE images.

CONCLUSIONS

CAPTURE provides a robust and fully automated solution for obtaining 4-D motion-resolved images in a free-breathing setting. With its unique tuning feature, CAPTURE can adapt to large intersubject and interscan variations. CAPTURE also enables better lesion delineation because of improved image sharpness, thereby increasing the visibility of small lesions.

Magnetic Resonance-based Motion Correction for Quantitative PET in Simultaneous PET-MR Imaging

Motion degrades image quality and quantitation of PET images, and is an obstacle to quantitative PET imaging. Simultaneous PET-MR offers a tool that can be used for correcting the motion in PET images by using anatomic information from MR imaging acquired concurrently. Motion correction can be performed by transforming a set of reconstructed PET images into the same frame or by incorporating the transformation into the system model and reconstructing the motion-corrected image. Several phantom and patient studies have validated that MR-based motion correction strategies have great promise for quantitative PET imaging in simultaneous PET-MR.

Concurrent Respiratory Motion Correction of Abdominal PET and Dynamic Contrast-Enhanced-MRI Using a Compressed Sensing Approach

We present an approach for concurrent reconstruction of respiratory motion-compensated abdominal dynamic contrast-enhanced (DCE)-MRI and PET data in an integrated PET/MR scanner. The MR and PET reconstructions share the same motion vector fields derived from radial MR data; the approach is robust to changes in respiratory pattern and does not increase the total acquisition time. Methods: PET and DCE-MRI data of 12 oncologic patients were simultaneously acquired for 6 min on an integrated PET/MR system after administration of ¹⁸F-FDG and gadoterate meglumine. Golden-angle radial MR data were continuously acquired simultaneously with PET data and sorted into multiple motion phases on the basis of a respiratory signal derived directly from the radial MR data. The resulting multidimensional dataset was reconstructed using a compressed sensing approach that exploits sparsity among respiratory phases. Motion vector fields obtained using the full 6-min (MC_6-min) and only the last 1 min (MC_1-min) of data were incorporated into the PET reconstruction to obtain motion-corrected PET images and in an MR iterative reconstruction algorithm to produce a series of motion-corrected DCE-MR images (moco_GRASP). The motion-correction methods (MC_6-min and MC_1-min) were evaluated by qualitative analysis of the MR images and quantitative analysis of SUV_max and SUV_mean, contrast, signal-to-noise ratio (SNR), and lesion volume in the PET images. Results: Motion-corrected MC_6-min PET images demonstrated 30%, 23%, 34%, and 18% increases in average SUV_max, SUV_mean, contrast, and SNR and an average 40% reduction in lesion volume with respect to the non-motion-corrected PET images. The changes in these figures of merit were smaller but still substantial for the MC_1-min protocol: 19%, 10%, 15%, and 9% increases in average SUV_max, SUV_mean, contrast, and SNR; and a 28% reduction in lesion volume. Moco_GRASP images were deemed of acceptable or better diagnostic image quality with respect to conventional breath-hold Cartesian volumetric interpolated breath-hold examination acquisitions. Conclusion: We presented a method that allows the simultaneous acquisition of respiratory motion-corrected diagnostic quality DCE-MRI and quantitatively accurate PET data in an integrated PET/MR scanner with negligible prolongation in acquisition time compared with routine PET/DCE-MRI protocols.

An MR-Based Model for Cardio-Respiratory Motion Compensation of Overlays in X-Ray Fluoroscopy

In X-ray fluoroscopy, static overlays are used to visualize soft tissue. We propose a system for cardiac and respiratory motion compensation of these overlays. It consists of a 3-D motion model created from real-time magnetic resonance (MR) imaging. Multiple sagittal slices are acquired and retrospectively stacked to consistent 3-D volumes. Slice stacking considers cardiac information derived from the ECG and respiratory information extracted from the images. Additionally, temporal smoothness of the stacking is enhanced. Motion is estimated from the MR volumes using deformable 3-D/3-D registration. The motion model itself is a linear direct correspondence model using the same surrogate signals as slice stacking. In X-ray fluoroscopy, only the surrogate signals need to be extracted to apply the motion model and animate the overlay in real time. For evaluation, points are manually annotated in oblique MR slices and in contrast-enhanced X-ray images. The 2-D Euclidean distance of these points is reduced from 3.85 to 2.75 mm in MR and from 3.0 to 1.8 mm in X-ray compared with the static baseline. Furthermore, the motion-compensated overlays are shown qualitatively as images and videos.

Advances in PET/MR instrumentation and image reconstruction

The combination of positron emission tomography (PET) and MRI has attracted the attention of researchers in the past approximately 20 years in small-animal imaging and more recently in clinical research. The combination of PET/MRI allows researchers to explore clinical and research questions in a wide number of fields, some of which are briefly mentioned here. An important number of groups have developed different concepts to tackle the problems that PET instrumentation poses to the exposition of electromagnetic fields. We have described most of these research developments in preclinical and clinical experiments, including the few commercial scanners available. From the software perspective, an important number of algorithms have been developed to address the attenuation correction issue and to exploit the possibility that MRI provides for motion correction and quantitative image reconstruction, especially parametric modelling of radiopharmaceutical kinetics. In this work, we give an overview of some exemplar applications of simultaneous PET/MRI, together with technological hardware and software developments.