contemporary_us <- us_states()
if (require(USAboundariesData, quietly = TRUE) && require(sf, quietly = TRUE)) {
historical_us <- us_states("1820-07-04")
contemporary_ne <- us_states(
states = c(
"Massachusetts",
"Vermont",
"Maine",
"New Hampshire",
"Rhode Island",
"Connecticut"
)
)
historical_ne <- us_states(
as.Date("1805-03-12"),
states = c(
"Massachusetts",
"Vermont",
"Maine",
"New Hampshire",
"Rhode Island",
"Connecticut"
),
resolution = "high"
)
plot(st_geometry(contemporary_us))
plot(st_geometry(historical_us))
plot(st_geometry(contemporary_ne))
plot(st_geometry(historical_ne))
}
Run the code above in your browser using DataLab