bonjour
je cherche à créer un fichier shape à partir de coordonnées RGF 93.
ci-dessous mon code:
#je pars d'un fichier shape (il contient un nuage de point)
# je partitionne le nuage de point
Code : Sélectionner tout - Visualiser dans une fenêtre à part shape <- readOGR (dsn=".",layer="pt2")
# je créé des barycentres des partitions
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2 X<- cbind(shape$XCOORD, shape$YCOORD) Y <-kmeans (X,8)
# je récupère les x y des barycentres
Code : Sélectionner tout - Visualiser dans une fenêtre à part points(Y$centers, pch = 13)
et à partir de là je coince...
Code : Sélectionner tout - Visualiser dans une fenêtre à part Z<-cbind(Y$centers)
Je voudrai transformer les x y en shape (epsg:2154). Si quelqu'un peut m'aider... merci d'avance.
fcka
Partager