MODIS arctic sites data request

MODIS_arctic_sites_download.R

# download NDVI time series
library(MODISTools)

df <- data.frame(
  site_name = 
    c("US-Bes",
      "US-Atq",
      "US-Ivo",
      "US-ICt",
      "DK-ZaH",
      "RU-Che",
      "RU-Cok",
      "RU-Sam",
      "CA-DL1",
      "CA-TVC"),
  lat = c(
    71.280881,
    70.4696228,
    68.48649,
    68.61667,
    74.4732,
    68.61689,
    70.82973,
    72.3733,
    64.8688553,
    68.74617),
  lon = c(
    -156.596467,
    -157.4089471,
    -155.75022,
    -149.3,
    20.5503,
    161.35089,
    147.48897,
    126.4978,
    -111.5747927,
    -133.50171
  )
  )


modis_data <- apply(df, 1, function(x){
  mt_subset(
    site_name = x['site_name'],
    lat = x['lat'],
    lon = x['lon'],
    product = "MOD13Q1",
    band = c("250m_16_days_EVI",
             "250m_16_days_NDVI",
             "250m_16_days_pixel_reliability"),
    km_lr = 1,
    km_ab = 1,
    internal = TRUE
  )  
})

data request