Усеченные полигоны плотности с помощью ggmap

У меня возникли проблемы с построением карт плотности с помощью R/ggmap. Мои данные выглядят так:

> head(W)
        date    lat     lon    dist
1 2010-01-01 31.942 -86.659 292.415
2 2010-01-10 32.970 -84.174  89.121
3 2010-01-17 31.000 -85.363 319.552
4 2010-01-17 31.457 -83.951 258.501
5 2010-01-17 31.073 -81.987 373.915
6 2010-01-17 33.210 -83.149 129.927

И я рисую это, используя это:

ggmap(atlanta.map, extent = "panel") +
  geom_point(data = W, aes(x = lon, y = lat)) +
  geom_density2d(data = W, aes(x = lon, y = lat)) +
  stat_density2d(data = W, aes(x = lon, y = lat, fill = ..level..,
    alpha = ..level..), size = 0.01, bins = 8, geom = 'polygon') +
  theme(axis.title = element_blank()) +
  scale_fill_gradient(low = "yellow", high = "red") +
  scale_alpha(range = c(.5, .75), guide = FALSE)

Но, хотя контуры выглядят нормально и НЕКОТОРЫЕ полигоны в порядке, многоугольники, пересекающие границы, разбиваются на две составляющие: правильная «гладкая» часть и прямолинейная часть, замыкающая многоугольник. Эти две части встречаются на границе карты.

Мои данные выходят за границы карты, поэтому у KDE достаточно информации для получения значимых оценок плотности вплоть до границ.

Кто-нибудь знает, в чем может быть проблема? И, что более важно, как я могу это исправить?

Спасибо, Эндрю.

введите здесь описание изображения


r ggplot2 ggmap
person DataWookie    schedule 10.12.2013    source источник

Ответы (1)


arrow_upward
3
arrow_downward

С некоторыми дальнейшими поисками в Google и охотой на SO я собрал решение.

ggmap(map, extent = "normal", maprange=FALSE) %+% W + aes(x = lon, y = lat) +
    geom_density2d() +
    stat_density2d(aes(fill = ..level.., alpha = ..level..),
                   size = 0.01, bins = 16, geom = 'polygon') +
    scale_fill_gradient(low = "green", high = "red") +
    scale_alpha(range = c(.00, .25), guide = FALSE) +
    coord_map(projection="mercator", 
              xlim=c(attr(map, "bb")$ll.lon, attr(map, "bb")$ur.lon),
              ylim=c(attr(map, "bb")$ll.lat, attr(map, "bb")$ur.lat)) +
    theme(legend.position = "none", axis.title = element_blank())
person DataWookie    schedule 10.12.2013
comment
ну и что исправить? - person rrs; 05.05.2016
comment
Здорово! Мне просто это тоже нужно. Это менялось с панели на нормальную и добавлялось coord_map(....) +, чтобы спецификация ограничений непосредственно из экстентов формы выполняла свою работу гладко. - person Peter.k; 05.02.2017