Imaging lung perfusion

From the first measurements of the distribution of pulmonary blood flow using radioactive tracers by West and colleagues (J Clin Invest 40: 1-12, 1961) allowing gravitational differences in pulmonary blood flow to be described, the imaging of pulmonary blood flow has made considerable progress. The researcher employing modern imaging techniques now has the choice of several techniques, including magnetic resonance imaging (MRI), computerized tomography (CT), positron emission tomography (PET), and single photon emission computed tomography (SPECT). These techniques differ in several important ways: the resolution of the measurement, the type of contrast or tag used to image flow, and the amount of ionizing radiation associated with each measurement. In addition, the techniques vary in what is actually measured, whether it is capillary perfusion such as with PET and SPECT, or larger vessel information in addition to capillary perfusion such as with MRI and CT. Combined, these issues affect quantification and interpretation of data as well as the type of experiments possible using different techniques. The goal of this review is to give an overview of the techniques most commonly in use for physiological experiments along with the issues unique to each technique.     

A computational model of hemodynamic parameters in cortical capillary networks

The analysis of hemodynamic parameters and functional reactivity of cerebral capillaries is still controversial. To assess the hemodynamic parameters in the cortical capillary network, a generic model was created using 2D voronoi tessellation in which each edge represents a capillary segment. This method is capable of creating an appropriate generic model of cerebral capillary network relating to each part of the brain cortex because the geometric model is able to vary the capillary density. The modeling presented here is based on morphometric parameters extracted from physiological data of the human cortex. The pertinent hemodynamic parameters were obtained by numerical simulation based on effective blood viscosity as a function of hematocrit and microvessel diameter, phase separation and plasma skimming effects. The hemodynamic parameters of capillary networks with two different densities (consistent with the variation of the morphometric data in the human cortical capillary network) were analyzed. The results show pertinent hemodynamic parameters for each model. The heterogeneity (coefficient variation) and the mean value of hematocrits, flow rates and velocities of the both network models were specified. The distributions of blood flow throughout the both models seem to confirm the hypothesis in which all capillaries in a cortical network are recruited at rest (normal condition). The results also demonstrate a discrepancy of the network resistance between two models, which are derived from the difference in the number density of capillary segments between the models.     

Inert gas clearance from tissue by co-currently and counter-currently arranged microvessels

To elucidate the clearance of dissolved inert gas from tissues, we have developed numerical models of gas transport in a cylindrical block of tissue supplied by one or two capillaries. With two capillaries, attention is given to the effects of co-current and counter-current flow on tissue gas clearance. Clearance by counter-current flow is compared with clearance by a single capillary or by two co-currently arranged capillaries. Effects of the blood velocity, solubility, and diffusivity of the gas in the tissue are investigated using parameters with physiological values. It is found that under the conditions investigated, almost identical clearances are achieved by a single capillary as by a co-current pair when the total flow per tissue volume in each unit is the same (i.e., flow velocity in the single capillary is twice that in each co-current vessel). For both co-current and counter-current arrangements, approximate linear relations exist between the tissue gas clearance rate and tissue blood perfusion rate. However, the counter-current arrangement of capillaries results in less-efficient clearance of the inert gas from tissues. Furthermore, this difference in efficiency increases at higher blood flow rates. At a given blood flow, the simple conduction-capacitance model, which has been used to estimate tissue blood perfusion rate from inert gas clearance, underestimates gas clearance rates predicted by the numerical models for single vessel or for two vessels with co-current flow. This difference is accounted for in discussion, which also considers the choice of parameters and possible effects of microvascular architecture on the interpretation of tissue inert gas clearance.     

A mixture theory model of fluid and solute transport in the microvasculature of normal and malignant tissues. I. Theory

In order to better understand the mechanisms governing transport of drugs, nanoparticle-based treatments, and therapeutic biomolecules, and the role of the various physiological parameters, a number of mathematical models have previously been proposed. The limitations of the existing transport models indicate the need for a comprehensive model that includes transport in the vessel lumen, the vessel wall, and the interstitial space and considers the effects of the solute concentration on fluid flow. In this study, a general model to describe the transient distribution of fluid and multiple solutes at the microvascular level was developed using mixture theory. The model captures the experimentally observed dependence of the hydraulic permeability coefficient of the capillary wall on the concentration of solutes present in the capillary wall and the surrounding tissue. Additionally, the model demonstrates that transport phenomena across the capillary wall and in the interstitium are related to the solute concentration as well as the hydrostatic pressure. The model is used in a companion paper to examine fluid and solute transport for the simplified case of an axisymmetric geometry with no solid deformation or interconversion of mass.     

Mathematical model for the estimation of hemodynamic and oxygenation variables by tissue spectroscopy

This article presents a quasistatic, compartmental model of tissue-level hemodynamics and oxygenation that leads to a set of formulas, which is suitable to calculate important physiological variables from the mean tissue concentration and saturation of hemoglobin, measured by tissue spectroscopy. Dimensioned quantities are represented relative to their baseline value in the equations (relative value = perturbed/baseline). All model parameters are non-dimensional. The model is based and extends on a number of previous works: previous models of similar aim and scope are consolidated, and every critical assumptions and approximations are treated explicitly; extensions include for example the incorporation of the Fahraeus-effect and the separate estimation of the volume changes of the arterial and the venous compartments. The information content of spectroscopic data alone is shown to be valuable, but limited: the relative venous volume, the oxygen extraction fraction and the relative cellulovascular coupling (defined as the ratio of blood flow and oxygen consumption) can be calculated from these data, if the alterations in arterial blood volume are negligible. The number of variables estimated by the derived formulas can be increased if local blood flow is measured simultaneously: in this case, the relative arterial and venous volume and resistance, the oxygen extraction fraction, and the relative oxygen consumption can be determined. Given that this model considers arterial blood pressure, saturation and hematocrit as its inputs, when measured, the model becomes applicable in such conditions as hyper- or hypotension, hypoxic hypoxia, hemodilution and hemorrhage, where these variables do change. The estimation of the changes in arterial resistance can be applied to estimate the extent of an autoregulatory response.     

A geometrical model of dermal capillary clearance

A new microscopic model is developed to describe the dermal capillary clearance process of skin permeants. The physiological structure is represented in terms of a doubly periodic array of absorbing capillaries. Convection-dominated transport in the blood flow within the capillaries is coupled with interstitial diffusion, the latter process being quantified via a slender-body-theory approach. Convection across the capillary wall and in the interstitial phase is treated as a perturbation which may be added to the diffusive transport. The model accounts for the finite permeability of the capillary wall as well as for the geometry of the capillary array, based on realistic values of physiological parameters. Calculated dermal concentration profiles for permeants having the size and lipophilicity of salicylic acid and glucose illustrate the power and general applicability of the model. Furthermore, validation of the model with published in vivo experimental results pertaining to human skin permeation of hydrocortisone is presented. The model offers the possibility for in-depth theoretical understanding and prediction of subsurface drug distribution in the human skin following topical application, as well as rates of capillary clearance into the systemic circulation. A simpler approach that treats the capillary bed as a homogeneously absorbing zone is also employed. The latter may be used in conjunction with the capillary exchange model to estimate measurable dermal transport and clearance parameters in a straightforward manner.     

Endothelialized Microfluidics for Studying Microvascular Interactions in Hematologic Diseases

Advances in microfabrication techniques have enabled the production of inexpensive and reproducible microfluidic systems for conducting biological and biochemical experiments at the micro- and nanoscales ^1,2. In addition, microfluidics have also been specifically used to quantitatively analyze hematologic and microvascular processes, because of their ability to easily control the dynamic fluidic environment and biological conditions^3-6. As such, researchers have more recently used microfluidic systems to study blood cell deformability, blood cell aggregation, microvascular blood flow, and blood cell-endothelial cell interactions^6-13.However, these microfluidic systems either did not include cultured endothelial cells or were larger than the sizescale relevant to microvascular pathologic processes. A microfluidic platform with cultured endothelial cells that accurately recapitulates the cellular, physical, and hemodynamic environment of the microcirculation is needed to further our understanding of the underlying biophysical pathophysiology of hematologic diseases that involve the microvasculature.Here, we report a method to create an "endothelialized" in vitro model of the microvasculature, using a simple, single mask microfabrication process in conjunction with standard endothelial cell culture techniques, to study pathologic biophysical microvascular interactions that occur in hematologic disease. This "microvasculature-on-a-chip" provides the researcher with a robust assay that tightly controls biological as well as biophysical conditions and is operated using a standard syringe pump and brightfield/fluorescence microscopy. Parameters such as microcirculatory hemodynamic conditions, endothelial cell type, blood cell type(s) and concentration(s), drug/inhibitory concentration etc., can all be easily controlled. As such, our microsystem provides a method to quantitatively investigate disease processes in which microvascular flow is impaired due to alterations in cell adhesion, aggregation, and deformability, a capability unavailable with existing assays.     

Kinetic modeling of [(18)F]FDG in skeletal muscle by PET: a four-compartment five-rate-constant model

Various modeling strategies have been developed to convert regional [(18)F]fluorodeoxyglucose ([(18)F]FDG) concentration measured by positron emission tomography (PET) to a measurement of physiological parameters. However, all the proposed models have been developed and tested mostly for brain studies. The purpose of the present study is to select the most accurate model for describing [(18)F]FDG kinetics in human skeletal muscle. The database consists of basal and hyperinsulinemic-euglycemic studies performed in normal subjects. PET data were first analyzed by an input-output modeling technique (often called spectral analysis). These results provided guidelines for developing a compartmental model. A new model with four compartments and five rate constants (5K model) emerged as the best. By accounting for plasma and extracellular and intracellular kinetics, this model allows, for the first time, PET assessment of the individual steps of [(18)F]FDG kinetics in human skeletal muscle, from plasma to extracellular space to transmembrane transport into the cell to intracellular phosphorylation. Insulin is shown to affect transport and phosphorylation but not extracellular kinetics, with the transport step becoming the main site of control. The 5K model also allows definition of the domain of validity of the classic three-compartment three- or four-rate-constant models. These models are candidates for an investigative tool to quantitatively assess insulin control on individual metabolic steps in human muscle in normal and physiopathological states.     

Positron emission tomography: a clinical and biological research tool

Medical and biological imaging has undergone a revolution in the past decade. Positron emission tomography (PET) has been developed to visualize biochemical and physiological phenomena in living humans and animals. For instance, blood flow, blood volume, glucose metabolism, amino acid metabolism, can be quantitatively estimated by means of PET with various radioactive tracers. This functional and molecular imaging technique has progressed rapidly from being a research technique in laboratories to a routine clinical imaging modality. The most widely used radiotracer in routine is 18F-fluorodeoxyglucose (18FDG), which is an analogue of glucose. Since glucose metabolism is increased many fold in malignant tumors, PET has a major role in the field of clinical oncology and recently in cardiology and neurology. PET is also a valuable tool to study cerebral or cardiac binding sites and to image the expression of reporter genes in small animals. In this review, we summarize the most recent developments in PET imaging with particular reference to the radiotracers available and their application.     

The identification of transport parameters from truncated physiological tracer curves

The primary data curves from the injection of tracers into the circulation are often obscured by the appearance of a recirculation hump. Previous techniques for parameter identification in the presence of such a complication have been extended from compartmental models to include circulatory models consisting of partial differential equations of the Taylor dispersion and Turner capacitance types. A comparative analysis of the efficiency of parameter identification by several computational strategies is presented.     

Tracer kinetic model of regional pulmonary function using positron emission tomography.

To determine the spatial distributions of pulmonary perfusion, shunt, and ventilation, we developed a compartmental model of regional (13)N-labeled molecular nitrogen ((13)NN) kinetics measured from positron emission tomography (PET) images. The model features a compartment for right heart and pulmonary vasculature and two compartments for each region of interest: 1) aerated alveolar units and 2) alveolar units with no gas content (shunting). The model was tested on PET data from normal animals (dogs and sheep) and from animals with experimentally injured lungs simulating acute respiratory distress syndrome. The analysis yielded estimates of regional perfusion, shunt fraction, and specific ventilation with excellent goodness-of-fit to the data (R(2) > 0.99). Model parameters were estimated to within 10% accuracy in the presence of exaggerated levels of experimental noise by using a Monte Carlo sensitivity analysis. Main advantages of the present model are that 1) it separates intraregional blood flow to aerated alveolar units from that shunting across nonaerated units and 2) it accounts and corrects for intraregional tracer removal by shunting blood when estimating ventilation from subsequent washout of tracer. The model was thus found to provide estimates of regional parameters of pulmonary function in sizes of lung regions that could potentially approach the intrinsic resolution for PET images of (13)NN in lung (approximately 7.0 mm for a multiring PET camera).     

Blood pressure and blood flow variation during postural change from sitting to standing: model development and validation.

Short-term cardiovascular responses to postural change from sitting to standing involve complex interactions between the autonomic nervous system, which regulates blood pressure, and cerebral autoregulation, which maintains cerebral perfusion. We present a mathematical model that can predict dynamic changes in beat-to-beat arterial blood pressure and middle cerebral artery blood flow velocity during postural change from sitting to standing. Our cardiovascular model utilizes 11 compartments to describe blood pressure, blood flow, compliance, and resistance in the heart and systemic circulation. To include dynamics due to the pulsatile nature of blood pressure and blood flow, resistances in the large systemic arteries are modeled using nonlinear functions of pressure. A physiologically based submodel is used to describe effects of gravity on venous blood pooling during postural change. Two types of control mechanisms are included: 1) autonomic regulation mediated by sympathetic and parasympathetic responses, which affect heart rate, cardiac contractility, resistance, and compliance, and 2) autoregulation mediated by responses to local changes in myogenic tone, metabolic demand, and CO(2) concentration, which affect cerebrovascular resistance. Finally, we formulate an inverse least-squares problem to estimate parameters and demonstrate that our mathematical model is in agreement with physiological data from a young subject during postural change from sitting to standing.     

A dynamic model of renal blood flow autoregulation

To test whether a mathematical model combining dynamic models of the tubuloglomerular feedback (TGF) mechanism and the myogenic mechanism was sufficient to explain dynamic autoregulation of renal blood flow, we compared model simulations with experimental data. To assess the dynamic characteristics of renal autoregulation, a broad band perturbation of the arterial pressure was employed in both the simulations and the experiments. Renal blood flow and tubular pressure were used as response variables in the comparison. To better approximate the situationin vivo where as large number of individual nephrons act in parallel, each simulation was performed with 125 parallel versions of the model. The key parameters of the 125 versions of the model were chosen randomly within the physiological range. None of the constituent models, i.e., the TGF and the myogenic, could alone reproduce the experimental observations. However, in combination they reproduced most of the features of the various transfer functions calculated from the experimental data. The major discrepancy was the presence of a bimodal distribution of the admittance phase in the simulations. This is not consistent with most of the experimental data, which shows a unimodal curve for the admittance phase. The ability of the model to reproduce the experimental data supports the hypothesis that dynamic autoregulation of renal blood flow is due to the combined action of TGF and the myogenic response.     

Extension of the theory of the multiple-indicator dilution technique to metabolic systems with an arbitrary number of rate constants

The multiple indicator dilution technique consists in the instantaneous injection of a mixture of tracers into the arterial perfusate flow of a catheterized or isolated perfused organ, followed by the analysis of the effluent perfusate. The theory of this technique, which has hitherto been developed for cases where metabolism of a tracer is confined to sequestration described by a single rate constant, is extended in this paper to include an arbitrary number of metabolic rate constants. Partial differential equations with constant coefficients describing the events in a single capillary are derived by applying conventional compartmental analysis to infinitesimally small sections of the capillary. Methods for solving such systems in the time as well as in the frequency domain are developed. From the solutions, the impulse response of the whole organ is evaluated assuming variable capillary and uniform large-vessel transit times. In addition, an efficient method using much less computer time was developed, based on the approximation of the distribution of the capillary transit times by a sum of exponentials. Evaluation of moments (recoveries and mean transit times) is also treated. The results are applied to an example from hepatic lactate metabolism.     

Modelling the impact of clot fragmentation on the microcirculation after thrombectomy

Many ischaemic stroke patients who have a mechanical removal of their clot (thrombectomy) do not get reperfusion of tissue despite the thrombus being removed. One hypothesis for this ‘no-reperfusion’ phenomenon is micro-emboli fragmenting off the large clot during thrombectomy and occluding smaller blood vessels downstream of the clot location. This is impossible to observe in-vivo and so we here develop an in-silico model based on in-vitro experiments to model the effect of micro-emboli on brain tissue. Through in-vitro experiments we obtain, under a variety of clot consistencies and thrombectomy techniques, micro-emboli distributions post-thrombectomy. Blood flow through the microcirculation is modelled for statistically accurate voxels of brain microvasculature including penetrating arterioles and capillary beds. A novel micro-emboli algorithm, informed by the experimental data, is used to simulate the impact of micro-emboli successively entering the penetrating arterioles and the capillary bed. Scaled-up blood flow parameters–permeability and coupling coefficients–are calculated under various conditions. We find that capillary beds are more susceptible to occlusions than the penetrating arterioles with a 4x greater drop in permeability per volume of vessel occluded. Individual microvascular geometries determine robustness to micro-emboli. Hard clot fragmentation leads to larger micro-emboli and larger drops in blood flow for a given number of micro-emboli. Thrombectomy technique has a large impact on clot fragmentation and hence occlusions in the microvasculature. As such, in-silico modelling of mechanical thrombectomy predicts that clot specific factors, interventional technique, and microvascular geometry strongly influence reperfusion of the brain. Micro-emboli are likely contributory to the phenomenon of no-reperfusion following successful removal of a major clot.     

Metabolic dynamics in skeletal muscle during acute reduction in blood flow and oxygen supply to mitochondria: in-silico studies using a multi-scale, top-down integrated model

Control mechanisms of cellular metabolism and energetics in skeletal muscle that may become evident in response to physiological stresses such as reduction in blood flow and oxygen supply to mitochondria can be quantitatively understood using a multi-scale computational model. The analysis of dynamic responses from such a model can provide insights into mechanisms of metabolic regulation that may not be evident from experimental studies. For the purpose, a physiologically-based, multi-scale computational model of skeletal muscle cellular metabolism and energetics was developed to describe dynamic responses of key chemical species and reaction fluxes to muscle ischemia. The model, which incorporates key transport and metabolic processes and subcellular compartmentalization, is based on dynamic mass balances of 30 chemical species in both capillary blood and tissue cells (cytosol and mitochondria) domains. The reaction fluxes in cytosol and mitochondria are expressed in terms of a general phenomenological Michaelis-Menten equation involving the compartmentalized energy controller ratios ATP/ADP and NADH/NAD(+). The unknown transport and reaction parameters in the model are estimated simultaneously by minimizing the differences between available in vivo experimental data on muscle ischemia and corresponding model outputs in coupled with the resting linear flux balance constraints using a robust, nonlinear, constrained-based, reduced gradient optimization algorithm. With the optimal parameter values, the model is able to simulate dynamic responses to reduced blood flow and oxygen supply to mitochondria associated with muscle ischemia of several key metabolite concentrations and metabolic fluxes in the subcellular cytosolic and mitochondrial compartments, some that can be measured and others that can not be measured with the current experimental techniques. The model can be applied to test complex hypotheses involving dynamic regulation of cellular metabolism and energetics in skeletal muscle during physiological stresses such as ischemia, hypoxia, and exercise.     

Effect of sepsis on skeletal muscle oxygen consumption and tissue oxygenation: interpreting capillary oxygen transport data using a mathematical model

Inherent in the inflammatory response to sepsis is abnormal microvascular perfusion. Maldistribution of capillary red blood cell (RBC) flow in rat skeletal muscle has been characterized by increased 1) stopped-flow capillaries, 2) capillary oxygen extraction, and 3) ratio of fast-flow to normal-flow capillaries. On the basis of experimental data for functional capillary density (FCD), RBC velocity, and hemoglobin O2 saturation during sepsis, a mathematical model was used to calculate tissue O2 consumption (Vo2), tissue Po2 (Pt) profiles, and O2 delivery by fast-flow capillaries, which could not be measured experimentally. The model describes coupled capillary and tissue O2 transport using realistic blood and tissue biophysics and three-dimensional arrays of heterogeneously spaced capillaries and was solved numerically using a previously validated scheme. While total blood flow was maintained, capillary flow distribution was varied from 60/30/10% (normal/fast/stopped) in control to 33/33/33% (normal/fast/stopped) in average sepsis (AS) and 25/25/50% (normal/fast/stopped) in extreme sepsis (ES). Simulations found approximately two- and fourfold increases in tissue Vo2 in AS and ES, respectively. Average (minimum) Pt decreased from 43 (40) mmHg in control to 34 (27) and 26 (15) mmHg in AS and ES, respectively, and clustering fast-flow capillaries (increased flow heterogeneity) reduced minimum Pt to 14.5 mmHg. Thus, although fast capillaries prevented tissue dysoxia, they did not prevent increased hypoxia as the degree of microvascular injury increased. The model predicts that decreased FCD, increased fast flow, and increased Vo2 in sepsis expose skeletal muscle to significant regions of hypoxia, which could affect local cellular and organ function.     

Analytical model for long-distance tracer-transport in plants

Recent investigations of long-distance transport in plants using non-invasive tracer techniques such as ¹¹C radiolabeling monitored by positron emission tomography (PET) combined with magnetic resonance imaging (MRI) revealed the need of dedicated methods to allow a quantitative data analysis and comparison of such experiments. A mechanistic compartmental tracer transport model is presented, defined by a linear system of partial differential equations (PDEs). This model simplifies the complexity of axial transport and lateral exchanges in the transport pathways of plants (e.g. the phloem) by simulating transport and reversible exchange within three compartments using just a few parameters which are considered to be constant in space and time. For this system of PDEs an analytical solution in Fourier-space was found allowing a fast and numerically precise evaluation. From the steady-state behavior of the model, the system loss (steadily fixed tracer along the transport conduits) was derived as an additional parameter that can be readily interpreted in a physiological way. The presented framework allows the model to be fitted to spatio-temporal tracer profiles including error and sensitivity analysis of the estimated parameters. This is demonstrated for PET data sets obtained from radish, sugar beet and maize plants.     

Kinetics of glucose transport and sequestration in lactating bovine mammary glands measured in vivo with a paired indicator/nutrient dilution technique.

To quantify kinetics of mammary glucose utilization in vivo, 24 paired glucose and extracellular indicator (p-aminohippuric acid) dilution curves across intact bovine mammary glands were obtained after bolus injections into the external iliac artery. Dilution curves were analyzed using a compartmental capillary, convolution integration model. Four candidate submodels of glucose transport and metabolism in capillary supply zones were fit to the glucose dilution curves and evaluated. Model I, with one extracellular compartment for glucose and first-order unidirectional uptake, failed, indicating that efflux of glucose from the intracellular space could not be ignored. Model II, with first-order exchanges between extracellular and intracellular compartments and sequestration from the latter, was overdefined because unidirectional clearance of glucose was at least five times the blood flow rate and 20 times the net clearance rate. Model III, combining extracellular and intracellular space into one compartment, was superior in its goodness-of-fit to curves and identifiability of parameters. Michaelis-Menten parameters of sequestration were not identifiable. Parameters of the optimal compartmental capillary, convolution integration model were applicable to both the dynamics of injected glucose dilution and the steady-state background arteriovenous difference of glucose. Glucose sequestration followed first-order kinetics between 0 and 7 mM extracellular glucose with an average rate constant of 0.006 s(-1) or a clearance of 44 ml/s. The ratio of intracellular to extracellular glucose distribution space was 0.34, which is considerably lower than the expected intracellular volume and suggests an intracellular occlusion compartment with which extracellular glucose rapidly exchanges.     

Mapping 3-D functional capillary geometry in rat skeletal muscle in vivo

We have developed a novel mapping software package to reconstruct microvascular networks in three dimensions (3-D) from in vivo video images for use in blood flow and O2 transport modeling. An intravital optical imaging system was used to collect video sequences of blood flow in microvessels at different depths in the tissue. Functional images of vessels were produced from the video sequences and were processed using automated edge tracking software to yield location and geometry data for construction of the 3-D network. The same video sequences were analyzed for hemodynamic and O2 saturation data from individual capillaries in the network. Simple user-driven commands allowed the connection of vessel segments at bifurcations, and semiautomated registration enabled the tracking of vessels across multiple focal planes and fields of view. The reconstructed networks can be rotated and manipulated in 3-D to verify vessel connections and continuity. Hemodynamic and O2 saturation measurements made in vivo can be indexed to corresponding vessels and visualized using colorized maps of the vascular geometry. Vessels in each reconstruction are saved as text-based files that can be easily imported into flow or O2 transport models with complete geometry, hemodynamic, and O2 transport conditions. The results of digital morphometric analysis of seven microvascular networks showed mean capillary diameters and overall capillary density consistent with previous findings using histology and corrosion cast techniques. The described mapping software is a valuable tool for the quantification of in vivo microvascular geometry, hemodynamics, and oxygenation, thus providing rich data sets for experiment-based computational models.