Using the normalized difference vegetation index, soil moisture, surface temperature and rainfall a non-linear multivariate optimization method has been developed that predicts crop yield at field and regional levels.