In this work, the aim was to evaluate the existence of spatial association of the municipal average official soybean yield (t ha-1) with agrometeorological data and vegetation indices. The information was observed by ten-day periods, in crop years 2010/2011, 2011/2012 and 2012/2013 in the State of Paraná. Local univariate spatial correlation (LISA index), as well as global bivariate correlation (L statistics) were calculated. With this study, we identified neighboring municipalities with high yield in the West as well as municipalities that are located with low-low yield Northwestern, showing positive spatial autocorrelation (IMG=1), significative (p-value < 0.05). In addition, there were differences between seeding times in different regions, and climate irregularity during flowering periods and grain filling in crop year 2011/2012 throughout the state, which caused a large drop in production in all municipalities of the state of Paraná. The analysis of local spatial association showed that in the three crop years, the Northwest region presented a significant low yield potential of soybean (p-value < 0.05). In addition, it was observed that the period from the 3rd ten-day period of October to the 2nd ten-day period of January was essential for the soybean cycle in the different regions of the state, since this period encompasses the critical phases of crop. Differences were also observed between the crop years studied, regarding the agrometeorological variables, which affected soybean yield mainly in the Western region of Paraná – Brazil.