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 лучше?
Большое спасибо
Это намного проще для нас, если вы даете нам образец набора данных, который мы можем работать. Можете ли вы вставить то, что выглядит 'dput (datainSUB)'? – nograpes
Почему вы обертываете его в pdf-вызов? Сначала запустите его на экране. – Spacedman