2017-05-08 7 views
0

を使用しながら、以下のようにIは、2つのデータフレームを持っている:私はドロップオフデータフレームで繰り返される撮像データフレームにピックアップ配位(経度及び緯度)を見つけることを試みる可視化問題リーフレット

PickUP <- data.frame(pickuplong = c(-73.93909 ,-73.94189 ,-73.93754,-73.91638,-73.92792 ,-73.88634), pickuplat =c(40.84408,40.83841,40.85311,40.84966,40.86284,40.85628)) 

Dropoff <- data.frame(pickuplong = c(-73.93351 ,-73.93909 ,-73.93909 ,-73.80747,-73.95722,-73.91880), pickuplat =c(40.76621,40.84408,40.85311,40.69951,40.68877,40.75917), Droplong =c(-73.91300,-73.96259 ,-73.94870,-73.93860,-73.93633, -73.90690), Droplat =c(40.77777,40.77488 ,40.78493,40.84463,40.75977,40.77013)) 

。私は以下のコードを持っているが、私はこの上でエラーを得た:

library(sp) 
library(rgdal) 
library(leaflet) 
library(mapview) 
library(dplyr) 
a <- semi_join(Dropoff , PickUP , by = c("pickuplong","pickuplat")) 
a$ID <- 1:nrow(a) 
Dropoff_p <- a[, c("ID", "Pickup_longitude", "Pickup_latitude")] 
Dropoff_d <- a[, c("ID", "Dropoff_longitude", "Dropoff_latitude")] 
coordinates(Dropoff_p) <- ~Pickup_longitude + Pickup_latitude 
coordinates(Dropoff_d) <- ~Dropoff_longitude + Dropoff_latitude 
proj4string(Dropoff_p) <- CRS("+init=epsg:4326") 
proj4string(Dropoff_d) <- CRS("+init=epsg:4326") 
map_p <- mapview(Dropoff_p, color = "red") 
map_d <- mapview(Dropoff_d, color = "blue") 
map_p + map_d 

私のエラーは次のとおりです。

エラー$ < -.data.frame(TMP、 "ID"、値=にC(1L、0L)): 交換は2行があり、データがラップアップ時に0エラーがあります。 接続

答えて

0

データフレームをサブセット化すると、あなたは同じ列名を使用する必要が開くことができません。 Dropoff_pDropoff_dcoordinates(Dropoff_p)、およびproj4string(Dropoff_d)の列名を変更して、スクリプトが機能します。

さらに、mapviewパッケージには新しいアップデートがあります。必要に応じてmapviewをバージョン2.0.1に更新することができます。 col.regions = "red"col.regions = "blue"を追加することもできます。これは、新しいバージョンのように思われるので、color引数はポイントのアウトラインのみを変更します。塗りの色を変更するには、col.regionsを使用します。

library(sp) 
library(rgdal) 
library(leaflet) 
library(mapview) 
library(dplyr) 

a <- semi_join(Dropoff , PickUP , by = c("pickuplong","pickuplat")) 
a$ID <- 1:nrow(a) 
Dropoff_p <- a[, c("ID", "pickuplong", "pickuplat")] 
Dropoff_d <- a[, c("ID", "Droplong", "Droplat")] 
coordinates(Dropoff_p) <- ~pickuplong + pickuplat 
coordinates(Dropoff_d) <- ~Droplong + Droplat 
proj4string(Dropoff_p) <- CRS("+init=epsg:4326") 
proj4string(Dropoff_d) <- CRS("+init=epsg:4326") 
map_p <- mapview(Dropoff_p, color = "red", col.regions = "red") 
map_d <- mapview(Dropoff_d, color = "blue", col.regions = "blue") 
map_p + map_d 
+0

私は、$ ID <2行目にエラーが発生しました - 1:あなたは、私は、長方形のような他の形状にマップ上の円のアイコンを変更する方法nrow(A) – Behzad

+0

を知っていますか。ありがとう – Behzad

+0

@Behzadシンボルを円から長方形に変更する簡単な方法があるかどうかわかりません。詳細については、 'mapview' https://environmentalinformatics-marburg.github.io/mapview/introduction.htmlまたは'リーフレット 'のドキュメントを参照してください。ところで、私のソリューションが少なくともあなたの元の質問を解決できるかどうか。それらを回答として受け入れてください。 – www