idshklein

pacman::p_load(tidyverse,sf,sfnetworks,nngeo,tidygraph)
options(scipen = 10000)
outer = matrix(c(0,0,1000,0,1000,1000,0,1000,0,0),ncol=2, byrow=TRUE)
make_grid <- function(street_length){st_polygon(list(outer)) %>% 
    st_make_grid(street_length) %>% 
    st_segments() }
x <- 200
map_df(c(2,4,5,8,10,20),~make_grid(1000/.x) %>% 
         st_segmentize(20) %>% 
         st_segments() %>% 

idshklein

pacman::p_load(tidyverse,jsonlite,lubridate,rayshader)
df_2019 <- read_json("https://data.gov.il/api/3/action/datastore_search?resource_id=5ed3ac9d-0709-46fb-aa42-1156682f3b85&limit=389&filters={%22OfficeLineId%22:%2224100%22}")
df_2020 <- read_json("https://data.gov.il/api/3/action/datastore_search?resource_id=6042d958-d86f-40bc-803b-2c5620d90394&limit=392&filters={%22OfficeLineId%22:%2224100%22}")
df_2021 <- read_json("https://data.gov.il/api/3/action/datastore_search?resource_id=09de377b-6462-48d5-8470-177c17016f6e&limit=405&filters={%22OfficeLineId%22:%2224100%22}")
df_2022 <- read_json("https://data.gov.il/api/3/action/datastore_search?resource_id=bfd22231-7a8b-4fc0-91bd-342397860ab5&limit=440&filters={%22OfficeLineId%22:%2224100%22}")
dfis <- function(df){df$result$records %>% map_df(~.x)}
dfis(df_2019) %>% 
  bind_rows(dfis(df_2020)) %>% 
  bind_rows(dfis(df_2021)) %>% 
  mutate_all(as.character) %>% 

idshklein

library(sf)
library(tidyverse)
library(sfnetworks)
library(tidygraph)
library(lwgeom)
library(mapview)
library(igraph)
library(qgisprocess)
library(nngeo)
library(dbscan)

idshklein

pacman::p_load(tidyverse,jsonlite,httr)
buses <- read_json("https://data.gov.il/api/3/action/datastore_search?resource_id=91d298ed-a260-4f93-9d50-d5e3c5b82ce1&limit=20000")
buses1 <- buses$result$records %>% map_df(~.x)
rishui <- read_json("https://data.gov.il/api/3/action/datastore_search?resource_id=8fd3aa6c-a64f-4c91-818e-fd35017cd19b&&filters={%22rishui_date%22:%222022-08-09T00:00:00%22}&limit=20000")
rishui1 <- rishui$result$records%>% map_df(~.x)
GET("https://open-bus-stride-api.hasadna.org.il/siri_rides/list?get_count=false&scheduled_start_time_from=2022-08-09T00%3A00%3A00%2B02%3A00&scheduled_start_time_to=2022-08-09T23%3A59%3A59%2B02%3A00&order_by=siri_route_id%20asc%2Cvehicle_ref%20desc&limit=200000",write_disk(siri <- tempfile(fileext = ".json")))
siri1 <- read_json(siri)%>% map_df(~.x)
rishui1$trips_count %>% sum()
GET("https://open-bus-stride-api.hasadna.org.il/gtfs_rides/list?get_count=false&limit=100000",write_disk(rides1 <- tempfile(fileext = ".json")))
GET("https://open-bus-stride-api.hasadna.

idshklein

pacman::p_load(tidyverse,sf,sfdep,sfnetworks,waywiser,tidygraph)
nc = st_read(system.file("shape/nc.shp", package="sf"))
# waywiser- only polygons
waywiser::ww_build_neighbors(nc)
# sfnetworks - neighboring points connected by edges
st_intersects(nc,nc) %>% 
  imap_dfr(~data.frame(from=.y,to=.x)) %>%
  filter(from!=to) %>% 
  as_tbl_graph() %>% 
  left_join(st_centroid(nc) %>% mutate(rn=row_number() %>% as.character()),by=c("name"="rn")) %>% 

idshklein

pacman::p_load(tidyverse,jsonlite,lubridate)
rishui <- read_json("https://data.gov.il/api/3/action/datastore_search?resource_id=8fd3aa6c-a64f-4c91-818e-fd35017cd19b&limit=300000")
rishui1 <- read_json("https://data.gov.il/api/3/action/datastore_search?resource_id=8fd3aa6c-a64f-4c91-818e-fd35017cd19b&offset=300000&limit=300000")
rishui2 <- read_json("https://data.gov.il/api/3/action/datastore_search?resource_id=8fd3aa6c-a64f-4c91-818e-fd35017cd19b&offset=600000&limit=300000")
rishui3 <- read_json("https://data.gov.il/api/3/action/datastore_search?resource_id=8fd3aa6c-a64f-4c91-818e-fd35017cd19b&offset=900000&limit=300000")
rishhi_all <- list(rishui,rishui1,rishui2,rishui3)
tbl <- rishhi_all %>% 
  map_df(~.x$result$records %>% map_df(~.x))

tbl %>% glimpse()

idshklein

library(tidyverse)
library(jsonlite)
library(lubridate)
library(hms)
res3 <- map(0:25,function(x){
  print(x)
  res2 <- read_json(paste0("https://open-bus-stride-api.hasadna.org.il/siri_rides/list?siri_route__line_refs=7005&order_by=siri_route_id%20asc%2Cvehicle_ref%20desc&offset=",x*100))
  df2 <- map_df(1:length(res2),~res2[[.x]])  
  df2
})

idshklein

library(tidyverse)
library(jsonlite)
library(lubridate)
library(sf)
res3 <- map(0:12,function(x){
  print(x)
  res2 <- read_json(paste0("https://open-bus-stride-api.hasadna.org.il/siri_vehicle_locations/list?lon__greater_or_equal=35.1998&lon__lower_or_equal=35.201&lat__greater_or_equal=31.754&lat__lower_or_equal=31.7558&offset=",1+x*100,"&recorded_at_time_from=2022-03-09T00:00:00.000Z&recorded_at_time_to=2022-03-20T00:00:00.000Z"))
  df2 <- map_df(1:length(res2),~res2[[.x]])  
  df2
})

idshklein

library(tidyverse)
library(sf)
library(sfnetworks)
library(tidygraph)
library(dbscan)
stations <- st_read("METRO_STAT.shp")
lines <- st_read("METRO_LINE.shp")
plot(stations$geometry)
clus <- st_coordinates(stations) %>% 
  dbscan(eps = 300,minPts = 1) %>% 

idshklein

library(tidyverse)
library(jsonlite)
library(lubridate)
library(hms)
library(sf)
library(ggspatial)
library(gtfstools)
library(sfnetworks)
library(ggimage)
library(imputeTS)