docs: def tile_coord(self, xtile, ytile, zoom):
          """
          This returns the NW-corner of the square. Use the function
          with xtile+1 and/or ytile+1 to get the other corners.
          With xtile+0.5 & ytile+0.5 it will return the center of the tile.
          http://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Tile_numbers_to_lon..2Flat._2
          """
          assert self.tile_srid == 3857, 'Custom tile projection not supported yet'
          n = 2.0 ** zoom
          lon_deg = xtile / n * 360.0 - 180.0
          lat_rad = math.atan(math.sinh(math.pi * (1 - 2 * ytile / n)))
          lat_deg = math.degrees(lat_rad)
          return (lon_deg, lat_deg)