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

```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
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

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
```