clear all
close all
clc

load('MAT\combi_vol');
load('MAT\rtk_vol');
load('MAT\sb25m_vol');
load('MAT\sb50m_vol');

figure
plot(rtk_vol(:,1),rtk_vol(:,2),'r^','MarkerFaceColor','r')
hold on
grid on
plot(sb25m_vol(:,1),sb25m_vol(:,2),'bo','MarkerFaceColor','w')
plot(sb50m_vol(:,1),sb50m_vol(:,2),'bo','MarkerFaceColor','b')
plot(combi_vol(:,1),combi_vol(:,2),'kx','MarkerSize',9)
hold on
grid on
startDate = datenum('10-01-2008');
endDate = datenum('12-31-2009');
xData = linspace(startDate,endDate,15);
set(gca,'XTick',xData);
datetick('x','mmm','keepticks')    
xlabel('datum')
ylabel('\DeltaV (m^3)')
title('Volumeverandering Galgeplaat suppletie t.o.v. T0')
legend('RTK','singlebeam 25 m','singlebeam 50 m','gecombineerde meting',4)
ylim([10 14]*10^4)
saveas(gcf,'PNG\kubering.png')
saveas(gcf,'EPS\kubering.eps','psc2')
