我有以下代码，可以创建纽约市县061和047的Map：

library(tidycensus)
library(tidyverse)
library(tigris)

ny_data <- get_acs(
  geography = "tract",
  variables = "B19013_001",
  state = "NY",
  county = c("061","047"),
  year = 2019,
  geometry = TRUE,
  cb = FALSE) %>%
  st_transform(26918) %>%
  erase_water(year = 2019)

ggplot(data=ny_data, aes(fill = estimate)) + 
  geom_sf(color = NA) + 
  scale_fill_viridis_c(option = "magma") +
  theme_minimal() +
  theme(panel.grid.major = element_blank(), 
        panel.grid.minor = element_blank(),
        axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))

NYC map1
但是，我想稍微放大一下这张Map。我知道我应该使用coord_sf函数。但是，当我将函数添加到ggplot时，我得到了意想不到的结果：

ggplot(data=ny_data, aes(fill = estimate)) + 
  geom_sf(color = NA) + 
  coord_sf(xlim = c(-74.05, -73.9), ylim = c(40.6, 40.85)) +
  scale_fill_viridis_c(option = "magma") +
  theme_minimal() +
  theme(panel.grid.major = element_blank(), 
        panel.grid.minor = element_blank(),
        axis.text.x = element_text(angle = 90, vjust = 0.5, hjust=1))

NYC map2
我怀疑这可能与坐标系统有关，但我不熟悉地理空间数据，也不知道如何修复它。这也可能与“st_transform”函数有关，该函数用于帮助我擦除水域。如果你能帮忙的话，我将不胜感激。
摘要：我试图使用coord_sf函数放大Map区域，经度和纬度的x和y限制不像预期的那样工作。