Haddar and Kress [9] extended the use of the conformal mapping approach [2, 8] to reconstruct the internal boundary curve Ti of a doubly connected domain from the Cauchy data on the external boundary of a harmonic function satisfying a homogeneous impedance boundary condition on Ti. However, the analysis of this scheme indicates non convergence of the proposed algorithm for small values of the impedance. In this paper, we modify the algorithm proposed in [9] in order to obtain a convergent and stable inversion process for small impedances. We illustrate the performance of the method through some numerical examples that also include the cases of variable impedances.