我已经在这方面工作了一段时间...我使用ncread(path,myvarname,[1 1 1],[]读取的每个netcdf文件都有684x447x72形式的netcdf数据684 447 72])。我需要在文件= 2的维度3中附加数据,以使变量现在为684x447x144,依此类推每个新文件。所以,例如,如果我有10个文件,那么最后的read_in_var = 684x447x720。如何在一次读取一个文件的循环中执行此操作?谢谢!!
答案 0 :(得分:1)
使用matrix indexing。 考虑到文件名是file1,file2,file3,... file10,以下内容将起作用:
read_in_var = zeros(684,447,720); %Pre-allocation
for k=1:10
read_in_var(:,:,72*(k-1)+1:72*k)) = ncread(['file',num2str(k),'.nc'], ...
myvarname,[1 1 1], [684 447 72]);
end