ggmapベースマップにラスタを追加：ggplot2のinset_raster（）にアルファ（透明）と塗りつぶしの色を設定します。

Question

ggplot2にGoogleマップのベースマップをオーバーレイするラスタをマップにプロットします。したがって、私はこのようなget_map()とinsert_raster()を使用：

library(ggplot2) 
library(ggmap) 

bm <- ggmap(get_map(location = "Bangkok", maptype = "hybrid")) 

bm + inset_raster(as.raster(r), xmin = r@extent[1], xmax = r@extent[2], 
        ymin = r@extent[3], ymax = r@extent[4]) 

はalphaを設定し、fill色を変更するには、いずれかの可能性はありますか？

結果は次のようになります。

Answer 1

：ところで周りのこの作品は、もともとこのリンクで育った

は、私たちことがわかっthis blog entryからの更なる情報について

下の元の記事を読みます直接に空間ポリゴンを使用できますggplot::geom_polygon()

あなただけ後述するように、何か

を視覚化する必要がある場合はSPEEDをプロット取り組む方法

、そのようなプロットは、ピクセルの数が多い非常に遅くなることがあります。したがって、ポリゴンに変換する前に、ピクセル数（ほとんどの場合、マップ内の情報量を減らすことはできません）を減らすことを検討してください。したがって、raster::aggregateを使用して、ピクセル数を合理的な量に減らすことができます。

この例では、ピクセル数が4（2 * 2、水平方向*垂直）の順にどのように減少するかを示しています。詳細は?raster::aggregateを参照してください。

r <- aggregate(r, fact = 2) 
# afterwards continue with rasterToPolygons(r)... 

オリジナルポスト：

しばらくして、私はこの問題を解決する方法を発見しました。 ラスタをポリゴンに変換する！この考えは基本的にMarc Needham's blog postの後に実装されました。

しかし、1つの欠点があります：ggplotは、必然的に直面する多数のポリゴンでは非常に遅くなります。ただし、png()（またはその他の）デバイスにプロットすることで、スピードアップできます。ここで

---

は、コード例です：

library(raster) 
library(ggplot2) 
library(ggmap) 

r <- raster(....) # any raster you want to plot 
rtp <- rasterToPolygons(r) 
rtp@data$id <- 1:nrow(rtp@data) # add id column for join 

rtpFort <- fortify(rtp, data = rtp@data) 
rtpFortMer <- merge(rtpFort, rtp@data, by.x = 'id', by.y = 'id') # join data 

bm <- ggmap(get_map(location = "Shanghai", maptype = "hybrid", zoom = 10)) 

bm + geom_polygon(data = rtpFortMer, 
        aes(x = long, y = lat, group = group, fill = layer), 
        alpha = 0.5, 
        size = 0) + ## size = 0 to remove the polygon outlines 
    scale_fill_gradientn(colours = topo.colors(255)) 

これは、このようなものになり：

Answer 2

はちょうどこのに自分自身を探していました。

library(ggplot2) 
library(ggmap) 

bm <- ggmap(get_map(location = "Bangkok", maptype = "hybrid")) 

bm <- bm + geom_raster(...) # insert your raster here 

bm <- bm + coord_cartesian() 

plot(bm) 

私は：

Error: geom_raster only works with Cartesian coordinates.

この問題の回避策は、次のように（）coord_cartesian使用することです：私は、ラスタとggmap出力をオーバーレイしようとしていた遭遇問題は、次のエラーでしたあなたのラスターrがどこから来ているのか分からない。これを行うには、ラスタrをデータフレームに変換し、geom_raster（）命令に従ってデータを追加するだけで、座標が緯度/経度（つまりマップと同じ）であることを確認してください。

質問に答えるには、geom_raster（）を使用してアルファと塗りつぶしを操作できます。

これが役に立ちます。 fortifyなしでさらに高速 https://groups.google.com/forum/embed/#!topic/ggplot2/nqzBX22MeAQ