inst/doc/VII_combine_raster_point.R

In rnpn: Interface to the National 'Phenology' Network 'API'

## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(echo = TRUE)
library(RColorBrewer)
library(rnpn)
library(raster)

## ----eval=FALSE---------------------------------------------------------------
#  npn_download_site_phenometrics(
#    request_source = 'Your Name Here',
#    years = '2013',
#    num_days_quality_filter = '30',
#    species_ids = '35',
#    phenophase_ids = '373',
#    download_path = 'cl_lilac_data_2013_SIxLeaf.csv',
#    six_leaf_layer = TRUE,
#    six_sub_model = 'lilac'
#  )

## ----eval=FALSE---------------------------------------------------------------
#  LilacLeaf2018<-npn_download_geospatial(
#    'si-x:lilac_leaf_ncep_historic',
#    '2018-01-01',
#  )
#  
#  
#  LilacLeaf2018Obs <-npn_download_individual_phenometrics(
#    request_source = 'Your Name Here',
#    years = '2018',
#    species_ids = '35',
#    phenophase_ids = '373'
#  )
#  
#  

## ----eval=FALSE---------------------------------------------------------------
#  
#  coords <- LilacLeaf2018Obs[ , c("longitude", "latitude")]
#  data <- as.data.frame(LilacLeaf2018Obs$first_yes_doy)
#  
#  crs <- CRS("+proj=utm +zone=18 +datum=WGS84 +units=m +no_defs
#                   +ellps=WGS84 +towgs84=0,0,0")
#  
#  LL_spdf <- SpatialPointsDataFrame(
#    coords = coords,
#    data = data,
#    proj4string = crs
#  )
#  

## ----eval=FALSE---------------------------------------------------------------
#  
#  my.palette <- brewer.pal(n=9,name="OrRd")
#  
#  plot(
#    LilacLeaf2018,
#    col = my.palette,
#    main="2018 Observed and Predicted Lilac Leaf Out"
#  )
#  
#  plot(
#    LL_spdf,
#    main="Lilac Obs",
#    pch = 21,
#    bg = my.palette,
#    col = 'black',
#    xlim=c(-125.0208,-66.47917),
#    ylim=c(24.0625 ,49.9375),
#    add = TRUE
#  )
#  
#  legend(
#    "bottomright",
#    legend=c("Cloned Lilac Leaf Out Observations"),
#    pch = 21,
#    bg = 'white',
#    col = 'black',
#    bty="n",
#    cex=.8
#  )
#  

## ---- echo=FALSE, out.width = "50%", fig.pos="h"------------------------------
knitr::include_graphics("figures/7-plot.png")