function [Y,context] = compute(X,parameter,context)

% EROSION - compute

SE = create_se(parameter);

Y = morph_erode(X,SE);