Hi,

> If I know the co-ordinates for all 4 corners of a non-north-up image, how
> can I calculate the 6 GDAL affine transform parameters?
>

There is a script gcps2wld.py in GDAL:

You should assign the corners (ground control points) to a raster by:

gdal_translate -of vrt -gcp 0 0 19.3 50.5 ... input.tif temp.vrt

and then

gcps2wld.py temp.vrt

will return a world file.

Best regards

Klokan

P.S.
Direct calculation in pure python is possible for example with:
http://code.google.com/p/oldmapsonline/source/browse/trunk/gcps2wld/python-no-dep/gcps2wld.py

with documetation and a Java version at:
http://blog.oldmapsonline.org/2008/06/georeferencing-images-by-control-points.html
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to