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

Fri Jan 27 18:44:05 CET 2012

Thanks a lot, everyone - I solved most of my challenges:

the last line of my code should read:
spplot(can2,zcol="groups",col.regions=classcolors,colorkey =
FALSE,lwd=.4,col='white',
  par.settings = list(axis.line = list(col='transparent')))

Thanks a lot for your help!
Dimitri

On Fri, Jan 27, 2012 at 12:00 PM, Dimitri Liakhovitski
<dimitri.liakhovitski at gmail.com> wrote:
> I found an important part of the solution - was able to get rid of the legend!
>
> 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,colorkey =
> FALSE,par.settings =
>    list(axis.line = list(col =  'transparent')))
>
> One last challenge:
> Is it possible to delete the black borders between small geographic
> units and only leave borders between provinces?
> (I need to leave only these borders:
> can1<-getData('GADM', country="CAN", level=1
> plot(can1)  )
> Or maybe one could change the color of the black borders between small
> geographic units?
>
> Thanks a lot!
> Dimitri
>
> On Fri, Jan 27, 2012 at 11:34 AM, Dimitri Liakhovitski
> <dimitri.liakhovitski at gmail.com> wrote:
>> 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
>
>
>
> --
> Dimitri Liakhovitski
> marketfusionanalytics.com

-- 
Dimitri Liakhovitski
marketfusionanalytics.com