https://stat.ethz.ch/pipermail/r-sig-geo/2023-June/029294.html
rw <- data.frame( ## Simplified neighborhood rectangle
Longitude = c(-8528150, -8528500, -8528500, -8528150),
Latitude = c( 4771475, 4771475, 4771880, 4771880)
)
prj <- "EPSG:3857"
library(vapour)
library(ximage) ## hypertidy/ximage
library(dsn) ## hypertidy/dsn
im <- gdal_raster_image(dsn::wms_virtualearth_street(), target_ext = c(range(rw$Longitude), range(rw$Latitude)) + c(-1, 1, -1, 1) * 5e4,
target_crs = prj, target_dim = dev.size("px"), resample = "cubic")
op <- par(mar = rep(0, 4))
ximage(im, asp = 1)
abline(v = rw$Longitude, h = rw$Latitude)
par(op)
library(sf)
rw %>%
st_as_sf(coords = c("Longitude", "Latitude"), dim = "XY") %>%
st_set_crs(3857) %>%
st_transform(crs = 4269) ## CRS 4269 is NAD83
Simple feature collection with 4 features and 0 fields
Geometry type: POINT
Dimension: XY
Bounding box: xmin: -76.61282 ymin: 39.34683 xmax: -76.60968 ymax: 39.34965
Geodetic CRS: NAD83
geometry
1 POINT (-76.60968 39.34683)
2 POINT (-76.61282 39.34683)
3 POINT (-76.61282 39.34965)
4 POINT (-76.60968 39.34965)
Thanks, Mike, for your work to help me with my problem.
Here's code to generate a map of my neighborhood, with dimensions in degrees:
Here's how to create the map I used to get the coordinates of the corners of the borders of my neighborhood:
The area you found is about 13 miles due north of my neighborhood:
![Screenshot from 2023-06-11 11-19-37](https://user-images.githubusercontent.com/1446154/244947623-08975092-ff58-45a0-86b9-a38c9d9100c8.png)
Thanks, again, for all your help.
-Kevin