1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54
| library (spatstat) # package nécessaire au fonctionnement de la fonction runifpoint
A=runifpoint(7,win=owin(c(0,470),c(0,670))) # génère n points aléatoirement dans un rectangle de côté 470*670
plot(A)# permet d'avoir un apercu de la conformation
X=A$x #permet de creer un vecteur X avec les valeurs x de chaque point généré
Y=A$y #permet de creer un vecteur Y avec les valeurs y de chaque point généré
# on déterminer les coordonnées pour chaque point A,B,C,D,E,F,G, en considérant A comme le point de départ
Xa=X[1]
Ya=Y[1]
Xb=X[2]
Yb=Y[2]
Xc=X[3]
Yc=Y[3]
Xd=X[4]
Yd=Y[4]
Xe=X[5]
Ye=Y[5]
Xf=X[6]
Yf=Y[6]
Xg=X[7]
Yg=Y[7]
#déterminer ensuite la longueur de tous les segments possible
AB=sqrt((Xb-Xa)^2+(Yb-Ya)^2)
AC=sqrt((Xc-Xa)^2+(Yc-Ya)^2)
AD=sqrt((Xd-Xa)^2+(Yd-Ya)^2)
AE=sqrt((Xe-Xa)^2+(Ye-Ya)^2)
AF=sqrt((Xf-Xa)^2+(Yf-Ya)^2)
AG=sqrt((Xg-Xa)^2+(Yg-Ya)^2)
BC=sqrt((Xc-Xb)^2+(Yc-Yb)^2)
BD=sqrt((Xd-Xb)^2+(Yd-Yb)^2)
BE=sqrt((Xe-Xd)^2+(Ye-Yb)^2)
BF=sqrt((Xf-Xb)^2+(Yf-Yb)^2)
BG=sqrt((Xg-Xb)^2+(Yg-Yb)^2)
CD=sqrt((Xd-Xc)^2+(Yd-Yc)^2)
CE=sqrt((Xe-Xc)^2+(Ye-Yc)^2)
CF=sqrt((Xe-Xc)^2+(Ye-Yc)^2)
CG=sqrt((Xe-Xc)^2+(Ye-Yc)^2)
DE=sqrt((Xe-Xd)^2+(Ye-Yd)^2)
DF=sqrt((Xf-Xd)^2+(Yf-Yd)^2)
DG=sqrt((Xg-Xd)^2+(Yg-Yd)^2)
GE=sqrt((Xe-Xg)^2+(Ye-Yg)^2)
GF=sqrt((Xf-Xg)^2+(Yf-Yg)^2)
FE=sqrt((Xe-Xf)^2+(Ye-Yf)^2) |