R Notebook

setwd("C:/Users/Sathiya Narayanan/Desktop/ANALYTICS/Maps")

The working directory was changed to C:/Users/Sathiya Narayanan/Desktop/ANALYTICS/Maps inside a notebook chunk. The working directory will be reset when the chunk is finished running. Use the knitr root.dir option in the setup chunk to change the the working directory for notebook chunks.

library(sp)

package <U+393C><U+3E31>sp<U+393C><U+3E32> was built under R version 3.4.3

ind1$NAME_1 = as.factor(ind1$NAME_1)
ind1$fake.data = runif(length(ind1$NAME_1))
spplot(ind1,"NAME_1",  col.regions=rgb(0,ind1$fake.data,0), colorkey=T, main="Indian States")

ind1$NAME_1

 [1] Andaman and Nicobar    Andhra Pradesh         Arunachal Pradesh     
 [4] Assam                  Bihar                  Chandigarh            
 [7] Chhattisgarh           Dadra and Nagar Haveli Daman and Diu         
[10] Goa                    Gujarat                Haryana               
[13] Himachal Pradesh       Jammu and Kashmir      Jharkhand             
[16] Karnataka              Kerala                 Lakshadweep           
[19] Madhya Pradesh         Maharashtra            Manipur               
[22] Meghalaya              Mizoram                Nagaland              
[25] NCT of Delhi           Odisha                 Puducherry            
[28] Punjab                 Rajasthan              Sikkim                
[31] Tamil Nadu             Tamil Nadu             Telangana             
[34] Tripura                Uttar Pradesh          Uttarakhand           
[37] West Bengal           
36 Levels: Andaman and Nicobar Andhra Pradesh Arunachal Pradesh ... West Bengal

wb1 = (ind1[ind1$NAME_1=="Tamil Nadu",])
spplot(wb1,"NAME_1", col.regions=rgb(0,0,1), main = "Tamil Nadu, India",scales=list(draw=T), colorkey =F)

kt1 = (ind1[ind1$NAME_1=="Karnataka",])
spplot(kt1,"NAME_1", col.regions=rgb(0,1,0), main = "Karnataka, India",scales=list(draw=T), colorkey =F)

ind2=readRDS("IND_adm2.rds")
ind3=readRDS("IND_adm3.rds")

wb2 = (ind2[ind2$NAME_1=="West Bengal",])
spplot(wb2,"NAME_1", main = "West Bengal Districts", colorkey =F)

wb2$NAME_2 = as.factor(wb2$NAME_2)
wb2$fake.data = runif(length(wb2$NAME_1)) 
spplot(wb2,"NAME_2",  col.regions=rgb(0,wb2$fake.data, 0), colorkey=T)

col = rainbow(length(levels(wb2$NAME_2)))
spplot(wb2,"NAME_2",  col.regions=col, colorkey=T)

col_no = as.factor(as.numeric(cut(wb2$fake.data, c(0,0.2,0.4,0.6,0.8,1))))
levels(col_no) = c("<20%", "20-40%", "40-60%","60-80%", ">80%")
wb2$col_no = col_no
myPalette = brewer.pal(5,"Greens")
spplot(wb2, "col_no", col=grey(.9), col.regions=myPalette, main="District Wise Data")

wb3 = (ind3[ind3$NAME_1=="West Bengal",])

wb3$NAME_3 = as.factor(wb3$NAME_3)
col = rainbow(length(levels(wb3$NAME_3)))
spplot(wb3,"NAME_3", main = "Taluk, District - West Bengal", colorkey=T,col.regions=col,scales=list(draw=T))

wb3$NAME_2

 [1] "Alipurduar"         "Bankura"            "Bankura"           
 [4] "Barddhaman"         "Barddhaman"         "Barddhaman"        
 [7] "Barddhaman"         "Barddhaman"         "Birbhum"           
[10] "Birbhum"            "Birbhum"            "Dakshin Dinajpur"  
[13] "Darjiling"          "Darjiling"          "Darjiling"         
[16] "Darjiling"          "Haora"              "Haora"             
[19] "Hugli"              "Hugli"              "Hugli"             
[22] "Hugli"              "Jalpaiguri"         "Koch Bihar"        
[25] "Koch Bihar"         "Koch Bihar"         "Koch Bihar"        
[28] "Koch Bihar"         "Kolkata"            "Maldah"            
[31] "Murshidabad"        "Murshidabad"        "Murshidabad"       
[34] "Murshidabad"        "Nadia"              "Nadia"             
[37] "Nadia"              "North 24 Parganas"  "North 24 Parganas" 
[40] "North 24 Parganas"  "North 24 Parganas"  "North 24 Parganas" 
[43] "Pashchim Medinipur" "Pashchim Medinipur" "Pashchim Medinipur"
[46] "Purba Medinipur"    "Purba Medinipur"    "Puruliya"          
[49] "South 24 Parganas"  "South 24 Parganas"  "South 24 Parganas" 
[52] "Uttar Dinajpur"     "Uttar Dinajpur"    

wb3 = (ind3[ind3$NAME_1=="West Bengal",])
n24pgns3 = (wb3[wb3$NAME_2=="North 24 Parganas",])
spplot(n24pgns3,"NAME_3", colorkey =F, scales=list(draw=T), main = "24 Pgns (N) West Bengal")

n24pgns3$NAME_3 = as.factor(n24pgns3$NAME_3)
n24pgns3$fake.data = runif(length(n24pgns3$NAME_3))
spplot(n24pgns3,"NAME_3", col.regions=rgb(0, n24pgns3$fake.data, 0), colorkey=T,scales=list(draw=T))

s24pgns3 = (wb3[wb3$NAME_2=="South 24 Parganas",])
spplot(s24pgns3,"NAME_3", colorkey =F, scales=list(draw=T), main = "24 Pgns (S) West Bengal")

s24pgns3$NAME_3 = as.factor(s24pgns3$NAME_3)
s24pgns3$fake.data = runif(length(s24pgns3$NAME_3))
spplot(s24pgns3,"NAME_3", col.regions=rgb(0, s24pgns3$fake.data, 0), colorkey=T,scales=list(draw=T),main = "24 Pgns (S) West Bengal")

mur3 = (wb3[wb3$NAME_2=="Murshidabad",])
spplot(mur3,"NAME_3", colorkey =F, scales=list(draw=T), main = "Murshidabad West Bengal")

mur3$NAME_3 = as.factor(mur3$NAME_3)
mur3$fake.data = runif(length(mur3$NAME_3))
spplot(mur3,"NAME_3", col.regions=rgb(0,0, mur3$fake.data), colorkey=T,scales=list(draw=T),main = "Murshidabad West Bengal")

n24pgns3 = (wb3[wb3$NAME_2=="North 24 Parganas",])
basirhat3 = (n24pgns3[n24pgns3$NAME_3=="Basirhat",])
spplot(basirhat3,"NAME_3", colorkey =F, scales=list(draw=T), main = "Basirhat,24 Pgns (N) West Bengal")

mur3 = (wb3[wb3$NAME_2=="Murshidabad",])
bahar3 = (mur3[mur3$NAME_3=="Baharampur",])
spplot(bahar3,"NAME_3", colorkey =F, scales=list(draw=T), main = "Baharampur, Murshidabad, West Bengal")

wb4 = (ind2[ind2$NAME_1=="West Bengal",])
wb4$NAME_2 = as.factor(wb4$NAME_2)
col = rainbow(length(levels(wb2$NAME_2)))
spplot(wb2,"NAME_2",  col.regions=col,scales=list(draw=T),ylim=c(23.5,25),xlim=c(87,89), colorkey=T)

library(ggmap)
TN2 = (ind2[ind2$NAME_1=="Tamil Nadu",])
nam = c("Chennai","Madurai","Hosur")
pos = geocode(nam)

Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=Chennai&sensor=false
Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=Madurai&sensor=false
Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=Hosur&sensor=false

tlat = pos$lat+0.05    # -- the city name will be above the marker
cities = data.frame(nam, pos$lon,pos$lat,tlat)
names(cities)[2] = "lon"
names(cities)[3] = "lat"
text1 = list("panel.text", cities$lon, cities$tlat, cities$nam,col="red", cex = 0.75)
mark1 = list("panel.points", cities$lon, cities$lat, col="blue")
text2 = list("panel.text",87.0,26.0,"GADM map", col = "dark green", cex = 1.2)
spplot(TN2, "NAME_1",
       sp.layout=list(text1,mark1, text2),
       main="Tamil Nadu Districts",
       colorkey=FALSE, scales=list(draw=TRUE))