3D optoacoustic numerical breast phantoms and simulated OAT measurement data (hemispherical shape, 4 lesions)

Companion dataset of the manuscript: Seonyeong Park, Umberto Villa, Fu Li, Refik Mert Cam, Alexander A. Oraevsky, Mark A. Anastasio, \"Stochastic three-dimensional numerical phantoms to enable computational studies in quantitative optoacoustic computed tomography of breast cancer,\" J. Biomed. Opt. 28(6) 066002 (20 June 2023) https://doi.org/10.1117/1.JBO.28.6.066002 This dataset contains 40 sets of three-dimensional (3D) numerical breast phantoms (NBPs) for use in virtual imaging trials (VITs) of optoacoustic tomography (OAT) and the corresponding simulated multi-wavelength optical fluence distributions, induced initial pressure distributions, and OAT measurement data. The NBPs are in hemispherical shapes, and four different-sized lesions, that are composed solely of a viable tumor cell region, were inserted into each of 40 NBPs. The NBPs correspond to one of the following four breast density types defined in Breast Imaging Reporting and Data System (BI-RAD®): A: Breast is almost entirely fatty; B: Breast has scattered areas of fibroglandular density; C: Breast is heterogeneously dense; D: Breast is extremely dense. Each NBP set consists of A tissue label map (anatomical NBP); Functional properties and chromophore concentrations maps (functional NBPs); Optical properties maps (optical absorption and scattering coefficients, scattering anisotropy, and refractive indexes) at multiple illumination wavelengths (optical NBPs); Acoustic properties (speed of sound, density, and acoustic attenuation) maps (acoustic NBPs); Simulated multi-wavelength optical fluence distributions; Simulated multi-wavelength induced initial pressure distributions; Simulated multi-wavelength acoustic measurements. The tissue label map (anatomical NBP) was created by use of our adaptation of tools from the Virtual Imaging Clinical Trials for Regulatory Evaluation (VICTRE) project at the U.S. Food and Drug Administration (FDA) and our python library to introduce blood vasculature under the skin layer. In each NBP, four different-sized (<10 mm in diameter) numerical lesion phantoms (NLPs) were inserted at locations randomly selected from those predicted by the VICTRE tools based on the duct and terminal duct lobular unit (TDLU) structures that are well-known sites for lesion formation. The considered tissue types and their unsigned 8-bit integer (uint8) labels are as belows. Background: 0 Fat: 1 Skin: 2 Glandular: 29 Nipple: 33 Ligament: 88 Terminal duct lobular unit (TDLU): 95 Duct: 125 Artery: 150 Vein: 225 Peripheral angiogenesis: 190 Viable tumor cell: 200 Necrotic core: 210 The functional, optical, and acoustic properties maps were produced employing our python libraries to assign corresponding properties to each breast tissue type. The optical fluence distributions were simulated using the MCX software [Fang2009], [Yu2018]. The induced initial pressure distributions were calculated via voxelwise multiplication of the optical absorption coefficient distributions and the simulated optical fluence distributions, assuming a Grüneisen Γ=1 as commonly done as constant for soft tissues. The acoustic measurements were simulated employing the k-wave software [Treeby2010]. Further details on 1) modifications and adaptations of the VICTRE NBPs for use in VITs of OAT; 2) virtual OAT imaging system and data acquisition are in the accompanying paper [Park2023]. The file naming convention of files in this dataset is {breast type}{seed number}{lesion presence}_{contained data}.mat, where: Breast type: A, B, C, or D; Seed number, i.e., a number randomly generated when producing each tissue label map; Lesion presence: absent (healthy, h) or present (l); Contained data: tissue label map (label), functional properties and chromophore concentrations maps (func), optical properties maps (opt), acoustic properties maps(acou), simulated optical fluence distributions (phi), simulated induced initial pressure distributions (p0), or simulated acoustic measurements (p). For example, the file name of a tissue label map of the type A breast with no lesion inserted (healthy breast) that was created using the seed number 123456 is A12345678h_label.mat. The actual data and metadata contained in each file are as below. {breast type}{seed number}{lesion presence}_label.mat Data label: 1360 x 1360 x 680 unit8 tissue label map that includes healthy tissues, viable tumor cell region, and necrotic core Metadata origin: 3 x 1 float32 values that specify x, y, and z coordinates of origin, (-85, -85, -85) mm voxel_size: float32 value that specifies voxel size, 0.125 mm tissue_type_label: 13 x 2 cell that specifies tissue type names and the corresponding uint8 values in label breast_type: Breast type, A: Breast is almost entirely fatty; B: Breast has scattered areas of fibroglandular density; C: Breast is heterogeneously dense; or D: Breast is extremely dense lesion_presence: Lesion presence, h: Lesion-absent numerical breast phantom (healthy) or l:...