Project

General

Profile

« Previous | Next » 

Revision 578abd5d

Added by Benoit Parmentier over 8 years ago

debugging code to find out issues for missing days for mosaicing for year 1999

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/01/2016  


  		




  
  17
  
    
#MODIFIED ON: 05/30/2016  


  		




  18
  18
  
    
#PROJECT: NCEAS INPLANT: Environment and Organisms                                                                           


  		




  19
  19
  
    

  		




  20
  20
  
    
#First source these files:


  		




  ......




  63
  63
  
    

  		




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


  		




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


  		




  66
  
  
    
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_05012016.R" #PARAM12


  		




  67
  
  
    
function_mosaicing <-"global_run_scalingup_mosaicing_05012016.R"


  		




  
  66
  
    
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_05302016.R" #PARAM12


  		




  
  67
  
    
function_mosaicing <-"global_run_scalingup_mosaicing_05302016.R"


  		




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


  		




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


  		




  70
  70
  
    

  		




  ......




  106
  106
  
    
max_mem = 1e+07 #param 9, arg 9


  		




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


  		




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


  		




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


  		




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


  		




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


  		




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


  		




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


  		




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


  		




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


  		




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


  		




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


  		




  ......




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


  		




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


  		




  120
  122
  
    
values_range <- "-100,100"


  		




  
  123
  
    
NA_flag_val <- -32768 #should be here


  		




  121
  124
  
    

  		




  122
  125
  
    
# var <- "TMAX" # variable being interpolated #param 1, arg 1


  		




  123
  126
  
    
# in_dir <- "/nobackupp6/aguzman4/climateLayers/out/" #PARAM2,arg 2


  		




  ......




  281
  284
  
    
#debug(debug_fun_test)


  		




  282
  285
  
    
#debug_fun_test(list_param_raster_prediction)


  		




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


  		




  284
  
  
    
#function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_04072016.R" #PARAM12


  		




  285
  
  
    
#function_mosaicing <-"global_run_scalingup_mosaicing_04072016.R"


  		




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


  		




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


  		




  
  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 


  		




  288
  291
  
    

  		




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


  		




  290
  293
  
    
  assessment_prediction_obj <- run_mosaicing_prediction_fun(i,list_param_run_mosaicing_prediction)


  		












Also available in:  Unified diff