A fast and automated method is proposed to get a spline-based model from voxel data. A hierarchical basis of C1C1 cubic polynomials provides the volumetric parametrization. The shape of the boundary is approximated using a simple template in 2D and 3D. Complex domains with voids can be handled in a single-patch framework.