我有一个netCDF文件,我打算从每年850 hpa的压力高度中提取点数据。但是,在使用raster :: rotate()校正x轴后,会出现一行缺失数据,如下图所示。
任何人都可以对此有所了解吗?
数据可用here。
library(maptools)
library(raster)
data(wrld_simpl)
r <- raster("z850.nc", band = 5)
r <- rotate(r) ## Corrects extent
plot(r)
plot(wrld_simpl, add = T)