product <- "MOD_NDVI_M"
year <- 2003
month <- 01
ndvi_03_01 <- NASA.neo.gfc.download(product, year, month)

NASA.stack <- stack(ndvi_03_01)
NASA.stack[NASA.stack==255] <- NA
NASA.stack <- NASA.stack/255
NASA.stack_03_01 <- NASA.stack - .1

temp <- raster(ndvi_03_01[1])
cols <- temp@legend@colortable
rm(temp)

plot(NASA.stack_03_01, ext = extent, col = cols, zlim = -.1,.9)

vals_03_01 <- extract(NASA.stack_03_01, shp)