Error caused by changes in default parameters of opgd in sesu_opgd

[SpatLyu]

library(terra)
library(tidyverse)
library(gdverse)
fvcpath = "https://github.com/SpatLyu/rdevdata/raw/main/FVC.tif"
fvc = terra::rast(paste0("/vsicurl/",fvcpath))
fvc
## class       : SpatRaster 
## dimensions  : 418, 568, 13  (nrow, ncol, nlyr)
## resolution  : 1000, 1000  (x, y)
## extent      : -92742.16, 475257.8, 3591385, 4009385  (xmin, xmax, ymin, ymax)
## coord. ref. : Asia_North_Albers_Equal_Area_Conic 
## source      : FVC.tif 
## names       :       fvc,   premax,   premin,   presum,    tmpmax,     tmpmin, ... 
## min values  : 0.1363270, 109.8619,  2.00000, 3783.904,  9.289694, -11.971293, ... 
## max values  : 0.9596695, 249.9284, 82.74928, 8549.112, 26.781870,   1.322163, ...

The original data resolution is 1000m, and then we construct the data under 2000-10000 m spatial units with 1000 spatial unit interval.

su = seq(1000,10000,by = 1000)
fvc1000 = tibble::as_tibble(terra::as.data.frame(fvc,na.rm = T))
fvc_other = 2:10 %>%
  purrr::map(\(.x) terra::aggregate(fvc,fact=.x ,fun="mean") %>%
               terra::as.data.frame(na.rm = T) %>%
               tibble::as_tibble())
fvc = c(list(fvc1000),fvc_other)
str(fvc)
## List of 10
##  $ : tibble [136,243 × 13] (S3: tbl_df/tbl/data.frame)
##   ..$ fvc   : num [1:136243] 0.198 0.193 0.192 0.189 0.208 ...
##   ..$ premax: num [1:136243] 163 161 160 159 164 ...
##   ..$ premin: num [1:136243] 7.95 6.8 5.24 5 9.98 ...
##   ..$ presum: num [1:136243] 3956 3892 3842 3808 4051 ...
##   ..$ tmpmax: num [1:136243] 20.8 20.7 20.9 21.1 20.6 ...
##   ..$ tmpmin: num [1:136243] -7.53 -7.55 -7.48 -7.39 -7.59 ...
##   ..$ tmpavg: num [1:136243] 8.05 8.02 8.15 8.35 7.97 ...
##   ..$ pop   : num [1:136243] 1.903 1.203 0.547 0.542 10.392 ...
##   ..$ ntl   : num [1:136243] 6.6 4.91 3.75 3.99 7.1 ...
##   ..$ lulc  : num [1:136243] 10 10 10 10 10 10 10 10 10 10 ...
##   ..$ elev  : num [1:136243] 1758 1754 1722 1672 1780 ...
##   ..$ slope : num [1:136243] 2.65 3.45 3.96 2.9 1.94 ...
##   ..$ aspect: num [1:136243] 176.4 169.6 138.5 110.9 99.5 ...
##  $ : tibble [33,722 × 13] (S3: tbl_df/tbl/data.frame)
##   ..$ fvc   : num [1:33722] 0.195 0.184 0.157 0.204 0.208 ...
##   ..$ premax: num [1:33722] 162 159 167 165 165 ...
##   ..$ premin: num [1:33722] 7.37 5.14 4.21 10.98 9.8 ...
##   ..$ presum: num [1:33722] 3935 3799 3975 4134 4089 ...
##   ..$ tmpmax: num [1:33722] 20.9 21.3 21.1 20.4 20.9 ...
##   ..$ tmpmin: num [1:33722] -7.39 -7.21 -7.3 -7.62 -7.33 ...
##   ..$ tmpavg: num [1:33722] 8.17 8.5 8.49 7.78 8.2 ...
##   ..$ pop   : num [1:33722] 18.69 0.91 8.94 10.31 6.46 ...
##   ..$ ntl   : num [1:33722] 6.15 4.32 2.14 4.64 6.79 ...
##   ..$ lulc  : num [1:33722] 10 10 10 10 10 10 10 10 10 10 ...
##   ..$ elev  : num [1:33722] 1720 1638 1662 1835 1717 ...
##   ..$ slope : num [1:33722] 3.89 2.63 3.06 3.37 3.93 ...
##   ..$ aspect: num [1:33722] 114 158 136 102 120 ...
##  $ : tibble [14,840 × 13] (S3: tbl_df/tbl/data.frame)
##   ..$ fvc   : num [1:14840] 0.205 0.197 0.175 0.163 0.21 ...
##   ..$ premax: num [1:14840] 165 161 160 162 165 ...
##   ..$ premin: num [1:14840] 10.35 6.02 4.97 4.48 11.04 ...
##   ..$ presum: num [1:14840] 4098 3908 3848 3910 4131 ...
##   ..$ tmpmax: num [1:14840] 20.7 21.2 21.5 21.3 20.4 ...
##   ..$ tmpmin: num [1:14840] -7.46 -7.19 -7.07 -7.24 -7.52 ...
##   ..$ tmpavg: num [1:14840] 8.04 8.44 8.67 8.52 7.84 ...
##   ..$ pop   : num [1:14840] 12.33 13.78 4.59 6.52 2.97 ...
##   ..$ ntl   : num [1:14840] 6.37 7.79 7.23 9.96 4.59 ...
##   ..$ lulc  : num [1:14840] 10 10 10 10 10 10 10 10 10 10 ...
##   ..$ elev  : num [1:14840] 1762 1654 1600 1649 1805 ...
##   ..$ slope : num [1:14840] 3.41 3.19 2.61 3.06 3.76 ...
##   ..$ aspect: num [1:14840] 126 130 179 201 146 ...
##  $ : tibble [8,268 × 13] (S3: tbl_df/tbl/data.frame)
##   ..$ fvc   : num [1:8268] 0.2 0.174 0.175 0.182 0.193 ...
##   ..$ premax: num [1:8268] 164 161 165 171 164 ...
##   ..$ premin: num [1:8268] 7.57 5.37 6.07 5.78 9.32 ...
##   ..$ presum: num [1:8268] 4022 3896 4026 4177 4072 ...
##   ..$ tmpmax: num [1:8268] 21 21.6 21.3 20.8 21 ...
##   ..$ tmpmin: num [1:8268] -7.22 -6.96 -7.15 -7.46 -7.15 ...
##   ..$ tmpavg: num [1:8268] 8.36 8.81 8.57 8.19 8.35 ...
##   ..$ pop   : num [1:8268] 5.82 15.87 20.4 8.66 1.55 ...
##   ..$ ntl   : num [1:8268] 8.33 8.39 13.18 2.69 11 ...
##   ..$ lulc  : num [1:8268] 10 10 10 10 10 10 10 10 10 10 ...
##   ..$ elev  : num [1:8268] 1684 1567 1642 1690 1693 ...
##   ..$ slope : num [1:8268] 3.43 2.13 3.48 3.22 3.56 ...
##   ..$ aspect: num [1:8268] 115 159 224 207 133 ...
##  $ : tibble [5,240 × 13] (S3: tbl_df/tbl/data.frame)
##   ..$ fvc   : num [1:5240] 0.188 0.162 0.168 0.186 0.189 ...
##   ..$ premax: num [1:5240] 163 162 168 174 164 ...
##   ..$ premin: num [1:5240] 6.86 5.23 4.15 5.99 7.86 ...
##   ..$ presum: num [1:5240] 3992 3922 4040 4254 4047 ...
##   ..$ tmpmax: num [1:5240] 21.2 21.7 21.2 20.8 21.2 ...
##   ..$ tmpmin: num [1:5240] -7.09 -6.9 -7.22 -7.42 -7 ...
##   ..$ tmpavg: num [1:5240] 8.54 8.92 8.53 8.21 8.58 ...
##   ..$ pop   : num [1:5240] 5.64 23.14 9.73 6.84 2.36 ...
##   ..$ ntl   : num [1:5240] 9.1 10.45 5.58 2.89 12.3 ...
##   ..$ lulc  : num [1:5240] 10 10 10 10 10 10 10 10 10 10 ...
##   ..$ elev  : num [1:5240] 1645 1539 1611 1677 1643 ...
##   ..$ slope : num [1:5240] 2.96 1.86 3.19 3.32 2.79 ...
##   ..$ aspect: num [1:5240] 122 174 192 213 132 ...
##  $ : tibble [3,607 × 13] (S3: tbl_df/tbl/data.frame)
##   ..$ fvc   : num [1:3607] 0.196 0.169 0.165 0.174 0.188 ...
##   ..$ premax: num [1:3607] 165 161 165 168 175 ...
##   ..$ premin: num [1:3607] 9.19 5.07 5.89 4.14 5.51 ...
##   ..$ presum: num [1:3607] 4081 3885 4035 4064 4281 ...
##   ..$ tmpmax: num [1:3607] 20.9 21.7 21.6 21.3 20.7 ...
##   ..$ tmpmin: num [1:3607] -7.2 -6.91 -6.99 -7.17 -7.39 ...
##   ..$ tmpavg: num [1:3607] 8.3 8.86 8.79 8.63 8.23 ...
##   ..$ pop   : num [1:3607] 2.69 11.89 27.15 12.59 4.31 ...
##   ..$ ntl   : num [1:3607] 8.82 9.36 12.72 6.77 2.09 ...
##   ..$ lulc  : num [1:3607] 10 10 10 10 10 10 10 10 10 10 ...
##   ..$ elev  : num [1:3607] 1705 1557 1577 1585 1680 ...
##   ..$ slope : num [1:3607] 3.37 1.92 2.69 2.89 3.33 ...
##   ..$ aspect: num [1:3607] 141 130 200 201 218 ...
##  $ : tibble [2,634 × 13] (S3: tbl_df/tbl/data.frame)
##   ..$ fvc   : num [1:2634] 0.172 0.159 0.177 0.179 0.208 ...
##   ..$ premax: num [1:2634] 161 163 166 170 165 ...
##   ..$ premin: num [1:2634] 5.53 4.98 4.19 4.17 8.3 ...
##   ..$ presum: num [1:2634] 3924 3969 4003 4115 4133 ...
##   ..$ tmpmax: num [1:2634] 21.6 21.8 21.2 21.3 21.1 ...
##   ..$ tmpmin: num [1:2634] -6.91 -6.9 -7.15 -7.17 -7.06 ...
##   ..$ tmpavg: num [1:2634] 8.84 8.96 8.52 8.64 8.46 ...
##   ..$ pop   : num [1:2634] 4.79 23.35 33.75 6.38 8.72 ...
##   ..$ ntl   : num [1:2634] 9.65 11.31 11.85 7.32 5.76 ...
##   ..$ lulc  : num [1:2634] 10 10 10 10 10 10 10 10 10 10 ...
##   ..$ elev  : num [1:2634] 1568 1528 1632 1585 1670 ...
##   ..$ slope : num [1:2634] 1.92 2.08 3.24 2.86 2.62 ...
##   ..$ aspect: num [1:2634] 129 181 169 222 164 ...
##  $ : tibble [2,002 × 13] (S3: tbl_df/tbl/data.frame)
##   ..$ fvc   : num [1:2002] 0.169 0.162 0.176 0.184 0.203 ...
##   ..$ premax: num [1:2002] 162 163 167 172 165 ...
##   ..$ premin: num [1:2002] 5.12 4.14 3.82 4.05 8.5 ...
##   ..$ presum: num [1:2002] 3957 3949 4022 4170 4145 ...
##   ..$ tmpmax: num [1:2002] 21.6 21.7 21.4 21.2 20.9 ...
##   ..$ tmpmin: num [1:2002] -6.92 -6.9 -7.08 -7.21 -7.09 ...
##   ..$ tmpavg: num [1:2002] 8.86 8.96 8.68 8.57 8.32 ...
##   ..$ pop   : num [1:2002] 5.81 15.22 26.95 6.31 11.41 ...
##   ..$ ntl   : num [1:2002] 9.42 11.06 12.13 6.12 3.6 ...
##   ..$ lulc  : num [1:2002] 10 10 10 10 10 10 10 10 10 10 ...
##   ..$ elev  : num [1:2002] 1560 1524 1584 1605 1708 ...
##   ..$ slope : num [1:2002] 2.16 2.37 2.9 3.02 2.77 ...
##   ..$ aspect: num [1:2002] 129 194 172 211 161 ...
##  $ : tibble [1,561 × 13] (S3: tbl_df/tbl/data.frame)
##   ..$ fvc   : num [1:1561] 0.175 0.169 0.179 0.196 0.198 ...
##   ..$ premax: num [1:1561] 163 163 168 164 166 ...
##   ..$ premin: num [1:1561] 5.42 3.72 3.66 8.46 5.79 ...
##   ..$ presum: num [1:1561] 4014 3950 4050 4134 4138 ...
##   ..$ tmpmax: num [1:1561] 21.5 21.7 21.4 20.8 21.4 ...
##   ..$ tmpmin: num [1:1561] -6.97 -6.91 -7.05 -7.08 -6.82 ...
##   ..$ tmpavg: num [1:1561] 8.78 8.95 8.76 8.28 8.75 ...
##   ..$ pop   : num [1:1561] 4.5 17.55 15.72 12.63 7.42 ...
##   ..$ ntl   : num [1:1561] 8.93 10.71 10.93 2.97 3.35 ...
##   ..$ lulc  : num [1:1561] 10 10 10 10 10 10 10 10 10 10 ...
##   ..$ elev  : num [1:1561] 1581 1524 1563 1723 1599 ...
##   ..$ slope : num [1:1561] 2.3 2.44 2.91 2.8 2.83 ...
##   ..$ aspect: num [1:1561] 137 191 175 177 150 ...
##  $ : tibble [1,253 × 13] (S3: tbl_df/tbl/data.frame)
##   ..$ fvc   : num [1:1253] 0.177 0.177 0.178 0.186 0.19 ...
##   ..$ premax: num [1:1253] 164 164 169 160 162 ...
##   ..$ premin: num [1:1253] 5.22 3.56 3.34 10.52 7.39 ...
##   ..$ presum: num [1:1253] 4046 3990 4098 4058 4069 ...
##   ..$ tmpmax: num [1:1253] 21.5 21.6 21.5 21.2 21 ...
##   ..$ tmpmin: num [1:1253] -6.98 -6.96 -7.07 -6.86 -6.98 ...
##   ..$ tmpavg: num [1:1253] 8.77 8.86 8.79 8.68 8.4 ...
##   ..$ pop   : num [1:1253] 6.1 18.38 8.83 9.96 8.91 ...
##   ..$ ntl   : num [1:1253] 7.901 11.324 9.294 0.611 2.963 ...
##   ..$ lulc  : num [1:1253] 10 10 10 10 10 10 10 10 10 10 ...
##   ..$ elev  : num [1:1253] 1581 1547 1552 1640 1692 ...
##   ..$ slope : num [1:1253] 2.41 2.83 3.04 2.29 2.85 ...
##   ..$ aspect: num [1:1253] 130 182 211 194 163 ...

discvar = names(select(fvc1000,-c(fvc,lulc)))
g1 = sesu_opgd(fvc ~ ., data = fvc,su = su,discvar = discvar,cores = 6)
## Error in `purrr::map2()` at gdverse/R/sesu_opgd.R:60:3:
## ℹ In index: 1.
## Caused by error in `checkForRemoteErrors()`:
## ! 6 nodes produced errors; first error: Can't extract columns past the end.
## ℹ Location 2 doesn't exist.
## ℹ There is only 1 column.
g1
## Error: object 'g1' not found
plot(g1)
## Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'plot': object 'g1' not found