speed up download for "infinite" basemap about contextily HOT 9 OPEN

Yes, but the application I developed uses matplotlib embedded in a wxpython application. So all the maps, rasters, vectors are displayed within a matplotlib canvas. I want to show the infinite basemaps as a background for the layers, so I had to embed the folium map in my matplotlib canvas, witch is impossible, I think. Thank you for your answer, best regards, Michele Zucchelli. Il dom 21 gen 2024, 15:37 Dani Arribas-Bel ***@***.***> ha scritto:

Ok, tomorrow I ' ll do it Il gio 25 gen 2024, 17:16 Martin Fleischmann ***@***.***> ha scritto:

Hi, I've tryed open a pull request but I cannot, 'cause I'm not a collaborator and I don't know how to become a collaborator. So I'm going to share my edits by email. Hi everyone, I have modified contextily on my own in order to speeding up download and reprojection. I would like to share this improvements with you: 1) Use always threads instead of processes to handle multiple connections. Simply change the following line in the function: "bounds2img": "preferred_backend = "threads" The threads are much faster to be initialized and for operations I/O like the tiles dowload python allows the GIL. 2)Create a session with the request package. The session keeps active the connection with the server, speeding up the tiles downlod time. session = requests.Session() # Configura il numero massimo di retry per le connessioni retries = Retry(total=5, backoff_factor=0.1, status_forcelist=[500, 502, 503, 504]) session.mount('http://', HTTPAdapter(max_retries=retries)) session.mount('https://', HTTPAdapter(max_retries=retries)) timeout = 10 #10 secondi def _retryer(tile_url, _, __): """ Attempt to get a tile from a URL. Arguments --------- tile_url : str Properly-formatted URL for a tile provider. Returns ------- request object containing the web response. """ try: response = session.get(tile_url, headers={"user-agent": USER_AGENT}, timeout = timeout) response.raise_for_status() # Solleva un'eccezione per codici di risposta 4xx o 5xx except requests.HTTPError as http_err: if response.status_code == 404: raise requests.HTTPError(f"Tile URL resulted in a 404 error. Check your tile URL:\n{tile_url}") from http_err else: raise # Rilancia l'eccezione per altri errori HTTP return response 3)I have implemented a reproject function that is faster, because doesn't create a temporary file but it is "on the fly": def warp_tiles(img, extent, t_crs, resampling=Resampling.bilinear): """ Riproietta un array NumPy da un CRS a un altro. Args: - img: Array NumPy dell'immagine. - t_crs: Sistema di coordinate destinazione. Returns: - img_reprojected: Array NumPy dell'immagine riproiettata. """ src_crs = 3857 # Ottieni le dimensioni dell'immagine e la trasformazione di partenza e di arrivo height, width, bands = img.shape # --- https://rasterio.readthedocs.io/en/latest/quickstart.html#opening-a-dataset-in-writing-mode minX, maxX, minY, maxY = extent x = np.linspace(minX, maxX, width) y = np.linspace(minY, maxY, height) resX = (x[-1] - x[0]) / width resY = (y[-1] - y[0]) / height transform = from_origin(x[0] - resX / 2, y[-1] + resY / 2, resX, resY)#trasformazione di partenza # riproietta le coordinate w, e, s, n = extent w, s = GIS.Convert_coordinates(w, s, src_crs, t_crs) e, n = GIS.Convert_coordinates(e, n, src_crs, t_crs) extent_reproj = (w, e, s, n) minX, maxX, minY, maxY = extent_reproj x = np.linspace(minX, maxX, width) y = np.linspace(minY, maxY, height) resX = (x[-1] - x[0]) / width resY = (y[-1] - y[0]) / height transform_t = from_origin(x[0] - resX / 2, y[-1] + resY / 2, resX, resY) #trasformazione di arrivo # Prepara l'array di output img = img.transpose(2, 0, 1) # Riproietta l'immagine img_reprojected = np.empty(shape=(bands, height, width), dtype=img.dtype) for i in range(bands): # Itera sui canali dell'immagine (per esempio, RGB) reproject( source=img[i, :, :], destination=img_reprojected[i, :, :], src_transform=transform, src_crs=src_crs, dst_transform=transform_t, dst_crs=t_crs, resampling=resampling) #restituisci l'immagine e l'estensione riproiettati return img_reprojected.transpose(1, 2, 0), extent_reproj I hope you enjoy this functions, let me know if you have any doubt or wheter I committed any error. Best regards, Michele Zucchelli. Il giorno gio 25 gen 2024 alle ore 17:54 Michele Zucchelli < ***@***.***> ha scritto: