Robust Fitting of Subdivision Surfaces for Smooth Shape Analysis