リーフレットR derivePolygons欠けている緯度が足りない

Question

マップ上に病気のイベントデータのサイトをプロットしようとしています。

私はデータをインポートするためにこれを使用します。

ByTown<-readOGR(dsn="C:/temp/lyme/Towns", layer="Towns", encoding = "UTF-8", verbose= FALSE) 

は、クラスを確認してください。

class(ByTown) 
#getting this result 
[1] "SpatialPolygonsDataFrame" 
attr(,"package") 
[1] "sp" 

それから私は、文字データへのすべての要素を変換し、私はまだSpatialPolygonsDataFrameを持っていることを確認してください

次に、オリジナルと同じタイトルケースにマージするデータをフォーマットします。

townCount$City<-str_to_title(townCount$City) 

それから私は、空間ポリゴンデータフレームにカウントデータをgeo_join：

ByTown<-geo_join(ByTown, townCount,"MCD_NAME", "City") 

その後、私はパレットを設定し、マッピングを実行します。

pal = colorQuantile("PuOr",ByTown$count, n=5) 
map<-leaflet(ByTown) %>% 
    addProviderTiles("CartoDB.Positron")%>% 
    addPolygons(fillColor = ~pal(count), 
      color = "#000000", 
      stroke = TRUE, 
      weight = 1, 
      smoothFactor = 0.5, 
      options(viewer = NULL)) 
map 

をそして、私はこのエラーを取得します：

Error in derivePolygons(data, lng, lat, missing(lng), missing(lat), "addPolygons") : 
    addPolygons must be called with both lng and lat, or with neither. 

私は座標スロットがあり、そこにデータがあります...私はエラーによって困惑し、オンラインで有用な回答を見つけることができません。ここでは、座標スロットの最初の多角形のヘッドれる：

head(nByTown@polygons[[1]]@Polygons[[1]]@coords) 

      [,1]  [,2] 
[1,] 1036519 916318.7 
[2,] 1036039 916355.8 
[3,] 1031757 916299.7 
[4,] 1027474 916244.5 
[5,] 1026709 916198.1 
[6,] 1026826 916248.3 

いずれも、すべてのは、この問題を持って、根本原因を特定し、それを修正しますか？

Answer 1

addPolygons関数では、経度の緯度が&の列を定義する必要があります。または、入力したデータから緯度を計算します。

lat/lon列を指定していない場合、それらの列はどれかを調べようとします。あなたのケースでは、それを見つけることができませんでした。これは主に、あなたのデータが緯度/経度形式でないことが原因です。

したがって、あなたはしばらくの間、これらの形状ファイルを使用している（未テスト）

nByTown_latlon <- spTransform(nByTown, CRS("+proj=longlat +datum=WGS84"))