-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
Hi Mitch,
At the moment, mpop can only write geotiffs with a given area (or
geotransform). We actually never needed to write unprojected data to
geotiff yet...
However, you can have a look in mpop/imageo/geo_image.py, there you
have a GeoImage method called geotiff_save. It's using gdal's python
api to write data to geotiffs. The gdal api is a bit obscure, but I
found that instead of using setGeoTransform, one could use setGCPs
(ground controll points).
I tried it a bit, but couldn't find any info on how to set up GCPs.
Here what I tried (after line 309):
logger.exception("Could not load geographic data, invalid area")
srs = osr.SpatialReference()
srs.ImportFromEPSG(4326)
lons, lats = area.get_lonlats()
dst_ds.SetGCPs((lats, lons), srs.ExportToWkt())
But I got some strange error. Feel free to try it out, or tell us if
you find more info on the web...
Best regards,
Martin
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v2.0.22 (GNU/Linux)
iQEcBAEBAgAGBQJVUaR4AAoJEBdvyODiyJI4n/cH/Axe5he1eZnhC3vPZeec0I+G
6JfgQdsH8LdoFtgvPotSM1gEF/IVuvk1bmRJ0amFrBMEOgcH71id4aPz8MzkMmQD
6ey4T3iVnL3E9aal5cvXzBgWvEYbZnP79I9vXLWEhF8O9v51jQxZVn+dii3jgYy3
7X1VPYEx7kXDREyEHJ9AmIrr6tWVKALo/D4sozjk4MTKW5962IW3DqAqB3gvRUOf
pMs6c4jNQdzh0lrImJJNdj4QUPUGmjQketjXqYPWPQxNHZUTlt57/Bzts/IkBdOm
0VRK3CWtKhaJWsAxYG9UtcivrlILXxm+rbt2Xx+zeY3TeKiGC3RKt5ogxxgD1AI=
=EFDr
-----END PGP SIGNATURE-----