% läs in bilden
f = double(imread('strand.png'));

% fx = df/dx, fy = df/dy
[fx,fy] = gradient(f);

% Genom att räkna ut skalärproduken av riktningsvektorn och gradienten får
% man ut riktningsderivatan
g = (fx * 2 + fy * -1) / sqrt(5);

% visa bilden (vi är bara intresserade av absolutvärdet)
imagesc(abs(g))
colormap gray