In this paper, the link between least-squares (LS) estimator and Laplacian smoothing method from the point of view of LS method is firstly disclosed. Further more, an M-estimator for mesh filtering is presented. At last, the M-estimator is extended to the re-weighted M-estimator to remove noise while preserving the surface feature efficiently. The re-weighted M-estimator is bilateral filtering in natural.