ambient-local-exporter/pkg/weather/enrich.go

90 lines
2.1 KiB
Go
Raw Normal View History

2025-01-04 22:51:55 +00:00
package weather
2025-01-12 20:30:37 +00:00
import (
"math"
"net/url"
"strconv"
2025-01-12 22:23:32 +00:00
"gitea.libretechconsulting.com/rmcguire/ambient-weather-local-exporter/pkg/ambient/config"
2025-01-12 20:30:37 +00:00
)
2025-01-04 22:51:55 +00:00
// Attempts to complete missing fields that may not
// be set by a specific provider, such as DewPoint and WindChill
2025-01-12 22:23:32 +00:00
func (u *WeatherUpdate) Enrich(weatherStations ...*config.WeatherStation) {
2025-01-05 22:36:39 +00:00
if u == nil {
return
}
2025-01-09 02:28:38 +00:00
if u.WindChillF == nil && u.TempOutdoorF != nil && u.WindSpeedMPH != nil {
wc := CalculateWindChill(*u.TempOutdoorF, *u.WindSpeedMPH)
u.WindChillF = &wc
2025-01-04 22:51:55 +00:00
}
2025-01-09 02:28:38 +00:00
if u.DewPointF == nil && (u.TempOutdoorF != nil && u.HumidityOudoor != nil) {
if *u.TempOutdoorF != 0 || *u.HumidityOudoor != 0 {
dp := CalculateDewPoint(*u.TempOutdoorF, float64(*u.HumidityOudoor))
u.DewPointF = &dp
}
2025-01-04 22:51:55 +00:00
}
2025-01-09 02:28:38 +00:00
if u.BaromAbsoluteIn == nil && u.BaromRelativeIn != nil {
2025-01-04 22:51:55 +00:00
u.BaromAbsoluteIn = u.BaromRelativeIn
}
}
2025-01-05 01:54:36 +00:00
func CalculateDewPoint(tempF, humidity float64) float64 {
2025-01-04 22:51:55 +00:00
// Convert temperature from Fahrenheit to Celsius
tempC := (tempF - 32) * 5 / 9
// Calculate the dew point using the Magnus-Tetens approximation
2025-01-05 01:54:36 +00:00
a := float64(17.27)
b := float64(237.7)
2025-01-04 22:51:55 +00:00
2025-01-05 01:54:36 +00:00
alpha := (a*tempC)/(b+tempC) + math.Log(humidity/100)
2025-01-04 22:51:55 +00:00
dewPointC := (b * alpha) / (a - alpha)
// Convert dew point back to Fahrenheit
dewPointF := (dewPointC * 9 / 5) + 32
return dewPointF
}
2025-01-05 01:54:36 +00:00
func CalculateWindChill(tempF float64, windSpeedMPH float64) float64 {
2025-01-04 22:51:55 +00:00
if windSpeedMPH <= 3 {
// Wind chill calculation doesn't apply for very low wind speeds
return tempF
}
// Formula for calculating wind chill
2025-01-05 01:54:36 +00:00
return 35.74 + 0.6215*tempF -
35.75*math.Pow(windSpeedMPH, 0.16) +
0.4275*tempF*math.Pow(windSpeedMPH, 0.16)
2025-01-04 22:51:55 +00:00
}
2025-01-12 20:30:37 +00:00
// Helper function to set values from fields
// typically from a WeatherUpdate
func SetURLVal(vals *url.Values, key string, value any) {
if value == nil {
return
}
switch v := value.(type) {
case *float64:
if v == nil {
return
}
str := strconv.FormatFloat(*v, 'f', 4, 64)
vals.Set(key, str)
case *int:
if v == nil {
return
}
str := strconv.FormatInt(int64(*v), 10)
vals.Set(key, str)
case *string:
if v == nil {
return
}
vals.Set(key, *v)
}
}