ABSTRACT
Background: Magnetic resonance elastography (MRE) is a non-invasive method to measure the viscoelastic properties of tissue and has been applied in multiple abdominal organs. However, abdominal MRE suffers from detrimental breathing motion causing misalignment of structures between repeated acquisitions for different MRE dimensions (e.g., motion encoding directions and wave phase offsets). This study investigated motion correction strategies to resolve all breathing motion on sagittal free-breathing MRE acquisitions in a phantom, in healthy volunteers and showed feasibility in patients. Methods: First, in silico experiments were performed on a static phantom dataset with simulated motion. Second, eight healthy volunteers underwent two sagittal MRE acquisitions in the pancreas and right kidney. The multi-frequency free-breathing spin-echo echo-planar-imaging (SE-EPI) MRE consisted of four frequencies (30, 40, 50, 60 Hz), eight wave-phase offsets, with 3 mm3 isotropic voxel size. Following data re-sorting in different number of motion states (4 till 12) based on respiratory waveform signal, three intensity-based registration methods (monomodal, multimodal, and phase correlation) and non-rigid local registration were compared. A ranking method was used to determine the best registration method, based on seven signal-to-noise and image quality measures. Repeatability was assessed for no motion correction (Original) and the best performing method (Best) using Bland-Altman analysis. Lastly, the best motion correction method was compared to no motion correction on patient MRE data [pancreatic ductal adenocarcinoma (PDAC, n=5) and metabolic dysfunction-associated steatotic liver disease (MASLD) (n=1)]. Results: In silico experiments showed a deviation of shear wave speed (SWS) with simulated motion to the ground truth, which was (partially) resolved using motion correction. In healthy volunteers ranking resulted in the best motion correction method of monomodal registration using nine motion states, while no motion correction was ranked last. Limits of agreement were (-0.18, 0.14), and (-0.25, 0.18) m/s for Best and Original, respectively. Using motion correction in patients resulted in a significant increase in SWS in the pancreas (Original: 1.39±0.10 and Best: 1.50±0.17 m/s). After motion correction PDAC had a mean SWS of 1.56±0.27 m/s (Original: 1.42±0.25 m/s). The fibrotic liver mean SWS was 2.07±0.20 m/s (Original: 2.12±0.18 m/s). Conclusions: Motion correction in sagittal free-breathing abdominal MRE results in improved data quality, inversion precision, repeatability, and is feasible in patients.
ABSTRACT
BACKGROUND: Magnetic resonance elastography (MRE) is used to non-invasively estimate biomechanical tissue properties via the imaging of propagating mechanical shear waves. Several factors including mechanical transducer design, MRI sequence design and viscoelastic reconstruction influence data quality and hence the reliability of the derived biomechanical properties. PURPOSE: To design and characterize a novel mechanical MRE transducer concept based on a rotational eccentric mass, coined the gravitational transducer. MATERIALS AND METHODS: Table top measurements were performed using accelerometers to characterize the frequency response of the new transducer concept at different driving frequencies (fâVIB) and different rotating masses. These were compared to a commercially available pneumatically driven MRE transducer. MR data were acquired on a 3T scanner using a fractionally encoded gradient echo MRE sequence in three healthy volunteers. Acceleration and displacement spectra were plotted in units of g and mm, respectively, and visually compared, emphasizing the ratio between the peaks at fâVIB and its 2nd harmonic, a known cause of error in the reconstruction of biomechanical properties as is explored in more detail in numerical simulations here. No formal statistical testing was performed in this proof-of-principle paper. RESULTS: The new transducer concept shows-as expected from theory-a quadratic or linear increase of acceleration amplitude with increase in fâVIB or mass, respectively. Furthermore, different versions of the transducer show markedly lower 2nd harmonic-to-fâVIB ratios compared to the commercially available pneumatically driven transducer. Displacement was constant over a range of fâVIB, in accordance with theory. Phantom and in vivo data show low nonlinearity and excellent data quality. CONCLUSION: The table top measurements are in concordance with the theory behind a transducer based on a rotational eccentric mass. The resulting constant displacement amplitude irrespective of fâVIB and low 2nd harmonic-to-fâVIB ratio result in low nonlinearity and high data fidelity in both phantom and in vivo examples.
