I (on Win7) was trying to get my shapefile data (Here are the data files.) to be displayed using leaflet
package. But without any success. I only get the background map tile but not my data. I am able to see my data in QGIS though. I got the following code and description from here.
library(rgdal)
shapeData <- readOGR(".",'myGIS')
ogrInfo(".",'myGIS')
leaflet() %>% addTiles() %>% setView(lng = -106.363590, lat=31.968483,zoom=11) %>% addPolygons(data=shapeData,weight=5,col = 'red') %>% addMarkers(lng = -106.363590,lat=31.968483,popup="Hi there")
This is what I see in my internet browser. I do not see the shapeData being projected on it:
The answer of Victorp works. However I would suggest to use:
shapeData <- spTransform(shapeData, CRS("+proj=longlat +datum=WGS84 +no_defs"))
This CRS corresponds to EPSG:4326. Leaflet actually takes care of transforming from EPSG:4326 to EPSG:3857 (known as "Google Mercator").
Now, the difference (GRS80 vs WGS84) is negligible (0,1 mm on one of the axis). And it seems that leaflet considers both to be the same Ellipsoid. However for test purposes would be advisable to have EPSG:4326.
If you love us? You can donate to us via Paypal or buy me a coffee so we can maintain and grow! Thank you!
Donate Us With