Я пытаюсь включить карту береговой линии, а также, возможно, и границы страны, поверх графика уровня с координатами широты и долготы.
Код, который у меня так далеко за уровень уровня:
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 ° с.ш. и 34 ° -31 ° в.д. на участке. Только береговые линии и границы страны, так что вы все еще можете видеть цветное затенение.
Спасибо.