X/^@szdZddlZddlmZmZddlmZmZmZm Z m Z m Z ddl m Z ddlmZmZddlmZmZmZddlmZdd lmZmZdd lmZmZdd lm Z m!Z!m"Z"dd l#m$Z$dd l%m&Z&Gddde Z'Gddde'Z(Gddde'Z)Gddde)Z*Gddde'Z+Gddde'Z,Gddde,Z-Gddde,Z.Gddde,Z/ie(d 6e)d!6e+d"6e-d#6e.d$6e/d%6e,d&6e*d'6e(d ej06e)d!ej06e+d"ej06e-d#ej06e.d$ej06e/d%ej06e,d&ej06Z1dS)(a| The OGRGeometry is a wrapper for using the OGR Geometry class (see http://www.gdal.org/classOGRGeometry.html). 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.999930378602485 29.999797886557641,-89.999930378602485 29.999797886557641) 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)a2b_hexb2a_hex)byrefc_char_pc_doublec_ubytec_void_p string_at)GDALBase)Envelope OGREnvelope) GDALException OGRIndexError SRSException) OGRGeomType)geomsrs)CoordTransformSpatialReference) hex_regex json_regex wkt_regex)six)rangec@sveZdZdZdddZddZddZd d Zed d Z d dZ ddZ ddZ ddZ ddZddZddZeddZddZdd ZeeeZed!d"Zed#d$Zed%d&Zed'd(Zed)d*Zed+d,Zed-d.Zed/d0Zed1d2Zd3d4Zd5d6Z eee Z!d7d8Z"d9d:Z#ee"e#Z$ed;d<Z%ed=d>Z&ed?d@Z'edAdBZ(e(Z)edCdDZ*edEdFZ+edGdHZ,edIdJZ-edKdLZ.dMdNZ/dOdPZ0dQdRdSZ1dTdUZ2dVdWZ3dXdYZ4dZd[Z5d\d]Z6d^d_Z7d`daZ8dbdcZ9dddeZ:ddfdgZ;edhdiZ<edjdkZ=dldmZ>dndoZ?dpdqZ@drdsZAdS)t OGRGeometryz'Generally encapsulates an OGR geometry.NcCsst|tj}|rQtj|rQtjt|jj}d}n|rt j|}t j|}|rH|j drt |j d}n|j djdkr t jt|j dj}t j|tt|j djqt jtt|j djdtt}q"|rft j|j}q"t|t jt|j}nt|tjrt jt|dttt|}nUt|trt j|j}n1t||jr |}ntdt||sAtdt|n||_|r\||_nt |j!j|_"dS) z>Initializes Geometry on either WKT or an OGR pointer as input.FsridtypeZ LINEARRINGwktNz4Invalid input type for OGR Geometry construction: %sz)Cannot create OGR Geometry from input: %s)# isinstancer string_typesrmatch memoryviewrupperencoderrgroupintcapiZ create_geomrnumZ import_wktrrZfrom_wktrZ from_jsonfrom_wkbbyteslenZptr_typer rstrptrr GEO_CLASSES geom_type __class__)selfZ geom_inputrZ str_instanceZwkt_mZjson_mgr2P/home/ubuntu/projects/ifolica/build/django/django/contrib/gis/gdal/geometries.py__init__As>! !.< 0   zOGRGeometry.__init__c Cs3ytj|jWnttfk r.YnXdS)zDeletes this Geometry.N)r&Z destroy_geomZ_ptrAttributeError TypeError)r0r2r2r3__del__zszOGRGeometry.__del__cCs4|j}|r|j}nd}t|j|fS)N)rrr)wkb)r0rr2r2r3 __getstate__s   zOGRGeometry.__getstate__cCs^|\}}tj|dttt|}|sHtdn||_||_dS)Nz.Invalid OGRGeometry loaded from pickled state.)r&r(rrr*r r,r)r0stater8rr,r2r2r3 __setstate__s  ' zOGRGeometry.__setstate__c Cs>|\}}}}td||||||||||f S)z3Constructs a Polygon from a bounding box (4-tuple).z,POLYGON((%s %s, %s %s, %s %s, %s %s, %s %s)))r)clsZbboxZx0Zy0Zx1y1r2r2r3 from_bboxszOGRGeometry.from_bboxcCs |j|S)z(Returns the union of the two geometries.)union)r0otherr2r2r3__or__szOGRGeometry.__or__cCs |j|S)z8Returns the intersection of this Geometry and the other.) intersection)r0r@r2r2r3__and__szOGRGeometry.__and__cCs |j|S)z2Return the difference this Geometry and the other.) difference)r0r@r2r2r3__sub__szOGRGeometry.__sub__cCs |j|S)z?Return the symmetric difference of this Geometry and the other.)sym_difference)r0r@r2r2r3__xor__szOGRGeometry.__xor__cCs$t|tr|j|SdSdS)z$Is this Geometry equal to the other?FN)rrequals)r0r@r2r2r3__eq__s zOGRGeometry.__eq__cCs ||k S)zTests for inequality.r2)r0r@r2r2r3__ne__szOGRGeometry.__ne__cCs|jS)z*WKT is used for the string representation.)r)r0r2r2r3__str__szOGRGeometry.__str__cCstj|jS)z6Returns 0 for points, 1 for lines, and 2 for surfaces.)r&Zget_dimsr,)r0r2r2r3 dimensionszOGRGeometry.dimensioncCstj|jS)z1Returns the coordinate dimension of the Geometry.)r&Z get_coord_dimr,)r0r2r2r3_get_coord_dimszOGRGeometry._get_coord_dimcCs2|dkrtdntj|j|dS)z/Sets the coordinate dimension of this Geometry.z(Geometry dimension must be either 2 or 3N)rNrO) ValueErrorr&Z set_coord_dimr,)r0dimr2r2r3_set_coord_dims zOGRGeometry._set_coord_dimcCstj|jS)z(The number of elements in this Geometry.)r&Zget_geom_countr,)r0r2r2r3 geom_countszOGRGeometry.geom_countcCstj|jS)z.Returns the number of Points in this Geometry.)r&Zget_point_countr,)r0r2r2r3 point_countszOGRGeometry.point_countcCs|jS)z7Alias for `point_count` (same name method in GEOS API.))rT)r0r2r2r3 num_pointsszOGRGeometry.num_pointscCs|jS)zAlais for `point_count`.)rT)r0r2r2r3 num_coordsszOGRGeometry.num_coordscCsttj|jS)z#Returns the Type for this Geometry.)rr&Z get_geom_typer,)r0r2r2r3r.szOGRGeometry.geom_typecCstj|jS)z"Returns the Name of this Geometry.)r&Z get_geom_namer,)r0r2r2r3 geom_nameszOGRGeometry.geom_namecCstj|jS)zIReturns the area for a LinearRing, Polygon, or MultiPolygon; 0 otherwise.)r&Zget_arear,)r0r2r2r3areaszOGRGeometry.areacCs"ttj|jttS)z'Returns the envelope for this Geometry.)r r&Z get_enveloper,rr )r0r2r2r3envelopeszOGRGeometry.envelopecCs |jjS)zDReturns the envelope as a 4-tuple, instead of as an Envelope object.)rYtuple)r0r2r2r3extentszOGRGeometry.extentc CsFy)tj|j}ttj|SWntk rAdSYnXdS)z0Returns the Spatial Reference for this Geometry.N)r&Z get_geom_srsr,rsrs_apiZ clone_srsr)r0srs_ptrr2r2r3_get_srss  zOGRGeometry._get_srscCst|tr|j}n\t|tjtjrLt|}|j}n+|dkrad}ntdt|tj |j|dS)z,Sets the SpatialReference for this geometry.Nz7Cannot assign spatial reference with object of type: %s) rrr,r integer_typesrr6rr&Z assign_srs)r0rr]srr2r2r3_set_srss     zOGRGeometry._set_srscCs|j}|r|jSdS)N)rr)r0rr2r2r3 _get_srids zOGRGeometry._get_sridcCs:t|tjs|dkr*||_n tddS)Nz!SRID must be set with an integer.)rrr_rr6)r0rr2r2r3 _set_srid"s zOGRGeometry._set_sridcCs#ddlm}||j|jS)z4Returns a GEOSGeometry object from this OGRGeometry.r) GEOSGeometry)Zdjango.contrib.gis.geosrdr8r)r0rdr2r2r3geos+szOGRGeometry.geoscCstj|jS)z/Returns the GML representation of the Geometry.)r&Zto_gmlr,)r0r2r2r3gml1szOGRGeometry.gmlcCst|jjS)z=Returns the hexadecimal representation of the WKB (a string).)rr8r")r0r2r2r3hex6szOGRGeometry.hexcCstj|jS)zF Returns the GeoJSON representation of this Geometry. )r&Zto_jsonr,)r0r2r2r3json;szOGRGeometry.jsoncCstj|jdS)z/Returns the KML representation of the Geometry.N)r&Zto_kmlr,)r0r2r2r3kmlCszOGRGeometry.kmlcCstj|jS)z#Returns the size of the WKB buffer.)r&Z get_wkbsizer,)r0r2r2r3wkb_sizeHszOGRGeometry.wkb_sizecCsftjdkrd}nd}|j}t|}tj|j|t|tj t ||S)z/Returns the WKB representation of the Geometry.littler) sys byteorderrjrr&Zto_wkbr,rrr!r )r0rnszbufr2r2r3r8Ms   zOGRGeometry.wkbcCstj|jttS)z/Returns the WKT representation of the Geometry.)r&Zto_wktr,rr)r0r2r2r3r[szOGRGeometry.wktcCs7|j}|r,|jr,d|j|jfS|jSdS)z0Returns the EWKT representation of the Geometry.z SRID=%s;%sN)rrr)r0rr2r2r3ewkt`s zOGRGeometry.ewktcCsttj|j|jS)zClones this OGR Geometry.)rr& clone_geomr,r)r0r2r2r3clonejszOGRGeometry.clonecCstj|jdS)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&Zgeom_close_ringsr,)r0r2r2r3 close_ringsnszOGRGeometry.close_ringsFcCs|r#|j}|j||St|trKtj|j|jnrt|trstj|j|jnJt|t j t j rt|}tj|j|jn t ddS)a Transforms this geometry to a different spatial reference system. May take a CoordTransform object, a SpatialReference object, string WKT or PROJ.4, and/or an integer SRID. By default nothing is returned and the geometry is transformed in-place. However, if the `clone` keyword is set, then a transformed clone of this geometry will be returned. zUTransform only accepts CoordTransform, SpatialReference, string, and integer objects.N) rs transformrrr&Zgeom_transformr,rZgeom_transform_torr_rr6)r0Z coord_transrsZkloner`r2r2r3ruws    zOGRGeometry.transformcCs1t|tstdn||j|jS)zA generalized function for topology operations, takes a GDAL function and the other geometry to perform the operation on.zrArCrErGrIrJrKpropertyrLrMrR coord_dimrSrTrUrVr.rWrXrYr[r^rarrbrcrrerfrgrhZgeojsonrirjr8rrqrsrtrurwrxrHryrzr{r|r}r~rrrrDrBrFr?r2r2r2r3r>sz 9                             rc@s^eZdZeddZeddZeddZeddZeZd S) PointcCstj|jdS)z(Returns the X coordinate for this Point.r)r&getxr,)r0r2r2r3xszPoint.xcCstj|jdS)z(Returns the Y coordinate for this Point.r)r&getyr,)r0r2r2r3yszPoint.ycCs&|jdkr"tj|jdSdS)z(Returns the Z coordinate for this Point.rOrN)rr&getzr,)r0r2r2r3zszPoint.zcCsH|jdkr|j|jfS|jdkrD|j|j|jfSdS)z Returns the tuple of this point.rNrON)rrrr)r0r2r2r3rZsz Point.tupleN) rrrrrrrrZcoordsr2r2r2r3rs rc@seZdZddZddZddZeddZeZd d Z ed d Z ed dZ eddZ dS) LineStringcCs|dkr||jkrttt}}}tj|j|t|t|t||j}|dkr|jfS|dkr|j|jfS|dkr|j|j|jfSntdt |dS)z%Returns the Point at the given index.rrlrNrOzindex out of range: %sN) rTrr&Z get_pointr,rrvaluerr+)r0indexrrrrQr2r2r3 __getitem__s.     zLineString.__getitem__ccs'x t|jD]}||VqWdS)z+Iterates over each point in the LineString.N)rrT)r0ir2r2r3__iter__szLineString.__iter__cCs|jS)z:The length returns the number of points in the LineString.)rT)r0r2r2r3__len__"szLineString.__len__cs)tfddttDS)z4Returns the tuple representation of this LineString.c3s|]}|VqdS)Nr2).0r)r0r2r3 )sz#LineString.tuple..)rZrr*)r0r2)r0r3rZ&szLineString.tuplecs&fddttDS)zp Internal routine that returns a sequence (list) corresponding with the given function. cs"g|]}j|qSr2)r,)rr)rvr0r2r3 1s z'LineString._listarr..)rr*)r0rvr2)rvr0r3_listarr,szLineString._listarrcCs|jtjS)z$Returns the X coordinates in a list.)rr&r)r0r2r2r3r3sz LineString.xcCs|jtjS)z$Returns the Y coordinates in a list.)rr&r)r0r2r2r3r8sz LineString.ycCs#|jdkr|jtjSdS)z$Returns the Z coordinates in a list.rON)rrr&r)r0r2r2r3r=sz LineString.zN) rrrrrrrrZrrrrrr2r2r2r3r s     rc@seZdZdS) LinearRingN)rrrr2r2r2r3rEs rc@seZdZddZddZddZeddZeZed d Z e Z ed d Z ed dZ dS)PolygoncCs|jS)z-The number of interior rings in this Polygon.)rS)r0r2r2r3rKszPolygon.__len__ccs'x t|jD]}||VqWdS)z*Iterates through each ring in the Polygon.N)rrS)r0rr2r2r3rOszPolygon.__iter__cCsZ|dks||jkr.td|n(ttjtj|j||jSdS)z%Gets the ring at the specified index.rzindex out of range: %sN)rSrrr&rr get_geom_refr,r)r0rr2r2r3rTszPolygon.__getitem__cCs|dS)z"Returns the shell of this Polygon.rr2)r0r2r2r3shell\sz Polygon.shellcs&tfddtjDS)z0Returns a tuple of LinearRing coordinate tuples.c3s|]}|jVqdS)N)rZ)rr)r0r2r3resz Polygon.tuple..)rZrrS)r0r2)r0r3rZbsz Polygon.tuplecs&tfddtjDS)z%The number of Points in this Polygon.c3s|]}|jVqdS)N)rT)rr)r0r2r3rlsz&Polygon.point_count..)sumrrS)r0r2)r0r3rThszPolygon.point_countcCs,ttd}tj|j|j|S)z/Returns the centroid (a Point) of this Polygon.r)rrr&Z get_centroidr,)r0pr2r2r3centroidnszPolygon.centroidN) rrrrrrrrZ exterior_ringrZrrTrr2r2r2r3rIs    rc@speZdZdZddZddZddZdd Zed d Z ed d Z e Z dS)GeometryCollectionzThe Geometry Collection class.cCsZ|dks||jkr.td|n(ttjtj|j||jSdS)z)Gets the Geometry at the specified index.rzindex out of range: %sN)rSrrr&rrrr,r)r0rr2r2r3r{szGeometryCollection.__getitem__ccs'x t|jD]}||VqWdS)zIterates over each Geometry.N)rrS)r0rr2r2r3rszGeometryCollection.__iter__cCs|jS)z5The number of geometries in this Geometry Collection.)rS)r0r2r2r3rszGeometryCollection.__len__cCst|trdt||jrKx=|D]}tj|j|jq(Wqtj|j|jnCt|tjrt|}tj|j|jn tddS)z-Add the geometry to this Geometry Collection.zMust add an OGRGeometry.N) rrr/r&Zadd_geomr,rrr )r0rr1tmpr2r2r3adds  zGeometryCollection.addcs&tfddtjDS)z1The number of Points in this Geometry Collection.c3s|]}|jVqdS)N)rT)rr)r0r2r3rsz1GeometryCollection.point_count..)rrrS)r0r2)r0r3rTszGeometryCollection.point_countcs&tfddtjDS)z;Returns a tuple representation of this Geometry Collection.c3s|]}|jVqdS)N)rZ)rr)r0r2r3rsz+GeometryCollection.tuple..)rZrrS)r0r2)r0r3rZszGeometryCollection.tupleN) rrrrrrrrrrTrZrr2r2r2r3rxs     rc@seZdZdS) MultiPointN)rrrr2r2r2r3rs rc@seZdZdS)MultiLineStringN)rrrr2r2r2r3rs rc@seZdZdS) MultiPolygonN)rrrr2r2r2r3rs rrlrNrOe)2rrmbinasciirrctypesrrrrrr Zdjango.contrib.gis.gdal.baser Z django.contrib.gis.gdal.enveloper r Zdjango.contrib.gis.gdal.errorr rrZ django.contrib.gis.gdal.geomtyperZ"django.contrib.gis.gdal.prototypesrr&rr\Zdjango.contrib.gis.gdal.srsrrZ!django.contrib.gis.geometry.regexrrrZ django.utilsrZdjango.utils.six.movesrrrrrrrrrrZwkb25bitr-r2r2r2r3(sJ .9//