A general method of surface profiling with phase-shifting interferometry techniques uses iterative linear regression to fit the sequence of interferograms to a physical model of the cavity. The physical model incorporates all important cavity influences, including environmentally induced rigid-body motion, phase shifter miscalibrations, multiple interference, geometry-induced spatial phase-shift variations, and their cross-couplings. By incorporating an initial estimate of the surface profile and iteratively solving for space- and time-dependent variables separately, convergence is robust and rapid. The technique has no restriction on surface shape or departure.