reproject doens't reproject AutoDim dimensions

[tiemvanderdeure]

The dim field of a lookup doesn't inherit the supertype of the dimension (XDim, etc), but defaults to Rasters.AutoDim(). reproject dispatches on this field, So reproject will silently not do anything and return the same values - but with a different crs associated!

using Rasters
import Rasters: XDim
xdim = X(Projected(1:10; crs = EPSG(4326)))
xdim isa XDim # true
lookup(xdim).dim isa XDim # false
lookup(reproject(xdim; crs = EPSG(3857))).data == lookup(xdim).data # true!

I get that when reprojecting we want all dimensions that are not x or y to keep the same values, and this is not a problem after constructing a Raster using the dimension - but maybe there's a way to change the dispatch on the internal reproject methods so this doesn't happen?

[rafaqz]

Arrg yes thats not cool. We should error on AutoDim. But if you format it first it should work?

[rafaqz]

Also nothing that is not X or Y should get a crs...

[tiemvanderdeure]

But if you format it first it should work?

Yes it does

Also nothing that is not X or Y should get a crs...

Looked at it again and this only affects Projected dimensions, and I don't think anything that is not X or Y is usually projected