A B-spline level set (BLS) based method is proposed for shape reconstruction in electrical impedance tomography (EIT). We assume that the conductivity distribution to be reconstructed is piecewise constant, transforming the image reconstruction problem into a shape reconstruction problem. The shape/interface of inclusions is implicitly represented by a level set function (LSF), which is modeled as a continuous parametric function expressed using B-spline functions. Starting from modeling the conductivity distribution with the B-spline based LSF, we show that the shape modeling allows us to compute the solution by restricting the minimization problem to the space spanned by the B-splines. As a consequence, the solution to the minimization problem is obtained in terms of the B-spline coefficients. We illustrate the behavior of this method using simulated as well as water tank data. In addition, robustness studies considering varying initial guesses, differing numbers of control points, and modeling errors caused by inhomogeneity are performed. Both simulation and experimental results show that the BLS-based approach offers clear improvements in preserving the sharp features of the inclusions in comparison to the recently published parametric level set method.