在R中将更改传单更改为不同的坐标系

时间:2017-04-17 18:39:43

标签: r leaflet tiles coordinate-systems map-projections

我正试图在R中的Leaflet地图中的荷兰Rijksdriehoekstelsel坐标中获取一张地图,该地图将在稍后的Shiny应用程序中显示。我将不得不稍后添加多边形,并将其X和Y坐标更改为WGS84将是可怕的。我从互联网上窃取了proj4def和分辨率,然而,我只是得到一个带有放大和缩小按钮的灰色屏幕。我搞砸了leafletOptions或tileLayer。我此时已尝试过多个urlTemplates,我真的不知道可能出错了什么。

如果我使用L.Proj.CRS.TMS我只是得到一个看起来更糟的白色屏幕。

library(leaflet)

leaflet(
 options =
    leafletOptions(
      crs=leafletCRS(
        crsClass="L.Proj.CRS",
        code='EPSG:28992',
        proj4def='+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812 +no_defs',
    resolutions = c(3440.640, 1720.320, 860.160, 430.080, 215.040, 107.520, 53.760, 26.880, 13.440, 6.720, 3.360, 1.680, 0.840, 0.420),
      origin=c(0,0))
  )) %>%
  addTiles(urlTemplate='http://geodata.nationaalgeoregister.nl/tms/1.0.0/brtachtergrondkaart/{z}/{x}/{y}.png', 
  tileOptions(minZoom=3,maxZoom=13,tms=TRUE,continuousWorld=TRUE)   ) %>%
  setView(52,5.1,13)

这是Javascript中的一个例子,它做了同样的事情: http://bl.ocks.org/webmappergists/raw/1ceac1353920b3810486/

0 个答案:

没有答案