An iterative formulation for the least squares approximation of 2-D signals