Anshul Singhvi
Anshul Singhvi
```julia # Prepare a test dataset rs = RasterStack((Raster(rand(X(1:10), Y(1:10)); missingval=NaN) for _ in 1:4)...; missingval = NaN) # Create areas of known zero, just for verification map(rs) do A...
Right now, some in memory operations work with a Raster that wraps a CategoricalArray. Here's an mwe: ```julia using Rasters, CategoricalArrays grain_order = ["clay", "silt", "sand"] grain_char = rand(grain_order, 6,...
It would be cool to have a tutorial that shows all the different things you can do (and have to be careful of) when building a raster from scratch. -...
Currently it's a bit inscrutable: ``` julia> rmin = Rasters.zonal(min, src_srtm; of = zion) ERROR: MethodError: no method matching min(::Base.SkipMissing{Raster{…}}) Closest candidates are: min(::Any, ::Missing) @ Base missing.jl:134 min(::Any, ::Any)...
From geocompjl: ```julia julia> california_raster1 = Rasters.rasterize( last, california_borders; # GI.MultiLineString fill = 1, res = 0.5, shape = :polygon ) ERROR: MethodError: no method matching Rasters.Edges(::GeoInterface.MultiLineStringTrait, ::GeoInterface.Wrappers.MultiLineString{…}, ::Tuple{…}; allocs::Rasters.Allocs{…})...
This would be nice to have! It would probably also be solved if Extent had a GI.RectangleTrait, but that's a more long term thing since it doesn't exist yet. I...
It's super easy to create dimension axes that have vector lookups: ```julia-repl julia> using CairoMakie, DimensionalData julia> points = Point2f.(1, 1:10) 10-element Vector{Point{2, Float32}}: [1.0, 1.0] [1.0, 2.0] [1.0, 3.0]...
``` julia> sstr = read(rs.SST) ERROR: required memory 235358208 is greater than system memory 38846464. Use `lazy=true` if you are loading dataset, and only call `read` on a subset after...
This is on the `cf` branch. I tried to load a raster where `grid_mapping => "crs"`, but there's no `crs` field in the dataset, only a `proj4` field. It also...