Le portail francophone de la géomatique

Merci beaucoup Mathieu, le milieu d'un segment fonctionne très bien.

J'ai un peu plus de mal en revanche pour la 2° question.
Je voudrais récupérer les points qui prolongent de 250 m. une piste d'aviation (fichier piste UTM.zip ci-joint) de part et d'autre de chacune des 2 extrémités de cette piste (Nord et Sud).

J'ai tapé le code suivant :
with tmp as
(
    SELECT 0 AS angle, --fixer l'angle, ici 0°
           250  AS rayon, --fixer la distance, ici 250 m.
           id,
           force_rhr(geometry) geom --on fixe la fin de ligne toujours du même côté, pas forcément utile selon la situation
    FROM   pisteutm
)


SELECT
    setsrid(
    make_point(
            X(st_line_interpolate_point(geom,1)) + rayon*sin(angle), --0 pour le départ de la ligne, 1 pour l'arrivée
            Y(st_line_interpolate_point(geom,1)) + rayon*cos(angle) --0 pour le départ de la ligne, 1 pour l'arrivée
            )
    ,32636) --fixer la projection
    AS geom,
       id
FROM tmp

En sortie j'ai saisi une couche appelée Extremite, mais elle n'apparaît pas dans mes couches lorsque j'exécute la requête, ce qui prouve que l'exécution de la requête ne se passe pas bien.

2 remarques :
- la projection est bien UTM Zone 36N
- l'angle de 0 degré signifie pour moi que l'on prolonge l'axe piste de 250 m. vers le Nord, ne faut-il pas mettre 180 degrés pour la prolonger vers le Sud? (ce qui entrainerait qu'il faudrait faire 2 requêtes séparées (une pour la prolongation Nord et une pour la prolongation Sud), à moins que l'on puisse faire les 2 prolongations en une seule fois (excusez-moi, je ne sais pas programmer en SQL).

Merci encore de vos conseils.

Il me semble que la réponse est dans la proposition de Mathieu : "--0 pour le départ de la ligne, 1 pour l'arrivée", ce qui donne :

SELECT
    setsrid(
    make_point(
            X(st_line_interpolate_point(geom,0)) + rayon*sin(angle), --0 pour le départ de la ligne, 1 pour l'arrivée
            Y(st_line_interpolate_point(geom,0)) + rayon*cos(angle) --0 pour le départ de la ligne, 1 pour l'arrivée
            )
    ,32636) --fixer la projection
    AS geom,
       id
FROM tmp

Bonjour,
Autrement, en version full boîte à outils Traitement:

: L'algorithme "Centroide" (si les entités sont constituées d'un seul segment) ou "Interpoler le point sur la ligneavec la formule length($geometry)/2 pour le paramètre distance

: une combinaison de Prolonger les lignes et Extraire les vertex spécifiques (valeurs 0,-1)

Ou les deux en 1, prolonger la ligne d'abord puis utiliser "Points le long de la géométrie" avec comme distance length($geometry)/2 (j'ai cru comprendre qu'on prolongeait de la même distance des deux côtés).

Bon j'avoue, c'est moins classe que le SQL mais c'est pour dire que la boîte à outils regorge d'outils pour régler un bon nombre de tâches apparemment complexes....

Re,

Pour un angle fixe : la réponse était dans mon premier post!
En relisant bien les commentaires, vous trouverez la ligne à modifier pour changer l'angle.
Dans mon exemple, l'angle était de 45°.
Voir ce qui est en gras dans le 1er bloc (bloc tmp)

Si l'angle est variable (entré en paramètre dans le script), c'est plus compliqué!
Il faut faire une fonction (mais il existe peut-être des choses + clic-bouton).
La question est comment est calculé l'angle.
Si l'angle est dans un champ (donc calculé au préalable), il suffit de remplacer le 45 par le nom du champ qui contient la valeur de l'angle.

J'espère que c'est assez clair.
Bonne soirée,

Dernière modification par Mathieu Denat (Tue 12 July 2022 22:29)

Rebonjour,

Je joins un petit schéma très simple pour essayer de comprendre pourquoi le paramétrage d'une valeur angulaire ne correspond pas à l'attendu (cf. mon dernier post).

Le segment AB est mon segment initial, je veux sortir le point C qui est à une distance D de B (appelé rayon dans le code SQL) et à un angle alpha (appelé angle dans le code SQL) de l'axe AB.

Est t'on bien d'accord sur les 2 points suivants ?
- l'angle alpha n'est pas un angle par rapport au Nord (axe des Y sur mon schéma) mais par rapport à l'axe du segment AB (Y' sur mon schéma)
- X(st_line_interpolate_point(geom,0)) + rayon*sin(angle), Y(st_line_interpolate_point(geom,0)) + rayon*cos(angle) sont bien mesurés dans le repère (B, X', Y') et non le repère (O (intersection des axes des X et des Y), X, Y).

Merci de votre avis, bonne soirée.

Bonjour Mathieu,

Merci beaucoup pour tes précieux conseils, ça a super bien marché.

J'ai fait un test d'un point situé à 30° et 10000 m. de l'extrémité Nord d'une piste d'aviation (orientée plein Nord, la chance ...).
J'ai tapé le code suivant :
with tmp as
(SELECT 30 *pi()/180 AS angle, --fix angle and convert it from rad to deg
           10000  AS rayon, --fix distance ATTENTION CARTESIAN DISTANCE
           id,
           force_rhr(geometry) geom -- we always fix the end of the line on the same side, not necessarily useful depending on the situation
    FROM   pisteutm – the runway
)
SELECT
    setsrid(
    make_point(
            X(st_line_interpolate_point(geom,1)) + rayon*sin(angle), --0 for start of the line, 1 for end
            Y(st_line_interpolate_point(geom,1)) + rayon*cos(angle) --0 for start of the line, 1 for end
            )
    ,32636) –fix EPSG
    AS geom,
       id
FROM tmp


Et ça m'a donné le résultat en PJ.
Nickel, encore merci à toi.

Avec plaisir pour le coup de pouce!

Je ne suis pas sûr d'avoir bien saisi la question, je reformule:
En déplaçant les points, on peut tomber sur un point avec une altitude différente.
Tu veux récupérer l'altitude du nouveau point, c'est bien ça?

Si c'est le cas, il faut utiliser un MNT ou récupérer la donnée d'une couche extérieure.
Ça revient à faire une jointure (spatiale).
Je suis à côté de la plaque ou c'est bien ça la question?