convert static plotting from base r to ggplot2 for episode 12 and 13 #190
Conversation
| plot(NDVI_HARV_stack, | ||
| zlim = c(.15, 1), | ||
| nc = 4) | ||
| # convert to a df for plotting in two steps, |
There was a problem hiding this comment.
These lines look identical to lines 217-225. Can they be removed or is there a difference I'm not seeing?
There was a problem hiding this comment.
Yes it is a bit heavy but because we changed NDVI_HARV_stack (divided by 10000), we need to recreate NDVI_HARV_stack _df again, isn't it?
|
Thanks @annefou, this looks really good. I left some notes - one place where there appears to be about 10 lines of duplicated code, and then a few places where I think we can shorten the ggplot calls. I'm happy to merge this as is (once the check completes) and make further fine-tuning adjustments later, so don't worry about following up on my review. I can test out the adjustments and see if they work. Thanks again! |
| # one single column called 'variable' | ||
| melt(id.vars = c('x','y','optional')) | ||
|
|
||
| # view a plot of all of the rasters |
There was a problem hiding this comment.
This also looks identical to lines 227-232
| > `par(mfrow=c(2, 2))` to create a 2x2 tiled layout. When you are done, be sure to | ||
| > reset your layout using: `par(mfrow=c(1, 1))`. | ||
| > | ||
| > `grid.arrange` from gridExtra R package. Make sure you load it by using: |
There was a problem hiding this comment.
I'm not clear on why we're introducing a new library in an exercise. I think that will be really, really challenging for people. Is this something that you can't do with facet_grid?
| Once the names for each band have been reassigned, we can render our plot with | ||
| the new labels using `names.attr=rasterNames`. | ||
| the new labels using a`labeller`. |
| ggtitle("Landsat NDVI - Julian Days \nHarvard Forest 2011") + | ||
| theme(plot.title = element_text(hjust = 0.5), | ||
| # remove x and y-axis title, text and ticks | ||
| axis.title.x=element_blank(), |
There was a problem hiding this comment.
I think we can do this by using theme_void() on line 207 (before the theme call to center the title). Then we can get rid of lines 209-214 (I think).
There was a problem hiding this comment.
Or possibly theme_blank(). https://gist.github.com/dgrtwo/38e70de658b48e166f90
There was a problem hiding this comment.
Yes if we use:
theme_void() + theme(plot.title = element_text(hjust = 0.5),
we can remove the rest.
|
@annefou - the check passed! I'm going to merge this now and we can make some minor adjustments as needed in a later PR. |
convert static plotting from base r to ggplot2 for episode 12 and 13
updated episode 12 and 13: ggplot2 is now used instead of base r and rasterVis as requested in #133