4

spplot() に領域の名前ラベルを追加したいと思います。

例:

load(url('http://gadm.org/data/rda/FRA_adm0.RData'))
FR <- gadm
FR <- spChFIDs(FR, paste("FR", rownames(FR), sep = "_"))
load(url('http://gadm.org/data/rda/CHE_adm0.RData'))
SW <- gadm
SW <- spChFIDs(SW, paste("SW", rownames(SW), sep = "_"))
load(url('http://gadm.org/data/rda/DEU_adm0.RData'))
GE <- gadm
GE <- spChFIDs(GE, paste("GE", rownames(GE), sep = "_"))

df <- rbind(FR, SW, GE)

## working
plot(df)
text(getSpPPolygonsLabptSlots(df), labels = c("FR", "SW", "GE"))


## not working
spplot(df[1-2,])
text((getSpPPolygonsLabptSlots(df), labels = c("FR", "SW"))

2枚目はラティスのせいで動かないかも!? ただし、spplot 機能が必要です。プロットのラベルを取得するにはどうすればよいですか?

4

2 に答える 2

6

テキストを追加する標準的な方法は ltext、格子関数を使用することですが、そこに指定された座標は常に絶対です。本質的に、テキストを追加した後に図を実際に再スケーリングすることはできません。例:

data(meuse.grid)
gridded(meuse.grid)=~x+y
meuse.grid$g = factor(sample(letters[1:5], 3103, replace=TRUE),levels=letters[1:10])
meuse.grid$f = factor(sample(letters[6:10], 3103, replace=TRUE),levels=letters[1:10])

spplot(meuse.grid, c("f","g"))
ltext(100,200,"Horror")

これらの図を生成します (スケーリング前とスケーリング後)

ホラー1 ホラー2

各パネル内の座標を使用して、カスタム パネル関数を使用できます。

myPanel <- function(x,y,xx,yy,labels,...){
    panel.xyplot(x,y,...)
    ltext(xx,yy,labels)
}

xyplot(1:10 ~ 1:10,data=quakes,panel=myPanel,
        xx=(1:5),yy=(1:5)+0.5,labels=letters[1:5])

(自分で実行して、どのように見えるかを確認してください)

このトリックは、splot 関数内でも使用できますが、使用するプロット関数を実際に確認する必要があります。splot のヘルプ ファイルには、可能なオプション (polygonsplot、gridplot、pointsplot) が記載されているため、それらのいずれかが意図したとおりに動作しているかどうかを確認する必要があります。上記のグリッドプロットを続けると、次のようになります。

myPanel <- function(x,y,z,subscripts,xx,yy,labels,...){
    panel.gridplot(x,y,z,subscripts,...)
    ltext(xx,yy,labels)
}
# I just chose some coordinates
spplot(meuse.grid, c("f","g"),panel=myPanel,xx=180000,yy=331000,label="Hooray")

これにより、再スケーリング可能な結果が得られ、各パネルにテキストが追加されます。

ここに画像の説明を入力 万歳

于 2011-02-28T17:02:43.653 に答える
5

ありがとう、ギャビン・シンプソン!

私はついに方法を見つけました。

将来的に他の人に役立つことを願って、私は自分のソリューションを投稿します:

sp.label <- function(x, label) {
    list("sp.text", coordinates(x), label)
}

ISO.sp.label <- function(x) {
    sp.label(x, row.names(x["ISO"]))
}

make.ISO.sp.label <- function(x) {
    do.call("list", ISO.sp.label(x))
}

spplot(df['ISO'], sp.layout = make.ISO.sp.label(df))
于 2011-02-28T16:53:39.273 に答える