利用ggplot2绘制风向玫瑰图

利用风资源软件画风向图根基上是千篇一律，假如画出新意的风向玫瑰图，可以利用R语言的ggplot2包，画图功效如下：

ggplot2（http://had.co.nz/ggplot2/）是CRAN上知名度很高的画图扩展包，也是R语言利用者最喜欢的画图东西。R语言是统计规模广为利用的开源编程语言，在海外更为风行。假如想相识R语言，可以去统计之都论坛，可能进修入门教程R in Action(http://ishare.iask.sina.com.cn/f/18406061.html)。

附:画图源代码（源数据是NRG Data Retriever导出的txt名目标原始测风数据）

library(ggplot2)
df.raw <- read.table(“winddata.txt”, sep = “\t”, skip = 170, header = T)
# 风向玫瑰图
brea <- c(-1, 11.25 + (22.5*0:16))
winddir <- cut(df.raw$CH7Avg, brea, labels = c(“N”, “NNE”, “NE”, “ENE”, “E”, “ESE”,
“SE”, “SSE”, “S”, “SSW”, “SW”, “WSW”, “W”, “WNW”, “NW”, “NNW”, “N2”))
levels(winddir)[17] = “N”
windspeed <- cut(df.raw$CH1Avg, breaks = c(0, 5, 10, 15, 20))
df.dir <- data.frame(winddir, windspeed)
rosebrea <- 1:5*max(table(winddir))/5
p3 <- ggplot(df.dir, aes(x = winddir, fill = windspeed))
p3 + labs(x = NULL, y = NULL) + geom_bar(aes(y = ..count..))+coord_polar(start = -pi/16) +
scale_y_continuous(breaks = rosebrea, labels = paste(round(rosebrea/length(winddir), 2)*100, “%”))
ggsave(“test.png”)