

Exploring vector & raster ORNL -- on Oak Ridge DAAC but not in Cloud Data

- [Tucker et al, _Nature_ 2023](https://doi.org/10.1038/s41586-022-05653-6)
- [Visualization tool](https://trees.pgc.umn.edu/app)
- [Data on ORNL DAAC](https://doi.org/10.3334/ORNLDAAC/2117)
- [github repo](https://github.com/ankitkariryaa/An-unexpectedly-large-count-of-trees-in-the-western-Sahara-and-Sahel/)

```{r}
library(earthdatalogin)
edl_set_token()
```


```{r}

```

3 terrabyte GPKG mosaic

```{r}
tree <- ("/vsicurl/https://daac.ornl.gov/daacdata/global_vegetation/Semi-Arid_Tree_Carbon_50cm/data/utm_28/SSAr2_32628_GE01-QB02-WV02-WV03_PV_001_027_mosaic.gpkg")
ex <- terra::vect(tree)
```



```{r}
tree_density <- "/vsicurl/https://daac.ornl.gov/daacdata/global_vegetation/Semi-Arid_Tree_Carbon_50cm/data/wgs_84_mosaic_100m_tree_carbon_density.tif"

r <- terra::rast(tree_density)

```

