Question

这是一个问题，看看是否有人在他们的旅行中看到过这样的事情。我正在处理大量的天气数据，我想根据wind barbs绘制风。

我已查看了包RadioSonde，但其plotwind()函数未执行我预期的工作。它确实有一个很好的数据类型示例data(ExampleSonde)

可以说我可以将TeachingDemos与my.symbols()结合使用来创建这些风倒钩。如果有人发现（或创造）一种绘制风倒钩的方法，我只是好奇。否则my.symbols()就是。

谢谢，

獾

Answer 1

另一种方法是使用grid图形创建风倒钩。

第一步是计算需要多少，以及需要什么类型的倒钩。正如here所述，我创建了三种类型，分别代表50,10和5节 - 我将速度向下舍入到最接近的5节。

wind_barb下面的函数会为每个风速生成一个新的grob。使用Integrating Grid Graphics Output with Base Graphics Output - Murrell (pg4)中的想法，您可以轻松绘制凹陷并通过旋转视口来表示风向。

一个例子

创建一些数据

set.seed(1)
dat <- data.frame(x=-2:2, y=-2:2, 
                  direction=sample(0:360, 5), 
                  speed=c(10, 15, 50, 75, 100))
#    x  y direction speed
# 1 -2 -2        95    10
# 2 -1 -1       133    15
# 3  0  0       205    50
# 4  1  1       325    75
# 5  2  2        72   100

剧情

library(gridBase)
library(grid)

with(dat, plot(x, y, ylim=c(-3, 3), xlim=c(-3, 3), pch=16))

vps <- baseViewports()
pushViewport(vps$inner, vps$figure, vps$plot)
# Plot
for (i in 1:nrow(dat)) {
    pushViewport(viewport(
                    x=unit(dat$x[i], "native"),
                    y=unit(dat$y[i], "native"), 
                    angle=dat$direction[i]))
        wind_barb(dat$speed[i])
    popViewport()
  }

popViewport(3)

哪个产生

wind_barb创建倒钩的功能（请简化我！）。您可以分别调整mlength和wblength参数来更改倒钩的高度和宽度。

wind_barb <- function(x, mlength=0.1, wblength=0.025) {

  # Calculate which / how many barbs
    # any triangles (50)
    fif <- floor(x /50)
    # and then look for longer lines for remaining speed (10)
    tn <- floor( (x - fif* 50)/10)
    # and then look for shorter lines for remaining speed (5)
    fv <- floor( (x - fif* 50 - tn* 10)/5)

    # Spacing & barb length
    yadj <- 0.5+mlength
    dist <- (yadj-0.5) / 10    
    xadj <- 0.5+wblength
    xfadj <- 0.5+wblength/2        

  # Create grobs
    main_grob <- linesGrob(0.5, c(0.5, yadj ))

    # 50 windspeed
    if(fif != 0) {
      fify <- c(yadj, yadj-dist*seq_len(2* fif) )
      fifx <- c(0.5, xadj)[rep(1:2, length=length(fify))]
      fif_grob <- pathGrob(fifx, fify, gp=gpar(fill="black"))
    } else {
      fif_grob <- NULL
      fify <- yadj+dist
    }

    # Ten windspeed
    if(tn != 0) {
    tny <- lapply(seq_len(tn) , function(x) min(fify) - dist*c(x, x-1))  
    tn_grob <- do.call(gList, 
                      mapply(function(x,y) 
                             linesGrob(x=x, y=y, gp=gpar(fill="black")),
                                      x=list(c(0.5, xadj)), y=tny, SIMPLIFY=FALSE))
    } else {
      tn_grob <- NULL
      tny <- fify
    }

    # Five windspeed
    if(fv != 0) {
    fvy <- lapply(seq_len(fv) , function(x) min(unlist(tny)) -dist* c(x, x-0.5))
    fv_grob <- do.call(gList, 
                        mapply(function(x,y) 
                              linesGrob(x=x, y=y, gp=gpar(fill="black")),
                                      x=list(c(0.5, xfadj)), y=fvy, SIMPLIFY=FALSE))
    } else {
      fv_grob <- NULL
    }    

    # Draw    
    #grid.newpage()
    grid.draw(gList(main_grob, fif_grob, tn_grob, fv_grob))
}

-------------------------------

来自下方的评论

绘制的风向是错误的。要拥有正确的气象风向，请使用angle = 360 - dat$direction[i]。见http://tornado.sfsu.edu/geosciences/classes/m430/Wind/WindDirection.html

在R中绘制风倒钩

1 个答案:

-------------------------------