function y = imsmooth(x, h) x = double(x); if nargin == 1, % smoothing filter. sum(sum(h)) = 1 h = 1/16*[1 2 1; 2 4 2; 1 2 1]; end y = conv2(x, h, 'same'); y = uint8(y);