njSSKJr SSKrSSKrSSKJrJrJr \Rr "SS\ 5r "SS\ 5r "SS \ 5r "S S \ 5rS rS!S jrSrSrSrSrSrSrSrSrSrSrSrSrS!SjrSrSrS"Sjr "SS5r!\"S :Xa#\RF"\ "\RH55 gg)#)print_functionN)gdalogrosrc$\rSrSrSrSrSrSrg)AffineTransformDecorator)z@A class providing some useful methods for affine TransformationscSUlSUlXlURSUlURSUlURS:aU=RS-slURSUlURSUlg)Nr)lrxlry geotransformscaleXscaleYulxuly)self transforms S/var/www/html/pdf-tiff/venv/lib/python3.13/site-packages/osgeo_utils/gdal_retile.py__init__!AffineTransformDecorator.__init__,sw%''* ''* ;;? KK2 K$$Q'$$Q'c/n/nURU-nURU-nURUR5 URUR5 URURU-5 URUR5 URURU-5 URURU-5 URUR5 URURU-5 X4/$N)rrappendrr)rwidthheightxlistylistwhs r pointsFor"AffineTransformDecorator.pointsFor7s KK%  KK&  TXX TXX TXX\" TXX TXX\" TXX\" TXX TXX\"~r)rrrrrrrN)__name__ __module__ __qualname____firstlineno____doc__rr%__static_attributes__rrrr)sJ (rrc*\rSrSrSrSrSrSrSrg) DataSetCacheHz A class for caching source tilesc.SUl/Ul0Ulg)N) cacheSizequeuedictrs rrDataSetCache.__init__Ks  rcXR;aURU$[R"U5nUc [S[-[ R S9 g[UR5UR:Xa(URRS5nURU URRU5 X RU'U$)NzError opening: %sfiler r) r5rOpenprint NameErrorsysstderrlenr4r3popr)rnameresulttoRemoves rgetDataSetCache.getPs 99 99T? "4 > % 1 C tzz?dnn ,zz~~a(H (# $ $ rcPURR5HnAM U?U?gr)r5valuesr4)rdatasets r__del__DataSetCache.__del___s&yy'')G* J Ir)r3r5r4N) r'r(r)r*r+rrErJr,r-rrr/r/Hs* rr/c$\rSrSrSrSrSrSrg) tile_infofz A class holding info how to tilec$XlX lX0lX@lSUlX:a-U=R[ X- X5- -S- X5- - 5- slSUlX$:a-U=R [ X$- XE- -S- XE- - 5- slXPlg)Nr )rr tileWidth tileHeight countTilesXint countTilesYoverlap)rxsizeysizerPrQrUs rrtile_info.__init__is  "$     "i&9:Q>9CVW!       #z';>%335!//2AAC &{'B'B'DEjj jj  & &q ) = = ? >hhjDGDG# "4::&E(66 **, GGEN' $$--/99;!9!9!9!9488 3t{{BCD UDHHtxx$74;;#FGHI rc U?U?gr)rerfr6s rrJmosaic_info.__del__s J  rc  URR5R5 URR5RXX45 /nSnURR5R 5nUcOUR U5 UcUR 5R5nOdUR 5R5n[USUS5[USUS5[USUS5[USUS54nMUcg[XS5[X6S5[X&S5[XFS54nURR5RS5 [X1- UR- S-5n [X$- UR- S-5n URRSXUR UR"/5n U R%XRSUSUR/5 ['SUR S-5HYn U R)U 5n UR*cM#U R-UR*5 U R/UR*5 M[ UGHnUR1S5nUR2R5U5n[7UR955nUR:UR<UR--UlUR@URBUR--Ul"[UR:U5n[UR>U5nURS:a-[UR@U5n[URDU5nO,[UR@U5n[URDU5n[UUR:- UR- S-5n[UUR@- UR- S-5n[UR<[UUR:- UR- S-55U- n[URB[UUR@- UR- S-55U- nUS::dUS::aGM[UU- UR- S-5n[UU- UR- S-5n[U R<[UU- UR- S-55U- n[U RB[UU- UR- S-55U- nUS::dUS::aGMUS:deUS:deUS:deUS:de['SUR S-5Hn UR)U 5nU R)U 5n URFbU RIURF5 U RKURLU S- 5 UROUUUUUUUR"5nUc[Q[RRT"55 U RWUUUUU5 M GM U $)z Find a gdal dataset representing a subset of a mosaic, based on a bounding box. Might overlap multiple tiles of the mosaic returns GDALDataset or None Nrr rarg?TEMP),rfrgrhSetSpatialFilterRectrirGetGeometryRef GetEnvelopeminmaxSetSpatialFilterrSrrrcCreaterlroSetGeoTransformryrmrsFillSetNoDataValuerjrerErrtr RasterXSizerr RasterYSizerrwSetRasterColorTableSetRasterColorInterpretationrx ReadRasterr<rGetLastErrorMsg WriteRaster)rminxminymaxxmaxyfeaturesenveloper featureEnv resultSizeX resultSizeYresultDSbandNrt_band featureNamesourceDSrtgw_ulxtgw_lrxtgw_ulytgw_lrysw_xoffsw_yoffsw_xsizesw_ysizetw_xofftw_yofftw_xsizetw_ysizes_banddatas r getDataSetmosaic_info.getDataSets $$&335 $$&;;DS))224CCEG OOG $"113??A$335AAC  1 x{3 1 x{3 1 x{3 1 x{3 $   qk " qk " qk " qk "   $$&77=4;$++5;< 4;$++5;< ??)) Kdjj$.."    $ Qa!MNAtzzA~.F++F3F{{& DKK(%%dkk2 / G!**1-Kzz~~k2H*8+C+C+EFCgg 4 4szz AACGgg 4 4szz AACG#''4(G#''4(G{{Qcggt,cggt,cggt,cggt,7SWW, :S@AG7SWW, :S@AGH((#w/@CJJ.NQT.T*UV  H((#w/@CJJ.NQT.T*UV 1}A 7T>T[[83>?G7T>T[[83>?GH((#w~.Ls.R*ST  H((#w~.Ls.R*ST 1}A a< <a< <a< <a< <4::>2!//7!//777&..tww733DGGFQJ4GH((NN<$..01""7GXxN'3g PrcAgrr-)rmemDSs r closeDataSetmosaic_info.closeDataSetAs rcN[SUR-5 [SURURUR4-5 [SUR UR 4-5 [SURURURUR4-5 g)Nz Filename: zFile Size: %dx%dx%dzPixel Size: %f x %fzUL:(%f,%f) LR:(%f,%f)) r<rdrVrWrlrrrrrrr6s rrZmosaic_info.reportEsv lT]]*+ #tzz4::tzz&JJK #t{{DKK&@@A '488TXXtxx*RRSr)rcrorlrerxrwrdrrrsrfrqrrrrrVrWN) r'r(r)r*r+rrJrrrZr,r-rrr]r]s#I0Jd CJTrr]c UR(a [S[UR5-SS9 [ URSUR SUR 5nURH|n[R"U5nUc g[UR55nURURUR5n[UR XUSUS5 AM~ UR(a [S5 U$)Nz*Building internal Index for %d tile(s) ... )end TileIndexrr finished)Verboser<r@NamescreateTileIndexTileIndexFieldNameTileIndexDriverTyprr;rrtr%rr addFeature)grf inputTilerrpointss rgetTileIndexFromFilesrLsyy :S\IsS$ ; 4 4dA>A-. /F %2>>A-. /F 77199  F c&k)Fzbmmbjj&@AGzbllRZZ&?@G]]FLLE!!&qVQG&qVD)7*"))+W,88@x(@L 'E8L 7771999Q U5\12-2 "  amm$Q*Q__.>q.ABCAtczH,HAx(((4 7rc[US5nUR5R5 UR5R5nUcGO[R R URS55nUR(aV[R R[R RURS555nUSS-U-nURU5 UR5nUR5nUH2n URUR5 URSU -5 M4 URS5 GM5UR5 g)Nr#rr rz%f )openrgrhrirrrrjrrrwriterr CsvDelimiterclose) rrrcsvfilerrrgeomcoordscoords rrrs8S!G NN!!# .."113 ? 77##G$4$4Q$78   bggoog.>.>q.ABCAtczH,H h%%'!!#E MM!.. ) MM$, '  d " MMOrc[RS:wa UR[RS5nURS5nUS:aSU-$USUS-S-XS-S-$)Nrrtmp_r )rrreplacerfind)rdposs r_createTempFileNamersg vv}##BFFC0 .. C Qw  Aa 6 )H1WY,? ??rc[RRU5(a[R"S5 UR R U5 [R"5 [RRU5(a[R"U5 UR RX!5 g)NCPLQuietErrorHandler) rrrrPushErrorHandlerDriverDeletePopErrorHandlerremoveRename)roldNamenewNames r_renameDatasetrso ww~~g 45    77>>' " " IIg HHOOG%rc [U5n URS-n URS-n [URX*--U SUR X;--SU /5n U(a0U R XE5n [URXvU SU S5 gURU RU R X\R--U RXLR--U R 5nUcgU R XE5n [URXvU SU S5 URc URnO URnU RU RSU R SU R/nURnURc)URRXUUXR 5nOURRSXEUU5nUc[#S[$R&S9 gUR)U5 UR+UR,5 [/SUS-5H]nUR1U5nUR2bUR5UR25 UR7UR8US- 5 M_ UR:b][/SUS-5HJnUR1U5nUR=UR:5 UR?UR:5 ML [@RB"UUSSURD5nUS:wa[#SU U4-[$R&S9 gURGU5 URcURI5 O:URRKU USUR 5nURI5 SnSn[MX U5 URN(aH[#US -[QU5-S -[QU5-S -[QU5-S -[QU5-5 gg) z Create an individual tile for the pyramids. The levelMosaicInfo object contains data about the mosaic at a given pyramid level. Returns None if successful and return 1 if there is an error rarr N'Creation failed, terminating gdal_tile.r9z$Reprojection failed for %s, error %d : |-->-))rrrrrrr%rrrBandTyperorl MemDriverrr CreateOptionsr<r>r?r SetProjectionrqryrmrwrrrxrsrrrReprojectImageResamplingMethodr FlushCache CreateCopyrrr)rlevelMosaicInforrrr tileNamerr temp_tilenamesxsyrrs_fhbtrrlt_fhbandrrestt_fhs rcreatePyramidTiler:s(1M   ! #B   ! #B "   ', .    ', .    Cu-1''&)VAYO  % % 6JJ..%**:L0Lcgg D | ]]5 )Fq##UfQiKzz  & & ZZGGSZZCGGQ CL  ! !E {{xx}VUBX{{!!"eUB? | 7cjjI&112a###D)    )  & &'9'9 :++O,>,>tax,HI $ )!UQY'D''-F KK.. /  ! !/"8"8 9(   dD$a6H6H IC ax 2mS5I I   & {{ ##M4AOOL  D1X.yy  'l  'l      %j   &k  rc [U5n URc URn O URn [URUR SUR SUR/5n U RX+R --U R SU R X;R--SU R/n U(a;[U 5n U RXE5n[URXvUSUS5 gURU RX+R --U R X;R--X[R--U RX+R --XKR --U R X;R--5nUcg[U 5n U RXE5n[URXvUSUS5 URnURc)URRXUUXR 5nOURRSXEUU 5nUc[#S[$R&S9 gUR)U 5 UR*b)UR-UR*R/55 [1UR2U5n[1UR4U5n[7SUS-5HnUR9U5nUR9U5nUR:bUR=UR:5 UR>b6URAUR>5 URCUR>5 URESSUUUUURF5nURISSUUUUUURF5 M URKU5 URcURM5 O:URROU USUR 5nURM5 SnSn[QX U5 URR(aH[#US-[UU5-S-[UU5-S -[UU5-S -[UU5-5 gg) a* Create add a vector feature representing the tile to the index, then recreate the Args: g (RetileGlobals): object with global script variables minfo (mosaic_info): mosiac_info object offsetX (int): The offset in the X direction. offsetY (int): The offset in the Y direction. width (int): The width of the tile. height (int): The height of the tile. tilename (str): The name of the tile. OGRDS (DataSource): The OGR DataSource object containing the tile index feature_only (bool): Whether to only generate features. Nrr r!r"r9r#r$r%r&)+rr'rorrrrrr%rrrrlr(rrr)r<r>r?rrr* ExportToWktrrrryrmrwrrsrrrrnrrr-r.rrr)rrrrrr rrrr1r5rrdec2rr4rlr6readXreadYr7rrrr9s rrrds$(1Mzz __ ZZ " ELL!UYY5<<@ C 'JJ&&   'JJ&&  L' 5.1''&)VAYO    'JJ&& 'JJ&&**)<< 'JJ&&);; 'JJ&&  D | $L 1D ^^E *Fq##UfQiK KKE{{xx}VUBX{{!!"eUB? | 7cjjI&|| 1<<3356   % (E   & )Ea###D)##D) 88   & &uxx 0 << # KK %  ! !%,, /  AueUE6??S1atUE6??S$ t{{ ##M4AOOL  D1X.yy  'l  'l      %j   &k  rc[R"U5nUc[S[RS9 gUR U5nUb6UR 5 URU5 U(a[SU-5 URU5nUc[SU-[RS9 gURSU[R5nUc[S[RS9 g[R"U[R5nUc[SU-[RS9 gURS 5 URU5S :wa[S U-[RS9 gU$) NzESRI Shapefile driver not foundr9r ztruncating index zCould not open datasource indexzCould not create LayerzCould not create FieldDefn for rzCould not create Field for )rrbr<r>r?r;DestroyDeleteDataSourceCreateDataSource CreateLayer wkbPolygon FieldDefn OFTStringSetWidth CreateField) rdsName fieldNamesrs driverName OGRDriver OGRDataSourceOGRLayer OGRFieldDefns rrrs5##J/I /cjjANN6*M ""6*  %. /..v6M *V3#**E((#s~~FH &SZZ8==CMM:L /);#**M#L)Q. +i7cjjI rc UR5n[R"UR55nUc[ S[ R S9 gURX5 SUSUSUSUSUSUSUSUSUSUS4 -n[R"XuR55nUc[ S[ R S9 gURU5 URU5 UR5 g) NzCould not create Featurer9r z*POLYGON ((%f %f,%f %f,%f %f,%f %f,%f %f ))rrarzCould not create Geometry) rgrFeature GetLayerDefnr<r>r?rCreateGeometryFromWktrSetGeometryDirectlyrrC) rrQlocationr!r"rR OGRFeaturewkt OGRGeometrys rrrs%%'HX2245J (szz:*5 6 a a a a a a a a a a : C++C1G1G1IJK ) ;"";/ :&rc$UR5 gr)rC)rQs rrrsrcUn[SURS-5HgnSUl[URU5n[ [ URS- 5[ URS- 5UUU5n [XX5nMi g)Nr r ra) ryLevelsrr]rdrMrSrVrWbuildPyramidLevel) rrcreatedTileIndexDSrPrQrUr~rr/levelOutputTileInfos r buildPyramidrcs Gq!((Q,' %enng>' %%) * %%) *     $A8KS(rc [[SURS-55n[[SURS-55n[ UR S[ U5-URURUR5nUHnUHnUS- URUR- -n US- URUR- -n URn URn X-UR:aURU - n X-UR:aURU - n [XX(Xs5n UR =(a ["R$R'U 5n[)UUU U U U U UU5 M M UR*b$[-X5UR*-n[/XU5 UR0b$[-X5UR0-n[3XU5 U$)zZ Build the pyramids at level N and returns an OGR dataset of the tile index at that level r TileResult_)rryrTrRrrrrrrrQrUrPrr rrrrrr:rrrrr)rr/rbrrrrrrrrr rrrrrs rr`r`.s%.::Q>? @F %.::Q>? @F  E "     EFz#..1D1L1LLGz#--0C0K0KKG)33F'11E!4!:!::+11G;"5"<"<<,33g="$7H88@x(@L  )B " *Q__< Ai0}} q(1==81W- LrcURnURU:a URn[[U55n[R R [R RUR55nUSSS:XaUSS[US5US4nS[U5-S-U4-n S[U5-S-U4-n UR(a[X5[U5-[R-US-S-U -S-U -n URU:ajX@l [R R[X5[U5-5(d+[R"[X5[U5-5 O[X5US-S-U -S-U -n URc XS- n U $U SUR-- n U $)z creates the tile file name r@r z%0i_.)rRrTr@rrrsplitextrrdrrrrrmkdir Extension) rrrrrrmaxim countDigitsparts yIndex_str xIndex_strfrmts rrris NNE ~~c%j/K GG  RWW--enn= >E Qx{cq!c%(m,eAh7[))C/F9>,q"83v;"FGGa/#f+=>A%a036CcIJV{{ a K akk!! Krc[S5 [R"5n[U5H.n[R"U5n[UR 5 M0 g)NzValid formats:r )r<rGetDriverCountry GetDriver ShortName)countrAdrivers r UsageFormatrzsH     !Eu& f rcVU(a[RO[Rn[SUS9 [SUS9 [SUS9 [SUS9 [SUS9 [SUS9 [SUS9 [S US9 [S US9 [S US9 [S US9 [S US9 U(aS$S$)Nz/Usage: gdal_retile.py [--help] [--help-general]r9z< [-v] [-q] [-co =]... [-of ]z( [-ps ]z! [-overlap ]zC [-ot {Byte|Int8|Int16|UInt16|UInt32|Int32|Float32|Float64|z0 CInt16|CInt32|CFloat32|CFloat64}]zB [-tileIndex [-tileIndexField ]]z1 [-csv [-csvDelim ]]zC [-s_srs ] [-pyramidOnly] -levels z6 [-r {near|bilinear|cubic|cubicspline|lanczos}]z% [-useDirForEachRow] [-resume]zA -targetDir []...rar)r>r?stdoutr<)isErrorfs rUsagers 3::A ;!D HqQ 41= -A6 OVWX <1E NUVW =AF OVWX BK 1: MTUV1QrcUc [5nUc[Rn[R"U5nUcgSnU[ U5:GaX#nUS:Xa [ SS9$US:XdUS:XaUS- nX#UlGOpUS:Xa`US- n[R"X#5Ul UR[R:Xa[S X#-[RS 9 gGO US :Xa$US- nURRX#5 GOUS :Xa S UlGOUS:Xa S UlGOUS:XaUS- nX#Ul[$R&R)UR"5(d)[SUR"-S-[RS 9 gUR"[ UR"5S- S[$R*:wa"UR"[$R*-UlGOUS:Xa0US- n[-X#5UlUS- n[-X#5UlGOUS:XaUS- n[-X#5UlGOUS:XaUS- nX#nUS:Xa[R4UlGO{US:Xa[R8UlGO^US:Xa[R:UlGOAUS:Xa[R<UlGO$US:Xa[R>UlGO[SU-5 gUS:XaBUS- n[-X#5Ul UR@S:a[SUR@-5 gGOUS:XaSUS- n[BRD"5Ul#URFRIX#5S:wa[SX#-5 gGOWUS:Xa S Ul%GOHUS :XaWUS- nX#Ul&[$R&ROURL5nUS(dU=RLS!- sl&OUS":XaUS- nX#Ul(OUS#:XaWUS- nX#Ul)[$R&ROURR5nUS(dU=RRS$- sl)OyUS%:XaUS- nX#Ul*OdUS&:XaS Ul+OVUS':XaS Ul,OHUSSS(:Xa$[S)U-[RS 9 [ S S9$URZRU5 US- nU[ U5:aGMURZ(d[S*5 [ S S9$UR.S:XdUR0S:Xa2[S+UR.UR04-[RS 9 gUR.UR2- S::dUR0UR2- S::a[S,[RS 9 gUR"c![S-[RS 9 [ S S9$URV(atURJ(dcUR"[]S5-[$R*-n[$R&R)U5(d[$R^"U5 UR@S:aSn[aXR@S-5Hn UR"[]U 5-[$R*-n[$R&R)U5(aMR[$R^"U5 [$R&R)U5(d[S.U-[RS 9 gUR(dM[S/U-[RS 9 M [Rb"UR5Ul2URdc/[S0UR-[RS 9 [g5$URdRi5n U Rk[Rl5Ul7S1U ;a[Rb"S25Ul8[sU5n U c[S3[RS 9 g[uURZSU 5n [wU RxU RzUR.UR0UR25n URFczU R|(ai[BRD"5Ul#URFRIU R|5S:wa&[S4U R|-[RS 9 gUR(a U R5 U R5 URJ(d[XU 5nU R5 OU nUR@S:a,[XXR.UR0UR25 UR(a [S55 g)6Nrr z--helpF)r}z-ofz-fz-otzUnknown GDAL data type: %sr9z-coz-vTz-qz -targetDirz TargetDir z does not existz-psz-overlapz-rnearbilinearcubic cubicsplinelanczoszUnknown resampling method: %sz-levelszInvalid number of levels : %dz-s_srszinvalid -s_srs: z -pyramidOnlyz -tileIndexz.shpz-tileIndexFieldz-csvz.csvz -csvDelimz-useDirForEachRowz-resumer&zUnrecognized command option: %szNo input files selected.zInvalid tile dimension %d,%dz'Overlap too big w.r.t tile height/widthz&Missing Directory for Tiles -targetDirzCannot create level dir: %szCreated level dir: %sz4Format driver %s not found, pick a supported driver. DCAP_CREATEr`zError building tile indexzinvalid projection FINISHED)C RetileGlobalsr>argvrGeneralCmdLineProcessorr@rFormatGetDataTypeByNamer' GDT_Unknownr<r?r)rrrrrrrrrS TileWidth TileHeightOverlapGRA_NearestNeighbourr, GRA_Bilinear GRA_CubicGRA_CubicSpline GRA_Lanczosr_rSpatialReferencerSetFromUserInputrrrkrrr rrrrrlryrbrrz GetMetadatarE DMD_EXTENSIONrmr(rr]rMrVrWrqrZrrCrc)argsrrrhargResamplingMethodStringrpleveldir startIndx levelIndxDriverMD tileIndexDSrrdsCreatedTileIndexs rmainrsy O |xx  ' ' -D | A c$i-g (?' ' E\SD[ FAwAH E\ FA//8AJzzT---2TW<3::N.E\ FA OO " "47 + D[AI D[AG L FA'AK77>>!++..lQ[[03DD3::V{{3q{{+a/12bff<kkBFF2 E\ FAdg,AK FAtw%>"':5%)%6%6"'72%)^^"'=8%)%9%9"'94%)%5%5"58NNO I  FA47|AHxx!|5@AH_ FA//1AL||,,TW5:(4723;N " AM L FA"gAOGG$$Q__5E86) % % FA#'7A F] FA GAMGG$$Q]]3E8 ' K  FA!WAN ' '!%A  I AH !W^ 3c9 K& & GGNN3  QO c$i-R 77 ()T""{{a1<<1, *akk1<<-H H {{QYY!#q||aii'?1'D 7cjjI{{ 6SZZHT"" !--;;Q'"&&0ww~~h'' HHX xx!| y((Q,7I{{S^3bff>(++3h>SZZPyyy-8szzJ8##AHH-AHxx BQXX M }xx##%H,,t112AKH$**51 '*K ) ;  K 0E 5;; Q[[!,, RB|| 0 0++- << ( ()9)9 :a ? (5+;+;;#** Myy  ==&q4(xx!|Q1;; aiiXyy j rc"\rSrSr/SQrSrSrg)riu)rrr)rrrrrr'rrmr(rrrr rrrr,r_rrrrc~SUlSUl/Ul/UlSUlSUlSUlSUlSUlSUl SUl SUl SUl SUl SUlSUlSUlSUlSUl[&R(UlSUlSUlS UlSUlSUlg) zOnly used for unit testsFrBrGTiffNrYMemory;r )rrr)rrrrrr'rrmr(rrrr rrrrrr,r_rrrrr6s rrRetileGlobals.__init__s       ",!"* $ 9 9   % r)r'r)r rrrmrrr_r(rrrrr,rrrrrrrrrrN)r'r(r)r* __slots__rr,r-rrrrusI8rr__main__)r )NN)% __future__rrr>osgeorrrTermProgress_nocbrobjectrr/rMr]rrrrrrrr:rrrrrcr`rrzrrrr'exitrr-rrrs>&  ! !v>6<00:FT&FTR6- AH..@ &l ^q h!H> T8v&R "C L88v zHHT#((^r