1. A method of computing image gradients comprising the steps of: identifying a set of voxels; neglecting all voxels within the set of voxels having a gradient length below a predetermined threshold length; and comparing remaining voxels to neighboring voxels to identify local maxima along the gradient.