if (FALSE) {
library(rgee)
library(tidyrgee)
ee_Initizialize()
modis_ic <- ee$ImageCollection("MODIS/006/MOD13Q1")
point_sample_buffered <- tidyrgee::bgd_msna |>
sample_n(3) |>
sf::st_as_sf(coords=c("_gps_reading_longitude",
"_gps_reading_latitude"), crs=4326) |>
sf::st_transform(crs=32646) |>
sf::st_buffer(dist = 500) |>
dplyr::select(`_uuid`)
modis_ic_tidy <- as_tidyee(modis_ic)
modis_monthly_baseline_mean <- modis_ic_tidy |>
select("NDVI") |>
filter(year %in% 2000:2015) |>
group_by(month) |>
summarise(stat="mean")
ndvi_monthly_mean_at_pt<- modis_monthly_baseline_mean |>
ee_extract(y = point_sample_buffered,
fun="mean",
scale = 500)
}
Run the code above in your browser using DataLab