【发布时间】：2020-10-23 07:55:17
【问题描述】：
使用的库：
library(sp)
library(sf)
library(ggplot2)
library(ggmap)
创建了一个名为“coordinate.data”的数据框，其中经度和纬度作为列名，气象站位置作为行名。
longitude <- c(-73.964482,-73.953678,-73.893522,-73.815856,-74.148499)
latitude <- c(40.767544,40.631762,40.872481,40.734335,40.604014)
coordinate.data <- data.frame(longitude,latitude)
rownames(coordinate.data) <- c("MANH","BKLN","BRON","QUEE","STAT")
然后我检索了新泽西县和纽约市行政区的 shapefile 数据，并删除了所有不必要的列，因此两个 shapefile 中只剩下几何字段。 NYC Boroughs shapefile 数据是从 NYC Open Data 下载的，而 NJ 县边界是从 NJGIN Open Data 下载的。
nj.shp <- st_read("~/Downloads/NJ/NJ_Counties.shp")
nj <- nj.shp[,-(1:21)]
nyc.shp <- st_read("~/Downloads/NY/NYC_Boroughs.shp")
nyc <- nyc.shp[,-(1:4)]
我将两个 shapefile 格式化为具有相同的投影（ESPG 代码 3857）并将它们组合成一个 shapefile 数据框，其中包含一个变量（几何）中的 26 个观测值（县/区）。
same.projection <- CRS("+init=EPSG:3857")
nj.data <- st_transform(nj,same.projection)
new.projection <- CRS("+init=EPSG:3857")
nyc.data <- st_transform(nyc,new.projection)
combined.data <- rbind(nj.data,nyc.data)
我现在正尝试在地图上绘制组合 shapefile（“combined.data”）以及气象站位置（“coordinate.data”）。 当我尝试这样做时，它不可避免地会运行并且 R 会关闭。如果我删除 geom_sf(...)，它会正确绘制站点并格式化所有内容，所以我认为问题出在这行代码上。
mesonet.map <-ggplot() +
ggtitle("NY Mesonet Site Locations") +
xlab("Longitude") +
ylab("Latitude") +
geom_point(data=coordinate.data,aes(x=longitude,y=latitude))+
geom_text(aes(x=longitude,y=latitude,label=rownames(coordinate.data)),size=3.25,nudge_y=0.02)+
geom_sf(data=combined.data,fill='darkgreen') +
mesonet.map + theme(
panel.background=element_rect(fill="lightblue",color="lightblue",size=0.5,linetype="solid"),
panel.grid.major=element_line(size=0.5,linetype='solid',color="white"),
panel.grid.minor=element_line(size=0.25,linetype='solid',color="white")
)
【问题讨论】：
-
除了一些小问题，您的代码在我的机器上运行良好。