Original image is directly processed by the existing image segmentation algorithms, which is easily affected by noise. A bi-dimensional ensemble empirical mode decomposition (BEEMD) method is introduced to improve the accuracy of MR image segmentation by distance regularized level set (DRLSE) method. The BEEMD method is the extension of one-dimensional noise assisted data analysis from ensemble empirical mode decomposition (EEMD). The key points of BEEMD are as follows. four-neighborhood optimization is used to find extermum; three-spline interpolation is used to obtain the envelope; amplitude standard of added white noise is restricted; a certain time of integration is used to avoid modality aliasing problem. The main steps of the proposed method are as follows. Firstly, the MR image is decomposed into a number of two-dimensional intrinsic mode functions (BIMF) by BEEMD method; different weighting coefficients are endued to BIMF for image reconstruction to enhance the segmentation target. Secondly, part of BIMF components are added into edge indicator function of DRLSE to recover the blurring boundary caused by Gauss smooth operation. Then DRLSE is used to segment the reconstructed MR image. High accuracy and robustness of proposed algorithm are obtained in both simulations and clinical MR images. However, compared with DRLSE, the proposed method is complex and time consuming because using BEEMD for preprocessing the segmentation image.