我想使用ggplot和facet_wrap生成12张Map，这些Map显示了1月至12月苏格兰大陆架的海洋温度。我得到了一个包含25，000个观测值的“.csv”文件，其具有2016-2018年的每月温度数据以及相关联的纬度和经度值。我曾尝试使用geom_rasterMap此数据，但收到以下错误消息

Error: cannot allocate vector of size 39.7 Gb
In addition: Warning messages:
1: In f(...) :
 Raster pixels are placed at uneven horizontal intervals and will be shifted. Consider using geom_tile() instead.
2: In f(...) :
 Raster pixels are placed at uneven vertical intervals and will be shifted. Consider using geom_tile() instead

当然，我也尝试过使用geom_tile，但当我运行代码时，我得到了一个没有颜色的空白Map，尽管我已经指定了填充和颜色。下面是一些类似于我的 Dataframe 的示例数据，请注意，我的latitude和longitude值的真实的数据不是由seq(...,...,0.1)均匀隔开的（抱歉，我不知道如何在不进行排序的情况下生成示例数据），因此您不会遇到与我相同的错误，代码也将正常工作

#In my dataframe the lat and long values are not equally spaced!
mapoc_temp = expand.grid(data.frame(Longitude= seq(-64.5,-62.4,0.1), 
                                    Latitude= seq(42.7,44.8,0.1),
                                    year = sample(c(2016,2017,2018), 22, replace = T),
                                    month = sample(month.abb, 22, replace = T)))
mapoc_temp$Temp = runif(nrow(mapoc_temp))

下面是geom_raster的代码，它给出了我提到的错误

library(mapdata)
library(ggplot2)

#map I'm using for ggplot
canada = map_data("worldHires", "Canada")

ggplot(mapoc_temp, aes(x=Longitude, y=Latitude)) + 
#try to map my temp on ocean
  geom_raster(aes(fill = Temp, x = Longitude), interpolate = TRUE) +
  geom_polygon(data = canada, aes(x=long, y=lat, group=group), colour="grey50", fill = 'grey55')+
  coord_sf(xlim=c(-64.5,-62.8), ylim=c(42.7,45)) +
#get months
facet_wrap(vars(month))

这是我尝试geom_tile时的代码，我再次注意到它在我的样本数据中有效，但在我的真实的数据中无效，我相信这与我的坐标之间的距离不相等有关。

ggplot(mapoc_temp, aes(x=Longitude, y=Latitude)) + 
  geom_tile(aes(fill = Temp, x = Longitude, y = Latitude), colour = mapoc_temp$Temp) +
  geom_polygon(data = canada, aes(x=long, y=lat, group=group), colour="grey50", fill = 'grey55')+
  coord_sf(xlim=c(-64.5,-62.8), ylim=c(42.7,45)) +
  facet_wrap(vars(month))

这是使用geom_tile获得的图片

这是我想要制作的大致内容，但显然是12张Map，因为我希望每个月都有一张分辨率更高的Map（如果可能，但在这一点上，我会选择任何用我的温度数据着色的Map）。

我有一种感觉，我必须沿着这类代码（我不久前发现的代码片段）做更多的事情，但我已经尝试过操纵它，没有成功。有什么建议吗？

#my CRS value!
latlong = "+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0"

#for transforming
planar ="+proj=utm +zone=20 +datum=NAD83 +units=km +no_defs +ellps=GRS80 +towgs84=0,0,0"

    out <- x%>%
      st_as_sf(coords=c("Longitude","Latitude"),crs=latlong)%>%
      st_transform(planar)%>%
      st_coordinates()%>%
      raster::rasterize(.,grid,field=xyz[,3],fun=rasterFun)%>%
      raster::projectRaster(.,crs=latlong)%>%
      raster::rasterToPolygons(.)%>% # this part is slow
      st_as_sf()%>%
      rename(MAP = layer)