3D numerical breast phantom: Normalization of optical fluence distribution for 3D functional OAT

Companion dataset of the manuscript: Seonyeong Park, Frank J. Brooks, Umberto Villa, Richard Su, Mark A. Anastasio, Alexander A. Oraevsky, "Normalization of optical fluence distribution for three-dimensional functional optoacoustic tomography of the breast," J. Biomed. Opt. 27(3) 036001 (16 March 2022) https://doi.org/10.1117/1.JBO.27.3.036001 This dataset contains an anatomically realistic numerical breast phantom (NBP) to validate optical fluence normalization methods for quantitative and functional optoacoustic tomography (OAT). Specifically, it includes 1) 3D maps of optical absorption coefficient; 2) 3D maps of simulated optical fluence distribution, and 3) 3D OAT reconstructed images at three wavelengths in the near-infrared range. The images were reconstructed from noisy synthetic data using filtered back-projection (FBP). The size of the reconstructed image is 480 × 480 × 240 voxels (120 × 120 × 60 mm3). To simulate the synthetic data, a NBP was created using a computational framework for virtual 3D OAT breast imaging trials, developed by the authors of the reference [Park2020]. Further details of the NBP; functional, optical, and acoustic property assignment; and simulation of optoacoustic signals are in the accompanying paper [Park2022]. The following files are included in this dataset: mu_a_w{757, 800, 850}.mat: Optical absorption coefficient distributions of a NBP at illumination wavelengths of 757 nm, 800 nm, and 850 nm; optical_fluence_w{757, 800, 850}.mat: Optical fluence distributions of a NBP that were simulated at illumination wavelengths of 757 nm, 800 nm, and 850 nm using the MCXLAB software [Fang2009], [Yu2018]; and RECON_NOISY_w{757, 800, 850}_FBP.mat: Images reconstructed from noisy synthetic measurements of a NBP, that were simulated at illumination wavelengths of 757 nm, 800 nm, and 850 nm, using FBP. These data were saved as MATLAB binary files version 5 (extension .mat). They can be imported in MATLAB using the load function. The code of optical fluence normalization is available from the GitHub repository: comp-imaging-sci/optical-fluence-normalization_3d-oat-breast.