У меня есть файл NETCDF с такими атрибутами, как долгота, широта, время, осадки. Данные охватывают определенную пространственную область. Это ежедневные данные с 1960 по 2100 год.
1) Я хотел бы разделить данные пространственно (например, широта [45,50] и долгота [-78, -85]) из основного домена
2) Из подмножества я хотел бы усреднить по всем сеткам и создать ежедневный временной ряд с одним столбцом, а затем записать его в файл .csv.
NB: мои данные содержат недостающие значения
Что-то в этом роде должно работать
library(raster)
b <- brick('file.nc')
be <- crop(b, extent(-85, -78, 45, 50))
a <- aggregate(be, dim(be)[2:1], na.rm=TRUE)
v <- values(a)
write.csv(v, 'precip.csv', row.names=FALSE)
Обычно, чтобы получить дату:
date <- getZ(be)
Or
date <- names(a)
date <- gsub('X', '', date)
А потом
v <- data.frame(date=date, prec=v)
write.csv(v, 'precip.csv', row.names=FALSE)
Доступна ли эта дата для чтения человеком, зависит от того, как она хранится в файле ncdf (т.е. соблюдаются ли определенные соглашения или нет).
Это сделает сокращение и усреднение, но записывает ответ в другой netcdf. Если вам действительно нужен CSV, вам нужно будет использовать эту часть решения выше.
cdo fldmean -sellonlatbox,-85,-78,45,50 in.nc out.nc
