[R] Recommended ways to draw US Census map on Open Street Map base map?

[Bert Gunter]

As you are using Rstudio stuff(ggplot, tidyverse,etc.), you might try the
RStudio support site here <https://community.rstudio.com/> (note that they
have rebranded under the banner "posit.co").

However, I would have also thought the R-sig-geo list was the best place to
look for such help.

Cheers,
Bert


On Tue, Jun 6, 2023 at 11:38 AM Kevin Zembower via R-help <
r-help using r-project.org> wrote:

> Hello, all,
>
> I asked a version of this question on the R-sig-geo list, but didn't get
> any response. I'm asking here in the hopes of a wider audience.
>
> I'm trying to draw US Census map data, fetched with tigris, on top of a
> base map fetched by the package OpenStreetMap. I'm hoping for the most
> straight-forward solution. I made significant progress with leaflet(),
> but didn't need the interactivity of the map. I just need a 2D, static
> map that I can print and include in a document.
>
> Here's some of what I've tried so far:
> ==========================================
> library(tidyverse)
> library(tigris)
> options(tigris_use_cache = TRUE)
> library(OpenStreetMap)
> library(ggplot2)
>
> ## Get an Open Street Map:
> rw_map <- openmap(nw, se,
>                    type = "osm",
>                    mergeTiles = TRUE) %>%
>      openproj(projection = "+proj=longlat +ellps=WGS84 +datum=WGS84
> +no_defs")
>
> ## Get an example census map:
> rw_tract <- tracts(state = "MD",
>                      county = "Baltimore city",
>                      year = "2020") %>%
>      filter(NAME == "2711.01")
>
> ## This works:
> autoplot.OpenStreetMap(rw_map)
>
> ## So does this:
> plot(rw_tract$geometry)
>
> ## These don't:
> autoplot.OpenStreetMap(rw_map) +
>      geom_sf(rw_tract$geometry)
>
> ggplot(map_data(rw_map), aes(long, lat))
>
>
> ggplot(aes(x="long", y="lat")) +
>      geom_sf(rw_map$geometry)
> =====================================================
>
> I think my problem in part is failing to fully understand the formats of
> the rw_map and rw_tract containers. rw_tract says it's a simple feature
> collection, but rw_map just gives me lists of the data.
>
> Can anyone help nudge me along in getting my rw_tract to be drawn on my
> rw_map? Any advice or guidance on putting together map data from
> different sources?
>
> And an over-arching question: Is moving in this direction, with ggplot2,
> the way you would recommend accomplishing this task? Is there a simpler,
> more straight-forward way of doing this?
>
> Thanks in advance for your help and efforts.
>
> -Kevin
>
> ______________________________________________
> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

	[[alternative HTML version deleted]]