データポイントをプロットし補間する

Question

データをマップにプロットして補間したいが、rでどのように行うのか分からない。これは私がこれまで行ってきたことです。

library(ggmap) 

gc <- geocode("Rakiraki,Fiji") 
Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=Rakiraki,Fiji&sensor=false 
map <- get_map(gc) 
Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=-17.399264,178.070532&zoom=10&size=640x640&scale=2&maptype=terrain&language=en-EN&sensor=false 
(bb <- attr(map, "bb")) 
    ll.lat ll.lon ur.lat ur.lon 
1 -17.81878 177.6318 -16.9801 178.5107 
(bbox <- bb2bbox(bb)) 
    left bottom  right  top 
177.63177 -17.81878 178.51067 -16.98010 
ggmap(map) 
+ geom_point(
+ + aes(x = lon, y = lat), 
+ + data = gc, colour = "red", size = 3 
Error: unexpected '=' in: 
"+ aes(x = lon, y = lat), 
+ data =" 
+) 
Error: unexpected ')' in "+)" 
data1 <-read.csv(file.choose(),header=T) 
data1 
    Longitude Latitude Wind.Speed 
1  177.8 -17.6  4.0 
2  177.5 -17.5  3.5 
3  178.0 -17.4  4.5 
4  178.1 -17.6  4.0 
5  178.2 -17.3  6.0 
6  178.3 -17.7  6.5 
7  178.3 -17.5  5.0 
8  178.4 -17.6  5.5 

どのように私はこれらのデータをプロットすることができますし、どのようにマップ上でそれを補間することができます。

私は、最終的なマップはこのような何かを見てみたいです。添付ファイルを見つけてくださいenter image description here

Answer 1

あなたの質問にはかなりの問題があります。かっこはgeom_pointの最後に閉じておらず、存在してはならない+のシンボルが存在します。

data1のデータを、つまり異なる場所の風速に追加するとします。このような何か作業をする必要があります：

ggmap(map) + 
    geom_point(aes(x = Longitude, 
       y = Latitude, 
       size = Wind.Speed), 
      data = data1, 
      colour = "red") 

それはあなたが「補間」おそらくstat_density2dのようなものは、あなたが考えていたものであることにより、何を意味するのか、正確には明らかではないのか？上記のコードにこれを追加します。

+ stat_density2d(data = data1, 
       aes(x = Longitude, y = Latitude))