Precise Non-invasive Liver Cirrhosis Quantification via Multimodal Stiffness-Structural Fusion Imaging

Summary

Traditional Magnetic Resonance Elastography separates stiffness and structural images, requiring manual ROI delineation and radiologist input. This study introduces a multimodal fusion method that combines both image types to produce stiffness-structural fusion images and liver stiffness distribution statistics, supporting accurate cirrhosis staging and enabling radiomics analysis based on liver stiffness distribution.

Abstract

Traditional Magnetic Resonance Elastography (MRE) presents stiffness and structural images separately, requiring manual ROI delineation for liver stiffness assessment. This introduces observer variability and may overlook significant lesions, failing to accurately represent whole-liver stiffness distribution patterns. The present study implemented multimodal image fusion integrating stiffness maps with structural images using: (1) acquisition on a 3.0T MRI with optimized T1-weighted imaging (TR/TE = 157.62/1.17 ms) and MRE (60 Hz vibrations); (2) computational preprocessing of DICOM data into 3D volumes; (3) precise co-registration and fusion calculation; and (4) quantitative liver stiffness distribution analysis across different cirrhosis stages. The stiffness-structural fusion imaging technique integrated anatomical structures with stiffness patterns in a single visualization, enabling simultaneous assessment of liver morphology and stiffness distribution. Comparative analysis revealed a clear differentiation between cirrhotic patients and healthy controls, with cirrhotic livers showing significantly higher proportions of advanced fibrosis. This technique eliminated variability due to manual ROI placement and established a standardized framework for comprehensive whole-liver assessment. The stiffness-structural fusion imaging technique provides comprehensive visualization of liver stiffness with precise anatomical correlation, potentially improving diagnostic accuracy compared to conventional MRE. This approach eliminates observer variability, enables standardized whole-liver assessment, and establishes a quantitative framework for monitoring disease progression and treatment response, while facilitating improved clinician-patient communication and offering a foundation for advanced radiomics in liver cirrhosis staging.

Introduction

Liver cirrhosis represents a significant global health challenge, with accurate staging being crucial for optimal patient management and treatment planning^1,2. Magnetic Resonance Elastography (MRE) has emerged as a valuable non-invasive technique for liver cirrhosis assessment, offering advantages over invasive biopsy procedures that carry inherent sampling error risks and patient discomfort^3,4,5. However, conventional MRE approaches present significant limitations in clinical practice that impact diagnostic accuracy and workflow efficiency.

Traditional MRE methodologies typically present stiffness maps and structural images as separate entities, requiring clinicians to mentally correlate anatomical features with stiffness patterns across different image series^5,6. This approach necessitates manual region of interest (ROI) placement for liver stiffness quantification, introducing considerable inter-observer variability and potential sampling bias^7,8. Furthermore, the conventional ROI-based approach fails to capture the heterogeneous nature of liver stiffness distribution, as highlighted in recent literature emphasizing the prognostic significance of parenchymal heterogeneity in chronic liver disease progression^9,10. The inability to comprehensively visualize and quantify whole-liver stiffness distribution represents a significant barrier to advanced radiomics analysis and precise disease characterization.

Alternative non-invasive techniques for liver cirrhosis assessment include ultrasound elastography and serum biomarker panels^11,12. While transient elastography offers point measurements of liver stiffness, it lacks the comprehensive spatial mapping capabilities of MRE ¹³. Serum biomarker panels provide indirect estimates of fibrosis but cannot visualize the spatial distribution of parenchymal changes¹¹. Recent comparative studies have consistently demonstrated superior diagnostic accuracy of MRE for intermediate and advanced fibrosis stages compared to these alternatives^13,14.

To address these limitations, we have developed a novel multimodal stiffness-structural fusion imaging technique that seamlessly integrates quantitative stiffness maps with high-resolution structural images. This approach enables simultaneous visualization of anatomical landmarks and stiffness patterns in a single integrated representation, eliminating the need for mental co-registration and reducing cognitive burden during interpretation^15,16. The technique provides comprehensive whole-liver stiffness distribution analysis, overcoming sampling limitations of conventional ROI-based methods.

This study utilized 3.0 Tesla MRI systems with elastography capabilities, standard patient preparation including 4-6 h fasting, and MATLAB with Image Processing Toolbox for computational analysis. These requirements are increasingly available at academic medical centers and specialized imaging facilities.

This protocol provides a detailed methodology for implementing the stiffness-structural fusion technique in clinical and research settings. The approach is particularly valuable for centers with access to MRE capabilities seeking to enhance diagnostic precision, improve clinical workflow efficiency, and establish quantitative frameworks for longitudinal monitoring of chronic liver disease. The technique is applicable across various etiologies of chronic liver disease and can be implemented with minimal additional post-processing time using the provided computational framework.

Protocol

The protocol was performed under institutional IRB approval. A patient from Dongzhimen Hospital, Beijing University of Chinese Medicine, underwent routine upper abdominal magnetic resonance imaging and MRE at YouAn Hospital with informed consent. This case was selected for the research methodology due to the confirmed progression from fibrosis to cirrhosis, highlighting the value of the approach in diagnosing the boundary between liver fibrosis and cirrhosis. The study also provides a quantitative comparison between this patient's liver stiffness distribution (LSD) and that of a healthy liver. All equipment and software tools used in this study are listed in the Table of Materials.

1. MRI scanning strategies

1. Prepare the patient by ensuring a fasting period of 4-6 h prior to examination to reduce physiological variability.
2. Position the patient on the 3.0 Tesla MRI system using a combination of a 16-channel abdominal array and a 32-channel posterior array coil.
3. Structural imaging acquisition
   1. Perform anatomical imaging using a breath-hold axial 2D T1-weighted dual-echo gradient-recalled echo sequence with the following parameters:
      NOTE: TR/TE = 157.62/1.17 ms; Flip angle = 60°; Field of view = 400 × 400 mm²; Acquisition matrix = 512 × 512; Slice thickness = 7 mm.
   2. Verify that the sequence provides high-resolution structural information (in-plane resolution of 0.78 × 0.78 mm²) with inherent fat-water separation capability.
4. MR Elastography acquisition
   1. Place the pneumatic driver over the right anterior liver to apply mechanical vibrations at 60 Hz.
   2. Perform liver stiffness mapping using a 2D echo-planar imaging (EPI)-based acquisition with the following parameters:
      NOTE: TR/TE = 1050.3/63.4 ms; Flip angle = 90°; Field of view = 420 × 420 mm²; Acquisition matrix = 256 × 256 (reconstructed to 256 × 256); Slice thickness = 10 mm; Bandwidth = 1953.1 Hz/pixel.
   3. Acquire eight phase offsets to compute wave propagation and generate quantitative stiffness maps (elastograms) with values expressed in kilopascals (kPa).
5. Image co-registration strategy
   1. Ensure accurate spatial alignment between structural images and elastograms by maintaining identical patient positioning and consistent breath-hold instructions.
   2. Perform all imaging in the axial plane with matching slice coverage of the liver.
   3. Plan the MRE sequence using the anatomical T1-weighted images as a reference to ensure optimal coverage of the liver parenchyma.
   4. Verify that both sequences utilize the same frame of reference to facilitate automated co-registration in the subsequent fusion process.

2. Data collection and preparation

1. Renaming the folder of every sequence
   NOTE: Add explicit names for each sequence, as DICOM data exported from the equipment lacks sequence identification. This facilitates subsequent analysis and processing.
   1. Execute the Description_Name function to add descriptive names to the folders for each sequence.
   2. Copy all DICOM data to a customized working directory.
   3. Navigate to the directory containing the data in MATLAB's working directory.
2. Quickly checking images of structural sequence
   NOTE: Examine structural images derived from T1 imaging, which share the same coordinate system with stiffness maps and provide clear liver contour visualization essential for subsequent stiffness distribution calculations.
   1. Generate a 3D volume matrix from the DICOM files of the MRI sequences using the following MATLAB code:
      f=dir('*.dcm');
      for i=1:length(f)
      V(:,:,i)=dicomread(fullfile(f(i).folder,f(i).name));
      ​End
   2. Visualize the image sequence using MATLAB's sliceViewer function for interactive examination:
      figure;
      H=sliceViewer(V);colormap(gray(1024));
      ​set(gcf, 'Toolbar', 'figure');
   3. Interact with the sliceViewer GUI. Use the scroll bar to browse through different slices. Verify that liver boundaries are clearly delineated with sharp contrast against surrounding tissues.
      NOTE: Acceptable images should show continuous liver parenchyma without motion artifacts or blurring (Figure 1). This interactive viewer allows a quick assessment of image quality and liver morphology.
   4. Find the icons for zooming in, zooming out, and returning to the global view in the upper-right corner of the GUI. Use the Data Tip icon to mark the coordinates of specific liver regions for further analysis.
   5. Observe that the default gray colormap displays values from low (black) to high (white). Right-click on the Color Bar in the pop-up menu to select different colormap options for optimal liver tissue visualization.
   6. If the contrast is not satisfactory, use the left mouse button to drag up and down in the middle of the figure to adjust the window level. Drag left and right to adjust the window width. The corresponding accurate contrast range will be displayed on the color bar.
      NOTE: Use these interactive controls to inspect MRI structural sequence characteristics across intensity and spatial dimensions. Assess liver anatomical contour integrity and image quality before proceeding to the next step.
3. Quickly checking MRE stiffness maps
   NOTE: Examine MRE stiffness maps for quantitative liver elasticity measurements, where different colors represent varying degrees of stiffness, which is essential for identifying fibrotic or cirrhotic areas.
   1. Generate a 3D volume matrix from the DICOM files of the MRE sequences using the following MATLAB code:
      f=dir('*.dcm');
      for i=1:length(f)
      V_LSD(:,:,i)=dicomread(fullfile(f(i).folder,f(i).name));
      ​End
   2. Visualize the MRE stiffness maps using MATLAB's sliceViewer function for interactive examination:
      figure;H=sliceViewer(V_LSD);colormap(jet(1024));set(gcf, 'Toolbar', 'figure');
   3. Interact with the sliceViewer GUI. Use the scroll bar at the bottom of the graphical user interface (GUI) to browse through different slices in the sequence (Figure 2). Note that this interactive viewer allows a quick assessment of stiffness distribution and liver elasticity patterns.
   4. Find the icons for zooming in, zooming out, and returning to the global view in the upper-right corner of the GUI. Use the Data Tip icon to mark the coordinates of specific liver regions with notable stiffness values for further analysis.
   5. Observe that the default jet colormap represents stiffness values from low (blue) to high (red). Right-click on the Color Bar in the pop-up menu to adjust the colormap range for optimal visualization of liver stiffness gradients.
   6. If the contrast is not satisfactory, use the left mouse button to drag up and down in the middle of the figure to adjust the window level. Drag left and right to adjust the window width. The corresponding accurate stiffness range will be displayed on the color bar.
   7. Confirm that stiffness maps display smooth color transitions within the liver parenchyma. Acceptable stiffness mapping should show minimal noise artifacts and clear differentiation between liver tissue (typically in the 1-8 kPa range) and background regions.
      NOTE: These interactive controls enable flexible inspection of MRE stiffness maps across both intensity space and sequence location, facilitating quick assessment of liver stiffness distribution and potential areas of fibrosis.

3. Stiffness-structural fusion

NOTE: Integrate both modalities since stiffness maps alone cannot accurately localize lesions without structural image correlation.

1. Define path variables for both structural images and stiffness maps, as well as scan files at approximately the same axial positions.
2. Select the Struct_Stiffness_Fusion function (I_fusion = Struct_Stiffness_Fusion(I_structural, I_stiffness)) to generate Stiffness-Structural Fusion Images in the current directory.
3. Locate the fusion image files in the file browser and open them using the default image viewer. Verify that anatomical structures and stiffness patterns are properly aligned.
   NOTE: Acceptable fusion images should demonstrate precise overlay without spatial misregistration between structural and elastography data. Following the steps above, aligned scan axial positions automatically output the fusion image sequence, which can be viewed using common RGB image viewers.

4. The liver stiffness distribution

NOTE: To complement the fusion image analysis shown in Figure 3, generate liver stiffness distribution statistics as demonstrated in Figure 4.

1. Repeat steps 1-3 to examine patients with different stages of liver cirrhosis.
2. Select the LSD_compute function ([LSD_distribution, LSD_stats] = LSD_compute(V_stiffness, liver_mask)) to calculate the Patient's LSD (Liver Stiffness Distribution) values across different stiffness stages¹⁵ (Figure 4).
   NOTE: The patient shown in Figure 3 is in late-stage liver fibrosis that has already progressed to cirrhosis, with the whole liver average stiffness reaching 5 kPa. Therefore, in the quantitative distribution in Figure 4, the majority of liver tissue stiffness values are above the F4 stage.

5. Final verification and data archiving

1. Verify that all fusion images have been successfully generated and saved in the working directory.
2. Confirm that LSD analysis results have been computed and documented for clinical interpretation.
3. Archive all processed data, including structural images, stiffness maps, fusion images, and quantitative LSD metrics, for future reference and longitudinal comparison.
   NOTE: The multimodal stiffness-structural fusion imaging protocol is now complete. The generated fusion images and quantitative LSD analysis provide comprehensive liver assessment data ready for clinical interpretation and patient communication.

Results

The visualization approach presented in this study offers comprehensive insights into liver structure and stiffness distribution, enabling precise assessment of liver cirrhosis. Through the GUI depicted in Figure 1, clinicians can better understand the liver's structural anatomy across multiple cross-sections, primarily evaluating its anatomical integrity. Examining Figure 2 reveals a critical limitation of isolated stiffness maps: without corresponding structur...

Discussion

Critical steps
Several critical steps in the protocol require particular attention to ensure the successful implementation of the multimodal stiffness-structural fusion imaging technique. First, the planning of the MRE sequence using anatomical T1-weighted images as reference ensures optimal coverage of liver parenchyma and facilitates subsequent automated co-registration¹⁵. Second, the preprocessing of DICOM data requires careful management of sequence identification. Since...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
MATLAB	GE Medical Systems	SIGNA Pioneer	3.0 Tesla
MRI system	Intelligent Entropy	LSD-Fusion V1.0	Beijing Intelligent Entropy Science & Technology Co Ltd. Modeling for CT/MRI fusion
Tools for Modeling	MathWorks	2022B	Computing and visualization