(e V ddlmZddlZddlZddlZddlZddlZddlZddlZddl Z ddl Z ddl Z ddl Z ddl Z ddlmZddlmZmZmZmZmZddlmZddlmZddlmZmZeZ ddlZddlm Z ddl!m"Z#d Z$d Z&d Z'd Z(ejRdZ*dZ+dZ,GddZ-ej\dZ/Gdde0Z1Gdde2Z3iZ4gdZ5ejlddZ7e7rejpjse7Z9ejejpjue9dD]Z;e ej~e>jdZAeAe*jde;zS e3jeAZDeDe4eDj<e5jeDje jZHdZId ZJGd!d"e2ZKGd#d$e2ZLGd%d&e2ZMGd'd(e0ZNdid)eOd*eOd+efd,ZPd-eQd+dfd.ZR djd/ZSdid0ZTd1ejd2ed+eeVfd3ZWd1ejd2ed+eeejeeOffd4ZYd5eejd2ed+eejfd6ZZd7ejd+e[fd8Z\ djd9ejd:eejd;eejd2eed+ejf d<Z]d=Z^ djd>ejd?eeVd2eed+ejfd@Z_dAeOd+eOfdBZ` djd>ejd2eed+ejfdCZad7ejd+eQfdDZbdEZcdkdIZddJeQdKeeeQeQfdLeOdFdGd2ef dMZedJeQdLeOdFdGd+eeeeQeQffdNZfdFdGd+eQfdOZgd+ejfdPZidldQeeOd+eeOeOeffdRZjd2edSeOdLeOd+efdTZkGdUdHe2ZlGdVdGe2ZmGdWdXe0ZnGdYdZe2ZodSeOdLeOd2ed+eemeelffd[ZpGd\d]e2Zqd^ZrdSeOdLeOd2ed+dfd_ZsdSeOdLeOd2ed+dfd`ZtGdadbe2ZuGdcdde2Zve jd fdQeeOd+eQfdeZxdmdQeeOd+eQfdfZyezdgk(r%e jexe jd hyy#e%$rd Z$YwxYw#e0$rdZAYwxYw#e1$rYfe0$re*jde;zYwxYw)n)divisionN)partial)AnyListNoReturnOptionalTuple)uuid4) ElementTree)gdalosr)ImageTFz$Id$) averagenearbilinearcubic cubicsplinelanczos antialiasmodemaxminmedq1q3)allgoogle openlayersleafletmapmlnone gdal2tilesc|jdr(tj|ddk7rtd|yt j |dy)z;Wrapper for os.makedirs() that can work with /vsi files too/vsiirzCannot create T)exist_okN) startswithr MkdirRecursive Exceptionosmakedirs)paths 8/usr/lib/python3/dist-packages/osgeo_utils/gdal2tiles.pyr*r*_sJ v   tU +q 0nTF34 4 1 D4(c|jdr;tj|}tytj|j St jj|S)z>Wrapper for os.path.isfile() that can work with /vsi files toor$F) r&r VSIStatLstatS_ISREGrr)r+isfile)r+stat_ress r,r2r2hsJ v==& <||HMM**ww~~d##r-ceZdZdZdZdZy)VSIFilez1Expose a simplistic file-like API for a /vsi filec ||_||_yN)filenamef)selfr8r9s r,__init__zVSIFile.__init__vs  r-ctj|dt||jt|k7rt d|j zy)NzError while writing into %s)r VSIFWriteLlenr9r(r8)r:contents r,writez VSIFile.writezs@ ??7As7|TVV <G L9DMMIJ J Mr-N)__name__ __module__ __qualname____doc__r;rAr-r,r5r5ss;Kr-r5c#^K|jdr_tj||}|td|d| t ||tj |dk7rtd|yt ||y#tj |dk7rtd|wxYww)zDWrapper for open() built-in method that can work with /vsi files toor$Nz Cannot open z in rz Cannot close )r&r VSIFOpenLr(r5 VSIFCloseLopen)r8rr9s r,my_openrKs6" NN8T * 9l8*D?@ @ <(A& &q!Q&-z :;;'8T""q!Q&-z :;;'s;B-B 6B-(B**B-c eZdZy)UnsupportedTileMatrixSetNrBrCrDrFr-r,rMrMr-rMcFeZdZd dZdZdZdZdZede ddfd Z y) TileMatrixSetreturnNcd|_d|_d|_d|_d|_d|_d|_d|_d|_yr7) identifiersrs topleft_x topleft_y matrix_width matrix_height tile_size resolution level_countr:s r,r;zTileMatrixSet.__init__sD !r-c|j|jz|z d|zz }t||jz ||zz }t||j|j |jz|jzz z ||zz }||fSN)r[rZintrVrWrY)r:xyzoverriden_tile_sizerestxtys r,GeorefCoordToTileCoordz$TileMatrixSet.GeorefCoordToTileCoordsoo.1DD1M !dnn$/B)BC D NN((4>>9DOOKL (( *  2v r-ct|jD]:}|j|jz|z d|zz }||kDs+t d|dz cS|jdz S)?Maximal scaledown zoom of the pyramid closest to the pixelSize.r`rr=)ranger\r[rZr)r: pixelSizereirfs r,ZoomForPixelSizezTileMatrixSet.ZoomForPixelSizesit''( %A//DNN25HHAqDQC31a!e}$ %!##r-c|j|jz|z d|zz }||z|jz}||z|j|j|jz|jzz z}||fS)FConverts pixel coordinates in given zoom level of pyramid to EPSG:3857r`)r[rZrVrWrY)r:pxpyzoomrerfmxmys r,PixelsToMeterszTileMatrixSet.PixelsToMetersstoo.1DD4P #X & #X NNT//$..@4??R R 2v r-c|j||z||z||\}}|j|dz|z|dz|z||\}}||||fS)z6Returns bounds of the given tile in georef coordinatesr=)rw) r:rgrhrtreminxminymaxxmaxys r, TileBoundszTileMatrixSet.TileBoundssw(( $ $ $ $    d (( !V* * !V* *    d dD$''r-jcd|vsJd|vsJd|vsJt|dtsJtj}|j t |ddk(sJ|j xs|j}d}|jr|j}n4|jr$|jtjzdz }t}||_t |d|_t#|dD]1\}}d|vsJt|dtsJ|d}t%|dk(sJd |vsJd |vsJd |vsJ|d}|d } |d } |d } | | k7r t'd |d dz|z } |r| |} }|dk(r=||_| |_| |_| |_d|vsJd|vsJ|d|_|d|_||j(k7s| |j*k7r t'dt5|j,d |zz | z d| zkDr t'd| |j.k7s)t'dt%|d|_|S)NrT supportedCRS tileMatrixr? topLeftCornerr`scaleDenominator tileWidth tileHeightr=zOnly square tiles supportedgמY2? matrixWidth matrixHeightz"All levels should have same origing:0yE>z2Only resolutions varying as power-of-two supportedz%All levels should have same tile size) isinstancelistr SpatialReferenceSetFromUserInputstrEPSGTreatsAsLatLongEPSGTreatsAsNorthingEasting IsProjectedGetLinearUnits IsGeographic GetSemiMajormathpirQrUrT enumerater?rMrVrWr[rZrXrYabsr\) r~rUswapaxis metersPerUnittmsrnrrrVrWrrr[s r,parsezTileMatrixSet.parsesq   """q   !L/4000""$##C.(9$:;q@@@**,Q0O0O0Q ?? ..0M    ,,.83>MoQ|_-&q7' MAz"j0 00j94@ @@&7M}%* **%3 33*, ,,:- --%a(I%a(I";/I#L1JJ&./LMM#$67'AMQJ'0)9 Av ) ) !+ ) $ 222%333#-m#< $.~$>! -cmm1K23WXXs~~a0:=> ARR2L -2?K' Pa o. r-rRN) rBrCrDr;rirorwr} staticmethoddictrrFr-r,rQrQs=   $("99/99r-rQ)mercatorgeodeticrasterr ztms_MapML_APSTILE.jsonz tms_*.jsonrbutf-8z Cannot parse au globalmaptiles.py Global Map Tiles as defined in Tile Map Service (TMS) Profiles ============================================================== Functions necessary for generation of global tiles used on the web. It contains classes implementing coordinate conversions for: - GlobalMercator (based on EPSG:3857) for Google Maps, Yahoo Maps, Bing Maps compatible tiles - GlobalGeodetic (based on EPSG:4326) for OpenLayers Base Map and Google Earth compatible tiles More info at: http://wiki.osgeo.org/wiki/Tile_Map_Service_Specification http://wiki.osgeo.org/wiki/WMS_Tiling_Client_Recommendation http://msdn.microsoft.com/en-us/library/bb259689.aspx http://code.google.com/apis/maps/documentation/overlays.html#Google_Maps_Coordinates Created by Klokan Petr Pridal on 2008-07-03. Google Summer of Code 2008, project GDAL2Tiles for OSGEO. In case you use this class in your product, translate it to another language or find it useful for your project please let me know. My email: klokan at klokan dot cz. I would like to know where it was used. Class is available under the open-source GDAL license (www.gdal.org). cpeZdZdZddeddfdZdZdZdZd Z d Z d Z d Z d Z dZdZdZdZdZy)GlobalMercatoram TMS Global Mercator Profile --------------------------- Functions necessary for generation of tiles in Spherical Mercator projection, EPSG:3857. Such tiles are compatible with Google Maps, Bing Maps, Yahoo Maps, UK Ordnance Survey OpenSpace API, ... and you can overlay them on top of base maps of those web mapping applications. Pixel and tile coordinates are in TMS notation (origin [0,0] in bottom-left). What coordinate conversions do we need for TMS Global Mercator tiles:: LatLon <-> Meters <-> Pixels <-> Tile WGS84 coordinates Spherical Mercator Pixels in pyramid Tiles in pyramid lat/lon XY in meters XY pixels Z zoom XYZ from TMS EPSG:4326 EPSG:387 .----. --------- -- TMS / \ <-> | | <-> /----/ <-> Google \ / | | /--------/ QuadTree ----- --------- /------------/ KML, public WebMapService Web Clients TileMapService What is the coordinate extent of Earth in EPSG:3857? [-20037508.342789244, -20037508.342789244, 20037508.342789244, 20037508.342789244] Constant 20037508.342789244 comes from the circumference of the Earth in meters, which is 40 thousand kilometers, the coordinate origin is in the middle of extent. In fact you can calculate the constant as: 2 * math.pi * 6378137 / 2.0 $ echo 180 85 | gdaltransform -s_srs EPSG:4326 -t_srs EPSG:3857 Polar areas with abs(latitude) bigger then 85.05112878 are clipped off. What are zoom level constants (pixels/meter) for pyramid with EPSG:3857? whole region is on top of pyramid (zoom=0) covered by 256x256 pixels tile, every lower zoom level resolution is always divided by two initialResolution = 20037508.342789244 * 2 / 256 = 156543.03392804062 What is the difference between TMS and Google Maps/QuadTree tile name convention? The tile raster itself is the same (equal extent, projection, pixel size), there is just different identification of the same raster tile. Tiles in TMS are counted from [0,0] in the bottom-left corner, id is XYZ. Google placed the origin [0,0] to the top-left corner, reference is XYZ. Microsoft is referencing tiles by a QuadTree name, defined on the website: http://msdn2.microsoft.com/en-us/library/bb259689.aspx The lat/lon coordinates are using WGS84 datum, yes? Yes, all lat/lon we are mentioning should use WGS84 Geodetic Datum. Well, the web clients like Google Maps are projecting those coordinates by Spherical Mercator, so in fact lat/lon coordinates on sphere are treated as if the were on the WGS84 ellipsoid. From MSDN documentation: To simplify the calculations, we use the spherical form of projection, not the ellipsoidal form. Since the projection is used only for map display, and not for displaying numeric coordinates, we don't need the extra precision of an ellipsoidal projection. The spherical projection causes approximately 0.33 percent scale distortion in the Y direction, which is not visually noticeable. How do I create a raster in EPSG:3857 and convert coordinates with PROJ.4? You can use standard GIS tools like gdalwarp, cs2cs or gdaltransform. All of the tools supports -t_srs 'epsg:3857'. For other GIS programs check the exact definition of the projection: More info at http://spatialreference.org/ref/user/google-projection/ The same projection is designated as EPSG:3857. WKT definition is in the official EPSG database. Proj4 Text: +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs Human readable WKT format of EPSG:3857: PROJCS["Google Maps Global Mercator", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]], PROJECTION["Mercator_1SP"], PARAMETER["central_meridian",0], PARAMETER["scale_factor",1], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]]] rZrRNc||_dtjzdz|jz |_dtjzdzdz |_y)z*Initialize the TMS Global Mercator pyramidr`iRa@N)rZrrinitialResolution originShift)r:rZs r,r;zGlobalMercator.__init__sA"!"TWWw!6!Gtww;036r-c||jzdz }tjtjd|ztjzdz tjdz z }||jzdz }||fS)zKConverts given lat/lon in WGS84 Datum to XY in Spherical Mercator EPSG:3857f@Zv@)rrlogtanr)r:latlonrurvs r,LatLonToMeterszGlobalMercator.LatLonToMeterssm4## #e + XXdhhSDGG3e;< =5 Q $"" "U *2v r-c ||jz dz}||jz dz}dtjz dtjtj|tjzdz ztjdz z z}||fS)zMConverts XY point from Spherical Mercator EPSG:3857 to lat/lon in WGS84 Datumrrr`r)rrratanexp)r:rurvrrs r,MetersToLatLonzGlobalMercator.MetersToLatLonsD$$$-D$$$- gg 499TXXcDGGme&;<==# M O Cxr-ct|j|}||z|jz }||z|jz }||fS)rq Resolutionr)r:rrrsrtrfrurvs r,rwzGlobalMercator.PixelsToMeterssCood# #X(( ( #X(( (2v r-ct|j|}||jz|z }||jz|z }||fS)zCConverts EPSG:3857 to pyramid pixel coordinates in given zoom levelr)r:rurvrtrfrrrss r,MetersToPixelszGlobalMercator.MetersToPixelssEood#4###s *4###s *2v r-cttj|t|jz dz }ttj|t|jz dz }||fS)z9Returns a tile covering region in given pixel coordinatesr=rarceilfloatrZr:rrrsrgrhs r, PixelsToTilezGlobalMercator.PixelsToTileY2dnn 556: ; 2dnn 556: ;2v r-c.|j|z}|||z fS)z7Move the origin of pixel coordinates to top-left cornerrZ)r:rrrsrtmapSizes r,PixelsToRasterzGlobalMercator.PixelsToRasters!..D(7R<r-cR|j|||\}}|j||S)z+Returns tile for given mercator coordinates)rr)r:rurvrtrrrss r, MetersToTilezGlobalMercator.MetersToTiles/$$RT2B  R((r-c|j||jz||jz|\}}|j|dz|jz|dz|jz|\}}||||fS)z9Returns bounds of the given tile in EPSG:3857 coordinatesr=)rwrZ)r:rgrhrtryrzr{r|s r,r}zGlobalMercator.TileBoundssu((dnn)>>QSWX d(( !Vt~~ %Q$..'@$ ddD$''r-c|j|||}|j|d|d\}}|j|d|d\}}||||fS)zHReturns bounds of the given tile in latitude/longitude using WGS84 datumrr=r`)r}r) r:rgrhrtboundsminLatminLonmaxLatmaxLons r,TileLatLonBoundszGlobalMercator.TileLatLonBoundssaR.,,VAYq B,,VAYq B//r-c&|jd|zz S)zDResolution (meters/pixel) for given zoom level (measured at Equator)r`)rr:rts r,rzGlobalMercator.Resolution s%%D11r-cttD]'}||j|kDstd|dz cStdz Srkrr=rl MAXZOOMLEVELrrr:rmrns r,rozGlobalMercator.ZoomForPixelSizeF|$ %A4??1--1a!e}$ %ar-c|d|zdz |z fS)z8Converts TMS tile coordinates to Google Tile coordinatesr`r=rF)r:rgrhrts r, GoogleTilezGlobalMercator.GoogleTilesAtGaK2%%%r-cd}d|zdz |z }t|ddD]4}d}d|dz z}||zdk7r|dz }||zdk7r|dz }|t|z }6|S)z3Converts TMS tile coordinates to Microsoft QuadTreer`r=r)rlr)r:rgrhrtquadKeyrndigitmasks r,QuadTreezGlobalMercator.QuadTreesgkR tQ# "AEQ Pixels <-> Tiles WGS84 coordinates Pixels in pyramid Tiles in pyramid lat/lon XY pixels Z zoom XYZ from TMS EPSG:4326 .----. ---- / \ <-> /--------/ <-> TMS \ / /--------------/ ----- /--------------------/ WMS, KML Web Clients, Google Earth TileMapService tmscompatiblerZrRNch||_|rd|jz |_yd|jz |_y)Nrr)rZresFact)r:rrZs r,r;zGlobalGeodetic.__init__Rs0" !4>>1DL!4>>1DLr-cN|jd|zz }d|z|z }d|z|z }||fS)zLConverts lon/lat to pixel coordinates in given zoom of the EPSG:4326 pyramidr`rrr)r:rrrtrfrrrss r,LonLatToPixelszGlobalGeodetic.LonLatToPixels^s;llQW$Ci3 3h# 2v r-cttj|t|jz dz }ttj|t|jz dz }||fS)zDReturns coordinates of the tile covering region in pixel coordinatesr=rrs r,rzGlobalGeodetic.PixelsToTilefrr-cR|j|||\}}|j||S)z@Returns the tile for zoom which covers given lon/lat coordinates)rr)r:rrrtrrrss r, LonLatToTilezGlobalGeodetic.LonLatToTilems/$$S#t4B  R((r-c&|jd|zz S)zAResolution (arc/pixel) for given zoom level (measured at Equator)r`rrs r,rzGlobalGeodetic.Resolutionss||ag%%r-cttD]'}||j|kDstd|dz cStdz Srrrs r,rozGlobalGeodetic.ZoomForPixelSizexrr-c|jd|zz }||jz|zdz ||jz|zdz |dz|jz|zdz |dz|jz|zdz fS)z Returns bounds of the given tiler`rrr=)rrZ)r:rgrhrtrfs r,r}zGlobalGeodetic.TileBoundss~llQW$  # % +  # % * !Vt~~ % +c 1 !Vt~~ % +b 0   r-cL|j|||}|d|d|d|dfS)z1Returns bounds of the given tile in the SWNE formr=rrr`)r})r:rgrhrtbs r,rzGlobalGeodetic.TileLatLonBoundss3 OOBD )!adAaD!A$''r-r)rBrCrDrErboolrar;rrrrror}rrFr-r,rr/sI D 2htn 2 2t 2) &   (r-rceZdZdZddZdZy)Zoomifyz_ Tiles compatible with the Zoomify viewer ---------------------------------------- cJ||_||_||f}tj||z tj||z f}g|_|jj |g|_|j j ||d|kDs|d|kDrtj|ddz tj|ddz f}tj|d|z tj|d|z f}|jj ||j j ||d|kDr|d|kDr|jj|j jt|j|_ g|_ d|jd<td|jdzD]Y}|jj |j|dz d|j|dz dz|j|dz z[y)z'Initialization of the Zoomify tile treerr=r`N) rZ tileformatrrtierSizeInTilesappend tierImageSizefloorreverser? numberOfTierstileCountUpToTierrl)r:widthheightrZr imagesizetilesrns r,r;zZoomify.__init__s#$FO 59,-tyy)9K/LM " ##E*  !!),lY&)A,*BIaL1$45tzz)A,QRBR7STI )A,23 )A,23E  ' ' .    % %i 0lY&)A,*B $$& ""$!!5!56"$$%q!q$,,q01 A  " " ) )$$QU+A.1E1Ea!e1LQ1OO((Q/0  r-c |||j|dzz|j|z}tjj dt j |dz z|d|d|d|jS)z0Returns filename for tile with given coordinatesrz TileGroup%.0fr-.)rrr)r+joinrrr)r:rbrcrd tileIndexs r, tilefilenamezZoomify.tilefilenamesnD003A6669O9OPQ9RR ww|| djjS9 91doo 6  r-N)rjpg)rBrCrDrEr;r rFr-r,rrs &P r-rc eZdZy) GDALErrorNrNrFr-r,rrrOr-rmessagedetailsrRctjjdtjjd|z|r"tjjd|ztjdy)Nz4Usage: gdal2tiles.py [options] input_file [output] zgdal2tiles.py: error: %s z %s r`)sysstderrrAexit)rrs r,exit_with_errorrsTJJMNJJ1G;< g-.HHQKr-cache_in_bytescxdt|dz dz ztjd<tj|y)Nz%d GDAL_CACHEMAX)rar)environr SetCacheMax)rs r, set_cache_maxrs4#'^d-BT-I)J"JBJJ^$r-c |sg}|||c|d<|d<|d<||d<d|vr||d<d|vrt|ddz |d<d|vrt|dd z|d<|gk(rd |d<|4d } tj|jtj|d <nEd } t j ||||d<d|||dfz|d <||||\|d<|d<|d<|d<|dk(r d|zdz|d<n| d|z|d<nd|d<|j} | s| rd} nd} d|z} | r| d|zz } |D]G\} } }|| | |\}}}}t j | ||}| d|| ||d|||||d| || |f zz } I| dz } | S)z6 Template for the KML. Returns filled string. rgrhtzrrZ minlodpixelsr` maxlodpixelsrFxml_escaped_titleT realtiley %d/%d/%d.kmlsouthwestnortheastrr= drawOrderz../../ra. %(xml_escaped_title)s a %(north).14f %(south).14f %(east).14f %(west).14f %(minlodpixels)d %(maxlodpixels)d %(drawOrder)d %(realtiley)d.%(tileformat)s %(north).14f %(south).14f %(east).14f %(west).14f a  %d/%d/%d.%s %.14f %.14f %.14f %.14f %d -1 %s%d/%d/%d.kml onRegion z )rar EscapeStringtitle CPLES_XML GDAL2TilesgetYTileurl)rgrhrtileextrZtileswneoptionschildrenargstilekmlr0scxcyczcsouthcwestcnorthceastytiles r, generate_kmlr@sL )+R&DJT DJ D$%[T!"4 #4q#89^T!"4 #4q#89^2~!^ z$($5$5gmmT^^$T !&//B@[$2b"d;>O5P$P !CKBPRTVCW@W tF|T']DL QwFQJ[ F[[ ++C  CC     45  <& B'/B';$vu##BG4  ,               +# # & P A Hr-c |j}|j}|j}|jdk(r`td|dzD]M}t j |j ||j |d}|dk7s>td||fzOy|jdk(rctr\|jdr tdtj|||ftj} t|D]7}tj|j |dzdd||| dddd|f<9t!j"| d } | j%||ft j&} t(j*j-|r,t!j.|} t!j0| | | } i} |j2d k(r!|j4rd | d <n|j6| d <| j8||j2fi| y|jdk(rtj:}n@|jdk(rtj<}n|jdk(rtj>}n|jdk(rtj@}n|jdk(rtjB}n|jdk(rtjD}n|jdk(rtjF}n|jdk(rtjH}n_|jdk(rtjJ}n?|jdk(rtjL}n|jdk(rtjN}|jQd|tS|z ddd|tS|z f|jQdt jT||dd}|dk7rtd||fzyy)z-Scales down query dataset to the tile datasetrr=rz+RegenerateOverview() failed on %s, error %drr$zDOutputting to /vsi file systems with antialias mode is not supportedNRGBAWEBPTlosslessqualityrrrrrrrrrrrrFrrFrFrFrz'ReprojectImage() failed on %s, error %d)+ RasterXSize RasterCount resamplingrlr RegenerateOverview GetRasterBandrnumpy_availabler&r(numpyzerosuint8 gdalarrayBandReadAsArrayr fromarrayresizeLANCZOSr)r+existsrJ composite tiledriver webp_lossless webp_qualitysaveGRA_NearestNeighbour GRA_Bilinear GRA_CubicGRA_CubicSpline GRA_LanczosGRA_ModeGRA_MaxGRA_MinGRA_MedGRA_Q1GRA_Q3SetGeoTransformrReprojectImage)dsquerydstiler3r  querysizerZ tilebandsrnrfarrayimim1im0paramsgdal_resamplings r,scale_query_to_tilersdsh##I""I""IY&q)a-( A))%%a(&*>*>q*A9CaxA\SVDWW     { *  " "6 *V   Y 9=u{{Ky! A&66%%a!e,aIyE!Q'N __UF +iiI. > 77>>, '**\*C//#sC0C    '$$%)z"$+$8$8y!w11!!'64O !8 9\3  l&>&>&KK  l&>&>&KK*$ F?#F#'LE>"F#'LG#$F#'LWH$%XF#'L55hGB"$"3"3FF"CBQ"GKD%"$"3"3FF"CBQ"GKD%99$$&*E4%?'335* -- (..0&2D2D2F  w LL^   " - -k: Fr-c|s|Stj|}|jd}||S|jD]q\}}tj}|j dd|i|j ||jd|j}|jd|stj|jS)NGDALWarpOptionsOptionnamer) r fromstringfinditems TreeBuilderstartdataendcloseinserttostringdecode) vrt_string warp_optionsvrt_rootr3keyvaluetbelems r,add_gdal_warp_options_to_stringrrs %%j1Hmm-.G"((* U  $ $ & FC=)  xxxzq$      ) 0 0 22r-warped_vrt_dataset nodata_valuesc |gk7sJ|jdd}t|ddd}tj|}|j ddj |Dcgc] }t |c}|r[|jrOtjdtd d 5}|j|jddd d d |S|Scc}w#1swY|SxYw) z^ Takes an array of NODATA values and forces them on the WarpedVRT file dataset passed xml:VRTrNO_DATAYES) INIT_DESTUNIFIED_SRC_NODATA NODATA_VALUES z/Modified warping result saved into 'tiles1.vrt' tiles1.vrtwN) GetMetadatarr OpenSetMetadataItemr rr~rrrJrA)rrr3rcorrected_datasetrnr9s r,update_no_data_valuesrs B  #// :1=J0)5IJ$ *-%%="Aa3q6"AB7?? FG , $ A GG%11)  ;;/ ) MH o6J%*//W*D <==     "BHH1 %- HH]BGGGFF=FF? 88:D OOE4 ==!56L  "BHH^R GGC1 FF> 88:D  "BHHX ,-GGCLFF8 88:D    ) 0 0 22r-c^|jdvr|jdd}t|}tj|}|r[|j rOt jdtdd5}|j|jddddd|S|S#1swY|SxYw)z Handles dataset with 1 or 3 bands, i.e. without alpha channel, in the case the nodata value has not been forced by options )r=rrrz9Modified -dstalpha warping result saved into 'tiles1.vrt'rrN) rIrrr rr~rrrJrA)rr3rr9s r,'update_alpha_value_for_non_alpha_inputsrs%%/'33I>qA 1*= !YYz2 w LLT UlC( FA*66yA!DE F   F s 2$B""B,c|jdj}|jtjzs|j dk(s|j dk(r|j dz S|j S)zG Return the number of data (non-alpha) bands of a gdal dataset r=rr`)rL GetMaskBand GetMaskFlagsr GMF_ALPHArI)r alphabands r, nb_data_bandsrsk%%a(446I    !DNN 2   ! #   ! #""Q&&   r-c|g}|jdk(r*|jrdg}|Sdt|jzg}|S)NrCz LOSSLESS=TruezQUALITY=)rXrYrrZ)r3coptss r,_get_creation_optionsrsL EV#  $%E L #g&:&:";;> )<[$/ 1GG$)88)33 )11!&/   *=   s AQ55Q>base_tz base_tiles output_folderc |dz }|dddz }|dddz }tj|||}tjj |t |t ||d|j } |jrtj| |jr-t| r"|jrtjdytjd} |j} tj| } |jdz} | j!dd|j"zd|j"z| }| j!d|j"|j"| }g}|D]}|d}|d}tj|||}tjj |t |t ||d|j }t|sttj$|tj&}|dzdk(rd}n |j"}|j(r'|j*d k(r|dzdk(rd}n$|j"}n|dzdk(r |j"}nd}|j,| dz k(rx| j/d|d}|j1t3d g|j"|j"zz}|j5dd|j"|j"|| g |}n|j,| k7r t7d |j9dd|j"|j"}|j5|||j"|j"|t;t=d| dz |j?||sytA|||| |jBdk7rM| j/| |dtE|| dz}tjF|tjH||jr+tjd|ddj |z|jJrtM||}|tOtjj |d|||fzd5}|jQtS||||j |j"|||Dcgc] }|d|d|fc}jUddddyyycc}w#1swYyxYw)zLGenerating an overview tile from no more than 4 underlying tiles(base tiles)r=rr +Tile generation skipped because of --resumeNr rr`rrwrz'Unexpected number of bands in base tilerrrz.aux.xmlz build from zoom z , tiles: %srvr%rr)+r.r/r)r+r rrr~rrr1r2r rr rr)rZrrxyzrrIrAddBand bytearrayr.r(r*rrlrrsrJrr/Unlinkr/r0rKrAr@r-) r<r=r>rr3 overview_tz overview_tx overview_tyoverview_ty_realr  mem_driverr  out_driverrlrirjusable_base_tiles base_tilebase_txbase_ty base_ty_realbase_tile_path dsquerytiletileposxtileposytmp_dsr base_dataaux_xmlr9r9ts r,create_overview_tilerXsA+KQ-"a'KQ-"a'K!**; WM77<< K K#]%A%AB L  \"~~&. ?? LLF G%%e,J++K%%k2J++a/I A '' ']-D-D)DiG   M # #]%<%F NN 00=3J3JJKD   ''''$+  !K  $ $ 1EF F** q-))=+B+B      # #  # # 5IM23     +A@,D |L[( &4I'4R  + == ! - KK  )'+>*AUUV]G4   !"k;@P%QQ    ###%44%//8BC1!A$!g.C fWo    &D!  s80S$(S:S$S$$S-ci}|j|\}}}}t||dz dD]E}|dz } t||dzD],} | dz } | |f} | | f} | |vrg|| <|| j| .G|dz }t||dzD]E} | dz } tjj |t |t | }t|Gt|jS)z6Group base tiles that belong to the same overview tiler=r) tminmaxrlrr)r+r rr*rvalues)r<r>roverview_to_basestminxtminytmaxxtmaxyrhrGrgrFrL overview_tilerE tiledirnames r,group_overview_base_tilesrcCs !.!6!6w!?E5%E519b) ?Ag uqy) ?B'KRI(+6M$5535!-0 m , 3 3I > ? ?A+KE519%Ag ggll=#k2BC DTU  !((* ++r-cd}t|jdz |jdz dD]<}|j|\}}}}|dt ||z zdt ||z zzz }>|S)Nrr=r)rltmaxztminzrZr)r tile_numberrr]r^r_r`s r,count_overview_tilesrh`sKM''!+]-@-@1-DbIK%2%:%:2%>"ueUC ..1s55=7I3IJJ K r-c d}tj|dtz}|jddddtdd j tz |jd d d dt dd j t z |jddddd|jdddd|jddddd|jd d!d"d#d$|jd%d&d'dd(|jd)dd*d+,|jd-d.dd/d0,|jd1d2dd3d4,|jd5d6dd7d8,|jd9d:d;d<=|jd>dd?d@,|jdAdBdCdDd;dEF|jdGdHdIdJgdIddKLtj|dMdN}|jdOdPdQddR|jdSdTdQdUdV|jdWdXdYdZ|j|tj|d[d\}|jd]d^d_dtd`d j tz |jdadbdcdd|jdedfdgdh|jdidjdkdl|jdmdndodp|j|tj|dqdr}|jdsdtddu|j|tj|dvdw}|jdxdytdzd{||jd}d~dd|j||jddddddddddd |S)z*Prepare the option parser for input (argv)z*Usage: %prog [options] input_file [output]z%prog )versionz-pz --profilerchoicezGTile cutting profile (%s) - default 'mercator' (Google Maps compatible)rv)desttypechoiceshelpz-rz --resamplingrJz*Resampling method (%s) - default 'average'z-sz--s_srsrSRSz;The spatial reference system used for the source input data)rlmetavarroz-zz--zoomrtz4Zoom levels to render (format:'2-5', '10-' or '10').)rlroz-ez--resumer1 store_truez)Resume mode. Generate only missing files.)rlactionroz-az --srcnodatarxNODATAz4Value in the input dataset considered as transparentz-dz--tmscompatiblerzfWhen using the geodetic profile, specifies the base resolution as 0.703125 or 2 tiles at zoom level 0.z--xyzrAzs r, process_argsrLs> _F%%4%0MGT  R  4y1} A 79<$ H aJ *  H: U  4yA~Q (()9)9*)EFqI G# ??j ($(G ! ??   gmmM J   gll= I%gz=IG }g --r-rc@|js$tjj||_d}d}t |dr|j r}t |j trc|j jdd}|d}t|}t|dk(r'|d}|r"t|}||krtd||fz|}||g|_|jr0|jjds|xjdz c_ |jrM|}|jdr|dd}|xjtjj|dzz c_ |jd k(rts t!d d |j"d k(ryt%j&|j" t!d |j(sC|j*dks|j*dkDr t!dt|j*|_|j,rt.j1dt|zt.j1d|t.j1d|t.j1dt%j2dz dz z|S)Nrtrr=rr`z%max zoom (%d) less than min zoom (%d)/rrz2'antialias' resampling algorithm is not available.z/Install PIL (Python Imaging Library) and numpy.rCzWEBP driver is not availabledz+webp_quality should be in the range [1-100]z Options: %szInput: zOutput: z Cache: %d MBr)r,r)r+rhasattrrtrrrzrar?r(r0endswithrJrMrrXr rrYrZr~rr GetCacheMax) r3rr>rfreminmaxzoom_minzoom_maxout_paths r,rrxs/ ==((4  E EwGLLZ c5R##C+!9H  v;! ayHH 5=#?5%.P E5>GL{{7;;//4 s {{    S !}H rww''1C77 [( @ = V#    2 2 3 ; : ;$$##q(G,@,@3,F MN#&w';';#r3rRNcd|_d|_d|_d|_d|_d|_d|_d|_d|_d|_ d|_ d|_ d|_ d|_ d|_d|_d|_d|_d|_d|_d|_d|_d|_|j.r|j.|_|j0|_|j0dk(rd|_nd|_|j4r't7|t9j:||_nt9j:|_t>j@jC|j<tEtGdz|_$d|_%d |j,z|_&d |_'||_||_||_(|jPjRd k(r|j,|_&n-|jPjRd k(r|j,d z|_&|jPjT\|_+|_,|jPjZ|_-y)z%Constructor function - initializationNrr|pngwebp)dirz.vrtTrFrrr`).r7r6warped_input_datasetout_srs nativezoomrZtsizerrr3out_gtr2r9ominxomaxxomaxyrrr> isepsg4326in_srsrrZryrXr1rxr*tempfilemkdtemptmp_dirr)r+r rr tmp_vrt_filename scaledqueryrk overviewqueryr3rJrtrfrer/)r:rr>r3s r,r;zGDAL2Tiles.__init__s  $(!     "       !    $--DN!,,    & DL!DL ;; ] ##++ >DL#++-DL " T\\3uw<&;P Q T^^+ #$* << " "f ,!^^DN \\ $ $ 2!^^a/DN!%!2!2 DJ<<##r-c&'tjtjj_tjd_jst djzj s t djr/tjjtj}n t djjr9tjd|j|j|j fz|st#djz|j dk(rt#djz|j%d j'rt#d d jz|j%d j(tj*k7rt#d d jzt-|j}jjr9tjd|j|j|j fzt/|j\__t5j0j_d_jj:dk7rfj0s t#ddt=|s t#ddj0j?j6j?k7s|jAdk7rvtC|j0j6_|r(tEj8|j_n&tGj8j_j8rmjjrWtjdj8jj8jj8j fzj8s|_tjdjIjJj8tMj8_'d_(tSjT}|jWd|jYtRjZj6r||j?j6j?k(rQd_(j\d_.j\r+jjrtjdj\d_.j8j__0j`dj`dfdk7r t#dj`d_1j`dj8jj`d zz_2j`d _3j`d j8jj`d zz _4jjrNtjd!tkjbd"jhjdjffzjj:d#k(r-tmjn$_8jpjr_:twtydtz_>tydtzD]}jpjjbjh|\}}jpjjdjf|\}}td|td|}}td|zd z |td|zd z |}}||||fj||<jxjpjj`d tj8jj8jztjnz _BjRjpjj`d _Etjj_Etjj_Bjjrtjd%tjpjjbjhtjpjjdjftjd&jztjd'jjpjjfzyyjj:d(k(rtjjjn$_Kjjr_:twtydtz_>tydtzD]}jjjbjh|\}}jjjdjf|\}}td|td|}}td|d zzd z |td|zd z |}}||||fj||<jxjjj`d tj8jj8jztjnz _BjRjjj`d _Etjj_Etjj_BjjrEtjd)jbjhjdjffzyyjj:dk(rd*} tdtttj| j8jtjnz tj| j8jtjnz _Pjjr"tjd+jzjd_Bj8j_Etjj_En(jjkDrd jjz z} jjd,vrd-} n1jjd.vrjj} nd/} tjjJ|j8j| zj8j| z| 0tjjJ_j8j__0j_Ptwtydjd z_>twtydjd z_Stydjd zD]}d1j|z zjnz} d\}}ttjj8j| z d z }ttjj8j| z d z }tj| j|<||||fj||<j\r:j2r.tSjj0|fd2} | _:yd3_:ytjj:}d_:twtyd|jd z_>tyd|jd zD]}|jjbjh|jn\}}|jjdjf|jn\}}td|td|}}t|jd|zzd z |t|jd|zzd z |}}||||fj||<jy|jj`d tj8jj8jztjnz jn_BjS|jj`d jn_Etjj_Etjj_Bjjrtjd4jbjhjdjffztjd&jztjd5jzyy)6z=Initialization of the input raster, reprojection if necessaryr zBThe '%s' driver was not found, is it available in this GDAL build?zCThe 'MEM' driver was not found, is it available in this GDAL build?zNo input file was specifiedz"Input file: (%dP x %dL - %d bands)z/It is not possible to open the input file '%s'.rz"Input file '%s' has no raster bandr=zFPlease convert this file to RGB/RGBA and run gdal2tiles on the result.zFrom paletted file you can create RGBA file (temp.vrt) by: gdal_translate -of vrt -expand rgba %s temp.vrt then run: gdal2tiles temp.vrtzCPlease convert this file to 8-bit and run gdal2tiles on the result.zsTo scale pixel values you can use: gdal_translate -of VRT -ot Byte -scale %s temp.vrt then run: gdal2tiles temp.vrtz(Preprocessed file:(%dP x %dL - %d bands)NrzInput file has unknown SRS.zEUse --s_srs EPSG:xyz (or similar) to provide source reference system.z|There is no georeference - neither affine transformation (worldfile) nor GCPs. You can generate only 'raster' profile tiles.zEither gdal2tiles with parameter -p 'raster' or use another GIS software for georeference e.g. gdal_transform -gcp / -a_ullr / -a_srs)r3z0Projected file: tiles.vrt (%dP x %dL - %d bands)rFrTzKML autotest OK!r`rrrznGeoreference of the raster contains rotation or skew. Such raster is not supported. Please use gdalwarp first.rz#Bounds (output srs): %f, %f, %f, %f rrzBounds (latlong): %s, %szMinZoomLevel: %dzMaxZoomLevel: %d (%f)rz Bounds (latlong): %f, %f, %f, %fcXtj|tjdz Sr_)rlog10)rbs r,log2z#GDAL2Tiles.open_input..log2a szz!}tzz!}44r-zNative zoom of the raster: %d)rrr)rrrrrrr)rr resampleAlgrc&d j|z z jdz} jd| jz|zz}| j|zz} jjr2 j | jz|zz }| j|zz }n1 j | jz|zz}| j|zz} js0j||dd\}}j||dd\}}||||fS)Nr`r=r) rerrZr3rArrrr) rbrcrd pixelsizexr'r)r(r&rr:s r,rastertileswnez-GDAL2Tiles.open_input..rastertileswne sdjj1n- A> ;;q>A,>,KKD$..:"==D||'' $ Q-?*-L L %(C C $ Q-?*-L L %(C C??&(&7&7e&DRa&H e&(&7&7e&DRa&H e $t33r-cyN)rrrrrFrbrcrds r,z'GDAL2Tiles.open_input.. r-zBounds (georef): %f, %f, %f, %fzMaxZoomLevel: %d)Zr AllRegisterrrXr7r6r(rrrr3r~rrrHrrIrrLGetRasterColorTabler|r}rrrrrrrrrrrrrrrrrr3rr rrrrr/rrrrrrroundrrZrrr2rrlrrZrrrrfrorrerrrrrrrrarrrrJ Translaterrrr\rirXrY)r:rtrsrs4326rr]r^r_r`roversample_factorrrrrrs` @r, open_inputzGDAL2Tiles.open_inputIs| ++DOO< ++E2 ||T//" ||U  ??*.))DOOTEUEU*VM9: : <<   LL4!--!--!--  ADOOS   $ $ ) @4??R S  & &q ) = = ? X&)-8   & &q ) 2 2dmm C U&)-8 ) E <<   LL:!--!--!-- (7}dll'S$ T_' T\\B %)! <<  8 +;;1[ $M2N\ ))+t||/I/I/KK))+q0,=!4;; -)0E119dll1D-1X114<<1D-((T\\-A-A F11==11==11==(((5D % U#..  ! !4#<#< ,D,E,EF &&(t$&&s'E'EF <Su$)5%#? R  @zz!!]];;KKN11==11== DNN+ , zz!!]];;DKKNK  TZZ8 TZZ4DJ||## 6zz4::tzz4::FG$ \\ ! !X - 5"  $ 9 9 E E"'"7!8   $ 9 9 E E"'"7!8 DO(||## DLAs23 @"99JJ B  u #99JJ B  u #1e}c!Umu"3#3#3ae#;a#?G%%2-15Ju%*5%#? R  @zz! 11KKN11==11== DNN+ , NN zz! 11$++a.$..Q  TZZ8 TZZ4DJ||## 5zz4::tzz4::FG /$**<= /$**<= $r-c t|j|jjdk(rU|jj |j |j\}}|jj |j|j\}}td|td|}}td|td|}}||||f|_ |jjdvr|jjdk(r|jjr3tt j"j%|jdskt't j"j%|jdd5}|j)|j+j-d d d d |jjd vr|jjr4tt j"j%|jd sMt't j"j%|jd d5}|j)|j/j-d d d d n|jjd k(rl|j |j}}|j|j}}td|td|}}td|td|}}||||f|_ n\|jjdk(r<|j |j}}|j|j}}||||f|_ nd |_ |jjdvr|jjr3tt j"j%|jdskt't j"j%|jdd5}|j)|j1j-d d d d |jj2s|j|jjr3tt j"j%|jdskt't j"j%|jdd5}|j)|j5j-d d d d |jjdvr|jj2r|jjdk7r|jjd k7s|jj6r|jjr3tt j"j%|jdskt't j"j%|jdd5}|j)|j9j-d d d d |j:ra|j<Sg}|j>|j@\}}} } tC|| dzD]4} tC|| dzD] } |jE| | |j@g"6|j:r|jjr3tt j"j%|jdst't j"j%|jdd5}|j)tGd d d |jH|jJ|j<|j|j-d d d d y y y y y #1swYxYw#1swYSxYw#1swYxYw#1swYxYw#1swYxYw#1swYy xYw)z Generation of main metadata files and HTML viewers (metadata related to particular tiles are generated during the tile processing). rgECUrgECU@r)rrzgooglemaps.htmlrrN)rrz leaflet.htmlrgVgV@r)rrzopenlayers.htmlztilemapresource.xml)rr z mapml.mapmlr=zdoc.kml)&r*r>r3rrrrrrrrrr9r~r1r2r)r+r rKrAgenerate_googlemapsr-generate_leafletgenerate_openlayersrAgenerate_tilemapresourcergenerate_mapmlr/r2rZrfrlrr@r1rZ) r:r&r'r(r)r9r4xminyminxmaxymaxrbrcs r,generate_metadatazGDAL2Tiles.generate_metadata s. ##$ <<  : ---66tzz4::NKE4--66tzz4::NKE4lE2C4E4Ek513ud3C4EeT2DI &&*;;LL((J6||**&GGLL!3!35FG3! T%7%79JKTL 8 8 : A A' JKL ||%%);;||**&GGLL!3!3^D3! T%7%7H$I 5 5 7 > >w GHII \\ ! !Z /**djj%D**djj%DeU+S->4EdE*Ct,<4EeT2DI \\ ! !X -**djj%D**djj%DeT2DIDI << ! !%: :<<&&f T//1BC/GGLL!3!35FGHGGD446==gFGH    %LL''bggll4+=+=?TUV T//1FG I557>>wGH I LL " "&6 6    $$0%%3t||7Q7QLL''bggll4+=+=}MNd&8&8-H$O ?ST++-44W=> ? 88 1H%)\\$**%= "D$d4* 8tTAX.8AOOQ4::$678 8xx||**&GGLL!3!3Y?3! T%7%7CT( $ $ $ $ $ $ $ ( %fWo 328OLLII8HH I I  ? ? sI?/^" /^/;/^</_ rr*r.r/r1r1r2rrr}rrrZ geo_queryrarrHrrrArr r rrXr/rfrrrrr+)!r:r]r^r_r`r5rlrktcountti tile_detailsrrgrbrhr?r rrr nativesizer!r"r#r$r%r&r'r(rxsizeysizeconfs! r,generate_base_tileszGDAL2Tiles.generate_base_tilesw s ||!! KK0 1 <<   LL  LLC D LLC D LL &*\\$**%="ueU  & &''!+ NN <<   LL-0C0CC D LL94 5c%%-((QUU]1C-CD  ZZuqy) "B'',,t'9'93r7CGLK [ ! "uqy"-e BE519-d a"++BDLLA!ww||&&GG$dll3   <<''LL"fl1K!KL<<&&6,+?||++ %RS<<'':5 00R %(xml_escaped_title)s %(srs)s r=rzB r`g\mEAg?z% )r r+r3r,r-r9rZr1r0rrrrrlrfrerr)r:r5r7rds r,rz#GDAL2Tiles.generate_tilemapresourceD s=$($5$5 LL  %  !DH99@W tF|T']DL NN[!\\\!\\--\,,..Y <<  : -%DK \\ ! !Z /%DK \\  ,,,,DK \\,,224DKDK   tzz4::>2 A||##x/]\*t23dkk!nD %%3]L)1kAqD.@!DE%%3]L)1hAoqAB# *  r-cNi}tj|jjtj|d<d|d<|jj dk7r*|dxxd|jj zz cc<d|d<nd|d<|j \|d <|d <|d <|d <|j|d <|j|d<|j|d<|j|d<|jj|d<|jj|d<d|d<d|d<d|z}|S)a Template for googlemaps.html implementing Overlay of tiles for 'mercator' profile. It returns filled string. Expected variables: title, googlemapskey, north, south, east, west, minzoom, maxzoom, tile_size, tileformat, publishurl r#z'https://maps.googleapis.com/maps/api/js googlemapsurlrz?key=rgooglemapsurl_hintz`r&r'r(r)minzoommaxzoomrZrrra // Beginning of https://github.com/gavinharriss/google-maps-v3-opacity-control/blob/master/CustomTileOverlay.js // with CustomTileOverlay.prototype.getTileUrl() method customized for gdal2tiles needs. /******************************************************************************* Copyright (c) 2010-2012. Gavin Harriss Site: http://www.gavinharriss.com/ Originally developed for: http://www.topomap.co.nz/ Licences: Creative Commons Attribution 3.0 New Zealand License http://creativecommons.org/licenses/by/3.0/nz/ ******************************************************************************/ CustomTileOverlay = function (map, opacity) { this.tileSize = new google.maps.Size(256, 256); // Change to tile size being used this.map = map; this.opacity = opacity; this.tiles = []; this.visible = false; this.initialized = false; this.self = this; } CustomTileOverlay.prototype = new google.maps.OverlayView(); CustomTileOverlay.prototype.getTile = function (p, z, ownerDocument) { // If tile already exists then use it for (var n = 0; n < this.tiles.length; n++) { if (this.tiles[n].id == 't_' + p.x + '_' + p.y + '_' + z) { return this.tiles[n]; } } // If tile doesn't exist then create it var tile = ownerDocument.createElement('div'); var tp = this.getTileUrlCoord(p, z); tile.id = 't_' + tp.x + '_' + tp.y + '_' + z tile.style.width = this.tileSize.width + 'px'; tile.style.height = this.tileSize.height + 'px'; tile.style.backgroundImage = 'url(' + this.getTileUrl(tp, z) + ')'; tile.style.backgroundRepeat = 'no-repeat'; if (!this.visible) { tile.style.display = 'none'; } this.tiles.push(tile) this.setObjectOpacity(tile); return tile; } // Save memory / speed up the display by deleting tiles out of view // Essential for use on iOS devices such as iPhone and iPod! CustomTileOverlay.prototype.deleteHiddenTiles = function (zoom) { var bounds = this.map.getBounds(); var tileNE = this.getTileUrlCoordFromLatLng(bounds.getNorthEast(), zoom); var tileSW = this.getTileUrlCoordFromLatLng(bounds.getSouthWest(), zoom); var minX = tileSW.x - 1; var maxX = tileNE.x + 1; var minY = tileSW.y - 1; var maxY = tileNE.y + 1; var tilesToKeep = []; var tilesLength = this.tiles.length; for (var i = 0; i < tilesLength; i++) { var idParts = this.tiles[i].id.split("_"); var tileX = Number(idParts[1]); var tileY = Number(idParts[2]); var tileZ = Number(idParts[3]); if (( (minX < maxX && (tileX >= minX && tileX <= maxX)) || (minX > maxX && ((tileX >= minX && tileX <= (Math.pow(2, zoom) - 1)) || (tileX >= 0 && tileX <= maxX))) // Lapped the earth! ) && (tileY >= minY && tileY <= maxY) && tileZ == zoom) { tilesToKeep.push(this.tiles[i]); } else { delete this.tiles[i]; } } this.tiles = tilesToKeep; }; CustomTileOverlay.prototype.pointToTile = function (point, z) { var projection = this.map.getProjection(); var worldCoordinate = projection.fromLatLngToPoint(point); var pixelCoordinate = new google.maps.Point(worldCoordinate.x * Math.pow(2, z), worldCoordinate.y * Math.pow(2, z)); var tileCoordinate = new google.maps.Point(Math.floor(pixelCoordinate.x / this.tileSize.width), Math.floor(pixelCoordinate.y / this.tileSize.height)); return tileCoordinate; } CustomTileOverlay.prototype.getTileUrlCoordFromLatLng = function (latlng, zoom) { return this.getTileUrlCoord(this.pointToTile(latlng, zoom), zoom) } CustomTileOverlay.prototype.getTileUrlCoord = function (coord, zoom) { var tileRange = 1 << zoom; var y = tileRange - coord.y - 1; var x = coord.x; if (x < 0 || x >= tileRange) { x = (x % tileRange + tileRange) % tileRange; } return new google.maps.Point(x, y); } // Modified for gdal2tiles needs CustomTileOverlay.prototype.getTileUrl = function (tile, zoom) { if ((zoom < mapMinZoom) || (zoom > mapMaxZoom)) { return "https://gdal.org/resources/gdal2tiles/none.png"; } var ymax = 1 << zoom; var y = ymax - tile.y -1; var tileBounds = new google.maps.LatLngBounds( fromMercatorPixelToLatLng( new google.maps.Point( (tile.x)*256, (y+1)*256 ) , zoom ), fromMercatorPixelToLatLng( new google.maps.Point( (tile.x+1)*256, (y)*256 ) , zoom ) ); if (mapBounds.intersects(tileBounds)) { return zoom+"/"+tile.x+"/"+tile.y+".png"; } else { return "https://gdal.org/resources/gdal2tiles/none.png"; } } CustomTileOverlay.prototype.initialize = function () { if (this.initialized) { return; } var self = this.self; this.map.overlayMapTypes.insertAt(0, self); this.initialized = true; } CustomTileOverlay.prototype.hide = function () { this.visible = false; var tileCount = this.tiles.length; for (var n = 0; n < tileCount; n++) { this.tiles[n].style.display = 'none'; } } CustomTileOverlay.prototype.show = function () { this.initialize(); this.visible = true; var tileCount = this.tiles.length; for (var n = 0; n < tileCount; n++) { this.tiles[n].style.display = ''; } } CustomTileOverlay.prototype.releaseTile = function (tile) { tile = null; } CustomTileOverlay.prototype.setOpacity = function (op) { this.opacity = op; var tileCount = this.tiles.length; for (var n = 0; n < tileCount; n++) { this.setObjectOpacity(this.tiles[n]); } } CustomTileOverlay.prototype.setObjectOpacity = function (obj) { if (this.opacity > 0) { if (typeof (obj.style.filter) == 'string') { obj.style.filter = 'alpha(opacity:' + this.opacity + ')'; } if (typeof (obj.style.KHTMLOpacity) == 'string') { obj.style.KHTMLOpacity = this.opacity / 100; } if (typeof (obj.style.MozOpacity) == 'string') { obj.style.MozOpacity = this.opacity / 100; } if (typeof (obj.style.opacity) == 'string') { obj.style.opacity = this.opacity / 100; } } } // End of https://github.com/gavinharriss/google-maps-v3-opacity-control/blob/master/CustomTileOverlay.js custom_tile_overlay_jsa) // Beginning of https://github.com/gavinharriss/google-maps-v3-opacity-control/blob/master/ExtDraggableObject.js /** * @name ExtDraggableObject * @version 1.0 * @author Gabriel Schneider * @copyright (c) 2009 Gabriel Schneider * @fileoverview This sets up a given DOM element to be draggable * around the page. */ /* * Licensed under the Apache License, Version 2.0 (the "License"); * you may not use this file except in compliance with the License. * You may obtain a copy of the License at * * http://www.apache.org/licenses/LICENSE-2.0 * * Unless required by applicable law or agreed to in writing, software * distributed under the License is distributed on an "AS IS" BASIS, * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. * See the License for the specific language governing permissions and * limitations under the License. */ /** * Sets up a DOM element to be draggable. The options available * within {@link ExtDraggableObjectOptions} are: top, left, container, * draggingCursor, draggableCursor, intervalX, intervalY, * toleranceX, toleranceY, restrictX, and restrictY. * @param {HTMLElement} src The element to make draggable * @param {ExtDraggableObjectOptions} [opts] options * @constructor */ function ExtDraggableObject(src, opt_drag) { var me = this; var event_ = (window["GEvent"]||google.maps.Event||google.maps.event); var opt_drag_=opt_drag||{}; var draggingCursor_ = opt_drag_.draggingCursor||"default"; var draggableCursor_ = opt_drag_.draggableCursor||"default"; var moving_ = false, preventDefault_; var currentX_, currentY_, formerY_, formerX_, formerMouseX_, formerMouseY_; var top_, left_; var mouseDownEvent_, mouseUpEvent_, mouseMoveEvent_; var originalX_, originalY_; var halfIntervalX_ = Math.round(opt_drag_.intervalX/2); var halfIntervalY_ = Math.round(opt_drag_.intervalY/2); var target_ = src.setCapture?src:document; if (typeof opt_drag_.intervalX !== "number") { opt_drag_.intervalX = 1; } if (typeof opt_drag_.intervalY !== "number") { opt_drag_.intervalY = 1; } if (typeof opt_drag_.toleranceX !== "number") { opt_drag_.toleranceX = Infinity; } if (typeof opt_drag_.toleranceY !== "number") { opt_drag_.toleranceY = Infinity; } mouseDownEvent_ = event_.addDomListener(src, "mousedown", mouseDown_); mouseUpEvent_ = event_.addDomListener(target_, "mouseup", mouseUp_); setCursor_(false); if (opt_drag_.container) { } src.style.position = "absolute"; opt_drag_.left = opt_drag_.left||src.offsetLeft; opt_drag_.top = opt_drag_.top||src.offsetTop; opt_drag_.interval = opt_drag_.interval||1; moveTo_(opt_drag_.left, opt_drag_.top, false); /** * Set the cursor for {@link src} based on whether or not * the element is currently being dragged. * @param {Boolean} a Is the element being dragged? * @private */ function setCursor_(a) { if(a) { src.style.cursor = draggingCursor_; } else { src.style.cursor = draggableCursor_; } } /** * Moves the element {@link src} to the given * location. * @param {Number} x The left position to move to. * @param {Number} y The top position to move to. * @param {Boolean} prevent Prevent moving? * @private */ function moveTo_(x, y, prevent) { var roundedIntervalX_, roundedIntervalY_; left_ = Math.round(x); top_ = Math.round(y); if (opt_drag_.intervalX>1) { roundedIntervalX_ = Math.round(left_%opt_drag_.intervalX); left_ = (roundedIntervalX_1) { roundedIntervalY_ = Math.round(top_%opt_drag_.intervalY); top_ = (roundedIntervalY_opt_drag_.toleranceX||(currentX_-(left_+src.offsetWidth))>opt_drag_.toleranceX)||((top_-currentY_)>opt_drag_.toleranceY||(currentY_-(top_+src.offsetHeight))>opt_drag_.toleranceY)) { left_ = originalX_; top_ = originalY_; } } if(!opt_drag_.restrictX&&!prevent) { src.style.left = left_ + "px"; } if(!opt_drag_.restrictY&&!prevent) { src.style.top = top_ + "px"; } } /** * Handles the mousemove event. * @param {event} ev The event data sent by the browser. * @private */ function mouseMove_(ev) { var e=ev||event; currentX_ = formerX_+((e.pageX||(e.clientX+document.body.scrollLeft+document.documentElement.scrollLeft))-formerMouseX_); currentY_ = formerY_+((e.pageY||(e.clientY+document.body.scrollTop+document.documentElement.scrollTop))-formerMouseY_); formerX_ = currentX_; formerY_ = currentY_; formerMouseX_ = e.pageX||(e.clientX+document.body.scrollLeft+document.documentElement.scrollLeft); formerMouseY_ = e.pageY||(e.clientY+document.body.scrollTop+document.documentElement.scrollTop); if (moving_) { moveTo_(currentX_,currentY_, preventDefault_); event_.trigger(me, "drag", {mouseX: formerMouseX_, mouseY: formerMouseY_, startLeft: originalX_, startTop: originalY_, event:e}); } } /** * Handles the mousedown event. * @param {event} ev The event data sent by the browser. * @private */ function mouseDown_(ev) { var e=ev||event; setCursor_(true); event_.trigger(me, "mousedown", e); if (src.style.position !== "absolute") { src.style.position = "absolute"; return; } formerMouseX_ = e.pageX||(e.clientX+document.body.scrollLeft+document.documentElement.scrollLeft); formerMouseY_ = e.pageY||(e.clientY+document.body.scrollTop+document.documentElement.scrollTop); originalX_ = src.offsetLeft; originalY_ = src.offsetTop; formerX_ = originalX_; formerY_ = originalY_; mouseMoveEvent_ = event_.addDomListener(target_, "mousemove", mouseMove_); if (src.setCapture) { src.setCapture(); } if (e.preventDefault) { e.preventDefault(); e.stopPropagation(); } else { e.cancelBubble=true; e.returnValue=false; } moving_ = true; event_.trigger(me, "dragstart", {mouseX: formerMouseX_, mouseY: formerMouseY_, startLeft: originalX_, startTop: originalY_, event:e}); } /** * Handles the mouseup event. * @param {event} ev The event data sent by the browser. * @private */ function mouseUp_(ev) { var e=ev||event; if (moving_) { setCursor_(false); event_.removeListener(mouseMoveEvent_); if (src.releaseCapture) { src.releaseCapture(); } moving_ = false; event_.trigger(me, "dragend", {mouseX: formerMouseX_, mouseY: formerMouseY_, startLeft: originalX_, startTop: originalY_, event:e}); } currentX_ = currentY_ = null; event_.trigger(me, "mouseup", e); } /** * Move the element {@link src} to the given location. * @param {Point} point An object with an x and y property * that represents the location to move to. */ me.moveTo = function(point) { moveTo_(point.x, point.y, false); }; /** * Move the element {@link src} by the given amount. * @param {Size} size An object with an x and y property * that represents distance to move the element. */ me.moveBy = function(size) { moveTo_(src.offsetLeft + size.width, src.offsetHeight + size.height, false); } /** * Sets the cursor for the dragging state. * @param {String} cursor The name of the cursor to use. */ me.setDraggingCursor = function(cursor) { draggingCursor_ = cursor; setCursor_(moving_); }; /** * Sets the cursor for the draggable state. * @param {String} cursor The name of the cursor to use. */ me.setDraggableCursor = function(cursor) { draggableCursor_ = cursor; setCursor_(moving_); }; /** * Returns the current left location. * @return {Number} */ me.left = function() { return left_; }; /** * Returns the current top location. * @return {Number} */ me.top = function() { return top_; }; /** * Returns the number of intervals the element has moved * along the X axis. Useful for scrollbar type * applications. * @return {Number} */ me.valueX = function() { var i = opt_drag_.intervalX||1; return Math.round(left_ / i); }; /** * Returns the number of intervals the element has moved * along the Y axis. Useful for scrollbar type * applications. * @return {Number} */ me.valueY = function() { var i = opt_drag_.intervalY||1; return Math.round(top_ / i); }; /** * Sets the left position of the draggable object based on * intervalX. * @param {Number} value The location to move to. */ me.setValueX = function(value) { moveTo_(value * opt_drag_.intervalX, top_, false); }; /** * Sets the top position of the draggable object based on * intervalY. * @param {Number} value The location to move to. */ me.setValueY = function(value) { moveTo_(left_, value * opt_drag_.intervalY, false); }; /** * Prevents the default movement behavior of the object. * The object can still be moved by other methods. */ me.preventDefaultMovement = function(prevent) { preventDefault_ = prevent; }; } /** * @name ExtDraggableObjectOptions * @class This class represents the optional parameter passed into constructor of * ExtDraggableObject. * @property {Number} [top] Top pixel * @property {Number} [left] Left pixel * @property {HTMLElement} [container] HTMLElement as container. * @property {String} [draggingCursor] Dragging Cursor * @property {String} [draggableCursor] Draggable Cursor * @property {Number} [intervalX] Interval in X direction * @property {Number} [intervalY] Interval in Y direction * @property {Number} [toleranceX] Tolerance X in pixel * @property {Number} [toleranceY] Tolerance Y in pixel * @property {Boolean} [restrictX] Whether to restrict move in X direction * @property {Boolean} [restrictY] Whether to restrict move in Y direction */ // End of https://github.com/gavinharriss/google-maps-v3-opacity-control/blob/master/ExtDraggableObject.js ext_draggable_object_jsav* %(xml_escaped_title)s %(googlemapsurl_hint)s
Generated by GDAL2Tiles, Copyright © 2008 Klokan Petr Pridal, GDAL & OSGeo GSoC
) r r+r3r,r-rr9rfrerZr1r0rr:r5r7s r,rzGDAL2Tiles.generate_googlemaps s[$($5$5 LL  %  !!J_ << ! !%; ;  !Wt||/E/E%E E !)+D% &s $ DH99@W tF|T']DL**Y**Y NN[!\\\!\\--\ LL22[v  $ x  % H q der  tr-ci}|jjjdd|d<tj|jjtj |d<|j \|d<|d<|d<|d<|d|dzd z |d <|d|dzd z |d <|j|d <|j|d <|j|d<|j|d<|j|d<|jj|d<|jjjdd|d<|jjrd|d<nd|d<d|z}|S)z Template for leaflet.html implementing overlay of tiles for 'mercator' profile. It returns filled string. Expected variables: title, north, south, east, west, minzoom, maxzoom, tile_size, tileformat, publishurl "z\"double_quote_escaped_titler#r&r'r(r)r centerlon centerlatrr beginzoomrZrrrrrr=aN %(xml_escaped_title)s
)r3r,replacer r+r-r9rfrerZr1r0rrArs r,rzGDAL2Tiles.generate_leafletsj-1\\-?-?-G-GU-S )*$($5$5 LL  %  !DH99@W tF|T']DL!']T']:cA[!&\DL8C?[**Y**Y JJ[ NN[!\\\!\\--\ LL22::3F[ <<  DKDKj VWk  ^r-c  i}tj|jjtj|d<|jj |d<|j |d<|j|d<|j|d<|j|d<|jj|d<|jj|d<|jjrd |d <nd |d <|j|d <|j|d <|j|d<|j |d<|j|jzdz |d<|j|j zdz |d<d|z}|jj"dk(s|jj"dk(r|d|zz }|jj"dk(r |d|zz }n|jj"dk(r|jj$rd|jz }nd|jz }t'|jdzDcgc] }|d|zz  }}ddj)d|Dzdz|d<|jjr'|jj$rd |d!<nd"|d!<d#|d$<n d%|d!<d&|d$<|d'|zz }n|jj"d(k(rd|j*z|j,dz}t'|jdzDcgc] }|d|zz  }}||j |d)<ddj)d*|Dzdz|d<d+|j|j|j|j fz|d,<|jjr$d-|j|j fz|d!<d#|d$<n#d-|j|jfz|d!<d&|d$<|d.|zz }nt.|jj"j0}t'|jdzDcgc] }|d|zz  }}||j |d)<ddj)d/|Dzdz|d<ddj)fd0t't3|Dzdz|d1<|jjr$d-j4j6fz|d!<d#|d$<n=d-j4j6j0jzz fz|d!<d&|d$<d+j4j6j8j0zjzz j4j:j0zjzzj6fz|d,<|d2|zz }|d3|zz }|jj"d4vrM|d|d5<|jj"dk(r#|jj$r |d5xxdz cc<|d6|zz }n|d7|zz }|jj"dk(r|d8z }n~|jj"dk7re|j<rY|j<j?r?|j<jAd9d:k(r!|d;|j<jCd9zz }|dz }|Scc}wcc}wcc}w)?z Template for openlayers.html, with the tiles as overlays, and base layers. It returns filled string. r#rrrrZrrrrsign_yrrrrrr`center_xcenter_yac %(xml_escaped_title)s
Generated by GDAL2Tiles    
)"r r+r3r,r-rrfrerZr1r0rrArrrrrrrlr rrrr[r?rVrWrYrXrrrr)r:r5r7base_resrnr3rs @r,rzGDAL2Tiles.generate_openlayers+s$$($5$5 LL  %  !,,..Y**Y**Y NN[!\\\!\\--\ LL22[ <<  DN DN W  W  W  W  JJ3q8Z JJ3q8Z! DE"  L <<  : -1E1E1S 'PQ(* AX <<  : - $% A,\\ ! !Z /||)) 4>>1 4>>149$**q.4IJq8ad?JKJchhD DDDsJ  ||<<--%0DN%1DN$2[!!-X$8[! :; AB\\ ! !X -T__- A>H49$**q.4IJq8ad?JKJ(4DNchhD DDDsJ  'B    E'D" #||!0DJJ 3K!KX$2[!!0DJJ 3K!KX$8[! 45 A@--.C~~H49$**q.4IJq8ad?JKJ(4DNchhD DDDsJ  (("3{#34    ||!0CMM3==3Q!QX$2[!!0MMMMCNNS]]$BB4"X%9[!&A   1 1CNN BS]] RR  0 03>> ACMM QQ E'D" # <= AD  9    <<  #; ; $YD ||##z1dll6P6P[!Q&! ( A + A <<  : - , ,A \\ ! !Z / KK++-KK006&@UW[WbWbWsWsX   <<  #; ; ; ;A    uKfKjKs+[[#["c |jjr|jj}ntjdd}t |dj j d}|jjdk(rd}n2|jjdk(rd}n|jj}|jd |}|jd |jjr|jjnd }|j|j\}}}}|jd t|}|jd ttj||j|j}|jdt|}|jdttj||j|j}|jdt|j}|jdt|j}|jdt|j}|jdt|j }|S)Nr ztemplate_tiles.mapmlrrrOSMTILErWGS84z${TILING_SCHEME}z${URL}z./z ${MINTILEX}z ${MINTILEY}z ${MAXTILEX}z ${MAXTILEY}z ${CURZOOM}z ${MINZOOM}z ${MAXZOOM}z ${TILEEXT})r3rr FindFilerJreadrrr%r0rZrerr.r/rfr1)r:templater7 tiling_schemer]r^r_r`s r,rzGDAL2Tiles.generate_mapmls << & &||22H}}V-CDH 4 % % ' . .w 7 <<  : -%M \\ ! !Z /#M LL00M II(- 8 IIhDLL4D4D 0 0$ O%)\\$**%="ueU IImSZ 0 II 3z225$**dllST  IImSZ 0 II 3z225$**dllST  IIlC O 4 IIlC O 4 IIlC O 4 IIlC $5 6r-c|jrS|jdk7rD|jdvr d|zdz |z St|j}|jd|zzdz |z S|S)z Calculates the y-tile number based on whether XYZ or TMS (default) system is used :param ty: The y-tile number :param tz: The z-tile number :return: The transformed y-tile number rr<r`r=)rArrrY)rhrr3rs r,r/zGDAL2Tiles.getYTileso ;;7??h6"::2 R'')C!!ArE)A-  r-r)r)rBrCrDrOptionsr;rrr r rr r rrrrrrrr/rFr-r,r.r.sH$3H$sH$WH$QUH$Vp>d zx_"U;Z8H+H%I_"B*BXA#AFU SU nK#KZhShT Br-r.ct|||}|j|j|j\}}||fSr7)r.rrr )rr>r3r"rrs r,worker_tile_detailsrIsHJ w?J  "","@"@"BM< , &&r-cNeZdZejfdeddfdZddZd deddfdZy) ProgressBar total_itemsrRNc.||_d|_||_y)Nr)rL nb_items_done progress_cbk)r:rLrOs r,r;zProgressBar.__init__s&(r-c*|jdddy)Nrr)rOr]s r,rzProgressBar.starts !R&r-nb_itemsc|xj|z c_t|j|jz }|j|ddy)Nr)rNrrLrO)r:rQprogresss r, log_progresszProgressBar.log_progresss@ h&++,t/?/?? (B-r-r)r=) rBrCrDr TermProgress_nocbrar;rrTrFr-r,rKrKs76:6L6L)C)QU) '.S..r-rKcjdk(rt}|j}|Sjdk(r#tj}|j}|Sjdk(rt j }|jd|jt jjr{jrot j }|jt j|jjt j||fd}|}|Sd}|Sd}|S)Nrrrrcd j|z z jdz} jd| jz|zz}| j|zz} jr5 jd| jz|zz }| j|zz }n1 j| jz|zz}| j|zz} j s0j ||dd\}}j ||dd\}}||||fS)Nr`r=rr)rerrZrArrr) rbrcrdrr'r)r(r&rr3rs r,rz%get_tile_swne..rastertileswnesH---12]5P5PQR5SS"//2-111J>?m55 BB;;%33A6m555 BC"M$;$;j$HHE&++a-2I2I.IJ.VV"M$;$;j$HHE$11"$"3"3D%"@!"DKD%"$"3"3D%"@!"DKD%dE4//r-cyrrFrs r,rzget_tile_swne..rr-)rrrrrr rrrrr/rrr) rr3r tile_swnerrrrrs `` @r,r0r0s&*$!#-- \ [ J &!'"7"78-- V U H $&&(t$&&s'E'EF   !9!9))+F  ) )#*H*H I  !9!9 :--fg>B 04'I  5I  r-c |jrtjdt|||\}}|jrtjd|js0|js$t t |}|j|D]8}t|||jr|jr)j:ttddrt` |jsIt|}|rr3r rbase_progress_barr2r,overview_progress_barr<base_tile_groupsr=s r,single_threaded_tilingrdss /0,ZPD, 34 ??7=='L(9:!#- {+w}}  * * , - {K.  ! ==$T*  KK4 5??(3E(:%%++-TZZ454WmTR* 5J *mT7 S??7==%224 55 MM"''//$--01r-c |jxsd}|jrtjdt |||\}}|jrtjd|js0|j s$t t|}|jtdtdt||z}|jtt|||D],} |jr|j rj.|j sIt|} | rrr3)rvr~rrrIrurKr?rrrimap_unorderedrr;rTrhrrlrerfrcrXr^r_r)r+r`r)rr>r3poolrvr rrarg_r,rbr<rcs r,multi_threaded_tilingrkEs'',1L /0,ZPD, 34 ??7=='L(9:!As3L 1\ ABCI   $'!-w}}  * * , - ==$T*  KK4 5??(3E(:%%++-TZZ454WmTR3sC(8$9\$IJK $$ $+"   %  5A??7==%224 55" MM"''//$--01r-ceZdZdZdZy) UseExceptionscxtj|_|jstjyyr7)r GetUseExceptionsold_used_exceptionsrmr]s r, __enter__zUseExceptions.__enter__}s-#'#8#8#: ''    (r-cF|jstjyyr7)rpr DontUseExceptionsr:rmrrs r,__exit__zUseExceptions.__exit__s''  " " $(r-N)rBrCrDrqrurFr-r,rmrm|s ! %r-rmceZdZdZdZdZy) DividedCachec||_yr7)rv)r:rvs r,r;zDividedCache.__init__s (r-ctj|_tdt j |j|j z }t|y)Ni)r rgdal_cache_maxrrrrvr)r:gdal_cache_max_per_processs r,rqzDividedCache.__enter__sH"..0%( D$7$7$:K:K$KL& " 01r-c.t|jyr7)rrzrts r,ruzDividedCache.__exit__sd))*r-N)rBrCrDr;rqrurFr-r,rwrws)2+r-rwc6tt|D]<}||dk(s |dzt|ks||dztj||dz<>d|vrddlm}ddlm}t5||jd5}| ddddddyt|jd |_ t|||jj| cdddcdddSt|| S#1swYnxYwdddy#1swYyxYw) Nz--configr`r=rwr)MPI)MPICommExecutor)rootT) unorderedr)rlr?r)rmpi4pyr~mpi4py.futuresrrm COMM_WORLDrryrhsubmainGet_size)rrrnr~rris r,mainrs3t9 2 7j QUSY%6&*1q5kBJJtAE{ #2$2 _ ocnn1E |   #*$((d"CD dCNN335HX    t.>??     s1.DC:DAC: D:D ?DDc tj|}|yt|dd|\}}}|r||_|jxsd}t 5|t ||||nd|dk(rt |||nQddl}t|dsd|_ ddl m } t|5| |5}t ||||dddddddddy#1swYxYw#1swYxYw#1swYyxYw)Nrr=r__spec__)Pool)r) r GeneralCmdLineProcessorrrvrmrkrd__main__rrmultiprocessingrrw) rri pool_sizerrr>r3rvrrs r,rrs   ' ' -D |)5 QR#3*&J w('',1L P   !*mWd K Q  ":}g F 8Z0$(! ,l+ PTL-I PT%j-$O P PP"  P P P PP" s=A C' C'C6C>C'C CC$ C''C0rr)rr7)rr r2r rRN)F)NrF)| __future__r contextlibglobjsonrrrr)r^r0rr threading functoolsrtypingrrrrr uuidr xml.etreer osgeor r rGrNPILrosgeo.gdal_array gdal_arrayrQrM ImportErrorrrr getLoggerrr*r2r5contextmanagerrKr(rMobjectrQrrrBr8r+r`r  tmsfilenamerJrCrloadsrr~errorrrrTrlocalr__doc__globalmaptilesrrrrrrrrarr@rsDatasetrrrrrrrrrrrrrrr;rXrcrhrrrrr r rr.rIrKr0rdrkrmrwrrrrBrrFr-r,rsN    77!  (O  M   < ()$ K K  # # y zFzz 1  4==!9 : ggooh'G tyyg|!DE, K&++-  4;;w/0A 9 LL;6 7  %%a(C "%s~~CNN+#,&ioo @ XVXv^(V^(B5 f5 p  S3%#%$%AEC LcL w4PU;B $<< $*1 $ 8C(( )8C= 89 $F,,-8? c""#.dllt"& C,,Cs++,C S)) *Cg  C  \\ CL30"&, ,;,g , \\ ,^63S63S63tDH /7/@ \\0 4<< C IXU UU38_%UU! U  Up, ,!$,5B, $uS#X  ,: #` x,,` F).tCy).U3WCT=U).X< <"%<69< <~GG:"4&"4J i J#J#ZF''$''29' ;Z( ()' .& .1h+2+2$'+229+2 +2\4242$'422942 42n%F%+6+"88e@tCy@@6 $s) PS H z CHHT#((T 23eLOn A (    LL;6 7  s<"Q(!Q6:R(Q32Q36RRR( R('R(