[This article was first published on The Devil is in the Data – The Lucid Manager, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

I am soon embarking on a trip around the world to speak at the World Water Congress in Tokyo, visit family and friends in the Netherlands and some tourism in San Francisco. One of the workshops I am involved with is about Fake News in the water industry.

Fake news and conspiracies are threatening western civilisation, and one of the craziest trends is the flat earth conspiracy. Some of you might ask whether my trip actually around the world or am I travelling across a flat disk?

This article discusses flat earth mathematics, and how to convert and visualise map projections in R. This article uses the code I published earlier in articles about creating flight maps and Pacific island hopping. You can view the code below or download the full version from GitHub.

## The Flat Earth

YouTube contains thousands of videos from people that claim the earth is flat and that science as we know it is a “globalist conspiracy”. While their claims are hilarious, people that believe our planet is flat are often passionate and even conduct simple experiments and engage in naive flat earth mathematics.

Adherents to the idea that the world is flat often propose Gleason’s Map (shown above) as their version of the correct representation of our planet. This map is not a map of the flat earth but a polar Azimuthal equidistant projection of the globe. Even though the map itself states that it is a projection of the world, flat earth believers nevertheless see it as a literal representation of the earth.

Projecting the spherical earth on a flat surface is a complex task which will always require compromise as it is impossible to draw the surface of the globe on a piece of paper truthfully.  The video below provides a practical introduction to map projections.

We can recreate Gleason’s map with ggplot, which incorporates the mapproj package to show maps in various projections. The Azimuthal equidistant method is popularised in the flag of the United Nations. Antarctica is in this projection displayed as a ring around the world. Flat earth evangelists believe that the South Pole a ring of ice that prevents us from proceeding beyond the disc.

library(tidyverse)
library(tidyverse)
world <- map_data("world")
worldmap <- ggplot(world) +
geom_path(aes(x = long, y = lat, group = group)) +
theme(axis.text = element_blank(),
axis.ticks = element_blank()) +
labs(x = "", y = "")
worldmap +
coord_map("azequidistant", orientation = c(90, 0, 270))


The coord_map function allows you to change projections. The orientation argument defines the centre point of the projection. The most common location is the North Pole (90, 0). The third value gives the clockwise rotation in degrees. Gleason’s map uses a rotation of 270 degrees.

## My Round-the-World Itinerary

My trip will take me from Australia to Japan, from where I go to the Netherlands. The last two legs will bring me back to Australia via San Francisco.

The itinerary is stored in a data frame, and the ggmap package geocodes the longitude and latitude of each of the locations on my trip. As the earth is a sphere, an intermediate point needs to be added for trips that pass the dateline, as I explained in my article about flight maps.

The itinerary is visualised using the same method as above but centring on Antarctica. The geosphere package helps to estimate the total travel distance, which is approximately 38,995 km, slightly less than a trip around the equator. This distance is the great circle distance, which is the shortest distance between two points on a sphere, adjusted for a spheroid (a flattened sphere).

The flight paths on this map are curved because of the inevitable distortions when projecting a sphere on a flat surface.

## Define itinerary
library(ggmap)
airports <- tibble(city = c("Melbourne", "Tokyo", "Amsterdam", "San Francisco", "Melbourne"))

## Find coordinates
## Where you receive an OVER_QUERY_LIMIT error, repeat the code
itinerary <- geocode(airports$city) %>% mutate(location = airports$city)

## Split travel past dateline
dl <- which(diff(itinerary$lon) > 180) dr <- ifelse(itinerary$lon[dl] < 0, -180, 180)
dateline <- tibble(lon = c(dr, -dr),
lat = rep(mean(itinerary$lat[dl:(dl + 1)]), 2), location = "dateline") itinerary <- rbind(itinerary[1:dl, ], dateline, itinerary[(dl + 1):nrow(itinerary), ]) itinerary ## Visualise worldmap + geom_point(data = itinerary, aes(lon, lat), colour = "red", size = 4) + geom_path(data = itinerary, aes(lon, lat), colour = "red", size = 1) + coord_map("azequidistant", orientation = c(-90, 0, 270)) ggsave("rtw.png", dpi = 150) ## Great Circle Distance library(geosphere) sapply(1:(nrow(itinerary) - 1), function(l) distVincentyEllipsoid(itinerary[l, 1:2], itinerary[(l + 1), 1:2]) / 1000) %>% sum()  ## Flat Earth Mathematics If the Gleason map were in reality a map of the flat earth, then the flight paths on the map would show as straight lines. I was unable to plot these lines on the map because ggplot corrects for the effects of projecting coordinates. The mapproj package contains the mapproject function to calculate the projected coordinates based on longitude and latitude. The output of this function is a grid with limits from $-\pi$ to $\pi$. On the polar Azimuthal equidistant projection, distances from any point to the centre of the map are not distorted. A line from the lon/lat 0,0 to the north pole has a projected distance of $\pi/2$, which in the spherical world is $\pi / 2 * 6378.137 = 10018.75$ km. In other words, we need to multiply the Euclidean distance with the radius of the Earth to estimate the length on a flat map. This last code snippet converts my itinerary to projected coordinates, calculates the Euclidean distance between the points and multiplies this with the Earth’s diameter. The outcome is a staggering 83,675 km, more than twice the length over a sphere. The main reason for this difference is the distortion of this projection. When taking Antartica as the centre of the projection, East-West lengths on the Northern hemisphere are massively distorted. When undertaking the same calculation with the North pole as the center, the total distance is 46,348 km, still a lot more than the actual length of the journey. The only leg with a similar distance to reality is from Brisbane to Tokyo because this route is almost a radial line without distortion. This article proves that the Gleason map is not a representation of a flat earth because aeroplanes would have to break the sound barrier to fly these distances in the time it takes to travel. Whichever way you project the globe on a flat map will lead to inevitable distortions. The Gleason map itself mentions that it is a projection. But these facts will not stop people believing in a flat earth, I am after all an engineer and thus part of the globalist science conspiracy. You can view the complete code on my GitHub repository. library(mapproj) coords <- mapproject(itinerary$lon, itinerary$lat, "azequidistant", orientation = c(-90, 0, 270)) coords <- tibble(x = coords$x, y = coords$y) sum(sqrt(diff(coords$x)^2 + diff(coords\$y)^2) * 6378.137)


The post Flat Earth Mathematics in the R Language appeared first on The Lucid Manager.