analysis/wind_flow.R

In khufkens/swiftr: Process bird (swift) GPS position tracker data and meta-data.

library(gdistance)
library(rWind)
library(swiftr)

r <- wind_direction("~/download.nc")
names(r) <- c("wind.speed","wind.direction")

loc <- data.frame(lat = seq(44,47,0.25))
loc$lon <- 77


# read position data, remove 0 values
df <- sr_read_pos("data-raw/Obs010617_050807_Tag13819 - Copy - Combined.pos")
df <- df[which(df$latitude != 0),]
bla <- df[88:92,]

r$wind.direction <- 180 + r$wind.direction

Conductance<-flow.dispersion(r)

AtoB<- shortestPath(Conductance,
                    c(bla$lon[1], bla$lat[1]),
                    c(bla$lon[3], bla$lat[3]), output="SpatialLines")
plot(r$wind.speed)
points(bla$lon, bla$lat)
lines(AtoB)

apply(loc, 1, function(l){
  lines(shortestPath(Conductance,
               c(bla$lon[1], bla$lat[1]),
               c(l['lon'], l['lat']), output="SpatialLines"))
})


# combine everything in an mapview visualization
m <- mapview(line) +
  mapview(test, zcol = "direction", popup = popupTable(test))

m1 <- mapview(r)

sync(m,m1)