[R-sig-Geo] Coloring Canada provinces in raster

[Dimitri Liakhovitski]

Vielen-vielen Dank, Rafael.
A factor!
That does exactly what I wanted! Thanks a million!

Still can't figure out how to get rid of the legend on the right. I
couldn't find the solution here:
http://127.0.0.1:14681/library/sp/html/spplot.html
Or, if it's impossible, maybe one could at least change the order in
the legend (to have 1 at the top and 6 at the bottom)?
Finally- anyone knows how to delete the borders between small
geographic units and only leave borders between provinces?

Thank you very much, everyone for your help!
Dimitri

The working code so far:

library(raster)
can2<-getData('GADM', country="CAN", level=2)
nrofunits<-length(can2$NAME_2)
groups<-c(rep(1:6,(nrofunits %/%
6)),rep(1,5))[order(c(rep(1:6,(nrofunits %/% 6)),rep(1,5)))]
# adding values (grouping values) to "data":
can2 at data[["groups"]]<-groups
can2 at data$groups<-as.factor(can2 at data$groups)
classcolors <- rainbow(6)
spplot(can2,zcol="groups",col.regions=classcolors)



On Fri, Jan 27, 2012 at 3:25 AM, Rafael Wüest <rafael.wueest at wsl.ch> wrote:
> Hi Dimitri
>
> you can use the col.regions argument in case groups is a factor:
>
> can2 at data$groups<-as.factor(can2 at data$groups)
> spplot(can2,zcol="groups",col.regions=c('red','blue','yellow','lightblue','green','orange'))
>
> Hope this does what you were looking for...
>
> Rafi
>
> On 26.01.2012, at 21:42, Dimitri Liakhovitski wrote:
>
>> Dear map-gurus,
>> I've made some progress based on one post I found, but now I am stuck.
>> I was able to group Canadian admin units into 6 groups so that each
>> group has the same color:
>>
>> # Getting the data:
>> library(raster)
>> can2<-getData('GADM', country="CAN", level=2)
>> class(can2)
>>
>> # Creating my (artificial) groupings of Canadian admin units:
>> nrofunits<-length(can2$NAME_2)
>> groups<-c(rep(1:6,(nrofunits %/%
>> 6)),rep(1,5))[order(c(rep(1:6,(nrofunits %/% 6)),rep(1,5)))]
>>
>> # adding values (grouping values) to "data" and calling this new
>> variables "groups"
>> can2 at data[["groups"]]<-groups
>>
>> # Plotting:
>> spplot(can2,"groups")
>>
>> Questions that I am still not sure how to address:
>> How can I use MY colors instead of the default colors? In fact, I
>> don't want heat map, I need 6 separate colors...
>> Can I remove the legend from the plot?
>> Can I eliminate borders between counties but have borders between
>> Canadian Provinces (like in the US example below)?
>> I.e., only borders from the plot for can1:
>> can1<-getData('GADM', country="CAN", level=1)
>> plot(can1)
>>
>> Thanks a lot!
>> Dimitri
>>
>>
>>
>> On Thu, Jan 26, 2012 at 1:14 PM, Dimitri Liakhovitski
>> <dimitri.liakhovitski at gmail.com> wrote:
>>> Dear R map gurus,
>>> I am very new to creating maps in R and have never used "raster"
>>> before, so very sorry if it's a simple question.
>>>
>>> I am trying to create a map of Canada with small administrative units
>>> (counties and the like) colored based on some criterion I set (e.g.,
>>> all units that belong to Group 1 get one color, to Group 2 get another
>>> color, etc). At the bottom of this e-mail is an example of what I am
>>> trying to achieve - but it's for the US (based on coloring US
>>> counties) and using "maps". I was not able to find small
>>> administrative divisions (like counties) for Canada in "maps", only
>>> cities.
>>> So, I was referred to "raster". I was able to read in the Canada data
>>> at the right administrative level:
>>>
>>> library(raster)
>>> can2<-getData('GADM', country="CAN", level=2)
>>> class(can2)
>>>
>>> # "NAME_2" in object "can2" contains exactly the level of
>>> administrative detail I need!
>>> can2$NAME_2
>>>
>>> # and plot and ssplot plot the map with those counties and colors them:
>>> plot(can2)     # No colors
>>> spplot(can2, "NAME_2")  # colored, but not sure based on what rule +
>>> creates a legend
>>>
>>> Question:
>>> Is it at all possible to color the units in can2$NAME_2 based on their
>>> groupings that I create (like in the US example below), and the colors
>>> I choose?
>>> Is it possible to remove the legend from the plot?
>>> Is it possible to have no borders between counties but have borders
>>> between Canadian Provinces (like in the US example below)?
>>>
>>> Thank you very much for your help!
>>> Dimitri Liakhovitski
>>>
>>>
>>> ### Example for creating a US map with colored counties using "maps":
>>>
>>> # 1. Creating an (ARTIFICIAL) criterion for coloring US counties in
>>> only few colors:
>>> library(maps)
>>> allcounties<-data.frame(county=map('county', plot=FALSE)$names)
>>> allcounties$group<-c(rep(1:6,513),rep(1,4))[order(c(rep(1:6,513),rep(1,4)))]
>>> ### My colors:
>>> mycolors <- rainbow(6)
>>> map('county',fill=TRUE,col=mycolors[allcounties$group],resolution=0,lty=0,bg=
>>> "transparent")
>>> map('state', lwd=1, add=TRUE)     # adding the state borders
>>
>>
>>
>> --
>> Dimitri Liakhovitski
>> marketfusionanalytics.com
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>
>
> --
> Rafael Wüest
> Swiss Federal Research Institute WSL
> Zürcherstrasse 111
> 8903 Birmensdorf
> Switzerland
>
> +41 44 7392126
> rafael.wueest at wsl.ch
> http://www.wsl.ch/info/mitarbeitende/wueest
>
>
>



-- 
Dimitri Liakhovitski
marketfusionanalytics.com