2014-02-10 2 views
0

У меня есть поле, которое географической привязкой я совмещен с функцией изображенияПостроение графиков Сетчатая поле

df <- datainSUB 

yr mo dy hr lon lat cell sst avg moavg 
1900 6 5 17 -73.5 -60.5 83 2.4 2.15 3.15 
1900 6 7 17 -74.5 -60.5 83 3.9 2.15 3.15 
1900 8 17 17 -70.5 -60.5 83 -0.9 2.15 0.60 
1900 8 18 17 -73.5 -60.5 83 2.1 2.15 0.60 
1900 9 20 17 -71.5 -60.5 83 0.2 2.15 2.20 
1900 9 21 17 -74.5 -61.5 83 1.6 2.15 2.20 

gridplot <- function(df){ 
    pdf(paste(df$mo,".pdf")) 
    # Compute the ordered x- and y-values 
    LON <- seq(-180, 180, by = space) 
    LAT <- seq(-90, 90, by = space) 

    # Build the matrix to be plotted 
    moavg <- matrix(NA, nrow=length(LON), ncol=length(LAT)) 
    moavg[cbind(match(round(df$lon, -1), LON), match(round(df$lat, -1), LAT))] <- df$moavg 

    # Plot the image 
    image(LON, LAT, moavg) 
    map(add=T,col="saddlebrown",interior = FALSE, database="world") 
    dev.off() 
} 

Я хочу, чтобы добавить цвет легенды к сюжету, но я не знаю, как сделать это. Может быть, ggplot лучше?

Большое спасибо

+0

Это намного проще для нас, если вы даете нам образец набора данных, который мы можем работать. Можете ли вы вставить то, что выглядит 'dput (datainSUB)'? – nograpes

+1

Почему вы обертываете его в pdf-вызов? Сначала запустите его на экране. – Spacedman

ответ

0

Добавьте следующую строку после построения ваши данные:

legend(x="topright", "your legend goes here", fill="saddlebrown") 
Смежные вопросы