# NOT RUN {
# Retrieve public earth temperature data for the last 30 days at Reefton Ews
# station, at a depth of 10cm
# Subtract 30 days from today's date to get the start date
last_month = paste(as.character(Sys.Date() - 30), 0)
reefton_earth = cf_query(cf_user(), cf_datatype(4, 3, 2), cf_station(),
start_date = last_month)
class(reefton_earth) # cfTemp object
# Plot the temperature data using the defaults
plot(reefton_earth)
# Enlarge the text and add the observations as points
library(ggplot2) # for element_text() and geom_point()
plot(reefton_earth, ggtheme = "bw", text = element_text(size = 16)) +
geom_point(size = 3, shape = 1)
# Save the plot as a png to the current working directory
library(ggplot2) # for ggsave()
ggsave("my_earthTemp_plot.png")
# }
Run the code above in your browser using DataLab