我正在尝试使用包含多边形的传单创建交互式绘图。 这些多边形的坐标不是经度/纬度,而是代表建筑平面图的笛卡尔x和y坐标。
对于单位数字的坐标,结果看起来很好:
library(leaflet)
library(sp)
Sr1 = Polygon(cbind(c(9, 8, 8, 9, 9),
c(3, 3, 4, 4, 3)))
Srs1 = Polygons(list(Sr1), "s1")
SpP = SpatialPolygons(list(Srs1), 1L)
leaflet() %>% addPolygons(data = SpP)
但是,当我将所有y坐标增加100时,就会出错:
library(leaflet)
library(sp)
Sr1 = Polygon(cbind(c(9, 8, 8, 9, 9),
c(103, 103, 104, 104, 103)))
Srs1 = Polygons(list(Sr1), "s1")
SpP = SpatialPolygons(list(Srs1), 1L)
leaflet() %>% addPolygons(data = SpP)
我猜测这与地图的默认投影设置有关。我尝试使用以下设置初始化地图:
leaflet(options(crsClass = "L.CRS.Simple"))
阅读this示例但未成功。