4D-MRI Technique Research Articles

Free-breathing time-resolved 4D MRI with improved T1-weighting contrast.

This work proposes MP-Grasp4D (magnetization-prepared golden-angle radial sparse parallel 4D) MRI, a free-breathing, inversion recovery (IR)-prepared, time-resolved 4D MRI technique with improved T1-weighted contrast. MP-Grasp4D MRI acquisition incorporates IR preparation into a radial gradient echo sequence. MP-Grasp4D employs a golden-angle navi-stack-of-stars sampling scheme, where imaging data of rotating radial stacks and navigator stacks (acquired at a consistent rotation angle) are alternately acquired. The navigator stacks are used to estimate a temporal basis for low-rank subspace-constrained reconstruction. This allows for the simultaneous capture of both IR-induced contrast changes and respiratory motion. One temporal frame of the imaging volume in MP-Grasp4D MRI is reconstructed from a single stack and an adjacent navigator stack on average, resulting in a nominal temporal resolution of 0.16 seconds per volume. Images corresponding to the optimal inversion time (TI) can be retrospectively selected for providing the best image contrast. Reader studies were conducted to assess the performance of MP-Grasp4D MRI in liver imaging across 30 subjects in comparison with standard Grasp4D MRI without IR preparation. MP-Grasp4D MRI received significantly higher scores (P < 0.05) than Grasp4D in all assessment categories. There was a moderate to almost perfect agreement (kappa coefficient from 0.42 to 0.9) between the two readers for image quality assessment. When the scan time is reduced, MP-Grasp4D MRI preserves image contrast and quality, demonstrating additional acceleration capability. MP-Grasp4D MRI improves T1-weighted contrast for free-breathing time-resolved 4D MRI and eliminates the need for explicit motion compensation. This method is expected to be valuable in different MRI applications such as MR-guided radiotherapy.

1248: A systematic review of 4D MRI techniques for abdominal radiotherapy treatment planning

1248: A systematic review of 4D MRI techniques for abdominal radiotherapy treatment planning

4D-MRI assisted stereotactic body radiation therapy for unresectable colorectal cancer liver metastases

This study evaluated the feasibilities and outcomes following four-dimensional magnetic resonance imaging (4D-MRI) assisted stereotactic body radiation therapy (SBRT) for unresectable colorectal liver metastases (CRLMs). From March 2018 to January 2022, we identified 76 unresectable CRLMs patients with 123 lesions who received 4D-MRI guided SBRT in our institution. 4D-MRI simulation with or without abdominal compression was conducted for all patients. The prescription dose was 50–65 Gy in 5–12 fractions. The image quality of computed tomography (CT) and MRI were compared using the Clarity Score. Clinical outcomes and toxicity profiles were evaluated. 4D-MRI improved the image quality compared with CT images (mean Clarity Score: 1.67 vs 2.88, P < 0.001). The abdominal compression reduced motions in cranial–caudal direction (P = 0.03) with two phase T2 weighted images assessing tumor motion. The median follow-up time was 12.5 months. For 98 lesions assessed for best response, the complete response, partial response and stable disease rate were 57.1 %, 30.6 % and 12.2 %, respectively. The local control (LC) rate at 1 year was 97.3 %. 46.1 % of patients experienced grade 1–2 toxicities and only 2.6 % patients experienced grade 3 hematologic toxicities. The 4D-MRI technique allowed accurate target delineation and motion tracking in unresectable CRLMs patients. Favorable LC rate and mild toxicities were achieved. This study provided evidence for using 4D-MRI assisted SBRT as an alternative treatment in unresectable CRLMs.

4D-MRI Guided Stereotactic Body Radiation Therapy for Unresectable Colorectal Liver Metastases

This study evaluated the feasibilities and outcomes following four-dimensional magnetic resonance imaging (4D-MRI) guided stereotactic body radiation therapy (SBRT) for unresectable colorectal liver metastases (CRLM). From March 2018 to January 2022, we identified 76 unresectable CRLM patients with 123 lesions who received 4D-MRI guided SBRT in our institution. 4D-MRI simulation with or without abdominal compression was conducted for all patients. The prescription dose was 50-65 Gy in 5-12 fractions. The image quality of computed tomography (CT) and MRI were compared using the Clarity Score. Clinical outcomes and toxicity profiles were evaluated. The 4D-MRI significantly improved the image quality compared with CT images (mean Clarity Score: 1.67 vs 2.88, P < 0.001). The abdominal compression significantly reduced motions in cranial-caudal direction (P = 0.03) with 2 phase T2 weighted images assessing tumor motion. The median follow-up time was 12.5 months. For 98 lesions assessed for best response, the complete response, partial response and stable disease rate were 57.1 %, 30.6 % and 12.2 %, respectively. The local control (LC) rate at 2 year was 97.3%. 46.1% of patients experienced grade 1-2 toxicities and only 2.6% patients experienced grade 3 hematologic toxicities. The 4D-MRI technique allowed precise target delineation and motion tracking in unresectable CRLM patients. High LC rate and mild toxicities were achieved. This study provided evidence for using 4D-MRI guided SBRT as an alternative treatment in unresectable CRLM.

Comparative study of sub-second temporal resolution 4D-MRI and 4D-CT for target motion assessment in a phantom model

To develop and investigate the feasibility of sub-second temporal resolution volumetric T1-weighted four-dimensional (4D-) MRI in comparison with 4D-CT for respiratory-correlated motion assessment using an MRI/CT-compatible phantom. Sub-second high temporal resolution (0.5 s) gradient-echo T1-weighted 4D-MRI was developed using a volumetric acquisition scheme with compressed sensing. An MRI/CT-compatible motion phantom (simulated liver tumor) with three sinusoidal movements of amplitudes and two respiratory patterns was introduced and imaged with 4D-MRI and 4D-CT to investigate the geometric accuracy of the target movement. The geometric accuracy, including centroid position, volume, similarity index of dice similarity coefficient (DSC), and Hausdorff distance (HD), was systematically evaluated. Proposed 4D-MRI achieved a similar geometric accuracy compared with 4D-CT regarding the centroid position, volume, and similarity index. The observed position differences of the absolute average centroid were within 0.08 cm in 4D-MRI and 0.03 cm in 4D-CT, less than the 1-pixel resolution for each modality. The observed volume difference in 4D-MRI/4D-CT was within 0.73 cm3 (4.5%)/0.29 cm3 (2.1%) for a large target and 0.06 cm3 (11.3%)/0.04 cm3 (11.6%) for a small target. The observed DSC values for 4D-MRI/4D-CT were at least 0.93/0.95 for the large target and 0.83/0.84 for the small target. The maximum HD values were 0.25 cm/0.31 cm for the large target and 0.21 cm/0.15 cm for the small target. Although 4D-CT potentially exhibit superior numerical accuracy in phantom studies, the proposed high temporal resolution 4D-MRI demonstrates sub-millimetre geometric accuracy comparable to that of 4D-CT. These findings suggest that the 4D-MRI technique is a viable option for characterizing motion and generating phase-dependent internal target volumes within the realm of radiotherapy.

Multi-contrast four-dimensional magnetic resonance imaging (MC-4D-MRI): Development and initial evaluation in liver tumor patients.

To develop a novel multi-contrast four-dimensional magnetic resonance imaging (MC-4D-MRI) technique that expands single image contrast 4D-MRI to a spectrum of native and synthetic image contrasts and to evaluate its feasibility in liver tumor patients. The MC-4D-MRI technique integrates multi-parametric MRI fusion, 4D-MRI, and deformable image registration (DIR) techniques. The fusion technique consists of native MRI as input, image pre-processing, fusion algorithm, adaptation, and fused multi-contrast MRI as output. Four-dimensional deformation vector fields (4D-DVF) were generated from an original T2/T1-w 4D-MRI by deforming end-of-inhalation (EOI) to nine other phase volumes via DIR. The 4D-DVF were applied to multi-contrast MRI to generate a spectrum of 4D-MRI in different image contrasts. The MC-4D-MRI technique was evaluated in five liver tumor patients on tumor contrast-to-noise ratio (CNR), internal target volume (ITV) contouring consistency, diaphragm motion range, and tumor motion trajectory; and in digital anthropomorphic phantoms on 4D-DIR introduced errors in tumor motion range, centroid location, extent, and volume. MC-4D-MRI consisting of 4D-MRIs in native image contrasts (T1-w, T2-w, and T2/T1-w) and synthetic image contrasts, such as tumor-enhanced contrast (TEC) were generated in five liver tumor patients. Patient tumor CNR increased from 2.6±1.8 in the T2/T1-w MRI, to -4.4±2.4, 6.6±3.0, and 9.6±3.9 in the T1-w, T2-w, and TEC MRI, respectively. Patient ITV inter-observer mean Dice similarity coefficient (mDSC) increased from 0.65±0.10 in the original T2/T1-w 4D-MRI, to 0.76±0.14, 0.77±0.12, and 0.86±0.05 in the T1-w, T2-w, and TEC 4D-MRI, respectively. Patient diaphragm motion range absolute differences between the three new 4D-MRIs and original T2/T1-w 4D-MRI were 1.2±1.3, 0.3±0.7, and 0.5±0.5mm, respectively. Patient tumor displacement phase-averaged absolute differences between the three 4D-MRIs and the original 4D-MRI were 0.72±0.33, 0.62±0.54, and 0.74±0.43mm in the superior-inferior (SI) direction, and 0.59±0.36, 0.51±0.30, and 0.50±0.24mm in the anterior-posterior (AP) direction, respectively. In the digital phantoms, phase-averaged absolute tumor centroid shift caused by the 4D-DIR were at or below 0.5mm in SI, AP, and left-right (LR) directions. We developed an MC-4D-MRI technique capable of expanding single image contrast 4D-MRI along a new dimension of image contrast. Initial evaluations in liver tumor patients showed enhancements in image contrast variety, tumor contrast, and ITV contouring consistencies using MC-4D-MRI. The technique might offer new perspectives on the image contrast of MRI and 4D-MRI in MR-guided radiotherapy.

Gastrointestinal 4D MRI with respiratory motion correction.

Gastrointestinal motion patterns such as peristalsis and segmental contractions can alter the shape and position of the stomach and intestines with respect to other irradiated organs during radiation therapy. Unfortunately, these deformations are concealed by conventional four-dimensional (4D)-MRI techniques, which were developed to visualize respiratory motion by binning acquired data into respiratory motion states without considering the phases of GI motion. We present a method to reconstruct breathing-compensated images showing the phases of periodic gastric motion and study the effect of this motion on regional anatomical structures. Sixty-seven DCE-MRI examinations were performed on patients undergoing MRI simulation for hepatocellular carcinoma using a golden-angle stack-of-stars sequence that collected 2000 radial spokes over 5min. The collected data were reconstructed using a method with integrated respiratory motion correction into a time series of 3D image volumes without visible breathing motion. From this series, a gastric motion signal was extracted by temporal filtering of time-intensity curves in the stomach. Using this motion signal, breathing-corrected back-projection images were sorted according to the gastric phase and reconstructed into 21 gastric motion state images showing the phases of gastric motion. Reconstructed image volumes showed gastric motion states clearly with no visible breathing motion or related artifacts. The mean frequency of the gastric motion signal was 3cycles/min with a standard deviation of 0.27cycles/min. Periodic gastrointestinal motion can be visualized without confounding respiratory motion using the presented GI 4D MRI technique. GI 4D MRIs may help define internal target volumes for treatment planning, aid in planning organ at risk volume definition, or support motion model development for gastrointestinal motion tracking algorithms for real-time MR-guided radiation therapy.

Respiratory-Correlated (RC) vs. Time-Resolved (TR) Four-Dimensional Magnetic Resonance Imaging (4DMRI) for Radiotherapy of Thoracic and Abdominal Cancer.

Recent technological and clinical advancements of both respiratory-correlated (RC) and time-resolved (TR) four-dimensional magnetic resonance imaging (4DMRI) techniques are reviewed in light of tumor/organ motion simulation, monitoring, and assessment in radiotherapy. For radiotherapy of thoracic and abdominal cancer, respiratory-induced tumor motion, and motion variation due to breathing irregularities are the major uncertainties in treatment. RC-4DMRI is developed to assess tumor motion for treatment planning, whereas TR-4DMRI is developed to assess both motion and motion variation for treatment planning, delivery and assessment. RC-4DMRI is reconstructed to provide one-breathing-cycle motion, similar to 4D computed tomography (4DCT), the current clinical standard, but with higher soft-tissue contrast, no ionizing radiation, and less binning artifacts due to the use of an internal respiratory surrogate. Recent studies have shown that its spatial resolution has reached or exceeded that of 4DCT and scanning time becomes clinically acceptable. TR-4DMRI is recently developed with an adequate spatiotemporal resolution to assess tumor motion and motion variations for treatment simulation, delivery and assessment. The super-resolution approach is most promising since it can image any organ/body motion, whereas RC-4D MRI are limited to resolve only respiration-induced motion and some TR-4DMRI approaches may more or less depend on RC-4DMRI. TR-4DMRI provides multi-breath motion data that are useful not only in MR-guided radiotherapy but also for building a patient-specific motion model to guide radiotherapy treatment using an non-MR-equipped linear accelerator. Based on 4DMRI motion data, motion-corrected dynamic contrast imaging and diffusion-weighted imaging have also been reported, aiming to facilitate tumor delineation for more accurate radiotherapy targeting. Both RC- and TR-4DMRI have been evaluated for potential clinical applications, such as delineation of tumor volumes, where sufficiently high spatial resolution and large field-of-view are required. The 4DMRI techniques are promising to play a role in motion assessment in radiotherapy treatment planning, delivery, assessment, and adaptation.

Gating Window Optimization for Respiratory Phase-Gated Radiotherapy Using the Ultra-Fast Volumetric 4D-MRI Generated 3D Probability Distribution Map

Respiratory phase-gated RT is clinically popular in reducing radiation dose to the surrounding normal tissues by managing the respiratory motion. Although the optimal setting is controversial, the gating windows are usually set at 30-50%, 30-60%, 40-60% or 30-70% (0% = inhale, 50% = exhale), and is usually based simply on the variation of the SI translational surrogate. In this study, we aim to optimize the gating window based on the volume of overlapping VOI with a probability of 90% (V90,i, i = respiratory phase) through an ultra-fast volumetric 4D-MRI technique. Free breathing volumetric abdominal MR images were acquired on 9 healthy subjects (34 ± 6 yrs) using a 1.5T MR-simulator (3D gradient echo sequence with 1.63 vol/sec, 56 axial slice/vol, voxel size = 2.7 x 2.7 x 4 mm3, 144 dynamics). Binary mask of each dynamic volumetric of the left kidney (LK) and right kidney (RK) was segmented. The images were binned into 10 respiratory phases based on liver position. For all phases, the binary mask of each volumetric time frame at phase i were summated prior to the normalization to the total number of dynamic volumes at phase i (resulting in a total of 10 probability distribution maps). For each phase i, V90, i and the SI translational range (RSI, i) of the liver position were evaluated. For each subject, the optimal gating window phase (Pw) and the gating window durations (Dw) was determined with two approaches based on (1) RSI (Pw = phases with minimum RSI, Dw = neighboring phases of Pw with RSI< 2mm); (2) dice similarity coefficient (DSC) (Pw = phase of maximum V90 (V90,max), DW = neighboring phases of Pw with DSC > 0.9). The mean RSI of the phase of V90,max for all VOIs, and the mean minimum RSI (RSI,min) were also reported. Comparison between RSI and DSC approaches was performed using signed rank test. Comparing between approaches, similar Pw at exhalation was obtained using RSI (50±11%) and DSC (LK:46±24%, p=1; RK:54±19%, p=0.69). The mean V90,max (LK=156±34mm3, RK=149±26mm3) was close to the mean organ volume (LK=160±34mm3; RK=154±27mm3). The mean RSI of the phase at V90,max (LK=2.48±1.92mm, p=0.38; RK=2.41±1.91mm, p=1) was similar to RSI,min (2.06±1.99mm). The choice of gating window using DSC was preferable because a highly repeatable volumetric coverage was obtained, while such volume information cannot be fully revealed by the RSI approach. A smaller Dw was obtained using RSI (14±5%) then DSC (LK:42±23%, p=0.016; RK:60±31%, p=0.016), though insignificant, might suggest a less effective gated-RT delivery based simply on SI displacement. Larger Dw was obtained based on DSC approach allows a more effective gated-RT delivery with the overlapping volume taken into consideration.

Probability-based 3D k-space sorting for motion robust 4D-MRI.

Current 4D-MRI techniques are prone to breathing-variation-induced motion artifacts. This study developed a novel method for motion-robust multi-cycle 4D-MRI using probability-based multi-cycle sorting to overcome this deficiency. The main cycles were first extracted from the breathing signal. 3D k-space data were then sorted using a result-driven method for each main cycle. The new method was tested on a 4D-extended cardiac-torso (XCAT) phantom with a patient and an artificially generated breathing curve. For comparison, the k-space data were sorted using conventional phase sorting to generate single-cycle 4D-MRI images. Signal-to-noise ratio (SNR) of tumor and liver, tumor volume consistency, and average intensity projection (AIP) accuracy were compared between the two methods. The original phantom images were used as references for the evaluation. The new method showed improved tumor-to-liver SNR and tumor volume consistency as compared to 3D k-space phase sorting in both the simulated artificial and real patient breathing signals. For the artificial breathing cycles, the average tumor-to-liver SNR and standard deviation (SD) of tumor volume were 2.53 and 3.80% for cycle 1, 2.24 and 6.16% for cycle 2 of probability-based sorting as compared to 1.47 and 21.83% obtained using the phase sorting method; for the patient breathing curve, values of 1.99 and 2.71%, 1.97 and 3.29%, 1.88 and 4.16% were observed for cycle 1, cycle 2 and cycle 3 of probability-based sorting, versus 1.44 and 7.20% for phase sorting method. Furthermore, the AIP accuracy was improved in the probability-based sorting approach when compared to phase sorting, with the average intensity difference per voxel reduced from 0.39 to 0.15 for the artificial curve, and from 0.46 to 0.21 for the patient curve. We demonstrated the feasibility of probability-based 3D k-space sorting for motion-robust multi-cycle 4D-MRI reconstruction with breathing variation induced motion artifact reduction compared with conventional 2D image sorting and 3D phase sorting methods. This new technique can potentially improve the accuracy of radiation treatment guidance for mobile targets.

Respiratory motion-resolved, self-gated 4D-MRI using rotating cartesian k-space (ROCK).

To propose and validate a respiratory motion resolved, self-gated (SG) 4D-MRI technique to assess patient-specific breathing motion of abdominal organs for radiation treatment planning. The proposed 4D-MRI technique was based on the balanced steady-state free-precession (bSSFP) technique and 3D k-space encoding. A novel rotating cartesian k-space (ROCK) reordering method was designed which incorporates repeatedly sampled k-space centerline as the SG motion surrogate and allows for retrospective k-space data binning into different respiratory positions based on the amplitude of the surrogate. The multiple respiratory-resolved 3D k-space data were subsequently reconstructed using a joint parallel imaging and compressed sensing method with spatial and temporal regularization. The proposed 4D-MRI technique was validated using a custom-made dynamic motion phantom and was tested in six healthy volunteers, in whom quantitative diaphragm and kidney motion measurements based on 4D-MRI images were compared with those based on 2D-CINE images. The 5-minute 4D-MRI scan offers high-quality volumetric images in 1.2 × 1.2 × 1.6 mm3 and eight respiratory positions, with good soft-tissue contrast. In phantom experiments with triangular motion waveform, the motion amplitude measurements based on 4D-MRI were 11.89% smaller than the ground truth, whereas a -12.5% difference was expected due to data binning effects. In healthy volunteers, the difference between the measurements based on 4D-MRI and the ones based on 2D-CINE were 6.2 ± 4.5% for the diaphragm, 8.2 ± 4.9% and 8.9 ± 5.1% for the right and left kidney. The proposed 4D-MRI technique could provide high-resolution, high-quality, respiratory motion-resolved 4D images with good soft-tissue contrast and are free of the "stitching" artifacts usually seen on 4D-CT and 4D-MRI based on resorting 2D-CINE. It could be used to visualize and quantify abdominal organ motion for MRI-based radiation treatment planning.

TH‐EF‐BRA‐11: Feasibility of Super‐Resolution Time‐Resolved 4DMRI for Multi‐Breath Volumetric Motion Simulation in Radiotherapy Planning

Purpose:To develop a novel super‐resolution time‐resolved 4DMRI technique to evaluate multi‐breath, irregular and complex organ motion without respiratory surrogate for radiotherapy planning.Methods:The super‐resolution time‐resolved (TR) 4DMRI approach combines a series of low‐resolution 3D cine MRI images acquired during free breathing (FB) with a high‐resolution breath‐hold (BH) 3DMRI via deformable image registration (DIR). Five volunteers participated in the study under an IRB‐approved protocol. The 3D cine images with voxel size of 5×5×5 mm3 at two volumes per second (2Hz) were acquired coronally using a T1 fast field echo sequence, half‐scan (0.8) acceleration, and SENSE (3) parallel imaging. Phase‐encoding was set in the lateral direction to minimize motion artifacts. The BH image with voxel size of 2×2×2 mm3 was acquired using the same sequence within 10 seconds. A demons‐based DIR program was employed to produce super‐resolution 2Hz 4DMRI. Registration quality was visually assessed using difference images between TR 4DMRI and 3D cine and quantitatively assessed using average voxel correlation. The fidelity of the 3D cine images was assessed using a gel phantom and a 1D motion platform by comparing mobile and static images.Results:Owing to voxel intensity similarity using the same MRI scanning sequence, accurate DIR between FB and BH images is achieved. The voxel correlations between 3D cine and TR 4DMRI are greater than 0.92 in all cases and the difference images illustrate minimal residual error with little systematic patterns. The 3D cine images of the mobile gel phantom preserve object geometry with minimal scanning artifacts.Conclusion:The super‐resolution time‐resolved 4DMRI technique has been achieved via DIR, providing a potential solution for multi‐breath motion assessment. Accurate DIR mapping has been achieved to map high‐resolution BH images to low‐resolution FB images, producing 2Hz volumetric high‐resolution 4DMRI. Further validation and improvement are still required prior to clinical applications.This study is in part supported by the NIH (U54CA137788/U54CA132378).

MO-FG-CAMPUS-JeP2-01: 4D-MRI with 3D Radial Sampling and Self-Gating-Based K-Space Sorting: Image Quality Improvement by Slab-Selective Excitation

Purpose: A recent 4D MRI technique based on 3D radial sampling and self-gating-based K-space sorting has shown promising results in characterizing respiratory motion. However due to continuous acquisition and potentially drastic k-space undersampling resultant images could suffer from low blood-to-tissue contrast and streaking artifacts. In this study 3D radial sampling with slab-selective excitation (SS) was proposed in attempt to enhance blood-to-tissue contrast by exploiting the in-flow effect and to suppress the excess signal from the peripheral structures particularly in the superior-inferior direction. The feasibility of improving image quality by using this approach was investigated through a comparison with the previously developed non-selective excitation (NS) approach. Methods: Two excitation approaches SS and NS were compared in 5 cancer patients (1 lung 1 liver 2 pancreas and 1 esophagus) at 3Tesla. Image artifact was assessed in all patients on a 4-point scale (0: poor; 3: excellent). Signal-tonoise ratio (SNR) of the blood vessel (aorta) at the center of field-of-view and its nearby tissue were measured in 3 of the 5 patients (1 liver 2 pancreas) and blood-to-tissue contrast-to-noise ratio (CNR) were then determined. Results: Compared with NS the image quality of SS was visually improved with overall higher signal in all patients (2.6±0.55 vs. 3.4±0.55). SS showed an approximately 2-fold increase of SNR in the blood (aorta: 16.39±1.95 vs. 32.19±7.93) and slight increase in the surrounding tissue (liver/pancreas: 16.91±1.82 vs. 22.31±3.03). As a result the blood-totissue CNR was dramatically higher in the SS method (1.20±1.20 vs. 9.87±6.67). Conclusion: The proposed 3D radial sampling with slabselective excitation allows for reduced image artifact and improved blood SNR and blood-to-tissue CNR. The success of this technique could potentially benefit patients with cancerous tumors that have invaded the surrounding blood vessels where radiation therapy is needed to remove tumor from those regions prior to surgical resection. This work is partially supported by NIH R03CA173273; and CTSI core voucher award.

Four Dimensional Magnetic Resonance Imaging With 3D Radial Sampling and Self-gating Based K-space Sorting: Early Clinical Experience on Pancreatic Cancer Patients

Dynamic magnetic resonance imaging (MRI) has been used to characterize internal organ motion, but real-time acquisition has been limited to 2 dimensions (2D). Methods have been developed to reconstruct four dimensional MRI (4D-MRI) based on rebinning time-stamped 2D images or 2D K-space data. These methods suffer from anisotropic resolution and rebinning artifacts. We developed a 4D-MRI technique using 3D radial sampling and self-gating based K-space sorting (SG-KS-4D-MRI) to overcome these limitations, and applied it to study the target motion of pancreatic cancer patients. Ten patients were imaged during free breathing using 4D-CT, cine 2D-MRI and the SG-KS-4D-MRI method, which is an 8-minute spoiled gradient recalled echo (GRE) sequence with 3D radial-sampling K-space projections and 1D projection based self-gating. Tumor volumes were drawn at the end of exhalation phases in the SG-KS-4D-MRI and 4D-CT for each patient, and then mapped to the other phases using deformable registration with B-spline algorithm. Tumor volumes and motion trajectories were compared. An isotropic resolution of 1.6 mm was achieved in the SG-KS-4D-MRI images, which showed superior soft tissue contrast to 4D-CT and were free of visible rebinning artifacts. The tumor motion trajectory cross-correlations between SG-KS-4D-MRI and cine 2D-MRI directions in superior-inferior (SI), anterior-posterior (AP) and medial-lateral (ML) directions were 0.93 ± 0.03, 0.83 ± 0.10, and 0.74 ± 0.18, respectively. The tumor motion trajectory cross-correlations between the SG-KS-4D-MRI and 4D-CT in SI, AP and ML were 0.91 ± 0.06, 0.72 ± 0.18, and 0.44 ± 0.24, respectively. The average motion amplitude differences (Ddiff) between SG-KS-4D-MRI and cine 2D-MRI in SI, AP, and ML directions were 0.83 ± 0.52, 0.40 ± 0.17, and 0.48 ± 0.24 mm. The Ddiff between SG-KS-4D-MRI and 4D-CT in SI, AP, and ML directions were 1.09 ± 0.43, 0.53 ± 0.24, and 0.47 ± 0.26 mm, respectively. The average standard deviations of the GTV volumes (GTV_σ) calculated from each of the ten breathing phases were 0.8 cm3 (SG-KS-4D-MRI) and 1.0 cm3 (4D-CT), P = .004, indicating a statistically significant difference between the two imaging modalities. A novel 4D-MRI imaging sequence, SG-KS-4D-MRI, was used to evaluate target motion for free breathing pancreatic cancer patients. Pancreatic motion trajectories agreed well with 4D-CT and cine 2D-MRI data in the SI and AP directions, but agreement was limited in the ML direction, where image pixel size and displacements are of similar size. Use of SG-KS-4D-MRI results in improved spatial resolution, voxel isotropy, and gives an average motion pattern over an extended period of time; SG-KS-4D-MRI also results in statistically more consistent tumor volume definition across the ten breathing phases.

T2-weighted four dimensional magnetic resonance imaging with result-driven phase sorting.

T2-weighted MRI provides excellent tumor-to-tissue contrast for target volume delineation in radiation therapy treatment planning. This study aims at developing a novel T2-weighted retrospective four dimensional magnetic resonance imaging (4D-MRI) phase sorting technique for imaging organ/tumor respiratory motion. A 2D fast T2-weighted half-Fourier acquisition single-shot turbo spin-echo MR sequence was used for image acquisition of 4D-MRI, with a frame rate of 2-3 frames/s. Respiratory motion was measured using an external breathing monitoring device. A phase sorting method was developed to sort the images by their corresponding respiratory phases. Besides, a result-driven strategy was applied to effectively utilize redundant images in the case when multiple images were allocated to a bin. This strategy, selecting the image with minimal amplitude error, will generate the most representative 4D-MRI. Since we are using a different image acquisition mode for 4D imaging (the sequential image acquisition scheme) with the conventionally used cine or helical image acquisition scheme, the 4D dataset sufficient condition was not obviously and directly predictable. An important challenge of the proposed technique was to determine the number of repeated scans (NR) required to obtain sufficient phase information at each slice position. To tackle this challenge, the authors first conducted computer simulations using real-time position management respiratory signals of the 29 cancer patients under an IRB-approved retrospective study to derive the relationships between NR and the following factors: number of slices (NS), number of 4D-MRI respiratory bins (NB), and starting phase at image acquisition (P0). To validate the authors' technique, 4D-MRI acquisition and reconstruction were simulated on a 4D digital extended cardiac-torso (XCAT) human phantom using simulation derived parameters. Twelve healthy volunteers were involved in an IRB-approved study to investigate the feasibility of this technique. 4D data acquisition completeness (Cp) increases as NR increases in an inverse-exponential fashion (Cp = 100 - 99 × exp(-0.18 × NR), when NB = 6, fitted using 29 patients' data). The NR required for 4D-MRI reconstruction (defined as achieving 95% completeness, Cp = 95%, NR = NR,95) is proportional to NB (NR,95 ∼ 2.86 × NB, r = 1.0), but independent of NS and P0. Simulated XCAT 4D-MRI showed a clear pattern of respiratory motion. Tumor motion trajectories measured on 4D-MRI were comparable to the average input signal, with a mean relative amplitude error of 2.7% ± 2.9%. Reconstructed 4D-MRI for healthy volunteers illustrated clear respiratory motion on three orthogonal planes, with minimal image artifacts. The artifacts were presumably caused by breathing irregularity and incompleteness of data acquisition (95% acquired only). The mean relative amplitude error between critical structure trajectory and average breathing curve for 12 healthy volunteers is 2.5 ± 0.3 mm in superior-inferior direction. A novel T2-weighted retrospective phase sorting 4D-MRI technique has been developed and successfully applied on digital phantom and healthy volunteers.

SU‐E‐J‐240: Development of a Novel 4D MRI Sequence for Real‐Time Liver Tumor Tracking During Radiotherapy

Purpose:To develop a Novel 4D MRI Technique that is feasible for realtime liver tumor tracking during radiotherapy.Methods:A volunteer underwent an abdominal 2D fast EPI coronal scan on a 3.0T MRI scanner (Siemens Inc., Germany). An optimal set of parameters was determined based on image quality and scan time. A total of 23 slices were scanned to cover the whole liver in the test scan. For each scan position, the 2D images were retrospectively sorted into multiple phases based on breathing signal extracted from the images. Consequently the 2D slices with same phase numbers were stacked to form one 3D image. Multiple phases of 3D images formed the 4D MRI sequence representing one breathing cycle.Results:The optimal set of scan parameters were: TR= 57ms, TE= 19ms, FOV read= 320mm and flip angle= 30°, which resulted in a total scan time of 14s for 200 frames (FMs) per slice and image resolution of (2.5mm,2.5mm,5.0mm) in three directions. Ten phases of 3D images were generated, each of which had 23 slices. Based on our test scan, only 100FMs were necessary for the phase sorting process which may lower the scan time to 7s/100FMs/slice. For example, only 5 slices/35s are necessary for a 4D MRI scan to cover liver tumor size ≤ 2cm leading to the possibility of tumor trajectory tracking every 35s during treatment.Conclusion:The novel 4D MRI technique we developed can reconstruct a 4D liver MRI sequence representing one breathing cycle (7s/ slice) without an external monitor. This technique can potentially be used for real‐time liver tumor tracking during radiotherapy.

TH-CD-207-09: Retrospective 4D-MRI with a Novel Image-Based Surrogate: A Sagittal-Coronal-Diaphragm Point of Intersection (SCD-PoI) Motion Tracking Method

Purpose: The unreliable stability of internal respiratory surrogates and inconvenience of external respiratory surrogates for current retrospective 4D-MRI techniques largely affects the image quality of 4D-MRI. This study aims at developing image-based surrogate, a sagittal-coronal-diaphragm point of intersection (SCD-PoI) motion tracking method for retrospective 4D-MRI reconstruction.Methods/materials: As a pre-estimate of respiratory motion pattern, single-slice sagittal cines (FIESTA) were acquired at a location near the dome of the diaphragm. Subsequently, multi-slice coronal cines (FIESTA) were acquired and used for 4D-MRI reconstruction with phase sorting. Diaphragm motion trajectories were measured from the point of intersection between sagittal MRI cine plane, coronal MRI cine plain and the diaphragm dome surface. This point is defined as sagittal-coronal-diaphragm point of intersection (SCD-PoI). We pre-estimate respiraotyr motion by tracking SCD-PoI on sagitall cine. Then coronal images were then re-binned to different phased bins according to SCD-PoI motion tracking on coronal cine. This 4D-MRI technique was evaluated on a 4D Digital Extended Cardiac-Torso (XCAT) human phantom with a hypothesized moving tumor, six healthy voluneteers and two cancer patients under an IRB-approved study. Region of interest (ROI: tumor for XCAT and patients, dome of left kidney for healthy volunteers) trajectories on 4D-MRI were measured and compared with the reference (input respiratory curve for XCAT and ROI trajectories extracted from reference single-slice MRI cine (FIESTA) for human subjects). Superior-inferior (SI) mean absolute amplitude difference (D) and cross-correlation coefficient (CC) were calculated. Results: 4D-MRI on XCAT demonstrated highly accurate motion information with a low D (1.13mm) and a high CC (0.98) in the SI direction. Minimal artifacts were observed in human participants’ 4D-MRI, and images were adequate to reveal the respiratory motion of organs and tumor (D=1.08±1.03mm; CC=0.96). Conclusion: A novel 4D-MRI technique with image-based respiratory surrogate has been developed and tested on a digital phantom and human subjects.

Four-dimensional MRI using three-dimensional radial sampling with respiratory self-gating to characterize temporal phase-resolved respiratory motion in the abdomen.

To develop a four-dimensional MRI (4D-MRI) technique to characterize the average respiratory tumor motion for abdominal radiotherapy planning. A continuous spoiled gradient echo sequence was implemented with 3D radial trajectory and 1D self-gating for respiratory motion detection. Data were retrospectively sorted into different respiratory phases based on their temporal locations within a respiratory cycle, and each phase was reconstructed by means of a self-calibrating CG-SENSE program. Motion phantom, healthy volunteer and patient studies were performed to validate the respiratory motion detected by the proposed method against that from a 2D real-time protocol. The proposed method successfully visualized the respiratory motion in phantom and human subjects. The 4D-MRI and real-time 2D-MRI yielded comparable superior-inferior (SI) motion amplitudes (intraclass correlation = 0.935) with up-to one pixel mean absolute differences in SI displacements over 10 phases and high cross-correlation between phase-resolved displacements (phantom: 0.985; human: 0.937-0.985). Comparable anterior-posterior and left-right displacements of the tumor or gold fiducial between 4D and real-time 2D-MRI were also observed in the two patients, and the hysteresis effect was shown in their 3D trajectories. We demonstrated the feasibility of the proposed 4D-MRI technique to characterize abdominal respiratory motion, which may provide valuable information for radiotherapy planning.

Retracted: Reducing motion artifacts in 4D MR images using principal component analysis (PCA) combined with linear polynomial fitting model.

We have previously developed a retrospective 4D‐MRI technique using body area as the respiratory surrogate, but generally, the reconstructed 4D MR images suffer from severe or mild artifacts mainly caused by irregular motion during image acquisition. Those image artifacts may potentially affect the accuracy of tumor target delineation or the shape representation of surrounding nontarget tissues and organs. So the purpose of this study is to propose an approach employing principal component analysis (PCA), combined with a linear polynomial fitting model, to remodel the displacement vector fields (DVFs) obtained from deformable image registration (DIR), with the main goal of reducing the motion artifacts in 4D MR images. Seven patients with hepatocellular carcinoma (2/7) or liver metastases (5/7) in the liver, as well as a patient with non‐small cell lung cancer (NSCLC), were enrolled in an IRB‐approved prospective study. Both CT and MR simulations were performed for each patient for treatment planning. Multiple‐slice, multiple‐phase, cine‐MRI images were acquired in the axial plane for 4D‐MRI reconstruction. Single‐slice 2D cine‐MR images were acquired across the center of the tumor in axial, coronal, and sagittal planes. For a 4D MR image dataset, the DVFs in three orthogonal direction (inferior–superior (SI), anterior–posterior (AP), and medial–lateral (ML)) relative to a specific reference phase were calculated using an in‐house DIR algorithm. The DVFs were preprocessed in three temporal and spatial dimensions using a polynomial fitting model, with the goal of correcting the potential registration errors introduced by three‐dimensional DIR. Then PCA was used to decompose each fitted DVF into a linear combination of three principal motion bases whose spanned subspaces combined with their projections had been validated to be sufficient to represent the regular respiratory motion. By wrapping the reference MR image using the remodeled DVFs, ‘synthetic’ MR images with reduced motion artifacts were generated at selected phase. Tumor motion trajectories derived from cine‐MRI, 4D CT, original 4D MRI, and ‘synthetic’ 4D MRI were analyzed in the SI, AP, and ML directions, respectively. Their correlation coefficient (CC) and difference (D) in motion amplitude were calculated for comparison. Of all the patients, the means and standard deviations (SDs) of CC comparing ‘synthetic’ 4D MRI and cine‐MRI were 0.98±0.01,0.98±0,01, and 0.99±0.01 in SI, AP, and ML directions, respectively. The mean±SD Ds were 0.59±0.09 mm,0.29±0.10 mm, and 0.15±0.05 mm in SI, AP and ML directions, respectively. The means and SDs of CC comparing ‘synthetic’ 4D MRI and 4D CT were 0.96±0.01,0.95±0.01, and 0.95±0.01 in SI, AP, and ML directions, respectively. The mean±SD Ds were 0.76±0.20 mm,0.33±0.14 mm, and 0.19±0.07 mm in SI, AP, and ML directions, respectively. The means and SDs of CC comparing ‘synthetic’ 4D MRI and original 4D MRI were 0.98±0.01,0.98±0.01, and 0.97±0.01 in SI, AP, and ML directions, respectively. The mean±SD Ds were 0.58±0.10 mm,0.30±0.09 mm, and 0.17±0.04 mm in SI, AP, and ML directions, respectively. In this study we have proposed an approach employing PCA combined with a linear polynomial fitting model to capture the regular respiratory motion from a 4D MR image dataset. And its potential usefulness in reducing motion artifacts and improving image quality has been demonstrated by the preliminary results in oncological patients.PACS numbers: 87.57.cp, 87.57.nj, 87.61.‐c

Four dimensional magnetic resonance imaging with retrospective k-space reordering: a feasibility study.

Current four dimensional magnetic resonance imaging (4D-MRI) techniques lack sufficient temporal/spatial resolution and consistent tumor contrast. To overcome these limitations, this study presents the development and initial evaluation of a new strategy for 4D-MRI which is based on retrospective k-space reordering. We simulated a k-space reordered 4D-MRI on a 4D digital extended cardiac-torso (XCAT) human phantom. A 2D echo planar imaging MRI sequence [frame rate (F) = 0.448 Hz; image resolution (R) = 256 × 256; number of k-space segments (NKS) = 4] with sequential image acquisition mode was assumed for the simulation. Image quality of the simulated "4D-MRI" acquired from the XCAT phantom was qualitatively evaluated, and tumor motion trajectories were compared to input signals. In particular, mean absolute amplitude differences (D) and cross correlation coefficients (CC) were calculated. Furthermore, to evaluate the data sufficient condition for the new 4D-MRI technique, a comprehensive simulation study was performed using 30 cancer patients' respiratory profiles to study the relationships between data completeness (Cp) and a number of impacting factors: the number of repeated scans (NR), number of slices (NS), number of respiratory phase bins (NP), NKS, F, R, and initial respiratory phase at image acquisition (P0). As a proof-of-concept, we implemented the proposed k-space reordering 4D-MRI technique on a T2-weighted fast spin echo MR sequence and tested it on a healthy volunteer. The simulated 4D-MRI acquired from the XCAT phantom matched closely to the original XCAT images. Tumor motion trajectories measured from the simulated 4D-MRI matched well with input signals (D = 0.83 and 0.83 mm, and CC = 0.998 and 0.992 in superior-inferior and anterior-posterior directions, respectively). The relationship between Cp and NR was found best represented by an exponential function (CP=1001-e(-0.18NR) , when NS = 30, NP = 6). At a CP value of 95%, the relative error in tumor volume was 0.66%, indicating that NR at a CP value of 95% (NR,95%) is sufficient. It was found that NR,95% is approximately linearly proportional to NP (r = 0.99), and nearly independent of all other factors. The 4D-MRI images of the healthy volunteer clearly demonstrated respiratory motion in the diaphragm region with minimal motion induced noise or aliasing. It is feasible to generate respiratory correlated 4D-MRI by retrospectively reordering k-space based on respiratory phase. This new technology may lead to the next generation 4D-MRI with high spatiotemporal resolution and optimal tumor contrast, holding great promises to improve the motion management in radiotherapy of mobile cancers.