In order to solve the problem that the gridless DOA estimation algorithms based on generalized finite rate of innovation (FRI) signal reconstruction model are not suitable for two-dimensional DOA estimation using planar array, a separable gridless DOA estimation algorithm exploiting bi-orthogonal sparse linear array (BSLA) structure is proposed in this paper, which is called 2D-SGFRI. The 2D-SGFRI algorithm firstly recovers the covariance data of the virtual array formed by BSLA through the matrix completion method, so as to obtain the complete covariance data vectors about two independent parameters respectively. Next, since the covariance data vector satisfies the constraints of annihilation filter equations, the generalized FRI signal reconstruction model can be utilized to retrieve DOA from the covariance data vector. Compared with the existing DOA estimation algorithms based on generalized FRI signal reconstruction model, the 2D-SGFRI algorithm can be can be effectively applied to twodimensional DOA estimation, and can obtain stable estimation results. At the same time, due to the reduction of the dimension of positive semidefinite matrix, the 2D-SGFRI algorithm can significantly reduce the computational complexity compared with the two-dimensional DOA estimation algorithms based on atomic norm minimization (ANM). A series of simulation experiments are shown to verify the effectiveness and superiority of 2D-SGFRI algorithm. INDEX TERMS Two-dimensional DOA estimation, bi-orthogonal sparse linear array (BSLA), finite rate of innovation (FRI), matrix completion, annihilation filter