def indices_to_coords(x, y, geomatrix) -> tuple:
	"""Returns geographic coordinates of pixel

	Calculates coordinates at x,y based on 
	geomatrix

	***

	Parameters
	----------
	x:int
		Column position of pixel within image
	y:int
		Row position of pixel within image
	geomatrix:geomatrix
		As returned by gdal.GetGeoTransform()

	Returns
	-------
	Tuple of geographic coordinates in
	form (x, y). Note that output will be in
	the projected coordinates of input geomatrix
	"""

	ulX = geomatrix[0]
	ulY = geomatrix[3]
	xDist = geomatrix[1]
	yDist = geomatrix[5]
	rtnX = geomatrix[2]
	rtnY = geomatrix[4]

	geoX = ulX + (xDist * x)
	geoY = ulY + (yDist * y)

	return (geoX, geoY)