Le portail francophone de la géomatique

Sophie_Géo · Mon 23 May 2016 21:21

Bonjour,

J'ai un fichier de noeuds et un fichier de polylignes, je suis sous spatialite (c'est le même SQL que pgsql).

Dans mon fichier polylignes, j'ai 2 champs, "NOEUDA" et "NOEUDB" qui sont vides.
Dans mon fichier noeuds, j'ai un champ "identifiant" unique pour chaque noeud.

Je souhaite récupérer grâce à une petite ligne de code l'identifiant du noeud à chaque extrémité de ma polyligne.

Je sais qu'il y a un st_intersect, st_startpoint, st_endpoint mais je ne sais pas comment les articuler pour remplir mes 2 champs dans ma table polylignes...

J'ai essayé ça, qui sélectionne bien mes points intersectant avec mes lignes mais comment remplir ensuite mes champs NOEUD A et NOEUD B? :
"select id from polylignes, noeuds where st_intersects(polylignes.geom,noeuds.geom)"

Merci pour votre aide précieuse

Sophie

SIGEAL · Tue 24 May 2016 07:45

Bonjour,

Une sous-requête peut faire l'affaire (pas très rapide, mais ça fonctionne) :

Code:

select id,
(select id_arret from noeuds n where st_within(st_startpoint(p.geometry), st_buffer(n.geometry, 1))) n_deb,
(select id_arret from noeuds n where st_within(st_endpoint(p.geometry), st_buffer(n.geometry, 1))) n_fin
from polylignes p;

À adapter pour mettre à jour la table polylignes.

Sophie_Géo · Fri 27 May 2016 16:25

Bonjour,

Merci, ça me mets déjà pas mal sur la piste. Cependant ça ne me donne pas les bons NOEUDA et NOEUDB aux extrémités des polylignes...
Et je ne vois pas du tout comment remplir ces données dans mes champs NOEUDA et NOEUDB.

Merci pour votre aide,

Sophie

SIGEAL · Fri 27 May 2016 17:22

Bonjour,

ça ne me donne pas les bons NOEUDA et NOEUDB aux extrémités des polylignes

Pouvez-vous être plus précise ?
Vos couches sont-elles topologiquement correctes (noeuds superposés aux extrémités de polylignes) ?
Vos polylignes sont-elle simples, ou des collections (multilinestring) ?

Pour la mise à jour, il faut utiliser update :

Code:

update polylignes set noeuda = 
(select a.id_arret from noeuds as a where st_within(st_startpoint(polylignes.geometry), st_buffer(noeuds.geometry, 1)))

Requête non testée, une recherche sur spatialite update join vous donnera plus d'informations.

Sophie_Géo · Mon 30 May 2016 10:39

Bonjour,

Oui je connais plutôt bien le problème.

Lors de l'import des shapes dans spatialite il me convertit les lignes en multilignes et les points en multipoint. C'est la raison pour laquelle je suppose que ça ne fonctionne pas super bien. D'ailleurs sur un ST_stratpoint seul, mes champs sont tous vides...

J'ai essayé de faire comme sur postgresql :

alter table matable alter column the_geom type geometry(Linestring, <srid>)using st_geometryN(the_geom, 1);

Mais spatialite ne tolère pas les commandes alter table et alter column. Je sais pas par quoi les remplacer?

Merci pour votre réponse.

Sophie

SIGEAL · Mon 30 May 2016 10:53

Bonjour,

Plutôt que de modifier le type de colonne, ce qui n'est effectivement pas possible sous SQLite, essayez d'extraire la première polyligne de la collection :

Code:

update polylignes set noeuda = 
(select a.id_arret from noeuds as a where st_within(st_startpoint(st_geometryN(polylignes.geometry, 1)), st_buffer(noeuds.geometry, 1)))

Il faudrait cependant vérifier préalablement la topologie de votre couche : Y a-t-il des polylignes avec plusieurs secitons ? Si oui pourquoi ?

Sophie_Géo · Mon 30 May 2016 11:07

Il n'y a pas de polylignes en plusieurs sections, c'est juste une 4 lignes droites parfaitement jointes par des points réparties un peu partout dans l'espace. J'ai crée un projet d'exemple tout simple.

Dernière modification par Sophie_Géo (Mon 30 May 2016 11:08)

SIGEAL · Mon 30 May 2016 14:05

Alors ça devrait fonctionner, sinon vous pouvez mettre vos données en pièce jointe.

Sophie_Géo · Mon 30 May 2016 15:13

Je dois avouer que cela m'as un peu gonflée donc je suis passée sur pgsql et j'ai correctement transformé mes multilignes en lignestring.
Du coup j'ai lancé votre 1ère requête avec le update et il me dit en erreur :

ERREUR: référence invalide d'une entrée de la clause FROM pour la table « noeuds »

SIGEAL · Mon 30 May 2016 15:40

Il y avait une coquille dans mon exemple :

Code:

update polylignes set noeuda = 
(select a.id_arret from noeuds as a where st_within(st_startpoint(st_geometryN(polylignes.geometry, 1)), st_buffer(a.geometry, 1)))

Sophie_Géo · Mon 30 May 2016 16:58

Bonjour,

Merci pour votre aide.

Nous avons cherché avec mon collègue et cette requête fonctionne aussi:

Code:

update polylignes set noeuda = ( select idnoeud from (select idnoeud,polylignes.gid from noeuds, polylignes where st_within(st_startpoint(polylignes.geom), st_buffer(noeuds.geom, 1))) as a where a.gid=polylignes.gid)

Merci à vous, vous nous avez été d'une grande aide!

Sophie

Dernière modification par Sophie_Géo (Mon 30 May 2016 16:58)

Sophie_Géo · Mon 30 May 2016 17:25

EDIT :

J'ai tout de même essayé votre requête et ça ne fonctionne pas.

Il me renvoie toujours l'erreur :
ERREUR: plus d'une ligne renvoyée par une sous-requête utilisée comme une expression

Il manque dans votre ligne une condition " = " à une table de référence qui point vers une ligne en particulier et non une matrice...

Du coup notre ligne de code focntionne.

Si ça peut être utile.

Merci

Ludivine · Sat 28 January 2017 08:49

Bonjour,

Je souhaiterai faire exactement ce qui est décrit ici, mais lorsque j'utilise le code (en modifiant les noms de couches/champs), ça ne fonctionne pas, j'ai une erreur de syntaxe...

Je n'ai pas l'habitude de passer par les lignes de code, aussi peut-être que je fais... n'importe quoi!... J'ai essayé de rechercher comment faire, mais je m'y perds un peu...

Pour utiliser le code, j'ouvre la console Python, et je colle directement le code en modifiant ce qu'il faut: chez moi ça donne ça:

update pipe set NODE1 = ( select DC_ID from (select DC_ID,pipe.gid from junction, pipe where st_within(st_startpoint(pipe.geom), st_buffer(junction.geom, 1))) as a where a.gid=pipe.gid)

puis je valide...

Quelqu'un pourrait-il m'aider un peu?

Merci
Ludivine

(j'utilise QGIS 2.18)

EDIT:
Bon, faute de réponses, j'ai fait autrement, par des créations de champs intermédiaires + jointures...
Toutefois, j'aurai bien aimé comprendre malgré tout comment faire ce type de manip... Quelqu'un pourrait-il m'éclairer ou m'indiquer ou trouver l'info (malgré mes recherches, je n'ai rien trouvé)?
Merci

Dernière modification par Ludivine (Fri 03 February 2017 09:53)

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 23 May 2016 21:21

QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

#2 Tue 24 May 2016 07:45

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

Code:

#3 Fri 27 May 2016 16:25

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

#4 Fri 27 May 2016 17:22

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

Code:

#5 Mon 30 May 2016 10:39

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

#6 Mon 30 May 2016 10:53

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

Code:

#7 Mon 30 May 2016 11:07

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

#8 Mon 30 May 2016 14:05

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

#9 Mon 30 May 2016 15:13

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

#10 Mon 30 May 2016 15:40

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

Code:

#11 Mon 30 May 2016 16:58

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

Code:

#12 Mon 30 May 2016 17:25

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

#13 Sat 28 January 2017 08:49

Re: QGIS: Récupérer les noms des Noeuds aux extrémités d'une polyligne

Pied de page des forums