Basic function to create a geotransform as used by GDAL.
geo_transform0(px, ul, sh = c(0, 0))
pixel resolution (XY, Y-negative)
grid offset, top-left corner
affine shear (XY)
vector of parameters xmin, xres, yskew, ymax, xskew, yres
geo_world0()
which uses the same parameters in a different order