Mo Kadbi

Technical Note: Characterization and correction of gradient nonlinearity induced distortion on a 1.0 T open bore MRâ SIM

PURPOSE Distortions in magnetic resonance imaging (MRI) compromise spatial fidelity, potentially impacting delineation and dose calculation. We characterized 2D and 3D large field of view (FOV), sequence-independent distortion at various positions in a 1.0 T high-field open MR simulator (MR-SIM) to implement correction maps for MRI treatment planning. METHODS A 36 × 43 × 2 cm(3) phantom with 255 known landmarks (∼1 mm(3)) was scanned using 1.0 T high-field open MR-SIM at isocenter in the transverse, sagittal, and coronal axes, and a 465 × 350 × 168 mm(3) 3D phantom was scanned by stepping in the superior-inferior direction in three overlapping positions to achieve a total 465 × 350 × 400 mm(3) sampled FOV yielding >13 800 landmarks (3D Gradient-Echo, TE/TR/α = 5.54 ms/30 ms/28°, voxel size = 1 × 1 × 2 mm(3)). A binary template (reference) was generated from a phantom schematic. An automated program converted MR images to binary via masking, thresholding, and testing for connectivity to identify landmarks. Distortion maps were generated by centroid mapping. Images were corrected via warping with inverse distortion maps, and temporal stability was assessed. RESULTS Over the sampled FOV, non-negligible residual gradient distortions existed as close as 9.5 cm from isocenter, with a maximum distortion of 7.4 mm as close as 23 cm from isocenter. Over six months, average gradient distortions were -0.07 ± 1.10 mm and 0.10 ± 1.10 mm in the x and y directions for the transverse plane, 0.03 ± 0.64 and -0.09 ± 0.70 mm in the sagittal plane, and 0.4 ± 1.16 and 0.04 ± 0.40 mm in the coronal plane. After implementing 3D correction maps, distortions were reduced to <1 pixel width (1 mm) for all voxels up to 25 cm from magnet isocenter. CONCLUSIONS Inherent distortion due to gradient nonlinearity was found to be non-negligible even with vendor corrections applied, and further corrections are required to obtain 1 mm accuracy for large FOVs. Statistical analysis of temporal stability shows that sequence independent distortion maps are consistent within six months of characterization.

Dosimetric and workflow evaluation of first commercial synthetic CT software for clinical use in pelvis

To evaluate a commercial synthetic CT (syn-CT) software for use in prostate radiotherapy. Twenty-five prostate patients underwent CT and MR simulation scans in treatment position on a 3T MR scanner. A commercially available MR protocol was used that included a T2w turbo spin-echo sequence for soft-tissue contrast and a dual echo 3D mDIXON fast field echo (FFE) sequence for generating syn-CT. A dual-echo 3D FFE B 0 map was used for patient-induced susceptibility distortion analysis and a new 3D balanced-FFE sequence was evaluated for identification of implanted gold fiducial markers and subsequent image-guidance during radiotherapy delivery. Tissues were classified as air, adipose, water, trabecular/spongy bone and compact/cortical bone and assigned bulk HU values. The accuracy of syn-CT for treatment planning was analyzed by transferring the structures and plan from planning CT to syn-CT and recalculating the dose. Accuracy of localization at the treatment machine was evaluated by comparing registration of kV radiographs to either digitally reconstructed radiographs (DRRs) generated from syn-CT or traditional DRRs generated from the planning CT. Similarly, accuracy of setup using CBCT and syn-CT was compared to that using the planning CT. Finally, a MR-only simulation workflow was established and end-to-end testing was completed on five patients undergoing MR-only simulation. Dosimetric comparison between the original CT and syn-CT plans was within 0.5% on average for all structures. The de-novo optimized plans on the syn-CT met institutional clinical objectives for target and normal structures. Patient-induced susceptibility distortion based on B 0 maps was within 1 mm and 0.5 mm in the body and prostate respectively. DRR and CBCT localization based on MR-localized fiducials showed a standard deviation of  <1 mm. End-to-end testing and MR simulation workflow was successfully validated. MRI derived synthetic CT can be successfully used for a MR-only planning and treatment for prostate radiotherapy.

Magnetic Resonance-Based Automatic Air Segmentation for Generation of Synthetic Computed Tomography Scans in the Head Region.

PURPOSE To incorporate a novel imaging sequence for robust air and tissue segmentation using ultrashort echo time (UTE) phase images and to implement an innovative synthetic CT (synCT) solution as a first step toward MR-only radiation therapy treatment planning for brain cancer. METHODS AND MATERIALS Ten brain cancer patients were scanned with a UTE/Dixon sequence and other clinical sequences on a 1.0 T open magnet with simulation capabilities. Bone-enhanced images were generated from a weighted combination of water/fat maps derived from Dixon images and inverted UTE images. Automated air segmentation was performed using unwrapped UTE phase maps. Segmentation accuracy was assessed by calculating segmentation errors (true-positive rate, false-positive rate, and Dice similarity indices using CT simulation (CT-SIM) as ground truth. The synCTs were generated using a voxel-based, weighted summation method incorporating T2, fluid attenuated inversion recovery (FLAIR), UTE1, and bone-enhanced images. Mean absolute error (MAE) characterized Hounsfield unit (HU) differences between synCT and CT-SIM. A dosimetry study was conducted, and differences were quantified using γ-analysis and dose-volume histogram analysis. RESULTS On average, true-positive rate and false-positive rate for the CT and MR-derived air masks were 80.8% ± 5.5% and 25.7% ± 6.9%, respectively. Dice similarity indices values were 0.78 ± 0.04 (range, 0.70-0.83). Full field of view MAE between synCT and CT-SIM was 147.5 ± 8.3 HU (range, 138.3-166.2 HU), with the largest errors occurring at bone-air interfaces (MAE 422.5 ± 33.4 HU for bone and 294.53 ± 90.56 HU for air). Gamma analysis revealed pass rates of 99.4% ± 0.04%, with acceptable treatment plan quality for the cohort. CONCLUSIONS A hybrid MRI phase/magnitude UTE image processing technique was introduced that significantly improved bone and air contrast in MRI. Segmented air masks and bone-enhanced images were integrated into our synCT pipeline for brain, and results agreed well with clinical CTs, thereby supporting MR-only radiation therapy treatment planning in the brain.

4D UTE flow: a phase-contrast MRI technique for assessment and visualization of stenotic flows.

Inaccuracy of conventional four‐dimensional (4D) flow MR imaging in the presence of random unsteady and turbulent blood flow distal to a narrowing has been an important challenge. Previous investigations have revealed that shorter echo times (TE) decrease the errors, leading to more accurate flow assessments.

Technical Note: Characterization and correction of gradient nonlinearity induced distortion on a 1.0 T open bore MR‐SIM

PURPOSE Distortions in magnetic resonance imaging (MRI) compromise spatial fidelity, potentially impacting delineation and dose calculation. We characterized 2D and 3D large field of view (FOV), sequence-independent distortion at various positions in a 1.0 T high-field open MR simulator (MR-SIM) to implement correction maps for MRI treatment planning. METHODS A 36 × 43 × 2 cm(3) phantom with 255 known landmarks (∼1 mm(3)) was scanned using 1.0 T high-field open MR-SIM at isocenter in the transverse, sagittal, and coronal axes, and a 465 × 350 × 168 mm(3) 3D phantom was scanned by stepping in the superior-inferior direction in three overlapping positions to achieve a total 465 × 350 × 400 mm(3) sampled FOV yielding >13 800 landmarks (3D Gradient-Echo, TE/TR/α = 5.54 ms/30 ms/28°, voxel size = 1 × 1 × 2 mm(3)). A binary template (reference) was generated from a phantom schematic. An automated program converted MR images to binary via masking, thresholding, and testing for connectivity to identify landmarks. Distortion maps were generated by centroid mapping. Images were corrected via warping with inverse distortion maps, and temporal stability was assessed. RESULTS Over the sampled FOV, non-negligible residual gradient distortions existed as close as 9.5 cm from isocenter, with a maximum distortion of 7.4 mm as close as 23 cm from isocenter. Over six months, average gradient distortions were -0.07 ± 1.10 mm and 0.10 ± 1.10 mm in the x and y directions for the transverse plane, 0.03 ± 0.64 and -0.09 ± 0.70 mm in the sagittal plane, and 0.4 ± 1.16 and 0.04 ± 0.40 mm in the coronal plane. After implementing 3D correction maps, distortions were reduced to <1 pixel width (1 mm) for all voxels up to 25 cm from magnet isocenter. CONCLUSIONS Inherent distortion due to gradient nonlinearity was found to be non-negligible even with vendor corrections applied, and further corrections are required to obtain 1 mm accuracy for large FOVs. Statistical analysis of temporal stability shows that sequence independent distortion maps are consistent within six months of characterization.

Four dimensional magnetic resonance imaging optimization and implementation for magnetic resonance imaging simulation

PURPOSE Precise radiation therapy for abdominal lesions is complicated by respiratory motion and suboptimal soft tissue contrast from 4-dimensional (4D) computed tomography, whereas 4D magnetic resonance imaging MRI (4DMRI) provides superior tissue discrimination. This work evaluates a novel 4DMRI algorithm for motion management in radiation therapy. METHODS AND MATERIALS Respiratory-triggered, T2-weighted, single-shot 4DMRI was evaluated for an open 1.0T magnetic resonance simulation platform. An in-house programmable platform was devised that translated objects for a variety of breathing patterns. Coronal 4DMRIs were acquired to evaluate the impact of number of phases on excursion and scan time. The impact of breathing period and regularity on scan time was assessed. A novel clinical 4D prototype phantom was scanned to characterize excursion and absolute volume differences between phase acquisitions. Optimized parameters were applied to abdominal 4DMRIs of 5 volunteers and 2 abdominal cancer patients on an institutional review board-approved protocol. Duty cycle, scan time, and waveform analysis were evaluated. Maximum intensity projection datasets were analyzed. RESULTS Two- to 5-fold acquisition time increase was measured for 10-phase versus 2-phase phantom experiments. Regular breathing patterns yielded higher duty cycles than irregular (48.5% and 35.9%, respectively, P < .001), whereas faster breathing rates yielded shorter 4DMRI acquisition times. Volumes of a hypodense target were underestimated 4% to 5% for 2 and 4 phases compared with 10 phases. Better agreement was obtained for 6- and 8-phase acquisitions (~3% different from 10 phase). Internal target volume centroids on minimum and maximum images across all phases were <2 mm different across all 10 phases, although slight target excursion variations (up to 4 mm) were observed. In humans, a strong negative association between breathing rate and acquisition time (Pearsons r = -0.68, P < .05) was observed. Eight-phase acquisition times ranged from 7 to 15 minutes, depending on the patient. CONCLUSION 4DMRI has been optimized and implemented. Irregular breathing patterns and slow breathing rate adversely impacted 4DMRI efficiency; thus, interventions such as biofeedback may be desirable.

A novel phase-corrected 3D cine ultra-short te (UTE) phase-contrast MRI technique

Phase-contrast (PC) MRI is a non-invasive technique to assess cardiovascular blood flow. However, this technique is not accurate for instance at the carotid bifurcation due to turbulent and disturbed blood flow in atherosclerotic disease. Flow quantification using conventional PC MRI distal to stenotic vessels suffers from intravoxel dephasing and flow artifacts. Previous studies have shown that short echo time (TE) potentially decreases the phase errors. In this work, a novel 3D cine UTE-PC imaging method is designed to measure the blood velocity in the carotid bifurcation using a UTE center-out radial trajectory and short TE time compared to standard PC MRI sequences. With a new phase error correction technique based on autocorrelation method, the proposed 3D cine UTE-PC has the potential to achieve high accuracy for quantification and visualization of velocity jet distal to a stenosis. Herein, we test the feasibility of the method in determining accurate flow waveforms in normal volunteers.

Assessment of flow and hemodynamics in the carotid artery using a reduced TE 4D flow spiral phase-contrast MRI

4D flow MRI is a powerful technique for quantitative flow assessment and visualization of complex flow patterns and hemodynamics of cardiovascular flows. This technique results in more anatomical information and comprehensive assessment of blood flow. However, conventional 4D PC MRI suffers from a few obstacles for clinical applications. The total scan time is long, especially in large volumes with high spatial resolutions. Inaccuracy of conventional Cartesian PC MRI in the setting of atherosclerosis and in general, disturbed and turbulent blood flow is another important challenge. This inaccuracy is the consequence of signal loss, intravoxel dephasing and flow-related artifact in the presence of disturbed and turbulent flow. Spiral k-space trajectory has valuable attributes which can help overcome some of the problems with 4D flow Cartesian acquisitions. Spiral trajectory benefits from shorter TE and reduces the flow-related artifacts. In addition, short spiral readouts with spiral interleaves can significantly reduce the total scan time, reducing the chances of patient motion which may also corrupt the data in the form of motion artifacts. In this paper, the accuracy of flow assessment and flow visualization with reduced TE 4D Spiral PC was investigated and good agreement was observed between the spiral and conventional technique. The systolic mean velocity, peak flow and the average flow in CCA and ICA of normal volunteers using 4D spiral PC MRI showed errors less than 10% compared to conventional 4D PC MRI. In addition, the scan time using spiral sequence was 3:31 min which is half of the time using conventional sequence.

Comparison of Cartesian, UTE radial, and Spiral Phase-Contrast MRI in Measurement of Blood Flow in Extracranial Carotid Arteries: Normal Subjects

Use of Phase contrast (PC) MRI in measurement of blood flow has significant clinical importance. In this paper, we compare the accuracy of the conventional approach to flow imaging to two de novo approaches in 3 normal subjects in the common, internal, and external carotid arteries and discuss and demonstrate the advantages and disadvantages of each method. The conventional PC sequence adopts a Cartesian read-out in k-space and requires longer acquisitions but exhibits flow artifacts in the setting of stenotic and disturbed flow. Spiral PC collects k-space data using spiral readout and is capable of reducing the TR and TE in order to minimize the total imaging time. Despite its efficiency in scan time, in the single shot mode, this technique suffers from off-resonance and inconsistent data artifacts. Use of multiple short spiral arms for providing k-space coverage resolves these issues. Ultra short TE (UTE) PC MRI is a novel technique which adopts a radial trajectory and provides improvements to the standard radial acquisition by reducing the echo time to less than 1 ms through combination of flow encoding and slice select gradients and by immediate sampling of the FID during readout. The ultra-short echo times, improves on intravoxel spin dephasing due to fluid mixing observed in imaging of disturbed flow and stenotic jets. Despite its capability of achieving the shortest TE, this method is hindered by longer acquisition times and phase corruption errors. We mitigate this by a novel 3-D acquisition which includes a phase correction step. All three approaches were found to be able to quantify the normal Carotid flow waveform with high accuracy.

Rapid flow quantification in iliac arteries with spiral phase-contrast MRI

Phase contrast MRI is a powerful tool for blood flow quantification. Conventional cartesian phase contrast sequences require lengthy acquisition on the order of several minutes. Spiral acquisition phase-contrast (PC) MRI is capable of reducing the TR and TE in order to minimize flow dependent artifacts and total imaging time. Despite this, in general, spiral phase contrast sequences suffer from off-resonance artifacts and inconsistent data artifacts. In this work, we show that short interleaved spiral readout trajectories have the capability to obtain high spatio-temporal resolution flow images in the common iliac artery distal to the aortoiliac bifurcation with little or no artifacts and with significant savings in image acquisition time over the Cartesian trajectory. To verify the accuracy, we compare our results with a Conventional cartesian trajectory.