r - ggmap を使用した切断密度ポリゴン

Question

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)

ただし、輪郭はきれいに見え、一部のポリゴンは問題ありませんが、境界と交差するポリゴンは、正しい「滑らかな」部分とポリゴンを閉じる直線部分の 2 つのコンポーネントに分割されます。これらの 2 つの部分は、マップの境界で交わります。

私のデータはマップの境界を越えて広がっているため、KDE ​​が境界までの密度の意味のある見積もりを導き出すための十分な情報があります。

問題が何であるかを知っている人はいますか？そして、もっと重要なのは、どうすれば修正できるのでしょうか?

ありがとう、アンドリュー。