310 likes | 418 Views
Discover how MR Elastography measures viscoelasticity of soft tissues for medical diagnoses and nondestructive testing. Learn about the MRE system, models, and analytical methods used for data processing and recovery.
E N D
Part 4: Viscoelastic Properties of Soft • Tissues in a Living Body Measured by MR • Elastography Gen Nakamura • Department of Mathematics, Hokkaido University, Japan • (Supported by Japan Science and Technology Agency) • Joint work with Yu, Jiang • ICMAT, Madrid, May 12, 2011
Magnetic Resonance Elastography, MRE • A newly developed non-destructive technique • (Muthupillai et al., Science, 269, 1854-1857, 1995, Mayo Clinic.) • Measure the viscoelasticity of soft tissues in a living body • Diagnosis: • the stage of liver fibrosis • early stage cancer: breast cancer, pancreatic cancer, prostate cancer, etc. • neurological diseases: Alzheimer’s disease, hydrocephalus, multiple sclerosis, etc. • Nondestructive testing (high frequency rheometer): • biological material, polymer material
MRE System in Hokkaido Univ. (JST Proj.) Japan Science and Technology Agency (JST) (1) External vibration system Electromagnetic vibrator (2) Pulse sequence with motion-sensitizing gradients (MSG) GFRP Bar 2~4 m (3) phantom Object External vibration system Micro-MRI Wave image Storage modulus (4) Inversion algorithm
MRE phantom: agarose or PAAm gel 100mm 65mm • --- time harmonic external vibration (3D vector) • --- frequency of external vibration (50~250Hz) • --- amplitude of external vibration (≤ 500 μm ) 70mm 10 mm hard soft
Viscoelastic wave in soft tissues • Time harmonic external vibration • Interior viscoelastic wave • --- amplitude of viscoelastic wave ( : real part, : imaginary part) viscoelastic body after some time
MRE measurements: phase image MRIsignal real part: R imaginary part: I 2D FFT magnitude image phase image MRE measurement
MRE measurements: phase image components in vertical direction (unit: )
Data analysis for MRE viscoelasticity of soft tissues or phantom interior wave displacement Step 1: modeling Step 3: recovery (inverse problem) Step 2: numerical simulation (forward problem) viscoelasticity models for soft tissues or phantom (PDE)
Viscoelasticity models for soft tissues • Time: • : bounded domain; • : Lipschitz continuous boundary; • Displacement: • General linear viscoelasticity model:
Viscoelasticity models for soft tissues • Stress tensor: • Density: • Small deformation (micro meter) ⇒ linear strain tensor
Constitutive equation • Voigt model: • Maxwell model: • Zener model:
Viscoelasticity tensors • full symmetries: • strong convexity (symmetric matrix ):
Time harmonic wave • Boundary: • : open subsets of with , Lipschitz continuous; • Time harmonic boundary input and initial condition: • Time harmonic wave (exponential decay): • Jiang, et. al., submitted to SIAM appl. math.. (isotropic, Voigt) • Rivera, Quar. Appl. Math., 3(4), 629-648, 1994. • Rivera, et. al., Comm. Math. Phys. 177(3), 583–602, 1996.
Time harmonic wave • Stationary model: • Sobolev spaces of fractional order 1/2 or 3/2 • an open subset with a boundary away from and • the set of distributions in the usual fractional Sobolev space compactly supported in • This can be naturally imbedded into
Constitutive equation (stationary case) • Voigt model: • Maxwell model: • Zener model:
Modified Stokes model • Isotropic+ nearly incompressible • Asymptotic analysis ⇒ modified Stokes model: • Jiang et. al., Asymptotic analysis for MRE, submitted to SIAP • H. Ammari et. al., Quar. Appl. Math., 2008: isotopic constant elasticity
Storage modulus andloss modulus • Storage ・ loss modulus () • Voigt model • Maxwell model • Zener • Angular frequency: • Shear modulus: • Shear viscosity: • Measured by rheometer
Modified Stokes model • 2D numerical simulation (Freefem++) • Plane strain assumption mm
Curl operator • Modified Stokes model: • Constants : Curl operator: filter of the pressure term
Pre-treatment: denoising • Mollifier(Murio, D. A.: Mollification and Space Marching) • Smooth function defined in the nbd of • : a bounded domain • : an extension of to • Function : a nonnegative function over such that and
Recovery of storage modulus • Constants: • Mollification: • Curl operator: • Numerical differentiation method • Numerical differentiation is an ill-posed problem • Numerical differentiation with Tikhonov regularization Unstable!!!
Recovery of storage modulus • Constants: • Mollification: • Curl operator: • Numerical Integration Method • : test region (2D or 3D) • : test function Unstable!!!
Recovery of storage modulus • Constants: • Mollification: • Curl operator: • Modified Integral Method • : test region (2D or 3D) • test region size: about one wavelength
Recovery from no noise simulated data Inclusion: small large outside • Exact value: 3.3 kPa 3.3 kPa 7.4 kPa • Mean value: 3.787 kPa 3.768 kPa 7.436 kPa • Stddev: 0.147 0.060 0.003 • Relative error: 0.1476 0.1418 0.00049
Recovery from noisy simulated data • 10% relative error Inclusion: small large outside • Exact value: 3.3 kPa 3.3 kPa 7.4 kPa • Mean value: 4.636 kPa 3.890 kPa 7.422 kPa • StdDev: 0.328 0.129 0.322 • Relative error: 0.4048 0.1788 0.00294
Recovery from experimental data Layered PAAm gel: hard (left) soft (right) Mean value: 31.100 kPa 10.762 kPa StdDev: 0.535 0.201 250 Hz 0.3 mm kPa cm
Recovery of storage modulus G’ Layered PAAm gel: hard (left) soft (right) Mean value: 31.100 (25.974) kPa 10.762 (8.988) kPa Standard deviation: 0.535 (6.982) 0.201 (4.407) modified method (old method (polynomial test function)) 250 Hz 0.3 mm kPa kPa cm cm
Recovery of storage modulus G’ Independent of frequencies (1 ~ 250 Hz) hard soft Rheometer: 32.5456 kPa 9.2472 kPa MRE, 250 Hz: 31.100 kPa 10.762 kPa Relative error: 0.0444 0.1638 Rheometer : ARES-2KFRT, TA Instruments Frequency: 0.1 ~ 10 Hz Strain mode: 5%