Anshul Singhvi

Zonal improvements without geometry lookup

11

- `spatialslices = false`: if true, run `f` on the slices of the datacube in `X, Y`. You can also pass a tuple of dimensions to specify which dims to...

RasterStack constructor with no rasters cause a stackoverflow

2

I forgot parens for namedtuple here ```julia-repl julia> precip_cellarea_stack = RasterStack(; precip = precip, cellarea = cellarea(precip)) ERROR: StackOverflowError: Stacktrace: [1] RasterStack(layers::@NamedTuple{}; resize::Rasters.NoKW, _layers::@NamedTuple{}, dims::Tuple{}, refdims::Tuple{}, missingval::@NamedTuple{}, metadata::DimensionalData.Dimensions.Lookups.NoMetadata, layermetadata::@NamedTuple{}, layerdims::@NamedTuple{},...

make zonal open stacks of length < 5

8

not sure if we want this size or not but that's easy to configure

[WIP] Add a `bylayer` keyword to `zonal`

5

Successor to #775. This provides a switch that allows you to iterate directly over the values of a rasterstack, instead of iterating over the layers and then the layer values....

load CDM CRS if they are available in string form

8

we still cannot load CF CRS natively, but most writers will also write the CRS in string form, which we can read. The key indicates the format. Fix #935 Partially...

Native resampling

1

It seems that with Interpolations.jl and Proj.jl, we can already get native resampling pretty easily using Images.jl's inverse warp functionality. This is only doable if we have a known desired...

Arbitrary metadata not read from rasters in GDAL

2

See the files in this zipfile: https://github.com/user-attachments/files/19419830/sub_leaf_type_2021.zip ```julia using ArchGDAL dataset = ArchGDAL.read("sub_leaf_type_2021/sub_leaf_type_2021.tif") ArchGDAL.metadata(dataset) ``` yields ```julia 2-element Vector{String}: "AREA_OR_POINT=Area" "DataType=Generic" ``` but if I open it as a Raster...

Use MetopDatasets as prototype for loading matrix-like data

5

It would be very cool to get https://github.com/eumetsat/MetopDatasets.jl as a backend for Rasters. It also looks like all the metopdatasets stuff has matrix lookups, so it would be a good...

Using spatial indices in polygon burning

1

Since GeometryOps will have the natural indexing structure soon, it may be useful to explore whether we should use that in Rasters as well. We can use it to accelerate...

experiment with best_manifold implementation here

2

seems we need this now, will move the code to GeometryOps / Proj / GFT / wherever later but want to at least get a working implementation here. this is...