def pixelToLatLng(raster_dataset, col, row):
  ds = raster_dataset
  gt = ds.GetGeoTransform()
  srs = osr.SpatialReference()
  srs.ImportFromWkt(ds.GetProjection())
  srsLatLong = srs.CloneGeogCS()
  ct = osr.CoordinateTransformation(srs,srsLatLong)
  ulon = col*gt[1]+gt[0]
  ulat = row*gt[5]+gt[3]
  # Transform the point into the GeoTransform space
  (lon,lat,holder) = ct.TransformPoint(ulon,ulat)
  return (lat, lon)
from http://zacharybears.com/using-python-to-translate-latlon-locations-to-pixels-on-a-geotiff/
dataset/ETHPy150Open trailbehind/DeepOSM/src/geo_util.py/pixelToLatLng