The aim of this paper is to propose a new regularized optimal control formulation for recovering an extended inclusion from boundary measurements. Our approach provides an optimal representation of the shape of the inclusion. It guarantees local Lipschitz stability for the reconstruction problem. Some numerical experiments are performed to demonstrate the validity and the limitations of the proposed reconstruction method.