%% load data
fid = fopen('fahrradwetterstatistik.dat');
fdata = cell2mat(textscan(fid, '%n%n%n','commentStyle', '//'));
fclose(fid)

clear fid


%% make overall pie chart
all = fdata(:,2:3);
figure
plotWeatherPie(all(:),1,1)

clear all


%% morning / evening pie chart
% find first entries for each day only
dif = diff(fdata(:,1));
firsts = fdata(logical([1;dif>0]),2:3);

fig = figure('Position',[144   204   748   364],...
             'PaperPositionMode','auto');

subplot(1,2,1)
plotWeatherPie(firsts(:,1),0,1)
text(0.5,1.1,'morning','Units','normalized','HorizontalAlignment','center')

subplot(1,2,2)
plotWeatherPie(firsts(:,2),1)
text(0.5,1.1,'evening','Units','normalized','HorizontalAlignment','center')

saveas(fig,'e:\website\images\fahrradwetterstatistik','png')

clear dif firsts fig


%% sequence plot
dif = diff(fdata(:,1));
firsts = fdata(logical([1;dif>0]),:);

dayes = datenum(num2str(firsts(:,1)),'yyyymmdd');
firstday = dayes(1);
dayes = dayes - dayes(1) + 1;

ndps = dayes(end)-29;
N = zeros(ndps,1);
percents = zeros(4,2,ndps);
for i = 1:ndps
    ind = find(dayes>=i & dayes<i+29);
    if isempty(ind)
        percents(:,:,i) = percents(:,:,i-1);
        N(i) = 0;
    else
        nums = histc(firsts(ind,2:3),-1:2);
        percents(:,:,i) = nums./repmat(sum(nums),4,1);
        N(i) = length(ind);
    end
end

fig = figure('Position',[144   204   748   364],...
             'PaperPositionMode','auto');

ticks = [1:round(ndps/10):ndps];
ny2007 = datenum('20070101','yyyymmdd')-firstday+1-14-.5;
ny2008 = datenum('20080101','yyyymmdd')-firstday+1-14-.5;
ny2009 = datenum('20090101','yyyymmdd')-firstday+1-14-.5;
ny2010 = datenum('20100101','yyyymmdd')-firstday+1-14-.5;

subplot(2,1,1)
plotWeatherArea(squeeze(percents(:,1,:))',N,1)
title('morning')
set(gca,'Layer','top')
set(gca,'XTick',ticks,'XTickLabel',{})
hold on
plot([ny2007,ny2007],[0,1],'w')
plot([ny2008,ny2008],[0,1],'w')
plot([ny2009,ny2009],[0,1],'w')
plot([ny2010,ny2010],[0,1],'w')
hold off

subplot(2,1,2)
plotWeatherArea(squeeze(percents(:,2,:))',N)
title('evening')
xlabel('date')
set(gca,'Layer','top')
set(gca,'XTick',ticks,'XTickLabel',datestr(firstday+(ticks-1)+14,'dd/mm'))
hold on
plot([ny2007,ny2007],[0,1],'w')
plot([ny2008,ny2008],[0,1],'w')
plot([ny2009,ny2009],[0,1],'w')
plot([ny2010,ny2010],[0,1],'w')
hold off
text(ny2007-0.5,1.06,'2006','Units','data','HorizontalAlignment','right')
text(ny2007+0.5,1.06,'2007','Units','data','HorizontalAlignment','left')
text(ny2008-0.5,1.06,'2007','Units','data','HorizontalAlignment','right')
text(ny2008+0.5,1.06,'2008','Units','data','HorizontalAlignment','left')
text(ny2009-0.5,1.06,'2008','Units','data','HorizontalAlignment','right')
text(ny2009+0.5,1.06,'2009','Units','data','HorizontalAlignment','left')
text(ny2010-0.5,1.06,'2009','Units','data','HorizontalAlignment','right')
text(ny2010+0.5,1.06,'2010','Units','data','HorizontalAlignment','left')
ylim([0,1])

saveas(fig,'e:\website\images\fahrradwetterzeitreihe','png')

clear dayes ind i nums fig N percents ndps ny2007 firstday ticks firsts dif