在现有的shapefile map ggplot r中组合osm数据(使用osmdata包下载)

# 安装并加载必要的包
install.packages("ggplot2")
install.packages("osmdata")
library(ggplot2)
library(osmdata)

# 下载OSM数据
bbox <- getbb("New York City")  # 获取纽约市的边界框
osm_data <- opq(bbox) %>%
  add_osm_feature(key = "building", value = "residential") %>%
  osmdata_sf()  # 下载住宅建筑数据

# 读取shapefile
shapefile_path <- "path_to_your_shapefile.shp"
shapefile_data <- st_read(shapefile_path)

# 绘制地图
ggplot() +
  geom_sf(data = shapefile_data, fill = "white", color = "black") +
  geom_sf(data = osm_data$osm_points, color = "blue", size = 2) +
  geom_sf(data = osm_data$osm_lines, color = "red") +
  geom_sf(data = osm_data$osm_polygons, fill = "green", alpha = 0.5) +
  theme_minimal() +
  labs(title = "Combining Shapefile and OSM Data in R")

shapefile_data <- st_transform(shapefile_data, crs = 4326)  # 转换为WGS84坐标系

ggplot() +
  geom_sf(data = osm_data$osm_polygons, fill = "green", alpha = 0.5) +
  geom_sf(data = shapefile_data, fill = "white", color = "black") +
  geom_sf(data = osm_data$osm_lines, color = "red") +
  geom_sf(data = osm_data$osm_points, color = "blue", size = 2) +
  theme_minimal() +
  labs(title = "Improved Visualization")