[WIP] add Rasters.jl benchmarks

.gitattributes

@@ -0,0 +1,2 @@
+# GitHub syntax highlighting
+pixi.lock linguist-language=YAML linguist-generated=true

.gitignore

@@ -1,3 +1,14 @@
+# R files
 .Rproj.user
 .Rhistory
 .Rdata
+# pixi environments
+.pixi
+*.egg-info
+# Julia environment
+rasters_jl/Manifest.toml
+# data
+data/LC08_L1TP_190024_20200418_20200822_02_T1.zip
+data/LC08_L1TP_190024_20200418_20200822_02_T1/*
+data/Mammals_Terrestrial/*
+data/stack.nc

README.md

@@ -21,6 +21,9 @@ You may also be interested in the [vector processing benchmarks](https://github.
 - [raster](https://github.com/rspatial/raster)
 - [exactextractr](https://github.com/isciences/exactextractr)
 
+**Juliia**:
+- [Rasters.jl](https://github.com/rafaqz/Rasters.jl)
+
 ## Reproduction
 1. Download raster data (851 MB) from [Google Drive](https://drive.google.com/uc?id=1lzglfQJqlQh9OWT_-czc5L0hQ1AhoR8M&export=download) or [Earth Explorer](https://earthexplorer.usgs.gov/) (original source, registration required) and then unzip to `data/`.
 2. Run all benchmarks using batch script (`run_benchmarks.sh`) or single benchmarks files.
@@ -40,6 +43,10 @@ Rscript stars/crop.R
 python3 rasterio/crop.py
 ```
 
+```
+julia rasters_jl/crop.jl
+```
+
 ## Dataset
 Landsat 8 satellite scene (10 bands, 30 m resolution, 7771 x 7871 pixels) was used for tests.
 

pixi.toml

@@ -0,0 +1,74 @@
+[project]
+name = "raster-benchmark"
+authors = ["Anshul Singhvi <[email protected]>"]
+channels = ["conda-forge"]
+description = "Benchmarking raster processing libraries"
+platforms = ["osx-64", "linux-64"]
+version = "0.1.0"
+
+[tasks.benchmarks]
+
+cmd = "echo 'Ran benchmarks'"
+depends-on = ["r-raster.benchmark", "r-stars.benchmark", "r-terra.benchmark", "r-exactextractr.benchmark", "py-rasterio.benchmark", "py-rasterstats.benchmark", "py-rioxarray.benchmark", "julia-rastersjl.benchmark"]
+
+
+[environments]
+r-raster = ["r", "r-sf", "r-raster"]
+r-stars = ["r", "r-sf", "r-stars"]
+r-terra = ["r", "r-sf", "r-terra"]
+r-exactextractr = ["r", "r-exactextractr", "r-terra"]
+
+py-rasterio = ["python", "py-rasterio"]
+py-rasterstats = ["python", "py-rasterstats"]
+py-rioxarray = ["python", "py-rioxarray"]
+
+julia-rasters = ["julia-rasters"]
+
+
+[feature.r.dependencies]
+r = ">=4.3,<5"
+r-tibble = ">=3.2.1,<4"
+
+[feature.r-sf.dependencies]
+r = ">=4.3.3,<5"
+r-sf = ">=1.0.15,<2"
+r-s2 = ">=1.1.4,<2"
+
+[feature.python.dependencies]
+python = ">=3.11.6,<4"
+xarray = ">=2024.9.0,<2026"
+fiona = ">=1.9.4,<2"
+pandas = ">=2.2.3,<3"
+geopandas = ">=1,<2"
+
+[feature.r-raster]
+dependencies = {r-raster = ">=3.6_26,<4"}
+# tasks = {benchmark = "cd raster; for path in \"${i}\"/*; do; echo \"$path\"; Rscript \"$path\"; done"}
+
+[feature.r-stars]
+dependencies = {r-stars = ">=0.6_6,<0.7", r-sf = ">=1.0_15,<2", r-s2 = ">=1.1.4,<2"}
+# tasks = {benchmark = "cd stars; for path in \"${i}\"/*; do; echo \"$path\"; Rscript \"$path\"; done"}
+
+[feature.r-terra]
+dependencies = {r-terra = ">=1.8.5,<2"}
+# tasks = {benchmark = "cd terra; for path in \"${i}\"/*; do; echo \"$path\"; Rscript \"$path\"; done"}
+
+[feature.r-exactextractr]
+dependencies = {r-exactextractr = ">=0.10,<0.11"}
+# tasks = {benchmark = "cd exactextractr; for path in \"${i}\"/*; do; echo \"$path\"; Rscript \"$path\"; done"}
+
+[feature.py-rasterio]
+dependencies = {rasterio = ">=1.3.7,<2"}
+# tasks = {benchmark = "cd rasterio; for path in \"${i}\"/*; do; echo \"$path\"; python3 \"$path\"; done"}
+
+[feature.py-rasterstats]
+dependencies = {rasterstats = ">=0.20.0,<0.21"}
+# tasks = {benchmark = "cd rasterstats; for path in \"${i}\"/*; do; echo \"$path\"; python3 \"$path\"; done"}
+
+[feature.py-rioxarray]
+dependencies = {rioxarray = ">=0.18.0,<0.19"}
+# tasks = {benchmark = "cd rioxarray; for path in \"${i}\"/*; do; echo \"$path\"; python3 \"$path\"; done"}
+
+[feature.julia-rasters]
+dependencies = {julia = ">=1.9.3,<2"}
+# tasks = {benchmark = "cd rasters_jl; for path in \"${i}\"/*; do; echo \"$path\"; BENCHMARKING=true julia \"$path\"; done"}