11 Jul
2022
11 Jul
'22
3:19 p.m.
Message posté par : ODJAC33
----------------------------------------
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.
----------------------------------------
Ce message est accompagné de fichiers, pour les télécharger, suivre le lien ci-dessous.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=354179#p354179
Pour y répondre : qgis_fr(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
11 Jul
11 Jul
5:33 p.m.
Message posté par : Sylvain M.
----------------------------------------
-----------------
Citation :
Comment faire pour récupérer le point qui est à 250m. au Sud de l'extrémité Sud de la
piste SVP?
-----------------
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 :
-----------------
Code :
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
-----------------
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=354182#p354182
Pour y répondre : qgis_fr(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
12 Jul
12 Jul
11:26 a.m.
Message posté par : SANTANNA
----------------------------------------
Bonjour,
Autrement, en version full boîte à outils Traitement:
-----------------
Citation :
1 point situé au milieu d'un segment
-----------------
: L'algorithme "Centroide" (si les entités sont constituées d'un seul
segment) ou "Interpoler le point sur la ligne
(https://docs.qgis.org/3.22/fr/docs/user_manual/processing_algs/qgis/vectorg…)
avec la formule length($geometry)/2 pour le paramètre distance
-----------------
Citation :
1 point situé à une distance et un angle paramétré à partir de l'extrémité d'un
segment
-----------------
: une combinaison de Prolonger les lignes
(https://docs.qgis.org/3.22/fr/docs/user_manual/processing_algs/qgis/vectorg…)
et Extraire les vertex spécifiques
(https://docs.qgis.org/3.22/fr/docs/user_manual/processing_algs/qgis/vectorg…)
(valeurs 0,-1)
Ou les deux en 1, prolonger la ligne d'abord puis utiliser "Points le long de la
géométrie
(https://docs.qgis.org/3.22/fr/docs/user_manual/processing_algs/qgis/vectorg…)
" 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....
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=354192#p354192
Pour y répondre : qgis_fr(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
10:26 p.m.
Message posté par : Mathieu Denat
----------------------------------------
Re,
-----------------
ODJAC33 écrit :
Mais ... comment faire pour sortir un point qui est à une certaine distance de
l'extrémité d'un segment, mais pas dans l'axe de ce segment? (donc à un angle
connu que l'on devra paramétrer dans l'algorithme).
SQL obligatoire ou géotraitements existant?
-----------------
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)
-----------------
Mathieu Denat écrit :
[...]
with tmp as
(
SELECT 45 AS angle, --fixer l'angle, ici 45°
10000 AS rayon, --fixer la distance, ici 1km
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 ligne
)
[...]
-----------------
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=354212#p354212
Pour y répondre : qgis_fr(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
13 Jul
13 Jul
6:17 p.m.
Message posté par : ODJAC33
----------------------------------------
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.
----------------------------------------
Ce message est accompagné de fichiers, pour les télécharger, suivre le lien ci-dessous.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=354245#p354245
Pour y répondre : qgis_fr(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
1186
days inactive
1188
days old
4 comments
1 participants
participants (1)
-
qgis_fr＠georezo.net