Revision b516681e
Added by Benoit Parmentier over 11 years ago
climate/research/oregon/interpolation/master_script_temp.R | ||
---|---|---|
80 | 80 |
#met_stations_outfiles_obj_file<-"met_stations_outfiles_obj_kriging_daily__365d_kriging_daily_mults10_lst_comb3_08062013.RData" |
81 | 81 |
|
82 | 82 |
var<-"TMAX" # variable being interpolated |
83 |
out_prefix<-"_365d_gam_cai_lst_comb3_08122013" #User defined output prefix
|
|
83 |
out_prefix<-"_365d_kriging_fus_lst_comb3_08122013" #User defined output prefix
|
|
84 | 84 |
out_suffix<-"_OR_08122013" #Regional suffix |
85 | 85 |
out_suffix_modis <-"_05302013" #pattern to find tiles produced previously |
86 | 86 |
|
87 | 87 |
#interpolation_method<-c("gam_fusion","gam_CAI","gam_daily") #other otpions to be added later |
88 |
interpolation_method<-c("gam_CAI") #other otpions to be added later |
|
88 |
#interpolation_method<-c("gam_CAI") #other otpions to be added later
|
|
89 | 89 |
#interpolation_method<-c("gam_fusion") #other otpions to be added later |
90 |
#interpolation_method<-c("kriging_fusion") #other otpions to be added later
|
|
90 |
interpolation_method<-c("kriging_fusion") #other otpions to be added later |
|
91 | 91 |
#interpolation_method<-c("gwr_fusion") #other otpions to be added later |
92 | 92 |
#interpolation_method<-c("gwr_CAI") #other otpions to be added later |
93 | 93 |
#interpolation_method<-c("kriging_CAI") |
... | ... | |
256 | 256 |
#LC1: Evergreen/deciduous needleleaf trees |
257 | 257 |
|
258 | 258 |
#Combination 3: for paper baseline=s(lat,lon)+s(elev) |
259 |
list_models<-c("y_var ~ s(lat,lon) + s(elev_s)", |
|
260 |
"y_var ~ s(lat,lon) + s(elev_s) + s(N_w)", |
|
261 |
"y_var ~ s(lat,lon) + s(elev_s) + s(E_w)", |
|
262 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST)", |
|
263 |
"y_var ~ s(lat,lon) + s(elev_s) + s(DISTOC)", |
|
264 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LC1)", |
|
265 |
"y_var ~ s(lat,lon) + s(elev_s) + s(CANHGHT)", |
|
266 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)", |
|
267 |
"y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,CANHGHT)") |
|
259 |
# list_models<-c("y_var ~ s(lat,lon) + s(elev_s)",
|
|
260 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(N_w)",
|
|
261 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(E_w)",
|
|
262 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST)",
|
|
263 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(DISTOC)",
|
|
264 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LC1)",
|
|
265 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(CANHGHT)",
|
|
266 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,LC1)",
|
|
267 |
# "y_var ~ s(lat,lon) + s(elev_s) + s(LST) + ti(LST,CANHGHT)")
|
|
268 | 268 |
|
269 | 269 |
#list_models<-c("y_var ~ lat*lon + elev_s") |
270 | 270 |
|
271 | 271 |
#list_models<-c("y_var ~ s(lat,lon) + s(elev_s)") |
272 | 272 |
|
273 |
#list_models<-c("y_var ~ lat*lon + elev_s",
|
|
274 |
# "y_var ~ lat*lon + elev_s + N_w",
|
|
275 |
# "y_var ~ lat*lon + elev_s + E_w",
|
|
276 |
# "y_var ~ lat*lon + elev_s + LST",
|
|
277 |
# "y_var ~ lat*lon + elev_s + DISTOC",
|
|
278 |
# "y_var ~ lat*lon + elev_s + LC1",
|
|
279 |
# "y_var ~ lat*lon + elev_s + CANHGHT",
|
|
280 |
# "y_var ~ lat*lon + elev_s + LST + I(LST*LC1)",
|
|
281 |
# "y_var ~ lat*lon + elev_s + LST + I(LST*CANHGHT)")
|
|
273 |
list_models<-c("y_var ~ lat*lon + elev_s", |
|
274 |
"y_var ~ lat*lon + elev_s + N_w",
|
|
275 |
"y_var ~ lat*lon + elev_s + E_w",
|
|
276 |
"y_var ~ lat*lon + elev_s + LST",
|
|
277 |
"y_var ~ lat*lon + elev_s + DISTOC",
|
|
278 |
"y_var ~ lat*lon + elev_s + LC1",
|
|
279 |
"y_var ~ lat*lon + elev_s + CANHGHT",
|
|
280 |
"y_var ~ lat*lon + elev_s + LST + I(LST*LC1)",
|
|
281 |
"y_var ~ lat*lon + elev_s + LST + I(LST*CANHGHT)")
|
|
282 | 282 |
|
283 | 283 |
#Default name of LST avg to be matched |
284 | 284 |
lst_avg<-c("mm_01","mm_02","mm_03","mm_04","mm_05","mm_06","mm_07","mm_08","mm_09","mm_10","mm_11","mm_12") |
Also available in: Unified diff
running kriging fusion comb3 OR with revised screening for paper analyses