Geodesic methods have been widely applied to image analysis. They are particularly efficient to extract a tubular structure, such as a blood vessel, given its two endpoints in a 2D or 3D medical image. We address here a more difficult problem: the extraction of a full vessel tree structure given a single initial root point, by growing a collection of keypoints or new initial source points, connected by minimal geodesic paths. In this article, those keypoints are iteratively added, using a new detection criteria, which utilize the weighted geodesic distances with respect to a radius-lifted Riemannian metric, the standard Euclidean curve length and a path score. Two main weaknesses of classical keypoints searching approach are that the weighted geodesic distance and the Euclidean path length do not take into account the orientation of the tubular structure or object boundaries, due to the use of an isotropic geodesic Riemannian metric, and suffer from a leakage problem. In contrast, we use an anisotropic geodesic Riemannian metric, and develop new criteria for selecting keypoints based on the path score and automatically stopping the tree growth. Experimental results demonstrate that our method can obtain the expected results, which can extract vessel structures at a finer scale, with increased accuracy.