Spatial trends represent an obstacle to genetic evaluation in maize breeding. Spatial analyses can correct spatial trends, which allow for an increase in selective accuracy. The objective of this study was to compare the spatial (SPA) and non-spatial (NSPA) models in diallel multi-environment trial analyses in maize breeding. The trials consisted of 78 inter-populational maize hybrids, tested in four environments (E1, E2, E3, and E4), with three replications, under a randomized complete block design. The SPA models accounted for autocorrelation among rows and columns by the inclusion of first-order autoregressive matrices (AR1 ⊗ AR1). Then, the rows and columns factors were included in the fixed and random parts of the model. Based on the Bayesian information criteria, the SPA models were used to analyze trials E3 and E4, while the NSPA model was used for analyzing trials E1 and E2. In the joint analysis, the compound symmetry structure for the genotypic effects presented the best fit. The likelihood ratio test showed that some effects changed regarding significance when the SPA and NSPA models were used. In addition, the heritability, selective accuracy, and selection gain were higher when the SPA models were used. This indicates the power of the SPA model in dealing with spatial trends. The SPA model exhibits higher reliability values and is recommended to be incorporated in the standard procedure of genetic evaluation in maize breeding. The analyses bring the parents 2, 10 and 12, as potential parents in this microregion.