création de shape

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:

---

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

---

# je partitionne le nuage de point

Code:

---

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

---

# je créé des barycentres des partitions

Code:

---

points(Y$centers, pch = 13)

---

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

Code:

---

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

bonjour,

j'ai trouvé :P

#transformation pt en shape (library maptools, sp, rgdal )

Code:

---

1 2 3

coordinates (pt) <- ~ V1 + V2 proj4string(pt)=CRS("++init=epsg:2154") pt.df <- SpatialPointsDataFrame(pt, data.frame(id=1:length(pt)))

---

#création shape

Code:

---

writeOGR(pt.df, dsn="." ,layer="t1",driver="ESRI Shapefile")

---

si une façon alternative existe, je suis preneur.
fcka