The determination of an interface in an elliptic differential equation from Dirichlet and Neumann boundary data is investigated. Uniqueness of the function characterizing the interface is proved. For numerical purposes the problem is formulated as an optimization problem involving constraints that are based on potential theory. An augmented Lagrangian formulation is used for the solution of the optimization problem. Its convergence is investigated and numerical experiments are described.