Multibaseline 2-D phase unwrapping (PU) is a critical step for the multibaseline synthetic aperture radar interferometry. Compared with the single-baseline PU, the multibaseline PU does not need to obey the phase continuity assumption, i.e., it is applicable to the terrain with the violent change. However, the performance of the multibaseline PU is directly related to noise level. In order to improve the noise robustness of the multibaseline PU, in this paper, we transplant the framework of the singlebaseline PU into the multibaseline PU and propose a two-stage programming approach, referred to as TSPA, which makes use of the gradient information of the interferogram similar to how the conventional single-baseline PU method does. Fortunately, although the proposed method belongs to the integer programming (usually, the integer programming is an NP-hard problem which is hard to solve), the constraint of the optimization model of the TSPA method is unimodular, so it can be efficiently solved. Furthermore, interestingly, some useful and important concepts of the single-baseline PU, for example, residue and branch cut, are also transplanted into the multibaseline PU in this paper, and we discuss the potential of extending most of the representative single-baseline PU methods into the multibaseline domain as well. Finally, the experiment results show the effectiveness and noise robustness of the TSPA multibaseline PU method.Index Terms-Multibaseline, robustness, single-baseline, synthetic aperture radar interferometry, unimodular, 2-D phase unwrapping (PU).