We present a phase unwrapping method using a maximum likelihood (ML) estimation technique together with frequency diversity information to reconstruct highly discontinuous ground elevation profiles. Frequency diversity can be obtained by considering the interferograms obtained by different couples of subband images.