Le portail francophone de la géomatique

Bonjour,

Une piste avec la même erreur : https://www.sigterritoires.fr/index.php … irtuelles/
Je n'ai pas vérifié dans le détail si le bug indiqué dans le lien correspond à ta requête. A creuser donc.

Dernière modification par p.jeremie (Fri 02 February 2018 19:48)

Bonjour,

Avec postgis, les fonctions make_line, make_point n'existent pas:
st_makePoint, st_makeLine, etc..

Nicolas

Il faut peut etre que je précise mon objectif. En fait, il me faut faire une carte des migrations professionnelles (carte de flux). J'ai récupéré une BD insee sur ma zone d'intéret (28 communes).
J'ai importé cette table dans postgres. 1 enregistrement = 1 individu.
Ma table se compose de plusieurs champs :
id unique
codgeo_res: identifiant de la commune de residence
x_resi: longitude de la commune de residence
y_resi : latitude de la commune de residence
codgeo_tr: identifiant de la commune de travail
x_travail: longitude de la commune de travail
y_travail :latitude de la commune de travail
codgeo_tr : identifiant de la commune de travail
ipondi : poids de l'individu (il s'agit d'une valeur ponderation affectée à chaque individu)
geom : geometrie de la commune de residence


L'idée serait donc de pouvoir generer une table pour laquelle 1 enregistrement correspondrait à un flux (commune A vers  commune B).
Il me faut donc generer les geometries utiles (ponctuels + fleches lineaires. Je dois  faire un SUM sur le champ ipondi et faire un group by sur le champ codgeo_residence...

J'espère etre clair. Si vous pensez pouvoir m'aider afin que je puisse traduire cet objectif en language sql fonctionnel...ce serait super...un grand merci .

Dernière modification par image95 (Sat 03 February 2018 13:23)

Merci à tous pour vos messages. Alors, je viens de lancer la requete dans pgadmin, apres avoir fait un cast STRING vers BIGINTnon seulement sur mes XY residence mais egalement pour les XY travail. Pour la premiere fois cette requete fonctionne.
PAr contre, comme le disait tumasgiu, cette requete fonctionne à condition que j'indique dans mon group by tous les champs presents dans mon select...mais le fait de faire celà modifir totalement le resultat en sortie non? Car moi je voulais faire un group by sur le codgeo de la residence, pas sur tous les champs !
Voici en pj, le resultat de ma requete, j'ai chargé la table temporaire créée sous qgis...euh bon, on voit bien qu il a créé la couche LIGNE et ca  semble ok pour les points residence...mais par contre les points travail semble totalement dans les choux...
Oui, bonne idée, je peux vous partager une extraction de mon fichier.

Si vouos pensez pouvoir m'eclairer? on avance quand meme...

Merci beaucoup pour votre aide.

Dernière modification par image95 (Sat 03 February 2018 19:13)

Bonsoir,

Revoyez la doc sur le GROUP BY: vous comprendrez pourquoi vous ne pouvez pas regrouper sur tous les champs:
Les champs qui ne sont pas dans le clause GROUP doivent etre utilisées dans des fonction agrégées, c'est a dire faisant un regroupement pour les valeurs des colonnes.

Si vous voulez garder toutes les valeurs pour une colonne non présente dans un GROUP, vous pouvez utiliser la fonction agrégée array_agg(col), qui produira un tableau contenant les valeurs pour la colonne "col".

Nicolas

Bonsoir,

Si je comprends bien, vous voulez représenter les flux de transports de ville a ville a partir de trajets d'individus.

On peut commencer par créer ces trajets pour chaque individu:

select e.gid, e.id, codgeo_res, libgeo_res,
  codgeo_tr, libgeo_tr, ipondi,
  st_makeline(st_makepoint(e.x_resi, e.y_resi),
              st_makepoint(e.x_tr, e.y_tr)) AS geom
from entrants e;

Ici, pas besoin de group by, on veut juste une nouvelle geom pour la table, avec les points de départ/arrivée des trajets.

Ensuite, pour "cumuler" ces trajets par ville et avoir une idée du poids, il doit y avoir pleins de méthodes différentes, selon que vous souhaitiez une représentation graphique, ou des stats, ou autre chose.

On peut par exemple "cumuler" les trajets identiques et les compter, pour un couple de villes donné. Ca permet d'avoir un poids plus important si un trajet est fait plus souvent (j'imagine que c'est ce que permet de calculer ipondi).

En faisait l'union des géométries linestring des trajets par couple de villes, on "elimine" les trajets identiques en les fusionnant, tout en gardant le compte du nombre de trajet (et somme ipondi)

On a alors tous les trajets entre villes avec des critères pour représenter différemment les trajets :

create table entrants_res as
with tmp as (
  select e.gid, e.id, codgeo_res, libgeo_res,
  codgeo_tr, libgeo_tr, ipondi,
  st_setSRID(st_makeline(st_makepoint(e.x_resi, e.y_resi),
              st_makepoint(e.x_tr, e.y_tr)), 2154) AS geom
  from entrants e
) SELECT codgeo_res, libgeo_res, codgeo_tr, libgeo_tr, sum(ipondi) as sumipondi,
  sum(ipondi)/1000 as linewidth,
  count(*) as num_tr, st_union(geom)::geometry(linestring, 2154) as geom
  FROM tmp
group by codgeo_res, libgeo_res, codgeo_tr, libgeo_tr;

Un exemple de représentation dans QGIS en PJ, avec largeur de la ligne issue du champ calculé "linewidth":
plus la ligne est large, plus ce trajet a été parcouru.

Nicolas

Nicolas

Cher Nicolas, je ne sais pas comment vous remercier. En fait, de mon coté, j'avais également réfléchi au st_union. Voilà ce que j'avais trouvé :

create table temp3 as select st_union(st_makeline(st_makepoint(entrants_v2.x_resi,entrants_v2.y_resi),st_makepoint(entrants_v2.x_tr,entrants_v2.y_tr))) as geometry,sum(entrants_v2.ipondi) as sumipondi, 
entrants_v2.codgeo_res as codgeo_res, entrants_v2.codgeo_tr as codgeo_tr

from entrants_v2
group by entrants_v2.codgeo_res, entrants_v2.codgeo_tr

voici une petite SS. Notre resultat semble similaire.Votre zone est plus restreinte car je vous avais transmis qu'une petite selection. Vous avez trouvé un bon style, en jouant à la fois sur la taille et la couleur.

Mais là, dans l'etat actuel des choses, la carte permet d'apprehender l'intensité du flux, mais pas le sens. Il y aurait il un moyen d orienter le flux (commune A vers B. Commune P vers F)...

De plus, Il s'agissait des flux entrants uniquement. Maintenant, j'imagine qu'il faut procéder de la meme facon pour les flux sortants? Je possède la table SORTANTS egalement.
Il doit y avoir moyen de faire les traitements entrants/sortants en une seule étape...en tout cas, je trouve cet exercice passionnant.

Dernière modification par image95 (Sat 03 February 2018 23:21)

Oui, on peut visualiser le sens: qgis permet de dessiner des fleches.

Un bel exemple de flux avec qgis: https://anitagraser.com/2016/12/18/deta … flow-maps/

Oui, il est possible de faire en une etape. On peut quasiment tout faire en une étape en SQL, mais souvent, découper en étapes logiques c'est pas mal.

Nicolas

OK super, merci...oui, je vais donner suite aux différents posts...

Bonjour,

Un travail similaire a faire...donc un très grand merci aux protagonistes de ces échanges pour leur réflexion et informations...j'y vois plus clair !

Si par contre il vous serez possible d'apporter une précision sur un point.

Dans cette requête on créer des points et des lignes, d’où cette expression

"....st_union(geom)::geometry(linestring, 2154)..."

Lorsque je l'inscrit dans ma requête il me sort une erreur a cause des points creer, existe il une alternative au st_union pour isoler et ne prendre que les lignes ??? 

Ma langue de Shakespeare  étant un peu rouillé, vous connaissez des liens en français permettant de travailler la lisibilité de ces flux une fois réalisé sous Qgis.

Encore merci

Merci tumasgiu, effectivement, je suis passé par un st_collectionextract, afin de ne récupérer que les lignes ce qui est la memee chose .

j'essaye de continuer le proccess de Nicolas, mais j'ai un peu de mal :

-  A quoi correspond  la fonction idtrajet%2
- La fonction st_curvetoline(format('CIRCULARSTRING(%s %s, %s %s, %s %s)',

qu'est ce que le 'CIRCULARSTRING(%s %s, %s %s, %s %s)' ???

Lorsque j'effectue cette requête le résultat renvoi parse error at position 18 within geometry

Si vous aviez la possibilité de détailler un peu plus cette requête ça serait sympa.

Bonne soirée

Dernière modification par Theos2000 (Mon 12 February 2018 19:51)

Merci pour toutes ces précisions et surtout l'asttuce de l'Id !

Pour l'instant ça bloque, je vais me contenter de lignes droites et approfondir ce "WKT"

Merci d'avoir pris le temps ! :-)