2016-11-02 7 views
3

リーフレットパッケージで多数のポリゴンを描画しようとしていますが、何がうまくいかないのか分かりません。R /リーフレット - 多数のポリゴンを描画

私が使用してシェープファイルは、ここで見つけることができます:https://www.data.gouv.fr/en/datasets/fond-de-carte-des-codes-postaux/

library(leaflet) 
library(rgdal) 
df <- readOGR("C:/Users/me/codes_postaux","codes_postaux_region") 
plot(df) 

output shapefile

シェープファイルは、私にはOKと思われると私が使用するコードはかなり単純です。しかし、私は出力としてのみ、ポリゴンはマップを取得します。私はかなり長い間、この問題に苦労してきました、誰かが私をここで助けることができたら本当に感謝します。

map <- leaflet(df) %>% 
    addProviderTiles("CartoDB.Positron")%>% 
    fitBounds(10,38,10,55) %>% 
    addPolygons(fillOpacity = 0.8, color = "Blue", weight = 1) 

map 

leaflet

+0

それは、このデータで動作します: 'DF < - ラスタ::のgetData( "GADM"、国= "FRA"、レベル= 4)'?あなたのデータかパッケージ関連の問題かどうか疑問に思っています。 – TimSalabim

答えて

4

[email protected]を見て、plot(df); axis(1); axis(2)の出力。シェイプファイルは特定のCRSを使用します。 SpatialPolygonsDataFrameを共通の​​(私はここからCRSコードを得ました:Leaflet for R: Raster Images)に変換する必要があります。

library(sp) 

pj <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs") 
df2 <- spTransform(df, pj) 

map2 <- leaflet(df2) %>% 
    addProviderTiles("CartoDB.Positron")%>% 
    fitBounds(10,38,10,55) %>% 
    addPolygons(fillOpacity = 0.8, color = "Blue", weight = 1) 
map2 

enter image description here

+0

ありがとうございます。これはうまくいきます。 CRS()で射影引数を見つけた方法を指定してください。 –

+0

@PierreDudek;私はそれが最も一般的なものであることを知っていたので、私はそれを使用しました(私は答えたときにそれを検索しませんでした)。しかし、今ではRのリーフレット記述でCRSを見つけました。CRSをCRSに書き直しました(その違いは単なる表現です)([R:Raster Imagesのリーフレット](https://rstudio.github.io/leaflet/)参照)。 raster.html)) – cuttlefish44

関連する問題