Re: [R-sig-Geo] European Union NUTS2, NUTS3 regions (SHP format)

[Tomas Zelinsky]

Thank you Barry for your help. It works now. I'm just still not sure how 
to use different widths of lines for different statistical levels. 
Hopefully I'll figure it out soon.

Thanks again.

Tomas




Barry Rowlingson  wrote / napísal(a):
2010/8/27 Tomas Zelinsky <tomas.zelin...@tuke.sk>:
  
Hello,

you probably know administrative maps of the EU and their availability (both
SHP and GDB) via
http://epp.eurostat.ec.europa.eu/portal/page/portal/gisco/popups/references/administrative_units_statistical_units_1
.

Basically there are three administrative levels in one map (NUTS 1, 2 and 3
levels).

1. Is there any possibility how to make NUTS3 level borderlines invisible if
my aim is to e.g. analyze NUTS2 level?

2. Or is it possible to make e.g. NUTS1 level borderlines thicker than the
other?
    

 Looking at NUTS_RG_60M_2006.SHP there's a STAT_LEVL_ variable which I
guess is the NUTS level code - this is probably documented in the
meta-data somewhere. So you can do:

 require(rgdal)
 s = readOGR(".","NUTS_RG_60M_2006")
 plot(s[s$STAT_LEVL_==0,])
 plot(s[s$STAT_LEVL_==1,])

 and so on. You might also want to extract just the polygons for a
particular level into a new object:

 n2 = s[s$STAT_LEVL_==2,]

 You can use lwd=X in plot for sp objects to change the line width.

Barry


_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo