API GEOS¶
Contexte¶
Qu’est-ce que GEOS ?¶
GEOS signifie Geometry Engine - Open Source. C’est une bibliothèque C++ dérivée de Java Topology Suite. GEOS implémente les prédicats et les opérateurs spatiaux de la spécification OpenGIS Simple Features for SQL. GEOS, dorénavant un projet OSGeo, a été initialement développée et maintenue par Refractions Research (Victoria, Canada).
Fonctionnalités¶
GeoDjango implémente une adaptation Python de haut niveau de la bibliothèque GEOS. Ses fonctionnalités comprennent notamment :
Une interface sous licence BSD des routines géométriques GEOS, implémentées purement en Python à l’aide de
ctypes
.Faiblement couplé à GeoDjango. Par exemple, les objets
GEOSGeometry
peuvent être utilisés en dehors d’un projet Django. En d’autres termes, il n’est pas nécessaire queDJANGO_SETTINGS_MODULE
ou une base de données soient définis.Mutabilité : les objets
GEOSGeometry
peuvent être modifiés.Multiplateforme et testé ; compatible avec les plateformes Windows, Linux, Solaris et Mac OS X.
Tutoriel¶
Cette section contient une brève introduction à l’utilisation des objets GEOSGeometry
.
Création d’un objet géométrique¶
Les objets GEOSGeometry
peuvent être créés de plusieurs façons. La première est de simplement créer une instance avec du contenu spatial. Les exemples suivants montrent comment créer le même objet géométrique avec du contenu WKT, HEX, WKB et GeoJSON :
>>> from django.contrib.gis.geos import GEOSGeometry
>>> pnt = GEOSGeometry('POINT(5 23)') # WKT
>>> pnt = GEOSGeometry('010100000000000000000014400000000000003740') # HEX
>>> pnt = GEOSGeometry(buffer('\x01\x01\x00\x00\x00\x00\x00\x00\x00\x00\x00\x14@\x00\x00\x00\x00\x00\x007@'))
>>> pnt = GEOSGeometry('{ "type": "Point", "coordinates": [ 5.000000, 23.000000 ] }') # GeoJSON
Un autre option est d’utiliser le constructeur correspondant au type d’objet géométrique que vous souhaitez créer. Par exemple un objet Point
peut être créé en passant des coordonnées X et Y à son constructeur :
>>> from django.contrib.gis.geos import Point
>>> pnt = Point(5, 23)
Tous ces constructeurs acceptent le paramètre nommé srid
. Par exemple :
>>> from django.contrib.gis.geos import GEOSGeometry, LineString, Point
>>> print(GEOSGeometry('POINT (0 0)', srid=4326))
SRID=4326;POINT (0.0000000000000000 0.0000000000000000)
>>> print(LineString((0, 0), (1, 1), srid=4326))
SRID=4326;LINESTRING (0.0000000000000000 0.0000000000000000, 1.0000000000000000 1.0000000000000000)
>>> print(Point(0, 0, srid=32140))
SRID=32140;POINT (0.0000000000000000 0.0000000000000000)
Pour terminer, il existe la méthode de fabrication fromfile()
qui renvoie un objet GEOSGeometry
à partir d’un fichier :
>>> from django.contrib.gis.geos import fromfile
>>> pnt = fromfile('/path/to/pnt.wkt')
>>> pnt = fromfile(open('/path/to/pnt.wkt'))
Les objets géométriques sont « pythonesques »¶
Les objets GEOSGeometry
sont « pythonesques », ce qui signifie qu’il est possible d’accéder, de modifier ou de passer en boucle leurs composants à l’aide des conventions Python standards. Par exemple, voici comment parcourir les coordonnées d’un Point
en boucle :
>>> pnt = Point(5, 23)
>>> [coord for coord in pnt]
[5.0, 23.0]
Avec n’importe quel objet géométrique, la propriété GEOSGeometry.coords
peut être utilisée pour obtenir les coordonnées géométriques sous forme de tuple Python :
>>> pnt.coords
(5.0, 23.0)
Il est possible d’obtenir et de définir les éléments des objets géométriques par des techniques d’indices Python standard. Cependant, les valeurs renvoyées dépendent du type d’objet géométrique. Par exemple, l’application d’un indice sur une ligne LineString
renvoie un tuple de coordonnées :
>>> from django.contrib.gis.geos import LineString
>>> line = LineString((0, 0), (0, 50), (50, 50), (50, 0), (0, 0))
>>> line[0]
(0.0, 0.0)
>>> line[-2]
(50.0, 0.0)
Alors que l’application d’un indice sur un Polygon
renvoie un anneau (un objet LinearRing
) correspondant à l’indice :
>>> from django.contrib.gis.geos import Polygon
>>> poly = Polygon( ((0.0, 0.0), (0.0, 50.0), (50.0, 50.0), (50.0, 0.0), (0.0, 0.0)) )
>>> poly[0]
<LinearRing object at 0x1044395b0>
>>> poly[0][-2] # second-to-last coordinate of external ring
(50.0, 0.0)
De plus, les coordonnées et les composants d’un objet géométrique peuvent être complétés ou modifiés de manière similaire à une liste Python :
>>> line[0] = (1.0, 1.0)
>>> line.pop()
(0.0, 0.0)
>>> line.append((1.0, 1.0))
>>> line.coords
((1.0, 1.0), (0.0, 50.0), (50.0, 50.0), (50.0, 0.0), (1.0, 1.0))
Les objets géométriques prennent en charge les opérateurs d’ensembles (set
) :
>>> from django.contrib.gis.geos import LineString
>>> ls1 = LineString((0, 0), (2, 2))
>>> ls2 = LineString((1, 1), (3, 3))
>>> print(ls1 | ls2) # equivalent to `ls1.union(ls2)`
MULTILINESTRING ((0 0, 1 1), (1 1, 2 2), (2 2, 3 3))
>>> print(ls1 & ls2) # equivalent to `ls1.intersection(ls2)`
LINESTRING (1 1, 2 2)
>>> print(ls1 - ls2) # equivalent to `ls1.difference(ls2)`
LINESTRING(0 0, 1 1)
>>> print(ls1 ^ ls2) # equivalent to `ls1.sym_difference(ls2)`
MULTILINESTRING ((0 0, 1 1), (2 2, 3 3))
L’opérateur d’égalité ne vérifie pas l’égalité spatiale
L’opérateur d’égalité de GEOSGeometry
utilise equals_exact()
, et non pas equals()
, c’est-à-dire que l’égalité dépend de coordonnées identiques et dans la même position :
>>> from django.contrib.gis.geos import LineString
>>> ls1 = LineString((0, 0), (1, 1))
>>> ls2 = LineString((1, 1), (0, 0))
>>> ls1.equals(ls2)
True
>>> ls1 == ls2
False
Objets géométriques¶
GEOSGeometry
¶
-
class
GEOSGeometry
(geo_input, srid=None)¶ Paramètres: - geo_input – Valeur d’entrée géométrique (chaîne ou tampon)
- srid (int) – identifiant de référence spatiale
Il s’agit de la classe de base de tous les objets géométriques GEOS. Elle s’initialise en fonction du paramètre geo_input
donné, puis en déduit la sous-classe géométrique appropriée (par ex. GEOSGeometry('POINT(1 1)')
crée un objet Point
).
Les formats de valeurs d’entrée suivants, ainsi que leur type Python correspondant, sont acceptés :
Format | Type en entrée |
---|---|
WKT / EWKT |
|
HEX / HEXEWKB |
|
WKB / EWKB | buffer |
GeoJSON (nécessite GDAL) |
|
Note
La nouvelle notation WKT 3D/4D avec un élément Z ou M intermédiaire (comme POINT Z (3, 4, 5)
) n’est prise en charge qu’à partir de GEOS 3.3.0.
Propriétés¶
-
GEOSGeometry.
coords
¶
Renvoie les coordonnées de l’objet géométrique sous forme de tuple.
-
GEOSGeometry.
dims
¶
Renvoie la dimension de l’objet géométrique :
0
pour lesPoint
s et lesMultiPoint
s1
pour lesLineString
s et lesMultiLineString
s2
pour lesPolygon
s et lesMultiPolygon
s-1
pour lesGeometryCollection
s videsla dimension maximale de ses éléments pour les
GeometryCollection
s non vides
-
GEOSGeometry.
empty
¶
Renvoie une valeur booléenne indiquant si l’ensemble des points de l’objet géométrique est vide.
-
GEOSGeometry.
geom_type
¶
Renvoie une chaîne correspondant au type de l’objet géométrique. Par exemple :
>>> pnt = GEOSGeometry('POINT(5 23)')
>>> pnt.geom_type
'Point'
-
GEOSGeometry.
geom_typeid
¶
Renvoie le numéro d’identification du type géométrique GEOS. Le tableau suivant expose la valeur correspondant à chaque type géométrique :
Objet géométrique |
ID |
---|---|
Point |
0 |
LineString |
1 |
LinearRing |
2 |
Polygon |
3 |
MultiPoint |
4 |
MultiLineString |
5 |
MultiPolygon |
6 |
GeometryCollection |
7 |
-
GEOSGeometry.
num_coords
¶
Renvoie le nombre de coordonnées de l’objet géométrique.
-
GEOSGeometry.
num_geom
¶
Renvoie le nombre d’objets géométriques dans cet objet. En d’autres termes, le résultat sera 1 pour toute autre chose que les collections d’objets géométriques.
-
GEOSGeometry.
hasz
¶
Renvoie une valeur booléenne indiquant si l’objet géométrique est en trois dimensions.
-
GEOSGeometry.
ring
¶
Renvoie une valeur booléenne indiquant si l’objet géométrique est un LinearRing
.
-
GEOSGeometry.
simple
¶
Renvoie une valeur booléenne indiquant si l’objet géométrique est « simple ». Un objet géométrique est simple si et seulement s’il n’a pas d’intersection avec lui-même (excepté aux points limites). Par exemple, un objet LineString
n’est pas simple si sa ligne se croise elle-même. Par définition les objets LinearRing
et Polygon
sont toujours simples car ils ne peuvent pas avoir d’intersections avec eux-mêmes.
-
GEOSGeometry.
valid
¶
Renvoie une valeur booléenne indiquant si l’objet géométrique est valide.
-
GEOSGeometry.
valid_reason
¶
Renvoie une chaîne décrivant la raison pour laquelle un objet géométrique n’est pas valable.
-
GEOSGeometry.
srid
¶
Propriété pouvant être utilisée pour récupérer ou définir le SRID associé à l’objet. Par exemple :
>>> pnt = Point(5, 23)
>>> print(pnt.srid)
None
>>> pnt.srid = 4326
>>> pnt.srid
4326
Propriétés en sortie¶
Les propriétés de cette section exportent l’objet GEOSGeometry
sous une forme différente. Cette forme peut être une chaîne, un tampon mémoire ou même un autre objet.
-
GEOSGeometry.
ewkt
¶
Renvoie le format WKT (Well-Known Text) étendu de l’objet géométrique. Cette représentation est spécifique à PostGIS et constitue un sur-ensemble du standard WKT OGC [1]. En bref, la représentation WKT est préfixée par le code SRID, par exemple SRID=4326;POINT(5 23)
.
Note
Le résultat de cette propriété n’inclut pas les informations 3dm, 3dz ou 4d que PostGIS prend en charge dans ses représentations EWKT.
-
GEOSGeometry.
hex
¶
Renvoie la représentation WKB de cet objet géométrique sous forme hexadécimale. Notez que la valeur SRID n’est pas incluse dans cette représentation car elle ne fait pas partie de la spécification OGC (utilisez plutôt la propriété GEOSGeometry.hexewkb
).
-
GEOSGeometry.
hexewkb
¶
Renvoie la représentation EWKB de cet objet géométrique sous forme hexadécimale. Il s’agit d’une extension de la spécification WKB incluant la valeur SRID correspondant à l’objet géométrique.
-
GEOSGeometry.
json
¶
Renvoie la représentation GeoJSON de l’objet géométrique. Notez que le résultat ne constitue pas une structure GeoJSON complète, mais seulement le contenu de la clé geometry
d’une structure GeoJSON. Voir aussi Sérialisation GeoJSON.
-
GEOSGeometry.
geojson
¶
Alias de GEOSGeometry.json
.
-
GEOSGeometry.
kml
¶
Renvoie une représentation KML (Keyhole Markup Language) de cet objet géométrique. Ceci ne devrait être utilisé que pour des objets géométriques ayant le SRID 4326 (WGS84), mais cette restriction n’est pas forcée.
-
GEOSGeometry.
ogr
¶
Renvoie un objet OGRGeometry
équivalent à l’objet géométrique GEOS.
Note
Nécessite GDAL.
-
GEOSGeometry.
wkb
¶
Renvoie la représentation WKB (Well-Known Binary) de cet objet géométrique sous forme de tampon Python (buffer). La valeur SRID n’est pas incluse. Si vous le souhaitez, utilisez plutôt la propriété GEOSGeometry.ewkb
.
-
GEOSGeometry.
ewkb
¶
Renvoie la représentation EWKB de cet objet géométrique sous forme de tampon Python (buffer). Il s’agit d’une extension de la spécification WKB qui inclut toute valeur SRID faisant partie de cet objet.
-
GEOSGeometry.
wkt
¶
Renvoie la représentation WKT (Well-Known Text) de cet objet géométrique (un standard OGC).
Méthodes de prédicats spatiaux¶
Toutes les méthodes de prédicats spatiaux ci-après acceptent en paramètre une autre instance GEOSGeometry
(other
) et renvoient une valeur booléenne.
-
GEOSGeometry.
contains
(other)¶
Renvoie True
si other.within(this)
renvoie True
.
-
GEOSGeometry.
crosses
(other)¶
Renvoie True
si la matrice d’intersection DE-9IM des deux objets géométriques est T*T******
(pour un point et une courbe, un point et une aire ou une ligne et une aire), ou 0********
(pour deux courbes).
-
GEOSGeometry.
disjoint
(other)¶
Renvoie True
si la matrice d’intersection DE-9IM des deux objets géométriques est FF*FF****
.
-
GEOSGeometry.
equals
(other)¶
Renvoie True
si la matrice d’intersection DE-9IM des deux objets géométriques est T*F**FFF*
.
-
GEOSGeometry.
equals_exact
(other, tolerance=0)¶
Renvoie True
si les deux objets géométriques sont exactement égaux, jusqu’à une certaine tolérance. La valeur de tolerance
doit être un nombre à virgule représentant la tolérance d’erreur dans la comparaison. Par exemple, poly1.equals_exact(poly2, 0.001)
compare l’égalité jusqu’à un millième d’unité.
-
GEOSGeometry.
intersects
(other)¶
Renvoie True
si GEOSGeometry.disjoint()
vaut False
.
-
GEOSGeometry.
overlaps
(other)¶
Renvoie True
si la matrice d’intersection DE-9IM des deux objets géométriques est T*T***T**
(pour deux points ou deux surfaces), ou 1*T***T**
(pour deux courbes).
-
GEOSGeometry.
relate_pattern
(other, pattern)¶
Renvoie True
si la matrice d’intersection DE-9IM des deux objets géométriques correspond au motif pattern
donné, une chaîne de neuf caractères dans l’alphabet {T
, F
, *
, 0
}.
-
GEOSGeometry.
touches
(other)¶
Renvoie True
si la matrice d’intersection DE-9IM des deux objets géométriques est FT*******
, F**T*****
ou F***T****
.
-
GEOSGeometry.
within
(other)¶
Renvoie True
si la matrice d’intersection DE-9IM des deux objets géométriques est T*F**F***
.
Méthodes topologiques¶
-
GEOSGeometry.
buffer
(width, quadsegs=8)¶
Renvoie un objet GEOSGeometry
représentant tous les points dont la distance à cet objet géométrique est plus petit ou égal à la longueur width
donnée. Le paramètre facultatif quadsegs
définit le nombre de segments utilisés pour faire l’approximation d’un quart de cercle (8 par défaut).
-
GEOSGeometry.
difference
(other)¶
Renvoie un objet GEOSGeometry
représentant les points composants cet objet géométrique qui ne composent pas l’autre objet.
-
GEOSGeometry.
interpolate
(distance)¶
-
GEOSGeometry.
interpolate_normalized
(distance)¶
Étant donné une distance (nombre à virgule), renvoie le point (ou le point le plus proche) de l’objet géométrique (LineString
ou MultiLineString
) à cette distance. La version normalisée accepte une distance sous forme d’un nombre à virgule entre 0 (origine) et 1 (point terminal).
Inverse de GEOSGeometry.project()
.
-
GEOSGeometry.
intersection
(other)¶
Renvoie un objet GEOSGeometry
représentant les points partagés entre cet objet géométrique et l’objet other
.
-
GEOSGeometry.
project
(point)¶
-
GEOSGeometry.
project_normalized
(point)¶
Renvoie la distance (nombre à virgule) à partir de l’origine de l’objet géométrique (LineString
ou MultiLineString
) jusqu’au point projeté sur la géométrie (c’est-à-dire le point de la ligne le plus proche du point donné). La version normalisée renvoie la distance sous forme d’un nombre à virgule entre 0 (origine) et 1 (point terminal).
Inverse de GEOSGeometry.interpolate()
.
-
GEOSGeometry.
relate
(other)¶
Renvoie la matrice d’intersection DE-9IM (une chaîne) représentant la relation topologique entre cet objet géométrique et l’objet en paramètre.
-
GEOSGeometry.
simplify
(tolerance=0.0, preserve_topology=False)¶
Renvoie un nouvel objet GEOSGeometry
, simplifié à la tolérance indiquée en utilisant l’algorithme Douglas-Peucker. Une plus haute valeur de tolérance produit moins de points dans le résultat. Si aucune tolérance n’est fournie, la valeur par défaut est 0.
Par défaut, cette fonction ne préserve pas la topologie. Par exemple, les objets Polygon
peuvent être divisés, être fusionnés pour former des lignes ou carrément disparaître. Des trous de Polygon
peuvent apparaître ou disparaître et des lignes peuvent se croiser. En définissant preserve_topology=True
, le résultat aura les mêmes dimensions et le même nombre de composants que l’objet de départ ; la fonction est cependant nettement plus lente.
-
GEOSGeometry.
sym_difference
(other)¶
Renvoie un objet GEOSGeometry
formé des points de cet objet géométrique qui ne sont pas dans l’objet en paramètre, ainsi que des points dans l’objet en paramètre qui ne sont pas dans cet objet.
-
GEOSGeometry.
union
(other)¶
Renvoie un objet GEOSGeometry
représentant tous les points de cet objet géométrique additionnés de ceux de l’objet en paramètre.
Propriétés topologiques¶
-
GEOSGeometry.
boundary
¶
Renvoie les limites sous forme d’une nouvelle instance d’objet géométrique.
-
GEOSGeometry.
centroid
¶
Renvoie un objet Point
représentant le centre géométrique de l’objet. Le point n’est pas assuré d’être à l’intérieur de l’objet géométrique.
-
GEOSGeometry.
convex_hull
¶
Renvoie le plus petit Polygon
contenant tous les points de l’objet géométrique.
-
GEOSGeometry.
envelope
¶
Renvoie un Polygon
qui représente l’enveloppe englobante de cet objet géométrique. Notez qu’un Point
peut aussi être renvoyé si l’objet en entrée est lui-même un point.
-
GEOSGeometry.
point_on_surface
¶
Calcule et renvoie un Point
assuré d’être à l’intérieur de l’objet géométrique.
Autres propriétés et méthodes¶
-
GEOSGeometry.
area
¶
Cette propriété renvoie l’aire de l’objet géométrique.
-
GEOSGeometry.
extent
¶
Cette propriété renvoie l’étendue de cet objet géométrique sous forme d’un tuple à 4 éléments, formé de (xmin, ymin, xmax, ymax)
.
-
GEOSGeometry.
clone
()¶
Cette méthode renvoie un objet GEOSGeometry
clone de l’original.
-
GEOSGeometry.
distance
(geom)¶
Renvoie la distance entre les points les plus proches de cet objet géométrique et l’objet geom
(un autre objet GEOSGeometry
).
Note
Les calculs de distance GEOS sont linéaires ; en d’autres termes, GEOS n’effectue pas de calculs sphériques même si le SRID définit un système de coordonnées géographique.
-
GEOSGeometry.
length
¶
Renvoie la longueur de cet objet géométrique (par ex. 0 pour un Point
, la longueur d’une ligne LineString
ou la circonférence d’un Polygon
).
-
GEOSGeometry.
prepared
¶
Renvoie un objet GEOS PreparedGeometry
relatif au contenu de cet objet géométrique. Les objets PreparedGeometry
sont optimisés pour les opérations contains, intersects, covers, crosses, disjoint, overlaps, touches et within. Référez-vous à la documentation Objets géométriques préparés pour plus d’informations.
-
GEOSGeometry.
srs
¶
Renvoie un objet SpatialReference
correspondant au SRID de l’objet géométrique ou .
Note
Nécessite GDAL.
-
GEOSGeometry.
transform
(ct, clone=False)¶
Transforme l’objet géométrique en fonction du paramètre de transformation de coordonnées donné (ct
), qui peut être un numéro de SRID, une chaîne de référence spatiale WKT, une chaîne PROJ.4, un objet SpatialReference
ou un objet CoordTransform
. Par défaut, l’objet géométrique est transformé lui-même et rien n’est renvoyé. Cependant, quand le paramètre clone
est défini, l’objet géométrique n’est pas modifié et c’est un clone transformé de l’objet géométrique qui est renvoyé.
Note
Exige GDAL. Génère une exception GEOSException
si GDAL n’est pas disponible ou si la valeur SRID de l’objet géométrique est None
ou plus petite que 0.
Point
¶
-
class
Point
(x, y, z=None, srid=None)¶ Les objets
Point
sont créés en utilisant des paramètres représentant les composantes de coordonnées du point ou une suite unique de coordonnées. Par exemple, les lignes suivantes sont équivalentes :>>> pnt = Point(5, 23) >>> pnt = Point([5, 23])
LineString
¶
-
class
LineString
(*args, **kwargs)¶ Les objets
LineString
sont créés en utilisant des paramètres représentant une suite de coordonnées ou une suite d’objetsPoint
. Par exemple, les lignes suivantes sont équivalentes>>> ls = LineString((0, 0), (1, 1)) >>> ls = LineString(Point(0, 0), Point(1, 1))
De plus, les objets
LineString
peuvent aussi être créés en passant une seule liste de coordonnées ou d’objetsPoint
:>>> ls = LineString( ((0, 0), (1, 1)) ) >>> ls = LineString( [Point(0, 0), Point(1, 1)] )
LinearRing
¶
-
class
LinearRing
(*args, **kwargs)¶ Les objets
LinearRing
sont construits de la même manière que les objetsLineString
. Cependant, les coordonnées doivent être fermées, c’est-à-dire que la première et la dernière coordonnée doivent être identiques. Par exemple :>>> ls = LinearRing((0, 0), (0, 1), (1, 1), (0, 0))
Remarquez que
(0, 0)
est à la fois la première et la dernière coordonnée ; si elles ne sont pas égales, une erreur est générée.
Polygon
¶
-
class
Polygon
(*args, **kwargs)¶ Les objets
Polygon
peuvent être créés en passant un ou plusieurs paramètres représentant les segments du polygone. Les paramètres doivent être soit des instances de segmentsLinearRing
, soit une suite d’éléments pouvant être utilisée pour construire un segmentLinearRing
:>>> ext_coords = ((0, 0), (0, 1), (1, 1), (1, 0), (0, 0)) >>> int_coords = ((0.4, 0.4), (0.4, 0.6), (0.6, 0.6), (0.6, 0.4), (0.4, 0.4)) >>> poly = Polygon(ext_coords, int_coords) >>> poly = Polygon(LinearRing(ext_coords), LinearRing(int_coords))
-
classmethod
from_bbox
(bbox)¶
Renvoie un objet polygone à partir du rectangle englobant donné, un tuple à 4 éléments comprenant
(xmin, ymin, xmax, ymax)
.-
num_interior_rings
¶
Renvoie le nombre de segments intérieurs de cet objet géométrique.
-
classmethod
Comparaison de polygones
Notez qu’il est possible de comparer des objets Polygon
directement avec <
or >
, mais comme la comparaison se fait à travers les segments LineString
du polygone, cela n’a pas grande signification (mais c’est cohérent et rapide). Vous pouvez toujours forcer la comparaison avec la propriété area
:
>>> if poly_1.area > poly_2.area:
>>> pass
Collections géométriques¶
MultiPoint
¶
MultiLineString
¶
-
class
MultiLineString
(*args, **kwargs)¶ Les objets
MultiLineString
peuvent être créés en passant un ou plusieurs objetsLineString
en paramètres ou une seule liste d’objetsLineString
:>>> ls1 = LineString((0, 0), (1, 1)) >>> ls2 = LineString((2, 2), (3, 3)) >>> mls = MultiLineString(ls1, ls2) >>> mls = MultiLineString([ls1, ls2])
-
merged
¶
Renvoie une ligne
LineString
représentant la fusion de tous les composants de cet objet-
MultiPolygon
¶
-
class
MultiPolygon
(*args, **kwargs)¶ Les objets
MultiPolygon
peuvent être créés en passant un ou plusieurs objetsPolygon
en paramètres ou une seule liste d’objetsPolygon
:>>> p1 = Polygon( ((0, 0), (0, 1), (1, 1), (0, 0)) ) >>> p2 = Polygon( ((1, 1), (1, 2), (2, 2), (1, 1)) ) >>> mp = MultiPolygon(p1, p2) >>> mp = MultiPolygon([p1, p2])
-
cascaded_union
¶
Renvoie un
Polygon
représentant l’union de tous les polygones composant cette collection. L’algorithme utilisé est nettement plus efficace (plus rapide) que de vouloir effectuer l’union individuelle des objets géométriques entre eux [2].-
GeometryCollection
¶
-
class
GeometryCollection
(*args, **kwargs)¶ Les objets
GeometryCollection
peuvent être créés en passant un ou plusieurs objetsGEOSGeometry
en paramètres ou une seule liste d’objetsGEOSGeometry
:>>> poly = Polygon( ((0, 0), (0, 1), (1, 1), (0, 0)) ) >>> gc = GeometryCollection(Point(0, 0), MultiPoint(Point(0, 0), Point(1, 1)), poly) >>> gc = GeometryCollection((Point(0, 0), MultiPoint(Point(0, 0), Point(1, 1)), poly))
Objets géométriques préparés¶
Pour obtenir un objet géométrique préparé, il suffit d’accéder à la propriété GEOSGeometry.prepared
. Avec une instance PreparedGeometry
, il est possible d’utiliser ses méthodes de prédicats spatiaux énumérées ci-dessous avec d’autres objets GEOSGeometry
. Une opération sur un objet géométrique préparé peut être beaucoup plus rapide ; plus l’objet géométrique est complexe, plus grand sera le gain de rapidité pour l’opération. Pour plus d’informations, consultez la page du wiki GEOS sur les objets géométriques préparés.
Par exemple :
>>> from django.contrib.gis.geos import Point, Polygon
>>> poly = Polygon.from_bbox((0, 0, 5, 5))
>>> prep_poly = poly.prepared
>>> prep_poly.contains(Point(2.5, 2.5))
True
PreparedGeometry
¶
-
class
PreparedGeometry
¶ Toutes les méthodes des objets
PreparedGeometry
acceptent un paramètreother
(autre) qui doit représenter une instanceGEOSGeometry
.-
contains
(other)¶
-
contains_properly
(other)¶
-
covers
(other)¶
-
crosses
(other)¶ Note
GEOS 3.3 est requis pour utiliser ce prédicat.
-
disjoint
(other)¶ Note
GEOS 3.3 est requis pour utiliser ce prédicat.
-
intersects
(other)¶
-
overlaps
(other)¶ Note
GEOS 3.3 est requis pour utiliser ce prédicat.
-
touches
(other)¶ Note
GEOS 3.3 est requis pour utiliser ce prédicat.
-
within
(other)¶ Note
GEOS 3.3 est requis pour utiliser ce prédicat.
-
Fabrication d’objets géométriques¶
-
fromfile
(file_h)¶ Paramètres: file_h (a Python file
object or a string path to the file) – Fichier d’entrée contenant des données spatialesType retourné: un objet GEOSGeometry
correspondant aux données spatiales du fichier
Exemple :
>>> from django.contrib.gis.geos import fromfile
>>> g = fromfile('/home/bob/geom.wkt')
-
fromstr
(string, srid=None)¶ Paramètres: Type retourné: un objet
GEOSGeometry
correspondant aux données spatiales de la chaîne
fromstr(string, srid)
est équivalente à GEOSGeometry(string, srid)
.
Exemple :
>>> from django.contrib.gis.geos import fromstr
>>> pnt = fromstr('POINT(-90.5 29.5)', srid=4326)
Objets d’entrées/sorties¶
Objets de lecture¶
Les classes d’entrées/sorties en lecture renvoient simplement une instance GEOSGeometry
à partir du contenu WKB ou WKT en entrée en fonction de leur méthode read(geom)
.
-
class
WKBReader
¶
Exemple :
>>> from django.contrib.gis.geos import WKBReader
>>> wkb_r = WKBReader()
>>> wkb_r.read('0101000000000000000000F03F000000000000F03F')
<Point object at 0x103a88910>
-
class
WKTReader
¶
Exemple :
>>> from django.contrib.gis.geos import WKTReader
>>> wkt_r = WKTReader()
>>> wkt_r.read('POINT(1 1)')
<Point object at 0x103a88b50>
Objets d’écriture¶
Tous les objets d’écriture ont une méthode write(geom)
qui renvoie la représentation WKB ou WKT de l’objet géométrique concerné. De plus, les objets WKBWriter
ont aussi des propriétés pouvant être utilisées pour changer l’ordre des octets et d’inclure la valeur SRID (autrement dit, du contenu EWKB).
-
class
WKBWriter
¶
WKBWriter
offre le plus de contrôle sur le contenu produit. Par défaut, il renvoie du contenu WKB conforme OGC lorsqu’on appelle sa méthode write
. Cependant, il possède des propriétés permettant la production de format EWKB, un sur-ensemble du standard WKB qui inclut des informations supplémentaires.
-
WKBWriter.
write
(geom)¶
Renvoie la représentation WKB de l’objet géométrique sous forme d’objet Python buffer
. Exemple :
>>> from django.contrib.gis.geos import Point, WKBWriter
>>> pnt = Point(1, 1)
>>> wkb_w = WKBWriter()
>>> wkb_w.write(pnt)
<read-only buffer for 0x103a898f0, size -1, offset 0 at 0x103a89930>
-
WKBWriter.
write_hex
(geom)¶
Renvoie la représentation WKB de l’objet géométrique sous forme hexadécimale. Exemple :
>>> from django.contrib.gis.geos import Point, WKBWriter
>>> pnt = Point(1, 1)
>>> wkb_w = WKBWriter()
>>> wkb_w.write_hex(pnt)
'0101000000000000000000F03F000000000000F03F'
-
WKBWriter.
byteorder
¶
Cette propriété peut être définie pour changer l’ordre des octets de la représentation géométrique.
Valeur |
Description |
---|---|
0 | Gros-boutiste (big endian, compatible avec les systèmes RISC) |
1 | Petit-boutiste (little endian, compatible avec les systèmes x86) |
Exemple :
>>> from django.contrib.gis.geos import Point, WKBWriter
>>> wkb_w = WKBWriter()
>>> pnt = Point(1, 1)
>>> wkb_w.write_hex(pnt)
'0101000000000000000000F03F000000000000F03F'
>>> wkb_w.byteorder = 0
'00000000013FF00000000000003FF0000000000000'
-
WKBWriter.
outdim
¶
Cette propriété peut être définie pour changer les dimensions du résultat de la représentation géométrique. En d’autres termes, si vous avez un objet 3D, définir cette valeur à 3 permet d’inclure la dimension Z dans le format WKB.
Valeur |
Description |
---|---|
2 | Par défaut, produit du contenu WKB en 2D. |
3 | Produit du contenu WKB en 3D. |
Exemple :
>>> from django.contrib.gis.geos import Point, WKBWriter
>>> wkb_w = WKBWriter()
>>> wkb_w.outdim
2
>>> pnt = Point(1, 1, 1)
>>> wkb_w.write_hex(pnt) # By default, no Z value included:
'0101000000000000000000F03F000000000000F03F'
>>> wkb_w.outdim = 3 # Tell writer to include Z values
>>> wkb_w.write_hex(pnt)
'0101000080000000000000F03F000000000000F03F000000000000F03F'
-
WKBWriter.
srid
¶
Définissez cette propriété avec une valeur booléenne pour indiquer si le code SRID de l’objet géométrique doit être inclus dans la représentation WKB. Exemple :
>>> from django.contrib.gis.geos import Point, WKBWriter
>>> wkb_w = WKBWriter()
>>> pnt = Point(1, 1, srid=4326)
>>> wkb_w.write_hex(pnt) # By default, no SRID included:
'0101000000000000000000F03F000000000000F03F'
>>> wkb_w.srid = True # Tell writer to include SRID
>>> wkb_w.write_hex(pnt)
'0101000020E6100000000000000000F03F000000000000F03F'
-
class
WKTWriter
¶
-
WKTWriter.
write
(geom)¶
Renvoie la représentation WKT de l’objet géométrique. Exemple :
>>> from django.contrib.gis.geos import Point, WKTWriter
>>> pnt = Point(1, 1)
>>> wkt_w = WKTWriter()
>>> wkt_w.write(pnt)
'POINT (1.0000000000000000 1.0000000000000000)'
Notes de bas de page
[1] | Voir PostGIS EWKB, EWKT and Canonical Forms, documentation de PostGIS au chap. 4.1.2. |
[2] | Pour plus d’informations, lisez l’article de blog de Paul Ramsey intitulé (Much) Faster Unions in PostGIS 1.4 et l’article de Martin Davis intitulé Fast polygon merging in JTS using Cascaded Union. |
Réglages¶
GEOS_LIBRARY_PATH
¶
Une chaîne indiquant l’emplacement de la bibliothèque GEOS C. Ce réglage n’est généralement nécessaire que si la bibliothèque GEOS C n’est pas à un emplacement standard (par ex. /home/bob/lib/libgeos_c.so
).
Note
Le réglage doit contenir le chemin complet vers la bibliothèque partagée en C ; en d’autres termes, il faut indiquer libgeos_c.so
, et non pas libgeos.so
.