Revision 7cd6d106
Added by Benoit Parmentier over 11 years ago
| climate/research/oregon/interpolation/GAM_fusion_function_multisampling.R | ||
|---|---|---|
| 8 | 8 |
# 5)runGAMFusion <- function(i,list_param) : daily step for fusion method, perform daily prediction |
| 9 | 9 |
# |
| 10 | 10 |
#AUTHOR: Benoit Parmentier |
| 11 |
#DATE: 04/20/2013
|
|
| 11 |
#DATE: 04/30/2013
|
|
| 12 | 12 |
#PROJECT: NCEAS INPLANT: Environment and Organisms --TASK#363-- |
| 13 | 13 |
|
| 14 | 14 |
##Comments and TODO: |
| ... | ... | |
| 258 | 258 |
#j indicate which month is predicted, var indicates TMIN or TMAX |
| 259 | 259 |
data_name<-paste(var,"_bias_LST_month_",j,"_",cname[k],"_",prop_month, |
| 260 | 260 |
"_",run_samp,sep="") |
| 261 |
raster_name<-paste("fusion_",data_name,out_prefix,".tif", sep="")
|
|
| 261 |
raster_name<-paste("fusion",data_name,out_prefix,".tif", sep="")
|
|
| 262 | 262 |
list_out_filename[[k]]<-raster_name |
| 263 | 263 |
} |
| 264 | 264 |
|
| ... | ... | |
| 472 | 472 |
#Saving kriged surface in raster images |
| 473 | 473 |
data_name<-paste("daily_delta_",y_var_name,"_",sampling_dat$date[i],"_",sampling_dat$prop[i],
|
| 474 | 474 |
"_",sampling_dat$run_samp[i],sep="") |
| 475 |
raster_name_delta<-paste("fusion_",var,"_",data_name,out_prefix,".tif", sep="")
|
|
| 475 |
raster_name_delta<-paste(interpolation_method,"_",var,"_",data_name,out_prefix,".tif", sep="")
|
|
| 476 | 476 |
writeRaster(daily_delta_rast, filename=raster_name_delta,overwrite=TRUE) #Writing the data in a raster file format...(IDRISI) |
| 477 | 477 |
|
| 478 | 478 |
#Now predict daily after having selected the relevant month |
Also available in: Unified diff
gam cai funciton test run