{ "id": 5467, "title": "R macro to visualize test flights (and great circle paths)", "url": "/competitions/flight2-milestone/discussion/5467", "author_name": "", "post_date": "2013-08-20T15:47:43.123Z", "votes": 2, "comment_count": 1, "views": 1316, "content": "

Below R-macro draws planes corresponding to test flights and all airports. It also draw s great circle path from current plane location to target airport.
Remark: currently it draws just first 50 test flights. Change parameter flightsToDraw if/when needed:

\n

# R-macro to plot planes and (all) airports corresponding to
# test flights. Draws great circle path from current plane
# location to target airport.

# Required libraries
library(maps)
library(geosphere)

\n

# Load in data sets
setwd("d:\\\\kaggle\\\\ge")
airports <- read.csv("airports.csv", header=TRUE)
testFlights <- read.csv("testFlights.csv", header=TRUE)

\n

# Calculate bounding box for test flights' longitudes
# and latitudes
xlim <- c(min(testFlights$CurrentLongitude),max(testFlights$CurrentLongitude))
ylim <- c(min(testFlights$CurrentLatitude),max(testFlights$CurrentLatitude))

\n

# Draw map using the bounding box
map("world", col="#ffffff", fill=TRUE, bg="white", lwd=0.05, xlim=xlim, ylim=ylim)

\n

# Plot desired number of test flights
flightsToDraw <- 50 # max value is nrow(testFlights)
# but drawing becomes too messy

for(i in 1:flightsToDraw){
  targetLongitude <- airports$longitude_degrees[airports$airport_icao_code==testFlights$ArrivalAirport[i]]
  targetLatitude <- airports$latitude_degrees[airports$airport_icao_code==testFlights$ArrivalAirport[i]]

 inter <- gcIntermediate(c(testFlights$CurrentLongitude[i],testFlights$CurrentLatitude[i]),
  c(targetLongitude, targetLatitude), n=50,  addStartEnd=TRUE)
  lines(inter,lty=2, col="#CC0000")
}

\n

# Plot planes
for(i in 1:flightsToDraw){
  inter <- gcIntermediate(c(testFlights$CurrentLongitude[i],testFlights$CurrentLatitude[i]),
c(testFlights$CurrentLongitude[i],testFlights$CurrentLatitude[i]), n=1, addStartEnd=TRUE)
  lines(inter, col="#0000FF",lwd=5)
}

\n

# Plot target airports
for(i in 1:nrow(airports)){
  inter <- gcIntermediate(c(airports$longitude_degrees[i],airports$latitude_degrees[i]),
c(airports$longitude_degrees[i],airports$latitude_degrees[i]), n=1,addStartEnd=TRUE)
  lines(inter, col="#00FF00",lwd=5)
}

\n

# Plot legend
legend(xlim[1] + (xlim[2]-xlim[1])*0.75,ylim[1] + (ylim[2]-ylim[1])*0.9,
  legend=c("Plane","Target airport"),text.col=c("#0000FF", "#00FF00"),ncol=1)

", "messages": [ { "id": "29065", "postDate": "08/20/2013 15:47:43", "content": "

Below R-macro draws planes corresponding to test flights and all airports. It also draw s great circle path from current plane location to target airport.
Remark: currently it draws just first 50 test flights. Change parameter flightsToDraw if/when needed:

\n

# R-macro to plot planes and (all) airports corresponding to
# test flights. Draws great circle path from current plane
# location to target airport.

# Required libraries
library(maps)
library(geosphere)

\n

# Load in data sets
setwd("d:\\\\kaggle\\\\ge")
airports <- read.csv("airports.csv", header=TRUE)
testFlights <- read.csv("testFlights.csv", header=TRUE)

\n

# Calculate bounding box for test flights' longitudes
# and latitudes
xlim <- c(min(testFlights$CurrentLongitude),max(testFlights$CurrentLongitude))
ylim <- c(min(testFlights$CurrentLatitude),max(testFlights$CurrentLatitude))

\n

# Draw map using the bounding box
map("world", col="#ffffff", fill=TRUE, bg="white", lwd=0.05, xlim=xlim, ylim=ylim)

\n

# Plot desired number of test flights
flightsToDraw <- 50 # max value is nrow(testFlights)
# but drawing becomes too messy

for(i in 1:flightsToDraw){
  targetLongitude <- airports$longitude_degrees[airports$airport_icao_code==testFlights$ArrivalAirport[i]]
  targetLatitude <- airports$latitude_degrees[airports$airport_icao_code==testFlights$ArrivalAirport[i]]

 inter <- gcIntermediate(c(testFlights$CurrentLongitude[i],testFlights$CurrentLatitude[i]),
  c(targetLongitude, targetLatitude), n=50,  addStartEnd=TRUE)
  lines(inter,lty=2, col="#CC0000")
}

\n

# Plot planes
for(i in 1:flightsToDraw){
  inter <- gcIntermediate(c(testFlights$CurrentLongitude[i],testFlights$CurrentLatitude[i]),
c(testFlights$CurrentLongitude[i],testFlights$CurrentLatitude[i]), n=1, addStartEnd=TRUE)
  lines(inter, col="#0000FF",lwd=5)
}

\n

# Plot target airports
for(i in 1:nrow(airports)){
  inter <- gcIntermediate(c(airports$longitude_degrees[i],airports$latitude_degrees[i]),
c(airports$longitude_degrees[i],airports$latitude_degrees[i]), n=1,addStartEnd=TRUE)
  lines(inter, col="#00FF00",lwd=5)
}

\n

# Plot legend
legend(xlim[1] + (xlim[2]-xlim[1])*0.75,ylim[1] + (ylim[2]-ylim[1])*0.9,
  legend=c("Plane","Target airport"),text.col=c("#0000FF", "#00FF00"),ncol=1)

", "rawMarkdown": "", "votes": null }, { "id": "29068", "postDate": "08/20/2013 16:35:14", "content": "

I have not checked if licenses of packages maps or geosphere are according to the competition rules (5 allowed licenses) - my guess is that they are not. So prepare to program (mostly gcIntermediate function) by yourself in case if you use above functionality as part of path estimations in your competition submissions.

", "rawMarkdown": "", "votes": null } ], "comments": [ { "id": 29068, "author_name": "sfin2015", "author_url": "", "post_date": "08/20/2013 16:35:14", "content": "

I have not checked if licenses of packages maps or geosphere are according to the competition rules (5 allowed licenses) - my guess is that they are not. So prepare to program (mostly gcIntermediate function) by yourself in case if you use above functionality as part of path estimations in your competition submissions.

", "votes": null, "replies": [] } ], "raw_markdown_by_id": { "29065": "", "29068": "" }, "source": "meta" }