# 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(time.stamp=winddata[,1])
neubuerg <- createMast(time.stamp=ts, set40, set30, set20)
neubuerg <- clean(mast=neubuerg)
# calculate energy
neubuerg.wb <- weibull(mast=neubuerg, v.set=1, print=FALSE)
neubuerg.e <- energy(wb=neubuerg.wb, print=FALSE)
neubuerg.e.2 <- energy(wb=neubuerg.wb, bins=NULL, print=FALSE) # no binning
# plot energy rose
plotEnergy(energy=neubuerg.e)
plotEnergy(energy=neubuerg.e.2)
# omit total amount
plotEnergy(energy=neubuerg.e, show.total=FALSE)
# change colours
plotEnergy(energy=neubuerg.e, col=gray(5:0 / 5.5))
# change text sizes
plotEnergy(energy=neubuerg.e, cex=0.8)
# manual definition of circles
plotEnergy(energy=neubuerg.e, circles=c(100, 500, 100))
# plot sectors in foreground
plotEnergy(energy=neubuerg.e, fg=TRUE)
# change position of axis labels
plotEnergy(energy=neubuerg.e, pos.axis=135)
# no legend
plotEnergy(energy=neubuerg.e, width.leg=0)
# freaky
plotEnergy(energy=neubuerg.e, border.leg=heat.colors(5), bty.leg="o",
cex.axis=0.5, cex.lab=2, cex.leg=0.5, circles=c(80, 800, 80),
col=rainbow(5), col.axis="green", col.border="orange",
col.circle="purple", col.cross="yellow", col.lab="pink",
col.leg="lightblue", fg=TRUE, lwd.border=2, lwd.circle=3, lwd.cross=4,
lty.circle="52168319", lty.cross="12223242", sec.space=0.6,
title.leg="* WiNd SpEeD *", x.intersp=2, y.intersp=5)
Run the code above in your browser using DataLab