[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