Revision 07434086
Added by Benoit Parmentier over 8 years ago
| climate/research/oregon/interpolation/global_run_scalingup_mosaicing.R | ||
|---|---|---|
| 418 | 418 |
"var_pred","list_models","use_autokrige","y_var_name","interpolation_method", |
| 419 | 419 |
"days_to_process","num_cores_tmp","NA_flag_val","file_format","out_dir_str", |
| 420 | 420 |
"out_suffix_str") |
| 421 |
|
|
| 421 |
#browser() |
|
| 422 | 422 |
list_create_accuracy_residuals_raster_obj <- lapply(1:length(day_to_mosaic),FUN=create_accuracy_residuals_raster, |
| 423 | 423 |
list_param=list_param_create_accuracy_residuals_raster) |
| 424 | 424 |
|
| ... | ... | |
| 426 | 426 |
#list_create_accuracy_residuals_raster_obj <- lapply(1:1,FUN=create_accuracy_residuals_raster, |
| 427 | 427 |
# list_param=list_param_create_accuracy_residuals_raster) |
| 428 | 428 |
|
| 429 |
#create_accuracy_residuals_raster_obj <- create_accuracy_metric_raster(1, list_param_create_accuracy_residuals_raster_obj)
|
|
| 429 |
#create_accuracy_residuals_raster_obj <- create_accuracy_residuals_raster(1, list_param_create_accuracy_residuals_raster_obj)
|
|
| 430 | 430 |
|
| 431 | 431 |
#note that three tiles did not produce a residuals surface!!! find out more later, join the output |
| 432 | 432 |
#to df_raste_tile to keep track of which one did not work... |
| ... | ... | |
| 436 | 436 |
#Plot as quick check |
| 437 | 437 |
#r1 <- raster(lf_mosaic[[1]][1]) |
| 438 | 438 |
#list_create_accuracy_residuals_raster_obj |
| 439 |
|
|
| 439 |
#browser() |
|
| 440 |
|
|
| 441 |
######################################### |
|
| 440 | 442 |
##Run for data_day_s |
| 441 | 443 |
## |
| 442 | 444 |
data_df <- data_day_s # data.frame table/spdf containing stations with residuals and variable |
| ... | ... | |
| 477 | 479 |
"var_pred","list_models","use_autokrige","y_var_name","interpolation_method", |
| 478 | 480 |
"days_to_process","num_cores_tmp","NA_flag_val","file_format","out_dir_str", |
| 479 | 481 |
"out_suffix_str") |
| 480 |
|
|
| 482 |
browser() |
|
| 481 | 483 |
list_create_accuracy_residuals_raster_obj <- lapply(1:length(day_to_mosaic),FUN=create_accuracy_residuals_raster, |
| 482 | 484 |
list_param=list_param_create_accuracy_residuals_raster) |
| 483 | 485 |
|
Also available in: Unified diff
mosaicing script, testing residuals generation function for raster surface