gkdZddlZddlZddlmZddlmZmZmZm Z m Z m Z ddl m Z ddlmZmZddlmZmZddlmZdd lmZdd lmZdd lmZmZdd lmZm Z m!Z!dd l"m#Z#ddl$m%Z%Gdde Z&Gdde&Z'Gdde&Z(Gdde(Z)Gdde&Z*Gdde&Z+Gdde+Z,Gdde+Z-Gdd e+Z.id!e'd"e(d#e*d$e,d%e-d&e.d'e+d(e)d)e'd*e(d+e*d,e,d-e-d.e.d/e+d0e'd1e(d2e*d3e,d4e-d5e.d6e+d!ej^ze'd"ej^ze(d#ej^ze*d$ej^ze,d%ej^ze-d&ej^ze.d'ej^ze+i Z0y)7a The OGRGeometry is a wrapper for using the OGR Geometry class (see https://gdal.org/api/ogrgeometry_cpp.html#_CPPv411OGRGeometry). OGRGeometry may be instantiated when reading geometries from OGR Data Sources (e.g. SHP files), or when given OGC WKT (a string). While the 'full' API is not present yet, the API is "pythonic" unlike the traditional and "next-generation" OGR Python bindings. One major advantage OGR Geometries have over their GEOS counterparts is support for spatial reference systems and their transformation. Example: >>> from django.contrib.gis.gdal import OGRGeometry, OGRGeomType, SpatialReference >>> wkt1, wkt2 = 'POINT(-90 30)', 'POLYGON((0 0, 5 0, 5 5, 0 5)' >>> pnt = OGRGeometry(wkt1) >>> print(pnt) POINT (-90 30) >>> mpnt = OGRGeometry(OGRGeomType('MultiPoint'), SpatialReference('WGS84')) >>> mpnt.add(wkt1) >>> mpnt.add(wkt1) >>> print(mpnt) MULTIPOINT (-90 30,-90 30) >>> print(mpnt.srs.name) WGS 84 >>> print(mpnt.srs.proj) +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs >>> mpnt.transform(SpatialReference('NAD27')) >>> print(mpnt.proj) +proj=longlat +ellps=clrk66 +datum=NAD27 +no_defs >>> print(mpnt) MULTIPOINT (-89.99993037860248 29.99979788655764,-89.99993037860248 29.99979788655764) The OGRGeomType class is to make it easy to specify an OGR geometry type: >>> from django.contrib.gis.gdal import OGRGeomType >>> gt1 = OGRGeomType(3) # Using an integer for the type >>> gt2 = OGRGeomType('Polygon') # Using a string >>> gt3 = OGRGeomType('POLYGON') # It's case-insensitive >>> print(gt1 == 3, gt1 == 'Polygon') # Equivalence works w/non-OGRGeomType objects True True N)b2a_hex)byrefc_char_pc_doublec_ubytec_void_p string_at)GDALBase)Envelope OGREnvelope) GDALException SRSException) OGRGeomType)geom)srs)CoordTransformSpatialReference) hex_regex json_regex wkt_regex)RemovedInDjango60Warning) force_bytesc0eZdZdZej ZdDdZdZdZ e dZ e dZ e dZe d Ze d Zd Zd Zd ZdZdZdZedZedZej4dZedZedZedZedZedZedZ edZ!edZ"edZ#edZ$edZ%dZ&ed Z'd!Z(d"Z)d#Z*ee)e*Z+d$Z,d%Z-ee,e-Z.d&Z/ed'Z0ed(Z1ed)Z2ed*Z3e3Z4ed+Z5ed,Z6ed-Z7ed.Z8ed/Z9d0Z:d1Z;dEd2Zd5Z?d6Z@d7ZAd8ZBd9ZCd:ZDd;ZEdDd<ZFed=ZGed>ZHd?ZId@ZJdAZKdBZLedCZMy)F OGRGeometryzEncapsulate an OGR geometry.Nc *t|t}|r5tj|r t t j |}d}|r]tj|}tj|}|r|drt|d}|djdk(rftjt|dj}tj|t!t#|dj%n$tj&t!t#|dj%dt!t)}n|r |j+|j%}nt|tjt|j}nt|tr|j-|}n`t|tr tj|j}n0t||j.r|}nt1dt3|z|st1d|z||_|r||_t8j;|j<jx}t?d |j<||_ y) z=Initialize Geometry on either WKT or an OGR pointer as input.Fsridtype LINEARRINGwktNz4Invalid input type for OGR Geometry construction: %sz)Cannot create OGR Geometry from input: %szUnsupported geometry type: )! isinstancestrrmatch memoryviewbytesfromhexrrintuppercapi create_geomrnum import_wktrrencodefrom_wktr _from_json _from_wkbptr_typer rptrr GEO_CLASSESget geom_type TypeError __class__)self geom_inputr str_instancewkt_mjson_mg geo_classs `/var/www/skyplay_api_hub/venv/lib/python3.12/site-packages/django/contrib/gis/gdal/geometries.py__init__zOGRGeometry.__init__Ds!*c2  IOOJ7#EMM*$=>J L OOJ/E%%j1F=eFm,C=&&(L8((U6])C)G)GHAOOAuXeEl6I6I6K-L'MN huU|':':'<=>eHJFWAOOJ$5$5$78J'$$[%<%@%@A  J /z*A  K 0  0A  DMM 2AFz"# ;jH  DH%););< __getstate__zOGRGeometry.__getstate__s/hh ''CCTXX##r@c|\}}tj|dttt |}|s t d||_||_y)Nz.Invalid OGRGeometry loaded from pickled state.)r(from_wkbrrlenr r1r)r7staterCrr1s r> __setstate__zOGRGeometry.__setstate__sGSmmCuXZ'8#c(C PQ Qr@cxtjt|dtt t |SrB)r(rGr$rrrH)clsr8s r>r/zOGRGeometry._from_wkbs,}} * tU8:%6J  r@c,tj|SrB)r( from_jsonr8s r>r.zOGRGeometry._from_jsons~~j))r@cd|\}}}}td|d|d|d|d|d|d|d|d|d|dS)z2Construct a Polygon from a bounding box (4-tuple).z POLYGON(( z, z)))r)rLbboxx0y0x1y1s r> from_bboxzOGRGeometry.from_bboxs8BB2r2r2r2r2 7  r@cPttjt|SrB)rr.rrOs r>rNzOGRGeometry.from_jsons;11+j2IJKKr@cJ|tjt|SrB)r(from_gmlr)rL gml_strings r>rZzOGRGeometry.from_gmls4==Z!89::r@c$|j|S)z'Return the union of the two geometries.)unionr7others r>__or__zOGRGeometry.__or__szz%  r@c$|j|S)z7Return the intersection of this Geometry and the other.) intersectionr^s r>__and__zOGRGeometry.__and__s  ''r@c$|j|S)z2Return the difference this Geometry and the other.) differencer^s r>__sub__zOGRGeometry.__sub__su%%r@c$|j|S)z?Return the symmetric difference of this Geometry and the other.)sym_differencer^s r>__xor__zOGRGeometry.__xor__s""5))r@cHt|txr|j|S)z$Is this Geometry equal to the other?)r requalsr^s r>__eq__zOGRGeometry.__eq__s%-D$++e2DDr@c|jS)z*WKT is used for the string representation.)rr7s r>__str__zOGRGeometry.__str__s xxr@c@tj|jS)z5Return 0 for points, 1 for lines, and 2 for surfaces.)r(get_dimsr1rns r> dimensionzOGRGeometry.dimension}}TXX&&r@c@tj|jS)z0Return the coordinate dimension of the Geometry.)r( get_coord_dimr1rns r> coord_dimzOGRGeometry.coord_dim!!$((++r@cd}tj|td|dvr tdt j |j |y)z.Set the coordinate dimension of this Geometry.z5coord_dim setter is deprecated. Use set_3d() instead.) stacklevel)ryz(Geometry dimension must be either 2 or 3N)warningswarnr ValueErrorr( set_coord_dimr1)r7dimmsgs r>rvzOGRGeometry.coord_dimsCF c3B f GH H 488S)r@c@tj|jS)z/Return the number of elements in this Geometry.)r(get_geom_countr1rns r> geom_countzOGRGeometry.geom_counts""488,,r@c@tj|jS)z-Return the number of Points in this Geometry.)r(get_point_countr1rns r> point_countzOGRGeometry.point_counts##DHH--r@c|jS)z7Alias for `point_count` (same name method in GEOS API.)rrns r> num_pointszOGRGeometry.num_pointsr@c|jS)zAlias for `point_count`.rrns r> num_coordszOGRGeometry.num_coordsrr@cRttj|jS)z"Return the Type for this Geometry.)rr( get_geom_typer1rns r>r4zOGRGeometry.geom_types4--dhh788r@c@tj|jS)z!Return the Name of this Geometry.)r( get_geom_namer1rns r> geom_namezOGRGeometry.geom_namerwr@c@tj|jS)zHReturn the area for a LinearRing, Polygon, or MultiPolygon; 0 otherwise.)r(get_arear1rns r>areazOGRGeometry.arearsr@c vttj|jt t S)z&Return the envelope for this Geometry.)r r( get_enveloper1rr rns r>envelopezOGRGeometry.envelopes'))$((E+-4HIJJr@c@tj|jSrB)r(is_emptyr1rns r>emptyzOGRGeometry.empty s}}TXX&&r@c.|jjS)zCReturn the envelope as a 4-tuple, instead of as an Envelope object.)rtuplerns r>extentzOGRGeometry.extents}}"""r@c@tj|jS)z.Return True if the geometry has Z coordinates.)r(is_3dr1rns r>rzOGRGeometry.is_3dszz$((##r@c|dur!tj|jdy|dur!tj|jdytd|d)z'Set if this geometry has Z coordinates.TFrz*Input to 'set_3d' must be a boolean, got ''.N)r(set_3dr1r~r7values r>rzOGRGeometry.set_3dsJ D= KK! $ e^ KK! $I%RTUV Vr@c@tj|jS)z.Return True if the geometry has M coordinates.)r( is_measuredr1rns r>rzOGRGeometry.is_measured#))r@c|dur!tj|jdy|dur!tj|jdytd|d)z'Set if this geometry has M coordinates.TrFrz0Input to 'set_measured' must be a boolean, got 'rN)r( set_measuredr1r~rs r>rzOGRGeometry.set_measured(sQ D=   dhh * e^   dhh *B5)2N r@c tj|j}tt j |S#t $rYywxYw)z/Return the Spatial Reference for this Geometry.N)r( get_geom_srsr1rsrs_api clone_srsr)r7srs_ptrs r>_get_srszOGRGeometry._get_srs6sC ''1G#G$5$5g$>? ?  s_set_srszOGRGeometry._set_srs>sq c+ ,ggG c3Z (!#&BffG [GIDQTIU  '*r@c8|j}|r |jSyrB)rrrDs r> _get_sridzOGRGeometry._get_sridTshh 88Or@cLt|ts|||_ytd)Nz!SRID must be set with an integer.)r r&rr5)r7rs r> _set_sridzOGRGeometry._set_sridZs# dC DLDH?@ @r@cDddlm}|j|jS)Nr GEOSGeometry)django.contrib.gis.geosrr/rCr7rs r> _geos_ptrzOGRGeometry._geos_ptrcs8%%dhh//r@cPddlm}||j|jS)z3Return a GEOSGeometry object from this OGRGeometry.rr)rrrrrs r>geoszOGRGeometry.geoshs 9DNN,dii88r@c@tj|jS)z.Return the GML representation of the Geometry.)r(to_gmlr1rns r>gmlzOGRGeometry.gmlos{{488$$r@cHt|jjS)zhexzOGRGeometry.hextstxx &&((r@c@tj|jS)zE Return the GeoJSON representation of this Geometry. )r(to_jsonr1rns r>jsonzOGRGeometry.jsonys ||DHH%%r@cBtj|jdS)z.Return the KML representation of the Geometry.N)r(to_kmlr1rns r>kmlzOGRGeometry.kmls{{488T**r@c@tj|jS)z"Return the size of the WKB buffer.)r( get_wkbsizer1rns r>wkb_sizezOGRGeometry.wkb_sizerr@c tjdk(rd}nd}|j}t|z}|jrt j nt j}||j|t|tt||S)z.Return the WKB representation of the Geometry.littlerr) sys byteorderrrrr( to_iso_wkbto_wkbr1rr#r )r7rszbufrs r>rCzOGRGeometry.wkbsn ==H $II ]]|%)$4$4$++txxE#J/)C,--r@c|jrtjntj}||jt t S)z.Return the WKT representation of the Geometry.)rr( to_iso_wktto_wktr1rr)r7rs r>rzOGRGeometry.wkts3 %)$4$4$++dhhhj 122r@c|j}|r(|jrd|jd|jS|jS)z/Return the EWKT representation of the Geometry.zSRID=;)rrrrDs r>ewktzOGRGeometry.ewkts2hh 388#&88TXX6 688Or@chttj|j|jS)zClone this OGR Geometry.)rr( clone_geomr1rrns r>clonezOGRGeometry.clones 4??4884dhh??r@cBtj|jy)z If there are any rings within this geometry that have not been closed, this routine will do so by adding the starting point at the end. N)r(geom_close_ringsr1rns r> close_ringszOGRGeometry.close_ringss dhh'r@c|r#|j}|j||St|tr+t j |j |j yt|tr+t j|j |j yt|ttfr6t|}t j|j |j ytd)aq Transform this geometry to a different spatial reference system. May take a CoordTransform object, a SpatialReference object, string WKT or PROJ, and/or an integer SRID. By default, return nothing and transform the geometry in-place. However, if the `clone` keyword is set, return a transformed clone of this geometry. zUTransform only accepts CoordTransform, SpatialReference, string, and integer objects.N) r transformr rr(geom_transformr1rgeom_transform_tor&r!r5)r7 coord_transrkloners r>rzOGRGeometry.transforms JJLE OOK (L k> 2   +// :  %5 6  " "488[__ =  c3Z 0!+.B  " "488RVV 4A r@crt|ts td||j|jS)zA generalized function for topology operations, takes a GDAL function and the other geometry to perform the operation on.z _topologyzOGRGeometry._topologys5%-N  DHHeii((r@cB|jtj|S)z7Return True if this geometry intersects with the other.)rr(ogr_intersectsr^s r> intersectszOGRGeometry.intersectss~~d11599r@cB|jtj|S)z8Return True if this geometry is equivalent to the other.)rr( ogr_equalsr^s r>rkzOGRGeometry.equals~~doou55r@cB|jtj|S)zBReturn True if this geometry and the other are spatially disjoint.)rr( ogr_disjointr^s r>disjointzOGRGeometry.disjoint~~d//77r@cB|jtj|S)z/Return True if this geometry touches the other.)rr( ogr_touchesr^s r>toucheszOGRGeometry.touches~~d..66r@cB|jtj|S)z/Return True if this geometry crosses the other.)rr( ogr_crossesr^s r>crosseszOGRGeometry.crossesrr@cB|jtj|S)z1Return True if this geometry is within the other.)rr( ogr_withinr^s r>withinzOGRGeometry.withinrr@cB|jtj|S)z0Return True if this geometry contains the other.)rr( ogr_containsr^s r>containszOGRGeometry.containsrr@cB|jtj|S)z0Return True if this geometry overlaps the other.)rr( ogr_overlapsr^s r>overlapszOGRGeometry.overlapsrr@ct|tr1t||j|j|jSt||j|jS)z?A helper routine for the OGR routines that generate geometries.)r rr1r)r7gen_funcr_s r>_geomgenzOGRGeometry._geomgensF e[ )x%))boundaryzOGRGeometry.boundarys}}T..//r@c@|jtjS)zk Return the smallest convex Polygon that contains all the points in this Geometry. )r r(geom_convex_hullrns r> convex_hullzOGRGeometry.convex_hulls }}T2233r@cB|jtj|S)z Return a new geometry consisting of the region which is the difference of this geometry and the other. )r r( geom_diffr^s r>rezOGRGeometry.differences }}T^^U33r@cB|jtj|S)zx Return a new geometry consisting of the region of intersection of this geometry and the other. )r r(geom_intersectionr^s r>rbzOGRGeometry.intersection#s }}T33U;;r@cB|jtj|S)zq Return a new geometry which is the symmetric difference of this geometry and the other. )r r( geom_sym_diffr^s r>rhzOGRGeometry.sym_difference*s }}T//77r@cB|jtj|S)z{ Return a new geometry consisting of the region which is the union of this geometry and the other. )r r( geom_unionr^s r>r]zOGRGeometry.union1s }}T__e44r@cttd}tj|j|j|S)z.Return the centroid (a Point) of this Polygon.Point)rrr( get_centroidr1)r7ps r>centroidzOGRGeometry.centroid8s0  G, - $((AEE*r@rB)F)N__name__ __module__ __qualname____doc__r( destroy_geom destructorr?rErJ classmethodr/ staticmethodr.rWrNrZr`rcrfrirlropropertyrrrvsetterrrrrr4rrrrrrrrrrrrrrrrrrrrgeojsonrrrCrrrrrrrrkrrrrrrr r rrerbrhr]rr@r>rr?se&""J?#D$  **  LL;; ! ( & *E '',, **--..    99,,''KK ''##$$W** +& 8X &C A Iy )D0 99 %%))&& G ++**.."33@(: ):6877688 =00444<85r@rceZdZfdZedZedZedZedZ edZ edZ e Z xZ S)rcxddlm}|jr|jj St |S)Nr)r)django.contrib.gisrrr _create_emptysuperr)r7rr6s r>rzPoint._geos_ptrCs,+-1ZZtzz'')PUW=N=PPr@cRtjtdjS)Npoint)r(r)rr*)rLs r>r,zPoint._create_emptyHs G 4 8 899r@cBtj|jdS)z'Return the X coordinate for this Point.r)r(getxr1rns r>xzPoint.xLyy1%%r@cBtj|jdS)z'Return the Y coordinate for this Point.r)r(getyr1rns r>yzPoint.yQr3r@c\|jr tj|jdSy)z'Return the Z coordinate for this Point.rN)rr(getzr1rns r>zzPoint.zVs% ::99TXXq) ) r@c\|jr tj|jdSy)z'Return the M coordinate for this Point.rN)rr(getmr1rns r>mzPoint.m\s'   99TXXq) ) r@cz|jr:|jr.|j|j|j|j fS|jr#|j|j|jfS|jr#|j|j|j fS|j|jfS)zReturn the tuple of this point.)rrr2r6r9r<rns r>rz Point.tuplebs ::$**664664664661 1 ::66466466) )   66466466) )vvtvv~r@)rrrrr#r,r%r2r6r9r<rcoords __classcell__)r6s@r>rrBsQ ::&&&&** ** Fr@rcreZdZdZdZedZeZdZedZ edZ edZ edZ y ) LineStringc d|cxkr|jkrWnnSttttf\}}}}tj|j|t |t |t |t ||j r:|jr.|j|j|j|jfS|j r#|j|j|jfS|jr#|j|j|jfS|j}|dk(r |jfS|dk(r|j|jfSytd|z)z$Return the Point at the given index.rrryz>Index out of range when accessing points of a line string: %s.N) rrr( get_pointr1rrrrrv IndexError)r7indexr2r6r9r<rs r> __getitem__zLineString.__getitem__qs!  ((( (!XZXZGJAq!Q NN488UE!HeAha%PQ( Szzd..ww!''99zzww00ww00..Caxz!))PSXX r@c|jS)z.Return the number of points in the LineString.rrns r>__len__zLineString.__len__sr@cPtfdttDS)z3Return the tuple representation of this LineString.c3(K|] }| ywrBr(.0ir7s r> z#LineString.tuple..s7T!W7s)rrangerHrns`r>rzLineString.tuples7eCI&6777r@crtt|Dcgc]}||j|c}Scc}w)zp Internal routine that returns a sequence (list) corresponding with the given function. )rOrHr1)r7rrMs r>_listarrzLineString._listarrs, ,1T+;r2z LineString.x}}TYY''r@c@|jtjS)z#Return the Y coordinates in a list.)rQr(r5rns r>r6z LineString.yrSr@cZ|jr|jtjSy)z#Return the Z coordinates in a list.N)rrQr(r8rns r>r9z LineString.zs# ::==+ + r@cZ|jr|jtjSy)z#Return the M coordinates in a list.N)rrQr(r;rns r>r<z LineString.ms%   ==+ + r@N) rrrrFrHr%rr>rQr2r6r9r<r(r@r>rArAps|* 88F=((((,, ,,r@rAc eZdZy) LinearRingNrrrr(r@r>rXrXr@rXcPeZdZdZdZedZeZedZeZ edZ y)Polygonc|jS)z4Return the number of interior rings in this Polygon.rrns r>rHzPolygon.__len__ r@cd|cxkr|jkrJnnGttjtj|j ||j Std|z)z$Get the ring at the specified index.rz9Index out of range when accessing rings of a polygon: %s.rrr(r get_geom_refr1rrDr7rEs r>rFzPolygon.__getitem__sZ  ' ' 1 1$((E BCTXX KeS r@c |dS)z!Return the shell of this Polygon.rr(rns r>shellz Polygon.shells Awr@cRtfdtjDS)z/Return a tuple of LinearRing coordinate tuples.c3<K|]}|jywrBrrKs r>rNz Polygon.tuple..CqT!W]]CrrOrrns`r>rz Polygon.tupleCE$//,BCCCr@cRtfdtjDS)z,Return the number of Points in this Polygon.c3<K|]}|jywrBrrKs r>rNz&Polygon.point_count..G147&&GrjsumrOrrns`r>rzPolygon.point_countGdoo0FGGGr@N) rrrrHrFr%re exterior_ringrr>rr(r@r>r\r\sV M DDF HHr@r\cFeZdZdZdZdZdZedZedZ e Z y)GeometryCollectionzThe Geometry Collection class.cd|cxkr|jkrJnnGttjtj|j ||j Std|z)z(Get the Geometry at the specified index.rz?Index out of range when accessing geometry in a collection: %s.rarcs r>rFzGeometryCollection.__getitem__s]  ' ' 1 1$((E BCTXX Q r@c|jS)zrHzGeometryCollection.__len__r_r@ct|trst||jr2|D],}tj|j |j .ytj|j |j yt|t r6t|}tj|j |j ytd)z-Add the geometry to this Geometry Collection.zMust add an OGRGeometry.N)r rr6r(add_geomr1r!r )r7rr<tmps r>addzGeometryCollection.adds dK ($/3AMM$((AEE23 dhh1 c "d#C MM$((CGG , :; ;r@cRtfdtjDS)z8Return the number of Points in this Geometry Collection.c3<K|]}|jywrBrrKs r>rNz1GeometryCollection.point_count..rorjrprns`r>rzGeometryCollection.point_countrrr@cRtfdtjDS)z:Return a tuple representation of this Geometry Collection.c3<K|]}|jywrBrhrKs r>rNz+GeometryCollection.tuple..rirjrkrns`r>rzGeometryCollection.tuplerlr@N) rrrr rFrHr{r%rrr>r(r@r>rurusG$  <HH DDFr@ruc eZdZy) MultiPointNrYr(r@r>rr rZr@rc eZdZy)MultiLineStringNrYr(r@r>rrrZr@rc eZdZy) MultiPolygonNrYr(r@r>rrrZr@rrryr{eiiiiiiii i i i i i i )1r rr|binasciirctypesrrrrrr django.contrib.gis.gdal.baser django.contrib.gis.gdal.enveloper r django.contrib.gis.gdal.errorr r django.contrib.gis.gdal.geomtyper"django.contrib.gis.gdal.prototypesrr(rrdjango.contrib.gis.gdal.srsrrdjango.contrib.gis.geometryrrrdjango.utils.deprecationrdjango.utils.encodingrrrrArXr\rurrrwkb25bitr2r(r@r>rs'R JJ1BE8;=HHH=-(F+K+\<,<,@  #Hk#HN,,` #  (  % uzwz    |  % * ' * / ,   %!" *#$ '*/,  e j g j o l 0; r@