In this paper, we will first introduce a novel multiscale representation (MSR) for shapes via level set motions and partial differential equations (PDEs). Based on the MSR, we will then design a surface inpainting algorithm to recover 3D geometry of blood vessels. Because of the nature of irregular morphology in vessels and organs, both phantom and real inpainting scenarios were tested using our new algorithm. Successful vessel recoveries are demonstrated with numerical estimation of the degree of arteriosclerosis and vessel occlusion. Key words. Level set motion, mean curvature, multiscale representation, partial differential equations, surface inpainting. AMS subject classifications. 92C55, 68U10, 65D18, 53C44.