Hi,
This is a great place to ask - and nice reproducible code!
You have a two issues as far as I can see. First it seems you have dropped
the "4" in "WGS84". And second you want to transform rw_pop, not the
output of leaflet(). The following works for me.
## Start
library(tidyverse)
library(tidycensus)
library(leaflet)
library(sf)
rw_blocks <- c(3000, 3001, 3002, 3005, 3006, 3007, 3008, 3009, 3010, 3011)
rw_pop <- get_decennial(
geography = "block",
variables = "P1_001N",
year = 2020,
state = "MD",
county = "Baltimore city",
geometry = TRUE
) %>%
filter(substr(GEOID, 6, 11) == "271101" &
substr(GEOID, 12, 15) %in% rw_blocks
) ## %>% st_transform('+proj=longlat +datum=WGS8')
(rw_pop_map <- rw_pop %>%
st_transform('+proj=longlat +datum=WGS84') %>%
leaflet() %>%
fitBounds(-76.616, 39.352, -76.610, 39.346) %>%
addTiles() %>%
addPolygons()
)
## End
On Wed, May 31, 2023 at 4:03 PM Kevin Zembower via R-sig-Geo <
r-sig-geo@r-project.org> wrote:
> Hello, all. Newbie to sf, tidycensus and the tidyverse here.
>
> First off, is this the appropriate list to ask this question? If not,
> let me know and I'll go away.
>
> I'm trying to map census blocks for my neighborhood to a base map. I'm
> using tidycensus to get the geometry of the census blocks, and leaflet
> to map them to the OSM base maps. Mostly, this is going really well, and
> I'm very pleased with the speed of development (I just started this
> morning) and results.
>
> However, I get this error:
>
> Warning message:
> sf layer has inconsistent datum (+proj=longlat +datum=NAD83 +no_defs).
> Need '+proj=longlat +datum=WGS84'
>
> I think I need to use st_transform, but can't get it to work.
>
> Here's a reproducible example, with some commented out lines of what
> I've tried:
>
> ## Reproducible example:
> library(tidyverse)
> library(tidycensus)
> library(leaflet)
> library(sf)
>
> rw_blocks <- c(3000, 3001, 3002, 3005, 3006, 3007, 3008, 3009, 3010, 3011)
>
> rw_pop <- get_decennial(
> geography = "block",
> variables = "P1_001N",
> year = 2020,
> state = "MD",
> county = "Baltimore city",
> geometry = TRUE
> ) %>%
> filter(substr(GEOID, 6, 11) == "271101" &
> substr(GEOID, 12, 15) %in% rw_blocks
> ) ## %>% st_transform('+proj=longlat +datum=WGS8')
>
> (rw_pop_map <- rw_pop %>%
> leaflet() %>%
> ## st_transform('+proj=longlat +datum=WGS8') %>%
> fitBounds(-76.616, 39.352, -76.610, 39.346) %>%
> addTiles() %>%
> addPolygons()
> )
> ## Error occurs when executing above block
>
> Can anyone offer me a hint as to how to resolve this error?
>
> Thanks so much for any advice and guidance.
>
> -Kevin
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
