# 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], v.std=winddata[,13])
ts <- formatTS(winddata[,1])
neubuerg <- createMast(time.stamp=ts, loc=NULL, desc=NULL,
set40, set30, set20)
neubuerg <- clean(neubuerg)
# create profile based on one height
neubuerg.wp <- profile(neubuerg, v.set=1, dir.set=1) # default alpha=0.2
neubuerg.wp.2 <- profile(neubuerg, v.set=1, dir.set=1, alpha=0.15)
# compare influence of alpha
neubuerg.wp
neubuerg.wp.2
# calculate profiles based on two heights
neubuerg.wp.3 <- profile(neubuerg, v.set=c(1,2), dir.set=1)
neubuerg.wp.4 <- profile(neubuerg, v.set=c(1,3), dir.set=1)
# change number of direction sectors
neubuerg.wp.5 <- profile(neubuerg, v.set=c(1,2), dir.set=1, num.sectors=8)
# compare calculated profiles
neubuerg.wp.3 # 40 m and 30 m dataset
neubuerg.wp.4 # 40 m and 20 m dataset
neubuerg.wp.5
Run the code above in your browser using DataLab