87 lines
1.9 KiB
Plaintext
87 lines
1.9 KiB
Plaintext
---
|
|
title: "Topic 8"
|
|
output: html_document
|
|
date: "`r Sys.Date()`"
|
|
---
|
|
|
|
```{r setup, include=FALSE}
|
|
knitr::opts_chunk$set(echo = TRUE)
|
|
```
|
|
|
|
# Download data
|
|
from OpenSky https://opensky-network.org/datasets/states/#flights/
|
|
```{r download}
|
|
library(dplyr)
|
|
library(lubridate)
|
|
library(readr)
|
|
library(utils)
|
|
|
|
# -------------------------------
|
|
# Step 1: Define dates
|
|
# -------------------------------
|
|
dates <- seq(as.Date("2022-09-01"), as.Date("2022-09-05"), by = "days")
|
|
|
|
tmp_dir <- tempdir()
|
|
all_flights <- list()
|
|
|
|
# -------------------------------
|
|
# Step 2: Download each day's CSV
|
|
# -------------------------------
|
|
for(d in dates) {
|
|
d <- as.Date(d)
|
|
|
|
# korrekter Dateiname
|
|
file_name <- paste0("flight_sample_", d, ".csv.gz")
|
|
url <- paste0("https://s3.opensky-network.org/data-samples/flights/", file_name)
|
|
|
|
message("Downloading: ", url)
|
|
|
|
dest <- file.path(tmp_dir, file_name)
|
|
|
|
tryCatch({
|
|
download.file(url, dest, mode = "wb")
|
|
|
|
# CSV lesen
|
|
df <- read_csv(dest, col_types = cols(.default = "c"))
|
|
all_flights[[length(all_flights)+1]] <- df
|
|
}, error = function(e) {
|
|
warning("Failed to download ", file_name)
|
|
})
|
|
}
|
|
|
|
# -------------------------------
|
|
# Step 3: Combine all days
|
|
# -------------------------------
|
|
flight_data <- bind_rows(all_flights)
|
|
|
|
# -------------------------------
|
|
# Step 4: Inspect data
|
|
# -------------------------------
|
|
head(flight_data)
|
|
colnames(flight_data)
|
|
str(flight_data)
|
|
|
|
if(interactive()) {
|
|
View(flight_data)
|
|
}
|
|
|
|
# -------------------------------
|
|
# Step 5: Filter flight BEL40V
|
|
# -------------------------------
|
|
flight_number <- "BEL40V"
|
|
|
|
if("callsign" %in% colnames(flight_data)) {
|
|
flight_BEL40V <- flight_data %>%
|
|
filter(callsign == flight_number)
|
|
|
|
# save
|
|
write_csv(flight_BEL40V, "flight_BEL40V.csv")
|
|
|
|
# view first row
|
|
head(flight_BEL40V)
|
|
} else {
|
|
warning("Column 'callsign' not found in flight_data!")
|
|
}
|
|
```
|
|
|