# 5)runGAMFusion <- function(i,list_param) : daily step for fusion method, perform daily prediction

#

#AUTHOR: Benoit Parmentier                                                                       

#PROJECT: NCEAS INPLANT: Environment and Organisms --TASK#363--   

##Comments and TODO:

#This script is meant to be for general processing tile by tile or region by region.

# Note that the functions are called from GAM_fusion_analysis_raster_prediction_mutlisampling.R.

# This will be expanded to other methods.

##################################################################################################

  clim_obj<-list(rast_clim_list,data_month,mod_list,list_formulas)

  names(clim_obj)<-c("clim","data_month","mod","formulas")

  return(clim_obj) 

}

  names(rast_clim_list)<-names(rast_bias_list)

  for (k in 1:nlayers(mod_rast)){

    clim_fus_rast<-LST-subset(mod_rast,k)

                     "_",run_samp,sep="")

    raster_name<-paste("fusion_",data_name,out_prefix,".tif", sep="")

    rast_clim_list[[k]]<-raster_name

  bias_rast<-interpolate(LST,fitbias) #interpolation using function from raster package

  #Saving kriged surface in raster images

                   "_",run_samp,sep="")

  raster_name_bias<-paste("fusion_",data_name,out_prefix,".tif", sep="")

  writeRaster(bias_rast, filename=raster_name_bias,overwrite=TRUE)  #Writing the data in a raster file format...(IDRISI)

  #now climatology layer

  clim_rast<-LST-bias_rast

                   "_",run_samp,sep="")

  raster_name_clim<-paste("fusion_",data_name,out_prefix,".tif", sep="")

  writeRaster(clim_rast, filename=raster_name_clim,overwrite=TRUE)  #Writing the data in a raster file format...(IDRISI)

  clim_obj<-list(rast_bias_list,rast_clim_list,data_month,mod_list,list_formulas)

  names(clim_obj)<-c("bias","clim","data_month","mod","formulas")

  return(clim_obj)

}

  daily_delta_rast<-interpolate(rast_clim_month,fitdelta) #Interpolation of the bias surface...

  #Saving kriged surface in raster images

                   "_",sampling_dat$run_samp[i],sep="")

  raster_name_delta<-paste("fusion_",var,"_",data_name,out_prefix,".tif", sep="")

  writeRaster(daily_delta_rast, filename=raster_name_delta,overwrite=TRUE)  #Writing the data in a raster file format...(IDRISI)

  obj_names<-c(y_var_name,"clim","delta","data_s","data_v",

               "sampling_dat",model_name)

  save(delta_obj,file= paste("delta_obj_",var,"_",sampling_dat$date[i],"_",sampling_dat$prop[i],

                                "_",sampling_dat$run_samp[i],out_prefix,".RData",sep=""))

  return(delta_obj)