Inlcude windspeed and pollution in pollutionrose

I would like to include the wind speed and the pollution variable together in one pollutionrose or windrose. The function windrose() or pollutionrose() only enables to include data series instead of wind speed.

The graph should show the pollution concentration by coluor, the wind direction by direction and the wind speed by radius like in the image.

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))

How it should looks like Image may be NSFW.
Clik here to view.