A fitting system is developed to fit subdivision surface with sharp feature from noising and dense triangular meshes of arbitrary topology. The system includes an improved mesh denoising method based on bilateral filtering of images, sharp features extraction, feature-preserving mesh simplification and topological optimization. An estimating method for Loop subdivision surface is introduced to predict how many subdivision iterations are necessary to meet a user-defined tolerance. The method of adaptive subdivision is proposed during the fitting process to handle local detailed surface features. Both experimental results and practical applications in engineering demonstrate that the system can effectively achieve a good quality of the fitting subdivision surface with nice details while using few facets in the approximation.