Hi Zia, I don't know whether it will work for you, but I was able to read in fairly boundary polygons (11k lines) and merge them using spRbind()
Not sure if it's exactly what you want as I'm pretty much a beginner over here on r-sig-geo... cheers, Sean On Sat, Jan 23, 2010 at 3:45 AM, Zia Ahmed <z...@cornell.edu> wrote: > Hi Matt: > With help of Jim Holtman (my R Guru!), I am able to modified your code for > selecting particular type of shape files from different folders and able > to combine them together. Thank you so much! > However, final polygon does not contain any filed (suppose county name, > FIPS code etc) of the shape files those I used for merging . How do I keep > these fields in my final polygon. Any idea? > Thanks > Zia > > # Rode for combing Code > library(rgdal) > setwd("D:/test") > files <- list.files(pattern="soilmu_a_ny.*.shp$", > recursive=TRUE,full.names=TRUE) > > # Get polygons from first file > > data.first <- readOGR(files[1],gsub("^.*/(.*).shp$", "\\1", files[1])) > polygons <- slot(data.first, "polygons") > > # add polygons from remaining files > > for (i in 2:length(files)) { > data.temp <- readOGR(files[i], gsub("^.*/(.*).shp$", "\\1",files[i])) > polygons <- c(slot(data.temp, "polygons"),polygons) > } > > # rename IDs of Polygons > > for (i in 1:length(polygons)) { > slot(polygons[[i]], "ID") <- paste(i) > } > > # ain't that spatial > > spatialPolygons <- SpatialPolygons(polygons) > spdf <- SpatialPolygonsDataFrame(spatialPolygons, > data.frame(NewID =1:length(polygons))) > > # output combined results > plot(spatialPolygons) > writeOGR(spdf, dsn="C:/test/combined.shp", layer="combined", driver="ESRI > Shapefile") > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo