
    nj                         S SK r S SKrS SKJrJrJr  S r\R                  4S jr\	S:X  a#  \R                  " \" \R                  5      5        gg)    N)gdalogrosrc                      [        S5        g)Nz~Usage: wcs_virtds_params.py [-lyr_name name] [-tileindex field_name] [-t_srs srsdef] -src_srs_name field_name ogr_ds_tileindex   )print     a/var/www/html/pdf-tiff/venv/lib/python3.13/site-packages/osgeo_utils/samples/wcs_virtds_params.pyUsager   (   s    	 	I r
   c                 	   [         R                  " U 5      n S nS nSnS nS n[        U 5      nSnXv:  a  X   S:X  a  XvS-
  :  a
  US-   nX   nOjX   S:X  a  XvS-
  :  a
  US-   nX   nOPX   S:X  a  XvS-
  :  a
  US-   nX   nO6X   S:X  a  XvS-
  :  a
  US-   nX   nOX   S   S:X  a
  [        5       $ Uc  X   nUS-   nXv:  a  M  Ub  Uc
  [        5       $ [         R                  " U5      nUc  [        S	U-  5      eUR                  5       S:X  a  UR                  S5      n	OUc  [        S
5      eUR                  U5      n	U	R                  5       R                  U5      S:  a  [        SU-  5      eU	R                  5       R                  U5      S:  a  [        SU-  5      eU	R                  5       n
Ub  [        R                  " 5       nUR                  U5      S:w  a  [        SU-  5      eU
bo  U
R                  5       n[        R                  " 5       n
U
R                  U5        U
R                  5       nUR                  5       nX:w  a  [        SU< SU< S35      eOU
nUc  [        S5      eUR!                  5         UR#                  S 5      nUR%                  S 5      nUR'                  5       nUS:w  d  Uc  [        SU-  5      eSnSnSn U	R)                  5       nUc  OUR+                  U5      n[,        R.                  R1                  U5      (       d=  [,        R.                  R3                  [,        R.                  R5                  U5      U5      n[6        R                  " U5      nUc  [        S	U-  5      e[6        R8                  " US U5      nUc  [        SU< SU< 35      eUR;                  5       nUUS   -  nUUS   -  nS nUS-   nGM  US:X  a  [        S5      eUU-  nUU-  nU	R=                  5       u  nnnn[?        UU-
  U-  S-   5      n[?        UU-
  [A        U5      -  S-   5      nU	RC                  5       nURE                  5       RC                  5       S:w  a;  [G        SU< SU< SU	RC                  5       < SU< S3	5        [G        S5        SU-  n OUn [G        S UU UUUUUUUUUUU[A        U5      4-  5        g)!Nlocation   z	-lyr_namez
-tileindexz-src_srs_namez-t_srsr   -zcannot open %sz-lyr_name should be specifiedz,%s field cannot be found in layer definitionzinvalid value for -t_srs : %sz9-t_srs overrides the layer SRS in an incompatible way : (z, )zcannot fetch source SRSEPSGz.cannot fetch source SRS as EPSG:XXXX code : %szcannot warp z to    ztileindex is emptyg      ?zESRI ShapefilezLAYER
      NAME "z[_tileindex"
      TYPE POLYGON
      STATUS OFF
      CONNECTIONTYPE OGR
      CONNECTION ",z'"
      PROJECTION
        "+init=epsg:z"
      END
    END z%s_tileindexa   LAYER
      NAME "%s"
      TYPE RASTER
      STATUS ON
      TILEINDEX "%s"
      TILEITEM "%s"
      TILESRS "%s"
      PROJECTION
        "+init=epsg:%s"
      END
      METADATA
       "wcs_label"       "%s"
       "wcs_rangeset_name"   "Range 1"  ### required to support DescribeCoverage request
       "wcs_rangeset_label"  "My Label" ### required to support DescribeCoverage request
       "wcs_extent"      "%f %f %f %f"
       "wcs_size"        "%d %d"
       "wcs_resolution"  "%f %f"
      END
    END)$r   GeneralCmdLineProcessorlenr   Open	ExceptionGetLayerCountGetLayerGetLayerByNameGetLayerDefnGetFieldIndexGetSpatialRefr   SpatialReferenceSetFromUserInputExportToProj4AutoIdentifyEPSGGetAuthorityNameGetAuthorityCodeExportToWktGetNextFeatureGetFieldospathisabsjoindirnamer   AutoCreateWarpedVRTGetGeoTransform	GetExtentintabsGetName	GetDriverr   )!argvogr_ds_namelyr_nametileitemtilesrssrsnamenArgciArgogr_dslyrlyr_srssrslyr_srs_proj4	srs_proj4authority_nameauthority_codedst_wktcounterxresyresfeatgdal_ds_namegdal_dswarped_vrt_dsgtxminxmaxyminymaxxsizeysize	layername	tileindexs!                                    r   mainrV   /   s   &&t,DKHHGGIED
,:$	)9!8DzHZ<'D19,<!8DzHZ?*tai/?!8DjGZ8#qy(8!8DjGZ]c!7N *Kax1 ,4 gowXXk"F~(;677"ooa 		788##H-
''1A5FQRR
''014FPQQ!G""$(A-;gEFF #113M**,G$$]3#113M))+I) -1 
 
{122))$/N))$/NooG>#9H7RSSGDD
!!#< }}X.ww}}\**77<<(DlSL))L)?,|;<<00$H lGLMM**,11A+) , !|,--GODGOD"}}T4t4$;$&,-E4$;#d)+c12EI!!#'77 +s{{}n
F	
 	b	"Y.				( I
'"	
$L r
   __main__)r)   sysosgeor   r   r   r   r5   rV   __name__exitr	   r
   r   <module>r\      sL   D 
 
     hh xv zHHT#((^ r
   