Total Variation-Based Reduction of Streak Artifacts, Ring Artifacts and Noise in 3D Reconstruction from Optical Projection Tomography

Optical projection tomography (OPT) is a computed tomography technique at optical frequencies for samples of 0.5–15 mm in size, which fills an important “imaging gap” between confocal microscopy (for smaller samples) and large-sample methods such as fluorescence molecular tomography or micro magnetic resonance imaging. OPT operates in either fluorescence or transmission mode. Two-dimensional (2D) projections are taken over 360° with a fixed rotational increment around the vertical axis. Standard 3D reconstruction from 2D OPT uses the filtered backprojection (FBP) algorithm based on the Radon transform. FBP approximates the inverse Radon transform using a ramp filter that spreads reconstructed pixels to neighbor pixels thus producing streak and other types of artifacts, as well as noise. Artifacts increase the variation of grayscale values in the reconstructed images. We present an algorithm that improves the quality of reconstruction even for a low number of projections by simultaneously minimizing the sum of absolute brightness changes in the reconstructed volume (the total variation) and the error between measured and reconstructed data. We demonstrate the efficiency of the method on real biological data acquired on a dedicated OPT device.