spatstat パッケージと ppp 関数を使用して、xy 位置に基づいて木の空間マップを作成しています。私はこれを行うことができました.各ポイントのサイズは木の直径に基づいています. 私ができなかったのは、ここでいくつかの助けを見つけたいと思っていることですが、これらの異なる円を表示し、プロット内に6つの樹種を含めることができる凡例を作成することです. 私が使用しているコードは次のとおりです。
df <- subset(plots, plots$spp == "DF") # Douglas-fir
dfx <- as.numeric(as.character(df$x))
dfy <- as.numeric(as.character(df$y))
dfd <- as.numeric(as.character(df$d))
dfp <- ppp(dfx, dfy, window = owin(c(0, 100), c(0, 100)),
unitname=c("metres","metres"), marks = dfd)
par(mar = c(2, 2, 2, 2))
plot(dfp, main = "", cex = 0.8, markscale = 0.04,
bg = rgb(0.1,0.9,0.3,0.5), fg = "black")
プロット内のすべての種に同様の構造があります (ニシツガとコードは同じです...)
ありがとう! </p>