Как создать еженедельную композицию из 5 дней подряд

Question

У меня есть файл NetCDF солености в воде Индонезии с 4 измерениями (долгота, широта, глубина и время). Как создать еженедельный композит из моих данных, загрузите данные здесь: https://onedrive.live.com/redir?resid=6FFDD661570C7D0A%21177 выходную карту здесь: https://onedrive.live.com/redir?resid=6FFDD661570C7D0A%21176

Я хотел бы преобразовать растр в вектор и применить, чтобы получить среднее значение, но у меня есть проблема, чтобы построить векторные данные, используя rasterVis

Answer 1

С твоим примером ни правда сложно:

# load needed librairies
library(rasterVis)

# open the data
salinity <- brick("data.nc", varname = "salinity")

salinity
# class       : RasterBrick 
# dimensions  : 61, 61, 3721, 5  (nrow, ncol, ncell, nlayers)
# resolution  : 0.08333333, 0.08333333  (x, y)
# extent      : 104.9583, 110.0417, -5.041667, 0.04166667  (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +datum=WGS84 
# data source : data.nc 
# names       : X252331200, X252417600, X252504000, X252590400, X252676800 
# z-value     : 252331200, 252417600, 252504000, 252590400, 252676800 
# varname     : salinity 
# level       : 1 

# Calculate the mean
m.salinity <- mean(salinity)

m.salinity

# class       : RasterLayer 
# dimensions  : 61, 61, 3721  (nrow, ncol, ncell)
# resolution  : 0.08333333, 0.08333333  (x, y)
# extent      : 104.9583, 110.0417, -5.041667, 0.04166667  (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +datum=WGS84 
# data source : in memory
# names       : layer 
# values      : 18.85652, 31.84299  (min, max)