Revision 25d7918b
Added by Benoit Parmentier over 8 years ago
| climate/research/oregon/interpolation/master_script_stage_7.R | ||
|---|---|---|
| 60 | 60 |
|
| 61 | 61 |
#script_path <- "/home/parmentier/Data/IPLANT_project/env_layers_scripts" |
| 62 | 62 |
script_path <- "/nobackupp8/bparmen1/env_layers_scripts" #path to script |
| 63 |
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_12192015.R" #PARAM12
|
|
| 64 |
function_mosaicing <-"global_run_scalingup_mosaicing_04052016.R"
|
|
| 63 |
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_04062016b.R" #PARAM12
|
|
| 64 |
function_mosaicing <-"global_run_scalingup_mosaicing_04062016.R"
|
|
| 65 | 65 |
source(file.path(script_path,function_mosaicing)) #source all functions used in this script |
| 66 | 66 |
source(file.path(script_path,function_mosaicing_functions)) #source all functions used in this script |
| 67 | 67 |
|
| ... | ... | |
| 117 | 117 |
create_out_dir_param <- TRUE #PARAM 12 |
| 118 | 118 |
|
| 119 | 119 |
#if daily mosaics NULL then mosaicas all days of the year #PARAM 13 |
| 120 |
day_to_mosaic <- c("19910101","19910102","19910103") #,"19920104","19920105") #PARAM9, two dates note in /tiles for now on NEX
|
|
| 121 |
|
|
| 120 |
#day_to_mosaic <- c("19910101","19910102","19910103") #,"19920104","19920105") #PARAM9, two dates note in /tiles for now on NEX
|
|
| 121 |
day_to_mosaic_range <- c("19910101","19910103") #if null run all year
|
|
| 122 | 122 |
#CRS_WGS84 <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +towgs84=0,0,0") #Station coords WGS84 #CONSTANT1
|
| 123 | 123 |
#CRS_locs_WGS84<-CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +towgs84=0,0,0") #Station coords WGS84
|
| 124 | 124 |
#proj_str<- CRS_WGS84 #PARAM 8 #check this parameter |
| ... | ... | |
| 166 | 166 |
#rasterOptions(maxmemory=1e+07,timer=TRUE) |
| 167 | 167 |
list_param_run_mosaicing_prediction <- list(in_dir,y_var_name,interpolation_method,region_name, |
| 168 | 168 |
mosaicing_method,out_suffix,out_suffix_str,metric_name,pred_mod_name,var_pred, |
| 169 |
create_out_dir_param,day_to_mosaic,proj_str,file_format,NA_value,num_cores, |
|
| 170 |
region_name,use_autokrige,infile_mask,df_assessment_files_name,mosaic_python,
|
|
| 169 |
create_out_dir_param,day_to_mosaic_range,proj_str,file_format,NA_value,num_cores,
|
|
| 170 |
use_autokrige,infile_mask,df_assessment_files_name,mosaic_python, |
|
| 171 | 171 |
python_bin,algorithm,match_extent,list_models) |
| 172 | 172 |
param_names <- c("in_dir","y_var_name","interpolation_method","region_name",
|
| 173 | 173 |
"mosaicing_method","out_suffix","out_suffix_str","metric_name","pred_mod_name","var_pred", |
| 174 |
"create_out_dir_param","day_to_mosaic","proj_str","file_format","NA_value","num_cores", |
|
| 175 |
"region_name","use_autokrige","infile_mask","df_assessment_files_name","mosaic_python",
|
|
| 174 |
"create_out_dir_param","day_to_mosaic_range","proj_str","file_format","NA_value","num_cores",
|
|
| 175 |
"use_autokrige","infile_mask","df_assessment_files_name","mosaic_python", |
|
| 176 | 176 |
"python_bin","algorithm","match_extent","list_models") |
| 177 | 177 |
names(list_param_run_mosaicing_prediction) <- param_names |
| 178 | 178 |
#list_param_run_mosaicing_prediction |
| ... | ... | |
| 180 | 180 |
#debug(debug_fun_test) |
| 181 | 181 |
#debug_fun_test(list_param_raster_prediction) |
| 182 | 182 |
i <- 1 #this select the first year of list_year_predicted |
| 183 |
function_mosaicing_functions <- "global_run_scalingup_mosaicing_function_04062016b.R" #PARAM12 |
|
| 184 |
function_mosaicing <-"global_run_scalingup_mosaicing_04062016.R" |
|
| 185 |
source(file.path(script_path,function_mosaicing)) #source all functions used in this script |
|
| 186 |
source(file.path(script_path,function_mosaicing_functions)) #source all functions used in this script |
|
| 187 |
|
|
| 183 | 188 |
if (stages_to_run[7]==7){
|
| 184 | 189 |
assessment_prediction_obj <- run_mosaicing_prediction_fun(i,list_param_run_mosaicing_prediction) |
| 185 | 190 |
#list_param_run_mosaicing_prediction |
Also available in: Unified diff
more changes to master script 7, stage 7 and testing