Compensation for non-uniform illumination and optical fluence attenuation in three-dimensional optoacoustic tomography of the breast

Optoacoustic tomography (OAT) is a promising modality for breast imaging that provides high resolution, detection sensitivity and diagnostic specificity for vascularized breast tumors. In OAT systems employing an arc- shaped illuminator, irregular overlaps of light beams can yield a non-uniform illumination throughout the entire volume of the breast. The imbalance in optical fluence leads to intensity loss in the reconstructed OAT images. Additionally, because optical fluence decreases with depth from breast skin surface, i.e., optical attenuation, deep breast tissues are diminished in the reconstructed images. For qualitative enhancement in 3D OAT imaging, we propose an image processing method to estimate, and compensate for, both the non-uniform incident optical fluence and the optical attenuation. We approximate the non-uniform illumination via maximum intensity extraction for polar angles in a spherical coordinate system. The location of the breast surface is estimated by detecting blood vessels nearest to the breast skin layer that appear with relatively high intensities in the reconstructed image. The breast depth is computed as the minimum distance between each voxel and the detected breast surface. The depth-dependent optical attenuation in the breast is estimated using the Beer- Lambert law down to the maximum penetration depth determined from an analysis of noise and artifacts in the reconstructed image. At each polar angle, the reciprocals of the estimated attenuation is used to compensate for the loss in intensity. The results are that previously invisible structures near the chest wall are revealed, and visible penetration depth was increased by 67% over the conventional, non-compensated volumes.