1 Mar
2023
1 Mar
'23
5:11 p.m.
Message posté par : SylvaDEES
----------------------------------------
Bonjour,
je cherche à réaliser une requête spatiale SQL permettant :
Dans un même schémas, de sélectionner toutes les
géométries (Points) de la table "Point_topo"
se situant à plus de 50 centimètres - des géométries (Points ou Polylignes) de 10
autres tables ( les 10 conditions doivent être remplies ).
L'objectif est d'effectuer un nettoyage-suppression des géomètries isolées de la
table "Point_topo",
via une vue dédiée, ou export au format shape.
Merci pour vos contributions.
Cordialement.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=358836#p358836
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
2 Mar
2 Mar
8:51 a.m.
Message posté par : Sylvain M.
----------------------------------------
Voila comment je procéderais :
1) faire un UNION des géométries des 10 tables en remplaçant le champ de geométrie (ligne
ou point) par un buffer de 50cm autour de cette géometrie :
-----------------
Code :
CREATE TABLE tables_union AS
(
SELECT ST_Buffer(geom, 0.5) as geom
FROM Table1
UNION
SELECT ST_Buffer(geom, 0.5) as geom
FROM Table2
UNION
..
SELECT ST_Buffer(geom, 0.5) as geom
FROM Table10
)
-----------------
1) Selectionner dans la table "Point_topo" les géométries qui n'intersectent
pas celles de "tables_union". Par exemple via la fonction ST_Disjoint
(https://postgis.net/docs/ST_Disjoint.html) (inverse de ST_Intersects
(https://postgis.net/docs/ST_Intersects.html) )
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=358851#p358851
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
5:38 p.m.
Message posté par : Sylvain M.
----------------------------------------
En effet, c'est mieux ! :D
(l’élève a encore du chemin ;) )
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=358869#p358869
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
3 Mar
3 Mar
4:20 p.m.
Message posté par : SylvaDEES
----------------------------------------
Bonjour Sylvain, bonjour Nicolas,
merci pour vos précieuses contributions. N'étant pas un maître du sql, j'ai tenté
d'interpréter au mieux vos directives,
et j'ai rédigé une ébauche de requête spatiale ( forcément en erreur ):
1ère étape : création de la table d'union de la
géomètrie des 10 tables concernées :
CREATE TABLE union_tables AS
(
SELECT gid, insee, the_geom
FROM table 1
UNION
SELECT gid, insee, the_geom
FROM table 2
UNION
SELECT gid, insee, the_geom
FROM table 3
CREATE INDEX union_tables_idx
ON union_tables
USING gist
(the_geom) ;
)
erreur au niveau de la rédaction du create index
2ème étape : requête spatiale st_dwithin sur la
table "Pont_topo" depuis la table "union_tables" :
SELECT *
from Point_topo a, union_tables b
WHERE ST_dwithin (a.the_geom, b.the_geom, 0.5)
Merci d'avance pour votre aide et vos correctifs.
Cordialement.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=358888#p358888
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
7 Mar
7 Mar
1:14 p.m.
Message posté par : SylvaDEES
----------------------------------------
Bonjour Sylvain,
merci pour ta réponse aux étapes 1 et 2, c'est beaucoup plus claire. Je vais tester
les deux requêtes.
Pour la 2éme étape, je confirme également que je ne comprends pas la fonction du
"null" dans la requête proposée par Nicolas.
Je reviens vers vous dès les tests réalisés.
Cordialement.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=358941#p358941
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
9 Mar
9 Mar
6:12 p.m.
Message posté par : SylvaDEES
----------------------------------------
Bonjour ,
merci pour votre précieuse contribution d'expert "sql". Après correction
puis test des deux requêtes étape 1 et étape 2 sur une base postgres test, puis en
base de production, tout fonctionne bien. Je comprends mieux, désormais,
l'utilisation des opérateurs "where exists" et "where not
exists" qui facilitent les requêtes spatiales.
La requête étape 2 a bien sélectionné tout les Point-topo étant à plus 50 cm des
objets/géométrie de la table UNION.
Bonne journée.
Géomatiquement.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=359013#p359013
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
2 Mar
2 Mar
9:23 a.m.
Message posté par : Nicolas Ribot
----------------------------------------
Bonjour
Un peu pareil que Sylvain, mais:
• faire un index spatial apres l'union
• ne pas faire le buffer: ca prend du temps et pas besoin ici: utiliser st_dwithin(geom,
geom, dist)
• ne pas utiliser st_disjoint, qui ne peut pas utiliser d'index spatial, mais plutot
la construction
select ...
from table t
where not exists (
select null
from t2
where st_dwithin(t.geom, t2.geom, 0.5)
)
Nicolas
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=358852#p358852
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
7:19 p.m.
Message posté par : Nicolas Ribot
----------------------------------------
:d
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=358877#p358877
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
3 Mar
3 Mar
5:15 p.m.
Message posté par : Sylvain M.
----------------------------------------
Pour la première étape, l'index spatial doit être calculé après la création de la
table.
Là, vous l'avez mis entre les parenthèses du CREATE TABLE.
Voici ce qui devrait marcher :
-----------------
Code :
CREATE TABLE union_tables AS
(
SELECT gid, insee, the_geom
FROM table 1
UNION
SELECT gid, insee, the_geom
FROM table 2
UNION
SELECT gid, insee, the_geom
FROM table 3
);
-----------------
-----------------
Code :
CREATE INDEX union_tables_idx
ON union_tables
USING gist
(the_geom) ;
-----------------
Moi, sous Spatialite, je fais mes index spatiaux comme ça : "SELECT
CreateSpatialIndex('ma_table', 'geom');", mais je pense que votre
écriture est bonne aussi. J'aurais juste mis "the_geom"
(union_tables_the_geom_idx) dans le nom de l'index, pour indiquer qu'il est lié à
ce champ.
Pour la 2e étape, ce doit être quelque chose comme ça (mais j'avoue que je ne maitrise
pas trop cette proposition, et du coup, si ça ne marche pas, on aura encore besoin de
l'aide de Nicolas ;) En l'occurrence, ce que je n'ai pas compris, c'est le
"null" dans la requête de Nicolas) :
-----------------
Code :
SELECT *
FROM Point_topo
WHERE NOT EXISTS (
SELECT null
FROM union_tables
WHERE st_dwithin(Point_topo.geom, union_tables.the_geom, 0.5)
)
-----------------
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=358893#p358893
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
8 Mar
8 Mar
9:01 a.m.
Message posté par : Nicolas Ribot
----------------------------------------
Bonjour,
Dans cette forme SQL 'where exists' ou 'where not exists', on n'est
pas interessé par les colonnes du select, mais juste de savoir s'il y a des lignes
(rows) retournées par le sous-select.
Dans ce cas, on peut écrire select null, ou select 1, ou tout autre chose.
"Exists" va renvoyer vrai s'il y a des lignes dans la sous requete.
Dans l'exemple, on se sert de l'index spatial pour trouver les lignes qui
matchent, puis le NOT exists va prendre le complémentaire de ces lignes, c'est à dire
les lignes qui ne matchent pas la condition (ici, points a plus de 50cm).
(https://www.postgresql.org/docs/15/functions-subquery.html#FUNCTIONS-SUBQUE…)
Nicolas
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=358957#p358957
Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
947
days inactive
955
days old
9 comments
1 participants
participants (1)
-
geobd＠georezo.net