function dim = decimate(image) 
% subsample 1/4 of the image pixels


s = size(image);
dim = image(1:2:s(1), 1:2:s(2));