I don't understand what you mean by "too small", and you haven't
provided a reproducible example.
The map will be as large as the default graphics window on your screen
will allow (leaving room for axes and labels as well).
If you are not happy with this, try something like:
x11(width=30, height=18)
map_world()
Otherwise, you will need to describe exactly what you mean by "really
too small to show the data effectively" and provide a reproducible example.
Ray Brownrigg
On 10/02/2013 5:06 a.m., Frank Rodgers wrote:
I am fairly new to R and am plotting flight data on a map. Everything is
working well except the size of the map is really too small to show the data
effectively and I can't seem to figure out how to make the output map
larger. Do I need to change the device characteristics or is it a map.???
call. Here is the code:
library(maps)
library(geosphere)
airports <- read.csv("airports.csv", header=TRUE)
flights <- read.csv("FSRflights.csv", header=TRUE, as.is=TRUE)
checkDateLine <- function(l){
n<-0
k<-length(l)
k<-k-1
for (j in 1:k){
n[j] <- l[j+1] - l[j]
}
n <- abs(n)
m<-max(n, rm.na=TRUE)
ifelse(m > 30, TRUE, FALSE)
}
clean.Inter <- function(p1, p2, n, addStartEnd){
inter <- gcIntermediate(p1, p2, n=n, addStartEnd=addStartEnd)
if (checkDateLine(inter[,1])){
m1 <- midPoint(p1, p2)
m1[,1] <- (m1[,1]+180)%%360 - 180
a1 <- antipode(m1)
l1 <- gcIntermediate(p1, a1, n=n, addStartEnd=addStartEnd)
l2 <- gcIntermediate(a1, p2, n=n, addStartEnd=addStartEnd)
l3 <- rbind(l1, l2)
l3
}
else{
inter
}
}
add_lines <- function(){
pal <- colorRampPalette(c("#00FF00", "#FF0000"))
colors <- pal(100)
fsub <- flights[flights$airline == "FSR",]
fsub <- flights[order(flights$cnt),]
maxcnt <- max(fsub$cnt)
for (j in 1:length(fsub$airline)) {
air1 <- airports[airports$iata == fsub[j,]$airport1,]
air2 <- airports[airports$iata == fsub[j,]$airport2,]
p1 <- c(air1[1,]$long, air1[1,]$lat)
p2 <- c(air2[1,]$long, air2[1,]$lat)
inter <- clean.Inter(p1,p2,n=100, addStartEnd=TRUE)
colindex <- round( (fsub[j,]$cnt / maxcnt) * length(colors))
lines(inter, col=colors[colindex], lwd=0.1)
}
}
map_usa <- function(){
xlim <- c(-171.738281, -56.601563)
ylim <- c(12.039321, 71.856229)
map("world", col="#191919", fill=TRUE, bg="#736F6E", lwd=0.05, xlim=xlim,
ylim=ylim)
add_lines()
}
map_world <- function(){
map("world", col="#191919", fill=TRUE, bg="#736F6E", lwd=.1)
add_lines()
}
map_world()
--
View this message in context:
http://r.789695.n4.nabble.com/R-maps-Help-tp4658030.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
