Read netCDFs in loop and export to table with different sizes

2 次查看（过去 30 天）
Siegmund Nuyts
Siegmund Nuyts 2022-9-23
评论： Siegmund Nuyts 2022-9-27
I am reading several netCDF files with a datetime and corresponding data.
All datafiles have a different lenght, so I thought to find the common range of datetime for all files and save the corresponding data in 1 table so that I have a table like:
[datetime sensor1 sensor2 sensor3 sensor...]
I already found a way to read the files and find the common date but struggle to get evertyhing into 1 table.
Does anyone has a suggestion?
addpath(genpath('./'))
path = 'F:\Waterlevel/';
netCDF = dir([path,strcat('**/Sensor_*')]);
for i = 1:length(netCDF)
baseFilename = netCDF(i).name;
fullFileName = fullfile(path, baseFilename);
fprintf(1, 'Now reading %s\n', fullFileName);
nt = ncread(fullFileName, 'datetime');
Times = datetime(nt, 'convertFrom', 'posixtime', 'Format', 'yyyy-MM-dd HH:mm:ss');
nlevel = ncread(fullFileName, 'level');
[commondates, whereina, whereinb] = intersect(Times,Times(1:length(i)));
firstcommon = commondates(1);
end

请先登录，再回答此问题。

采纳的回答

Eric Sofen
Eric Sofen 2022-9-23
编辑：Eric Sofen 2022-9-26
I would suggest using timetable and the synchronize method.
I assume nlevel is the sensor data? If so, build up a cell array of timetables in the loop, then synchronize them.
Edit: fixed concatenation code.
addpath(genpath('./'))
path = 'F:\Waterlevel/';
netCDF = dir([path,'**/Sensor_*']);
alldata = {};
for i = 1:length(netCDF)
baseFilename = netCDF(i).name;
fullFileName = fullfile(path, baseFilename);
fprintf(1, 'Now reading %s\n', fullFileName);
nt = ncread(fullFileName, 'datetime');
Times = datetime(nt, 'convertFrom', 'posixtime', 'Format', 'yyyy-MM-dd HH:mm:ss');
nlevel = ncread(fullFileName, 'level');
alldata = [alldata, {timetable(Times,nlevel)}];
end
finaldata = synchronize(alldata{:},'intersection');
  4 个评论
显示 2更早的评论
Eric Sofen
Eric Sofen 2022-9-26
Whoops! alldata should be built up as a cell array of tables:
alldata = [alldata, {timetable(Times, nlevel)}];

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 NetCDF 的更多信息

标签

产品


版本

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by