[R-sig-Geo] writeRaster - netcdf - Rotation

[Philipp Reiter]

Hey there,

I’ve got a problem concerning exporting a raster*-Object as nc-File:
When exporting rasterLayers or rasterBricks as nc-Files with writeRaster from the raster-Package, easting and northing in the resulting nc-File are switched, which results in a rotated plot. I can correct this afterwards, but it should be possible to be correct right away?!

Example code for a rasterLayer:
tmp <- raster(readGDAL(‘example.asc’,p4s="+init=epsg:31467"))
writeRaster(tmp, filename=’export.nc’, format="CDF", overwrite=TRUE)
nc <- open.ncdf("export.nc")
x <- get.var.ncdf(nc,”easting”)
y <- get.var.ncdf(nc,”northing”)

R version 2.15.3 (2013-03-01)
Platform: x86_64-w64-mingw32/x64 (64-bit)
ncdf_1.6.6    raster_2.1-16 rgdal_0.8-5   sp_1.0-6


I hope someone can give me a hint! Thanks in advance!
Philipp