在污染指数中包含风速和污染

Question

我想将风速和污染变量一起包含在一个污染like或windrose中。函数windrose()或pollutionrose()只能包含数据序列，而不能包含风速。

图中应以色度表示污染浓度，以风向表示风向，以半径表示风速。

pollutionRose(data, ws="ws", wd="wd", pollutant = "PM1", 
              breaks = c(0, 10, 20, 30, 40, 50, 60, 70))

windRose(data, ws = "ws", wd = "wd", pollutant = "PM1", 
         breaks = c(0, 10, 20, 30, 40, 50, 60, 70))

它应该是什么样子

​

Answer 1

您可以使用openAir中的polarplot函数生成一个图，该图显示了浓度随风速和风向的变化。我在下面为您提供了一个带有输出的代码示例，其中可以用您的输入文件替换Mock_Data。

library(openair)
library(truncnorm)

#Mock data set providing three sites with 1,000 measurements of wind direction
#wind speed and concentrations of a pollutant, in this case, PM10
Site <- c(rep("Site_A",1000), rep("Site_B", 1000), rep("Site_C", 1000))
ws <- rtruncnorm(n=3000, a=0, b=7, mean=2.5, sd=1.75)
wd <- rtruncnorm(n=3000, a=0, b=359, mean=220, sd=50)
PM10 <- rtruncnorm(n=3000, a=0.00, b=7.50, mean=2.5, sd=0.85)  
Mock_Data <- data.frame(Site, ws, wd, PM10)     

#Polar plot function
polarPlot(Mock_Data,pollutant= "PM10",
      type = "Site", col = "jet",
      key.position= "right",key.header= "mean PM10 (ug/m3)", key.footer=NULL)

给出这个输出-

现在，这里要强调的是，使用这个模拟数据集并不是复制该领域实际发生的事情的最佳方式。在那里，出现了流行的风向，风速随上述流行的风向而变化，并且存在污染源，这意味着您可以指出您可能正在调查的污染物浓度较高的特定地区。尽管如此，这只是一个正在工作的函数的示例，并且使用您自己的数据，应该复制您附加的图像。