Я не могу понять, что происходит с этими двумя файлами, я думаю, что это проблема проецирования, или, может быть, мне нужен новый многоугольник.
data:
library(sf)
provinces <- st_read("https://code.highcharts.com/mapdata/countries/ca/ca-all.geo.json")
provinces <- st_transform(provinces, 4326)
pts <- tribble(~Name, ~lat, ~lng,
"Kelowna",49.9,-119.483333,
"Beaumont",53.352436,-113.415313,
"Edmonton",53.55,-113.5,
"Paradise",47.526367,-52.871569,
"Toronto",43.651653,-79.382898)
pts <- st_as_sf(pts, coords = c("lng", "lat"), crs = 4326 )
показываеттаким образом, полигоны явно слишком малы (это должны быть провинции Канады):
leaflet(provinces) %>%
addProviderTiles(providers$CartoDB.Positron) %>%
addPolygons() %>%
addMarkers(data = pts)
Структура показывает, что bbox-ы разные, нопроекции одинаковы:
> print(provinces)
Simple feature collection with 14 features and 22 fields
geometry type: MULTIPOLYGON
dimension: XY
bbox: xmin: -96.01125 ymin: 40.00688 xmax: -95.8891 ymax: 40.08523
epsg (SRID): 4326
proj4string: +proj=longlat +datum=WGS84 +no_defs
> print(pts)
Simple feature collection with 5 features and 1 field
geometry type: POINT
dimension: XY
bbox: xmin: -119.4833 ymin: 43.65165 xmax: -52.87157 ymax: 53.55
epsg (SRID): 4326
proj4string: +proj=longlat +datum=WGS84 +no_defs
# A tibble: 5 x 2
Как мне получить эти два на одной странице?