Revision 96cf6457
Added by Benoit Parmentier over 8 years ago
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: 08/16/2016
|
|
17 |
#MODIFIED ON: 08/22/2016
|
|
18 | 18 |
#PROJECT: NCEAS INPLANT: Environment and Organisms |
19 | 19 |
|
20 | 20 |
#First source these files: |
... | ... | |
87 | 87 |
script_path <- "/nobackupp8/bparmen1/env_layers_scripts" #path to script |
88 | 88 |
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_08222016.R" |
89 | 89 |
#function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_07052016.R" #PARAM12 |
90 |
function_mosaicing <-"global_run_scalingup_mosaicing_08162016.R"
|
|
90 |
function_mosaicing <-"global_run_scalingup_mosaicing_08222016.R"
|
|
91 | 91 |
source(file.path(script_path,function_mosaicing)) #source all functions used in this script |
92 | 92 |
source(file.path(script_path,function_mosaicing_functions)) #source all functions used in this script |
93 | 93 |
|
... | ... | |
129 | 129 |
num_cores <- 6 #PARAM 8, arg 8 |
130 | 130 |
max_mem = 1e+07 #param 9, arg 9 |
131 | 131 |
mosaicing_method <- "use_edge_weights" #PARAM10, arg 10 |
132 |
metric_name <- "rmse" # "mae", "r" for MAE, R etc.; can also be ns or nv? #PARAM 11, arg 11 |
|
133 |
#metric_name <- "n"
|
|
132 |
#metric_name <- "rmse" # "mae", "r" for MAE, R etc.; can also be ns or nv? #PARAM 11, arg 11
|
|
133 |
metric_name <- "n" |
|
134 | 134 |
#metric_name <- "mae" |
135 | 135 |
|
136 | 136 |
#day_start <- "19990701" #PARAM 12 arg 12 |
... | ... | |
143 | 143 |
#df_assessment_files_name <- "/nobackupp6/aguzman4/climateLayers/out/reg5/assessment/output_reg5_1991/df_assessment_files_reg5_1991_reg5_1991.txt" # data.frame with all files used in assessmnet, PARAM 15 |
144 | 144 |
#df_assessment_files_name <- "/nobackupp6/aguzman4/climateLayers/out/reg5/assessment/output_reg5_1985/df_assessment_files_reg5_1985_reg5_1985.txt" |
145 | 145 |
algorithm <- "python" #PARAM 16 #if R use mosaic function for R, if python use modified gdalmerge script from Alberto Guzmann |
146 |
layers_option <- c("var_pred") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred |
|
146 |
#layers_option <- c("var_pred") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred
|
|
147 | 147 |
#layers_option <- c("ac_training") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred |
148 | 148 |
#layers_option <- c("res_training") # #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred |
149 | 149 |
#layers_option <- c("res_testing") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred |
150 |
#layers_option <- c("ac_testing") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred
|
|
150 |
layers_option <- c("ac_testing") #arg 17 ,param 17, options are:#res_training, res_testing,ac_training, ac_testing, var_pred |
|
151 | 151 |
|
152 | 152 |
tmp_files <- FALSE #arg 18, param 18, keep temp files if TRUE |
153 | 153 |
data_type <- "Int16" #, param 19, use int32 for output layers mosaiced |
154 | 154 |
scaling <- 100 #, param 20, if null use 1 |
155 |
values_range <- "-100,100" |
|
155 |
scaling <- 1 #use this if predicting n rather than other variables |
|
156 |
#values_range <- "-100,100" |
|
157 |
values_range <- "0,32767" #this is for n variable |
|
156 | 158 |
NA_flag_val <- -32768 #should be here |
157 | 159 |
|
158 | 160 |
|
Also available in: Unified diff
more testing of stage 7 with observation maps of n