我有一个使用netcdf文件中的raster
包生成的光栅文件,我需要提取fillvalue,比例因子和偏移量。使用包ncdf
,我可以使用下一个代码提取这些属性:
fillvalue <- att.get.ncdf(ncdf,"tp","_FillValue")
scale <- att.get.ncdf(ncdf,"tp","scale_factor")
offset <- att.get.ncdf(ncdf,"tp","add_offset")
但是使用raster
包我不知道如何提取它们,是否可能?或者我必须使用ncdf
?
答案 0 :(得分:1)
我没有要测试的ncdf文件,但是如果您使用raster
阅读文件并使用raster
保存,我相信以下内容应该有效。如果没有尝试使用readGDAL
包中的rgdal
来阅读ncdf图片。
library(raster)
img <- raster(file_path)
scale <- gain(img)
offset <- offs(img)
fillvalue <- NAvalue(img)