R функция addPolygons

Question

Я пытаюсь добавить на карту несколько полигонов, но не знаю, почему мой код не работает. Я использовал тот же код раньше для других файлов .shp, и он отлично работает. Я получил пространственный файл отсюда: http://library.queensu.ca/geo/health_regions . Он находится в Ontario-> Local Health Integrated Networks -> 2006 -> Boundary Files.

  lhin<- readOGR("lhin_divisions/HR035B08.shp",layer = "HR035B08")
  leaflet() %>% addTiles()%>%setView(lng=-81,lat=45,zoom=6)
  %>%addPolygons(data=lhin,stroke=TRUE, weight=1, col='red')

Я чувствую, что мне не хватает чего-то простого, но я не могу понять, что?

Полученная карта должна иметь указанные здесь границы:http://library.queensu.ca/madgic/free/health_regions/2005/images/lhin_ont.pdf

Accepted Answer

эта работа. Вы можете изменить цвет плитки, я просто использовал черный и белый, так как не знал, что ищу. Шейп-файл был в формате utm без проекции. Поэтому я установил проекцию и конвертировал в лонглат. proj4string(lhin) <- CRS("+proj=utm +ellps=WGS84 +datum=WGS84") lhin <- spTransform(lhin, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")) leaflet() %>% addProviderTiles("Stamen.Toner") %>% addPolygons(data=lhin,stroke=TRUE, weight=1, color = ~"red")

R функция addPolygons

Ответы

Ответ 1 · рейтинг 2 · MLavoie

Похожие материалы