我正在尝试在带有纬度和经度坐标的水平图顶部包含海岸线地图,可能还有国家边界。
到目前为止,我对 levelplot 的代码是:
library(lattice)
library(ncdf4)
ncin <- nc_open("testplot.nc")
tp <- ncvar_get(ncin,"tp")
lon <- ncvar_get(ncin,"longitude")
lat <- ncvar_get(ncin,"latitude")
grid <- expand.grid(lon=lon, lat=lat)
levelplot(tp ~ lon * lat, data=grid)
我得到的图像是这样的:
我怎样才能在情节中包含地图?
更新:我想要地块上位于 37°-31°N 和 34°-31°E 之间的地球切片。 只有海岸线和国家边界,所以您仍然可以看到颜色阴影。
谢谢你。