1041

A dynamic digital phantom with realistic vasculature and perfusion based on MR histology
Chengyue Wu1, David A. Hormuth2, Federico Pineda3, Gregory S. Karczmar3, Robert D. Moser2,4, and Thomas E. Yankeelov1,2,5,6
1Department of Biomedical Engineering, University of Texas at Austin, Austin, TX, United States, 2Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, TX, United States, 3Department of Radiology, University of Chicago, Chicago, IL, United States, 4Department of Mechanical Engineering, University of Texas at Austin, Austin, TX, United States, 5Department of Diagnostic Medicine, University of Texas at Austin, Austin, TX, United States, 6Department of Oncology, University of Texas at Austin, Austin, TX, United States

Synopsis

Digital phantoms are valuable tools for developing or optimizing new imaging techniques, devices, and analyses. In this contribution, we seek to develop a dynamic digital phantom which contains a detailed representation of vascular structure, tissue properties, and perfusion based on high-resolution MRI data of a rat kidney (courtesy of the Duke Center for In Vivo Microscopy). This dynamic digital phantom can be used to simulate perfusion and diffusion MRI techniques, and systematically evaluate new magnetic resonance imaging acquisition reconstruction/image processing techniques.

Introduction:

Digital phantoms are important research tools for developing new medical imaging techniques and devices1. By providing the “ground truth” of a particular object of interest and its associated biophysical properties, digital phantoms enable the systematic evaluation and optimization of new imaging techniques. Many previous studies have developed digital phantoms focusing on large scale anatomy, motion, or fixed tissue properties. For example, the 4D XCAT phantom2 models cardiac-torso anatomy and body movements. The VICTRE Phantom3 recapitulates the tissues and large vessels of the breast, to enable in silico replication of existing clinical features, and demonstrates that computational modeling can play a central role in the regulatory assessment of imaging products. However, relatively little work has been done on generating digital phantoms with both realistic anatomical (e.g., microvasculature structure) and functional (e.g., perfusion and interstitial transport) microenvironment characteristics. These features are fundamentally important in the development of realistic simulations of dynamic contrast enhanced (DCE-) MRI and diffusion weighted (DW-) MRI measurements. In this contribution, we present preliminary efforts at developing a digital phantom which contains a detailed representation of vascular structure, tissue properties, and perfusion based on MR microscopy imaging of a murine kidney. This dynamic digital phantom could be a useful tool for evaluating new DCE- and DW-MRI techniques or associated image processing methods.

Methods:

Data:
The digital phantom is based on the MR histology data of excised rat kidneys published by Xie et al.4. The dataset used in this study (available via the Duke Center for In Vivo Microscopy5), consists of contrast-enhanced T1-weighted and T2*-weighted images of a male Sprague-Dawley rat (52 weeks). All images are acquired with a matrix size of 1024 × 1024 × 512 slices, yielding an isotropic resolution of 31-microns. Acquisition details can be found in Xie's paper4.
Generation of the geometric model:
The kidney volume is segmented by applying a k-means clustering algorithm on the T1-weighted images. The edge of the segmented mask is then smoothed with a Gaussian filter with a standard deviation of 1. A semi-automatic algorithm is used to extract the vascular from the T2*-weighted images, after normalizing the global intensity to suppress background noise. The algorithm is initialized with a manually identified starting point within the renal artery, from which vessel centerlines are automatically tracked to determine the local vascular radius and orientation along the main trunk; this provides the first generation of the vascular tree. Manually selected branching points are identified along the main trunk, and the algorithm is repeated along these branches to identify subsequent generations in the image set. The end result is the hierarchical, tree-like structure of vasculature.
Generation of the dynamic model:
We apply an image-based computational fluid dynamics model6 to compute the dynamics of blood flow and tracer propagation on the vascular geometry detected in the previous step. Literature values of vascular permeability7, Lp = 1.0 × 10-9 g-1·cm2·s and diffusivity8, ADC = 2.69 × 10-3 mm2/s are assigned on the vascular and extravascular tissues, respectively. Vascular inlet and outlet pressures9,10 are prescribed as Pin = 105 mmHg and Pout = 15 mmHg, respectively. Our computational fluid dynamics model then returns estimates of 1) steady-state flow fields, including blood flow through the vasculature, and interstitial flow in the extravascular renal tissue, modeled by a 1D-3D coupled fluid system6, and 2) the spatio-temporal evolution of the field of contrast agent concentration, modeled by an advection-diffusion equation.

Results:

Figure 1 shows the vascular geometry of the kidney in the digital phantom. The left panel (a) provides an overall view of the entire vascular tree and kidney. Panels (b) to (g) shows the six individually reconstructed renal arterial trees. In this case, the major branches of renal arteries and three generations of further arterioles are completely constructed based on the MR histology image, preserving the realistic geometry and volume.
The calculated blood and interstitial flow on the same digital phantom are presented in Figure 2. The solution within the vascular domain provides the blood pressure, blood flow rate, and flow direction. With these calculated parameters, we can further compute the bolus-arrival time through the vasculature, with the time of the bolus arrival at the roots of the renal arteries considered as the initial time point; see panels (a) – (c). Within the interstitial space, the interstitial pressure and flow velocity fields are solved; see panels (d) – (f). This calculation provides spatially-resolved estimates of renal arterial perfusion and interstitial transport.
Based on the computed steady-state flow, the delivery of contrast agent over time can be estimated using the advection-diffusion equation. As shown in Figure 3, the tracer propagates through vasculature, and is delivered through the whole kidney tissue via convection and diffusion over time.

Discussion and Conclusions:

We have provided a preliminary report on a digital phantom which provides a detailed characterization of vasculature, realistic tissue properties, and physical flow dynamics. This dynamic digital phantom can be used to simulate perfusion and diffusion MRI techniques, and systematically evaluate new MRI acquisition reconstruction and image processing techniques.

Acknowledgements

NCI U01CA142565, U01CA174706, and R01 CA218700. CPRIT RR160005. T.E.Y. is a CPRIT Scholar in Cancer Research. MR histology images courtesy of Duke Center for In Vivo Microscopy.

References

[1] Hsu CM. Three-dimensional computerized breast phantom based on empirical data. Dissertation Abstracts International. 2010;72(02).

[2] Segars WP, Sturgeon G, Mendonca S, Grimes J, Tsui BM. 4D XCAT phantom for multimodality imaging research. Medical physics. 2010;37(9):4902-15.

[3] Badano A, Badal A, Glick S, Graff CG, Samuelson F, Sharma D, Zeng R. In silico imaging clinical trials for regulatory evaluation: initial considerations for VICTRE, a demonstration study. In Medical Imaging 2017: Physics of Medical Imaging, vol. 10132, p. 1013220.

[4] Xie L, Cianciolo RE, Hulette B, Lee HW, Qi Y, Cofer G, Johnson GA. Magnetic resonance histology of age-related nephropathy in the Sprague Dawley rat. Toxicologic pathology. 2012;40(5):764-78.

[5] Supplemental Material for Magnetic Resonance Histology of Age-related Nephropathy in Sprague-Dawley Rats. http://www.civm.duhs.duke.edu/lx201107/index.html

[6] Wu C, Hormuth DA, Oliver TA, Pineda F, Karczmar GS, Moser RD, Yankeelov TE. Patient-specific characterization of breast tumor-associated flow using image-guided computational fluid dynamics. Abstract #4230. May 2019, ISMRM 27th annual meeting, Montreal, QC, Canada.

[7] Jain RK. Normalizing tumor microenvironment to treat cancer: bench to bedside to biomarkers. Journal of Clinical Oncology. 2013;31(17):2205.

[8] Yang D, Ye Q, Williams DS, Hitchens TK, Ho C. Normal and transplanted rat kidneys: diffusion MR imaging at 7 T. Radiology. 2004;231(3):702-9.

[9] Winston JA, Safirstein RO. Reduced renal blood flow in early cisplatin-induced acute renal failure in the rat. American Journal of Physiology-Renal Physiology. 1985;249(4):F490-6.

[10] Shore AC. Capillaroscopy and the measurement of capillary pressure. British journal of clinical pharmacology. 2000;50(6):501-13.

Figures

Figure 1. Geometric model of the digital phantom. Panel (a) presents the 3D rendering of the whole vasculature, with the kidney shown as semi-translucent volume. Panels (b) – (g) show the individual renal arterial trees, respectively corresponding to the labels (blue) in the panel (a).

Figure 2. Steady state flow of the digital phantom. Panels (a) – (c) present MIPs (maximum intensity projection) of the calculated bolus-arrival time (BAT) in the axial, coronal, and sagittal orientations. The BAT maps indicate the direction of blood flow. Panel (d) presents the calculated interstitial pressure in the extravascular space of the kidney tissue. Panels (e) and (f) show the computed interstitial flow field, where panel (e) emphasizes the distribution of the flow velocity magnitude by color, and the panel (f) visualizes the streamlines to show flow directions.

Figure 3. Dynamic model of the digital phantom. Panels (a) – (e) shows the spatial distribution of contrast agent concentration at 5s, 10s, 15s and 20s, and 25s after the initial delivery of the bolus, respectively.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
1041