A reconstruction model based on convex-nonconvex finite element total variation (CNC-FETV) regularization is proposed to avoid the biased estimation of L1 regularization in diffuse optical tomography. First, the finite element method was used to divide the computational domain into a finite number of triangles, after which a continuous piecewise polynomial function was used to approximate the absorption coefficient value in each triangle. Then, the derived difference matrix was assembled element by element to obtain a representation of the FETV regularization. Subsequently, the CNC-FETV regularization was obtained by the construction method based on convex-nonconvex sparse regularization. Results theoretically proved that the nonconvex regularization term could maintain the overall convexity of the objective function under certain conditions. Finally, the alternating direction multiplier method was used to solve the proposed model. Numerical experiments show that compared with the Tikhonov and FETV regularization models, the proposed CNC-FETV regularization model has superior performances in both numerical criteria and visual effects for diffusion optical tomography reconstructions.