The aim of this paper is to achieve the 3D reconstruction of blood vessels from a limited number of 2D transversal cuts obtained from scanners. This is motivated by the fact that data can be missing. The difficulty of this work is to connect the blood vessels between some widely spaced cuts in order to produce the graph corresponding to the network of vessels. We identify the vessels on each transversal cut as a mass to be transported along a graph which allows to determine the bifurcation points of vessels. Specifically, we are interested in mass transfer problem Monge-Kantorovich [10] to model an optimized graph associated to the network of vessels. At this stage, we are able to reconstruct a 3D level set function by using the 2D level set functions given by the transversal cuts and the graph information. When the whole scanners data are available, a global reconstruction is proposed in a simple manner, without using the mass transfer problem.