# load and prepare data
data(winddata)
set40 <- createSet(height=40, v.avg=winddata[,2], dir.avg=winddata[,14])
set30 <- createSet(height=30, v.avg=winddata[,6], dir.avg=winddata[,16])
set20 <- createSet(height=20, v.avg=winddata[,10])
ts <- formatTS(winddata[,1])
neubuerg <- createMast(time.stamp=ts, loc=NULL, desc=NULL,
set40, set30, set20)
neubuerg <- clean(neubuerg)
# plot v vs. dir
plotPolar(neubuerg)
plotPolar(neubuerg, v.set=3, dir.set=2)
# customize plot
plotPolar(neubuerg, col="black", pch=1)
plotPolar(neubuerg, col=2, pch="*")
plotPolar(neubuerg, pch="+", cex=1.2, cex.pts=0.8)
plotPolar(neubuerg, pch="+", cex=0.8, cex.pts=1.2)
Run the code above in your browser using DataLab