Le portail francophone de la géomatique

MathieuB · Mon 28 June 2021 12:02

Bonjour à toutes et tous,

peux présent sur le forum depuis trop longtemp, me revoilà avec une incohérence que je n'arrive pas à comprendre, dans spatialite, pour les besoins d'un étudiant :

Je dispose de deux tables, une de point d'observation d'espèce (109418 lignes), une de mailles carrée (2921 enregistrements)
Je souhaite compter le nombre de taxons observés par maille mais le résultat est incohérent.

Une requête toute basique qui liste les mailles en intersection avec les données me retourne l'ensemble des mailles (2921 lignes) au lieu des 642 attendues...
Est-ce que je passe à coté d'une subtilité de spatialite ?

Code:

select  distinct carroyage.id
FROM carroyage, donnees
WHERE st_intersects(carroyage.geom, donnees.geom)

La syntaxe JOIN ON me retourne bien sûr le même résultat.

Merci d'avance pour vos éclairages éventuels...

tumasgiu · Mon 28 June 2021 13:36

Salut

y'a t'il dans données des lignes avec des géométries qui valent NULL ?

Si oui, st_intersects va alors renvoyer -1 qui apparemment serait évalué comme la valeur vrai
dans la clause WHERE, et donc tous les carreaux vont avoir au moins un match.

http://www.gaia-gis.it/gaia-sins/spatia … atest.html

Dernière modification par tumasgiu (Mon 28 June 2021 13:39)

MathieuB · Mon 28 June 2021 14:24

Merci Tumas !!!

C'était bien ça. Je n'y aurai jamais pensé ! En tout cas pas avant quelques jours de repos qui se font attendre :-)

J'ai donc complété la clause WHERE avec des conditions NOT NULL sur les deux géométries.
Et j'ai un résultat conforme à ce que j'obtiens avec PostGIS.

Code:

select  distinct carroyage.id
FROM carroyage JOIN donnees ON st_intersects(carroyage.geom, donnees.geom)
WHERE carroyage.geom IS NOT NULL AND donnees.geom IS NOT NULL

tumasgiu · Mon 28 June 2021 15:53

Le repos c'est important en effet

Je pense que la requête aurait même pu être simplifiée comme ceci :

Code:

select  distinct carroyage.id
FROM carroyage, donnees
WHERE st_intersects(carroyage.geom, donnees.geom) > 0

MathieuB · Tue 29 June 2021 09:53

Oui c'est sûr !

Je n'ai pas pris le temps de regarder la doc de sqlite.
Qui mentionne ceci :

2.1. Boolean Datatype
SQLite does not have a separate Boolean storage class. Instead, Boolean values are stored as integers 0 (false) and 1 (true).

il me faudra creuser mais cje trouve curieux que

Code:

SELECT -1 is true
ou 
SELECT -500 is true

retournent 1.
Comme si en fait toute valeur entière, positive ou négative autre que 0 soit évaluée comme "true"

Mais c'est surtout la doc de la fonction st_intersects de spatialite que j'aurais du lire !

The return type is Integer, with a return value of 1 for TRUE, 0 for FALSE, and –1 for UNKNOWN corresponding to a function invocation on NULL arguments;
convenience predicate: TRUE if the intersection of g1 and g2 is not empty

Merci encore Tumas.

tumasgiu · Tue 29 June 2021 11:43

il me faudra creuser mais cje trouve curieux que
Code:

SELECT -1 is true
ou
SELECT -500 is true

retournent 1.
Comme si en fait toute valeur entière, positive ou négative autre que 0 soit évaluée comme "true"

Je pense que ca vient du C (sqlite est écrit en C) ou le type booléen n'existe pas et que faux = 0 et vrai = tout le reste.

Ce que je me demande, c'est pourquoi les dev de spatialite ont décidé qu'une fonction dont au moins l'un des arguments
est NULL ne devait pas renvoyer NULL (je ne sais pas si c'est dans le standard SQL, mais c'est tout de même
une convention que j'ai tout le temps rencontré dans les autres sgbd).
Cela a peut être à voir avec le fait que spatialite utilise geos comme backend, et que les fonctions geos renvoie des entiers,
c'était peut être plus simple à développer comme çà.

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Mon 28 June 2021 12:02

Incohérence de résultat d'intersection Spatialite

Code:

#2 Mon 28 June 2021 13:36

Re: Incohérence de résultat d'intersection Spatialite

#3 Mon 28 June 2021 14:24

Re: Incohérence de résultat d'intersection Spatialite

Code:

#4 Mon 28 June 2021 15:53

Re: Incohérence de résultat d'intersection Spatialite

Code:

#5 Tue 29 June 2021 09:53

Re: Incohérence de résultat d'intersection Spatialite

Code:

#6 Tue 29 June 2021 11:43

Re: Incohérence de résultat d'intersection Spatialite

Pied de page des forums