Creating a map using .shp and .csv

40 views
Skip to first unread message

Samaira Batool

unread,
May 13, 2014, 4:26:15 AM5/13/14
to ggp...@googlegroups.com
I have been struggling with a map that I have been trying to produce using library(maptools) and library(ggplot2) in R. I have a shapefile which contains names of 8 Towns in Faisalabad  and a csv file.  I have tried the following commands but the problem is that I can't merge the shspefile with my csv file:


library(ggplot2)
 setwd("C:/Users/Sam/Desktop/sam")
 data<-read.table("C:/Users/Sam/Desktop/sam/sidra.csv",header=TRUE,sep=",")
library(maptools)
FSDMap <- readShapeSpatial("Export_Output.shp")
plot(FSDMap)
str(FSDMap)

Tandianwala_Town <- data.frame(FSDMap@polygons[[1]]@Polygons[[1]]@coords)
colnames(Tandianwala_Town) <- c("Lat", "Long")
Tandianwala_Town$Town <- rep(x="Tandianwala", times=nrow(Tandianwala_Town)) 

Sumindri_Town <- data.frame(FSDMap@polygons[[2]]@Polygons[[1]]@coords)
colnames(Sumindri_Town) <- c("Lat", "Long")
Sumindri_Town$Town <- rep(x="Sumindri", times=nrow(Sumindri_Town))

Iqbal_Town  <- data.frame(FSDMap@polygons[[3]]@Polygons[[1]]@coords)
colnames(Iqbal_Town ) <- c("Lat", "Long")
Iqbal_Town $Town <- rep(x="Iqbal", times=nrow(Iqbal_Town )) tand 

Jinnah_Town <- data.frame(FSDMap@polygons[[4]]@Polygons[[1]]@coords)
colnames(Jinnah_Town) <- c("Lat", "Long")
Jinnah_Town$Town <- rep(x="Jinnah", times=nrow(Jinnah_Town))

Madina_Town <- data.frame(FSDMap@polygons[[5]]@Polygons[[1]]@coords)
colnames(Madina_Town) <- c("Lat", "Long")
Madina_Town$Town <- rep(x="Madina", times=nrow(Madina_Town)) 

Jarawanla_Town <- data.frame(FSDMap@polygons[[6]]@Polygons[[1]]@coords)
colnames(Jarawanla_Town) <- c("Lat", "Long")
Jarawanla_Town$Town <- rep(x="Jarawanla", times=nrow(Jarawanla_Town))

Layyalpur_Town <- data.frame(FSDMap@polygons[[7]]@Polygons[[1]]@coords)
colnames(Layyalpur_Town) <- c("Lat", "Long")
Layyalpur_Town$Town <- rep(x="Layyalpur", times=nrow(Layyalpur_Town))

Chak_Jhumra_Town <- data.frame(FSDMap@polygons[[8]]@Polygons[[1]]@coords)
colnames(Chak_Jhumra_Town) <- c("Lat", "Long")
Chak_Jhumra_Town$Town <- rep(x="Chak_Jhumra", times=nrow(Chak_Jhumra_Town))

FSDCoords <- 
         rbind(
                 Tandianwala_Town
               , Sumindri_Town
               , Iqbal_Town
               , Jinnah_Town
               ,Madina_Town
               , Jarawanla_Town
               , Layyalpur_Town
               , Chak_Jhumra_Town
               )

          
library(ggplot2)

FSDPlot <- 
         ggplot(data=FSDCoords, mapping = aes(x=Lat, y=Long, group=Town)) +
         geom_polygon(aes(fill = Town))
print(FSDPlot)
I would very much appreciate if somebody could help me out with this
Thank you very much in advance.
fshapefile.zip
sidra.csv
Reply all
Reply to author
Forward
0 new messages