График интенсивности непрерывного с geom_tile в ggplot

Question

Я пытаюсь построить непрерывную переменную в пространстве. Я видел этот пример, который получает тот же результат, который мне нужен:

library("MASS")
library("ggplot2")
library(reshape2) 

DB<-melt(volcano)
ggplot(DB, aes(x=Var1, y=Var2, fill=value)) +geom_tile()

Итак, я сделал с моими данными:

library(repmis)
url<-"https://www.dropbox.com/s/4m5qk32wjgrjq40/dato.RDATA"
source_data(url)

library(ggplot2)
ggplot(dato,aes(y=variable,x=y,fill=value))+geom_tile()

Это прекрасно. Но мои "х" и "у" - это километры (восток и север) от точки в пространстве. Я преобразовал их в широту и долготу. Но теперь мой сюжет не работает!

ggplot(dato,aes(y=lat,x=long,fill=value))+geom_tile()

Я не понимаю почему. В любом случае, при построении моих данных в виде точек результат очень похож:

ggplot(dato,aes(y=lat,x=long,fill=value))+geom_point()
ggplot(dato,aes(y=variable,x=y,fill=value))+geom_point()

Answer 1

Вы можете немного обмануть и использовать geom_point с квадратной формой:

#devtools::install_github("sjmgarnier/viridis")
library(viridis)
library(ggplot2)
library(ggthemes)
library(scales)
library(grid)

gg <- ggplot(dato)
gg <- gg + geom_point(aes(x=long, y=lat, color=value), shape=15, size=5)
gg <- gg + coord_equal()
gg <- gg + scale_color_viridis(na.value="#FFFFFF00")
gg <- gg + theme_map()
gg <- gg + theme(legend.position="right")
gg

Я не проецировал пары lat/long и просто использовал coord_equal, Вы должны использовать правильную проекцию для отображаемого региона.

И теперь у вас есть любопытный вопрос о том, что это за горячие точки вокруг Милана:-)

gmap <- get_map(location=c(9.051062, 45.38804, 9.277473, 45.53438),
                 source="stamen", maptype="toner", crop=TRUE)
gg <- ggmap(gmap)
gg <- gg + geom_point(data=dato, aes(x=long, y=lat, color=value), shape=15, size=5, alpha=0.25)
gg <- gg + coord_map()
gg <- gg + scale_color_viridis(na.value="#FFFFFF00")
gg <- gg + theme_map()
gg <- gg + theme(legend.position="right")
gg