% t = jarkus_transects('areacode',7,'id',[7003575:7004000]);
% 
% t_m = jarkus_merge(t);
% t_mi = jarkus_interpolatenans(t_m, 'dim', 2);
% 
figure;
pcolor(t.x, t.y, squeeze(t_mi.altitude(:,:)));
shading interp; axis equal; hold on;

stride = 1;

load('D:\xbeach\egmond\GIS2\elevations.mat');

x = info.xllcorner+[0:stride:info.ncols-1]*info.cellsize;
y = info.yllcorner+[0:stride:info.nrows-1]*info.cellsize;
[X Y] = meshgrid(x, y);

A = A / 100;
pcolor(X, Y, A(1:stride:end, 1:stride:end));
shading flat; axis equal; hold on;

% load('D:\xbeach\egmond\GIS2\houses.mat');
% 
% x = info.xllcorner+[0:stride:info.ncols-1]*info.cellsize;
% y = info.yllcorner+[0:stride:info.nrows-1]*info.cellsize;
% [X Y] = meshgrid(x, y);
% 
% A = A / 100;
% pcolor(X, Y, A(1:stride:end, 1:stride:end));
% shading flat; axis equal; hold on;
% 
% load('D:\xbeach\egmond\GIS2\roads.mat');
% 
% x = info.xllcorner+[0:stride:info.ncols-1]*info.cellsize;
% y = info.yllcorner+[0:stride:info.nrows-1]*info.cellsize;
% [X Y] = meshgrid(x, y);
% 
% A = A / 100;
% pcolor(X, Y, A(1:stride:end, 1:stride:end));
% shading flat; axis equal; hold on;

colorbar;