在污染上升中包括風速和污染

Question

我想將風速和污染變量一起包含在一個污染玫瑰或風玫瑰中。 函數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)

給出這個輸出—— 現在，在這里要強調的是，使用這個模擬數據集並不是復制現場實際發生的事情的最佳方式。 在那里，盛行風向出現，風速因所述盛行風向而異，並且存在污染源，這意味着您可以指出您可能正在調查的污染物濃度較高的特定區域。 盡管如此，這只是工作中的功能的一個示例，並且在您自己的數據存在的情況下，應該復制您附加的圖像。