Folleto R derivar Polígonos que faltan lat faltan largos
Estoy tratando de trazar el sitio de algunos datos de eventos de enfermedades en un mapa.
Lo uso para importar los datos:
ByTown<-readOGR(dsn="C:/temp/lyme/Towns", layer="Towns", encoding = "UTF-8", verbose= FALSE)
revisa la clase:
class(ByTown)
#getting this result
[1] "SpatialPolygonsDataFrame"
attr(,"package")
[1] "sp"
Luego convierto todos los factores en datos de caracteres y verifico que todavía tengo unSpatialPolygonsDataFrame
utilizandoclass
de nuevo, lo que hago:
Luego formateo los datos que deseo combinar en el mismo caso de título que el original:
townCount$City<-str_to_title(townCount$City)
Luego geo_join los datos de recuento al marco de datos de polígono espacial:
ByTown<-geo_join(ByTown, townCount,"MCD_NAME", "City")
Luego configuro la paleta y ejecuto el mapeo:
pal = colorQuantile("PuOr",ByTown$count, n=5 )
map<-leaflet(ByTown) %>%
addProviderTiles("CartoDB.Positron")%>%
addPolygons(fillColor = ~pal(count),
color = "#000000",
stroke = TRUE,
weight = 1,
smoothFactor = 0.5,
options(viewer = NULL))
map
Y me sale este error:
Error in derivePolygons(data, lng, lat, missing(lng), missing(lat), "addPolygons") :
addPolygons must be called with both lng and lat, or with neither.
He buscado en las ranuras de coordenadas y hay datos allí ... Estoy desconcertado por el error y no encuentro ninguna respuesta útil en línea. Aquí está la cabeza del primer polígono en la ranura de coordenadas:
head(nByTown@polygons[[1]]@Polygons[[1]]@coords )
[,1] [,2]
[1,] 1036519 916318.7
[2,] 1036039 916355.8
[3,] 1031757 916299.7
[4,] 1027474 916244.5
[5,] 1026709 916198.1
[6,] 1026826 916248.3
¿Alguien tiene este problema, identifica la causa raíz y lo arregla?