I have come to a place in some personal-project code where I need to use skimage.transform.warp and I keep tabbing over to the documentation, sighing, and tabbing away.

I want to map from a small input to a large output with interpolation in the output and I’m going to be able to figure it out but it’s going to feel like doing my taxes the whole time.

@vruba Depending on the nature of the transform, this may be possible with gdalwarp. I'm pretty well versed in its esoteric features, so consider me nerd sniped if you want to send something along.

gdalwarp directly supports an affine transform, and I think PROJ supports a homography. A geolocation array opens up lots of possibilities depending on what can be placed in DST_SRS.

@geowurster So I sort of talked myself out of GDAL at the beginning of thinking about this, but maybe it’s been the obvious choice all along.

I’m just trying to reproject PACE OCI L1B data, which does the weather satellite thing of having bands for longitude and latitude values, so it’s essentially one GCP per pixel.

@geowurster The thing is that it’s so oblique at the edges that the terrain correction baked into the GCPs (or whatever the better name is) matters quite a bit. So the common sense thing would be to discard most of them, but I’d rather not do that. Therefore it’s order 1e6 GCPs, which feels like a lot to ask even GDAL, but I haven’t tested!