Magnetic resonance angiography (MRA) is an imaging technique to show the blood vessels but suppress signals from all the other tissues. There are two approaches to acquire an image of MRA. One is done in two-dimension (2D) by projecting the three-dimensional (3D) vessels onto 2D plane. The other is to directly obtain the complete 3D information. The advantage of 3D MRA is that one can view the data from arbitrary direction. However, the scan time is usually very long for 3D MRA. When the scan time is limited, we must use 2D MRA and the depth information is sacrificed. The depth information can be partially recovered using stereoscopic angiography, i.e., acquiring two 2D images from different viewing angles. The modality used more frequently for stereoscopic angiography is digital subtraction X-ray angiography (DSA). For X-ray DSA, the pixel intensity is directly related to the integration of the attenuation value in the path of the X-ray. We can derive the vessel shape by solving the integrals from the two views. On the other hands, there are many imaging parameters in MRA (such as T1, T2, and proton density). Therefore, it is difficult to obtain the relation between the vessel shape and the pixel intensity. For this reason, we attempt to reconstruct the shape of the vessels simply from the edge information of the two views. We assume that the shape of the vessel on every cross-section is an ellipse. Then, we develop an algorithm to estimate the parameters of the ellipse from the boundaries of the projective images. A 3D MRA data set was used to test the capability of our algorithm. From these data, we make two projective images. The two projections are 30 degree(s) apart. We employed our algorithm to estimate all the ellipses and reconstruct the 3D model of the vessels. Comparing the boundaries of the original projective images with the boundaries of the reconstructed 3D model, the average error is 0.471 pixels.