Ramin S. Sahebjavaher

Rapid acquisition of multifrequency, multislice and multidirectional MR elastography data with a fractionally encoded gradient echo sequence

In MR elastography (MRE), periodic tissue motion is phase encoded using motion‐encoding gradients synchronized to an externally applied periodic mechanical excitation. Conventional methods result in extended scan time for quality phase images, thus limiting the broad application of MRE in the clinic. For practical scan times, researchers have been relying on one‐dimensional or two‐dimensional motion‐encoding, low‐phase sampling and a limited number of slices, and artifact‐prone, single‐shot, echo planar imaging (EPI) readout. Here, we introduce a rapid multislice pulse sequence capable of three‐dimensional motion encoding that is also suitable for simultaneously encoding motion with multiple frequency components. This sequence is based on a gradient‐recalled echo (GRE) sequence and exploits the principles of fractional encoding. This GRE MRE pulse sequence was validated as capable of acquiring full three‐dimensional motion encoding of isotropic voxels in a large volume within less than a minute. This sequence is suitable for monofrequency and multifrequency MRE experiments. In homogeneous paraffin phantoms, the eXpresso sequence yielded similar storage modulus values as those obtained with conventional methods, although with markedly reduced variances (7.11 ± 0.26 kPa for GRE MRE versus 7.16 ± 1.33 kPa for the conventional spin‐echo EPI sequence). The GRE MRE sequence obtained better phase‐to‐noise ratios than the equivalent spin‐echo EPI sequence (matched for identical acquisition time) in both paraffin phantoms and in vivo data in the liver (59.62 ± 11.89 versus 27.86 ± 3.81, 61.49 ± 14.16 versus 24.78 ± 2.48 and 58.23 ± 10.39 versus 23.48 ± 2.91 in the X, Y and Z components, respectively, in the case of liver experiments). Phase‐to‐noise ratios were similar between GRE MRE used in monofrequency or multifrequency experiments (75.39 ± 14.93 versus 86.13 ± 18.25 at 28 Hz, 71.52 ± 24.74 versus 86.96 ± 30.53 at 56 Hz and 95.60 ± 36.96 versus 61.35 ± 26.25 at 84Hz, respectively). Copyright

Transperineal prostate MR elastography: initial in vivo results.

This article presents a new approach to magnetic resonance elastography of the prostate using transperineal mechanical excitation. This approach is validated using a prostate elasticity phantom and in vivo studies of healthy volunteers. It is demonstrated that the transperineal approach can generate shear wave amplitudes on the order of 6–30 μm in the mid‐gland region. The driver was implemented using an electromagnetic actuator with a hydraulic transmission system. The magnetic resonance elastography acquisition time has been reduced significantly by using a “second harmonic” approach. Displacement fields are processed using the established three‐dimensional local frequency estimation algorithm. The three‐dimensional curl‐based direct inversion was used to calculate the local wavelength. The traveling wave expansion algorithm was used to reconstruct the wave damping image for one case. Using the proposed method, it was possible to resolve lesions of 0.5 cc in the phantom study. Repeatability experiments were performed and analyzed. The results from this study indicate that transperineal magnetic resonance elastography—without an endorectal coil—is a suitable candidate for a patient study involving multiparametric magnetic resonance imaging of prostate cancer, where magnetic resonance elastography may provide additional information for improved diagnosis and image‐based surveillance. Magn Reson Med, 2013.

Travelling Wave Expansion: A Model Fitting Approach to the Inverse Problem of Elasticity Reconstruction

In this paper, a novel approach to the problem of elasticity reconstruction is introduced. In this approach, the solution of the wave equation is expanded as a sum of waves travelling in different directions sharing a common wave number. In particular, the solutions for the scalar and vector potentials which are related to the dilatational and shear components of the displacement respectively are expanded as sums of travelling waves. This solution is then used as a model and fitted to the measured displacements. The value of the shear wave number which yields the best fit is then used to find the elasticity at each spatial point. The main advantage of this method over direct inversion methods is that, instead of taking the derivatives of noisy measurement data, the derivatives are taken on the analytical model. This improves the results of the inversion. The dilatational and shear components of the displacement can also be computed as a byproduct of the method, without taking any derivatives. Experimental results show the effectiveness of this technique in magnetic resonance elastography. Comparisons are made with other state-of-the-art techniques.

Sparsity regularization in dynamic elastography

We consider the inverse problem of continuum mechanics with the tissue deformation described by a mixed displacement-pressure finite element formulation. The mixed formulation is used to model nearly incompressible materials by simultaneously solving for both elasticity and pressure distributions. To improve numerical conditioning, a common solution to this problem is to use regularization to constrain the solutions of the inverse problem. We present a sparsity regularization technique that uses the discrete cosine transform to transform the elasticity and pressure fields to a sparse domain in which a smaller number of unknowns is required to represent the original field. We evaluate the approach by solving the dynamic elastography problem for synthetic data using such a mixed finite element technique, assuming time harmonic motion, and linear, isotropic and elastic behavior for the tissue. We compare our simulation results to those obtained using the more common Tikhonov regularization. We show that the sparsity regularization is less dependent on boundary conditions, less influenced by noise, requires no parameter tuning and is computationally faster. The algorithm has been tested on magnetic resonance elastography data captured from a CIRS elastography phantom with similar results as the simulation.

Curl-Based Finite Element Reconstruction of the Shear Modulus Without Assuming Local Homogeneity: Time Harmonic Case

In elasticity imaging, the shear modulus is obtained from measured tissue displacement data by solving an inverse problem based on the wave equation describing the tissue motion. In most inversion approaches, the wave equation is simplified using local homogeneity and incompressibility assumptions. This causes a loss of accuracy and therefore imaging artifacts in the resulting elasticity images. In this paper we present a new curl-based finite element method inversion technique that does not rely upon these simplifying assumptions. As done in previous research, we use the curl operator to eliminate the dilatational term in the wave equation, but we do not make the assumption of local homogeneity. We evaluate our approach using simulation data from a virtual tissue phantom assuming time harmonic motion and linear, isotropic, elastic behavior of the tissue. We show that our reconstruction results are superior to those obtained using previous curl-based methods with homogeneity assumption. We also show that with our approach, in the 2-D case, multi-frequency measurements provide better results than single-frequency measurements. Experimental results from magnetic resonance elastography of a CIRS elastography phantom confirm our simulation results and further demonstrate, in a quantitative and repeatable manner, that our method is accurate and robust.

Prostate MR elastography with transperineal electromagnetic actuation and a fast fractionally encoded steady-state gradient echo sequence

Our aim is to develop a clinically viable, fast‐acquisition, prostate MR elastography (MRE) system with transperineal excitation. We developed a new actively shielded electromagnetic transducer, designed to enable quick deployment and positioning within the scanner. The shielding of the transducer was optimized using simulations. We also employed a new rapid pulse sequence that encodes the three‐dimensional displacement field in the prostate gland using a fractionally encoded steady‐state gradient echo sequence, thereby shortening the acquisition time to a clinically acceptable 8–10 min. The methods were tested in two phantoms and seven human subjects (six volunteers and one patient with prostate cancer). The MRE acquisition time for 24 slices, with an isotropic resolution of 2 mm and eight phase offsets, was 8 min, and the total scan, including positioning and set‐up, was performed in 15–20 min. The phantom study demonstrated that the transducer does not interfere with the acquisition process and that it generates displacement amplitudes that exceed 100 µm even at frequencies as high as 300 Hz. In the in vivo human study, average wave amplitudes of 30 µm (46 µm at the apex) were routinely achieved within the prostate gland at 70 Hz. No pain or discomfort was reported. Results in a single patient suggest that MRE can identify cancer tumors, although this result is preliminary. The proposed methods allow the integration of prostate MRE with other multiparametric MRI methods. The results of this study clearly motivate the clinical evaluation of transperineal MRE in patients. Copyright

Registration of Whole-Mount Histology and Volumetric Imaging of the Prostate Using Particle Filtering

Registration of histological slices to volumetric imaging of the prostate is an important task that can be used to optimize imaging for cancer detection. Such registration is challenging due to physical changes of the specimen during excision and fixation, and misalignment of the histological slices during preparation and digital scanning. In this work, we consider a multi-slice to volume registration method in which a stack of sparse, unaligned 2-D whole-mount histological slices is registered to a 3-D volumetric imaging of the prostate. We propose a particle filtering framework to contend with the high dimensionality of the search space and multimodal nature of the optimization. Such framework allows modeling of the uncertainty in the pose of the slices and in the imaged information, in order to derive optimal registration parameters in a Bayesian approach. Intensity-, region-, and point-based similarity metrics were incorporated into the registration algorithm to account for different imaging modalities. We demonstrate and evaluate our method on a diverse set of data that includes a synthetic volume, ex vivo and in vivo magnetic resonance imaging, and in vivo ultrasound.

Permanent magnet desktop magnetic resonance imaging system with microfabricated multiturn gradient coils for microflow imaging in capillary tubes

A prototype for a desktop high-resolution magnetic resonance imaging (MRI) velocimetry instrument to characterize flow fields in a capillary tube is demonstrated. This inexpensive compact system is achieved with a 0.6 T permanent magnetic configuration (Larmor frequency of 25 MHz) and temperature compensation using off-the-shelf NdFeB permanent magnets. A triaxial gradient module with microfabricated copper coils using a lithographic fabrication process has been developed. This gradient module is capable of generating fast-switching gradients (<100 micros) with amplitudes up to 1.7 T/m using custom made current amplifiers, and was optimized for microflow imaging. The radio frequency probe is integrated with the gradient module and is driven by custom electronics. A two-dimensional (2D) cross-sectional static image of the inside of a capillary tube with an inner diameter of 1.67 mm is acquired at an in-plane spatial resolution of better than 40 microm. Time-of-flight flow measurements were also obtained using this MRI system to measure the velocity profile of water flowing at average velocities of above 50 mm/s. The flow profile for slower flow velocities was obtained using phase-encoded techniques, which provides quantitative velocity information in 2D.

Model-based registration of ex vivo and in vivo MRI of the prostate using elastography

Registration of histopathology to in vivo magnetic resonance imaging (MRI) of the prostate is an important task that can be used to optimize in vivo imaging for cancer detection. Such registration is challenging due to the change in volume and deformation of the prostate during excision and fixation. One approach towards this problem involves the use of an ex vivo MRI of the excised prostate specimen, followed by in vivo to ex vivo MRI registration of the prostate. We propose a novel registration method that uses a patient-specific biomechanical model acquired using magnetic resonance elastography to deform the in vivo volume and match it to the surface of the ex vivo specimen. The forces that drive the deformations are derived from a region-based energy, with the elastic potential used for regularization. The incorporation of elastography data into the registration framework allows inhomogeneous elasticity to be assigned to the in vivo volume. We show that such inhomogeneity improves the registration results by providing a physical regularization of the deformation map. The method is demonstrated and evaluated on six clinical cases.

Real-Time quantitative elasticity imaging of deep tissue using free-hand conventional ultrasound

In this article an ultrasound elastography technology is reported which provides quantitative images of tissue elasticity from deep soft tissue. The technique is analogous to Magnetic Resonance Elastography in the use of external mechanical vibrations which can penetrate deep tissue. Multifrequency steady-state mechanical vibrations are applied to the tissue at the skin and tissue displacements are measured by a conventional ultrasound system. Absolute values of tissue elasticity are computed in real-time for each frequency and displayed to the physician. The quantitative elasticity images produced by the technology are validated with magnetic resonance elastography images as the gold standard on standard elasticity phantoms. Preliminary in-vivo data from healthy volunteers are presented which show the potential of the technology for clinical use. The system is currently being used in different clinical studies to image kidney fibrosis, liver fibrosis, and prostate cancer.