This paper presents a variational model for the denoising of wrapped phase images. By enforcing the required Pythagorean trigonometric identity between the real and imaginary components of the signal, this model improves the signal-to-noise ratio of the restored signal. To preserve phase map discontinuities, the model is based on total variation. The existence and uniqueness of the model’s solution are demonstrated using standard techniques. In addition, the convergence of a rapid fixed-point method to determine the numerical solution is demonstrated. Experiments on both synthetic and actual patterns validate the model’s performance.