GeoDjango propose actuellement les moteurs de base de données spatiales suivants :
\ndjango.contrib.gis.db.backends.postgis
django.contrib.gis.db.backends.mysql
django.contrib.gis.db.backends.oracle
django.contrib.gis.db.backends.spatialite
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 :
\n\n\nActuellement, MySQL n’implémente pas ces fonctions [
\nContains,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.
RasterField n’est actuellement implémenté que pour le moteur PostGIS. Les requêtes spatiales sont disponibles pour les champs matriciels, mais pas les fonctions de base de données spatiales ni les agrégats.
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\n>>> z = Zipcode(code=77096, poly='POLYGON(( 10 10, 10 20, 20 20, 20 15, 10 10))')\n>>> z.save()\nLes objets GEOSGeometry peuvent aussi être utilisés pour enregistrer des modèles géométriques :
>>> from django.contrib.gis.geos import GEOSGeometry\n>>> poly = GEOSGeometry('POLYGON(( 10 10, 10 20, 20 20, 20 15, 10 10))')\n>>> z = Zipcode(code=77096, poly=poly)\n>>> z.save()\nDe 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'\n>>> z = Zipcode(code=78212, poly=poly_3084)\n>>> z.save()\n>>> from django.db import connection\n>>> print(connection.queries[-1]['sql']) # printing the last SQL statement executed (requires DEBUG=True)\nINSERT INTO "geoapp_zipcode" ("code", "poly") VALUES (78212, ST_Transform(ST_GeomFromWKB('\\\\001 ... ', 3084), 4326))\nAinsi, 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]. 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.
Lors de la création de modèles matriciels, le champ matriciel convertira implicitement la donnée d’entrée en un objet GDALRaster par évaluation différée. Ce champ accepte donc les mêmes données que le constructeur GDALRaster.
Voici un exemple de la manière dont on peut créer un objet matriciel à partir d’un fichier matriciel volcano.tif (en supposant que le modèle Elevation existe) :
>>> from elevation.models import Elevation\n>>> dem = Elevation(name='Volcano', rast='/path/to/raster/volcano.tif')\n>>> dem.save()\nLes objets GDALRaster peuvent aussi être utilisés pour enregistrer des modèles matriciels :
>>> from django.contrib.gis.gdal import GDALRaster\n>>> rast = GDALRaster({'width': 10, 'height': 10, 'name': 'Canyon', 'srid': 4326,\n... 'scale': [0.1, -0.1], 'bands': [{"data": range(100)}]})\n>>> dem = Elevation(name='Canyon', rast=rast)\n>>> dem.save()\nNotez que c’est l’équivalent de :
\n>>> dem = Elevation.objects.create(\n... name='Canyon',\n... rast={'width': 10, 'height': 10, 'name': 'Canyon', 'srid': 4326,\n... 'scale': [0.1, -0.1], 'bands': [{"data": range(100)}]},\n... )\nLes types de requêtes de GeoDjango peuvent être utilisés avec n’importe quelle méthode de gestionnaire comme filter(), exclude(), etc. Cependant, les types de requêtes spécifiques à GeoDjango ne sont disponibles que pour les champs spatiaux.
Les filtres pour les champs « normaux » (par ex. CharField) peuvent s’enchaîner avec des filtres sur des champs géographiques. Les requêtes géographiques acceptent des entrées géométriques et matricielles des deux côtés et les types d’entrées peuvent être librement mélangés.
La structure générale des requêtes géographiques est décrite ci-dessous. Une référence complète peut être consultée dans la référence sur les recherches spatiales.
\nLes requêtes géographiques avec des objets géométriques prennent la forme générale suivante (modèle Zipcode d’exemple utilisé dans API de modèle de GeoDjango) :
>>> qs = Zipcode.objects.filter(<field>__<lookup_type>=<parameter>)\n>>> qs = Zipcode.objects.exclude(...)\nPar exemple :
\n>>> qs = Zipcode.objects.filter(poly__contains=pnt)\n>>> qs = Elevation.objects.filter(poly__contains=rst)\nDans ce cas, poly est le champ géographique, contains est le type de recherche spatiale, pnt est le paramètre (qui pourrait être un objet GEOSGeometry ou une chaîne GeoJSON , WKT ou HEXEWKB) et rst est un objet GDALRaster.
La syntaxe des requêtes matricielles est semblable à celle pour les géométries. La seule différence est qu’un indice de bande peut être indiqué en supplément. Si aucun indice de bande n’est indiqué, c’est la première bande (indice 0) qui est utilisée par défaut. Dans ce cas, la syntaxe est identique à la syntaxe des requêtes géométriques.
Pour indiquer l’indice de bande, un paramètre supplémentaire peut être ajouté des deux côtés de la requête. Du côté gauche, c’est la syntaxe de double soulignement qui est utilisée pour transmettre l’indice de bande. Du côté droit, un tuple contenant la matrice et l’indice de bande peut être utilisé.
\nCela aboutit à la forme générale suivante pour les requêtes impliquant des matrices (modèle Elevation d’exemple utilisé dans API de modèle de GeoDjango)
>>> qs = Elevation.objects.filter(<field>__<lookup_type>=<parameter>)\n>>> qs = Elevation.objects.filter(<field>__<band_index>__<lookup_type>=<parameter>)\n>>> qs = Elevation.objects.filter(<field>__<lookup_type>=(<raster_input, <band_index>)\nPar exemple :
\n>>> qs = Elevation.objects.filter(rast__contains=geom)\n>>> qs = Elevation.objects.filter(rast__contains=rst)\n>>> qs = Elevation.objects.filter(rast__1__contains=geom)\n>>> qs = Elevation.objects.filter(rast__contains=(rst, 1))\n>>> qs = Elevation.objects.filter(rast__1__contains=(rst, 1))\nDu côté gauche de l’exemple, rast est le champ matriciel géographique et contains est le type de requête spatiale. Du côté droit, geom est l’entrée géométrique et rst est un objet GDALRaster. L’indice de bande vaut 0 par défaut dans les deux premières requêtes et est défini à 1 dans les suivantes.
Bien que toutes les requêtes spatiales puissent être utilisées avec des objets matriciels des deux côtés, certains opérateurs ne savent pas gérer des entrées matricielles. Pour les cas où l’opérateur s’attend à une entrée géométrique, la matrice est automatiquement convertie en géométrie. Il est important d’avoir cela à l’esprit lors de l’interprétation des résultats de requête.
\nLe type de prise en charge matricielle est énuméré pour chaque requête dans la table de compatibilité. Les requêtes impliquant des matrices ne sont actuellement disponibles que pour le moteur PostGIS.
\nLes 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.
\nDisponibilité : PostGIS, Oracle, SpatiaLite, PGRaster (natif)
\nLes recherches de distance suivantes sont disponibles :
\n\n\nLes recherches de distance acceptent un paramètre sous forme de tuple comprenant :
\nUn objet géométrique ou matriciel 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.
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\n\nclass SouthTexasCity(models.Model):\n name = models.CharField(max_length=30)\n # A projected coordinate system (only valid for South Texas!)\n # is used, units are in meters.\n point = models.PointField(srid=32140)\nLes requêtes de distance suivantes pourraient être effectuées :
\n>>> from django.contrib.gis.geos import GEOSGeometry\n>>> from django.contrib.gis.measure import D # ``D`` is a shortcut for ``Distance``\n>>> from geoapp.models import SouthTexasCity\n# Distances will be calculated from this point, which does not have to be projected.\n>>> pnt = GEOSGeometry('POINT(-96.876369 29.905320)', srid=4326)\n# If numeric parameter, units of field (meters in this case) are assumed.\n>>> qs = SouthTexasCity.objects.filter(point__distance_lte=(pnt, 7000))\n# Find all Cities within 7 km, > 20 miles away, and > 100 chains away (an obscure unit)\n>>> qs = SouthTexasCity.objects.filter(point__distance_lte=(pnt, D(km=7)))\n>>> qs = SouthTexasCity.objects.filter(point__distance_gte=(pnt, D(mi=20)))\n>>> qs = SouthTexasCity.objects.filter(point__distance_gte=(pnt, D(chain=100)))\nLes requêtes matricielles fonctionnent de la même manière en remplaçant simplement le champ géométrique point par un champ marticiel, ou l’objet pnt par un objet matriciel, ou les deux. Pour indiquer l’indice de bande d’une entrée matricielle dans le côté droit, un tuple ternaire peut être passé à la requête comme ceci :
>>> qs = SouthTexasCity.objects.filter(point__distance_gte=(rst, 2, D(km=7)))\nOù la bande avec indice 2 (troisième bande) de la matrice rst serait utilisée pour la requête.
Le tableau suivant fournit un résumé des requêtes spatiales disponibles pour chaque moteur de base de données spatial. Les requêtes PostGIS Raster (PGRaster) sont divisées dans les trois catégories décrites dans les détails de requêtes matricielles: prise en charge native N, prise en charge native bilatérale B et prise en charge par conversion géométrique C.
Type de recherche | \nPostGIS | \nOracle | \nMySQL [6] | \nSpatiaLite | \nPGRaster | \n
|---|---|---|---|---|---|
| \n | X | \n\n | X | \nX | \nN | \n
| \n | X | \n\n | X | \nX | \nN | \n
| \n | X | \n\n | X | \nX | \nN | \n
| \n | X | \nX | \nX | \nX | \nB | \n
| \n | X | \n\n | \n | \n | B | \n
| \n | X | \nX | \n\n | \n | B | \n
| \n | X | \nX | \n\n | \n | B | \n
| \n | X | \n\n | \n | X | \nC | \n
| \n | X | \nX | \nX | \nX | \nB | \n
| \n | X | \nX | \n\n | X | \nN | \n
| \n | X | \nX | \n\n | X | \nN | \n
| \n | X | \nX | \n\n | X | \nN | \n
| \n | X | \nX | \n\n | X | \nN | \n
| \n | X | \nX | \n\n | X | \nB | \n
| \n | X | \nX | \nX | \nX | \nC | \n
| \n | X | \nX | \nX | \nX | \nB | \n
| \n | X | \nX | \nX | \nX | \nB | \n
| \n | X | \nX | \nX (≥ 5.7.5) | \nX (LWGEOM) | \n\n |
| \n | X | \nX | \nX | \nX | \nB | \n
| \n | X | \nX | \n\n | X | \nC | \n
| \n | X | \nX | \nX | \nX | \nB | \n
| \n | X | \nX | \nX | \nX | \nB | \n
| \n | X | \nX | \nX | \nX | \nB | \n
| \n | X | \n\n | \n | \n | C | \n
| \n | X | \n\n | \n | \n | C | \n
| \n | X | \n\n | \n | \n | B | \n
| \n | X | \n\n | \n | \n | B | \n
| \n | X | \n\n | \n | \n | C | \n
| \n | X | \n\n | \n | \n | C | \n
| \n | X | \n\n | \n | \n | C | \n
| \n | X | \n\n | \n | \n | C | \n
Le tableau suivant indique en résumé quelles sont les fonctions géométriques de base de données disponibles pour chaque moteur spatial.
\nFonction | \nPostGIS | \nOracle | \nMySQL | \nSpatiaLite | \n
|---|---|---|---|---|
| \n | X | \nX | \nX | \nX | \n
| \n | X | \n\n | X (≥ 5.7.5) | \nX | \n
| \n | X | \nX | \n\n | X | \n
| \n | X | \n\n | \n | X | \n
| \n | X | \n\n | \n | X | \n
| \n | X | \n\n | \n | X (LWGEOM) | \n
| \n | X | \nX | \n\n | \n |
| \n | X | \nX | \nX | \nX | \n
| \n | X | \nX | \nX | \nX | \n
| \n | X | \nX | \nX | \nX | \n
| \n | X | \n\n | X | \nX | \n
| \n | X | \n\n | \n | X | \n
| \n | X | \n\n | \n | \n |
| \n | X | \n\n | X (≥ 5.7.5) | \nX (LWGEOM) | \n
| \n | X | \nX | \nX | \nX | \n
| \n | X | \nX | \nX (≥ 5.7.5) | \nX (LWGEOM) | \n
| \n | X | \nX | \nX | \nX | \n
| \n | X | \n\n | \n | X | \n
| \n | X | \n\n | \n | X (LWGEOM) | \n
| \n | X | \n\n | \n | \n |
| \n | X | \nX | \nX | \nX | \n
| \n | X | \nX | \nX | \nX | \n
| \n | X | \nX | \n\n | X | \n
| \n | X | \nX | \n\n | X | \n
| \n | X | \nX | \n\n | X | \n
| \n | X | \n\n | \n | X | \n
| \n | X | \n\n | \n | X | \n
| \n | X | \nX | \nX | \nX | \n
| \n | X | \nX | \n\n | X | \n
| \n | X | \n\n | \n | X | \n
| \n | X | \nX | \nX | \nX | \n
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.
\nMéthode d’agrégation | \nPostGIS | \nOracle | \nSpatiaLite | \n
|---|---|---|---|
| \n | X | \n\n | X | \n
| \n | X | \nX | \nX | \n
| \n | X | \n\n | \n |
| \n | X | \n\n | X | \n
| \n | X | \nX | \nX | \n
Notes de bas de page
\n\n