Three-dimensional freehand ultrasound imaging produces a set of irregularly spaced B-scans, which are typically reconstructed on a regular grid for visualization and data analysis. Most standard reconstruction algorithms are designed to minimize computational requirements and do not exploit the underlying shape of the data. We investigate whether an approximation with splines holds any promise as a better reconstruction method. A radial basis function approximation method is implemented and compared with three standard methods. While the radial basis approach is computationally expensive, it produces accurate reconstructions without the kind of visible artefacts common with the standard methods. The other potential advantages of radial basis functions, such as the direct computation of derivatives, make further investigation worthwhile.