风玫瑰图制作，求指点

zhjl

请求大神指点，如何利用风向和风速数据，用什么软件工具可以制作如图所示的风玫瑰图（蒲福氏风级——风玫瑰图）。

MLSHWJZ

http://r.photo.store.qq.com/psb?/ee3f036b-72a7-4dd3-b6b5-2590ba916caa/XRg9xCDWlSVYoVs60FQpd2DidAwHaoEFI67iC28gzlk!/o/dPsAAAAAAAAA&ek=1&kp=1&pt=0&bo=*wKAArcDGgMDAPc!&su=1121381521&tm=1486206000&sce=0-12-12&rf=2-9，
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")

大海牛

excel就可以

liuyw

ncl可以画的

duguannan1985

谢谢大神提示！

zhjl

大海牛 发表于 2017-2-4 09:30
excel就可以

能具体说一下吗？

zhjl

能具体说说吗？谢谢！

MLSHWJZ

                               
登录/注册后可看大图

R语言可以做

MLSHWJZ

效果如图。

zhjl

MLSHWJZ 发表于 2017-2-4 19:11
http://r.photo.store.qq.com/psb?/ee3f036b-72a7-4dd3-b6b5-2590ba916caa/XRg9xCDWlSVYoVs60FQpd2DidAwHao ...

你好，你的程序很好，但是和我要实现的内容不是很一致，麻烦可以详细探讨嘛？