.Rhistory

read.table("airport_coordinates.txt",sep=":")
read.table("airport_coordinates.txt",sep=":")
flights <- read.table("flightdata_zrh")
flights <- read.table("flightdata_zrh.csv")
flights <- read.table("flightdata_zrh.csv",sep=";")
flights <- read.table("flightdata_zrh.csv",sep=";")
flights <- read.table("flightdata_zrh.csv",sep=";")
flights <- read.table("flightdata_zrh.csv",sep=";")
airports <- read.table("airport_coordinates.txt",sep=":")
View(flights)
flights <- read.table("flightdata_zrh.csv",sep=";",header=TRUE)
airports <- read.table("airport_coordinates.txt",sep=":",header=TRUE)
flights %>% left_join(aiports, by=c("origin_destination"="iata"))
library(tidyverse)
flights %>% left_join(aiports, by=c("origin_destination"="iata"))
flights %>% left_join(airports, by=c("origin_destination"="iata"))
flights_airports <- flights %>% left_join(airports, by=c("origin_destination"="iata"))
View(flights_airports)
flights_airports2 <- flights %>% left_join(airports, by=c("origin_destination"="icao"))
View(flights)
View(airports)
flights_airports2 %>% filter(is.na(icao))
flights_airports2 <- flights %>% left_join(airports, by=c("origin_destination"="icao"))
summary(flights_airports2)
summary(flights_airports)
flights_airports %>% group_by(airports) %>% summarize(n=n())
flights_airports %>%
group_by(origin_destination) %>%
summarize(n=n())
airports_n <- flights_airports %>%
group_by(origin_destination) %>%
summarize(n=n())
View(airports_n)
airports_n <- flights_airports %>%
group_by(origin_destination) %>%
summarize(n=n())
View(airports_n)
#airport database : http://www.partow.net/miscellaneous/airportdatabase/
airports <- read.table("airports.csv",sep=";",header=TRUE)
#airport database : http://www.partow.net/miscellaneous/airportdatabase/
airports <- read.csv("airports.csv",header=TRUE)
#airport database : http://www.partow.net/miscellaneous/airportdatabase/
airports <- read.csv("airports.csv",sep=";",header=TRUE)
View(airports)
#airport database : http://www.partow.net/miscellaneous/airportdatabase/
airports <- read.csv("airports.csv",sep=",",header=TRUE)
#airport database : http://www.partow.net/miscellaneous/airportdatabase/
airports <- read.csv("airports.csv",sep=";")
View(airports)
#airport database : http://www.partow.net/miscellaneous/airportdatabase/
airports <- read.csv("airports.csv",sep=",")
View(airports)
#airport database : http://www.partow.net/miscellaneous/airportdatabase/
airports <- read.csv("airports.csv")
#airport database : http://www.partow.net/miscellaneous/airportdatabase/
airports <- read.csv("airports.csv")
View(airports)
# airport data for 151657 flights not available?
flights_airports <- flights %>% left_join(airports, by=c("origin_destination"="iata_code"))
summary(fligths_airports)
summary(fligths_airports)
summary(flights_airports)
airports %>%
group_by(iata_code) %>%
count()
airports %>%
group_by(iata_code) %>%
count() %>%
count()
iata_red <-airports %>%
group_by(iata_code) %>%
count()
View(iata_red)
View(flights_airports)
flights <- read.table("flightdata_zrh.csv",sep=";",header=TRUE) %>% mutate(id=nrow())
flights <- read.table("flightdata_zrh.csv",sep=";",header=TRUE) %>%
mutate(id=row_number())
# airport data for 71 flights not available? redundat data for ca.14'000 flights?
flights_airports <- flights %>% left_join(airports, by=c("origin_destination"="iata_code"))
flights_airports %>% group_by(id) %>% summarize(n=n())
# airport data for 71 flights not available? redundat data for ca.14'000 flights?
flights_airports <- flights %>% left_join(airports, by=c("origin_destination"="iata_code"))
flights <- read.table("flightdata_zrh.csv",sep=";",header=TRUE) %>%
mutate(rowid=row_number())
# airport data for 71 flights not available? redundat data for ca.14'000 flights?
flights_airports <- flights %>% left_join(airports, by=c("origin_destination"="iata_code"))
flights_airports %>%
group_by(rowid) %>%
summarize(n=n())
flights_airports %>%
group_by(rowid) %>%
summarize(n=n()) %>%
arrange(desc(n))
duplicates <- flights_airports %>%
group_by(rowid) %>%
mutate(n=n()) %>%
arrange(desc(n))
View(duplicates)
View(flights_airports)
#airport database : http://ourairports.com/ -> open aiport data http://ourairports.com/data/
airports <- read.csv("airports.csv") %>% filter(iata_code!="small_airport")
# airport data for 71 flights not available? redundat data for ca.14'000 flights?
flights_airports <- flights %>% left_join(airports, by=c("origin_destination"="iata_code"))
airports <- read.csv("airports.csv") %>%
filter(iata_code!="small_airport")
#airport database : http://ourairports.com/ -> open aiport data http://ourairports.com/data/
airports <- read.csv("airports.csv") %>%
filter(type!="small_airport")
# airport data for 71 flights not available? redundat data for ca.14'000 flights?
flights_airports <- flights %>% left_join(airports, by=c("origin_destination"="iata_code"))
airports_n <- flights_airports %>%
group_by(origin_destination) %>%
summarize(n=n())
iata_red <-airports %>%
group_by(iata_code) %>%
count()
View(iata_red)
duplicates <- flights_airports %>%
group_by(rowid) %>%
mutate(n=n()) %>%
arrange(desc(n))
View(duplicates)
airports <- read.csv("airports.csv") %>%
#filter small airports as ther appears to be a minor airport with the same iatacode as CDG (paris)
filter(type!%in% c("small_airport","closed"))
#airport database : http://ourairports.com/ -> open aiport data http://ourairports.com/data/
airports <- read.csv("airports.csv") %>%
#filter small airports as ther appears to be a minor airport with the same iatacode as CDG (paris)
filter(type %in% !c("small_airport","closed"))
#airport database : http://ourairports.com/ -> open aiport data http://ourairports.com/data/
airports <- read.csv("airports.csv") %>%
#filter small airports as ther appears to be a minor airport with the same iatacode as CDG (paris)
filter(!type %in% c("small_airport","closed"))
# airport data for 71 flights not available? redundat data for ca.14'000 flights?
flights_airports <- flights %>%
left_join(airports, by=c("origin_destination"="iata_code"))
#redundant flights? double airports?
duplicates <- flights_airports %>%
group_by(rowid) %>%
mutate(n=n()) %>%
arrange(desc(n))
duplicates <- flights_airports %>%
group_by(rowid) %>%
mutate(n=n()) %>%
arrange(desc(n))
View(duplicates)
airports <- read.csv("airports.csv") %>%
#filter small airports as ther appears to be a minor airport with the same iatacode as CDG (paris) or were closed and replaced by new airports (e.g. Hongkong)
filter(!type %in% c("small_airport","closed") | !is.na(latitude_deg))
flights_airports <- flights %>%
left_join(airports, by=c("origin_destination"="iata_code"))
#airport database : http://ourairports.com/ -> open aiport data http://ourairports.com/data/
airports <- read.csv("airports.csv") %>%
#filter small airports as ther appears to be a minor airport with the same iatacode as CDG (paris) or were closed and replaced by new airports (e.g. Hongkong)
filter(!type %in% c("small_airport","closed"))
#airport database : http://ourairports.com/ -> open aiport data http://ourairports.com/data/
airports <- read.csv("airports.csv") %>%
#filter small airports as ther appears to be a minor airport with the same iatacode as CDG (paris) or were closed and replaced by new airports (e.g. Hongkong)
filter(!type %in% c("small_airport","closed")) %>%
filter(is.na(!latitude_deg))
#airport database : http://ourairports.com/ -> open aiport data http://ourairports.com/data/
airports <- read.csv("airports.csv") %>%
#filter small airports as ther appears to be a minor airport with the same iatacode as CDG (paris) or were closed and replaced by new airports (e.g. Hongkong)
filter(!type %in% c("small_airport","closed")) %>%
filter(!is.na(latitude_deg))
# airport data for 71 flights not available? redundat data for ca.14'000 flights?
flights_airports <- flights %>%
left_join(airports, by=c("origin_destination"="iata_code"))
airports1 <- read.csv("airports.csv") %>%
#filter small airports as ther appears to be a minor airport with the same iatacode as CDG (paris) or were closed and replaced by new airports (e.g. Hongkong)
filter(!type %in% c("small_airport","closed"))
airports1 <- read.csv("airports.csv") %>%
#filter small airports as ther appears to be a minor airport with the same iatacode as CDG (paris) or were closed and replaced by new airports (e.g. Hongkong)
filter(!type %in% c("small_airport","closed")) %>%
filter(!is.na(latitude_deg))
View(airports1)
#airport database : http://ourairports.com/ -> open aiport data http://ourairports.com/data/
airports1 <- read.csv("airports.csv") %>%
#filter small airports : one small airports has the same iatacode as CDG (paris), others were closed and replaced by new airports (e.g. Hongkong) and are redundant
filter(!type %in% c("small_airport","closed")) %>%
filter(!is.na(latitude_deg))
airports <- read.csv("airports.csv") %>%
#filter small airports : one small airports has the same iatacode as CDG (paris), others were closed and replaced by new airports (e.g. Hongkong) and are redundant
filter(!type %in% c("small_airport","closed")) %>%
filter(!is.na(latitude_deg))
# airport data for 71 flights not available? redundat data for ca.14'000 flights?
flights_airports <- flights %>%
left_join(airports, by=c("origin_destination"="iata_code"))
View(duplicates)
airports <- read.csv("airports.csv") %>%
#filter small airports : one small airports has the same iatacode as CDG (paris), others were closed and replaced by new airports (e.g. Hongkong) and are redundant
filter(!type %in% c("small_airport","closed")) %>%
filter(latitude_deg==0 & longitude_deg==0)
View(airports)
airports <- read.csv("airports.csv") %>%
#filter small airports : one small airports has the same iatacode as CDG (paris), others were closed and replaced by new airports (e.g. Hongkong) and are redundant
filter(!type %in% c("small_airport","closed")) %>%
#filter airports with missing geocoordinates
filter(latitude_deg!=0 & longitude_deg!=0)
# airport data for 71 flights not available? redundat data for ca.14'000 flights?
flights_airports3 <- flights %>%
left_join(airports, by=c("origin_destination"="iata_code"))
duplicates <- flights_airports3 %>%
group_by(rowid) %>%
mutate(n=n()) %>%
arrange(desc(n))
View(duplicates)
airports <- read.csv("airports.csv") %>%
#filter small airports : one small airports has the same iatacode as CDG (paris), others were closed and replaced by new airports (e.g. Hongkong) and are redundant
filter(!type %in% c("small_airport","closed")) %>%
#filter airports with missing geocoordinates
filter(latitude_deg!=0 & longitude_deg!=0) %>%
#filter duplicate OR Tambo International Airport
filter(id!=31055)
# airport data for 71 flights not available? redundat data for ca.14'000 flights?
flights_airports <- flights %>%
left_join(airports, by=c("origin_destination"="iata_code"))
summarize(flights_airports)
summarise(flights_airports)
summary(flights_airports)
flights_airports$effective_time
hms(flights_airports$effective_time[1])
lubridate::hms(flights_airports$effective_time[1])
flights_airports$effective_time[1]
flights_airports$effective_time[1,]
flights_airports <- flights %>%
left_join(airports, by=c("origin_destination"="iata_code")) %>%
#generate time difference - planed vs effective
mutate(timediff=lubridate::hms(flights_airports$planed_time)-lubridate::hms(flights_airports$effective_time))
View(flights_airports)
# airport data for 71 flights not available?
flights_airports <- flights %>%
left_join(airports, by=c("origin_destination"="iata_code")) %>%
#generate time difference - planed vs effective
mutate(difference=difftime(effective_time,planed_time,units="minutes") )
difftime(flights$effective_time[1]-flights$planed_time[1])
difftime(flights$effective_time[1],flights$planed_time[1])
flights$effective_time[1]
flights_airports <- flights %>%
left_join(airports, by=c("origin_destination"="iata_code")) %>%
#generate time difference - planed vs effective
mutate(effective_time1=lubridate::hms(effective_time),
planed_time1=lubridate::hms(planed_time))
difftime(flights$effective_time1[1],flights$planed_time1[1])
flights$effective_time1[1]
difftime(flights_airports$effective_time1[1],flights_airports$planed_time1[1])
as.period(flights_airports$effective_time1[1]-flights_airports$planed_time1[1])
rm(duplicates)
rm(flights)
rm(flights_airports2)
flight_sf <- readRDS("L:/STAT/08_DS/06_Diffusion/OGD/data_publication/twist_zrh/flight_sf.RDS")
View(flight_sf)
plot(flight_sf)
plot(flight_sf)
plot(flight_sf$geometry)
plot(flight_sf)
View(flight_sf)
flight_sf %>% filter(!is.na(geometry)) %>% plot()
library(tidyverse)
flight_sf %>%
filter(!is.na(geometry)) %>%
plot()
flight_sf %>%
filter(!is.na(geometry)) %>%
plot(geometry)
ggplot(flight_sf)+
geom_sf()