0

R を使用して、ウェイポイントとその補間 (N と呼ばれる値に基づく) を確認できる補間マップを取得しようとしています。背景には、道路と小道が見える白黒の地図が必要でした。複雑なことではないと確信していますが、ここでどこを間違えているのかわかりません。

私のデータは次のようになります。

  ALLdata
  ID  lat      long     N
1 213 52.36483 20.77724 5
2 214 52.36483 20.77724 5
3 215 52.36483 20.77724 5
4 216 52.20367 20.97026 6
5 217 52.20267 20.96939 7
6 218 52.20371 20.96483 0
7 219 52.20347 20.96595 1
8 220 52.20260 20.96884 2

これが私のコードです:

ALLdata <- read.csv("C:/Users/Me/Desktop/sample_animals.csv",header=TRUE,sep=";")
head(ALLdata)  
setwd("C:/Users/Me/Desktop/dataFieldseason2016/Graphics") #save graphs and maps

# grab a map. get_map creates a raster object
library(ggmap)
rwanda1 <- get_map(location = c(lon = 20.7, lat = 52.3),
               zoom = 9,
               maptype = "toner",
               source = "stamen")

g1 <- ggmap(rwanda1)
g1

# plot map and animals data
# use coord_map with default mercator projection
g1 + 
  geom_tile(data = ALLdata, aes(x = long, y = lat, z = N, fill = N), alpha = 0.8) +
  stat_contour(data = ALLdata, aes(x = long, y = lat, z = N)) +
  ggtitle("animals presence") +
  xlab("Longitude") +
  ylab("Latitude") +
  scale_fill_continuous(name = "animals (N)",
                    low = "white", high = "blue") +
  theme(plot.title = element_text(size = 25, face = "bold"),
    legend.title = element_text(size = 15),
    axis.text = element_text(size = 15),
    axis.title.x = element_text(size = 20, vjust = -0.5),
    axis.title.y = element_text(size = 20, vjust = 0.2),
        legend.text = element_text(size = 10)) +
  coord_map()

私はこれを得た: ここに画像の説明を入力

以下はうまくいかなかったものです。ウェイポイントをこの補間クラウド上のポイントとして表示したいと思います。最高はより暗く、最低はより明るいまたは白です。 ここに画像の説明を入力

ポイントまたは補間を確認できません。

4

0 に答える 0