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]

Merci pour vos réponses.

Je risque de revenir avec d'autres questions

Et bonne journée à tous et à toutes.

[bouda53]

Pour calculer les angles formés dans mes parcours j'ai trouvé ça :"Calculer le produit vectoriel de ces deux vecteurs, ça donne le cosinus. Puis le déterminant, ça te donne le sinus, et pof, t'as ton angle".

Ok ça c'est compris, mon soucis pour le moment est :

a)-Comment calculer avec R le produit vectoriel de mes deux vecteurs Xi et Yi ?
b)-Comment calculer avec R le déterminant ?

Merci d'avance.

[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.

[alhazred]

"Calculer le produit vectoriel de ces deux vecteurs, ça donne le cosinus. Puis le déterminant, ça te donne le sinus, et pof, t'as ton angle".

Ce ne serait pas plutôt le produit scalaire ?