In this paper we generalize the class of STARMA models to include contemporaneous spatial effects. Under the assumption of stationarity, invertibility and structural identifiability, we derive maximum likelihood estimators of the model coefficients by means of the Kalman filter recursions and the scoring algorithm, presenting some simulation results