Abstract. Fuzzy boundaries of anatomical structures in medical images make segmentation a challenging task. We present a new segmentation method that addresses the fuzzy boundaries problem. Our method maps the lengths of 3D rays cast from a seed point to the unit sphere, estimates the fuzzy boundaries location by thresholding the gradient magnitude of the rays lengths, and derives the true boundaries by Laplacian interpolation on the sphere. Its advantages are that it does not require a global shape prior or curvature based constraints, that it has an automatic stopping criteria, and that it is robust to anatomical variability, noise, and parameters values settings. Our experimental evaluation on 23 segmentations of kidneys and on 16 segmentations of abdominal aortic aneurysms (AAA) from CT scans yielded an average volume overlap error of 12.6% with respect to the ground-truth. These results are comparable to those of other segmentation methods without their underlying assumptions.