Project

General

Profile

« Previous | Next » 

Revision 1cef03c8

Added by Benoit Parmentier over 8 years ago

mosaicing script master, testing all options for mosaicing

View differences:

climate/research/oregon/interpolation/master_script_stage_7.R
14 14 

  		




  15
  15
  
    
#AUTHOR: Benoit Parmentier                                                                        


  		




  16
  16
  
    
#CREATED ON: 01/01/2016  


  		




  17
  
  
    
#MODIFIED ON: 06/08/2016  


  		




  
  17
  
    
#MODIFIED ON: 06/09/2016  


  		




  18
  18
  
    
#PROJECT: NCEAS INPLANT: Environment and Organisms                                                                           


  		




  19
  19
  
    

  		




  20
  20
  
    
#First source these files:


  		




  ......




  65
  65
  
    

  		




  66
  66
  
    
#script_path <- "/home/parmentier/Data/IPLANT_project/env_layers_scripts"


  		




  67
  67
  
    
script_path <- "/nobackupp8/bparmen1/env_layers_scripts" #path to script


  		




  68
  
  
    
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_05312016.R" #PARAM12


  		




  69
  
  
    
function_mosaicing <-"global_run_scalingup_mosaicing_06082016.R"


  		




  
  68
  
    
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_06112016.R" #PARAM12


  		




  
  69
  
    
function_mosaicing <-"global_run_scalingup_mosaicing_06112016d.R"


  		




  70
  70
  
    
source(file.path(script_path,function_mosaicing)) #source all functions used in this script 


  		




  71
  71
  
    
source(file.path(script_path,function_mosaicing_functions)) #source all functions used in this script 


  		




  72
  72
  
    

  		




  ......




  109
  109
  
    
mosaicing_method <- "use_edge_weights" #PARAM10, arg 10


  		




  110
  110
  
    
metric_name <- "rmse" # "mae", "r" for MAE, R etc.; can also be ns or nv? #PARAM 11, arg 11


  		




  111
  111
  
    
#metric_name <- "n"


  		




  
  112
  
    
#metric_name <- "mae"


  		




  
  113
  
    

  		




  112
  114
  
    
#day_start <- "19990701" #PARAM 12 arg 12


  		




  113
  115
  
    
#day_end <- "19990703" #PARAM 13 arg 13


  		




  114
  
  
    
day_start <- "19990108" #PARAM 12 arg 12


  		




  
  116
  
    
day_start <- "19990107" #PARAM 12 arg 12


  		




  115
  117
  
    
day_end <- "19990108" #PARAM 13 arg 13


  		




  116
  118
  
    
infile_mask <- "/nobackupp8/bparmen1/NEX_data/regions_input_files/r_mask_LST_reg4.tif" #PARAM 14, arg 14


  		




  117
  119
  
    
df_assessment_files_name <- "/nobackupp6/aguzman4/climateLayers/out/reg4/assessment/output_reg4_1999/df_assessment_files_reg4_1999_reg4_1999.txt"  # data.frame with all files used in assessmnet, PARAM 15


  		




  118
  120
  
    
algorithm <- "python" #PARAM 16 #if R use mosaic function for R, if python use modified gdalmerge script from Alberto Guzmann


  		




  119
  121
  
    
#layers_option <- c("var_pred") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred


  		




  120
  
  
    
#layers_option <- c("ac_testing") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred


  		




  121
  
  
    
#layers_option <- c("res_testing") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred


  		




  122
  
  
    
#layers_option <- c("res_training") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred


  		




  123
  122
  
    
layers_option <- c("ac_training") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred


  		




  
  123
  
    
# #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred


  		




  
  124
  
    
#layers_option <- c("res_training") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred


  		




  
  125
  
    
#layers_option <- c("ac_testing") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred


  		




  124
  126
  
    

  		




  125
  127
  
    
tmp_files <- TRUE #arg 18, param 18, keep temp files if TRUE


  		




  126
  128
  
    
data_type <- "Int16" #, param 19, use int32 for output layers mosaiced


  		












Also available in:  Unified diff