code: def pixelToMap(pX, pY, geoTransform):
      
      pX = np.asarray(pX, dtype=float)
      pY = np.asarray(pY, dtype=float)
      pX += 0.5
      pY += 0.5
      mX, mY = applyGeoTransform(pX, pY, geoTransform)
      return mX, mY
docs: Convert pixel coordinates to map coordinates based on geotransform
      
      Accepts float or NumPy arrays
  
      GDAL model used here - upper left corner of upper left pixel for mX, mY (and in GeoTransform)