In this paper, we present a shape recovery technique in 2D and 3D with specific applications in visualizing and measuring anatomical shapes from medical images. This algorithm models extremely corrugated structures like the brain, is topologically adaptable, is robust, and runs in O(NlogN) time where N is the total number of points in the domain. Our two-stage technique is based on the level set shape recovery scheme introduced in [11, 12, 4] and the fast marching method in [19] for computing solutions to static Hamilton-Jacobi equations.