ここで、ランドスケープ GIS レイヤーは、特に種やサンプルの分布を表示するのに非常に魅力的であることがわかりました。Rまたはその他のリソースで到達できるかどうかを知りたいですか?
GIS レイヤーは、この記事の図 1 で使用されました ( http://onlinelibrary.wiley.com/doi/10.1111/j.1469-8137.2010.03479.x/full )。
この図 1 の画像は次のとおりです。
R-sig-geo でこの質問を提出しようとしたことがあります。しかし、私は失敗しました。ここでいくつかのヘルプ/指示が得られることを期待しています。
方向性についてどうもありがとう。
幸運をお祈りしています、
このファイルをダウンロードして R で読み込み、正しい地理座標を持つように構成して、オーバープロットが簡単に機能するようにし、正しい配色で画像を表示することなどは非常に可能です。ただし、必要なすべてのデータの取得を自動化するのはそれほど簡単ではありません。
各ピクセルのRGB値の正しいセットをプロットできるように、GIFファイルのカラーテーブルが必要です（情報はファイルにありますが、Rで直接取得できるかどうかはわかりません。確認します-確かに GDAL を使用することもできますが、自動化された方法でこれらの値を抽出することは、利用可能なさまざまなツールに依存します)。
更新:ラスター パッケージが色情報を正しく取得し、それをプロットすることが判明しました。以下を参照してください。
また、地理空間情報、つまり参照ピクセルの座標 (ピクセルの左上隅など) とスケール (ピクセルの地理的な幅と高さ) も必要ですが、この情報はファイルに保存されません。また、ファイルの座標系はファイルに含まれておらず、画像データで明示的に提供されていない可能性が非常に高いです。
色と座標系がファイルに保存されている場合は、すべて簡単で、次のようなもので十分です。
（これは一度だけうまくいきましたが、その後のリクエストはサーバーによってブロックされていると思うので、ファイルを一度だけダウンロードしてみてください）。
u <- "http://onlinelibrary.wiley.com/store/10.1111/j.1469-8137.2010.03479.x/asset/image_n/NPH_3479_f1.gif?v=1&t=gskxvi17&s=0f13fa9dae78bd6837aeee594065c6ca112864d2"
imfile <- paste(tempfile(), ".gif", sep = "")
download.file(u, imfile, mode = "wb")
library(raster) ## rgdal also required for this file format
library(rgdal)
im <- raster(imfile)
plot(im)
これは問題ないように見えますが、「現実世界」の座標系がないことがわかります。これは、ピクセル 1 から X 次元の数値までの軸です (Y についても同じです)。
axis(1, pos = 2)
そのため、画像の適切な参照座標を見つけるために手動で作業する必要があります。ここでの推測は問題なく機能しますが、それでも推測にすぎず、一見単純に見えるものに対して多くの苦痛を生み出すことになる可能性があります。
ポイントをインタラクティブにプロットするだけで十分な場合は、および、および関連するプロット関数とlocator組み合わせて使用できます。pointslinestext
Feng、Googleドキュメントを正しく読んだら、追加のパラメータスタイルと要素を使用してラベルと表示された機能を変更できます。これらのカスタムパラメータをRgoogleMapsパッケージに含めませんでしたが、path引数を介して任意の追加パラメータを簡単に渡すことができます。GetMapのヘルプファイルを注意深く読むと、次の例に気付くでしょう。
＃次の例は、地方道路が明るい緑に変更され、住宅地が黒に変更されたブルックリンの地図を示しています。##実行されません：GetMap（center ='Brooklyn'、zoom = 12、maptype = "roadmap" 、path = "＆style = feature：road.local | element：geometry | hue：0x00ff00 | saturation：100＆style = feature：landscape | element：geometry | lightness：-100"、sensor ='false'、destfile = "MyTile4.png "、RETURNIMAGE = FALSE）;
お役に立てれば、
Markus Loecher
この画像のようなデータだけが必要な場合は、 と のツールを使用して、画像に直接アクセスするためのパッケージがありspますrgdal。gmapこの例は、dismoパッケージでの使用に近いです。
library(dismo)
e <- extent(-7, 5, 38, 44)
gm <- gmap(e, type = "terrain")
plot(gm)
「longlat」で範囲を指定している間、画像はネイティブ (Google) メルカトルに戻ることに注意してください。
print(gm)
?gmap独自のデータを変換してイメージのプロジェクションに一致させるためのその他のオプション、または で設定されたより広範な関数raster、rgdalおよびその他のオプションについては、を参照してくださいsp。推奨される可能性のある画像プロバイダーは他にもあり、提供されたパッケージの R スイートにはかなりの数のオプションがあります。