Discussion :

[bouda53]

Bonjour à vous,

j'aimerais savoir comment comptabiliser le nombre de croisements dans un itinéraire, j'ai les coordonnées de position (x,y) ,en utilisant une boucle.

[A. D.]

Bonjour,

Personnellement je n'ai pas l'habitude de traiter de telles données et je dois dire que j'ai un peu de mal à comprendre votre problématique...

• A quoi ressemble vos données?
• Qu'appelez-vous un "itinéraire"? Un ensemble de coordonnées (x,y), ie. votre jeu de données?
• Comment se définit un "croisement"?

Cordialement,

A.D.

[bouda53]

Les X et les Y sont des coordonnées géo (de position).
Exemple :

Xi=(-22,7777 -23,2974 -23,7238 -23,7238 -24,5113 -24,7471 -24,7471 -24,4075 -22,881 -23,9479 -23,483 -23,8066 -23,6016 -23,1155)

Yi=(-24,4335 -27,6925 -30,1107 -30,1107 -29,2514 -29,2514 -29,2514 -29,2514 -26,8242 -24,3981 -22,5708 -19,2992 -16,0103 -12,7461)

Je veux calculer le nombre de croisements (d'intersections) dans le parcours (graphique (Xi,Yi) ), et calculer le cumul des angles formés dans ce parcours.

Merci.
J’espère que c'est plus clair.

[mathieu_r]

bonjour,

cela est possible avec la librairie spatstat.

Utiliser la classe psp planar segment pattern pour créer une ligne.

Vous aurez deux lignes.

Ensuite utiliser crossing.psp

Si deux lignes dans un même psp, utiliser self.crossing.psp

[bouda53]

Bonjour
si j'ai bien compris dans le ce cas là, je doit découper mes vecteur Xi et Yi
pour obtenir deux point uniquement par groupe (pour former les différent segments)

mais la je suis face à un second souci
la création d'une boucle pour automatiser la segmentation
EX: Xi= 1 2 3 4 5 6
Yi= 1 2 3 4 5 6

afin d'obtenir (1,1) (2,2) (3,3) (4,4) (5,5) (6,6)

sachant que le 1 er segment (la ligne "graphiquement") est former par le point de départ (1,1) et le point d'arrivé (2,2)
le 2 eme segment départ (2,2) arrivée (3,3) ....etc

merci d'avance

[lilly74]

Bonsoir,

pour ce qui est de l'automatisation de la création des segments, voici mon idée:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
x <- 1:6
y <- 1:6
 
u <- cbind(x,y)
u
 
l <- dim(u)[1]
 
sgmts <- cbind(u, rbind(u[2:l,], u[1,]))
sgmts

On obtient:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
> sgmts
     x y x y
[1,] 1 1 2 2
[2,] 2 2 3 3
[3,] 3 3 4 4
[4,] 4 4 5 5
[5,] 5 5 6 6
[6,] 6 6 1 1

J'espère que cela vous donnera une piste.

Bonne continuation

[bouda53]

Je vous présente ci-joint un tableau recapitulatif des points X et Y d'un parcours pour lequel j'aimerais pouvoir calculer le nombre de croisements.

Merci d'avance de l'aide.

[mathieu_r]

Bonjour,

Normalement, la fonction de spatstat appelée crossing permet de détecter les intersections entre lignes.