Accurate segmentation of cerebral venous vasculature from gradient echo data is of central importance in several areas of neuroimaging such as for the susceptibility-based assessment of brain oxygenation or planning of electrode placement in deep brain stimulation. In this study, a vein segmentation algorithm for single- and multi-echo gradient echo data is proposed. First, susceptibility maps, true susceptibility-weighted images, and, in the multi-echo case, R2* maps were generated from the gradient echo data. These maps were filtered with an inverted Hamming filter to suppress background contrast as well as artifacts from field inhomogeneities at the brain boundaries. A shearlet-based scale-wise representation was generated to calculate a vesselness function and to generate segmentations based on local thresholding. The accuracy of the proposed algorithm was evaluated for different echo times and image resolutions using a manually generated reference segmentation and two vein segmentation algorithms (Frangi vesselness-based, recursive vesselness filter) as a reference with the Dice and Cohen's coefficients as well as the modified Hausdorff distance. The Frangi-based and recursive vesselness filter methods were significantly outperformed with regard to all error metrics. Applying the algorithm, susceptibility differences likely related to differences in blood oxygenation between superficial and deep venous territories could be demonstrated.