#Analyses, figures, tables and data are also produced in the script.

#AUTHOR: Benoit Parmentier 

#CREATED ON: 04/14/2015  

#Version: 6

#PROJECT: Environmental Layers project     

#COMMENTS: analyses run for reg4 1991 for test of mosaicing using 1500x4500km and other tiles

      #browser()

      lf_accuracy_testing_raster <- generate_ac_assessment_layers_by_tile(lf,layers_option,df,df_tile_processed,metric_name,

                                                    var_pred,list_models,use_autokrige,pred_mod_name,

                                                    y_var_name,interpolation_method,

      #undebug(mosaicFiles)

      #can also loop through methods!!!

                                mosaic_method="use_edge_weights",

                                num_cores=num_cores,

                                r_mask_raster_name=infile_mask,

                                scaling=scaling,

                                values_range=values_range)

      ##Took 12-13 minutes for 28 tiles and one date...!!! 

      if(tmp_files==F){ #if false...delete all files with "_tmp"

        ##now delete temporary files...

      }

      list_mosaic_obj[[i]] <- mosaic_obj

    }

      days_to_process <- day_to_mosaic[i]

      #browser()

      lf_accuracy_training_raster<- generate_ac_assessment_layers_by_tile(lf,layers_option,df,df_tile_processed,metric_name,

                                                    var_pred,list_models,use_autokrige,pred_mod_name,

                                                    y_var_name,interpolation_method,

      #undebug(mosaicFiles)

      #can also loop through methods!!!

      #browser()

                                mosaic_method="use_edge_weights",

                                num_cores=num_cores,

                                r_mask_raster_name=infile_mask,

                                data_type=data_type,

                                scaling=scaling,

                                values_range=values_range)

      list_mosaic_obj[[i]] <- mosaic_obj

    }

      days_to_process <- day_to_mosaic[i]

      #browser()

      lf_accuracy_residuals_testing_raster <- generate_ac_assessment_layers_by_tile(lf,layers_option,df,df_tile_processed,metric_name,

                                                    var_pred,list_models,use_autokrige,pred_mod_name,

                                                    y_var_name,interpolation_method,

      mosaic_method <- "use_edge_weights" #this is distance from edge

      out_suffix_tmp <- paste(interpolation_method,"kriged_residuals","data_day_v",day_to_mosaic[i],out_suffix,sep="_")

      #lf_tmp<-list.files(pattern="*kriged_residuals.*.tif",full.names=T)

      #lf_accuracy_residuals_raster[[i]]

      #debug(mosaicFiles)

      mosaic_obj <- mosaicFiles(lf_tmp,

                                data_type=data_type,

                                scaling=scaling,

                                values_range=values_range)

      list_mosaic_obj[[i]] <- mosaic_obj

    }      

      days_to_process <- day_to_mosaic[i]

      #browser()

      lf_accuracy_residuals_training_raster <- generate_ac_assessment_layers_by_tile(lf,layers_option,df,df_tile_processed,metric_name,

                                                    var_pred,list_models,use_autokrige,pred_mod_name,

                                                    y_var_name,interpolation_method,

                                data_type=data_type,

                                scaling=scaling,

                                values_range=values_range)

      list_mosaic_obj[[i]] <- mosaic_obj

      #Took 11 to 12 minues for one day and 28 tiles in region 4

    }

    ##End of mosaicing function for region predictions