Anshul Singhvi
Anshul Singhvi
I specifically want to iterate over namedtuples here.
Let's just make it a kwarg, I can see situations where both options (none missing vs all missing vs some specific layers missing) are necessary.
Huh. I guess that's a GeometryOps issue then, but I'm surprised all the numbers came out NaN. Would it make sense for indexing by a NaN extent to return a...
Yeah I put that there for empty geometries since they wouldn't plot otherwise (and I think reconstruct might also have failed at some point?)
> Question is does that error, or does it return a zero length array. I would say it should return a zero length array? If it was a type issue,...
Also, holy shit, polygonize just works :O
`zonal` works but `polygonize` does not roundtrip through it, not sure what's going on there. ```julia julia> using StatsBase; Rasters.zonal(countmap, grain; of = GO.polygonize(grain)) 3-element Vector{Union{Missing, Dict{CategoricalValue{String, UInt32}, Int64}}}: Dict{CategoricalValue{String,...
Aaaah I didn't realize DimensionalArrays has a CategoricalArrays extension! This makes much more sense now :D
A basic version of this will go into geocompjl, but I'm not sure where the rest should lie.
I think you could hack around this using categorical arrays for now, you'd basically rasterize the category in as an integer, then rebuild the raster with `rebuild(raster, data = CategoricalArray(parent(raster),...