nj=SSKrSSKrSSKJr SSKJrJrJrJr SSK r SSK J r SSK Jr SSKJrJr SSKJr SSKJrJr SS KJrJrJr SS \S \S \\\S \\\\4S\\\S\S\S\\\\4S\\\\4S\S\S\S\\4Sjjr"SS\5r\R@4Sjr!\"S:Xa#\RF"\!"\R@55 gg)N)Number)OptionalSequenceTupleUnion)gdal) PathLikeOrStr)GDALArgumentParser GDALScript))GDALTypeCodeAndNumericTypeCodeFromDataSet)OptionalProgressCallbackget_progress_callback)PathOrDS get_bandsopen_dssrcfiledstfilesrcwinskip band_numsdelim skip_nodata src_nodata dst_nodatareturn_np_arrayspre_allocate_np_arraysprogress_callbackreturnc v Sn [U 5n [U5n U c[SUS35e[X5n[ U5nU R 5nUcSSU R U R4n[U 5unnUb [US5nOU (aSnO[RnU(aU[R:XdU[R:Xa$SU-[ U5-RU5S-nO#SU-[ U5-RU5S-n[!US5S :a_[!US 5S :aM[!U R US -5S :a.[!U RUS -5S :aS U-S -U-S-nOSU-S-U-S-n[#U[$5(aU/U-nOUc['SU55nSU;aSnUb[(R*"UUS9n[#U[$5(aU/U-nUb SU;dUcSnUb[(R*"UUS9nU=(a USLnU(+=(a USLnU=(d Un[#U[,5(aUunnOU=nnUunnnn[ U5nUU- U-nUU- U-n UU -n!Sn"Sn#Sn$U (aNU (aU!OSn%[(R."U%5n&[(R."U%5n'[(R."U%U4US9n(Sn)[1UUU-U5GHn*U (aUOSn%[(R."U%U4US9n+[3U5HAun,n-U-R5UU*US 5n.U (a U.SU+U,'M+[(R6"U+U.SS9n+MC [1SUU5GHXn/U"S - n"U (a#U"U!- n0[9U0U$-5n1U1U#:a U1n#U "U05 U+SS2U/4n2U(a.[(R:"UU25(aU(aMhU(aUn2U/U-n3USU3S-US --U*S-US--n4US U3S-US--U*S-US --n5U(a9W[=U25-n6W[?U45[?U55U64-n7URAU75 U (a\U (aU4W&U)'U5W'U)'U2W(U)'OE[(R6"W&U45n&[(R6"W'U55n'[(R6"W(U2/SS9n(U)S - n)GM[ GM U (aCU(aSO U(aUOUn8U)U":waW&SU)n&W'SU)n'W(SU)2SS24n(W&W'W(RC5U84n U $)a translates a raster file (or dataset) into xyz format skip - how many rows/cols to skip each iteration srcwin (xoff, yoff, xsize, ysize) - Selects a subwindow from the source image for copying based on pixel/line location. band_nums - selected input bands to process, None to process all. delim - the delimiter to use between values in a line skip_nodata - Exclude the output lines with nodata value (as determined by srcnodata) src_nodata - The nodata value of the dataset (for skipping or replacing) default (`None`) - Use the dataset NoDataValue; `Sequence`/`Number` - use the given nodata value (per band or per dataset). dst_nodata - Replace source nodata with a given nodata. Has an effect only if not setting `-skipnodata` default(`None`) - use srcnodata, no replacement; `Sequence`/`Number` - replace the `srcnodata` with the given nodata value (per band or per dataset). srcfile - The source dataset filename or dataset object dstfile - The output dataset filename; for dstfile=None - if return_np_arrays=False then output will be printed to stdout return_np_arrays - return numpy arrays of the result, otherwise returns None pre_allocate_np_arrays - pre-allocated result arrays. Should be faster unless skip_nodata and the input is very sparse thus most data points will be skipped. progress_callback - progress callback function. use None for quiet or Ellipsis for using the default callback NzCould not open .rwtz%d z%gz%.10gz%sz%.3fc3@# UHoR5v M g7fN)GetNoDataValue).0bands P/var/www/html/pdf-tiff/venv/lib/python3.13/site-packages/osgeo_utils/gdal2xyz.py gdal2xyz..sBED--//Es)dtyped)axisg?)"rr ExceptionrlenGetGeoTransform RasterXSize RasterYSizer opensysstdoutr GDT_Int32 GDT_UInt32rstripabs isinstancerlistnpasarrayremptyrange enumerate ReadAsArrayappendint array_equaltuplefloatwrite transpose)9rrrrrrrrrrrrresultdsbands band_countgtdtnp_dtdst_fh band_formatfrmtreplace_nodataprocess_nodatax_skipy_skipx_offy_offx_sizey_size bands_countnXBlocksnYBlocks progress_end progress_curr progress_prevprogress_partssize all_geo_x all_geo_yall_dataidxydatai_bndr+ band_datax_i progress_fracprogressx_i_dataxgeo_xgeo_yband_strlinenodatas9 r,gdal2xyzr|2sTHF-.?@  B z/'!455 b $EUJ   B~Q79"=IBgt$   2#8 5LCJ6>>uELK 5LCJ6>>uELK 1J BqE S BNNRU*+c1BNNRU*+c1U?W,u4t;DE>F*U2T9D*f%% \J.  BEBB  z ZZ %8 *f%% \J.  2 8J ZZ %8 ::T#9K%oCJd,BN 2NN$!!#) E5&&e*K6)H6)Hh&LMMN5|1HHTN HHTN 88T:.e< C 5%&.& 14{!xxve4$U+KE4((61=I%'lU yyyq9 ,FF+C Q M - < }~=>m+$,M%m4AsF|H"..X"F"F#)He AqEQW1--SBqE0AAEqEQW1--SBqE0AAE&x8uU|U5\8DD T")%*IcN%*IcN$,HSM " )U ;I " )U ;I!yyH:AFH 1HCG,2`$*J - !$3I!$3Ia(HIx'9'9';VC Mc^^\rSrSrU4SjrS\4SjrS\4U4SjjrS\4Sjr Sr Sr U=r $) GDAL2XYZcf>[TU]5 SUl[R"S5Ulg)Nz(Translates a raster file into xyz formata\ The gdal2xyz utility can be used to translate a raster file into xyz format. It can be used as an alternative to gdal_translate of=xyz, But supporting other options, for example: * Select more then one band; * Skip or replace nodata value; * Return the output as numpy arrays.)super__init__titletextwrapdedent description)self __class__s r,rGDAL2XYZ.__init__s+ ? #?? 4 r}rc URnURSSSSSS9 URSSS [S S S 9 URS SSSS[SSS9 URSSSSSS9 URSSSSSSS 9 URS!S"S#S$SS%S9 URS&S'S([S)S*S+9 URS,S-[S)S.S+9 URS/S0[S1S29 URS3S4[S5S29 U$)6Nz-skipr store_truer%z-How many rows/cols to skip in each iteration.)destactiondefaulthelpz-srcwin)xoffyoffxsizeysizerr4zRSelects a subwindow from the source image for copying based on pixel/line location)metavarrtypenargsr-b-band--bandrr++a Select bands from the input spectral bands for output. Bands are numbered from 1 in the order spectral bands are specified. Multiple -b switches may be used. When no -b switch is used, the first band will be used.In order to use all input bands set -allbands or -b 0..)rrrrrz -allbandsz --allbandsallbandszSelect all input bands.)rrrz-csvr, store_constz*Use comma instead of space as a delimiter.)rconstrrrz -skipnodataz --skipnodataz -skip_nodatarzHExclude the output lines with nodata value (as determined by srcnodata). -srcnodata -nodatavaluer*zThe nodata value of the dataset (for skipping or replacing) Default (None) - Use the dataset nodata value; Sequence/Number - Use the given nodata value (per band or per dataset).)rrrr -dstnodatarzReplace source nodata with a given nodata. Has an effect only if not setting -skipnodata. Default(None) - Use srcnodata, no replacement; Sequence/Number - Replace the srcnodata with the given nodata value (per band or per dataset).r src_datasetzsThe source dataset name. It can be either file name, URL of data source or subdataset name for multi-dataset files.)rrrr dst_datasetzThe destination file name.)parser add_argumentrMrJstr)rargvrs r, get_parserGDAL2XYZ.get_parsers @    6e      F     *     =      [     V    )    !M    !-   r}c@>SnSn[U5n/nSn/n/n/n Xd:GaXn U S;aXdS- :Xa[SU S35eUS- nU"X5(d[SU S 35eXd:aFU"X5(a7UR[X55 US- nXd:aU"X5(aM7GO.U S ;aXdS- :Xa[SU S35eUS- nU"X5(d[SU S 35eXd:aFU"X5(a7UR[ X55 US- nXd:aU"X5(aM7OU S ;aXdS- :Xa[SU S35eUS- nU"X5(d[SU S 35eXd:aFU"X5(a7U R[ X55 US- nXd:aU"X5(aM7OUS- nURU 5 Xd:aGM[ [ U]U5n U(aX{S 'U(aXS'U (aXS'U $)Nc<[U5 g![a gf=fNTF)rJ ValueErrorss r,checkInt GDAL2XYZ.parse..checkInt\s# A    c<[U5 g![a gf=fr)rMrrs r, checkFloat"GDAL2XYZ.parse..checkFloatcs# a  rr)rrrr%zMissing argument following z: zArgument following z should be a integer)rrz should be a floatrrrr)r6r5rIrJrMrrparse) rrrrcountnew_argvirrrargkwargsrs r,rGDAL2XYZ.parse[sW  D     i'C// >#&A#b$IJJQ((#&9#>R$STTiHTW$5$5$$S\2FAiHTW$5$566 >#&A#b$IJJQ!$'**#&9#>P$QRRiJtw$7$7%%eDGn5FAiJtw$7$7& >#&A#b$IJJQ!$'**#&9#>P$QRRiJtw$7$7%%eDGn5FAiJtw$7$7Q$?iBx,X6 "+;  #-< #-<  r}czURS5(aSUS'OURS5(dSUS'US U$)Nrrr%)getrrs r,augment_kwargsGDAL2XYZ.augment_kwargssA ::j ! !"&F; K(("#F;  :  r}c [S0UD6$)N)r|rs r,doit GDAL2XYZ.doits!&!!r})rr) __name__ __module__ __qualname____firstlineno__rr rdictrrr__static_attributes__ __classcell__)rs@r,rrs=  b"4bH@T@D""r}rc4[5RU5$r()rmain)rs r,rrs :??4  r}__main__) NNr%NrFNNFT.)$r;rnumbersrtypingrrrrnumpyrCosgeorosgeo_utils.auxiliary.baser #osgeo_utils.auxiliary.gdal_argparser r osgeo_utils.auxiliary.numpy_utilr osgeo_utils.auxiliary.progressr rosgeo_utils.auxiliary.utilrrrrJrboolr|rrrrexitrr}r,rsk@ 334NVDC "&*&')-4848"#'25t t t Xc] #t Xc]" # t  & t  ttx/01tx/01tt!t0te_tn}"z}"@hh! zHHT#((^r}