9

In the igraph package for R, I am struggling to plot a social network using latitude/longitude coordinates as the layout of the graph.

Imagine this simple example: a network with 4 nodes of which you know the geographical location and the connections:

df<-data.frame("from" = c("Bob", "Klaus", "Edith", "Liu"), "to"= c("Edith", "Edith", "Bob", "Klaus"))

Here you have the meta-data for the nodes, so that Bob lives in NY, Klaus in Berlin, Edith in Paris and Liu in Bejing:

meta <- data.frame("name"=c("Bob", "Klaus", "Edith", "Liu"), "lon"=c(-74.00714, 13.37699, 2.34120, 116.40708),  "lat"=c(40.71455, 52.51607, 48.85693, 39.90469))

We make g the igraph object...

g <- graph.data.frame(df, directed=T, vertices=meta)

...and we define our layout as the longitude/latitude coordinates

lo <- layout.norm(as.matrix(meta[,2:3]))
plot.igraph(g, layout=lo)

If you run this example with these (real) geo-coordinates, you'll see that it is "relatively" accurate in the sense that the locations are right relative to each other. However, if I plot a lot of coordinates like this, the worlds cartesian map looks "stretched out".

Is there a way I can really plot my nodes on a world map so that the coordinates are 100% right and I get to see the connections between my nodes? I really want to keep using the igraph package as it offers a lot of functionalities I might need later on when I want to analyze the links between nodes.

josliber
  • 43,891
  • 12
  • 98
  • 133
willemr
  • 93
  • 1
  • 4
  • have you considered making a data frame from the igraph nodes and then plotting them over a map with the ggplot subsystem? – hrbrmstr Jun 04 '15 at 16:00
  • To @hrbrmstr's suggestion, this creates the world map, but to overlay or integrate edges and igraph capabilities is beyond me. require(maps) countries <- map_data("world") countries <- subset(countries, region != "Antarctica") ggplot(meta, aes(x=lon, y=lat)) + geom_point(shape = 22, fill = "red", size = 6) + geom_polygon(data = countries, aes(x=long, y=lat, group=group), col= "grey30", fill = "transparent", lwd=0) – lawyeR Jun 04 '15 at 16:46
  • 2
    You need to call here `igraph::plot()` with `rescale=FALSE` parameter, otherwise it rescales the coordinates along both dimensions into a `[-1; 1]` interval. – deeenes Jun 04 '15 at 19:40
  • Also, you would probably map the latitudes and longitudes into the XY plane using some kind of a map projection; the Mercator projection (see https://www.wikiwand.com/en/Mercator_projection ) is usually okay for most practical purposes. – Tamás Jun 04 '15 at 19:57
  • How about plotting great circles? http://dbsgeo.com/arc.js/ Needs some Javascript chops until someone makes a plugin for the leaflet R package... – Spacedman Jun 04 '15 at 21:56
  • @hrbrmstr and @lawyeR: I have indeed considered that, but the problem is connecting my nodes on the map without connecting them all by using the `+ geom_line()` command. @deeenes: if I include the `rescale=F` parameter, it gives me a blank graph... – willemr Jun 05 '15 at 06:32

1 Answers1

7

One element of the solution is no doubt the rescale = FALSE parameter to igraph::plot() as I suggested in comment. OP asked why ey gets an empty plot with this? It's because the plotting area is still limited to the [-1; 1] interval along both x and y axes. This is default of igraph::plot(). So we need to give xlim = c(-180, 180) and ylim = c(-90, 90) parameters. This already gives a correct positioning. However, if our aim is to produce a figure with the map of the world, maybe it is the best to write the igraph plot onto a cairo SVG device. Then we will be able to place the map behind the graph in any SVG editor (e.g. Inkscape is a great solution), and we are still free to scale and edit the graph and the labels. For doing this, some other igraph.plotting parameters necessary to be set, but this is already about proportions and aesthetics. Here is the code I used to produce the SVG output:

#!/usr/bin/Rscript

require(igraph)
require(Cairo)

df <- data.frame("from" = c("Bob", "Klaus", "Edith", "Liu"), 
    "to" = c("Edith", "Edith", "Bob", "Klaus"))

meta <- data.frame("name" = c("Bob", "Klaus", "Edith", "Liu"), 
    "lon" = c(-74.00714, 13.37699, 2.34120, 116.40708), 
    "lat" = c(40.71455, 52.51607, 48.85693, 39.90469))

g <- graph.data.frame(df, directed = TRUE, vertices = meta)

lo <- layout.norm(as.matrix(meta[,2:3]))

dpi = 1.0
Cairo(file = 'map-graph.svg', type = "svg", 
    units = "in", 
    width = 4 / dpi, 
    height = 2 / dpi, 
    dpi = dpi)

plot.igraph(g, 
    layout = lo, 
    xlim = c(-180, 180), 
    ylim = c(-90, 90), 
    rescale = FALSE, 
    edge.curved = TRUE, 
    edge.arrow.size = 10 / dpi, 
    edge.arrow.width = 0.5 / dpi, 
    vertex.label.dist = 50 / dpi, 
    vertex.label.degree = 90 / dpi, 
    vertex.size = 200 / dpi, 
    vertex.label.cex = 21 / dpi,
    vertex.frame.color = NA, 
    vertex.label.color = '#FFFF00', 
    edge.color = '#FFFFFF',
    vertex.label.family = 'sans-serif',
    edge.width = 16 / dpi)

dev.off()

When the SVG produced by igraph looks fine, we can open it in Inkscape. Then import (Ctrl+i) the map in case it is pixmap; or open if it is vector graphics (e.g. PDF, SVG). Manually scale and position the map to set the same scale as the graph in the SVG (i.e. until the points get to their right place) – for proportional scaling, hold the Ctrl in Inkscape. Here is the result of this method:

enter image description here

(The map image is made available for non-commercial public use by Wikimedia Commons).

I think igraph is capable to produce figures like this, but this is not the primary aim of this software, so it has its limitations. At some point you might consider to use some Geographic Information System (GIS) software which is designed exactly to do things like this. I have no experience with those, but qgis probably worths to take a look at.

deeenes
  • 4,148
  • 5
  • 43
  • 59
  • 1
    This is perfect. Great solution. The great circles suggestion was really useful as well (http://flowingdata.com/2011/05/11/how-to-map-connections-with-great-circles/) to visually do what I needed. – willemr Jun 08 '15 at 09:10