I like to write a r-function to create 100 shape files in a loop. I
have three data files. Each file has 100 data fields (like sim1, sim2,
sim3....sim100) with a unique ID fields (Thana_ID). First I want to
create 100 data frame from these files and then merge or join each of
this 100 files with a polygon files individually to create 100 ESRI
shape files. I am trying to this following way. But I need do this in a
loop. Any suggestion will be appreciated.
Thanks
Zia
#--------------------------------------------------
library(sp)
library(spdep)
library(gstat)
library(maptools)
library(maps)
library(rgdal)
######
setwd("H:/GWPR")
gas<-read.csv("thana_sis_logGAS_intake_246.csv",header=TRUE)
was<-read.csv("thana_sis_logWAS_intake_246.csv",header=TRUE)
total<-read.csv("thana_sis_logTotal_intake_246.csv",header=TRUE)
names(thana)
names(was)
names(gas)
names(total
# Create 100 data-frames by selecting one column from each file.
sim1<-cbind(THANA_ID=thana$THANA_ID, logWAS=was$sim1, logGAs=gas$sim1,
logTotal=total$sim1)
sim2<-cbind(THANA_ID=thana$THANA_ID, logWAS=was$sim2, logGAs=gas$sim2,
logTotal=total$sim2)
.
.
.
sim100<-cbind(THANA_ID=thana$THANA_ID, logWAS=was$sim100,
logGAs=gas$sim100, logTotal=total$sim100)
#Join each data framewith a shape file and then write #ESRIShape files.
# Total number of shapefiles willbe100.
sim1.thana<-readShapePoly("bd_thana_246.shp")
sim1.thana@data <- merge(thana@data,sim1,by.x="THANA_ID",
by.y="THANA_ID", all.x=T, sort=F)
writeOGR(sim1.thana,dsn="/GWPR/Shape_Files", layer="sim1.thana",
driver="ESRI Shapefile")
sim2.thana<-readShapePoly("bd_thana_246.shp")
sim2.thana@data <- merge(thana@data,sim2,by.x="THANA_ID",
by.y="THANA_ID", all.x=T, sort=F)
writeOGR(sim2.thana,dsn="/GWPR/Shape_Files", layer="sim2.thana",
driver="ESRI Shapefile")
.
.
.
sim100.thana<-readShapePoly("bd_thana_246.shp")
sim100.thana@data <- merge(thana@data,sim100,by.x="THANA_ID",
by.y="THANA_ID", all.x=T, sort=F)
writeOGR(sim100.thana,dsn="/GWPR/Shape_Files", layer="sim100.thana",
driver="ESRI Shapefile")
--
---------------------
Zia Ahmed, PhD
Research Associate
Department of Crop and Soil Sciences
1002 Bradfield Hall, Cornell University
Ithaca, NY 14853-4203
t. 607.255.9387
f. 607.255.3207
email [email protected]
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
