Tomokazu Numano1,2, Daiki Ito2,3,4, Kazuyuki Mizuhara5, Toshikatsu Washio2, Tetsushi Habe3,4, Toshiki Maeno3, Masaki Misawa2, and Naotaka Nitta2
1Radiological Sciences, Tokyo Metropolitan University, Tokyo, Japan, 2National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Japan, 3Tokyo Metropolitan University, Tokyo, Japan, 4Keio University Hospital, Tokyo, Japan, 5Tokyo Denki University, Tokyo, Japan
Synopsis
We developed a new technique for dynamic MR elastography (MRE) using a MR magnitude image. A general MRE uses a MR phase image as a wave image. Proposed technique was used a MR magnitude image instead of MR phase image as a wave image by using a special vibration. Proposed technique (special vibration MRE) performance was comparable to that of a continuous vibration MRE. Since special vibration MRE need an only few second vibration, it could dramatically eliminated the patients' vibration-related discomfort.
INTRODUCTION
Magnetic resonance elastography (MRE) is a non-invasive technique for
measuring the mechanical properties of soft tissues. Assuming that biologic
tissues are homogeneous, isotropic and linearly elastic, μ is calculated as follows: μ=ρλ2f2 where λ is the wavelength of shear wave (m), f is the frequency of the
shear wave (Hz) and ρ is the soft
tissue density (assumed to be 1000 kg/m3). Since the frequency and
the density were already known, if the wavelength is imaged, the shear modulus
can be calculated. For general MRE, continuous harmonic mechanical excitation
is applied to the tissue. In order to detect the propagating shear wave in the
tissue, a motion encoding gradient is used to encode the shear wave information
onto the MR phase image (wave image). When the accumulated MR phase sifts are
out of the [-π, π] range, the true phase values will be wrapped back to this
range, creating discontinuities in the MR phase image and subsequently in the
estimated motion. Thus the phase needs to be unwrapped before the motion can be
estimated accurately. However, the out of the range of [-π, π] is also caused
by significantly different magnetic susceptibility. If the MR phase shifts and
the significantly different magnetic susceptibility region are overlapped,
there is a possibility that an error occurs in the phase unwrapping process.
Against this background, we developed a novel MRE technique which used only MR
magnitude images instead MR phase images. The MR magnitude image MRE (MI-MRE)
was used a special vibration technique instead a continuous vibration. This
study shows how the observed pattern of MR magnitude signal variations in
MI-MRE related to the underlying shear wave motion, and how these variations
can be effected by the special vibrations.METHODS
Our previous study reported a new method
for MRE using the GRE-MultiEcho-MRE sequence without motion encoding gradient
(MEG)1. The GRE-MultiEcho-MRE uses a series of echoes acquired as a
train following after a single excitation pulse. The multiple symmetrical
gradient-echoes in the GRE-MultiEcho-MRE are acquired by symmetrical bipolar
readout gradient (GR). This GR has a function comparable to MEG (MEG-like
effect). In the all MRE, the vibration power was giving continuous throughout
the whole acquisition. For a new attempt, we changed the vibration of the
GRE-MultiEcho-MRE from continuous to the special vibration (Fig.1). Specifically,
arbitrary Gaussian distributions of vibration power were formed on half an
acquisition time. Interestingly, it was observed that the propagating wave
pattern was visualized on MR magnitude image. Glaser et al. (2003) reported a
phenomenon that the propagating wave pattern was visualized on MR magnitude
image as intravoxel phase dispersion (IVPD)2. We found the
phenomenon by using the special vibration was an entirely different from IVPD.
All MRE experiments were performed on a
clinical MR imager (Achieva 3.0 T, Philips) using a flex-M coil (gel-phantom)
and a torso coil (volunteer lower back). A self-made waveform generation system
(LabVIEW, USB-6221, National Instruments) was used to generate the vibration
waveform. This system is capable of generating sinusoidal waveforms with
arbitrary frequencies and phases. For the special vibration, the waveform
generation system can control arbitrary Gaussian distributions of output
voltage at during acquisition. Power amplifier (XTi 1000, Crown) and a
pneumatic pressure generator (Subwoofer TIT320C-4 12”, Dayton Audio) units were
used to supply vibrations to a vibration pad. The vibration pad was designed
using a three-dimensional printer (3D touch; 3D system) in order to adapt to
the gel-phantom and the lower back region. The GRE-MultiEcho-MRE sequence
parameters were TR, 40ms; 1st TE, 2.2ms; dTE, 3.3ms@150Hz, 5ms@100Hz, 10ms@50Hz;
flip angle, 20degree; scan matrix, 512×256; image matrix, 512×512; vibration
frequency, 50, 100, 150Hz; vibration phase offset, 4; total acquisition time,
82s; MEG-like effect direction, L-R. All elastograms were produced by Local
Frequency Estimate (LFE) algorithm freeware (MRE/Wave, MAYO CLINIC).RESULTS and DISCUSSION
Figure 2 shows the MI-MRE of gel-phantom
by using 100 Hz vibration. In Fig.2a, continuous vibration, the shear wave information
was visualized on MR phase image. Since the shear wave information was never
visualized on MR magnitude image, the maximum power was well adjusted for
prevent IVPD. In Fig.2b, special vibration, the shear wave information was
visualized on both images. Since the maximum vibration power of special
vibration (peak of Gaussian distribution) was set to the power of continuous
vibration, the IVPD was prevented on this gel-phantom experiment.
Figure 3 shows spatial frequency analyses of
each vibration frequency and variation of special vibration. Figure 4 shows the
difference between the k-space data of continuous vibration and the one of
special vibration. When the vibration frequency increases, the wavelength is
shortened, and high spatial frequency information is required. Since the gap of
the difference domain of k-space (white arrow) was broadening dependent on
vibration frequency, we considered that the difference domain of k-space (blue
arrow) has key component of MI-MRE. Figure 5 shows the difference domain of
k-space (blue arrow) and the width of special vibration (blue belt). In
Amplitude Convex 05 had only convolved shear wave information in the
low-spatial frequency domain, blue arrows and blue belt was not overlapped on
100, 150 Hz, the image contrast of MI-MRE decreased (Fig.3e).Acknowledgements
This work was supported by JSPS KAKENHI Grant NumberJP19K09579, Japan.References
- Numano T, Mizuhara K, Hata
J, et al. A simple method for MR elastography: a gradient-echo type multi-echo
sequence. Magn Reson Imaging 2015;33(1):31-7
- Glaser KJ, Felmlee JP, Manduca A, et al. Shear stiffness estimation
using intravoxel phase dispersion in magnetic resonance elastography. Magn
Reson Med. 2003 Dec;50(6):1256-65.