File: //lib/python3/dist-packages/django/contrib/gis/gdal/__pycache__/geometries.cpython-310.pyc
o
�aP] � @ s� d Z ddlZddlmZ ddlmZmZmZmZm Z m
Z
ddl m Z ddl
mZmZ ddlmZmZ ddlmZ dd lmZmZ dd
lmZm Z dd l m Z m Z m Z dd l!m"Z" G d
d� de �Z#G dd� de#�Z$G dd� de#�Z%G dd� de%�Z&G dd� de#�Z'G dd� de#�Z(G dd� de(�Z)G dd � d e(�Z*G d d � d e(�Z+d e$d e%d!e'd"e)d#e*d$e+d%e(d&e&d ej, e$d ej, e%d!ej, e'd"ej, e)d#ej, e*d$ej, e+d%ej, e(iZ-dS )'a�
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.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)�b2a_hex)�byref�c_char_p�c_double�c_ubyte�c_void_p� string_at)�GDALBase)�Envelope� OGREnvelope)�
GDALException� SRSException)� OGRGeomType)�geom�srs)�CoordTransform�SpatialReference)� hex_regex�
json_regex� wkt_regex)� force_bytesc @ s� e Zd ZdZejZd|dd�Zdd� Zdd� Z e
d d
� �Z e d d � �Z
e
d
d� �Ze dd� �Ze
dd� �Zdd� Zdd� Zdd� Zdd� Zdd � Zd d � Zed d � �Zd!d"� Zd#d$� Zeee�Zed%d&� �Z ed'd(� �Z ed)d*� �Z ed+d,� �Z ed-d.� �Z ed/d0� �Z!ed1d2� �Z"ed3d4� �Z#ed5d6� �Z$ed7d8� �Z%d9d:� Z&d;d<� Z'ee&e'�Z(d=d>� Z)d?d@� Z*ee)e*�Z+dAdB� Z,edCdD� �Z-edEdF� �Z.edGdH� �Z/edIdJ� �Z0e0Z1edKdL� �Z2edMdN� �Z3edOdP� �Z4edQdR� �Z5edSdT� �Z6dUdV� Z7dWdX� Z8d}dZd[�Z9d\d]� Z:d^d_� Z;d`da� Z<dbdc� Z=ddde� Z>dfdg� Z?dhdi� Z@djdk� ZAdldm� ZBd|dndo�ZCedpdq� �ZDedrds� �ZEdtdu� ZFdvdw� ZGdxdy� ZHdzd{� ZIdS )~� OGRGeometryz Encapsulate an OGR geometry.Nc C sr t |t�}|rt�|�rtt�|��}d}|rxt�|�}t�|�}|ra|d r-t |d �}|d �
� dkrNt � t
|d �j�}t �|tt|d �� ��� nRt �tt|d �� ��dtt� ��}n?|rk| �|�� �}n5t
|� t � t
|�j�}n(t |t�r�| �|�}n t |t
�r�t � |j�}nt || j�r�|}ntdt|� ��|s�td| ��|| _|r�|| _t | j j | _ dS ) z=Initialize Geometry on either WKT or an OGR pointer as input.F�srid�type�
LINEARRING�wktNz4Invalid input type for OGR Geometry construction: %sz)Cannot create OGR Geometry from input: %s) �
isinstance�strr �match�
memoryview�bytes�fromhexr r �int�upper�capi� create_geomr �num�
import_wktr r �encode�from_wktr �
_from_json� _from_wkb�ptr_typer r �ptrr � GEO_CLASSES� geom_type� __class__)�self�
geom_inputr � str_instance�wkt_m�json_m�g� r7 �D/usr/lib/python3/dist-packages/django/contrib/gis/gdal/geometries.py�__init__? s>
&
zOGRGeometry.__init__c C s$ | j }|r |j}nd }t| j�|fS �N)r r r �wkb�r1 r r7 r7 r8 � __getstate__x s
zOGRGeometry.__getstate__c C s> |\}}t �|d tt� �t|��}|std��|| _|| _d S )Nz.Invalid OGRGeometry loaded from pickled state.)r$ �from_wkbr r �lenr r- r )r1 �stater; r r- r7 r7 r8 � __setstate__� s
zOGRGeometry.__setstate__c C s t �t|�d tt� �t|��S r: )r$ r> r r r r? )�clsr2 r7 r7 r8 r+ � s zOGRGeometry._from_wkbc C s
t �| �S r: )r$ � from_json�r2 r7 r7 r8 r* � �
zOGRGeometry._from_jsonc C s, |\}}}}t d||||||||||f
�S )z2Construct a Polygon from a bounding box (4-tuple).z,POLYGON((%s %s, %s %s, %s %s, %s %s, %s %s)))r )rB �bbox�x0�y0�x1�y1r7 r7 r8 � from_bbox� s �zOGRGeometry.from_bboxc C s t t �t| ���S r: )r r* r rD r7 r7 r8 rC � � zOGRGeometry.from_jsonc C s | t �t|���S r: )r$ �from_gmlr )rB �
gml_stringr7 r7 r8 rM � rL zOGRGeometry.from_gmlc C �
| � |�S )z'Return the union of the two geometries.)�union�r1 �otherr7 r7 r8 �__or__� rE zOGRGeometry.__or__c C rO )z7Return the intersection of this Geometry and the other.)� intersectionrQ r7 r7 r8 �__and__� rE zOGRGeometry.__and__c C rO )z2Return the difference this Geometry and the other.)�
differencerQ r7 r7 r8 �__sub__� rE zOGRGeometry.__sub__c C rO )z?Return the symmetric difference of this Geometry and the other.)�sym_differencerQ r7 r7 r8 �__xor__� rE zOGRGeometry.__xor__c C s t |t�o | �|�S )z$Is this Geometry equal to the other?)r r �equalsrQ r7 r7 r8 �__eq__� � zOGRGeometry.__eq__c C � | j S )z*WKT is used for the string representation.)r �r1 r7 r7 r8 �__str__� � zOGRGeometry.__str__c C � t �| j�S )z5Return 0 for points, 1 for lines, and 2 for surfaces.)r$ �get_dimsr- r^ r7 r7 r8 � dimension� � zOGRGeometry.dimensionc C ra )z0Return the coordinate dimension of the Geometry.)r$ �
get_coord_dimr- r^ r7 r7 r8 �_get_coord_dim� � zOGRGeometry._get_coord_dimc C s" |dvrt d��t�| j|� dS )z.Set the coordinate dimension of this Geometry.)� � z(Geometry dimension must be either 2 or 3N)�
ValueErrorr$ �
set_coord_dimr- )r1 �dimr7 r7 r8 �_set_coord_dim� s zOGRGeometry._set_coord_dimc C ra )z/Return the number of elements in this Geometry.)r$ �get_geom_countr- r^ r7 r7 r8 �
geom_count� rd zOGRGeometry.geom_countc C ra )z-Return the number of Points in this Geometry.)r$ �get_point_countr- r^ r7 r7 r8 � point_count� rd zOGRGeometry.point_countc C r] )z7Alias for `point_count` (same name method in GEOS API.)�rq r^ r7 r7 r8 �
num_points� � zOGRGeometry.num_pointsc C r] )zAlias for `point_count`.rr r^ r7 r7 r8 �
num_coords� rt zOGRGeometry.num_coordsc C s t t�| j��S )z"Return the Type for this Geometry.)r r$ �
get_geom_typer- r^ r7 r7 r8 r/ � s zOGRGeometry.geom_typec C ra )z!Return the Name of this Geometry.)r$ �
get_geom_namer- r^ r7 r7 r8 � geom_name� rd zOGRGeometry.geom_namec C ra )zHReturn the area for a LinearRing, Polygon, or MultiPolygon; 0 otherwise.)r$ �get_arear- r^ r7 r7 r8 �area� rd zOGRGeometry.areac C s t t�| jtt� ���S )z&Return the envelope for this Geometry.)r
r$ � get_enveloper- r r r^ r7 r7 r8 �envelope� s zOGRGeometry.envelopec C ra r: )r$ �is_emptyr- r^ r7 r7 r8 �empty� rg zOGRGeometry.emptyc C s | j jS )zCReturn the envelope as a 4-tuple, instead of as an Envelope object.)r| �tupler^ r7 r7 r8 �extent� � zOGRGeometry.extentc C s2 zt �| j�}tt�|��W S ty Y dS w )z/Return the Spatial Reference for this Geometry.N)r$ � get_geom_srsr- r �srs_api� clone_srsr
)r1 �srs_ptrr7 r7 r8 �_get_srs s �zOGRGeometry._get_srsc C s` t |t�r |j}n t |ttf�rt|�}|j}n|du r d}ntdt|� ��t�| j|� dS )z+Set the SpatialReference for this geometry.Nz7Cannot assign spatial reference with object of type: %s) r r r- r" r � TypeErrorr r$ �
assign_srs)r1 r r� �srr7 r7 r8 �_set_srs s
zOGRGeometry._set_srsc C s | j }|r|jS d S r: )r r r< r7 r7 r8 � _get_srid s zOGRGeometry._get_sridc C s$ t |t�s |d u r|| _d S td��)Nz!SRID must be set with an integer.)r r" r r� )r1 r r7 r7 r8 � _set_srid% s
zOGRGeometry._set_sridc C s ddl m} |�| j�S )Nr �� GEOSGeometry)�django.contrib.gis.geosr� r+ r; �r1 r� r7 r7 r8 � _geos_ptr. s zOGRGeometry._geos_ptrc C s ddl m} || �� | j�S )z3Return a GEOSGeometry object from this OGRGeometry.r r� )r� r� r� r r� r7 r7 r8 �geos2 s zOGRGeometry.geosc C ra )z.Return the GML representation of the Geometry.)r$ �to_gmlr- r^ r7 r7 r8 �gml8 rd zOGRGeometry.gmlc C s t | j��� S )z<Return the hexadecimal representation of the WKB (a string).)r r; r# r^ r7 r7 r8 �hex= � zOGRGeometry.hexc C ra )zE
Return the GeoJSON representation of this Geometry.
)r$ �to_jsonr- r^ r7 r7 r8 �jsonB s zOGRGeometry.jsonc C � t �| jd�S )z.Return the KML representation of the Geometry.N)r$ �to_kmlr- r^ r7 r7 r8 �kmlJ r� zOGRGeometry.kmlc C ra )z"Return the size of the WKB buffer.)r$ � get_wkbsizer- r^ r7 r7 r8 �wkb_sizeO rd zOGRGeometry.wkb_sizec C sF t jdkrd}nd}| j}t| � }t�| j|t|�� tt ||��S )z.Return the WKB representation of the Geometry.�little� r )
�sys� byteorderr� r r$ �to_wkbr- r r r )r1 r� �sz�bufr7 r7 r8 r; T s
zOGRGeometry.wkbc C s t �| jtt� ��S )z.Return the WKT representation of the Geometry.)r$ �to_wktr- r r r^ r7 r7 r8 r b s zOGRGeometry.wktc C s& | j }|r|jrd|j| jf S | jS )z/Return the EWKT representation of the Geometry.z
SRID=%s;%s)r r r r< r7 r7 r8 �ewktg s
zOGRGeometry.ewktc C s t t�| j�| j�S )zClone this OGR Geometry.)r r$ �
clone_geomr- r r^ r7 r7 r8 �cloneq r\ zOGRGeometry.clonec C s t �| j� dS )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_ringsr- r^ r7 r7 r8 � close_ringsu s zOGRGeometry.close_ringsFc C s� |r
| � � }|�|� |S t|t�r t�| j|j� dS t|t�r+t�| j|j� dS t|t t
f�r@t|�}t�| j|j� dS t d��)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 r r$ �geom_transformr- r �geom_transform_tor" r r� )r1 � coord_transr� �kloner� r7 r7 r8 r� ~ s
zOGRGeometry.transformc C s t |t�s td��|| j|j�S )z�A generalized function for topology operations, takes a GDAL function and
the other geometry to perform the operation on.z<Must use another OGRGeometry object for topology operations!)r r r� r- )r1 �funcrR r7 r7 r8 � _topology� s
zOGRGeometry._topologyc C � | � tj|�S )z7Return True if this geometry intersects with the other.)r� r$ �ogr_intersectsrQ r7 r7 r8 �
intersects� � zOGRGeometry.intersectsc C r� )z8Return True if this geometry is equivalent to the other.)r� r$ �
ogr_equalsrQ r7 r7 r8 rZ � r� zOGRGeometry.equalsc C r� )zBReturn True if this geometry and the other are spatially disjoint.)r� r$ � ogr_disjointrQ r7 r7 r8 �disjoint� r� zOGRGeometry.disjointc C r� )z/Return True if this geometry touches the other.)r� r$ � ogr_touchesrQ r7 r7 r8 �touches� r� zOGRGeometry.touchesc C r� )z/Return True if this geometry crosses the other.)r� r$ � ogr_crossesrQ r7 r7 r8 �crosses� r� zOGRGeometry.crossesc C r� )z1Return True if this geometry is within the other.)r� r$ �
ogr_withinrQ r7 r7 r8 �within� r� zOGRGeometry.withinc C r� )z0Return True if this geometry contains the other.)r� r$ � ogr_containsrQ r7 r7 r8 �contains� r� zOGRGeometry.containsc C r� )z0Return True if this geometry overlaps the other.)r� r$ � ogr_overlapsrQ r7 r7 r8 �overlaps� r� zOGRGeometry.overlapsc C s2 t |t�rt|| j|j�| j�S t|| j�| j�S )z?A helper routine for the OGR routines that generate geometries.)r r r- r )r1 �gen_funcrR r7 r7 r8 �_geomgen� s
zOGRGeometry._geomgenc C � | � tj�S )z%Return the boundary of this geometry.)r� r$ � get_boundaryr^ r7 r7 r8 �boundary� rd zOGRGeometry.boundaryc C r� )zk
Return the smallest convex Polygon that contains all the points in
this Geometry.
)r� r$ �geom_convex_hullr^ r7 r7 r8 � convex_hull� s zOGRGeometry.convex_hullc C r� )z�
Return a new geometry consisting of the region which is the difference
of this geometry and the other.
)r� r$ � geom_diffrQ r7 r7 r8 rV � � zOGRGeometry.differencec C r� )zx
Return a new geometry consisting of the region of intersection of this
geometry and the other.
)r� r$ �geom_intersectionrQ r7 r7 r8 rT � r� zOGRGeometry.intersectionc C r� )zq
Return a new geometry which is the symmetric difference of this
geometry and the other.
)r� r$ �
geom_sym_diffrQ r7 r7 r8 rX � r� zOGRGeometry.sym_differencec C r� )z{
Return a new geometry consisting of the region which is the union of
this geometry and the other.
)r� r$ �
geom_unionrQ r7 r7 r8 rP � r� zOGRGeometry.unionr: )F)J�__name__�
__module__� __qualname__�__doc__r$ � destroy_geom�
destructorr9 r= rA � classmethodr+ � staticmethodr* rK rC rM rS rU rW rY r[ r_ �propertyrc rf rm � coord_dimro rq rs ru r/ rx rz r| r~ r� r� r� r r� r� r r� r� r� r� r� �geojsonr� r� r; r r� r� r� r� r� r� rZ r� r� r� r� r� r� r� r� r� rV rT rX rP r7 r7 r7 r8 r ; s�
9
r c s\ e Zd Z� fdd�Zedd� �Zedd� �Zedd� �Zed d
� �Z ed d � �Z
e
Z � Z S )
�Pointc s& ddl m} | jr|j�� S t� �� S )Nr )r� )�django.contrib.gisr� r~ r� �
_create_empty�superr� )r1 r� �r0 r7 r8 r� � s zPoint._geos_ptrc C s t �td�j�S )N�point)r$ r% r r& )rB r7 r7 r8 r� � s zPoint._create_emptyc C r� )z'Return the X coordinate for this Point.r )r$ �getxr- r^ r7 r7 r8 �x r� zPoint.xc C r� )z'Return the Y coordinate for this Point.r )r$ �getyr- r^ r7 r7 r8 �y r� zPoint.yc C s | j dkr t�| jd�S dS )z'Return the Z coordinate for this Point.ri r N)r� r$ �getzr- r^ r7 r7 r8 �z
s
�zPoint.zc C s4 | j dkr | j| jfS | j dkr| j| j| jfS dS )z Return the tuple of this point.rh ri N)r� r� r� r� r^ r7 r7 r8 r s
�z Point.tuple)
r� r� r� r� r� r� r� r� r� r� r �coords�
__classcell__r7 r7 r� r8 r� � s
r� c @ sX e Zd Zdd� Zdd� Zedd� �ZeZdd� Zed d
� �Z ed d � �Z
ed
d� �Z dS )�
LineStringc C s� d| kr | j k rKn n>t� t� t� }}}t�| j|t|�t|�t|�� | j}|dkr3|jfS |dkr=|j|jfS |dkrI|j|j|jfS dS td| ��)z$Return the Point at the given index.r r� rh ri z>Index out of range when accessing points of a line string: %s.N) rq r r$ � get_pointr- r r� �value�
IndexError)r1 �indexr� r� r� rl r7 r7 r8 � __getitem__ s � zLineString.__getitem__c C r] )z.Return the number of points in the LineString.rr r^ r7 r7 r8 �__len__+ r` zLineString.__len__c s t � fdd�tt� ��D ��S )z3Return the tuple representation of this LineString.c 3 s � | ]}� | V qd S r: r7 ��.0�ir^ r7 r8 � <genexpr>2 s � z#LineString.tuple.<locals>.<genexpr>)r �ranger? r^ r7 r^ r8 r / s zLineString.tuplec s � �fdd�t t���D �S )zp
Internal routine that returns a sequence (list) corresponding with
the given function.
c s g | ]}� �j |��qS r7 )r- r� �r� r1 r7 r8 �
<listcomp>: s z'LineString._listarr.<locals>.<listcomp>)r� r? )r1 r� r7 r� r8 �_listarr5 s zLineString._listarrc C r� )z#Return the X coordinates in a list.)r� r$ r� r^ r7 r7 r8 r� < rd z LineString.xc C r� )z#Return the Y coordinates in a list.)r� r$ r� r^ r7 r7 r8 r� A rd z LineString.yc C s | j dkr | �tj�S dS )z#Return the Z coordinates in a list.ri N)r� r� r$ r� r^ r7 r7 r8 r� F s
�z LineString.zN) r� r� r� r� r� r� r r� r� r� r� r� r7 r7 r7 r8 r� s
r� c @ � e Zd ZdS )�
LinearRingN�r� r� r� r7 r7 r7 r8 r� N � r� c @ sT e Zd Zdd� Zdd� Zedd� �ZeZedd� �ZeZ ed d
� �Z
ed d � �Z d
S )�Polygonc C r] )z4Return the number of interior rings in this Polygon.�ro r^ r7 r7 r8 r� T r` zPolygon.__len__c C �B d| kr | j k rn ntt�t�| j|��| j�S td| ��)z$Get the ring at the specified index.r z9Index out of range when accessing rings of a polygon: %s.�ro r r$ r� � get_geom_refr- r r� �r1 r� r7 r7 r8 r� X � zPolygon.__getitem__c C s | d S )z!Return the shell of this Polygon.r r7 r^ r7 r7 r8 �shell` r� z
Polygon.shellc � t � fdd�t� j�D ��S )z/Return a tuple of LinearRing coordinate tuples.c 3 � � | ]}� | j V qd S r: �r r� r^ r7 r8 r� i � � z Polygon.tuple.<locals>.<genexpr>�r r� ro r^ r7 r^ r8 r f � z
Polygon.tuplec r )z,Return the number of Points in this Polygon.c 3 r r: rr r� r^ r7 r8 r� p r z&Polygon.point_count.<locals>.<genexpr>��sumr� ro r^ r7 r^ r8 rq l � zPolygon.point_countc C s t td��}t�| j|j� |S )z.Return the centroid (a Point) of this Polygon.r� )r r r$ � get_centroidr- )r1 �pr7 r7 r8 �centroidr s zPolygon.centroidN) r� r� r� r� r� r� r �
exterior_ringr r� rq r r7 r7 r7 r8 r� R s
r� c @ sD e Zd ZdZdd� Zdd� Zdd� Zedd � �Zed
d � �Z e Z
d S )
�GeometryCollectionz The Geometry Collection class.c C r� )z(Get the Geometry at the specified index.r z?Index out of range when accessing geometry in a collection: %s.r� r r7 r7 r8 r� r z GeometryCollection.__getitem__c C r] )z<Return the number of geometries in this Geometry Collection.r� r^ r7 r7 r8 r� � r` zGeometryCollection.__len__c C sv t |t�r$t || j�r|D ]
}t�| j|j� q
dS t�| j|j� dS t |t�r7t|�}t�| j|j� dS td��)z-Add the geometry to this Geometry Collection.zMust add an OGRGeometry.N)r r r0 r$ �add_geomr- r r )r1 r r6 �tmpr7 r7 r8 �add� s
�
zGeometryCollection.addc r )z8Return the number of Points in this Geometry Collection.c 3 r r: rr r� r^ r7 r8 r� � r z1GeometryCollection.point_count.<locals>.<genexpr>r r^ r7 r^ r8 rq � r z GeometryCollection.point_countc r )z:Return a tuple representation of this Geometry Collection.c 3 r r: r r� r^ r7 r8 r� � r z+GeometryCollection.tuple.<locals>.<genexpr>r r^ r7 r^ r8 r � r zGeometryCollection.tupleN) r� r� r� r� r� r� r r� rq r r� r7 r7 r7 r8 r | s
r c @ r� )�
MultiPointNr� r7 r7 r7 r8 r � r� r c @ r� )�MultiLineStringNr� r7 r7 r7 r8 r � r� r c @ r� )� MultiPolygonNr� r7 r7 r7 r8 r � r� r r� rh ri � � � � �e ).r� r� �binasciir �ctypesr r r r r r � 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.prototypesr r$ r r� �django.contrib.gis.gdal.srsr r �django.contrib.gis.geometryr r r �django.utils.encodingr r r� r� r� r� r r r r �wkb25bitr. r7 r7 r7 r8 �<module> sP ( >$4**
�