2016-03-23 27 views
1

Bezogen auf meine vorherige Frage über calculate atan2 from two raster object in R?. Es ist möglich, ein Vektordiagramm zu erstellen, in dem die Geschwindigkeit (Neigung) und Richtung (Aspekt) mit einer Reihe von Pfeilen (zB Köcher in Matlab) basierend auf dieser Gleichung angezeigt wird:Wie erstellt man Vektorplot aus Geschwindigkeitsdatensatz?

für Geschwindigkeit: ws < - sqrt (u^2 + v^2)

für Richtung: wd < - (180/pi) * (atan2 (u, v))

my-Datensatz Probe kann

u <- raster(flname, varname = 'U') 
v <- raster(flname, varname = 'V') 
uv <- stack(u,v) 

I verwendet rasterVis downloaded here sein Paket zur Abbildung der Handlung, aber das Ergebnis ist nicht geeignet.

library(rasterVis) 
vectorplot(uv) 

Here the result

+0

Die Datei scheint beschädigt hochgeladen. Vielleicht ist etwas falsch von meinem Ende, nicht sicher. Bitte überprüfen. Die Daten, die ich sehe, sehen wie folgt aus: '~~V¼ €, a¼Ôw,]^k]^k]^k]^k]^k¼Ë» ù: ƒ½Dî "» oŽŽ =: þ? <Þ9j]^k]^k]^k]^k]^k]^k]^k½ ~ Lù; ØØ]^k]^k]^k]^k]^k]^k]^k]^k]^k]^k]^k¾ * 1: ' – JCKaz

Antwort

1

Von der Hilfeseite von vectorplot:

Wenn 'Isfield =' dxy '' 'Objekt' muss ein Raster sein * mit zwei Schichten die horizontale darstellt und die vertikale Komponenten, .

So ist die Lösung:

vectorplot(uv, isField = 'dXY') 
+0

Danke @Oscar Perpiñán, –