看到 CRAN 手册,我认为不可能完成这项工作。我的解决方法是手动创建一个覆盖指南针的多边形。这是一项乏味的工作,但这也许是现在要走的路。第一步是创建一个多边形。extent()
为您提供经度和纬度的最小值和最大值。我使用 xmax 和 ymax 为多边形指定五个点。(我使用了值并找到了最佳值。)我们需要为多边形分配一个正确的投影。如果您ny.map
在 R 控制台中键入,您将看到coord. ref.
. 这是您需要的投影。tm_shape()
然后,我使用和将多边形添加到地图中tm_fill()
。
lon <- c(extent(ny.map)[2]-2100, extent(ny.map)[2]-250, extent(ny.map)[2]-250, extent(ny.map)[2]-2100, extent(ny.map)[2]-2100)
lat <- c(extent(ny.map)[4]-250, extent(ny.map)[4]-250, extent(ny.map)[4]-2100, extent(ny.map)[4]-2100,extent(ny.map)[4]-250)
foo <- SpatialPolygons(list(Polygons(list(Polygon(cbind(lon, lat))), ID = 1)),
proj4string = CRS("+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"))
tm_shape(shp = ny.map) +
tm_raster() +
tm_shape(shp = foo) +
tm_fill("red", alpha = 0.4) +
tm_compass(position = c("right", "top"),
type = "4star",
show.labels = 2) +
tm_layout(title = "New York",
title.bg.color = "white",
title.bg.alpha = 0.5)