## R 4.3.0
## sf 1.0.12
## stars 0.6.1
## terra 1.7.29
## raster 3.6.21
## exactextractr 0.9.1
## Python 3.8.10
## rasterio 1.3.6
## rasterstats 0.18.0
## rioxarray 0.14.1

Extract values by points

68 275 points

Downsample

30 to 90 m

Calculate NDVI

Load rasters

Write multilayer

10 layers with masks, LZW compression

Crop

by rectangular extent

Zonal statistics

mean values for 1 940 buffers

See also the benchmark from the {exactextractr} package: https://github.com/isciences/exactextractr#performance