In this paper, we introduce an iterative scheme for phase demodulation of interferograms with nonuniformly spaced phase shifts. Our proposal consists of two stages: first, the phase map is obtained through a least squares fitting; second, the phase steps are retrieved using a statistical robust estimator. In particular, we use Tukey’s biweighted M-estimator because it can cope with both noisy data and outliers in comparison with the ordinary least squares estimator. Furthermore, we provide the frequency description of the algorithm and the phase demodulation allowing us to analyze the procedure and estimation according to the frequency transfer function (FTF) formalism for phase-shifting algorithms. Results show that our method can accurately retrieve the phase map and phase shifts, and it converges by the 10th iteration.