This paper presents an iterative method for the unbiased identification of linear Multiple-Input Multiple-Output (MIMO) discrete two-dimensional (2D) systems. The system discussed here has Auto-Regressive Moving-Average model with exogenous inputs (ARMAX model). The proposed algorithm functions on the basis of the traditional Iterative Generalized Least Squares (IGLS) method. In summary, this paper proposes a two-dimensional Multiple-Input Multiple-Output Iterative Generalized Least Squares (2DMIGLS) algorithm to estimate the unknown parameters of the ARMAX model. Finally, simulation results show the efficiency and accuracy of the presented algorithm in estimating the unknown parameters of the model in the presence of colored noise.