!\i3 SrSSKJr SSKJrJr SSKrSSKrSSKrSSK r SSK r SSK r SSK r SSKrSSKJr SSKJr SSKJrJr SSKJr SS KJr SS KJr SS KJr \R:"\5rSSSSSS\R@S S S S SS4 Sjr!g)zRaster stacking tool.)Iterable) ExitStackcontextmanagerN)disjoint_bounds) Resampling) RasterioError StackError) DatasetWriter)windows)Affine) subdivideF@c [US[[R45(a[R nO[ S5nUnSn[5nU"US5nURnURnURnURSnUnURSnUcUVs/sHnSPM nnURS5nSSS5 /n/nSn[[!X55GH ununnU"U5n U R"n!UcUU R$- nO)[U[&5(aUS- nOU[)U5- nU (a[+WU RSS9nU!R,(d [/S5eU!R0S:a [/S5eU!R2S:a [/S 5eU R4un"n#n$n%SSS5 UR7W"W$/5 UR7W#W%/5 GM U(aUun&n'n(n)O/[+U5[+U5[9U5[9U54un&n'n(n)U(dWnO<[U[:R<5(aX"4nO[)U5S:Xa USUS4nU(a[>R@"U&US- 5US-n&[>RB"U(US- 5US-n([>R@"U'US- 5US-n'[>RB"U)US- 5US-n)['[EU(U&- US- 55n*['[EU)U'- US- 55n+[FRH"U&U)5[FRJ"USUS*5-n,UbUn[LROS U5 UbUn[LROS U5 S n-WGb[PRR"W[PRT5(a@[PRV"U5n.U.R*Us=:*=(a U.R8:*Os n-O[XRZ"U5(ap[PR\"U5n.U.R*Us=:*=(a U.R8:*Os n-[PR^"U5n/U-[PR`"U/U5-n-OS n-U-(d[bRd"SUSUS35 WnO[LROS5 Sn[fRh"SSU*U+5n0U b[U [j5(aU nOgWn1U1Rl"S%0U =(d 0D6 U,U1S'U+U1S'U*U1S'UU1S'WU1S'UbUU1S'[R "U S40U1D6nURoU5 U S-[PRp"W5RrU-- n2U*U+-U2:aU0/n3O;[>R@"[>Rt"U255n4[wU0U4U45n3OU0/n3[LROSU35 U3GHn5[fR4"U5U,5un&n'n(n)[PRx"UU5RzU5R|4WS9n6U-(aU6RU5 Sn7[[!X55GHSun8unnU"U5n [U&U'U(U)4U R45(a&[LROSU U&U'U(U)45 SSS5 MXWU R:wa[SU35e[fR"U&U'U(U)U R"5RES5n9Uc U RnO[U[&5(aU/n[)U5U5RzU5R|4n:U RU:U9S S UUS9n;SSS5 [U[&5(aU6U7SS2SS24n[PRRU>5n?[PR"US!S"9 GMV U(dGM[fR"U&U'U(U)U,5RES5n@URU6U@S#9 GM Uc3U (a[PRRW6US S$9n6W6U,4sSSS5 $W(aURSU5 UR5 SSS5 gs snf![a SnGNf=f!,(df  GN=f!,(df  GN=f!,(df  GNE=f!,(df  g=f)&am Copy valid pixels from input files to an output file. All files must have the same data type, and coordinate reference system. Rotated, flipped, or upside-down rasters cannot be stacked. Geospatial bounds and resolution of a new output file in the units of the input file coordinate reference system may be provided and are otherwise taken from the source datasets. Parameters ---------- sources : list A sequence of dataset objects opened in 'r' mode or Path-like objects. bounds: tuple, optional Bounds of the output image (left, bottom, right, top). If not set, bounds are determined from bounds of input rasters. res: tuple, optional Output resolution in units of coordinate reference system. If not set, a source resolution will be used. If a single value is passed, output pixels will be square. use_highest_res: bool, optional. Default: False. If True, the highest resolution of all sources will be used. If False, the first source's resolution will be used. masked: bool, optional. Default: False. If True, return a masked array. Note: nodata is always set in the case of file output. nodata: float, optional nodata value to use in output file. If not set, uses the nodata value in the first input raster. dtype: numpy.dtype or string dtype to use in outputfile. If not set, uses the dtype value in the first input raster. indexes : list of ints or a single int, optional bands to read and stack. output_count: int, optional If using callable it may be useful to have additional bands in the output in addition to the indexes specified for read. resampling : Resampling, optional Resampling algorithm used when reading input files. Default: `Resampling.nearest`. target_aligned_pixels : bool, optional Whether to adjust output image bounds so that pixel coordinates are integer multiples of pixel size, matching the ``-tap`` options of GDAL utilities. Default: False. mem_limit : int, optional Process stack output in chunks of mem_limit MB in size. dst_path : str or PathLike, optional Path of output dataset. dst_kwds : dict, optional Dictionary of creation options and other parameters that will be overlaid on the profile of the output dataset. Returns ------- tuple Two elements: dest: numpy.ndarray Contents of all input rasters in single array out_transform: affine.Affine() Information for mapping pixel coordinates in `dest` to another coordinate system Raises ------ StackError When sources cannot be stacked due to incompatibility between them or limitations of the tool. rc3 # Uv g!f=f7f)N)objs Q/var/www/html/kml_chatgpt/mouzaenv/lib/python3.13/site-packages/rasterio/stack.py nullcontextstack..nullcontextus  s  Nc[U[R5(aU$[R"USS-USS--5$)Nrr) isinstancenumbersNumbermathsqrt)xs rstack..sC%a88'('>!YYqtqy1Q419'<='>)keyz3Rotated, non-rectilinear rasters cannot be stacked.zHRasters with negative pixel width ("flipped" rasters) cannot be stacked.zMRasters with negative pixel height ("upside down" rasters) cannot be stacked.z Set dtype: %szSet nodataval: %rFTz)Ignoring nodata value. The nodata value, z8, cannot safely be represented in the chosen data type, zs. Consider overriding it using the --nodata option for better results. Falling back to first source's nodata value.zSet nodataval to 0 transformheightwidthcountdtypenodatawg.Az Chunks=%r)r'z"Skipping source: src=%r, bounds=%rzCRS mismatch with source: ) out_shapewindow boundlessmaskedindexes resamplingrunsafe)casting)r,)copyr)PrstrosPathLikerasterioopenrrprofilecrsres nodatavalsdtypescolormap ValueError enumeratezipr#r&intlenminis_rectilinearr aeboundsextendmaxrrrfloorceilroundr translationscaleloggerdebugnp issubdtypeintegeriinfocmathisfinitefinfomin_scalar_typecan_castwarningswarnr Windowr update enter_contextr'itemsizerr zerosr$r%fillrr from_boundsr/readrisnaniscloseshapemagetmaskcopytowrite masked_equalwrite_colormapclose)AsourcesrHr;r(r'r/ output_countr0target_aligned_pixels mem_limituse_highest_resr.dst_pathdst_kwdsdataset_openerrdst exit_stackfirst first_profile first_crsbest_resfirst_nodataval nodatavaldtsfirst_colormapxsysidataset src_indexessrc src_transformleftbottomrighttopdst_wdst_sdst_edst_n output_width output_heightoutput_transforminrangeinfo nodata_dt dout_window out_profile max_pixelschunksnchunkdestdst_idxidx src_window temp_shapetemp_srcregion region_masktemp temp_mask dst_windowsA rstackrs n'!*sBKK011!     % C  GAJ '5!MMM IyyH#..q1O'IaB)01A41 &!&!2("  )23w3H)I %A%(C # & CII-L S11 A%L C $44L"" > H%33$M!??Q&$b!??Q&$g,/::(feSE)H IItUm $ IIvsm $M*JR )/ &E5%),R#b'3r7CG)K &E5%C W^^ , ,*C X]q63q6"C JJus1v~.Q7EIIec!fn-A6EJJus1v~.Q7EIIec!fn-A6E5%%-3q6!9:; E55=CF":;< !--eU;fll FSVG?    B LL" -  I LL,i 8  }}R,,xx|((i;;488;>>),,88B 9g% g g% g" g%% g3)"__doc__collections.abcr contextlibrrloggingr5rrVr[rnumpyrRr7rasterio.coordsrrasterio.enumsrrasterio.errorsrr rasterio.ior r rasterio.transformr rasterio.windowsr getLogger__name__rPnearestrrr!rrs$0 +%5%%&   8 $    !!   Pr!