Project

General

Profile

« Previous | Next » 

Revision 0cdba0f3

Added by Benoit Parmentier over 8 years ago

master script mosaicing debugging of accuracy layers

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: 05/30/2016  


  		




  
  17
  
    
#MODIFIED ON: 05/31/2016  


  		




  18
  18
  
    
#PROJECT: NCEAS INPLANT: Environment and Organisms                                                                           


  		




  19
  19
  
    

  		




  20
  20
  
    
#First source these files:


  		




  ......




  22
  22
  
    
#source /nobackupp6/aguzman4/climateLayers/sharedModules2/etc/environ.sh  


  		




  23
  23
  
    
#MODULEPATH=$MODULEPATH:/nex/modules/files


  		




  24
  24
  
    
#module load pythonkits/gdal_1.10.0_python_2.7.3_nex


  		




  25
  
  
    

  		




  
  25
  
    
#Set permissions:


  		




  
  26
  
    
#setfacl -Rm u:aguzman4:rwx /nobackupp6/aguzman4/climateLayers/LST_tempSpline/


  		




  
  27
  
    
#


  		




  26
  28
  
    
## TODO:


  		




  27
  29
  
    
# 


  		




  28
  30
  
    
## Comments: dealing with flags and datatype in output mosaics


  		




  ......




  63
  65
  
    

  		




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


  		




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


  		




  66
  
  
    
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_05302016.R" #PARAM12


  		




  67
  
  
    
function_mosaicing <-"global_run_scalingup_mosaicing_05302016.R"


  		




  
  68
  
    
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_05312016.R" #PARAM12


  		




  
  69
  
    
function_mosaicing <-"global_run_scalingup_mosaicing_05312016.R"


  		




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


  		




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


  		




  70
  72
  
    

  		




  ......




  108
  110
  
    
metric_name <- "rmse" #RMSE, MAE etc. #PARAM 11, arg 11


  		




  109
  111
  
    
#day_start <- "19990701" #PARAM 12 arg 12


  		




  110
  112
  
    
#day_end <- "19990703" #PARAM 13 arg 13


  		




  111
  
  
    
day_start <- "19990109" #PARAM 12 arg 12


  		




  112
  
  
    
day_end <- "19990111" #PARAM 13 arg 13


  		




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


  		




  
  114
  
    
day_end <- "19990108" #PARAM 13 arg 13


  		




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


  		




  114
  116
  
    
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


  		




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


  		




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


  		




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


  		




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


  		




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


  		




  118
  120
  
    

  		




  119
  
  
    
tmp_files <- FALSE #arg 18, param 18, keep temp files if TRUE


  		




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


  		




  120
  122
  
    
data_type <- "Int16" #, param 19, use int32 for output layers mosaiced


  		




  121
  123
  
    
scaling <- 100 #, param 20, if null use 1


  		




  122
  124
  
    
values_range <- "-100,100"


  		




  ......




  284
  286
  
    
#debug(debug_fun_test)


  		




  285
  287
  
    
#debug_fun_test(list_param_raster_prediction)


  		




  286
  288
  
    
i <- 1 #this select the first year of list_year_predicted


  		




  287
  
  
    
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_05302016.R" #PARAM12


  		




  288
  
  
    
function_mosaicing <-"global_run_scalingup_mosaicing_05302016.R"


  		




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


  		




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


  		




  
  289
  
    
#function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_05312016.R" #PARAM12


  		




  
  290
  
    
#function_mosaicing <-"global_run_scalingup_mosaicing_05312016.R"


  		




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


  		




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


  		




  291
  293
  
    

  		




  292
  294
  
    
if (stages_to_run[7]==7){


  		




  293
  295
  
    
  assessment_prediction_obj <- run_mosaicing_prediction_fun(i,list_param_run_mosaicing_prediction)


  		












Also available in:  Unified diff