njbN$&SSKrSSKrSSKJrJr SSKJr SSKJrJ r J r J r SSK r SSKJrJrJr SSKJr SSKJrJr SSKJr SS KJrJr SS KJr SS KJrJ r J!r!J"r"J#r#J$r$ SS K%J&r&J'r'J(r(J)r)J*r* "S S\5r+"SS\5r,\\ \RZ\+\ 4r.SSSSSSSSSSS\R^S4 S\&S\S\S\.S\\S\\0S\\1S\\ \1S\2S\2S\\3S\2S \S!\S"\2S#\ \ Rh\ Rh\ Rh44 S$jjr5SS\R^S4S\&S\S\S\\0S\\ \1S%\\,4 S&jjr6S'\7S(\S)\S*\2S+\2S,\24 S-jr8"S.S/\5r9\Rt4S0jr;\\rSrSr\"5r\"5r\"5rSrg)LocationInfoSRS?N) __name__ __module__ __qualname____firstlineno__r PixelLine SameAsDS_SRSSameAsDS_SRS_GeogCS__static_attributes__r"`/var/www/html/pdf-tiff/venv/lib/python3.13/site-packages/osgeo_utils/samples/gdallocationinfo.pyr r ?sI6L&r+r ch\rSrSr\"5r\"5r\"5r\"5r\"5r \"5r Sr g)LocationInfoOutputEr"N) r#r$r%r&r PixelLineValPixelLineValVerboseValOnlyXMLLifOnlyQuietr*r"r+r,r.r.Es-6L&fG &CfG FEr+r.FTgfilename_or_dsxysrs axis_order open_optionsovr_idx band_numsinline_xy_replacementreturn_ovr_pixel_linetransform_round_digitsallow_xy_outside_extent pixel_offset line_offset quiet_modereturnc [XS9n[U5(aUOSnUc[SUS35e[U[R 5(dU/n[U[R 5(dU/n[ U5[ U5:wa#[S[ U5S[ U535e[ U5n[Rn[U[R5(a[U[R5(d,[R"UUS9n[R"UUS9nSnUc[RnU[R:wGa>U[R:waUR5nSn[U[R 5(aUnO:U[R":XaUR%5nO ['X4S 9n[)UU5nUbVU(d"UR+5nUR+5nSn[-UX5 U bUR/XS 9 UR/XS 9 UR15n[2R4"U5nUc [S 5eUS US U--USU--USUSU--USU--p!SnUR6UR8nn[;UXvS9nUS R<US R>nnU(a UU- UU- 4OSunn[ U5n U U4n![AU5u n"[RB"U!U"S9n#U(+=(d U (+n$U$(a[RD"US :5(dZ[RD"UU:5(d<[RD"US :5(d[RD"UU:5(a&Sn%U (d [U%5eU(d [GU%5 US :XaUn&OU (aUU-nUn&OUU-n&US :XaUn'OU (aUU-nUn'OUU-n'S =n(n)[AU5un*n+[RB"U)U(/U+S9n,[I[KU&U'55HMun-un.n/[IU5H5un0n1[MU1S U.S- U/S- S S U,U*USS5 S :XdM'U,S S U#U0U-'M7 MO [OU5un2n3n4U2(a0[I[KU3U455Hun0un5n6U#U0U5-U6-U#U0'M XU#4$)Nr;zCould not open .zlen(x)=z should be the same as len(y)=)dtypeT)r:)outzFailed InvGeoTransform()r)r<)rLrL)shaperJz-Passed coordinates are not in dataset extent!g?)(rr Exception isinstancer__args__lennpfloat64ndarrayarrayr r'r( GetSpatialRefr CoordinateTransformationr) CloneGeogCSrrcopyrroundGetGeoTransformr InvGeoTransform RasterXSize RasterYSizerXSizeYSizeremptyanyprint enumeratezipr r)7r6r7r8r9r:r;r<r=r>r?r@rArBrC resample_algrDdsfilename point_countrJds_srsct points_srs geotransforminv_geotransformxsizeysizebands ovr_xsize ovr_ysize pixel_fact line_fact bnd_countrQnp_dtyperesults check_outsidemsgpixels_qlines_q buf_xsize buf_ysizebuf_typetypecodebuf_objidxpixellinebnd_idxband is_scaledscalesoffsetsscaleoffsets7 r,gdallocationinforSs$  ;B!-n!=!=~2H z/(1566 a++ , , C a++ , , C 1vQ'#a&)GAxPQQa&K JJE a $ $Jq"**,E,E HHQe $ HHQe $ $ {'' o''' /.. .%%'FB#s;;<</===!'!3!3!5J!(!DJ":v6~,AA,0) Q*)5GG2G:GG2G:))+ // =  #67 7 Q "21"5"9 9>2>>5E b) 5E 8>>58>>yI29U I-.vJ E I  $EB2FHhhhU(3G"NA*A&AM q1u U ++rvva!e}}qEz@R@R=&C.  #JQ  Zz>A~  Yi-I B2FHhhh 9-X>G'Hg(>?]eT&u-MGT!CK3J   )0 1  %#.@("8!>Ivw(1#fg2F(G $G_eV&w/%7&@GG )H =r+ output_modec Uc[RnU[R[R4;a[ SUS35eU[R:HnU[R [R [R4;n U (+n [XS9n [X5n[SU UUUUU US.UD6upn U RU RpU[R:HnU (GaRU[R:Xa [S5 [![#X55GHununnU(dUS:dUU :d US:dUU:a[SUSUS 35 M=U(a[%U5[%U5nnU[R:Xa[S US US 35 [!U5HunnU UUnU[R :Xa [U5 M/U[R :Xa[USUSU35 MXU[R:XdMn[S US35 [SU35 M GM U $)NzSorry, output mode z5 is not implemented yet. you may use the c++ version.rG)r6r7r8r=rjr>rDzReport:rz Location  z is off this file!z Location: (zP,zL)z Band :z Value: r")r.r5r3r4rRr2r0r1rrrrarbr GRIORA_NearestNeighbourrgrhriint)r6r7r8r;r=rjrkwargsrD print_moder>rkr|rsrtis_nearest_neighbourrrrrband_numvals r,gdallocationinfo_utilrs7(.. )--/A/I/IJJ!+.c d   2 8 88J""''..!J !+N  ;Bb,I$  !3   MA'>>2>>5'4+G+GG ,@@ @ ) "+CI"6 C%AIE>!85= %$/ABC'"%e*c$i4E"4"H"HHM%4&;<)29)=%GX!'*3/C"&8&@&@@c $(:(G(GGqau56$(:(N(NNz34 C512*>#70 Nr+rl longitudelatitudecoordtype_georefprint_xy print_valuesc[R"U[R5nUc[SU-5eU(aUnUnOn[R "5n U R UR55 U R5n [R"X5n U RX5upxn UR5n [R"U 5up[USUSU--USU--5n[USUSU--USU--5nU(a[SUU4-5 US:d&UUR:dUS:dUUR :a [S 5eUR#UUSS5nU(aF[%UR&5S:Xa[USS5 U$UHn[USS5 M U$) zk val_at_coord is a simplified version of gdallocationinfo. It accepts a single point and has less options. zCannot open %srrLrMrNrOrPz x=%d, y=%dz,Passed coordinates are not in dataset extent)r Open GA_ReadOnlyrRr SpatialReference ImportFromWkt GetProjectionr\r[TransformPointr_r`rrgrarb ReadAsArrayrUrQ)rlrrrrrrkXYr9 srsLatLongroheightrqsuccessrrr7r8resrs r, val_at_coordr)s 8T-- .B z(8344  ""$ "**,-__&  ) )* :**9?v%%'L"&"6"6|"DW Q "21"5"9 9Q  #a&)  Jc!fQi  Jr+cJ^\rSrSrU4SjrS\4SjrS\4SjrSr Sr U=r $)GDALLocationInfoi]ct>[TU]5 SUl[R"S5UlSUlg)NzRaster query toolz The gdallocationinfo utility provide a mechanism to query information about a pixel given its location in one of a variety of coordinate systems. Several reporting options are provided.)super__init__titletextwrapdedent descriptioninteractive_mode)self __class__s r,rGDALLocationInfo.__init__^s6 ( #?? b  !%r+rEc URnUR5nURSSSSS9 URSSSSS9 URS S SS S9 URS S SSS9 URSSSSS9 UR5nURSSS[SS9 URSSSSS9 URSSSSS9 URSSSSS9 URS S!S"S#S9 URS$S%SS&S9 URS'S(S)[S*S+S,9 URS-S.S/[S0S9 URS1S2S3S4/[SS5S69 URS7S8S9S:S*S;9 URSS?9 URS@SASB[ SCSD9 U$)ENz-xmlxml store_truezGThe output report will be XML formatted for convenient post processing.)destactionhelpz-lifonlylifonlyzThe only output is filenames production from the LocationInfo request against the database (i.e. for identifying impacted file from VRT).z-valonlyvalonlyzXThe only output is the pixel values of the selected pixel on each of the selected bands.z-plbplbzmThe output will be a series of lines in the form of Pixel,Line,band0,band1... for each of the selected bands.z-quietquietzNo output will be printed.z-l_srsr9srs_defz1The coordinate system of the input x, y location.)rmetavartyperz-geolocgeoloczIIndicates input x,y points are in the georeferencing system of the image.z -llgeolocllgeoloczmIndicates input x,y points are in the long, lat (geographic) based on the georeferencing system of the image.z-wgs84wgs84z/Indicates input x,y points are WGS84 long, lat.z-extent_strictrA store_falsezzIf set, input points outside the raster extent will raise an exception, otherwise a warning will be issued (unless quiet).z-interprjzXIf set, a Bilinear interpolation would be used, otherwise the NearestNeighbour sampling.z-br=r+zXSelects a band to query. Multiple bands can be listed. By default all bands are queried.)rrrnargsrz -overviewr<overview_levelzQuery the (overview_level)th overview (overview_level=1 is the 1st overview), instead of the base band. Note that the x,y location (if the coordinate system is pixel/line) must still be given with respect to the base band.z -axis_orderr:gis authorityzSX, Y Axis order: Traditional GIS, Authority complaint or otherwise utility default.)rchoicesrdefaultrz-oor;z NAME=VALUEz&Dataset open option (format specific).)rrrrr6rlz'The source GDAL raster datasource name.)rrrxyzx y*zseries of X Y pairs of location of target pixel. By default the coordinate system is pixel/line unless -l_srs, -wgs84 or -geoloc supplied.)rrrr)parseradd_mutually_exclusive_group add_argumentstrrfloat)rargvrgroups r, get_parserGDALLocationInfo.get_parserhs335  Z   J    k   .    -  335  D    \   ?    B    * A    k    k    $M    K(f     9    :    G   r+c[US5S:*UlUR(dD[R"USSSS25US'[R"USSSS25US'US(a[R US'O-US (a[R US'OUS (aS US'US n[U[5(aO[U[5(aUR5S :HnONUc%USSL=(a US[R :gnO&[U[5(aUnO[SU35e[U5nX!S 'US(a[RUS'OUS(a[R US'OmUS(a[R"US'OOUS(a[R$US'O1US(a[R&US'O[R(US'US(a[*R,O[*R.US'US US US US US US US US US U$)NrrLrrMr7r8rr9rrir:rzUnknown axis order rrrrrrrj)rUrrVrYr r(r)rSrrlowerboolrRrr.r3r4r2r0r5r1r GRIORA_Bilinearr)rrr: gis_orders r,augment_kwargsGDALLocationInfo.augment_kwargss9 #F4L 1Q 6$$((6$<1#56F3K((6$<1#56F3K ( +88F5M J +??F5M G_ F5ML) j/ 2 2 *c**&,,.%7 #"5M55M_%A%AAJ--& "5j\ BCC6yAJ)| %=$6$:$:F= ! I $6$>$>F= ! I $6$>$>F= ! E]$6$C$CF= ! G_$6$<$sysrenumrrnumbersrtypingrrrr numpyrVosgeor r r osgeo.gdal_arrayr osgeo_utils.auxiliary.array_utilrrosgeo_utils.auxiliary.baser#osgeo_utils.auxiliary.gdal_argparserr osgeo_utils.auxiliary.numpy_utilrosgeo_utils.auxiliary.osr_utilrrrrrrosgeo_utils.auxiliary.utilrrrrrr r.r[CoordinateTransformationOrSRSrrrrrrXrrrrrrrr#exitr"r+r,rs~@ 33&&.I3NV!d! !) # & & ?@!*.,0#'!)-"'"'.2$(22!MMMM ' M ) M 4. Mc]M &M M M%UOM"MMM !M" 2::rzz2:: -.#Mh$()-2204CCCC4. C  & C,-CL1111 1  1  1h]3z]3@hh) zHHT#((^r+