In view of several potential applications in multivariable two-dimensional (2-D) systems theory, a practical 2-D matrix Pade approximation is introduced by using a generalized inverse of the matrices. The approximants are expressed in the form of the 2-D Thiele-type continued fractions and are computed by an efficient recursive algorithm. As it's an application, the state-space realization problem of the 2-D filters is discussed.