---
title: "data_analysis"
author: "Xiaohao Yang"
vignette: >
  %\VignetteIndexEntry{data analysis}
  %\VignetteEngine{knitr::rmarkdown}
  %\VignetteEncoding{UTF-8}
---

```{r eval=FALSE}
bbox <- c(-83.087174,42.333373,-83.042542,42.358748)
ndvi <- greenSD::get_esa_wc(bbox = bbox, datatype = 'ndvi', year = 2021)
seg <- greenSD::ndvi_to_sem(ndvi$NDVI_p50, threshold = c(0.2, 0.6))
```

```{r eval=FALSE}
sample_data <- terra::rast(system.file("extdata", "detroit_gs.tif", package = "greenSD"))
pwgf <- greenSD::compute_exposure(sample_data, pop_year = 2020, radius = 500, height = F)
```

```{r eval=FALSE}
green <- greenSD::get_tile_green(bbox = c(-83.087174,42.333373,-83.042542,42.358748), 
                                 provider = "esri",
                                 zoom = 16)
r <- terra::ifel(green$green == 0, NA, 1)
m <- greenSD::compute_morphology(r, grid_size = 400)
```