提交: bd1f9f14a411b9cdebee9aa7a71c01798e046c11 - myprojects-hg-reborn (hg)

Commit MetaInfo

• 修订版: bd1f9f14a411b9cdebee9aa7a71c01798e046c11 (tree)
• 时间: 2014-02-27 17:29:53
• 作者: Lorenzo Isella <lorenzo.isella@gmai...>
• Commiter: Lorenzo Isella

Log Message

I added a code to superimpose a network on a geographhical map.

更改概述

• add: R-codes/plot_geo_network.R

差异

diff -r 96a16c42c924 -r bd1f9f14a411 R-codes/plot_geo_network.R
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/R-codes/plot_geo_network.R	Thu Feb 27 09:29:53 2014 +0100
@@ -0,0 +1,53 @@
+rm(list=ls())
+ 
+library(igraph)
+library(maps)
+ 
+ 
+t <- map("italy", plot=FALSE)
+ 
+ 
+pos1 <- c(16,40.3)
+pos2 <- c(12,44)
+pos3 <- c(13,42)
+ 
+pdf("italy_raw.pdf")
+plot(t$x, t$y, "l")
+points(pos1[1],pos1[2], bg="red", col="red", pch=21)
+points(pos2[1],pos2[2], bg="blue", col="blue", pch=21)
+points(pos3[1],pos3[2], bg="brown", col="brown", pch=21)
+ 
+ 
+dev.off()
+ 
+ 
+ 
+el <- matrix( c("foo", "bar", "bar", "foobar"), nc=2, byrow=TRUE)
+g <- graph.edgelist(el)
+ 
+g <- as.undirected(g)
+ 
+l <- matrix(ncol=2, nrow=3)
+ 
+l[1,] <- pos1
+l[2,] <- pos2
+l[3,] <- pos3
+ 
+ 
+ 
+pdf("italy.pdf")
+plot(t$x, t$y, "l", xaxt='n', yaxt='n',ann=FALSE,  axes=F)
+## plot(g, layout=l,
+##       vertex.label.dist=0.5,vertex.label=NA, vertex.size=5)
+
+autocurve.edges (g, start = 0.5)
+
+plot(g, add=TRUE, rescale=FALSE, layout=l,
+     vertex.label.dist=0.5,vertex.label=NA, vertex.size=15
+     )
+
+ 
+dev.off()
+
+
+print("So far so good")