Discussion :

[fcka0001]

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)

Code :

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

shape <- readOGR (dsn=".",layer="pt2")

# je partitionne le nuage de point

Code :

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

1
2
X<- cbind(shape$XCOORD, shape$YCOORD)
Y <-kmeans (X,8)

# je créé des barycentres des partitions

Code :

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

points(Y$centers, pch = 13)

# je récupère les x y des barycentres

Code :

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

Z<-cbind(Y$centers)

et à partir de là je coince...
Je voudrai transformer les x y en shape (epsg:2154). Si quelqu'un peut m'aider... merci d'avance.

fcka