This is a demonstration of a 2D phase unwrapping routine I wrote in MATLAB to extract topographical measurements from interferometry data generated by a custom-built Fizeau interferometer. The algorithm implements the Goldstein branch cut method. It finds local errors in the wrapped phase due to aliasing, singularities, or noise that make the unwrapping procedure path-dependent, and routes the unwrapping path around those points.