drawing order in spplot and Line width in sp maps

Sat, Mar 30, 2013 1:43 PM

On Sat, Mar 30, 2013 at 8:35 AM, Edzer Pebesma <edzer.pebesma at uni-muenster.de> wrote:

Thank you for making these wonderful packages!

This did the trick! I would never ever guess what was missing...

However this breaks colors for _sp.polygons_ if several were used :(
E.g. in my example I had

,----[ Only first color is effective with set_Polypath(FALSE) ]
| l.wshd <- list("sp.polygons", bdry, col=c("magenta", "blue"), lwd=1.5)
`----

So I had to have a logical block in panel function

,----[ snippet from a function supplied as panel ]
|          set_Polypath(FALSE)
|          panel.polygonsplot(...)
|          set_Polypath(TRUE)
|          sp.panel.layout(list(l.wshd, l.streams, l.struc))
`----

I guess it is hard to hit 2 birds with 1 stone...

I'll check it out once it propagates to a binary package @CRAN.

Also I feel like that documentation is misleading a bit regarding _cuts_
saying that it is for SpatialPointsDataFrame only. It appears that it does
change # of breaks for polygon just fine.


P.S. I thought I found an answer to my question about padding when I came
across [1]. However it looks like axis.padding is not for z values...

On 03/30/2013 06:12 AM, Mikhail Titov wrote:

Hi, All!

I know I'm excavating an old thread[1] (and top posting) but that is
exactly what bugs me.

I am surprised to see that I can't pass lwd directly to spplot's main
plot and not extra layouts. I can pass lwd to sp.polygons but not to
spplot directly :( Here is an incomplete snippet showing what I've tried

,----[ How do I change line width of shp polygons supplied to spplot? ]
| spplot(shp, names(d), col="grey80", lwd=5, cuts=5,
| #        par.settings = simpleTheme(lwd=3),
|       par.settings = list(superpose.polygon = list(lwd=5)),
|        panel = function(...) {
| #          panel.polygonsplot(..., lwd=3)
|          panel.polygonsplot(...)
|          sp.panel.layout(list(l.wshd, l.streams, l.struc))
|          if (panel.number()==1) {
|            sp.panel.layout(list(l.north, l.scale.bar, l.scale.text))
|          }
|        },
|        colorkey = list(
|          labels=list(
|            at = labelat,
|            labels = labeltext
|            )
|          ),
|        names.attr=droplevels(unique(rest$case)),
|        col.regions=rev(greens(14, .6))[-1])#rev(heat.colors(100)))
`----

I see that lwd for my layouts apply well, e.g., for l.wshd & l.streams
in the example above. However I can't seem to figure out how to make shp
lines thinner. I do have several variables to plot on a few panels.

Did I miss something?

P.S. Not quite related question but how does levelplot (?) estimates
color scale labels? If I leave colorkey alone, I get misplaced labels
provided cuts=5. For now, I do something like

,----[ Dirty fix for ticks & label alignment ]
| dd.rng <- range(unlist(dd), finite=TRUE)
| delta <- diff(dd.rng)
| labelat <- seq(dd.rng[1] - delta*.07, dd.rng[2] + delta*.07, length.out=7)
| labeltext <- format(exp(labelat), digits=2)
`----

I feel like there is almost a standard function for that but I can't
recall its name. Would you mind reminding if there is any?

Footnotes: 
[1]  https://r-forge.r-project.org/scm/viewvc.php/pkg/R/levelplot.R?view=markup&root=lattice

Mikhail

drawing order in spplot and Line width in sp maps

Thread (3 messages)