# 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)
# calculate energy
neubuerg.wb <- weibull(neubuerg, 1)
neubuerg.e <- energy(neubuerg.wb)
neubuerg.e.2 <- energy(neubuerg.wb, bins=NULL) # no binning
# plot energy rose
plotEnergy(neubuerg.e)
plotEnergy(neubuerg.e.2)
# omit total amount
plotEnergy(neubuerg.e, show.total=FALSE)
# change colour and text size
plotEnergy(neubuerg.e, col=gray(4:1 / 5), cex=0.9)
plotEnergy(neubuerg.e.2, col="limegreen", cex=1.3)
Run the code above in your browser using DataLab