Transform from pixel-space to world-space.

get_coordinates.py

#!/usr/bin/env python

import argparse
import rasterio as rio
import rasterio.warp


def get_coordinates(srcpath, pixel):

    with rio.open(srcpath) as src:
        print("Coordinate reference system %s" % src.crs['init'])
        src_crs = src.crs

        # Transform from pixel coordinates to world coordinates in the CRS of the image
        xy = src.ul(*pixel)

    # To transform to another coordinate reference system
    # Note the need to use the image's CRS as an intermediary. There might
    # be a way to skip this step.
    dst_crs = {'init': 'epsg:4326'}
    print rasterio.warp.transform(src_crs, dst_crs, *zip(xy))


if __name__ == '__main__':

    parser = argparse.ArgumentParser()
    parser.add_argument("srcpath")
    parser.add_argument("--pixel", nargs=2, type=int)
    args = parser.parse_args()

    get_coordinates(args.srcpath, args.pixel)