This paper presents a fast algorithm （BS2 Algorithm） for fitting C 1 surfaces to scat- tered data points. By using energy minimization, the bivariate spline space S 2 1（△ m,n （2） ） is introduced to construct a Cl-continuous piecewise quadratic surface through a set of irregularly 3D points. Moreover, a multilevel method is also presented. Some experimental results show that the accuracy is satisfactory. Furthermore, the BS2 Algorithm is more suitable for fitting surfaces if the given data points have some measurement errors.
