(e=ddlZddlZddlmZddlmZmZmZmZddl Z ddl m Z ddl mZddlmZmZddlmZddlmZmZdd lmZmZmZ dd ed ed eeed eeeefdeeedededeeeefdeeeefdedededeefdZGddeZej@fdZ!e"dk(r#ejFe!ej@yy)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 ` d} t| } t|} | td|dt| |}t |}| j }|dd| j | jf}t| \}}| t|d}n| rd}ntj}|r|tjk(s|tjk(r$d|zt |zj|dz}n#d|zt |zj|dz}t!|dd kr\t!|d d krKt!| j |d zd kr-t!| j|d zd krd |zd z|zdz}nd|zdz|zdz}t#|t$r|g|z}n|t'd|D}d|vrd}|t)j*||}t#|t$r|g|z}|d|vs|d}|t)j*||}|xr|du}| xr|du}|xs|}t#|t,r|\}}n|x}}|\}}}}t |}||z |z}||z |z} || z}!d}"d}#d}$| rI| r|!nd}%t)j.|%}&t)j.|%}'t)j.|%|f|}(d})t1|||z|D]}*| r|nd}%t)j.|%|f|}+t3|D]<\},}-|-j5||*|d }.| r |.d|+|,<%t)j6|+|.d}+>t1d||D],}/|"d z }"| r"|"|!z }0t9|0|$z}1|1|#kDr |1}#| |0|+dd|/f}2|rt)j:||2r|rQ|r|}2|/|z}3|d|3dz|d zz|*dz|dzz}4|d |3dz|dzz|*dz|d zz}5|r9t=|2z}6t?|4t?|5|6fz}7|jA|7| rW| r|4&|)<|5'|)<|2(|)<nEt)j6&|4}&t)j6'|5}'t)j6(|2gd}(|)d z })/| r8|rdn|r|n|}8|)|"k7r&d|)}&'d|)}'(d|)ddf}(&'(jC|8f} | S)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<K|]}|jywN)GetNoDataValue).0bands 6/usr/lib/python3/dist-packages/osgeo_utils/gdal2xyz.py zgdal2xyz..sBD$--/Bs)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{2sYHF-.?@  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%&.& 1.4{!xxve4$U+ :KE4((61=I%'lU yyyq9  :FF+# C Q M - < }~=>m+$,M%m4AsF|H"..X"F#)He AqEQW1--SBqE0AAEqEQW1--SBqE0AAE&x8uU|U5\8DD T")%*IcN%*IcN$,HSM " )U ;I " )U ;I!yyH:AFH 1HCG# .`$*J - !$3I!$3Ia(HIx'9'9';VC McLeZdZfdZdefdZdeffd ZdefdZdZ xZ S)GDAL2XYZcdt|d|_tjd|_y)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,rzGDAL2XYZ.__init__s+ ? #?? 4 r|rc |j}|jddddd|jddd td d |jd ddddtdd|jddddd|jdddddd |jd!d"d#d$dd%|jd&d'd(td)d*+|jd,d-td)d.+|jd/d0td12|jd3d4td52|S)6Nz-skipr store_truer%z-How many rows/cols to skip in each iteration.)destactiondefaulthelpz-srcwin)xoffyoffxsizeysizerr3zRSelects 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_argumentrLrIstr)rargvrs r, get_parserzGDAL2XYZ.get_parsers @    6e      F     *     =      [     V    )    !M    !-   r|cd}d}t|}g}d}g}g}g} ||kr||} | dvr~||dz k(rtd| d|dz }|||std| d ||krT|||rH|jt|||dz }||kr |||r4n| d vrz||dz k(rtd| d|dz }|||std| d ||kr|||r|jt |||dz }||kr|||r3n| d vrz||dz k(rtd| d|dz }|||std| d ||krU|||rJ| jt |||dz }||kr#|||r3n|dz }|j| ||krt t ||} |r|| d <|r|| d<| r| | d<| S)Nc: t|y#t$rYywxYwNTF)rI ValueErrorss r,checkIntz GDAL2XYZ.parse..checkInt\s# A    c: t|y#t$rYywxYwr)rLrrs r, checkFloatz"GDAL2XYZ.parse..checkFloatcs# a  rr)rrrr%zMissing argument following z: zArgument following z should be a integer)rrz should be a floatrrrr)r5r4rHrIrLrr~parse) rrrrcountnew_argvirrrargkwargsrs r,rzGDAL2XYZ.parse[so  D     %iq'C// >#&A#b$IJJQQ(#&9#>R$STT%iHT!W$5$$Sa\2FA%iHT!W$566 >#&A#b$IJJQ!$q'*#&9#>P$QRR%iJtAw$7%%eDGn5FA%iJtAw$7& >#&A#b$IJJQ!$q'*#&9#>P$QRR%iJtAw$7%%eDGn5FA%iJtAw$7Q$?%iBx,X6 "+F;  #-F< #-F<  r|cf|jdrd|d<n|jdsd|d<|d=|S)Nrrr%)getrrs r,augment_kwargszGDAL2XYZ.augment_kwargss; ::j !"&F; K("#F;  :  r|c tdi|S)N)r{rs r,doitz GDAL2XYZ.doits!&!!r|) __name__ __module__ __qualname__rr rdictrrr __classcell__)rs@r,r~r~s8  b"4bH@T@D"r|r~c4tj|Sr()r~main)rs r,rrs :??4  r|__main__) NNr%NrFNNFT.)$r:rnumbersrtypingrrrrnumpyrBosgeorosgeo_utils.auxiliary.baser #osgeo_utils.auxiliary.gdal_argparser r osgeo_utils.auxiliary.numpy_utilr osgeo_utils.auxiliary.progressr rosgeo_utils.auxiliary.utilrrrrIrboolr{r~rrrexitrr|r,rsm@ 334NVDC "&*&')-4848"#'25t t t Xc] #t Xc]" # t  & t  ttx/01tx/01tt!t0te_tn}"z}"@hh! z CHHT#((^r|