In forensic craniofacial reconstruction, facial features of an unknown individual are estimated from an unidentified skull, based on a mixture of experimentally obtained guidelines on the relationship between soft tissues and the underlying skeleton. In this paper, we investigate the possibility of using full 3D cross-sectional CT images for establishing a reference database of densely sampled distances between the external surfaces of the skull and head for automated craniofacial reconstruction. For each CT image in the reference database, the hard tissue (skull) and soft tissue (head) volumes are automatically segmented and transformed into signed distance transform (sDT) images, representing for each voxel in this image the Euclidean distance to the closest point on the skull and head surface, respectively, distances being positive (negative) for voxels inside (outside) the skull/head. Multiple craniofacial reconstructions are obtained by first warping each reference skull sDT maps to the target skull sDT using a B-spline based free form deformation algorithm and subsequently applying these warps to the reference head sDT maps. A single reconstruction of the target head surface is defined as the zero level set of the arithmetic average of all warped reference head sDT maps, but other reconstructions are possible, biasing the result to subject specific attributes (age, BMI, gender). Both qualitative and quantitative tests (measuring the similarity between the 3D reconstructed and corresponding original head surface) on a small (N = 20) database are presented to proof the validity of the concept.