Some information: I must not re-sample the image, the pixels are rotated, and the image opens in Arcmap just fine.
now, there are 15 GCPs on the raster (which make the pixels rotate)
1.When I first try to use
gdallocationinfo I get:
ERROR 6: No translation for an empty SRS to PROJ.4 format is known.
gdal_edit -ro -a_srs epsg:4326 "Image", this creates a .aux.xml file with the same name as my image, gdal now reads them both together, image still works on Arcmap.
then I use
gdallocationinfo again and get:
ERROR 1: Cannot get geotransform
gcps2wld "Image" to get the geotransformation parameters, and insert it to the .aux.xml file via python GDAL
now Arcmap can't read the image.
and when I use
gdalinfo, I get the following:
Corner Coordinates: Upper Left ( 0.0000360, -0.0000277) ( 0d 0' 0.13"E, 0d 0' 0.10"S)
Lower Left ( -0.145, 887589.718) ( 0d 8'43.59"W,Invalid angle)
Upper Right ( -0.068,-1152505.040) ( 0d 4' 3.91"W,Invalid angle)
Lower Right ( -0.213, -264915.321) ( 0d12'47.63"W,Invalid angle)
Center ( -0.107, -132457.661) ( 0d 6'23.75"W,Invalid angle)
and if I Ignore that and use
gdallocationinfo, coordinates that are supposed to be in the image give me errors for being outside the image borders.
I tried using
gdal_edit -a_ullr to set the corner coordinates, but that just norths up the image (no rotation).
how do I set up the image right (without resampling) so I can get the right pixels in the right coordinates?