On Tue, Jan 10, 2012 at 5:36 AM, Knut-Frode Dagestad <knutfrodesop...@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, warmer...@pobox.com light and sound - activate the windows | http://pobox.com/~warmerdam and watch the world go round - Rush | Geospatial Software Developer _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev