Le portail francophone de la géomatique

bonjour

avant d'aller plus loin je m'aperçois que j'ai un réel souci de compréhension

si je considère les requêtes suivantes :

select x(coord),y(coord) from table where date between '2014-02-08 14:00:00' and '2014-02-08 14:00:01' and coord::geography && st_setsrid('BOX(-90 -90,90 90)::geography'::box2d,4326);

select x(coord),y(coord) from table where date between '2014-02-08 14:00:00' and '2014-02-08 14:00:01' and coord::geography && st_setsrid('BOX(0 -90,0 90)::geography'::box2d,4326);

select x(coord),y(coord) from table where date between '2014-02-08 14:00:00' and '2014-02-08 14:00:01' and coord::geography && st_setsrid('BOX(-90 -90,0 90)::geography'::box2d,4326);

et même en réduisant xmin et xmax à 0
select x(coord),y(coord) from table where date between '2014-02-08 14:00:00' and '2014-02-08 14:00:01' and coord::geography && st_setsrid('BOX(0 -90,0 90)::geography'::box2d,4326);

j'ai toujours le même résultat à savoir 50 lignes avec des valeurs en x  qui ne semblent pas contraintes par les bornes xmin etxmax comme par exemple :
      x      |     y     
-------------+------------
  -49.616782 | -28.158047
...
112.294614 |  -3.002758
   37.745654 |  -14.76219
-128.763314 | -25.480079
   15.080735 |   1.323699
...
  -48.793458 | -30.385019
(50 lignes)

si par contre je touche à la latitude le nombre de valeurs change

est-ce dû au cast des coordonnées (de type geometry) dans le type geography

merci bien

jean-Marie

Bonjour,

En fait, le problème vient de la bbox permettant de choisir les bons points:
a 90° (ou -90°) de latitude, on est au pole, sur un point ou la notion de longitude n'existe plus.
Les coordonnées (-90 -90) et (90 90) (ou toutes leurs combinaisons ayant 90 et -90 de latitude) ne permettent pas de définir un objet avec 4 sommets (polygone), mais juste le méridien (ligne). Ce n'est pas un objet valide pour l'index spatial de Postgis utilisé avec le type geography.

Concernant le problème de recherche, il se ramène à un cas avec le type geography:
vous voulez trouvez tous les objets situés sur une certaine surface. Pour définir cette surface, la distance entre les points, en longitude, doit etre inférieure à 180°, par convention dans Postgis.
Un extent (-20 20) est valide: les deux points seront joints par le plus petit arc sur le globe.
Un extent (-160 160) sera joint par le plus petit arc (40° de longitude) et non par un arc de 320°.

Du coup, pour chercher sur un grand extent (non representable avec une geography), il faut chercher les objets ne se trouvant pas dans le petit extent complémentaire (je devrais faire un schéma... )

Dans votre cas, il faudrait analyser l'extent de recherche et adapter la requete de recherche en fonction de cet extent.
C'est possible avec un fonction pl/pgsql.

Nicolas

Le plus simple ne serait-il pas de couper l'extent en deux : [-160;0] et [0;160], avec un OR entre les deux conditions ?

Bonjour,

Oui effectivement, c'est le plus simple.

Du coup ca doit pouvoir s'écrire en sql directement:
Un case sur l'extent de recherche: si largeur > 180 °: découper en deux extents,
sinon, garder l'extent initial et fabriquer un deuxieme extent "vide"

Dans la sous requete, utiliser tjs la meme condition, sur les deux extents générés plus haut:
geog && extent1 OR geog && extent2

Nico