Image Reconstruction for Diffuse Optical Tomography Based on Radiative Transfer Equation

Algorithm 1

Reconstruction algorithm based on the regularizing Levenberg-Marquardt method.

Input: Set the initial guess of ; The regularization parameter , .

Output: Approximate minimizer .

  (1)  for    do

  (2)  for    do

      (i) Compute the Fréchet derivative , and  

  (3)  end for

      (i) Compute , and

      (ii) Update by solving the linearization problem

      (

      (iii) Check the stopping criterion.

  (4) end for