Revision 4d0fc33b
Added by Benoit Parmentier about 8 years ago
climate/research/oregon/interpolation/global_product_assessment_part0_functions.R | ||
---|---|---|
204 | 204 |
tile_coord <- names(list_spdf)[i] |
205 | 205 |
r_ref <- list_r_ref[[i]] |
206 | 206 |
|
207 |
tile_id <- paste0("tile_",i) |
|
207 | 208 |
df_tmp <- subset(df_missing,date==date_val,select=tile_coord) |
208 | 209 |
#for each row (date) |
209 | 210 |
val <- df_tmp[[tile_coord]] |
... | ... | |
216 | 217 |
tile_spdf$predicted <- val |
217 | 218 |
tile_spdf$tile_coord <- tile_coord |
218 | 219 |
tile_spdf$overlap <- 1 |
220 |
tile_spdf$tile_id <- tile_id |
|
219 | 221 |
|
220 | 222 |
#r <- rasterize(tile_spdf,r_ref,"predicted") |
221 | 223 |
#r <- rasterize(tile_spdf,r_ref,col_name) |
... | ... | |
231 | 233 |
r <- init(r, fun=set1f, overwrite=TRUE) |
232 | 234 |
} |
233 | 235 |
|
234 |
return(r) |
|
236 |
out_dir <- "." #can set this up later |
|
237 |
out_suffix_str <- paste0(col_name,"") # can set this parameter later |
|
238 |
raster_name <- file.path(out_dir,paste("r_",tile_id,"_",tile_coord,"_",out_suffix_str,".tif",sep="")) |
|
239 |
#raster_name <- |
|
240 |
writeRaster(r, NAflag=NA_flag_val,filename=raster_name,overwrite=TRUE) #unweighted mean |
|
241 |
|
|
242 |
return(raster_name) |
|
235 | 243 |
} |
236 | 244 |
|
237 | 245 |
############################ END OF SCRIPT ################################## |
Also available in: Unified diff
fixing generate raster function from tile