Revision 2411ddd5
Added by Benoit Parmentier over 8 years ago
| climate/research/oregon/interpolation/global_product_assessment_part0.R | ||
|---|---|---|
| 9 | 9 |
# |
| 10 | 10 |
#AUTHOR: Benoit Parmentier |
| 11 | 11 |
#CREATED ON: 10/27/2016 |
| 12 |
#MODIFIED ON: 11/06/2016
|
|
| 12 |
#MODIFIED ON: 11/07/2016
|
|
| 13 | 13 |
#Version: 1 |
| 14 | 14 |
#PROJECT: Environmental Layers project |
| 15 | 15 |
#COMMENTS: |
| ... | ... | |
| 385 | 385 |
list_r_ref=list_r_ref, |
| 386 | 386 |
col_name="overlap", |
| 387 | 387 |
date_val=df_missing$date[1]) |
| 388 |
##check that everything is correct: |
|
| 389 |
plot(r_mask) |
|
| 390 |
plot(list_predicted,add=T) |
|
| 391 |
plot(spdf_tiles_test,add=T,border="green",usePolypath = FALSE) #added usePolypath following error on brige and NEX |
|
| 392 |
|
|
| 393 |
list_predicted <- mclapply(1:6, |
|
| 394 |
FUN=rasterize_tile_day, |
|
| 395 |
list_spdf=shps_tiles, |
|
| 396 |
df_missing=df_missing, |
|
| 397 |
list_r_ref=list_r_ref, |
|
| 398 |
col_name = "overlap", |
|
| 399 |
date_val=df_missing$date[1], |
|
| 400 |
mc.preschedule=FALSE, |
|
| 401 |
mc.cores = num_cores) |
|
| 388 |
#list_predicted <- mclapply(1:6, |
|
| 389 |
# FUN=rasterize_tile_day, |
|
| 390 |
# list_spdf=shps_tiles, |
|
| 391 |
# df_missing=df_missing, |
|
| 392 |
# list_r_ref=list_r_ref, |
|
| 393 |
# col_name = "overlap", |
|
| 394 |
# date_val=df_missing$date[1], |
|
| 395 |
# mc.preschedule=FALSE, |
|
| 396 |
# mc.cores = num_cores) |
|
| 402 | 397 |
|
| 403 | 398 |
list_predicted <- mclapply(1:length(shps_tiles), |
| 404 | 399 |
FUN=rasterize_tile_day, |
| ... | ... | |
| 409 | 404 |
date_val=df_missing$date[1], |
| 410 | 405 |
mc.preschedule=FALSE, |
| 411 | 406 |
mc.cores = num_cores) |
| 407 |
|
|
| 408 |
##check that everything is correct: |
|
| 409 |
plot(r_mask) |
|
| 410 |
plot(raster(list_predicted[[1]]),add=T) |
|
| 411 |
plot(spdf_tiles_test,add=T,border="green",usePolypath = FALSE) #added usePolypath following error on brige and NEX |
|
| 412 |
|
|
| 412 | 413 |
|
| 413 | 414 |
### Make a list of file |
| 414 | 415 |
out_suffix_str_tmp <- paste0(region_name,"_",out_suffix) |
| ... | ... | |
| 435 | 436 |
out_mosaic_name=out_mosaic_name_predicted_m, |
| 436 | 437 |
raster_ref_name = rast_ref_name) ##if NA, not take into account |
| 437 | 438 |
r_overlap_raster_name <- mosaic_overlap_tiles_obj$out_mosaic_name |
| 438 |
cmd_str1 <- mosaic_weights_obj$cmd_str
|
|
| 439 |
cmd_str1 <- mosaic_overlap_tiles_obj$cmd_str
|
|
| 439 | 440 |
|
| 440 | 441 |
r_overlap <- raster(r_overlap_raster_name) |
| 441 | 442 |
r_mask <- raster(infile_mask) |
| ... | ... | |
| 443 | 444 |
out_mosaic_name_overlap_masked <- file.path(out_dir_str,paste("r_overlap_sum_masked_",out_suffix_str_tmp,".tif",sep=""))
|
| 444 | 445 |
|
| 445 | 446 |
r_overlap_m <- mask(r_overlap,r_mask,filename=out_mosaic_name_overlap_masked) |
| 447 |
plot(r_overlap_m) |
|
| 448 |
#plot(spdf_tiles_test,add=T,border="green",usePolypath = FALSE) #added usePolypath following error on brige and NEX |
|
| 446 | 449 |
|
| 450 |
|
|
| 451 |
r_table <- ratify(r_overlap_m) # build the Raster Attibute table |
|
| 452 |
rat <- levels(r_table)[[1]]#get the values of the unique cell frot the attribute table |
|
| 453 |
#rat$legend <- paste0("tile_",1:26)
|
|
| 454 |
tb_freq <- as.data.frame(freq(r_table)) |
|
| 455 |
rat$legend <- tb_freq$value |
|
| 456 |
|
|
| 457 |
levels(r_table) <- rat |
|
| 458 |
#my_col=c('blue','red','green')
|
|
| 459 |
my_col <- rainbow(length(tb_freq$value)) |
|
| 460 |
|
|
| 461 |
plot(r_table,col=my_col,legend=F,box=F,axes=F) |
|
| 462 |
legend(x='topright', legend =rat$legend,fill = my_col,cex=0.8) |
|
| 463 |
|
|
| 464 |
r <- raster(matrix(runif(20),5,4)) |
|
| 465 |
r[r>.5] <- NA |
|
| 466 |
g <- as(r, 'SpatialGridDataFrame') |
|
| 467 |
p <- as(r, 'SpatialPixels') |
|
| 468 |
p_spdf <- as(r_overlap_m,"SpatialPointsDataFrame") |
|
| 469 |
plot(r) |
|
| 470 |
points(p) |
|
| 447 | 471 |
#http://stackoverflow.com/questions/19586945/how-to-legend-a-raster-using-directly-the-raster-attribute-table-and-displaying |
| 448 | 472 |
# |
| 449 | 473 |
#http://gis.stackexchange.com/questions/148398/how-does-spatial-polygon-over-polygon-work-to-when-aggregating-values-in-r |
Also available in: Unified diff
adding plotting of overlap raster and testing rasterize day function