0

我正在尝试在带有纬度和经度坐标的水平图顶部包含海岸线地图,可能还有国家边界。

到目前为止,我对 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 之间的地球切片。 只有海岸线和国家边界,所以您仍然可以看到颜色阴影。

谢谢你。

4

0 回答 0