Add tests for unitful `cellarea`
Does what it says on the tin. Unfortunately it looks like Unitful can't convert to Float64 (see https://github.com/PainterQubits/Unitful.jl/issues/742) so going to Proj is where this fails.
It also occurs to me that we could try to promote and rebuild e.g. bignumbers or Unitful quantities when going into reproject. That would allow unitful axes to pass through reproject unchanged. But that's a different PR.
Hmm with units, probably we need to check they are the same units as the crs and otherwise convert. You can Unitful.ustrip to a specific unit, like ustrip(u"m", x)
(We will need a Unitful.jl extension)
I think if the crs units are "US survey foot" and the lookup has meter units we can just error. I guess a dumb first round check is is it in "meters" and then do the ustrip above, or just error. Maybe "kilometers" too. I don't know where to find the spec for the allowed units...
That sounds like a way bigger effort :D, but doable at some point via Proj. It might be nice for reproject, actually, since we could discover the CRS's units (maybe stash that in lookup metadata) and this would all optimize away at compile-time. Then indexing with arbitrary units could just work TM via Unitful's automatic conversion.
Most of the unit finding work would happen in Proj.