Based on remote sensing, the study analyzed the dynamics of a wheat crop using regression analysis to estimate production in relation to the NDMI, NDVI, MSAVI and NBR indices. The wheat crop considered in the study, was located in the area of Sacalaz, Timis county, Romania. Seven series of satellite images were taken, between February and June 2021. Based on the spectral information, the NDMI, NDVI, MSAVI and NBR indices were calculated and series of 595 values were obtained for each index. Very strong correlations were recorded between NBR and NDMI (r=0.999***), between MSAVI and NDVI (r=0.994***), between NBR and MSAVI (r=0.919**), between NDVI and NDMI ( r=0.911) and between NBR and NDVI (r=0.909**). Polynomial equations described the variation of NBR in relation to NDMI (R2=0.987), in relation to MSAVI (R2=0.794) and in relation to NDVI (R2=0.795). The variation of indices in relation to time, during the study period, was described by polynomial equations of the 2nd and 3rd degree, under statistical accuracy conditions (R2=0.854 for NDMI; R2=0.919 for NDVI; R2=0.956 for MSAVI, and R2=0.873 for NBR). Through the regression analysis, the wheat production was predicted, based on the calculated indices, and in different combinations of indices. The highest level of accuracy in the prediction of wheat production was recorded in the case of using the NDMI and NBR indices (p<0.001, RMSPE=0.24847), followed by the analysis variant in which the NDVI and MSAVI indices were considered (p=0.0002, RMSEP=1.59084) and the analysis variant in which NDVI and NBR indices were used (p=0.00081, RMSEP=5.20218). Graphical models, in 3D format and in the form of isoquants, have described the variation of wheat production in relation to the NDMI, NDVI, MSAVI and NBR indices.