我试图在传单（或地图视图）中显示多个图层，其中之一是 EPSG:27700 中的栅格。我设法充分覆盖这些图层的唯一方法是通过默认的经纬度投影，这意味着栅格的重新投影及其插值。我无法在这个项目中进行插值，所以我需要在 EPSG:27700 上工作。
如何向未投影的栅格显示附加图层？我尝试过使用 CRS.Simple，因为我想在简单的笛卡尔计划中显示所有内容，但没有成功。我不介意失去美丽的背景瓷砖。但无论我尝试什么，我都无法让我的多边形（也包括 EPSG27700）图层（或任何 sp 对象）与我的未插值栅格正确显示。我希望下面的 MWE 能够有效地说明我的问题：
library("raster")
library("leaflet")
library("eurostat")
library("sf")
## get UKK spdf projected on british grid EPSG27700
europe <- get_eurostat_geospatial(resolution = 10, nuts_level = 1, year = 2021)
UK_spdf <- as_Spatial(europe[grepl("UK", europe$id),])
UK_spdf <- spTransform(UK_spdf, crs("+init=epsg:27700 +units=km +datum=WGS84"))
## build a dummy raster projected on EPSG:27700
r <- rasterize(UK_spdf, raster(UK_spdf, ncols = 100, nrows = 200))
## the two layers overlay well in default plots
plot(r) ; plot(UK_spdf, add=TRUE)
## raster can …