raster icon indicating copy to clipboard operation
raster copied to clipboard

Error using `compression` argument while writing a netCDF file (`writeRaster`)

Open elgabbas opened this issue 1 year ago • 0 comments

Hello, I would like to export a raster object into netCDF file. However, I receive an error when using compression argument.

# raster version: 3.6-26     # ncdf4 version: 1.22
require(raster); require(ncdf4)

Map <- raster::raster(
  nrows = 404, ncols = 390,  xmn = 2630000, xmx = 6530000, 
  ymn = 1380000, ymx = 5420000, crs = "epsg:3035", vals = 1:157560)

raster::writeRaster(Map, filename = "Map1.nc", format = "CDF")
file.size("Map1.nc") # 637396
raster::writeRaster(Map, filename = "Map2.nc", format = "CDF", compression = 5)

A smaller file was created, but with the following error:

# Error in R_nc4_open: NetCDF: Unknown file format
# Error in ncdf4::nc_open(filename, readunlim = FALSE, suppress_dimvals = TRUE) : 
# Error in nc_open trying to open file D:\BioDT_IAS\Map2.nc (return_on_error= FALSE )

file.size("Map2.nc") # 201879

Setting options(HDF5_USE_FILE_LOCKING=FALSE), as suggested here, did not help.

I tried using terra:: but found an issue with the CRS information.

Any suggestions on how to avoid this issue?

elgabbas avatar Mar 14 '24 14:03 elgabbas