Le portail francophone de la géomatique

En enlevant la première "(" et la dernière ")", est-ce que ça marche ?

Ça dépend de la source de données. Ça peut être "geometry" dans certaines sources. Il faudrait spécifier a. ou t. si il s'agit de couches géométriques, sinon c'est ambigu.

Dernière modification par carteq (Thu 02 April 2020 21:54)

En fait, ici c'est soit "a", ou "t" qui contient la géométrie, ou les deux. Il faut indiquer de quelle table on la prend pour les opérations. Tous les champs sont référencés a. ou .t dans la requête, mais jamais "geom". Ça ne fait pas de sens, surtout si "geom" se trouve dans chaque table. Il faudrait nommer les champs qui contiendront les 2 opérations géométriques aussi. Il faudrait aussi inclure le champ "geom" dans le Select puisqu'il fait partie du Group BY.

Tous les champs hors du Group BY doivent être agrégés également, si il y en a.

Dernière modification par carteq (Thu 02 April 2020 23:20)

@Julien, selon un post précédent, PLM veut faire ceci :

Du coup ma solution ne marchera pas.
Il faut donc que PLM améliore sa formule.

A mon avis, il faut que tu commences par cela :

Il faudrait aussi inclure le champ "geom" dans le Select puisqu'il fait partie du Group BY

Un exemple, en couche virtuelle :

Select
st_collect(st_line_interpolate_point(st_linemerge(geometry) ,(a.dist/st_length(geometry)))) as geom2 from LIGNES as a

"geometry" est implicite : pas visible dans la table d'attribut.

Si tu es en SQLite, c'est "geom" je crois. À vérifier.

"Pour les géopackages il faut faire des conversions explicites de la géométrie (qui est stockée sous un BLOB différent) pour pouvoir utiliser les fonctions spatiales.
Ou alors lancer la requête suivant avant les requêtes spatiales :


Code:
select EnableGpkgAmphibiousMode()"

Dernière modification par carteq (Fri 03 April 2020 00:30)

Autre exemple, plus proche :
Select a.CL_CHE_PLA, b.testing,
st_line_interpolate_point(st_linemerge(a.geometry), 1-(a.ID_CHEMIN/st_length(a.geometry)))  as geometry from LaCouche as a JOIN LaTable as b on a.ident = b.ident

Bonjour @carteq, j'ai changé de stratégie,

SELECT ST_Line_Interpolate_Point( geometry , 0.5),* FROM Troncon 

je vais utiliser cette formule mais je veux remplacer la valeur 0.5 par une valeur contenue dans un champs A de la couche tronçon comme l'écrire ?

Ensuite est-ce que cette formule peut fonctionner sachant que chacune de mes entités dans Troncon ont une une valeur de champ A différente.

Commençons par le SQL, essaie ceci :

SELECT ST_Line_Interpolate_Point( a.geometry , b.A)  FROM Troncon as a, Troncon as b where a.ID = b.ID

(non testé)

Normal, c'est juste un select. Si vous faites une couche de points vide, disons "pts" :

Insert Into pts (geometry) select  ST_Line_Interpolate_Point( a.geometry , b.danom) as geometry  FROM Troncon as a, itv2 as b where a.id_troncon = b.idtronco

- non testé

Je vais devoir tester.

Bonjour,
J'ai testé ce week-end différentes choses suivant vos conseils éclairés. Mais je suis également reparti de ma formule qui fonctionnait l'année dernière. Et je n'arrive une fois de plus à aucun résultat.
Je redétaille mon problème :


   
J'ai un shape de polyligne (en pièce jointe) avec des identifiants unique contenu dans le champ (id)
    J'ai un CSV (en pièce jointe) dans lequel les données en colonne sont entre autres :
   

       
Colonne (idtroncon)=(id de la shape) : Identifiant unique qui renvoi à l'identifiant de la polyligne,
        Colonne (danom) : Un pourcentage, correspondant à la distance du point à créer depuis le début de la ligne.
        Autres colonnes  : Informations à intégrer dans le futur shape de points
   
   


Je dois donc créer des points sur les polylignes à partir des informations inscrites dans le CSV.
Sachant qu'il y a plusieurs points à créer sur les même lignes à des distances différentes.

J'espère que les pièces jointes seront d'une plus grande aide.

Merci de votre aide. Je désespère.

J'ai fait un test avec un nombre limité de champs, en Spatialite.

Insert Into pts_finaux2 (geometry, distanomcc, codanom, idtroncon, danom )
select  ST_Line_Interpolate_Point( a.geometry ,  (b.danom * st_length(a.geometry))) as geometry, b.distanomcc, b.codanom, b.idtroncon, b.danom
FROM troncon as a,
itv2 as b
where a.id = b.idtroncon

Ça fonctionne.

Dernière modification par carteq (Mon 06 April 2020 21:47)