Hi Oscar, Sorry to bother you again. Please consid

Oh, yes, it's true. Then try this: <div class="snippet-clipboard

Ok, then you have to use as.layer with <code class="n

Thank you, it now works. <a target="_blank" rel="noopener noreferrer

Legend for vectorplot doesn't appear when associated with levelplot about rastervis HOT 7 CLOSED

oscarperpinan commented on September 28, 2024

Legend for vectorplot doesn't appear when associated with levelplot

from rastervis.

Comments (7)

oscarperpinan commented on September 28, 2024 1

Oh, yes, it's true.

Then try this:

pvc <- vectorplot(r1, region = FALSE,
                  key.arrow = list(size = 5, label = 'm/s'))

plv <- levelplot(r1, margin = FALSE)

## reuse the legend of vectorplot in levelplot
plv$legend$top <- pvc$legend$top

plv + pvc

from rastervis.

oscarperpinan commented on September 28, 2024

Hi,

You can use the region argument to display the background. For example:

vectorplot(r1, region = TRUE, key.arrow = list(size = 5, label = 'm/s'))

or with a different Raster:

vectorplot(r1, region = r2, key.arrow = list(size = 5, label = 'm/s'))

from rastervis.

pascaloettli commented on September 28, 2024

Hi,

This only applies if the rasters have the same extent, which is not my case. And with your solution, I think I cannot tune neither the order of the layers nor the colorkey for the region.

Best

from rastervis.

oscarperpinan commented on September 28, 2024

Ok, then you have to use as.layer with under = TRUE:

pvc <- vectorplot(r1, region = FALSE,
                  key.arrow = list(size = 5, label = 'm/s'))

plv <- levelplot(r1, margin = FALSE)

pvc + as.layer(plv, under = TRUE)

from rastervis.

oscarperpinan commented on September 28, 2024

No problem!

Cheers.

from rastervis.

pascaloettli commented on September 28, 2024

Oscar,

Sorry, too fast. No colorkey for levelplot appears.

library(rasterVis)

proj <- CRS('+proj=longlat +datum=WGS84')

df <- expand.grid(x=seq(-2, 2, .01), y=seq(-2, 2, .01))
df$z <- with(df, (3*x^2 + y)*exp(-x^2-y^2))
r1 <- rasterFromXYZ(df, crs=proj)
df$z <- with(df, x*exp(-x^2-y^2))
r2 <- rasterFromXYZ(df, crs=proj)

pvc <- vectorplot(r1, region = FALSE, key.arrow = list(size = 5, label = 'm/s'))

plv <- levelplot(r1, margin = FALSE)

pvc + as.layer(plv, under = TRUE)