#MODIFIED ON: 12/17/2015            

function_mosaicing <-"global_run_scalingup_mosaicing_function_12172015.R"

#Take into account multiple region in some cases!!!  

reg_lf_mosaic <- vector("list",length=length(region_names))

for(k in 1:length(region_names)){

  in_dir_tiles_tmp <- file.path(in_dir_tiles, region_names[k])

  reg_lf_mosaic[[k]] <- lapply(1:length(day_to_mosaic),FUN=function(i){list.files(path=file.path(in_dir_tiles_tmp),    

           pattern=paste(".*.",day_to_mosaic[i],".*.tif$",sep=""),full.names=T,recursive=F)})

}

#use reg4 to test the code for now, redo later for any regions!!!

lf_mosaic <- reg_lf_mosaic[[2]] #list of files to mosaic for date 1, there a 28 files for reg4, South America

num_cores_tmp <- 7

#interpolation_method #c("gam_CAI") #PARAM3, need to select reg!!

df_tile_processed$path_NEX <- as.character(df_tile_processed$path_NEX)

df_tile_processed$reg <- basename(dirname(df_tile_processed$path_NEX))

##By regions, selected earlier

for(k in 1:length(region_names)){

  df_tile_processed_reg <- subset(df_tile_processed,reg==region_names[k])#use reg4

  #i<-1 #loop by days/date to process!!

  #test on the first day 

  list_param_create_accuracy_residuals_raster <- list(lf_day_tiles,data_df,df_tile_processed_reg,

  names(list_param_create_accuracy_residuals_raster) <- c("lf_day_tiles","data_df","df_tile_processed_reg",

  #list_create_accuracy_residuals_raster_obj <- lapply(1:length(day_to_mosaic),FUN=create_accuracy_residuals_raster_obj,

  #                                list_param=list_param_create_accuracy_residuals_raster_obj)

  #undebug(create_accuracy_residuals_raster)

  list_create_accuracy_residuals_raster_obj <- lapply(1:1,FUN=create_accuracy_residuals_raster,

  #create_accuracy_residuals_raster_obj <- create_accuracy_metric_raster(1, list_param_create_accuracy_residuals_raster_obj)

}

lf_accuracy_residuals_raster <- lapply(1:length(list_raster_created_obj),FUN=function(i){unlist(list_raster_created_obj[[i]]$list_raster_name)}) 

#Plot as quick check

r1 <- raster(lf_mosaic[[1]][1]) 

list_create_accuracy_residuals_raster_obj

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

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

  mosaic_edge_obj_residuals <- mosaicFiles(lf_accuracy_residuals_raster[[i]],

                                        mosaic_method="use_edge_weights",

                                        num_cores=num_cores,

                                        r_mask_raster_name=infile_mask,

                                        python_bin=python_bin,

                                        mosaic_python=mosaic_python,

                                        algorithm=algorithm,

                                        df_points=NULL,

                                        NA_flag=NA_flag_val,

                                        file_format=file_format,

                                        out_suffix=out_suffix_tmp,

                                        out_dir=out_dir)

  list_mosaic_obj[[i]] <- list(prediction=mosaic_edge_obj_prediction,accuracy=mosaic_edge_obj_accuracy)