[gdal-dev] Longitude and latitude raster bands

[Frank Warmerdam]

On Tue, Jan 10, 2012 at 5:36 AM, Knut-Frode Dagestad
<knutfrodesoppel at hotmail.com> wrote:
> 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.

Knut-Fronde,

I am not aware of a preexisting solution for either operation.

Best regards,
-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Software Developer