[gdal-dev] Longitude and latitude raster bands

[Knut-Frode Dagestad]

Hi,

For datasets with only GCPs for geolocation, is there any convenient and 
efficient method to extract latitude and longitude values as arrays of 
same size as the raster bands, using the Python API?

And similarly: is there an efficient method to return full arrays of the 
azimuth direction, e.g. defined as 0 for pixels where north is "up" (as 
it would always be in a lon-lat or Mercator projection).
This angle is relevant e.g. when one want to know the sensor look 
direction for satellite images.

Best regards from Knut-Frode