Getting geometry of raster file

[maxott]

Most likely a newbie question but how can I get the geometry or boundary of a geotiff file?

[visr]

What would you like the result to look like. Is this sufficient?

using ArchGDAL
const AG = ArchGDAL

p = raw"c:\Users\visser_mn\.julia\dev\ArchGDAL\test\data\utmsmall.tif"
ds = AG.read(p)

geotransform = AG.getgeotransform(ds)
xmin = geotransform[1]
ymax = geotransform[4]
xmax = xmin + geotransform[2] * AG.width(ds)
ymin = ymax + geotransform[6] * AG.height(ds)

Although I think in GeoArrays.jl it is done a bit more properly, by converting this to an AffineMap from CoordinateTransformations.jl. I think that should handle rotation as well, if needed.

Perhaps something similar would be good to add here as well, as I imagine this is a common use case.