Approximating surfaces by the moving least absolute deviations method