使用R的raster包，我从文件中获得brick，带有以下ncdump标题（我显示一个小示例文件，实际文件要大得多）：
dimensions:
lon = 2 ;
lat = 3 ;
time = UNLIMITED ; // (125000 currently)
variables:
float lon(lon) ;
lon:standard_name = "longitude" ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
lon:axis = "X" ;
float lat(lat) ;
lat:standard_name = "latitude" ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
lat:axis = "Y" ;
double time(time) ;
time:standard_name = "time" ;
time:long_name = "Time" ;
time:units = "seconds since 2001-1-1 00:00:00" ;
time:calendar = "standard" ;
time:axis = "T" ;
short por(time, lat, lon) ;
por:_FillValue = 0s ;
por:missing_value = 0s ;
在R：
class : RasterBrick
dimensions : 3, 2, 6, 125000 (nrow, ncol, ncell, nlayers)
resolution : 0.008999825, 0.009000778 (x, y)
extent : 6.4955, 6.5135, 44.0955, 44.1225 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
data source : /home/clima-archive/afantini/chym/chym_output/test.nc
names : X0, X3600, X7200, X10800, X14400, X18000, X21600, X25200, X28800, X32400, X36000, X39600, X43200, X46800, X50400, ...
z-value : 0, 449996400 (min, max)
varname : por
但是，为了更快地访问和更高压缩，已经交换了两个文件维度，因此对于我们需要的那种使用，分块更好。所以文件将是这样的（link to download the 1MB file）：
dimensions:
lon = UNLIMITED ; // (2 currently)
lat = 3 ;
time = 125000 ;
variables:
float lon(lon) ;
lon:standard_name = "longitude" ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
lon:axis = "X" ;
float lat(lat) ;
lat:standard_name = "latitude" ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
lat:axis = "Y" ;
double time(time) ;
time:standard_name = "time" ;
time:long_name = "Time" ;
time:units = "seconds since 2001-1-1 00:00:00" ;
time:calendar = "standard" ;
time:axis = "T" ;
short por(lon, lat, time) ;
por:_FillValue = 0s ;
por:missing_value = 0s ;
在R：
class : RasterBrick
dimensions : 3, 125000, 375000, 2 (nrow, ncol, ncell, nlayers)
resolution : 3600, 0.009000778 (x, y)
extent : -1800, 449998200, 44.0955, 44.1225 (xmin, xmax, ymin, ymax)
coord. ref. : NA
data source : /home/clima-archive/afantini/chym/chym_output/test_swapped.nc
names : X6.5, X6.50899982452393
degrees_east: 6.5, 6.50899982452393
varname : por
正如您所看到的那样，文件打开就好像列数为125000.我想将列数与层数交换，而不读取所有数据。我想从光栅手册中我应该使用layer或lvar，因为：
图层：整数。要在多层文件中使用的图层（变量）， 或从RasterStack / Brick或提取的层 SpatialPixelsDataFrame或SpatialGridDataFrame。一个空的 如果'layer = 0'
，则返回RasterLayer（无关联值）.......
'lvar'：整数＆gt; 0（默认= 3）。选择'级别变量' （第三维变量），如果文件有4个维度 （例如深度而不是时间）
但这似乎不起作用，例如layer="time"，因为它什么都没改变。
我该怎么做？
答案 0 :(得分：3)
如果您不打算在打开/阅读后重塑，我认为您可以使用ncdf4库读取变量中的数据然后转置它。类似的东西：
nc <- nc_open(*your_nc_file*)
data <- ncvar_get(nc, por) # "por" is the name of your variable, right ?
data_new <- aperm(data, c(1,3,2)) # "transpose" the matrix
可能的问题可能是data_new不再是raster*对象，但您可以轻松地从中重新创建一个{{1}}对象。
HTH，
洛伦佐