Intensity inhomogeneities in images cause problems in gray-value based image segmentation since the varying intensity often dominates over gray-value differences of the image structures. In this paper we propose a novel biconvex variational model that includes the intensity inhomogeneities to tackle this task. We combine a total variation approach for multi class segmentation with a multiplicative model to handle the inhomogeneities. In our model we assume that the image intensity is the product of a smoothly varying part and a component which resembles important image structures such as edges. Therefore, we penalize in addition to the total variation of the label assignment matrix a quadratic difference term to cope with the smoothly varying factor. A critical point of the resulting biconvex functional is computed by a modified proximal alternating linearized minimization method (PALM). We show that the assumptions for the convergence of the algorithm are fulfilled. Various numerical examples demonstrate the very good performance of our method. Particular attention is paid to the segmentation of 3D FIB tomographical images serving as a motivation for our work.