Just training and fiddling with spatial analytics
Markdown with R Commander
Ajay Ohri
2015-03-15
> library(maptools)
> library(raster)
> adm <- getData('GADM', country='IND', level=2)
> mahadm=adm[adm$NAME_1=="Maharashtra",]
> head(mahadm,20)
PID ID_0 ISO NAME_0 ID_1 NAME_1 ID_2 NAME_2 NL_NAME_2 VARNAME_2 TYPE_2 ENGTYPE_2
306 17478 105 IND India 21 Maharashtra 306 Ahmednagar Ahmadnagar District District
307 17479 105 IND India 21 Maharashtra 307 Akola District District
308 17480 105 IND India 21 Maharashtra 308 Amravati Amaravati, Amraoti, Amaraoti District District
309 17481 105 IND India 21 Maharashtra 309 Aurangabad District District
310 17482 105 IND India 21 Maharashtra 310 Bhandara District District
311 17483 105 IND India 21 Maharashtra 311 Bid Bir|Beed|Bhir|Bidh District District
312 17484 105 IND India 21 Maharashtra 312 Buldana District District
313 17485 105 IND India 21 Maharashtra 313 Chandrapur Chanda District District
314 17486 105 IND India 21 Maharashtra 314 Dhule Dhulia, West Khandesh District District
315 17487 105 IND India 21 Maharashtra 315 Garhchiroli District District
316 17488 105 IND India 21 Maharashtra 316 Gondiya District District
317 17489 105 IND India 21 Maharashtra 317 Greater Bombay District District
318 17490 105 IND India 21 Maharashtra 318 Hingoli District District
319 17491 105 IND India 21 Maharashtra 319 Jalgaon East Khandesh District District
320 17492 105 IND India 21 Maharashtra 320 Jalna District District
321 17493 105 IND India 21 Maharashtra 321 Kolhapur District District
322 17494 105 IND India 21 Maharashtra 322 Latur Kulaba, Kolaba, Kolabad District District
323 17495 105 IND India 21 Maharashtra 323 Nagpur District District
324 17496 105 IND India 21 Maharashtra 324 Nanded Nander District District
325 17497 105 IND India 21 Maharashtra 325 Nandurbar District District
> mahadm$pop=as.factor(sample(1:10,34,T))
> mahadm$pop2=as.factor(sample(1:10,34,T))
> mahadm$pop3=as.factor(sample(1:10,34,T))
> head(mahadm,20)
PID ID_0 ISO NAME_0 ID_1 NAME_1 ID_2 NAME_2 NL_NAME_2 VARNAME_2 TYPE_2 ENGTYPE_2
306 17478 105 IND India 21 Maharashtra 306 Ahmednagar Ahmadnagar District District
307 17479 105 IND India 21 Maharashtra 307 Akola District District
308 17480 105 IND India 21 Maharashtra 308 Amravati Amaravati, Amraoti, Amaraoti District District
309 17481 105 IND India 21 Maharashtra 309 Aurangabad District District
310 17482 105 IND India 21 Maharashtra 310 Bhandara District District
311 17483 105 IND India 21 Maharashtra 311 Bid Bir|Beed|Bhir|Bidh District District
312 17484 105 IND India 21 Maharashtra 312 Buldana District District
313 17485 105 IND India 21 Maharashtra 313 Chandrapur Chanda District District
314 17486 105 IND India 21 Maharashtra 314 Dhule Dhulia, West Khandesh District District
315 17487 105 IND India 21 Maharashtra 315 Garhchiroli District District
316 17488 105 IND India 21 Maharashtra 316 Gondiya District District
317 17489 105 IND India 21 Maharashtra 317 Greater Bombay District District
318 17490 105 IND India 21 Maharashtra 318 Hingoli District District
319 17491 105 IND India 21 Maharashtra 319 Jalgaon East Khandesh District District
320 17492 105 IND India 21 Maharashtra 320 Jalna District District
321 17493 105 IND India 21 Maharashtra 321 Kolhapur District District
322 17494 105 IND India 21 Maharashtra 322 Latur Kulaba, Kolaba, Kolabad District District
323 17495 105 IND India 21 Maharashtra 323 Nagpur District District
324 17496 105 IND India 21 Maharashtra 324 Nanded Nander District District
325 17497 105 IND India 21 Maharashtra 325 Nandurbar District District
pop pop2 pop3
306 5 2 6
307 8 5 1
308 2 4 4
309 2 6 10
310 8 7 5
311 9 10 1
312 8 9 4
313 4 9 3
314 1 8 7
315 9 5 7
316 7 8 4
317 9 8 6
318 1 10 3
319 1 8 6
320 8 6 10
321 8 4 8
322 4 9 2
323 5 3 5
324 4 8 10
325 1 5 1
> par(mfrow=c(3,1))
> plot(mahadm,col=mahadm$pop)
> plot(mahadm,col=mahadm$pop2)
> plot(mahadm,col=mahadm$pop3)
Like this:
Like Loading...