3

ライブラリ（RgoogleMaps）を使用して、地図上の測定位置（ポイント）をプロットします。ポイントごとに異なる機器があり、機器ごとに別々の色のポイントを取得することに成功しました。

theplot <- PlotOnStaticMap(lat=sitecoord$lat, lon=sitecoord$lon, 
                           cex=.7, pch=20, 
                           col=sitecoord$equipmentType, 
                           MyMap=Map, NEWMAP=FALSE)

結果のマッププロットに凡例を追加して、どの機器が青い点で表され、どの機器が赤で表されているかなどを確認するにはどうすればよいですか？

アップデート：

@Rguyの非常に優れた推奨事項を使用します。私はなんとか伝説を取り入れることができました。他の人のために、ここに私のテストコードがあります（いいえ、私はアイスランドで測定していません。例として使用しました）：

library(RgoogleMaps)
library(RColorBrewer)

Equipment <- c("AA","AA","BB","CC")
lat <- c(63.90,66.20,64.80,64.50)
lon <- c(-22.40,-14.20,-18.60,-15.00)
tblDataPoints <- data.frame(Equipment,lat,lon)

My.Pal <- brewer.pal(3, "Reds")
tblDataPoints$colorz <- My.Pal[tblDataPoints$Equipment]

plot.new()
bb <- qbbox(lat=range(tblDataPoints$lat), lon=range(tblDataPoints$lon))
m <- c(mean(tblDataPoints$lat), mean(tblDataPoints$lon))
zoom <- min(MaxZoom(latrange=bb$latR,lonrange=bb$lonR))
Map <- GetMap.bbox(bb$lonR, bb$latR, zoom=zoom, maptype="roadmap", NEWMAP=TRUE)
tmp <- PlotOnStaticMap(lat=lat, lon=lon, cex=.7, pch=20, col=tblDataPoints$colorz, MyMap=Map, NEWMAP=FALSE)

tblLgd <- unique(tblDataPoints[,c("Equipment","colorz")])
row.names(tblLgd) <- NULL

legend("topright", legend = tblLgd$Equipment, fill = tblLgd$colorz, bg = "white")
4

1 に答える 1

1

私は前にこれをしました。関数で発生している問題の再現可能な例を作成した場合は、legendそれについて話し合うことができます。それまでは、漠然とした説明があります。

1.RColorBrewerを使用して口蓋を作成します。例えば：

library(RColorBrewer)
My.pal <- brewer.pal(9, "reds")

2.何らかの方法で、各ポイントに色を割り当てます。私の場合、WTビンの列とベクトルがあったので、重みをビニングし、対応するエントリをmy.palそのポイントの色にすることで、ポイントごとの色を生成しました。この例ではbinz、口蓋の赤の色合いが9つしかないため、ベクトルのビンは9つ未満であることに注意してください。

colorz <-My.Pal [cut（datas $ WT、labels = FALSE）]

3.色の引数を渡して、マップ上にプロットします。

PlotOnStaticMap（MyMap、lat = datas $ LAT、lon = datas $ LON、col = colorz）

4.最後に、凡例を作成し、それをマップに追加します。

legend（ "bottomleft"、legend = legend.txt、fill = My.pal、title = "I AM"、bg = "white"）

あなたがそれをすべて理解することを願っています！

于 2011-10-17T19:34:53.950 に答える