From bca00e7f79bdc33cd48de29278e60683146a291d Mon Sep 17 00:00:00 2001 From: Anshul Singhvi Date: Tue, 22 Oct 2024 17:42:34 -0700 Subject: [PATCH] Update zonal.jl --- src/methods/zonal.jl | 10 +++++++++- 1 file changed, 9 insertions(+), 1 deletion(-) diff --git a/src/methods/zonal.jl b/src/methods/zonal.jl index 47638668..8da0d679 100644 --- a/src/methods/zonal.jl +++ b/src/methods/zonal.jl @@ -4,9 +4,17 @@ Calculate zonal statistics for the the zone of a `Raster` or `RasterStack` covered by the `of` object/s. +Since `f` is the first argument, you can use Julia's `do`-block syntax with `zonal`, +like so: +```julia +zonal(raster; of = geom) do r # r is an iterator if skipmissing=true, or a raster if skipmissing=false + return r # or do some operation +end +``` + # Arguments -- `f`: any function that reduces an iterable to a single value, such as `sum` or `Statistics.mean` +- `f`: any function that reduces an iterable to a single value, such as `sum` or `Statistics.mean` - but it can even be `identity`. - `x`: A `Raster` or `RasterStack` - `of`: A `DimTuple`, `Extent`, $OBJ_ARGUMENT