!\iIBSrSSKJrJr SSKrSSKrSSKrSSKrSSKrSSK r SSK r SSK r SSK Jr SSKJrJrJrJr SSKJr SSK Jr SSKJr SS KJr \R6"\5rS rS rS r S r!Sr"Sr#\\\ \!\"\#S.r$SSSSSSS\RJSSSSSSS4Sjr&g)z5Copy valid pixels from input files to an output file.) ExitStackcontextmanagerN) Resampling) MergeErrorRasterioDeprecationWarning RasterioError WindowError) DatasetWriter)windows)Affine) subdividec [R"USS9n[R"X5S9 [R"X%US9 [R"XUSS9 g)z"Returns the first available pixel.booldtypeoutunsafewherecastingN)np empty_like logical_not logical_andcopyto merged_datanew_data merged_masknew_maskkwargsmasks Q/var/www/html/kml_chatgpt/mouzaenv/lib/python3.13/site-packages/rasterio/merge.py copy_firstr%s< ==F 3DNN8&NN;$/IIk4Bc [R"USS9n[R"X5S9 [R"XUSS9 g)z!Returns the last available pixel.rrrrrN)rrrrrs r$ copy_lastr(%s. ==F 3DNN8&IIk4Br&c *[R"USS9n[R"X#US9 [R"XUS9 [R"XXSS9 [R"X5S9 [R "X%US9 [R "XUSS9 g)z Returns the minimum value pixel.rrrrrrrrN)rr logical_orrminimumrrrs r$copy_minr-,f ==F 3DMM+T2NN4"JJ{+8TNN8&NN;$/IIk4Br&c *[R"USS9n[R"X#US9 [R"XUS9 [R"XXSS9 [R"X5S9 [R "X%US9 [R "XUSS9 g)z Returns the maximum value pixel.rrrrr*rN)rrr+rmaximumrrrs r$copy_maxr17r.r&c *[R"USS9n[R"X#US9 [R"XUS9 [R"XXSS9 [R"X5S9 [R "X%US9 [R "XUSS9 g)z$Returns the sum of all pixel values.rrrrr*rNrrr+raddrrrs r$copy_sumr5Bsf ==F 3DMM+T2NN4"FF;kxPNN8&NN;$/IIk4Br&c *[R"USS9n[R"X#US9 [R"XUS9 [R"XXSS9 [R"X5S9 [R "X%US9 [R "XUSS9 g)z"Returns the count of valid pixels.rrrrr*rNr3rs r$ copy_countr7Msf ==F 3DMM+T2NN4"FF;+8LNN8&NN;$/IIktX>r&)firstlastminmaxsumcountr8F@cUb"S[5 U [;a [U nOI[ U 5(aU nO6[ SR U [[R5555e[US[[R45(a[RnO[S5nUnSn[!5nU"US5nUR"nUR$nUR&nUR(SnUnUR*SnUc UR,nO#[U[.5(aUnO [1U5nUR3S5nSSS5 U(dWnU(a Uunnn n!GO /n"/n#[5U5Hun$n%U"U%5n&U&R6n'U (a[9WU&R&SS9nU'R:(d [=S 5eU'R>S:a [=S 5eU'R@S:a [=S 5eU&RBun(n)n*n+SSS5 U"REW(W*/5 U#REW)W+/5 M [9U"5[9U#5[GU"5[GU#54unnn n!U(dWnO<[U[HRJ5(aX"4nO[1U5S:Xa USUS4nU (a[LRN"UUS- 5US-n[LRP"U US- 5US-n [LRN"UUS- 5US-n[LRP"U!US- 5US-n![/[SU U- US- 55n,[/[SU!U- US- 55n-[TRV"UU!5[TRX"USUS*5-n.UbUnR]S U5 UbUn[ZR]S U5 Sn/WGb[^R`"W[^Rb5(a@[^Rd"U5n0U0R8Us=:*=(a U0RF:*Os n/O[fRh"U5(ap[^Rj"U5n0U0R8Us=:*=(a U0RF:*Os n/[^Rl"U5n1U/[^Rn"U1U5-n/OSn/U/(d[R"SUSUS35 WnO[ZR]S5 Sn[pRr"SSU,U-5n2Ub[U[t5(aUnOgWn3U3Rv"S%0U=(d 0D6 U.U3S'U-U3S'U,U3S'UU3S'WU3S'UbUU3S'[R"US40U3D6nURyU5 U S-[^Rz"W5R|U-- n4U,U--U4:aU2/n5O;[LRN"[LR~"U455n6[U2U6U65n5OU2/n5Sn7U5GH`n8[pRB"U8U.5unnn n![^R"UU8RU8R4WS9n9U/(aU9RU5 [pRB"U8U.5n:[pR6"U8U.5n;Sn<[5U5GHun=n%U"U%5n&WU&R$:wa[SU%35eU7"U&RBU:U;5n>[pR"/U>QU&R6P76n?[pR"/U>QU;P76n@U<"W@5n@U@R5unAnBU9SS2UAUB4nC[fR"U5(a[^R"WC5nDOQ[^R`"WCRz[^Rb5(d[^R"WCU5nDOWCU:HnDU&RWW@RU@R4USU?US!9nEU"WCUEWDUERU=U@RU@RS"9 SSS5 GM U(dGM0[pR"/U:QU.P76nFU<"UF5nFURU9UFS#9 GMc Uc3U (a[^RRW9USS$9n9W9U.4sSSS5 $W(aURSU5 UR5 SSS5 g![ a SnGNf=f!,(df  GN=f!,(df  GN"=f![ [4a. [ZRS U&U&RB5 SSS5 GMf=f!,(df  GM=f!,(df  g=f)&aCopy valid pixels from input files to an output file. All files must have the same number of bands, data type, and coordinate reference system. Rotated, flipped, or upside-down rasters cannot be merged. Input files are merged in their listed order using the reverse painter's algorithm (default) or another method. If the output file exists, its values will be overwritten by input values. 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 first input file. 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. nodata: float, optional nodata value to use in output file. If not set, uses the nodata value in the first input raster. masked: bool, optional. Default: False. If True, return a masked array. Note: nodata is always set in the case of file output. dtype: numpy.dtype or string dtype to use in outputfile. If not set, uses the dtype value in the first input raster. precision: int, optional This parameters is unused, deprecated in rasterio 1.3.0, and will be removed in version 2.0.0. indexes : list of ints or a single int, optional bands to read and merge 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`. method : str or callable pre-defined method: first: reverse painting last: paint valid new on top of existing min: pixel-wise min of existing and new max: pixel-wise max of existing and new or custom callable with signature: merged_data : array_like array to update with new_data new_data : array_like data to merge same shape as merged_data merged_mask, new_mask : array_like boolean masks where merged/new data pixels are invalid same shape as merged_data index: int index of the current dataset within the merged dataset collection roff: int row offset in base array coff: int column offset in base array 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 merge 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 ------ MergeError When sources cannot be merged due to incompatibility between them or limitations of the tool. NzLThe precision parameter is unused, deprecated, and will be removed in 2.0.0.z0Unknown method {}, must be one of {} or callablerc3 # Uv g!f=f7f)N)objs r$ nullcontextmerge..nullcontexts  s  c[U[R5(aU$[R"USS-USS--5$)NrrE) isinstancenumbersNumbermathsqrt)xs r$merge..!sE)!W^^<<+,+B!%1Q419qtqy+@!A+Br&)keyz2Rotated, non-rectilinear rasters cannot be merged.zGRasters with negative pixel width ("flipped" rasters) cannot be merged.zLRasters with negative pixel height ("upside down" rasters) cannot be merged.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 transformheightwidthr=rnodatawg.AcF[USUS5n[USUS5nX4:a[eURS:a0[USUS5n[USUS5nXV:a[eO/[USUS5n[USUS5nXe:a[eX5XF4$)zBased on gdal_merge.py.rrGrE)r;r: ValueErrore)bounds1bounds2rQint_wint_eint_sint_ns r$_intersect_bounds merge.._intersect_boundss GAJ/E GAJ/E~  {{QGAJ 3GAJ 3>$$"GAJ 3GAJ 3>$$- -r&rcH[R"URS-5n[R"URS-5n[R"URS-5n[R"UR S-5n[ R"X!XC5$)a Equivalent to rounding both offsets and lengths. This method computes offsets, width, and height that are useful for compositing arrays into larger arrays and datasets without seams. It is used by Rasterio's merge tool and is based on the logic in gdal_merge.py. Returns ------- Window g?g?)rKfloorrow_offcol_offrRrSr Window)windowrdrerRrSs r$ win_alignmerge..win_alignsq**V^^c%9:**V^^c%9:FMMC$78 6<<##56~~gFFr&zCRS mismatch with source: z"Skipping source: src=%r, bounds=%r) out_shapeindexesmaskedrg resampling)indexroffcoff)rg)copyrA)Uwarningswarnr MERGE_METHODScallablerXformatlistkeysrHstrosPathLikerasterioopenrrprofilecrsres nodatavalsdtypesr=intlencolormap enumeraterQr:is_rectilinearrarYboundsextendr;rIrJrKrcceilroundr translationscaleloggerdebugr issubdtypeintegeriinfocmathisfinitefinfomin_scalar_typecan_castr rfr update enter_contextritemsizerLr zerosrRrSfillr from_boundsr infotoslicesisnanisclosereadr#rdrewritema masked_equalwrite_colormapclose)GsourcesrrrTr precisionrk output_countrmmethodtarget_aligned_pixels mem_limituse_highest_resrldst_pathdst_kwdsrdataset_openerrCdst exit_stackr8 first_profile first_crsbest_resfirst_nodataval nodatavaldt src_countfirst_colormapdst_wdst_sdst_edst_nxsysidatasetsrc src_transformleftbottomrighttop output_width output_heightoutput_transforminranger nodata_dt dout_window out_profile max_pixelschunksnr`chunkdest chunk_boundschunk_transformrhidxiboundsswcwrowscolsregion region_maskdatadwsG r$mergerbsl Z & v& &   > E E]//12   '!*sBKK011!     % C  GAJ '5!MMM IyyH#..q1O'IaB!KK GS))# L  &!&!2!(($L )/ &E5%BB'0 7#G,$'MMM&#&$GG!B$)77(P%*(e%*(j03zz,D&%7-: 4-( 63-(?1B*-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>%1ABL%//7GHO G$!*' 2 W#G,!CGG++.H ,RSS !"3JJ o#%00I'I3==I$00K'K?K#2B!#JD$!!T4-0F{{9--&(hhv&6 ]]6<<DD&(jj&C &, &9 88#,bii"B '#!#- $D# !ZZZZK-,!3`s((I,I8HIr] $r *_b ;uu))$ )F))I L""1n5 IIKQ $ &!% &#( '@-,p' 4! @#szz!!-,!-,[ s f1A8d9d&%>f1#B e 0Sf1f 'Ae5Crs;0 % &%&   8 $CCCCC?           !!    !_r&