#Different options to explore mosaicing are tested. This script only contains functions.

#AUTHOR: Benoit Parmentier 

#CREATED ON: 04/14/2015  

#Version: 2

#PROJECT: Environmental Layers project     

#COMMENTS: first commit of function script to test mosaicing using 1500x4500km and other tiles

  #NAflag,file_format,out_suffix etc...

  file_format <- list_param$file_format

  out_dir_str <- list_param$out_dir

  ####### START SCRIPT #####

  #use mclapply  

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

  names(list_param_raster_df_centroids) <- c("df_centroids","metric_name","NA_flag_val","file_format","out_dir","out_suffix")

  #undebug(create_raster_df_centroids_fun)

    max_val <- cellStats(r_dist,max)

    r <- abs(r_dist - max_val)/ (max_val - min_val)

  }

  extension_str <- extension(lf[i])

  raster_name_tmp <- gsub(extension_str,"",basename(lf[i]))

  raster_name <- file.path(out_dir_str,paste(raster_name_tmp,"_",method,"_weights_",out_suffix_str,file_format,sep=""))

    #names(list_param_create_weights) <- c("lf","df_points","r_feature","method","out_dir_str") 

    #num_cores <- 11

    #undebug(create_weights_fun)

    #This is the function creating the weights by tile. Distance from the centroids needs to be change from distance to

    #the edges...can use rows and columsn to set edges to 1 and 0 for the others.

      #input_file <- filename_list_mosaics_weights_m

      module_path <- mosaic_python #this should be a parameter for the function...

      #python /nobackupp6/aguzman4/climateLayers/sharedCode/gdal_merge_sum.py 

      #--config GDAL_CACHEMAX=1500 --overwrite=TRUE -o  outputname.tif --optfile input.txt

      #r_weights_sum_raster_name <- mosaic_python_merge(module_path=mosaic_python,