79 個の変数を含む netcdf ファイルがあり、変数ごとに 365 個のバンドがあります (1 年の各日に 1 つ)。1 つの変数 (つまり、365 レイヤーのラスター) のすべてのバンドを読み取りたいと考えています。「raster」パッケージでは、これは次のように機能します。
dailyvalues <- brick('GLOBAL_2010_day.nc', varname ='WDEP_PREC')
結果は、365 層の RasterBrick です。
> dailyvalues
class : RasterBrick
dimensions : 180, 360, 64800, 365 (nrow, ncol, ncell, nlayers)
resolution : 1, 1 (x, y)
extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +no_defs
source : GLOBAL_2010_day.nc
names : X2010.01.01, X2010.01.02, X2010.01.03, X2010.01.04, X2010.01.05, X2010.01.06, X2010.01.07, X2010.01.08, X2010.01.09, X2010.01.10, X2010.01.11, X2010.01.12, X2010.01.13, X2010.01.14, X2010.01.15, ...
Date : 2010-01-01, 2010-12-31 (min, max)
varname : WDEP_PREC
しかし、「terra」パッケージでこれを行うことに成功していません。私は試した
dailyvalues <- rast('GLOBAL_2010_day.nc')
> dailyvalues
class : SpatRaster
dimensions : 180, 360, 28835 (nrow, ncol, nlyr)
resolution : 1, 1 (x, y)
extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs
sources : GLOBAL_2010_day.nc:WDEP_PREC (365 layers)
GLOBAL_2010_day.nc:WDEP_SOX (365 layers)
GLOBAL_2010_day.nc:WDEP_OXN (365 layers)
... and 76 more source(s)
varnames : WDEP_PREC (WDEP_PREC)
WDEP_SOX (WDEP_SOX)
WDEP_OXN (WDEP_OXN)
...
names : WDEP_PREC_1, WDEP_PREC_2, WDEP_PREC_3, WDEP_PREC_4, WDEP_PREC_5, WDEP_PREC_6, ...
unit : mm, mm, mm, mm, mm, mm, ...
time : 2010-01-01 18:00:00 to 2010-12-31 12:00:00
結果の SpatRaster には 79 個の「ソース」がありますが、1 つの「ソース」を抽出する場合に使用する構文は何ですか? rast 関数に varname = 'WDEP_PREC' を追加しても機能しません。私はdailyvalues$を試しました...しかし、それは単一のレイヤーのみを呼び出します(「名前」の下にリストされています)。