Steps taken to break apart a timetable by years and plot in stackedplot
- The timetable is broken apart by year using arrayfun and each subtable is stored in a cell array TTbyYear.
- The dates of each subtable are shifted to 1952, preserving the month and day. Note that 1952 is a leapyear. It's important that all dates share the same year and since some years in your data may be leapyears, the shared year should be a leapyear.
- The subtables are joined back together using outerjoin where each year has its own column and the variable names in the new table TTy are the original years.
- The new table TTy is sent to stackedplot().
Demo
% Create demo timetable
TT = timetable(round(rand(3652,1)*100)/10-5,...
'StartTime',datetime(1956,1,1,'Format','dd.MM.yyyy'),...
'TimeStep',days(1),...
'VariableNames',{'DailyAvgTemp'});
TT.Properties.DimensionNames{1} = 'Datum';
head(TT) % show sample
% Break apart timetable by year
TTyears = year(TT.Datum);
unqYears = unique(TTyears);
TTbyYear = arrayfun(@(y){TT(TTyears==y,:)},unqYears);
% Join the tables back into a timetable with 1 column for each year
baseYear = 1956; % should be a leap year
TTy = timetable();
for i = 1:numel(TTbyYear)
% Change year to base year
TTbyYear{i}.Datum = datetime(baseYear,month(TTbyYear{i}.Datum),day(TTbyYear{i}.Datum));
% join tables
TTy = outerjoin(TTy, TTbyYear{i});
end
% Change variable names to years
TTy.Properties.VariableNames = compose('%d',unqYears);
% Change format of datetime values to hide the year
TTy.Datum.Format = 'dd.MM';
head(TTy) % Show sample
% Plot
stackedplot(TTy)
