tb_date <- subset(tb,date==date_processed & pred_mod==pred_mod_name)

  tb_date$tile_coord <- as.character(tb_date$tile_coord)

  df_centroids <- merge(df_centroids,tb_date,by="tile_coord")

  #use mclapply  

  #list_raster_name <- vector("list",length=length(lf))

  list_param_raster_df_centroids <- list(df_centroids,metric_name,NA_flag_val,file_format,out_dir_str,out_suffix_str)

  #undebug(create_raster_df_centroids_fun)

  #test_lf <- lapply(1,FUN=create_raster_df_centroids_fun,list_param=list_param_raster_df_centroids)                           

  raster_created_obj <- list(list_raster_name,df_centroids)

  names(raster_created_obj) <- c("list_raster_name","df_centroids")

    #debug(create_accuracy_metric_raster)

    #list_raster_created_obj <- lapply(1:1,FUN=create_accuracy_metric_raster,

    #                                  list_param=list_param_accuracy_metric_raster)

    #Extract list of files for rmse and date 1 (19920101), there should be 28 raster images

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