EhVDSrSSKJrJrJr SSKrSSKrSSKrSSK r SSK r SSK r SSK r SSK r SSKrSSKrSSKrSSKJrJr SSKJr SSKJrJr SSKJr SSKJr SS KJ r SS K!J"r" S S K#J$r$ SS K%J&r& /SQr'S\RP"5RR-r*\ RV"5r,\"\,SS9r-Sr.Sr/\ R`"\/5 S"Sjr1S#Sjr2Sr3Sr4S\"Rj4Sjr6\"Rj4Sjr7Sr8Sr9S$Sjr:Sr;SrS!r?g)&z1Tools for downloading map tiles from coordinates.)absolute_importdivisionprint_functionN)ImageUnidentifiedImageError)Memory)Paralleldelayed) from_origin) MemoryFile) WarpedVRT) Resampling) providers) TileProvider) bounds2raster bounds2img warp_tileswarp_img_transformhowmany set_cache_dirz contextily-)verbosec.U[Rlg)a Set a cache directory to use in the current python session. By default, contextily caches downloaded tiles per python session, but will afterwards delete the cache directory. By setting it to a custom path, you can avoid this, and re-use the same cache a next time by again setting the cache dir to that directory. Parameters ---------- path : str Path to the cache directory. N)memory store_backendlocation)paths bC:\Users\julio\OneDrive\Documentos\Trabajo\Ideas Frescas\venv\Lib\site-packages\contextily/tile.pyrr,s%)F!c6[R"[SS9 g)NT) ignore_errors)shutilrmtreetmpdirrr _clear_cacher&=s MM&-rc U(d[X5up[X#5up#[UUUUUUSU U S9 upU RupnU unnnn[R"UUU5n[R"UUU5nUSUS- U- nUSUS- U- n[ USUS- - USUS- -UU5n[ R"USSUUU[U RR5SUS 9 n[U5H"nURU S S 2S S 2U4US -5 M$ S S S 5 X4$!,(df  X4$=f) a Take bounding box and zoom, and write tiles into a raster file in the Spherical Mercator CRS (EPSG:3857) Parameters ---------- w : float West edge s : float South edge e : float East edge n : float North edge zoom : int Level of detail path : str Path to raster file to be written source : xyzservices.TileProvider object or str [Optional. Default: OpenStreetMap Humanitarian web tiles] The tile source: web tile provider or path to local file. The web tile provider can be in the form of a :class:`xyzservices.TileProvider` object or a URL. The placeholders for the XYZ in the URL need to be `{x}`, `{y}`, `{z}`, respectively. For local file paths, the file is read with `rasterio` and all bands are loaded into the basemap. IMPORTANT: tiles are assumed to be in the Spherical Mercator projection (EPSG:3857), unless the `crs` keyword is specified. ll : Boolean [Optional. Default: False] If True, `w`, `s`, `e`, `n` are assumed to be lon/lat as opposed to Spherical Mercator. wait : int [Optional. Default: 0] if the tile API is rate-limited, the number of seconds to wait between a failed request and the next try max_retries: int [Optional. Default: 2] total number of rejected requests allowed before contextily will stop trying to fetch more tiles from a rate-limited API. n_connections: int [Optional. Default: 1] Number of connections for downloading tiles in parallel. Be careful not to overload the tile server and to check the tile provider's terms of use before increasing this value. E.g., OpenStreetMap has a max. value of 2 (https://operations.osmfoundation.org/policies/tiles/). If allowed to download in parallel, a recommended value for n_connections is 16, and should never be larger than 64. use_cache: bool [Optional. Default: True] If False, caching of the downloaded tiles will be disabled. This can be useful in resource constrained environments, especially when using n_connections > 1, or when a tile provider's terms of use don't allow caching. Returns ------- img : ndarray Image as a 3D array of RGB values extent : tuple Bounding box [minX, maxX, minY, maxY] of the returned image T)zoomsourcell n_connections use_cacherwGTiffz epsg:3857driverheightwidthcountdtypecrs transformNr) _sm2llrshapenplinspacer rioopenstrr6namerangewrite)r/senrr(r)r*wait max_retriesr+r,ZexthbminXmaxXminYmaxYxyresXresYr8rasterbands rrrDs`N a|a|      # FAggGA! D$d D$"A D$"A bEAaDLA D bEAaDLA DAaD4!8OQrUTAX-=tTJI   !'',,   !HD LL1a:q 1  6M   6Ms 02D-- D=c ^^^U(d[X5up[X#5up#[U5n US:Hn U (a [XX#5nU (aXK- n[XLU S9n[ [ R "XX#U/55nUVs/sH1oRURURURS9PM3 nnU S:d[U [5(d [S5eU S:XdU (dSOSnU (a[R[ 5O[ m[#U US9"UUU4S jU55n[%UU5unnUunnnn[ R&"UU5unn[ R&"UU5unnUUUU4nUU4$s snf) a Take bounding box and zoom and return an image with all the tiles that compose the map and its Spherical Mercator extent. Parameters ---------- w : float West edge s : float South edge e : float East edge n : float North edge zoom : int Level of detail source : xyzservices.TileProvider object or str [Optional. Default: OpenStreetMap Humanitarian web tiles] The tile source: web tile provider or path to local file. The web tile provider can be in the form of a :class:`xyzservices.TileProvider` object or a URL. The placeholders for the XYZ in the URL need to be `{x}`, `{y}`, `{z}`, respectively. For local file paths, the file is read with `rasterio` and all bands are loaded into the basemap. IMPORTANT: tiles are assumed to be in the Spherical Mercator projection (EPSG:3857), unless the `crs` keyword is specified. ll : Boolean [Optional. Default: False] If True, `w`, `s`, `e`, `n` are assumed to be lon/lat as opposed to Spherical Mercator. wait : int [Optional. Default: 0] if the tile API is rate-limited, the number of seconds to wait between a failed request and the next try max_retries: int [Optional. Default: 2] total number of rejected requests allowed before contextily will stop trying to fetch more tiles from a rate-limited API. n_connections: int [Optional. Default: 1] Number of connections for downloading tiles in parallel. Be careful not to overload the tile server and to check the tile provider's terms of use before increasing this value. E.g., OpenStreetMap has a max. value of 2 (https://operations.osmfoundation.org/policies/tiles/). If allowed to download in parallel, a recommended value for n_connections is 16, and should never be larger than 64. use_cache: bool [Optional. Default: True] If False, caching of the downloaded tiles will be disabled. This can be useful in resource constrained environments, especially when using n_connections > 1, or when a tile provider's terms of use don't allow caching. zoom_adjust : int or None [Optional. Default: None] The amount to adjust a chosen zoom level if it is chosen automatically. Values outside of -1 to 1 are not recommended as they can lead to slow execution. Returns ------- img : ndarray Image as a 3D array of RGB values extent : tuple Bounding box [minX, maxX, minY, maxY] of the returned image auto)rW)rPrQzrz/n_connections must be a positive integer value.threads processes)n_jobspreferc3J># UHn[T5"UTT5v M g7fN)r ).0tile_url fetch_tile_fnrGrFs r bounds2img..s&FLU x{;;Is #)r9_process_source_calculate_zoom_validate_zoomlistmttiles build_urlrPrQrX isinstanceint ValueErrorrcache _fetch_tiler _merge_tilesxy)r/rCrDrEr(r)r*rFrGr+r, zoom_adjustprovider auto_zoomritile tile_urlspreferred_backendarraysmergedextentwestsoutheastnorthleftbottomrighttopras `` @rrrssR a|a|v&HIqQ*  $y 9D !tf- .ENSTed##dff$&&#AeITq =# > >JLL $q( 2;FLL- M ]3D EFLUFF"%0NFF%D%u55u%LD&tU#JE3 5&# %F 6>+Us8F cUc[RRnU$[U[5(a [ USSS9nU$[U[ 5(d [S5eSU;a [S5eUnU$)Nurl)r attributionr@zAThe 'url' needs to be a xyzservices.TileProvider object or stringz2The 'url' dict should at least contain a 'url' key) r OpenStreetMapHOTrkr?rdict TypeErrorrm)r)rss rrdrd(s ~**.. O FC FG O % % O   f MNN Orc[XU5nU$r^)_retryer)r`rFrGarrays rroro8s X[ 1E Lrz EPSG:4326cURupEnUupxp[R"XxU5n [R"XU5n U SU S- U- n U SU S- U- n[U SU S- - U SUS- -X5n[ UR SSS5USX#5unnnUR URURUR4nUR SSS5U4$)a Reproject (warp) a Web Mercator basemap into any CRS on-the-fly NOTE: this method works well with contextily's `bounds2img` approach to raster dimensions (h, w, b) Parameters ---------- img : ndarray Image as a 3D array (h, w, b) of RGB values (e.g. as returned from `contextily.bounds2img`) extent : tuple Bounding box [minX, maxX, minY, maxY] of the returned image, expressed in Web Mercator (`EPSG:3857`) t_crs : str/CRS [Optional. Default='EPSG:4326'] Target CRS, expressed in any format permitted by rasterio. Defaults to WGS84 (lon/lat) resampling : [Optional. Default=Resampling.bilinear] Resampling method for executing warping, expressed as a `rasterio.enums.Resampling` method Returns ------- img : ndarray Image as a 3D array (h, w, b) of RGB values (e.g. as returned from `contextily.bounds2img`) ext : tuple Bounding box [minX, maxX, minY, maxY] of the returned (warped) image r-rr.rz EPSG:3857) r:r;r<r _warper transposerrrr)imgrzt_crs resamplingrJr/rKrLrMrNrOrPrQrRrSr8w_imgbounds_s rrr=s@iiGA!#D D"A D"A bEAaDLA D bEAaDLA DAaD4!8OQrUTAX-=tJI aA ;E61[[&,, vzz AF ??1a #V ++rc([XX#U5upVnXW4$)a+ Reproject (warp) an `img` with a given `transform` and `s_crs` into a different `t_crs` NOTE: this method works well with rasterio's `.read()` approach to raster's dimensions (b, h, w) Parameters ---------- img : ndarray Image as a 3D array (b, h, w) of RGB values (e.g. as returned from rasterio's `.read()` method) transform : affine.Affine Transform of the input image as expressed by `rasterio` and the `affine` package s_crs : str/CRS Source CRS in which `img` is passed, expressed in any format permitted by rasterio. t_crs : str/CRS Target CRS, expressed in any format permitted by rasterio. resampling : [Optional. Default=Resampling.bilinear] Resampling method for executing warping, expressed as a `rasterio.enums.Resampling` method Returns ------- w_img : ndarray Warped image as a 3D array (b, h, w) of RGB values (e.g. as returned from rasterio's `.read()` method) w_transform : affine.Affine Transform of the input image as expressed by `rasterio` and the `affine` package )r)rr8s_crsrrrr w_transforms rrrns!F$CE*MEk  rc URupVn[5nURSUUU[URR 5UUS9n U R U5 SSS5 UR5n [XUS9n U R5nU Rn U RnSSS5 SSS5 SSS5 UW U4$!,(df  No=f!,(df  N4=f!,(df  N==f!,(df  NF=f)zK Warp an image. Returns the warped image and updated bounds and transform. r0r1N)r7r) r:r r>r?r6r@rBr readrr8) rr8rrrrKrJr/memfilemrastervrtrs rrrsiiGA!  \\ciinn%  MM#  \\^w7*EhhjMM F $  !!#  FE^ sS3C; CC;7 C*)C+C*3C; C C; C' #C** C8 4C;; D c [R"US[0S9nUR5 [R "UR 5n[R"U5RS5n[R"U5nUR5 UsSSS5 $!,(df  g=f![R[4a WRS:Xa%[R"SR!U55eUS:a)["R$"U5 US-n['XU5ng[R"S URS UR(S UR*35ef=f) a  Retry a url many times in attempt to get a tile and read the image Arguments --------- tile_url : str string that is the target of the web request. Should be a properly-formatted url for a tile provider. wait : int if the tile API is rate-limited, the number of seconds to wait between a failed request and the next try max_retries : int total number of rejected requests allowed before contextily will stop trying to fetch more tiles from a rate-limited API. Returns ------- array of the tile z user-agent)headersRGBANiz@Tile URL resulted in a 404 error. Double-check your tile url: {}rrz;Connection reset by peer too many times. Last message was: z Error: z for url: )requestsget USER_AGENTraise_for_statusioBytesIOcontentrr>convertr;asarrayclose HTTPErrorr status_codeformattimesleeprreasonr)r`rFrGrequest image_streamimagers rrrs9(\,,x, 1KL  " ZZ (LJJ|,44V>E>Q>Q=RS33:>>2B*W[[M*[\\\s2A B5 A B$ B5$ B2.B52B55AE1c U(d[X5up[X#5up#US:Xa [XX#5n[[[R "XX#U/555nU(a[ SXG4-5 U$)a Number of tiles required for a given bounding box and a zoom level Parameters ---------- w : float West edge s : float South edge e : float East edge n : float North edge zoom : int Level of detail verbose : Boolean [Optional. Default=True] If True, print short message with number of tiles and zoom. ll : Boolean [Optional. Default: False] If True, `w`, `s`, `e`, `n` are assumed to be lon/lat as opposed to Spherical Mercator. rWz0Using zoom level %i, this will download %i tiles)r9relenrgrhriprint)r/rCrDrEr(rr*ris rrrsf. a|a| v~qQ* RXXaA4&12 3E @D=PQ LrcURn[R"URURUR S5n[R"UR URUR S5nUR SURUS5- UR SURUS5- 4URUS5URUS544nU$)a8 Convert XY bounding box into the window of the tile raster Parameters ---------- bb : tuple (left, bottom, right, top) in the CRS of `rtr` rtr : RasterReader Open rasterio raster from which the window will be extracted Returns ------- window : tuple ((row_start, row_stop), (col_start, col_stop)) rrr.) rr;r<rrr:rr searchsorted)bbrtrrbbxiyiwindows rbb2wdwrs **C SXXsyy#))A, 7B SZZ#))A, 7B 11. . ! rrRSu?U0UV A A!78F Mrc Sn[RU-nX- S-nX- S-nS[R- S[R"[R"U[R-S- 55-[RS- - -nXE4$)ad Transform Spherical Mercator coordinates point into lon/lat NOTE: Translated from the JS implementation in http://dotnetfollower.com/wordpress/2011/07/javascript-how-to-convert-mercator-sphere-coordinates-to-latitude-and-longitude/ ... Arguments --------- x : float Easting y : float Northing Returns ------- ll : tuple lon/lat coordinates g@TXAgf@g@)r;piarctanexp)rPrQrMajorshiftlonlats rr9r9sy(F EEFNE )e C )e C "%%-3266#+2E+F!GG"%%RU+U VC 8Orc[R"X /5SSS2n[R"X/5SSS2n[R"U6n[R"U6n[R"[R"SU- 55n[R"[R"SU- 55n [R "X/5n [ U 5$)aAutomatically choose a zoom level given a desired number of tiles. .. note:: all values are interpreted as latitude / longitutde. Parameters ---------- w : float The western bbox edge. s : float The southern bbox edge. e : float The eastern bbox edge. n : float The northern bbox edge. Returns ------- zoom : int The zoom level to use in order to download this number of tiles. Nr-g@)r;sortsubtractceillog2minrl) r/rCrDrE lon_range lat_range lon_length lat_lengthzoom_lonzoom_latr(s rrere7s,"%I"%Ii(Ji(Jwwrwwy:567Hwwrwwy:567H 668& 'D t9rcrURSS5nSU;aURS5nSnOSnSnX0s=::aU::aU$ U(aSOSnS RX`5nU(aUS RX45- nOUS - nU(a)X:aU(a[R"U5 U$US - n[ U5e) ad Validate the zoom level and if needed raise informative error message. Returns the validated zoom. Parameters ---------- zoom : int The specified or calculated zoom level provider : dict auto : bool Indicating if zoom was specified or calculated (to have specific error message for each case). Returns ------- int Validated zoom level. min_zoomrmax_zoomTFinferred specifiedzDThe {0} zoom level of {1} is not valid for the current tile providerz (valid zooms: {0} - {1})..z This can indicate that the extent of your figure is wrong (e.g. too small extent, or in the wrong coordinate reference system))rrwarningswarnrm)r(rsrWrrmax_zoom_knownmodemsgs rrfrfZs(||J*HX<< +#8# $:;D P W W  C +228FF s  ?~ MM# O  I  S/rc[R"UVs/sHo"RUR4PM sn5nX3R SS9- nUSR upVnUS-R SS9up[R"XY-Xh-U4[RS9n [XA5H&upU upXX-US-U-2X-U S-U-2SS24'M( [R"UVs/sHn[R"U5PM sn5n[ USS2S45[ USS2S45[ USS2S45[ USS2S454unnnnU UUUU44$s snfs snf)a Merge a set of tiles into a single array. Parameters --------- tiles : list of mercantile.Tile objects The tiles to merge. arrays : list of numpy arrays The corresponding arrays (image pixels) of the tiles. This list has the same length and order as the `tiles` argument. Returns ------- img : np.ndarray Merged arrays. extent : tuple Bounding box [west, south, east, north] of the returned image in long/lat. r)axisr)r6Nr.r) r;rrPrQrr:maxzerosuint8ziprhr)rirxttile_xysindicesrJr/dn_xn_yrindarrrPrQrr{r|r}r~s rrprpsk*xxU3U##qssU34H1--GQiooGA!!   a (HC ((AGQWa( 9C(;>AEQUaK !a%1!4a 78)XXU3Uryy|U3 4F F1a4L F1a4L F1a4L F1a4L  D%u udE* **54$4s !E!' E&)rWNFrr.rT)rWNFrr.rTN)TF)T)@__doc__ __future__rrruuid mercantilerhratexitrrr"tempfilernumpyr;rasterior=PILrrjoblibr_Memoryr r rasterio.transformr rasterio.ior rasterio.vrtr rasterio.enumsrrr xyzservicesr__all__uuid4hexrmkdtempr$rrr&registerrrrdrobilinearrrrrrrr9rerfrpr%rrrs77@@  -$$*""%$ TZZ\-- -      #)".    pp   nb   #.*:M:M.,bAK@S@S$N"4,\\ F68 F2j/+r