[R-sig-Geo] Mapping my own polygon?

Kevin Zembower kev|n @end|ng |rom zembower@org
Sat Jun 10 21:27:25 CEST 2023


In my continuing work on reporting on US Census data for my 
neighborhood, I'd like to draw a map of the boundaries of it. I was 
successful in creating and printing an OSM basemap, with the US Census 
blocks that make up my neighborhood on it.

Now, I'd like to create my own polygon, of the boundaries of my 
neighborhood, because the census blocks don't line up exactly with the 
neighborhood boundaries. I need help creating a polygon that I can 
submit to read_osm() that will correctly return an OSM map of my area.

Here's what I've tried so far:
## Reproducible simple example:
library(tidyverse)
library(sf)

rw <- data.frame( ## Simplified neighborhood rectangle
     Longitude = c(-8528150, -8528500, -8528500, -8528150),
     Latitude  = c( 4771475,  4771475,  4771880,  4771880)
)

rw ## Returns (as expected):

 > rw
   Longitude Latitude
1  -8528150  4771475
2  -8528500  4771475
3  -8528500  4771880
4  -8528150  4771880
 >

rw %>%
     st_as_sf(coords = c("Longitude", "Latitude"), dim = "XY") %>%
     st_set_crs(6487) %>% ## CRS 6487 is NAD83 (2011) Maryland in meters
     st_transform(crs = 4269) ## CRS 4269 is NAD83

## Returns:

Simple feature collection with 4 features and 0 fields
Geometry type: POINT
Dimension:     XY
Bounding box:  xmin: 171.777 ymin: 24.65904 xmax: 171.7818 ymax: 24.66314
Geodetic CRS:  NAD83
                    geometry
1 POINT (171.7818 24.66192)
2 POINT (171.7806 24.65904)
3  POINT (171.777 24.66026)
4 POINT (171.7781 24.66314)
 >

I expected the POINTS to look like:
1 POINT (-76.61246 39.35010)

Can anyone suggest what I'm doing wrong? Thanks so much in advance. I've 
worked on it all day today, without making much progress.

-Kevin



More information about the R-sig-Geo mailing list