Rafael Schouten
Rafael Schouten
With tiled batches we can just keep filling the same small ector of points that we pass to Proj to tranform. We can swap over `resample` at the same time....
Actually we just let DiskArrays handle that currently
Ah right we need a version for Irregular and Regular
No regular is by definition gapless, there is no way to represent gaps. Actually only Explicit lookups have gaps and not Irregular, Irregular cant represent gaps either
Its probably already the fastest that exists by an order of magnitude, so I'm happy to call it off wherever. But generally switching from GDAL to Proj for these things...
Good idea, minimise conversions where possible
You can use the `crs` keyword to `write` (we do want to retain the crs in DataAPI metadata as a general GeoInterface pattern, but no one has implemented it yet)
Yeah, the real problem is `testgeometry` doesn't test anywhere near enough things
How do we get rid of ArchGDAL?
But how do we deal with Proj/Wkt, a lot of different text strings can be wgs84