- Notifications
You must be signed in to change notification settings - Fork6
geanders/weathermetrics
Folders and files
| Name | Name | Last commit message | Last commit date | |
|---|---|---|---|---|
Repository files navigation
Theweathermetrics package provides the following functions to calculate or convert between several weather metrics:
| Weather variable | Function | Input and / or output metric choices |
|---|---|---|
| Temperature | convert_temperature | "kelvin","celsius","fahrenheit" |
| Wind speed | convert_wind_speed | "knots","mph","mps","ftps","kmph" |
| Precipitation | convert_precip | "inches","mm","cm" |
| Dew point temperature | humidity.to.dewpoint | "celsius","fahrenheit" |
| Relative humidity | dewpoint.to.humidity | "celsius","fahrenheit" |
| Heat index | heat.index | "celsius","fahrenheit" |
Algorithms for heat index and wind speed are adapted for R from the algorithms used by the United States National Weather Service'sonline heat index calculator (accessed December 18, 2015) and the National Oceanic and Atmospheric Administration'sonline wind speed conversion (accessed February 22, 2016).
This package includes a function,convert_temperature, that allows you to convert between common temperature measures including degrees Celsius, Fahrenheit, and Kelvin. The following examples show the use of this function for three example datasets:
To convert between degrees Celsius, Fahrenheit, and Kelvin, use theconvert_temperature function:
#Convert from degrees Celsius to degress Fahrenheitdata(lyon)lyon$TemperatureF<- convert_temperature(lyon$TemperatureC,old_metric="celsius",new_metric="fahrenheit")lyon$DewpointF<- convert_temperature(lyon$DewpointC,old_metric="celsius",new_metric="fahrenheit")lyon#> Date TemperatureC DewpointC TemperatureF DewpointF#> 1 2000-06-18 22 13 71.6 55.4#> 2 2000-06-19 24 13 75.2 55.4#> 3 2000-06-20 24 9 75.2 48.2#> 4 2000-06-21 21 10 69.8 50.0#> 5 2000-06-22 23 13 73.4 55.4#> 6 2000-06-23 20 8 68.0 46.4#> 7 2000-06-24 16 11 60.8 51.8#Convert from degrees Fahrenheit to degrees Celsiusdata(norfolk)norfolk$TemperatureC<- convert_temperature(norfolk$TemperatureF,old_metric="f",new_metric="c")norfolk$DewpointC<- convert_temperature(norfolk$DewpointF,old_metric="f",new_metric="c")norfolk#> Date TemperatureF DewpointF TemperatureC DewpointC#> 1 2006-03-12 70 56 21.11 13.33#> 2 2006-03-13 74 59 23.33 15.00#> 3 2006-03-14 63 40 17.22 4.44#> 4 2006-03-15 50 19 10.00 -7.22#> 5 2006-03-16 47 32 8.33 0.00#> 6 2006-03-17 45 39 7.22 3.89#> 7 2006-03-18 44 16 6.67 -8.89#Convert from degrees Kelvin to degrees Celsiusdata(angeles)angeles$TemperatureC<- convert_temperature(angeles$TemperatureK,old_metric="kelvin",new_metric="celsius")angeles$DewpointC<- convert_temperature(angeles$DewpointK,old_metric="kelvin",new_metric="celsius")angeles#> date TemperatureK DewpointK TemperatureC DewpointC#> 1 2010-12-19 288.71 287.04 15.56 13.89#> 2 2010-12-20 287.60 285.37 14.45 12.22#> 3 2010-12-21 286.48 284.82 13.33 11.67#> 4 2010-12-22 287.04 284.82 13.89 11.67#> 5 2010-12-23 286.48 283.15 13.33 10.00#> 6 2010-12-24 287.60 282.04 14.45 8.89#> 7 2010-12-25 285.93 283.15 12.78 10.00
You can specify whether air temperature and dew point temperature inputs are in degrees Celsius, Fahrenheit, or Kelvin using theold_metric andnew_metric options (possible values are'celsius','fahrenheit','kelvin', or'c','f', and'k' for the same). The input forold_metric should be the temperature measure that you want to convert from, and the input fornew_metric should be the temperature measure you wish to convert to.
Theconvert_temperature function is a wrapper function for a variety of individual temperature conversion functions, including:celsius.to.fahrenheit,fahrenheit.to.celsius,celsius.to.kelvin,kelvin.to.celsius,fahrenheit.to.kelvin, andkelvin.to.fahrenheit functions, which you can use individually if you would like.
Theweathermetrics package includes two functions for converting between air temperature, dew point temperature, and relative humidity:dewpoint.to.humidity andhumidity.to.dewpoint.
For example, thelyon data set includes daily values of both air temperature (lyon$TemperatureC) and dew point temperature (lyon$DewpointC). Since this data set includes both air temperature and dew point temperature, you can calculate relative humidity using thedewpoint.to.humidity function:
data(lyon)lyon$RH<- dewpoint.to.humidity(t=lyon$TemperatureC,dp=lyon$DewpointC,temperature.metric="celsius")lyon#> Date TemperatureC DewpointC RH#> 1 2000-06-18 22 13 56.78893#> 2 2000-06-19 24 13 50.28890#> 3 2000-06-20 24 9 38.56789#> 4 2000-06-21 21 10 49.54916#> 5 2000-06-22 23 13 53.43049#> 6 2000-06-23 20 8 46.07966#> 7 2000-06-24 16 11 72.40572
You should specify whether air temperature and dew point temperature inputs are in degrees Fahrenheit or Celsius using thetemperature.metric argument (possible values are'fahrenheit' and'celsius'). If input values for temperature and dew point temperature are in different metrics (i.e., one is in degrees Fahrenheit and the other in degrees Celsius), you must convert one of the inputs using eithercelsius.to.fahrenheit orfahrenheit.to.celsius before you can input the values to thedewpoint.to.humidity function.
As an example of calculating dew point temperature, thenewhaven data set gives daily values of air temperature in degrees Fahrenheit (newhaven$TemperatureF) and relative humidity in % (newhaven$Relative.Humidity). Since this data set includes values for both temperature and relative humidity, you can calculate dew point temperature using thehumidity.to.dewpoint function:
data(newhaven)newhaven$DP<- humidity.to.dewpoint(t=newhaven$TemperatureF,rh=newhaven$Relative.Humidity,temperature.metric="fahrenheit")newhaven#> Date TemperatureF Relative.Humidity DP#> 1 2008-10-19 46 57 31.47#> 2 2008-10-20 48 55 32.45#> 3 2008-10-21 48 68 37.84#> 4 2008-10-22 46 60 32.75#> 5 2008-10-23 44 64 32.48#> 6 2008-10-24 44 67 33.63#> 7 2008-10-25 56 84 51.18
Relative humidity must be input as %, and you must specify the metric of air temperature using thetemperature.metric argument (possible values:'fahrenheit' or'celsius'). The dew point temperature will be calculated using the same metric as the air temperature input to the function. If you wish to get dew point temperature in a different metric than air temperature, you can use theconvert_temperature function. For example:
data(newhaven)newhaven$DP<- humidity.to.dewpoint(t=newhaven$TemperatureF,rh=newhaven$Relative.Humidity,temperature.metric="fahrenheit")newhaven$DP_C<- convert_temperature(newhaven$DP,old_metric="f",new_metric="c")newhaven#> Date TemperatureF Relative.Humidity DP DP_C#> 1 2008-10-19 46 57 31.47 -0.29#> 2 2008-10-20 48 55 32.45 0.25#> 3 2008-10-21 48 68 37.84 3.24#> 4 2008-10-22 46 60 32.75 0.42#> 5 2008-10-23 44 64 32.48 0.27#> 6 2008-10-24 44 67 33.63 0.91#> 7 2008-10-25 56 84 51.18 10.66
Calculations between air temperature, relative humidity, and dew point temperature are based on algorithms used by the United States National Weather Service'sonline heat index calculator (accessed December 18, 2015). These are approximations rather than exact conversions.
Theweathermetrics package includes a function,heat.index, that allows you to calculate a vector of heat index values from vectors of air temperature and either dew point temperature or relative humidity. For example, thesuffolk data set gives daily values of air temperature in degrees Fahrenheit (suffolk$TemperatureF) and relative humidity in % (suffolk$Relative.Humidity) for Suffolk, VA, for the week of July 12, 1998. To calculate daily heat index values for this data set, use theheat.index function:
data(suffolk)suffolk$HI<- heat.index(t=suffolk$TemperatureF,rh=suffolk$Relative.Humidity,temperature.metric="fahrenheit",output.metric="fahrenheit")suffolk#> Date TemperatureF Relative.Humidity HI#> 1 1998-07-12 72 69 72#> 2 1998-07-13 73 66 73#> 3 1998-07-14 74 74 75#> 4 1998-07-15 78 86 80#> 5 1998-07-16 78 100 81#> 6 1998-07-17 80 98 89#> 7 1998-07-18 81 78 86
You must specify whether the air temperature input to the function is in degrees Celsius or Fahrenheit using thetemperature.metric option (possible values:'fahrenheit' or'celsius'). You can choose which metric for heat index output using using theoutput.metric option (the default is to give heat index in the same metric as the air temperature values input to the function).
As another example, thelyon data set gives daily values of air temperature (lyon$TemperatureC) and dew point temperature (lyon$DewpointC), both in degrees Celsius. You can use this data to calculate daily heat index values in degrees Fahrenheit using:
data(lyon)lyon$HI_F<- heat.index(t=lyon$TemperatureC,dp=lyon$DewpointC,temperature.metric="celsius",output.metric="fahrenheit")lyon#> Date TemperatureC DewpointC HI_F#> 1 2000-06-18 22 13 71#> 2 2000-06-19 24 13 75#> 3 2000-06-20 24 9 74#> 4 2000-06-21 21 10 69#> 5 2000-06-22 23 13 73#> 6 2000-06-23 20 8 67#> 7 2000-06-24 16 11 60
When calculating heat index from air temperature and dew point temperature, both must be input in the same metric (either degrees Fahrenheit or degrees Celsius). If this is not the case, you can useconvert_temperature to convert one of the metrics before usingheat.index.
The algorithm for calculating heat index is adapted for R from the algorithms used by the United States National Weather Service'sonline heat index calculator (accessed December 18, 2015). Therefore, results should agree with results from the US National Weather Service online calculator. However, heat index is sometimes calculated using a simpler algorithm. Therefore, heat index values from the function will sometimes differ by one or two degrees compared to other heat index calculators or charts.
Theweathermetrics package includes a function,convert_wind_speed, that allows you to convert wind speed values between several common units of wind speed: knots ('knots'), miles per hour('mph'), meters per second ('mps'), feet per second ('ftps'), and kilometers per hour ('kmph'). The following code shows examples of applying this function to several sample datasets:
data(beijing)beijing$knots<- convert_wind_speed(beijing$kmph,old_metric="kmph",new_metric="knots")beijing#> date TemperatureF MPH mps ftps kmph knots#> 1 2016-01-10 24 4 1.8 5.9 6.5 3.5#> 2 2016-01-11 21 9 4.0 13.2 14.5 7.8#> 3 2016-01-12 18 11 4.9 16.2 17.8 9.6#> 4 2016-01-13 26 11 4.9 16.2 17.8 9.6#> 5 2016-01-14 26 4 1.8 5.9 6.5 3.5#> 6 2016-01-15 28 3 1.3 4.4 4.8 2.6#> 7 2016-01-16 27 5 2.2 7.3 8.0 4.3data(foco)foco$mph<- convert_wind_speed(foco$knots,old_metric="knots",new_metric="mph",round=0)foco$mph<- convert_wind_speed(foco$knots,old_metric="knots",new_metric="mps",round=NULL)foco$kmph<- convert_wind_speed(foco$mph,old_metric="mph",new_metric="kmph")foco#> date TemperatureF knots mph kmph#> 1 2015-10-11 65 4.3 2.212111 3.6#> 2 2015-10-12 59 9.6 4.938666 8.0#> 3 2015-10-13 64 2.6 1.337555 2.2#> 4 2015-10-14 63 2.6 1.337555 2.2#> 5 2015-10-15 54 5.2 2.675111 4.3#> 6 2015-10-16 51 5.2 2.675111 4.3#> 7 2015-10-17 59 3.5 1.800555 2.9
You must specify the unit of wind speed that you wish to convert from using theold_metric option, and the unit of wind speed you wish to convert to using thenew_metric option (possible values:'knots','mph','mps','ftps', or'kmph'). The unit forold_metric cannot be the same as the unit fornew_metric. You can specify the number of decimal places you wish to round to using theround argument. The default value forround is 1, consistant with the algorithms used by the National Oceanic and Atmospheric Administration'sonline wind speed conversion (accessed February 22, 2016). Whenround is set toNULL, the output value will not be rounded.
Theweathermetrics package includes a function,convert_precip, that allows you to convert a vector of precipitation measurement values between inches ('inches'), millimeters ('mm'), and centimeters ('cm'). For example:
data(breck)breck$Precip.mm<- convert_precip(breck$Precip.in,old_metric="inches",new_metric="mm",round=2)breck#> date Precip.in Precip.mm#> 1 2015-06-28 0.04 1.02#> 2 2015-06-29 0.00 0.00#> 3 2015-06-30 0.07 1.78#> 4 2015-07-01 0.00 0.00#> 5 2015-07-02 0.05 1.27#> 6 2015-07-03 0.01 0.25#> 7 2015-07-04 0.01 0.25data(loveland)loveland$Precip.in<- convert_precip(loveland$Precip.mm,old_metric="mm",new_metric="inches",round=NULL)loveland$Precip.cm<- convert_precip(loveland$Precip.mm,old_metric="mm",new_metric="cm",round=3)loveland#> date Precip.mm Precip.in Precip.cm#> 1 2013-09-08 0.00 0.0000000 0.000#> 2 2013-09-09 6.86 0.2700787 0.686#> 3 2013-09-10 12.95 0.5098425 1.295#> 4 2013-09-11 13.21 0.5200787 1.321#> 5 2013-09-12 19.56 0.7700787 1.956#> 6 2013-09-13 4.57 0.1799213 0.457#> 7 2013-09-14 6.60 0.2598425 0.660
You must specify the unit of precipitation measure that you wish to convert from using theold_metric option, and the unit of precipitation measure you wish to convert to using thenew_metric option (possible values:'inches','mm', and'cm'). You can specify the number of decimal places you wish to round to using theround argument. The default value forround is 2. Whenround is set toNULL, the output value will not be rounded.
Calculations between inches and metric units for precipitation measures use the algorithms used by the United States National Weather Service'sMeteorological Conversions (accessed March 20, 2016).
When any of the functions in this package encounter a missing value (NA) within any of the input vectors, the output weather metric for that observation will also be set asNA. For example:
df<-data.frame(T= c(NA,90,85),DP= c(80,NA,70))df$RH<- dewpoint.to.humidity(t=df$T,dp=df$DP,temperature.metric="fahrenheit")df#> T DP RH#> 1 NA 80 NA#> 2 90 NA NA#> 3 85 70 60.88112
Certain values of dew point temperature or relative humidity are impossible. Relative humidity cannot be lower than 0% or higher than 100%. Dew point temperature cannot be higher than air temperature (except in the case of supersaturation) . When any of these functions encounter an impossible weather metric in an input vector, it returnsNA as the output weather metric for that observation. For example:
df<-data.frame(T= c(90,90,85),DP= c(80,95,70))df$heat.index<- heat.index(t=df$T,dp=df$DP,temperature.metric='fahrenheit')#> Warning in dewpoint.to.humidity(t = t, dp = dp, temperature.metric =#> temperature.metric): For some observations, dew point temperature was#> higher than temperature. Since dew point temperature cannot be higher than#> air temperature, relative humidty for these observations was set to 'NA'.df#> T DP heat.index#> 1 90 80 108#> 2 90 95 NA#> 3 85 70 90
Additionally, the function returns a warning to alert the user that the input data includes impossible values for some observations.
All functions have defaults for rounding that are consistent with the algorithms used by the United States National Weather Service's online converters. For several of the functions, you may also specify that outputs are rounded to a different number of digits using theround option. For example:
data(suffolk)suffolk$TempC<- convert_temperature(suffolk$TemperatureF,old_metric="f",new_metric="c",round=5)suffolk$HI<- heat.index(t=suffolk$TemperatureF,rh=suffolk$Relative.Humidity,round=3)suffolk#> Date TemperatureF Relative.Humidity TempC HI#> 1 1998-07-12 72 69 22.22222 72.143#> 2 1998-07-13 73 66 22.77778 73.102#> 3 1998-07-14 74 74 23.33333 74.578#> 4 1998-07-15 78 86 25.55556 80.426#> 5 1998-07-16 78 100 25.55556 80.605#> 6 1998-07-17 80 98 26.66667 88.688#> 7 1998-07-18 81 78 27.22222 86.051
For conversions other than heat index, cite this package as:
G. Brooke Anderson, Roger D. Peng, and Joshua M. Ferreri. 2016.weathermetrics: Functions to Convert Between Weather Metrics. R package version 1.2.2.
To cite this package when calculating the heat index, use:
Anderson GB, Bell ML, Peng RD. 2013. Methods to calculate the heat index as an exposure metric in environmental health research. Environmental Health Perspectives 121(10):1111-1119.
About
Resources
Uh oh!
There was an error while loading.Please reload this page.
Stars
Watchers
Forks
Packages0
Contributors3
Uh oh!
There was an error while loading.Please reload this page.