Wie interpoliere Region auf Vektorplot?

Question

Weiter von meiner vorherigen Frage 'create vectorplot from velocity dataset'. Ich habe noch 2 Frage, wie man die Figur wie folgt aussehen lässt:

1) Wie die Region interpolieren? Ich habe versucht, interpolate = TRUE verwendet, aber nicht funktioniert.
2) Wie das Pfeilsymbol mit der gleichen Länge definieren (es ist auf den Pfeil zeigt nur bedeuten, die Geschwindigkeitsrichtung)

Hier habe meine Daten uv.nc und Syntax ich geschrieben:

library (raster) 
flname <- 'uv.nc' 
u <- raster(flname, varname = 'U') 
v <- raster(flname, varname = 'V') 
uv <- stack(u,v) 
s <- sqrt(u^2 + v^2) 

library(rasterVis) 
jet <- colorRampPalette(c('#00007F', 'blue', '#007FFF', 'cyan','#7FFF7F', 'yellow', '#FF7F00', 'red', '#7F0000')) 
range = seq(0, 0.5, 0.05) 
vectorplot(uv, isField = 'dXY', interpolate = TRUE, col.regions = jet, region=s, length=0.05) 

Answer 1

(erste Frage) Das interpolate Argument benötigt die Funktion, die von levelplot aufgerufen wird (intern von vectorplot verwendet, um den Hintergrund zu rendern). Dies funktioniert jedoch nicht direkt mit der aktuellen Version rasterVis. Sie können diesen Trick versuchen:

levelplot(s, 
      panel = panel.levelplot.raster, 
      interpolate = TRUE, 
      margin = FALSE) + 
vectorplot(uv, isField = 'dXY', region = FALSE) 

(Zweite Frage) Die Länge der Pfeile durch die Daten bestimmt wird, da Sie dXY = TRUE verwenden. Daher sollten Sie Ihre Daten ändern, um Vektoren mit der gleichen Größe zu erhalten.

uv0 <- uv/s 
vectorplot(uv0, isField = 'dXY', region = s)