如何根据地区多年气温，降水，ndvi栅格数据，拟合历年的ndvi预测栅格数据

是先以ndvi为因变量，气温降水为自变量，建立多元线性回归，然后根据当年气候数据，求当年残差，最后求ndvi预测值么
网上只找到有关残差的matlab代码，数据量太大了，总是不成功