def _latlon_to_tile(lat, lon, zoom): n = 2 ** zoom x = n * (lon + 180) / 360. y = n * (1 - (np.arcsinh(np.tan(np.deg2rad(lat))) / np.pi)) / 2. return int(x), int(y)