Geometric validation of self-gating k-space-sorted 4D-MRI vs 4D-CT using a respiratory motion phantom

Scoping review of magnetic resonance motion imaging phantoms

To review and analyze the currently available MRI motion phantoms. Publications were collected from the Toronto Metropolitan University Library, PubMed, and IEEE Xplore. Phantoms were categorized based on the motions they generated: linear/cartesian, cardiac-dilative, lung-dilative, rotational, deformation or rolling. Metrics were extracted from each publication to assess the motion mechanisms, construction methods, as well as phantom validation. A total of 60 publications were reviewed, identifying 48 unique motion phantoms. Translational movement was the most common movement (used in 38% of phantoms), followed by cardiac-dilative (27%) movement and rotational movement (23%). The average degrees of freedom for all phantoms were determined to be 1.42. Motion phantom publications lack quantification of their impact on signal-to-noise ratio through standardized testing. At present, there is a lack of phantoms that are designed for multi-role as many currently have few degrees of freedom.

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

PURPOSE

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.

METHODS AND MATERIALS

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.

RESULTS

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.5 mm, 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.43 mm in the superior-inferior (SI) direction, and 0.59 ± 0.36, 0.51 ± 0.30, and 0.50 ± 0.24 mm 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.5 mm in SI, AP, and left-right (LR) directions.

CONCLUSION

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.

Evaluating differences in respiratory motion estimates during radiotherapy: a single planning 4DMRI versus daily 4DMRI

Background

In radiotherapy, respiratory-induced tumor motion is typically measured using a single four-dimensional computed tomography acquisition (4DCT). Irregular breathing leads to inaccurate motion estimates, potentially resulting in undertreatment of the tumor and unnecessary dose to healthy tissue. The aim of the research was to determine if a daily pre-treatment 4DMRI-strategy led to a significantly improved motion estimate compared to single planning 4DMRI (with or without outlier rejection).

Methods

4DMRI data sets from 10 healthy volunteers were acquired. The first acquisition simulated a planning MRI, the respiratory motion estimate (constructed from the respiratory signal, i.e. the 1D navigator) was compared to the respiratory signal in the subsequent scans (simulating 5–29 treatment fractions). The same procedure was performed using the first acquisition of each day as an estimate for the subsequent acquisitions that day (2 per day, 4–20 per volunteer), simulating a daily MRI strategy. This was done for three outlier strategies: no outlier rejection (NoOR); excluding 5% of the respiratory signal whilst minimizing the range (Min95) and excluding the datapoints outside the mean end-inhalation and end-exhalation positions (MeanIE).

Results

The planning MRI median motion estimates were 27 mm for NoOR, 18 mm for Min95, and 13 mm for MeanIE. The daily MRI median motion estimates were 29 mm for NoOR, 19 mm for Min95 and 15 mm for MeanIE. The percentage of time outside the motion estimate were for the planning MRI: 2%, 10% and 32% for NoOR, Min95 and MeanIE respectively. These values were reduced with the daily MRI strategy: 0%, 6% and 17%. Applying Min95 accounted for a 30% decrease in motion estimate compared to NoOR.

Conclusion

A daily MRI improved the estimation of respiratory motion as compared to a single 4D (planning) MRI significantly. Combining the Min95 technique with a daily 4DMRI resulted in a decrease of inclusion time of 6% with a 30% decrease of motion. Outlier rejection alone on a planning MRI often led to underestimation of the movement and could potentially lead to an underdosage.

Trial registration: protocol W15_373#16.007

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-021-01915-1.

Design, Implementation, and Validation of a Pulsatile Heart Phantom Pump

The developments in Computed Tomography (CT) and Magnetic Resonance allow visualization of blood flow in vivo using these techniques. However, validation tests are needed to determine a gold standard. For the validation tests, controllable systems that can generate pulsatile flow are needed. In this study, we aimed to develop an affordable pulsatile pump and an artificial circulatory system to simulate the blood flow for validation purposes. Initially, the prerequisites for the phantom were pulsating flow output equal to that of the human cardiac pulse pattern; the flow pattern of the mimicked cardiac output should be equal to that of a human, a variable stroke volume (40-120 ml/beat), and a variable heart rate (60-170 bpm). The developed phantom setup was tested with CT scanner. A washout profile was created based on the image intensity of the selected slice. The test was successful for a heart rate of 70 bpm and a stroke volume of 68 ml, but the system failed to work at various heartbeats and stroke volumes. This was due to the problems with software of the microcontroller. As conclusion in this study, we present a proof of concept for a pulsatile heart phantom pump that can be used in validation tests.

A modern review of the uncertainties in volumetric imaging of respiratory-induced target motion in lung radiotherapy

Radiotherapy has become a critical component for the treatment of all stages and types of lung cancer, often times being the primary gateway to a cure. However, given that radiation can cause harmful side effects depending on how much surrounding healthy tissue is exposed, treatment of the lung can be particularly challenging due to the presence of moving targets. Careful implementation of every step in the radiotherapy process is absolutely integral for attaining optimal clinical outcomes. With the advent and now widespread use of stereotactic body radiation therapy (SBRT), where extremely large doses are delivered, accurate, and precise dose targeting is especially vital to achieve an optimal risk to benefit ratio. This has largely become possible due to the rapid development of image-guided technology. Although imaging is critical to the success of radiotherapy, it can often be plagued with uncertainties due to respiratory-induced target motion. There has and continues to be an immense research effort aimed at acknowledging and addressing these uncertainties to further our abilities to more precisely target radiation treatment. Thus, the goal of this article is to provide a detailed review of the prevailing uncertainties that remain to be investigated across the different imaging modalities, as well as to highlight the more modern solutions to imaging motion and their role in addressing the current challenges.

A novel amplitude binning strategy to handle irregular breathing during 4DMRI acquisition: improved imaging for radiotherapy purposes

Background

For radiotherapy of abdominal cancer, four-dimensional magnetic resonance imaging (4DMRI) is desirable for tumor definition and the assessment of tumor and organ motion. However, irregular breathing gives rise to image artifacts. We developed a outlier rejection strategy resulting in a 4DMRI with reduced image artifacts in the presence of irregular breathing.

Methods

We obtained 2D T2-weighted single-shot turbo spin echo images, with an interleaved 1D navigator acquisition to obtain the respiratory signal during free breathing imaging in 2 patients and 12 healthy volunteers. Prior to binning, upper and lower inclusion thresholds were chosen such that 95% of the acquired images were included, while minimizing the distance between the thresholds (inclusion range (IR)). We compared our strategy (Min95) with three commonly applied strategies: phase binning with all images included (Phase), amplitude binning with all images included (MaxIE), and amplitude binning with the thresholds set as the mean end-inhale and mean end-exhale diaphragm positions (MeanIE). We compared 4DMRI quality based on:

Data included (DI); percentage of images remaining after outlier rejection.

Reconstruction completeness (RC); percentage of bin-slice combinations containing at least one image after binning.

Intra-bin variation (IBV); interquartile range of the diaphragm position within the bin-slice combination, averaged over three central slices and ten respiratory bins.

IR.

Image smoothness (S); quantified by fitting a parabola to the diaphragm profile in a sagittal plane of the reconstructed 4DMRI.

A two-sided Wilcoxon’s signed-rank test was used to test for significance in differences between the Min95 strategy and the Phase, MaxIE, and MeanIE strategies.

Results

Based on the fourteen subjects, the Min95 binning strategy outperformed the other strategies with a mean RC of 95.5%, mean IBV of 1.6 mm, mean IR of 15.1 mm and a mean S of 0.90. The Phase strategy showed a poor mean IBV of 6.2 mm and the MaxIE strategy showed a poor mean RC of 85.6%, resulting in image artifacts (mean S of 0.76). The MeanIE strategy demonstrated a mean DI of 85.6%.

Conclusions

Our Min95 reconstruction strategy resulted in a 4DMRI with less artifacts and more precise diaphragm position reconstruction compared to the other strategies.

Trial registration

Volunteers: protocol W15_373#16.007; patients: protocol NL47713.018.14

Electronic supplementary material

The online version of this article (10.1186/s13014-019-1279-z) contains supplementary material, which is available to authorized users.

Novel 4D-MRI of tumor infiltrating vasculature: characterizing tumor and vessel volume motion for selective boost volume definition in pancreatic radiotherapy

BACKGROUND

Pancreatic ductal adenocarcinoma has dismal prognosis. Most patients receive radiation therapy (RT), which is complicated by respiration induced organ motion in upper abdomen. The purpose of this study is to report our early clinical experience in a novel self-gated k-space sorted four-dimensional magnetic resonance imaging (4D-MRI) with slab-selective (SS) excitation to highlight tumor infiltrating blood vessels for pancreatic RT.

METHODS

Ten consecutive patients with borderline resectable or locally advanced pancreatic cancer were recruited to the study. Non-contrast 4D-MRI with and without slab-selective excitation and 4D-CT with delay contrast were performed on all patients. Vessel-tissue CNR were calculated for aorta and critical vessels (superior mesenteric artery or superior mesenteric vein) encompassed by tumor. Respiratory motion trajectories for tumor, as well as involved vessels were analyzed on SS-4D-MRI. Intra-class cross correlation (ICC) between tumor volume and involved vessels were calculated.

RESULTS

Among all 4D imaging modalities evaluated, SS-4D-MRI sampling trajectory results in images with highest vessel-tissue CNR comparing to non-slab-selective 4D-MRI and 4D-CT for all patients studied. Average (±standard deviation) CNR for involved vessels are 13.1 ± 8.4 and 3.2 ± 2.7 for SS-4D-MRI and 4D-CT, respectively. The ICC factors comparing tumor and involved vessels motion trajectories are 0.93 ± 0.10, 0.65 ± 0.31 and 0.77 ± 0.23 for superior-inferior, anterior-posterior and medial-lateral directions respectively.

CONCLUSIONS

A novel 4D-MRI sequence based on 3D-radial sampling and slab-selective excitation has been assessed for pancreatic cancer patients. The non-contrast 4D-MRI images showed significantly better contrast to noise ratio for the vessels that limit tumor resectability compared to 4D-CT with delayed contrast. The sequence has great potential in accurately defining both the tumor and boost volume margins for pancreas RT with simultaneous integrated boost.

Computed tomography trachea volumetry in patients with scleroderma: Association with clinical and functional findings

BACKGROUND

In scleroderma, excessive collagen production can alter tracheal geometry, and computed tomography (CT) volumetry of this structure may aid in detecting possible abnormalities. The objectives of this study were to quantify the morphological abnormalities in the tracheas of ​​patients with scleroderma and to correlate these findings with data on clinical and pulmonary function.

METHODS

This was a cross-sectional study in which 28 adults with scleroderma and 27 controls matched by age, gender and body mass index underwent chest CT with posterior segmentation and skeletonization of the images. In addition, all participants underwent pulmonary function tests and clinical evaluation, including the modified Rodnan skin score (mRSS).

RESULTS

Most patients (71.4%) had interstitial lung disease on CT. Compared to controls, patients with scleroderma showed higher values ​​in the parameters measured by CT trachea volumetry, including area, eccentricity, major diameter, minor diameter, and tortuosity. The tracheal area and equivalent diameter were negatively correlated with the ratio between forced expiratory flow and forced inspiratory flow at 50% of forced vital capacity (FEF50%/FIF50%) (r = -0.44, p = 0.03 and r = -0.46, p = 0.02, respectively). The tracheal tortuosity was negatively correlated with peak expiratory flow (r = -0.51, p = 0.008). The mRSS showed a positive correlation with eccentricity (r = 0.62, p < 0.001) and tracheal tortuosity (r = 0.51, p = 0.007), while the presence of anti-topoisomerase I antibody (ATA) showed a positive correlation with tracheal tortuosity (r = 0.45, p = 0.03).

CONCLUSIONS

In a sample composed predominantly of scleroderma patients with associated interstitial lung disease, there were abnormalities in tracheal geometry, including greater eccentricity, diameter and tortuosity. In these patients, abnormalities in the geometry of the trachea were associated with functional markers of obstruction. In addition, tracheal tortuosity was correlated with cutaneous involvement and the presence of ATA.

Design and evaluation of an MRI-compatible linear motion stage

PURPOSE

To develop and evaluate a tool for accurate, reproducible, and programmable motion control of imaging phantoms for use in motion sensitive magnetic resonance imaging (MRI) appli cations.

METHODS

In this paper, the authors introduce a compact linear motion stage that is made of nonmagnetic material and is actuated with an ultrasonic motor. The stage can be positioned at arbitrary positions and orientations inside the scanner bore to move, push, or pull arbitrary phantoms. Using optical trackers, measuring microscopes, and navigators, the accuracy of the stage in motion control was evaluated. Also, the effect of the stage on image signal-to-noise ratio (SNR), artifacts, and B0 field homogeneity was evaluated.

RESULTS

The error of the stage in reaching fixed positions was 0.025 ± 0.021 mm. In execution of dynamic motion profiles, the worst-case normalized root mean squared error was below 7% (for frequencies below 0.33 Hz). Experiments demonstrated that the stage did not introduce artifacts nor did it degrade the image SNR. The effect of the stage on the B0 field was less than 2 ppm.

CONCLUSIONS

The results of the experiments indicate that the proposed system is MRI-compatible and can create reliable and reproducible motion that may be used for validation and assessment of motion related MRI applications.