замаскировать область океана в файле формы, используя ggplot2 - PullRequest
Новая
       32

замаскировать область океана в файле формы, используя ggplot2

0 голосов
/ 13 июня 2018

Я рисую распределение температуры с помощью ggplo2 на основе файла netcdf и файла формы мира

, но я не могу замаскировать океанскую область, мне нужна картинка, раскрашенная только на суше иоставьте область океана синим, geom_polygon может заполнить только площадь суши, как сделать противоположное заполнение?

thx

вот оригинальный код; 

library(RNetCDF)
library(ggplot2)
library(maptools)
library(rgdal)
library(sp)
setwd("F:\\Programing\\Rnetcdf")
f<-open.nc("Soong.nc")

basemap<-readShapeSpatial("world.shp")
head(basemap)
shore <- fortify(basemap)
#shore$long<-shore$long+180
head(shore)

var.inq.nc(f,"TMP_P0_L109_GLL0")
apt<-var.get.nc(f,"APTMP_P0_L103_GLL0")

lon<-var.get.nc(f,"lon_0")
lat<-var.get.nc(f,"lat_0")

newlon<-rep(lon,each=length(lat))
newlat<-rep(lat,length(lon))
newapt<-as.vector(t(apt))

newlon[newlon>180]<-newlon[newlon>180]-360

draw<-as.data.frame(cbind(newlat,newlon,newapt))

finaldraw<-ggplot(data=draw,aes(x=newlon,y=newlat,fill=newapt))+
  geom_raster(interpolate=F)+
  scale_fill_gradientn(colours = rev(rainbow(7)), na.value = NA) +
  geom_polygon(data = shore, aes(x=long, y = lat, group = group), color = "black",fill=NA) +
  theme_bw()

, и это изображение кода выше

enter image description here

...