A statistical algorithm for the reconstruction from time sequence echocardiographic images is proposed in this paper.The ability to jointly restore the images and reconstruct the 3D images without blurring the boundary is the main innovation of this algorithm. First, a Bayesian model based on MAP-MRF is used to reconstruct 3D volume, and extended to deal with the images acquired by rotation scanning method. Then, the spatiotemporal nature of ultrasound images is taken into account for the parameter of energy function, which makes this statistical model anisotropic. Hence not only can this method reconstruct 3D ultrasound images, but also remove the speckle noise anisotropically. Finally, we illustrate the experiments of our method on the synthetic and medical images and compare it with the isotropic reconstruction method.