A swarm intelligence based method is presented in this paper to solve the problem of parameterization of given set of planar and spatial data points, which may be uniformly or non-uniformly spaced, for least squares fitting with Bezier curves. The basic framework of the proposed method relies on fireworks algorithm. Performance of the proposed method is tested through several example curves of varying complexities.