2D tumor perfusion numerical phantom for dynamic contrast enhanced photoacoustic tomography virtual imaging studies

This dataset provides anatomically realistic spatiotemporal maps of anatomy and photoacoustically induced pressure distribution for virtual imaging studies of dynamic photoacoustic tomography of small animal models. Physical dimensions of the 2D-plus-time object are 3.72 cm x 3.72 cm and a 648 seconds time horizon (with a breathing cycle of 5s). The object is rendered on a spatiotemporal Cartesian grid with 186 x 186 pixels and 1,296 time frames. The anatomical information of the spatiotemporal object is derived from a two-dimensional (2D) abdominal slice extracted from an anatomically realistic whole-body murine numerical phantom (MOBY, Segars et Al, 2004). MOBY comprises detailed three-dimensional (3D) anatomical structures featuring tens of organs and models physiologically realistic cardiac and respiratory motions. A spherical lesion of diameter 7 mm was inserted in the abdominal region at a depth of 5 mm to mimic a subcutaneous tumor flank. To simulate the dynamic imaging of a tumor perfusion model, the concentration of an optical contrast agent in the lesion and blood-filled organs was computed using a two compartment Toft's perfusion model. In this model, the contrast agent was introduced after 90 seconds and injected over a 30 second window following a slow bolus pattern. Optical properties (optical absorption coefficient and scattering coefficient) corresponding to a 800 nm laser pulse wavelength were then assigned to each tissue type based on physiologically plausible chromophore concentrations (oxygenated hemoglobin, deoxygenated hemoglobin, water, fat, contrast agent). Contrast agent optical properties were chosen to mimic those of indocyanine green (ICG), those absorption spectrum has a peak near the chosen wavelength.. Each organ/tissue was then assigned a plausible value of the optical absorption coefficient in the near-infrared range corresponding to the diffusion of an optical contrast agent. The optical scattering coefficient was then assumed to be constant throughout the object. For each time frame, the photoacoustically induced pressure distribution was simulated by solving the diffusion approximation to the radiative transport equations with uniform illumination. The data set consists of the following files: induced_pressure.mat: Induced pressure map on abdominal slice of murine phantom anatomy.mat: Anatomy map on abdominal slice of murine phantom This dataset was saved as a MATLAB binary file version 5 (extension .mat). It can be imported in MATLAB using the load function. For more information about the MOBY phantom and how to obtain it, please see https://otc.duke.edu/technologies/4d-mouse-whole-body-moby-phantom-version-2-0/