A method for three-dimensional imaging of a sample (302) comprises: receiving (102) interference patterns (208) acquired using light-detecting elements (212), wherein each interference pattern (208) is formed by scattered light from the sample (302) and non-scattered light from a light source (206; 306), wherein the interference patterns (208) are acquired using different angles between the sample (302) and the light source (206; 306); performing digital holographic reconstruction applying an iterative algorithm to change a three-dimensional scattering potential of the sample (302) to improve a difference between the received interference patterns (208) and predicted interference patterns based on the three-dimensional scattering potential; wherein the iterative algorithm reduces a sum of a data fidelity term and a non-differentiable regularization term and wherein the iterative algorithm includes a forward-backward splitting method alternating between forward gradient descent (108) on the data fidelity term and backward gradient descent (110) on the regularization term.