2013-07-19 2 views
0

У меня есть растровый файл, сгенерированный с пакетом raster из файла netcdf, и мне нужно извлечь значение fillvalue, масштабный коэффициент и смещение. С пакетом ncdf я могу извлечь эти атрибуты со следующим кодом:Извлечение масштабного коэффициента и смещения из растра, сделанного с помощью netcdf в R

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?

ответ

1

У меня нет файла ncdf для проверки, но если вы читаете файл с помощью raster и сохранили его с помощью raster, я считаю, что нижеследующее должно работать. Если вы не попытаетесь использовать readGDAL из пакета rgdal, чтобы прочитать изображение ncdf.

library(raster)  
img <- raster(file_path) 
scale <- gain(img) 
offset <- offs(img) 
fillvalue <- NAvalue(img) 

 Смежные вопросы

  • Нет связанных вопросов^_^