- en
- Langue : fr
API de base de données GeoDjango¶
Moteurs de base de données spatiales¶
GeoDjango propose actuellement les moteurs de base de données spatiales suivants :
- django.contrib.gis.db.backends.postgis
- django.contrib.gis.db.backends.mysql
- django.contrib.gis.db.backends.oracle
- django.contrib.gis.db.backends.spatialite
Limitations spatiales de MySQL¶
Les extensions spatiales de MySQL ne prennent en charge que les opérations de rectangle englobant (ce que MySQL appelle « minimum bounding rectangle » ou MBR). Spécifiquement, MySQL ne se conforme pas au standard OGC :
Actuellement, MySQL n’implémente pas ces fonctions [Contains, Crosses, Disjoint, Intersects, Overlaps, Touches, Within] en accord avec la spécification. Celles qui sont implémentées renvoient le même résultat que les fonctions correspondantes qui agissent sur le MBR.
En d’autres termes, même si les requêtes spatiales telles que contains sont disponibles dans GeoDjango en utilisant MySQL, les résultats renvoyés sont tout à fait équivalentes à ce que renverrait bbcontains avec un autre moteur spatial.
Avertissement
Les vrais index spatiaux (R-trees) ne sont pris en charge qu’avec les tables MyISAM de MySQL. [5] En d’autres termes, lorsque vous utilisez les extensions spatiales de MySQL, vous devez choisir entre des requêtes spatiales rapides et l’intégrité de vos données (les tables MyISAM ne gérant pas les transactions ni les contraintes de clé étrangère).
Création et enregistrement de modèles géographiques¶
Voici un exemple de la manière dont il faut créer un objet géométrique (en supposant que le modèle Zipcode existe) :
>>> from zipcode.models import Zipcode
>>> z = Zipcode(code=77096, poly='POLYGON(( 10 10, 10 20, 20 20, 20 15, 10 10))')
>>> z.save()
Les objets GEOSGeometry peuvent aussi être utilisés pour enregistrer des modèles géométriques :
>>> from django.contrib.gis.geos import GEOSGeometry
>>> poly = GEOSGeometry('POLYGON(( 10 10, 10 20, 20 20, 20 15, 10 10))')
>>> z = Zipcode(code=77096, poly=poly)
>>> z.save()
De plus, si l’objet GEOSGeometry est dans un autre système de coordonnées (sa valeur SRID est différente) que celui du champ, il sera implicitement converti dans le SRID du champ du modèle, en utilisant la procédure de transformation de la base de données spatiale :
>>> poly_3084 = GEOSGeometry('POLYGON(( 10 10, 10 20, 20 20, 20 15, 10 10))', srid=3084) # SRID 3084 is 'NAD83(HARN) / Texas Centric Lambert Conformal'
>>> z = Zipcode(code=78212, poly=poly_3084)
>>> z.save()
>>> from django.db import connection
>>> print(connection.queries[-1]['sql']) # printing the last SQL statement executed (requires DEBUG=True)
INSERT INTO "geoapp_zipcode" ("code", "poly") VALUES (78212, ST_Transform(ST_GeomFromWKB('\\001 ... ', 3084), 4326))
Ainsi, les paramètres de géométrie peuvent être transmis en utilisant un objet GEOSGeometry, une syntaxe WKT (Well Known Text [1]), HEXEWKB (spécifique à PostGIS, un objet géométrique WKB en hexadécimal [2]), et GeoJSON [3] (nécessite GDAL). En fait, si la donnée n’est pas un objet GEOSGeometry, le champ géométrique essaie de créer une instance GEOSGeometry à partir des données fournies.
Pour plus d’informations sur la création d’objets GEOSGeometry, consultez le tutoriel GEOSl.
Recherches spatiales¶
Les types de recherches GeoDjango peuvent être utilisés avec toute méthode de gestionnaire comme filter(), exclude(), etc. Cependant, les types de recherches spécifiques à GeoDjango ne sont disponibles que pour les champs géométriques. Les filtres sur les champs « normaux » (par ex. CharField) peuvent être enchaînés avec d’autres sur les champs géographiques. Ainsi, les requêtes géographiques ont la forme générale suivante (en prenant comme exemple le modèle Zipcode utilisé dans l’API de modèle de GeoDjango) :
>>> qs = Zipcode.objects.filter(<field>__<lookup_type>=<parameter>)
>>> qs = Zipcode.objects.exclude(...)
Par exemple :
>>> qs = Zipcode.objects.filter(poly__contains=pnt)
Dans ce cas, poly est le champ géographique, contains est le type de recherche spatiale et pnt est le paramètre (qui pourrait être un objet GEOSGeometry ou une chaîne GeoJSON , WKT ou HEXEWKB).
Une référence complète peut être consultée dans la référence sur les recherches spatiales.
Note
GeoDjango construit le code SQL spatial à l’aide des objets GeoQuerySet, une sous-classe de QuerySet. L’instance GeoManager liée à votre modèle est ce qui permet d’utiliser GeoQuerySet.
Requêtes de distance¶
Introduction¶
Les calculs de distance sur des données spatiales sont complexes, parce que malheureusement, la Terre n’est pas plate. Certaines requêtes de distance avec des champs dans un système de coordonnées géographiques doivent parfois être exprimées différemment en raison de limites dans PostGIS. Pour plus de détails, consultez la section Choix d’un SRID dans la documentation API de modèle de GeoDjango.
Recherches de distance¶
Disponibilité : PostGIS, Oracle, SpatiaLite
Les recherches de distance suivantes sont disponibles :
Note
Pour des mensurations, au lieu d’interroger sur les distances, utilisez la méthode GeoQuerySet.distance().
Les recherches de distance acceptent un paramètre sous forme de tuple comprenant :
Un objet géométrique sur lequel baser les calculs ; et
Un nombre ou un objet Distance contenant la distance.
Si un objet Distance est utilisé, il peut être exprimé dans n’importe quelle unité (le code SQL généré utilisera des unités converties dans celles du champ) ; sinon, les paramètres numériques sont supposés être dans la même unité que celle du champ.
Note
Dans PostGIS, ST_Distance_Sphere ne restreint plus les types géométriques sur lesquels les requêtes de distance géographiques sont appliquées [4]. Cependant, ces requêtes peuvent être très lentes, car les distances de grand cercle doivent être calculées à la volée pour chaque ligne de la requête. La raison en est que l’index spatial des champs géométriques traditionnels ne peut pas être exploité.
Pour de bien meilleures performances concernant les requêtes de distance WGS84, évaluez la possibilité d’utiliser plutôt des colonnes de type géographique dans votre base de données car leur index spatial peut être exploité dans les requêtes de distance. Vous pouvez indiquer à GeoDjango de créer des colonnes géographiques en définissant geography=True dans la définition du champ.
Par exemple, si l’on dispose d’un modèle SouthTexasCity (exemple tiré des tests de distance de GeoDjango) avec un système de coordonnées projeté valide pour les villes du sud du Texas :
from django.contrib.gis.db import models
class SouthTexasCity(models.Model):
name = models.CharField(max_length=30)
# A projected coordinate system (only valid for South Texas!)
# is used, units are in meters.
point = models.PointField(srid=32140)
objects = models.GeoManager()
Les requêtes de distance suivantes pourraient être effectuées :
>>> from django.contrib.gis.geos import GEOSGeometry
>>> from django.contrib.gis.measure import D # ``D`` is a shortcut for ``Distance``
>>> from geoapp.models import SouthTexasCity
# Distances will be calculated from this point, which does not have to be projected.
>>> pnt = GEOSGeometry('POINT(-96.876369 29.905320)', srid=4326)
# If numeric parameter, units of field (meters in this case) are assumed.
>>> qs = SouthTexasCity.objects.filter(point__distance_lte=(pnt, 7000))
# Find all Cities within 7 km, > 20 miles away, and > 100 chains away (an obscure unit)
>>> qs = SouthTexasCity.objects.filter(point__distance_lte=(pnt, D(km=7)))
>>> qs = SouthTexasCity.objects.filter(point__distance_gte=(pnt, D(mi=20)))
>>> qs = SouthTexasCity.objects.filter(point__distance_gte=(pnt, D(chain=100)))
Tableaux de compatibilité¶
Recherches spatiales¶
Le tableau suivant indique en résumé quelles sont les recherches spatiales disponibles pour chaque moteur de base de données spatiale.
Type de recherche |
PostGIS | Oracle | MySQL [6] | SpatiaLite |
---|---|---|---|---|
bbcontains | X | X | X | |
bboverlaps | X | X | X | |
contained | X | X | X | |
contains | X | X | X | X |
contains_properly | X | |||
coveredby | X | X | ||
covers | X | X | ||
crosses | X | X | ||
disjoint | X | X | X | X |
distance_gt | X | X | X | |
distance_gte | X | X | X | |
distance_lt | X | X | X | |
distance_lte | X | X | X | |
dwithin | X | X | ||
equals | X | X | X | X |
exact | X | X | X | X |
intersects | X | X | X | X |
overlaps | X | X | X | X |
relate | X | X | X | |
same_as | X | X | X | X |
touches | X | X | X | X |
within | X | X | X | X |
left | X | |||
right | X | |||
overlaps_left | X | |||
overlaps_right | X | |||
overlaps_above | X | |||
overlaps_below | X | |||
strictly_above | X | |||
strictly_below | X |
Méthodes GeoQuerySet¶
Le tableau suivant indique en résumé quelles méthodes GeoQuerySet sont disponibles pour chae moteur de base de données spatiale. Notez que MySQL n’en prend en charge aucune et qu’il est donc exclu du tableau.
Méthode |
PostGIS | Oracle | SpatiaLite |
---|---|---|---|
GeoQuerySet.area() | X | X | X |
GeoQuerySet.centroid() | X | X | X |
GeoQuerySet.difference() | X | X | X |
GeoQuerySet.distance() | X | X | X |
GeoQuerySet.envelope() | X | X | |
GeoQuerySet.force_rhr() | X | ||
GeoQuerySet.geohash() | X | ||
GeoQuerySet.geojson() | X | X | |
GeoQuerySet.gml() | X | X | X |
GeoQuerySet.intersection() | X | X | X |
GeoQuerySet.kml() | X | X | |
GeoQuerySet.length() | X | X | X |
GeoQuerySet.mem_size() | X | ||
GeoQuerySet.num_geom() | X | X | X |
GeoQuerySet.num_points() | X | X | X |
GeoQuerySet.perimeter() | X | X | |
GeoQuerySet.point_on_surface() | X | X | X |
GeoQuerySet.reverse_geom() | X | X | |
GeoQuerySet.scale() | X | X | |
GeoQuerySet.snap_to_grid() | X | ||
GeoQuerySet.svg() | X | X | |
GeoQuerySet.sym_difference() | X | X | X |
GeoQuerySet.transform() | X | X | X |
GeoQuerySet.translate() | X | X | |
GeoQuerySet.union() | X | X | X |
Fonctions d’agrégation¶
Le tableau suivant indique en résumé quelles fonctions d’agrégation géographiques sont disponibles pour chaque moteur de base de données spatiale. Notez que MySQL n’en prend en charge aucune et qu’il est donc exclu du tableau.
Méthode d’agrégation |
PostGIS | Oracle | SpatiaLite |
---|---|---|---|
Collect | X | (à partir de la v3.0) |
|
Extent | X | X | (à partir de la v3.0) |
Extent3D | X | ||
MakeLine | X | ||
Union | X | X | X |
Notes de bas de page
[1] | Voir Open Geospatial Consortium, Inc., OpenGIS Simple Feature Specification For SQL, Document 99-049 (5 mai 1999), au Ch. 3.2.5, p. 3-11 (SQL Textual Representation of Geometry). |
[2] | Voir PostGIS EWKB, EWKT and Canonical Forms, documentation de PostGIS au chap. 4.1.2. |
[3] | Voir Howard Butler, Martin Daly, Allan Doyle, Tim Schaub, & Christopher Schmidt, The GeoJSON Format Specification, révision 1.0 (16 juin 2008). |
[4] | Voir la documentation de PostGIS pour ST_distance_sphere. |
[5] | Voir Création d’index spatiaux dans le manuel de référence de MySQL :
|
[6] | Référez-vous à la section Limitations spatiales de MySQL pour plus de détails. |