joshualerickson / explorergee Goto Github PK

Hi,

Been exploring the package and really enjoy the ease with which data can be pulled from gee. Found a possible bug in the band() function for the bands SWIR1 / SWIR2 and temporal = 'year_month' . Below is a quick example. All other combinations appear to work .

library(exploreRGEE)
library(rgee)
library(ggplot2)

ee_Initialize()

huc <- exploreRGEE::huc

ts <- get_landsat(huc, method = 'ts', startDate = '2019-01-01', endDate = '2020-12-31', c.low = 1, c.high = 12)

nothing returned for SWIR1 or SWIR2 for year_month

ts %>% band(band = 'SWIR2', scale = 250, ggplot = TRUE, temporal = 'year_month')

other temporal periods work, e.g.

ts %>% band(band = 'SWIR2', scale = 250, ggplot = TRUE, temporal = 'monthly')

other bands work, e.g.

ts %>% band(band = 'Green', scale = 250, ggplot = TRUE, temporal = 'year_month')

sessionInfo()

R version 4.0.4 (2021-02-15)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 19042)

Matrix products: default

locale:
[1] LC_COLLATE=English_Australia.1252 LC_CTYPE=English_Australia.1252 LC_MONETARY=English_Australia.1252
[4] LC_NUMERIC=C LC_TIME=English_Australia.1252

attached base packages:
[1] stats graphics grDevices utils datasets methods base

other attached packages:
[1] ggplot2_3.3.5.9000 exploreRGEE_0.1 rgee_1.1.3

loaded via a namespace (and not attached):
[1] httr_1.4.2 tidyr_1.2.0 splines_4.0.4 jsonlite_1.8.0 shiny_1.7.1
[6] assertthat_0.2.1 askpass_1.1 sp_1.4-6 pillar_1.7.0 lattice_0.20-41
[11] glue_1.6.2 reticulate_1.24 digest_0.6.29 promises_1.2.0.1 colorspace_2.0-3
[16] htmltools_0.5.2 httpuv_1.6.5 Matrix_1.4-0 pkgconfig_2.0.3 httpcode_0.3.0
[21] purrr_0.3.4 xtable_1.8-4 scales_1.1.1 processx_3.5.0 later_1.3.0
[26] tibble_3.1.6 openssl_2.0.0 proxy_0.4-26 googledrive_2.0.0 mgcv_1.8-39
[31] farver_2.1.0 generics_0.1.2 ellipsis_0.3.2 withr_2.5.0 geojsonio_0.9.4
[36] lazyeval_0.2.2 cli_3.1.0 magrittr_2.0.2 crayon_1.5.1 mime_0.12
[41] maptools_1.1-3 ps_1.6.0 fs_1.5.2 fansi_1.0.3 nlme_3.1-155
[46] foreign_0.8-82 class_7.3-20 geojsonsf_2.0.2 tools_4.0.4 gargle_1.2.0
[51] lifecycle_1.0.1 geojson_0.3.4 stringr_1.4.0 V8_4.1.0 munsell_0.5.0
[56] compiler_4.0.4 e1071_1.7-9 jqr_1.2.3 rlang_1.0.2 classInt_0.4-3
[61] units_0.8-0 grid_4.0.4 rstudioapi_0.13 htmlwidgets_1.5.4 crosstalk_1.2.0
[66] labeling_0.4.2 gtable_0.3.0 DBI_1.1.2 curl_4.3 R6_2.5.1
[71] lubridate_1.8.0 knitr_1.37 dplyr_1.0.8 fastmap_1.1.0 rgeos_0.5-9
[76] utf8_1.2.2 KernSmooth_2.23-20 stringi_1.7.6 crul_1.2.0 Rcpp_1.0.8.3
[81] vctrs_0.4.0 sf_1.0-7 png_0.1-7 leaflet_2.1.1 tidyselect_1.1.2
[86] xfun_0.30

Convert the existing utils_filters.R file to one with modified style from map_date_to_bandname_ic and year_filter2 functions.

Need to fix or somehow adapt 'stdDev' to be able to reduce multiple bands. Right now, it works if param = 'NDVI' one band only, if not then it won't reduce. I think this has to do with using data$reduce(ee$Reducer$stdDev()) in data_stat().

I want to add value to list to ee.Image.

MODDAY_ext_lst <- select(MODDAY_ext, LST_Day_1km)
print(MODDAY_ext_lst)
class(MODDAY_ext_lst)

Result:

class(MODDAY_ext_lst)
[1] "data.frame"

typeof(MODDAY_ext_lst)
[1] "list"

cal_value <- na_interpolation(MODDAY_ext_lst, option = "linear", maxgap = Inf)
class(cal_value)
typeof(cal_value)

Result:

class(MODDAY_ext)
[1] "data.frame"

typeof(MODDAY_ext)
[1] "list"

MODDAY_ext["LST_Day_1km"] <- cal_value["LST_Day_1km"]
print(MODDAY_ext)
class(MODDAY_ext)
typeof(MODDAY_ext)

Result:

class(MODDAY_ext)
[1] "data.frame"

typeof(MODDAY_ext)
[1] "list"

I want to add the created parameter(MODDAY_ext) to 'ee.Image' file. How can I do? Thank you.

need to add TRMMm and TRMMh to banding class for getting Date column.

add @details to get_landsat() in roxygen

Need to change the use of data to image so it's not confusing for the user.