Я пытаюсь экспортировать некоторые файлы NetCDF проекта CMIP5 в GeoTiff, как указано здесь , но полученные изображения имеют смещение.Вот код, который я использую:
library(ncdf4)
library(raster)
FILENAME_IN = "thetao_Omon_HadGEM2-ES_rcp45_r1i1p1_200512-201511.nc"
nc <- nc_open(FILENAME_IN)
summary(nc)
print(nc)
print(names(nc$var))
varsal <- brick(FILENAME_IN, varname="thetao", level=2, stopIfNotEqualSpaced=FALSE)
summary(varsal)
print(varsal)
r <- stack(varsal)
print(r)
names(r)
extent(r) <- extent(0, 360, -90, 90)
rr <- rotate(r)
print(rr)
names(rr)
writeRaster(rr, filename="thetao_Omon_HadGEM2-ES_rcp45_r1i1p1_lvl2.tif", format="GTiff",overwrite=TRUE)
Просмотр изображения:
