
    nj                         S SK r S SKJr  S r\ R                  4S jr\S:X  a#  \ R                  " \" \ R                  5      5        gg)    N)ogrc                  0    [        S5        [        S5        g)NzHUsage: vec_tr_spat.py [-spat xmin ymin xmax ymax] infile outfile [layer]    )print     [/var/www/html/pdf-tiff/venv/lib/python3.13/site-packages/osgeo_utils/samples/vec_tr_spat.pyUsager   &   s    	
TU	"Ir	   c                 \   S nS nS nSnU[        U 5      :  a  X   nUS:X  aF  [        XS-      5      n[        XS-      5      n[        XS-      5      n[        XS-      5      n	US-   nOUc  UnOUc  UnOUc  UnO
[        5       $ US-   nU[        U 5      :  a  M  Uc
  [        5       $ [        R                  " USS9n
Ub  U
R                  U5      nOU
R                  S5      nUR                  5       n[        R                  " S5      nUR                  U5        UR                  U5      nUR                  UR                  5       UR                  5       UR                  5       S	9nUR                  5       n[!        U5       H  nUR#                  U5      n[        R$                  " UR                  5       UR'                  5       5      nUR)                  UR+                  5       5        UR-                  UR/                  5       5        UR1                  U5        M     UR3                  WWWW	5        [        R4                  " [        R6                  5      nUR9                  Xg5        UR9                  X5        UR9                  X5        UR9                  Xi5        UR9                  Xg5        [        R4                  " [        R:                  5      nUR=                  U5        UR?                  5       nUb  URA                  5       nURC                  U5      S:w  aT  [        RD                  " UR                  5       S
9nURG                  U5        URI                  U5        URK                  5         URK                  5         UR?                  5       nUb  M  URK                  5         U
RK                  5         g)N   z-spatr         r   )updatezESRI Shapefile)	geom_typesrs)feature_def)&lenintr   r   OpenGetLayerByNameGetLayerGetLayerDefnGetDriverByNameDeleteDataSourceCreateDataSourceCreateLayerGetNameGetGeomTypeGetSpatialRefGetFieldCountrangeGetFieldDefn	FieldDefnGetTypeSetWidthGetWidthSetPrecisionGetPrecisionCreateFieldSetSpatialFilterRectGeometrywkbLinearRingAddPoint_2D
wkbPolygonAddGeometryGetNextFeatureGetGeometryRef	IntersectFeatureSetFromCreateFeatureDestroy)argvinfileoutfile
layer_nameiargs_minxs_minys_maxxs_maxyin_dsin_layerin_defn
shp_drivershp_ds	shp_layerin_field_count	fld_indexsrc_fdfdring	filt_geomin_featgeomout_feats                            r
   mainrQ   ,   s/   FGJ	A
c$i-g'>!e%F!e%F!e%F!e%FAA^F_GJ 7NE- c$i-0 w
 HHVA&E''
3>>!$##%G
 $$%56J(((1F""W%8%8%:@V@V@X # I **,N>*	%%i0]]6>>+V^^-=>
FOO%&
++-.b! + !!&&&&A
 <<))*DV$V$V$V$V$S^^,I$
 %%'G

%%'>>)$){{y/E/E/GHHW%##H-))+ 
" NN	MMOr	   __main__)sysosgeor   r   r8   rQ   __name__exitr   r	   r
   <module>rW      sE   >  
 hh m` zHHT#((^ r	   