forked from impact-initiatives-geospatial/surveyGEER
-
Notifications
You must be signed in to change notification settings - Fork 0
/
reprex_20220929.R
57 lines (43 loc) · 1.75 KB
/
reprex_20220929.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
library(rgee)
library(tidyrgee)
ee_Initialize()
## 1 reprex posted
# works fine with rgee
gldas <- ee$ImageCollection("NASA/GLDAS/V021/NOAH/G025/T3H")
sm <- gldas$select("SoilMoi0_10cm_inst")
sm_filt <- sm$filterDate("2004-01-01","2009-12-31")
sm_2004_2009_mean <- sm$filterDate("2004-01-01","2009-12-31")$mean()
sm_2004_2009_mean$bandNames()$getInfo()
# tidyversion throws error
sm_tidy <- as_tidyee(sm)
sm_2004_2009_mean_tidy<- sm_tidy|>
filter(year %in% c(2004:2009)) |>
summarise(stat="mean")
sm_2004_2009_mean_tidy$ee_ob$bandNames()$getInfo()
# 2 is it reducer? - no this works
gldas <- ee$ImageCollection("NASA/GLDAS/V021/NOAH/G025/T3H")
sm <- gldas$select("SoilMoi0_10cm_inst")
sm_2004_2009_filtered <- sm$filterDate("2004-01-01","2009-12-31")
ee_reducer <- tidyrgee:::stat_to_reducer_full("mean")
sm_2004_2009_reduced <- ee_reducer(sm_2004_2009_filtered)
sm_2004_2009_reduced$bandNames()$getInfo()
# 3 reprex attempt -1b - this works so not the problem
gldas <- ee$ImageCollection("NASA/GLDAS/V021/NOAH/G025/T3H")
sm <- gldas$select("SoilMoi0_10cm_inst")
sm_tidy <- as_tidyee(sm)
sm_tidy$ee_ob$first()$bandNames()$getInfo()
# 4
#debugonce(tidyrgee:::summarise_pixels)
# 5 - try another IC -- chirps
chirps <- ee$ImageCollection("UCSB-CHG/CHIRPS/DAILY")
chirps <- chirps$select("precipitation")
chirps_filt <- chirps$filterDate("2004-01-01","2009-12-31")
chirps_2004_2009_mean <- chirps$filterDate("2004-01-01","2009-12-31")$mean()
chirps_2004_2009_mean$bandNames()$getInfo()
# Here it works with CHIRPS, but it does seem to
# have a considerable lag compared to above
chirps_tidy <- as_tidyee(chirps)
chirps_2004_2009_mean_tidy<- chirps_tidy|>
filter(year %in% c(2004:2009)) |>
summarise(stat="mean")
chirps_2004_2009_mean_tidy$ee_ob$bandNames()$getInfo()