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

[fcka0001]

bonjour,

j'ai trouvé

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

Code :

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

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 :

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

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

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